Testing the Toxicofera: comparative transcriptomics casts doubt on the single, early evolution of the reptile venom system

Proteomic Analysis of Heloderma horridum horridum Venom: Assessment to Its Transcriptome and Newfound Proteins

Heloderma horridum horridum, a venomous reptile native to America, has a venom with potential applications in treating type II diabetes. In this work, H. h. horridum venom was extracted, lyophilized, and characterized using enzymatic assays for hyaluronidase, phospholipase, and protease. Proteomic analysis of the venom was conducted employing bottom-up/shotgun approaches, SDS-PAGE, high-pH reversed-phase chromatography, and fractionation of tryptic peptides using nano-LC-MS/MS. The proteins found in H. h. horridum venom were reviewed according to the classification of the transcriptome previously reported. The proteomic approach identified 101 enzymes, 36 other proteins, 15 protein inhibitors, 11 host defense proteins, and 1 toxin, including novel venom components such as calcium-binding proteins, phospholipase A2 inhibitors, serpins, cathepsin, subtilases, carboxypeptidase-like, aminopeptidases, glycoside hydrolases, thioredoxin transferases, acid ceramidase-like, enolase, multicopper oxidases, phosphoglucose isomerase (PGI), fructose-1,6-bisphosphatase class 1, pentraxin-related, peptidylglycine α-hydroxylating monooxygenase/peptidyl-hydroxyglycine α-amidating lyase, carbonic anhydrase, acetylcholinesterase, dipeptidylpeptidase, and lysozymes. These findings contribute to understanding the venomous nature of H. h. horridum and highlight its potential as a source of bioactive compounds. Data are available via PRoteomeXchange with the identifier PXD052417.

Mandibular gland proteomics of the Mexican alligator lizard, Abronia graminea, and the red-lipped arboreal alligator lizard, Abronia lythrochila

A useful approach to deepen our knowledge about the origin and evolution of venom systems in Reptilia has been exploring the vast biodiversity of this clade of vertebrates in search of orally produced proteins with toxic actions, as well as their corresponding delivery systems. The occurrence of toxins in anguimorph lizards has been demonstrated experimentally or inferred from reports of the toxic effects of the oral secretions of taxa within the Varanidae and Helodermatidae families. In the present study, we have focused on two alligator lizards of the Anguidae family, the Mexican alligator lizard, Abronia graminea, and the red-lipped arboreal alligator lizard, A. lythrochila. In addition, the fine morphology of teeth of the latter species is described. The presence of a conserved set of proteins, including B-type natriuretic peptides, cysteine-rich secretory proteins, group III phospholipase A2, and kallikrein, in submandibular gland extracts was demonstrated for both Abronia species. These proteins belong to toxin families found in oral gland secretions of venomous reptile species. This finding, along with previous demonstration of toxin-producing taxa in both paleo- and neoanguimorpha clades, provides further support for the existence of a handful of conserved toxin families in oral secretions across the 100+ million years of Anguimorpha cladogenesis.

Red-on-Yellow Queen: Bio-Layer Interferometry Reveals Functional Diversity Within Micrurus Venoms and Toxin Resistance in Prey Species

Snakes in the family Elapidae largely produce venoms rich in three-finger toxins (3FTx) that bind to the α 1 subunit of nicotinic acetylcholine receptors (nAChRs), impeding ion channel activity. These neurotoxins immobilize the prey by disrupting muscle contraction. Coral snakes of the genus Micrurus are specialist predators who produce many 3FTx, making them an interesting system for examining the coevolution of these toxins and their targets in prey animals. We used a bio-layer interferometry technique to measure the binding interaction between 15 Micrurus venoms and 12 taxon-specific mimotopes designed to resemble the orthosteric binding region of the muscular nAChR subunit. We found that Micrurus venoms vary greatly in their potency on this assay and that this variation follows phylogenetic patterns rather than previously reported patterns of venom composition. The long-tailed Micrurus tend to have greater binding to nAChR orthosteric sites than their short-tailed relatives and we conclude this is the likely ancestral state. The repeated loss of this activity may be due to the evolution of 3FTx that bind to other regions of the nAChR. We also observed variations in the potency of the venoms depending on the taxon of the target mimotope. Rather than a pattern of prey-specificity, we found that mimotopes modeled after snake nAChRs are less susceptible to Micrurus venoms and that this resistance is partly due to a characteristic tryptophan → serine mutation within the orthosteric site in all snake mimotopes. This resistance may be part of a Red Queen arms race between coral snakes and their prey.

Distinct regulatory networks control toxin gene expression in elapid and viperid snakes

BACKGROUND

Venom systems are ideal models to study genetic regulatory mechanisms that underpin evolutionary novelty. Snake venom glands are thought to share a common origin, but there are major distinctions between venom toxins from the medically significant snake families Elapidae and Viperidae, and toxin gene regulatory investigations in elapid snakes have been limited. Here, we used high-throughput RNA-sequencing to profile gene expression and microRNAs between active (milked) and resting (unmilked) venom glands in an elapid (Eastern Brown Snake, Pseudonaja textilis), in addition to comparative genomics, to identify cis- and trans-acting regulation of venom production in an elapid in comparison to viperids (Crotalus viridis and C. tigris).

RESULTS

Although there is conservation in high-level mechanistic pathways regulating venom production (unfolded protein response, Notch signaling and cholesterol homeostasis), there are differences in the regulation of histone methylation enzymes, transcription factors, and microRNAs in venom glands from these two snake families. Histone methyltransferases and transcription factor (TF) specificity protein 1 (Sp1) were highly upregulated in the milked elapid venom gland in comparison to the viperids, whereas nuclear factor I (NFI) TFs were upregulated after viperid venom milking. Sp1 and NFI cis-regulatory elements were common to toxin gene promoter regions, but many unique elements were also present between elapid and viperid toxins. The presence of Sp1 binding sites across multiple elapid toxin gene promoter regions that have been experimentally determined to regulate expression, in addition to upregulation of Sp1 after venom milking, suggests this transcription factor is involved in elapid toxin expression. microRNA profiles were distinctive between milked and unmilked venom glands for both snake families, and microRNAs were predicted to target a diversity of toxin transcripts in the elapid P. textilis venom gland, but only snake venom metalloproteinase transcripts in the viperid C. viridis venom gland. These results suggest differences in toxin gene posttranscriptional regulation between the elapid P. textilis and viperid C. viridis.

CONCLUSIONS

Our comparative transcriptomic and genomic analyses between toxin genes and isoforms in elapid and viperid snakes suggests independent toxin regulation between these two snake families, demonstrating multiple different regulatory mechanisms underpin a venomous phenotype.

Domain loss enabled evolution of novel functions in the snake three-finger toxin gene superfamily

Three-finger toxins (3FTXs) are a functionally diverse family of toxins, apparently unique to venoms of caenophidian snakes. Although the ancestral function of 3FTXs is antagonism of nicotinic acetylcholine receptors, redundancy conferred by the accumulation of duplicate genes has facilitated extensive neofunctionalization, such that derived members of the family interact with a range of targets. 3FTXs are members of the LY6/UPAR family, but their non-toxin ancestor remains unknown. Combining traditional phylogenetic approaches, manual synteny analysis, and machine learning techniques (including AlphaFold2 and ProtT5), we have reconstructed a detailed evolutionary history of 3FTXs. We identify their immediate ancestor as a non-secretory LY6, unique to squamate reptiles, and propose that changes in molecular ecology resulting from loss of a membrane-anchoring domain and changes in gene expression, paved the way for the evolution of one of the most important families of snake toxins.

Expansion and collapse of VEGF diversity in major clades of the animal kingdom

Together with the platelet-derived growth factors (PDGFs), the vascular endothelial growth factors (VEGFs) form the PDGF/VEGF subgroup among cystine knot growth factors. The evolutionary relationships within this subgroup have not been examined thoroughly to date. Here, we comprehensively analyze the PDGF/VEGF growth factors throughout all animal phyla and propose a phylogenetic tree. Vertebrate whole-genome duplications play a role in expanding PDGF/VEGF diversity, but several limited duplications are necessary to account for the temporal pattern of emergence. The phylogenetically oldest PDGF/VEGF-like growth factor likely featured a C-terminus with a BR3P signature, a hallmark of the modern-day lymphangiogenic growth factors VEGF-C and VEGF-D. Some younger VEGF genes, such as VEGFB and PGF, appeared completely absent in important vertebrate clades such as birds and amphibia, respectively. In contrast, individual PDGF/VEGF gene duplications frequently occurred in fish on top of the known fish-specific whole-genome duplications. The lack of precise counterparts for human genes poses limitations but also offers opportunities for research using organisms that diverge considerably from humans. Sources for the graphical abstract: 326 MYA and older [1]; 72-240 MYA [2]; 235-65 MYA [3].

Proteomic analysis of the mandibular glands from the Chinese crocodile lizard, Shinisaurus crocodilurus - Another venomous lizard?

Based on its phylogenetic relationship to monitor lizards (Varanidae), Gila monsters (Heloderma spp.), and the earless monitor Lanthanotus borneesis, the Chinese crocodile lizard, Shinisaurus crocodilurus, has been assigned to the Toxicofera clade, which comprises venomous reptiles. However, no data about composition and biological activities of its oral secretion have been reported. In the present study, a proteomic analysis of the mandibular gland of S. crocodilurus and, for comparison, of the herbivorous Solomon Island skink Corucia zebrata, was performed. Scanning electron microscopy (SEM) of the teeth from S. crocodilurus revealed a sharp ridge on the anterior surface, but no grooves, whereas those of C. zebrata possess a flattened crown with a pointed cusp. Proteomic analysis of their gland extracts provided no evidence of venom-derived peptides or proteins, strongly supporting the non-venomous character of these lizards. Data are available via ProteomeXchange with identifier PXD039424.

Snake Venomics: Fundamentals, Recent Updates, and a Look to the Next Decade

Venomic research, powered by techniques adapted from proteomics, transcriptomics, and genomics, seeks to unravel the diversity and complexity of venom through which knowledge can be applied in the treatment of envenoming, biodiscovery, and conservation. Snake venom proteomics is most extensively studied, but the methods varied widely, creating a massive amount of information which complicates data comparison and interpretation. Advancement in mass spectrometry technology, accompanied by growing databases and sophisticated bioinformatic tools, has overcome earlier limitations of protein identification. The progress, however, remains challenged by limited accessibility to samples, non-standardized quantitative methods, and biased interpretation of -omic data. Next-generation sequencing (NGS) technologies enable high-throughput venom-gland transcriptomics and genomics, complementing venom proteomics by providing deeper insights into the structural diversity, differential expression, regulation and functional interaction of the toxin genes. Venomic tissue sampling is, however, difficult due to strict regulations on wildlife use and transfer of biological materials in some countries. Limited resources for techniques and funding are among other pertinent issues that impede the progress of venomics, particularly in less developed regions and for neglected species. Genuine collaboration between international researchers, due recognition of regional experts by global organizations (e.g., WHO), and improved distribution of research support, should be embraced.

Co-option of the same ancestral gene family gave rise to mammalian and reptilian toxins

Background

Evolution can occur with surprising predictability when organisms face similar ecological challenges. For most traits, it is difficult to ascertain whether this occurs due to constraints imposed by the number of possible phenotypic solutions or because of parallel responses by shared genetic and regulatory architecture. Exceptionally, oral venoms are a tractable model of trait evolution, being largely composed of proteinaceous toxins that have evolved in many tetrapods, ranging from reptiles to mammals. Given the diversity of venomous lineages, they are believed to have evolved convergently, even though biochemically similar toxins occur in all taxa.

Results

Here, we investigate whether ancestral genes harbouring similar biochemical activity may have primed venom evolution, focusing on the origins of kallikrein-like serine proteases that form the core of most vertebrate oral venoms. Using syntenic relationships between genes flanking known toxins, we traced the origin of kallikreins to a single locus containing one or more nearby paralogous kallikrein-like clusters. Additionally, phylogenetic analysis of vertebrate serine proteases revealed that kallikrein-like toxins in mammals and reptiles are genetically distinct from non-toxin ones.

Conclusions

Given the shared regulatory and genetic machinery, these findings suggest that tetrapod venoms evolved by co-option of proteins that were likely already present in saliva. We term such genes ‘toxipotent’—in the case of salivary kallikreins they already had potent vasodilatory activity that was weaponized by venomous lineages. Furthermore, the ubiquitous distribution of kallikreins across vertebrates suggests that the evolution of envenomation may be more common than previously recognized, blurring the line between venomous and non-venomous animals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01191-1.

Utilisation of compounds from venoms in drug discovery

Difficult drug targets are becoming the normal course of business in drug discovery, sometimes due to large interacting surfaces or only small differences in selectivity regions. For these, a different approach is merited: compounds lying somewhere between the small molecule and the large antibody in terms of many properties including stability, biodistribution and pharmacokinetics. Venoms have evolved over millions of years to be complex mixtures of stable molecules derived from other somatic molecules, the stability comes from the pressure to be ready for delivery at a moment's notice. Snakes, spiders, scorpions, jellyfish, wasps, fish and even mammals have evolved independent venom systems with complex mixtures in their chemical arsenal. These venom-derived molecules have been proven to be useful tools, such as for the development of antihypotensive angiotensin converting enzyme (ACE) inhibitors and have also made successful drugs such as Byetta® (Exenatide), Integrilin® (Eptifibatide) and Echistatin. Only a small percentage of the available chemical space from venoms has been investigated so far and this is growing. In a new era of biological therapeutics, venom peptides present opportunities for larger target engagement surface with greater stability than antibodies or human peptides. There are challenges for oral absorption and target engagement, but there are venom structures that overcome these and thus provide substrate for engineering novel molecules that combine all desired properties. Venom researchers are characterising new venoms, species, and functions all the time, these provide great substrate for solving the challenges presented by today's difficult targets.

Old World Vipers-A Review about Snake Venom Proteomics of Viperinae and Their Variations

Fine-tuned by millions of years of evolution, snake venoms have frightened but also fascinated humanity and nowadays they constitute potential resources for drug development, therapeutics and antivenoms. The continuous progress of mass spectrometry techniques and latest advances in proteomics workflows enabled toxinologists to decipher venoms by modern omics technologies, so-called ‘venomics’. A tremendous upsurge reporting on snake venom proteomes could be observed. Within this review we focus on the highly venomous and widely distributed subfamily of Viperinae (Serpentes: Viperidae). A detailed public literature database search was performed (2003–2020) and we extensively reviewed all compositional venom studies of the so-called Old-World Vipers. In total, 54 studies resulted in 89 venom proteomes. The Viperinae venoms are dominated by four major, four secondary, six minor and several rare toxin families and peptides, respectively. The multitude of different venomics approaches complicates the comparison of venom composition datasets and therefore we differentiated between non-quantitative and three groups of quantitative workflows. The resulting direct comparisons within these groups show remarkable differences on the intra- and interspecies level across genera with a focus on regional differences. In summary, the present compilation is the first comprehensive up-to-date database on Viperinae venom proteomes and differentiating between analytical methods and workflows.

Tracking the recruitment and evolution of snake toxins using the evolutionary context provided by the Bothrops jararaca genome

Venom is a key adaptive innovation in snakes, and how nonvenom genes were co-opted to become part of the toxin arsenal is a significant evolutionary question. While this process has been investigated through the phylogenetic reconstruction of toxin sequences, evidence provided by the genomic context of toxin genes remains less explored. To investigate the process of toxin recruitment, we sequenced the genome of Bothrops jararaca, a clinically relevant pitviper. In addition to producing a road map with canonical structures of genes encoding 12 toxin families, we inferred most of the ancestral genes for their loci. We found evidence that 1) snake venom metalloproteinases (SVMPs) and phospholipases A2 (PLA2) have expanded in genomic proximity to their nonvenomous ancestors; 2) serine proteinases arose by co-opting a local gene that also gave rise to lizard gilatoxins and then expanded; 3) the bradykinin-potentiating peptides originated from a C-type natriuretic peptide gene backbone; and 4) VEGF-F was co-opted from a PGF-like gene and not from VEGF-A. We evaluated two scenarios for the original recruitment of nontoxin genes for snake venom: 1) in locus ancestral gene duplication and 2) in locus ancestral gene direct co-option. The first explains the origins of two important toxins (SVMP and PLA2), while the second explains the emergence of a greater number of venom components. Overall, our results support the idea of a locally assembled venom arsenal in which the most clinically relevant toxin families expanded through posterior gene duplications, regardless of whether they originated by duplication or gene co-option.

Recruitment of toxin-like proteins with ancestral venom function supports endoparasitic lifestyles of Myxozoa

Cnidarians are the oldest lineage of venomous animals and use nematocysts to discharge toxins. Whether venom toxins have been recruited to support parasitic lifestyles in the Endocnidozoa (Myxozoa + Polypodium) is, however, unknown. To examine this issue we variously employed transcriptomic, proteomic, associated molecular phylogenies, and localisation studies on representative primitive and derived myxozoans (Malacosporea and Myxosporea, respectively), Polypodium hydriforme, and the free-living staurozoan Calvadosia cruxmelitensis. Our transcriptomics and proteomics analyses provide evidence for expression and translation of venom toxin homologs in myxozoans. Phylogenetic placement of Kunitz type serine protease inhibitors and phospholipase A2 enzymes reveals modification of toxins inherited from ancestral free-living cnidarian toxins, and that venom diversity is reduced in myxozoans concordant with their reduced genome sizes. Various phylogenetic analyses of the Kunitz-type toxin family in Endocnidozoa suggested lineage-specific gene duplications, which offers a possible mechanism for enhancing toxin diversification. Toxin localisation in the malacosporean Buddenbrockia plumatellae substantiates toxin translation and thus illustrates a repurposing of toxin function for endoparasite development and interactions with hosts, rather than for prey capture or defence. Whether myxozoan venom candidates are expressed in transmission stages (e.g. in nematocysts or secretory vesicles) requires further investigation.

The Curious Case of the "Neurotoxic Skink": Scientific Literature Points to the Absence of Venom in Scincidae

In contrast to the clearly documented evolution of venom in many animal lineages, the origin of reptilian venom is highly debated. Historically, venom has been theorised to have evolved independently in snakes and lizards. However, some of the recent works have argued for the common origin of venom in “Toxicofera” reptiles, which include the order Serpentes (all snakes), and Anguimorpha and Iguania lizards. Nevertheless, in both these contrasting hypotheses, the lizards of the family Scincidae are considered to be harmless and devoid of toxic venoms. Interestingly, an unusual clinical case claiming neurotoxic envenoming by a scincid lizard was recently reported in Southern India. Considering its potentially significant medicolegal, conservation and evolutionary implications, we have summarised the scientific evidence that questions the validity of this clinical report. We argue that the symptoms documented in the patient are likely to have resulted from krait envenomation, which is far too frequent in these regions.

Convergent evolution of pain-inducing defensive venom components in spitting cobras

Convergent evolution provides insights into the selective drivers underlying evolutionary change. Snake venoms, with a direct genetic basis and clearly defined functional phenotype, provide a model system for exploring the repeated evolution of adaptations. While snakes use venom primarily for predation, and venom composition often reflects diet specificity, three lineages of cobras have independently evolved the ability to spit venom at adversaries. Using gene, protein, and functional analyses, we show that the three spitting lineages possess venoms characterized by an up-regulation of phospholipase A2 (PLA2) toxins, which potentiate the action of preexisting venom cytotoxins to activate mammalian sensory neurons and cause enhanced pain. These repeated independent changes provide a fascinating example of convergent evolution across multiple phenotypic levels driven by selection for defense.

In-depth transcriptome reveals the potential biotechnological application of Bothrops jararaca venom gland

Background:

Lack of complete genomic data of Bothrops jararaca impedes molecular biology research focusing on biotechnological applications of venom gland components. Identification of full-length coding regions of genes is crucial for the correct molecular cloning design.

Methods:

RNA was extracted from the venom gland of one adult female specimen of Bothrops jararaca. Deep sequencing of the mRNA library was performed using Illumina NextSeq 500 platform. De novo assembly of B. jararaca transcriptome was done using Trinity. Annotation was performed using Blast2GO. All predicted proteins after clustering step were blasted against non-redundant protein database of NCBI using BLASTP. Metabolic pathways present in the transcriptome were annotated using the KAAS-KEGG Automatic Annotation Server. Toxins were identified in the B. jararaca predicted proteome using BLASTP against all protein sequences obtained from Animal Toxin Annotation Project from Uniprot KB/Swiss-Pro database. Figures and data visualization were performed using ggplot2 package in R language environment.

Results:

We described the in-depth transcriptome analysis of B. jararaca venom gland, in which 76,765 de novo assembled isoforms, 96,044 transcribed genes and 41,196 unique proteins were identified. The most abundant transcript was the zinc metalloproteinase-disintegrin-like jararhagin. Moreover, we identified 78 distinct functional classes of proteins, including toxins, inhibitors and tumor suppressors. Other venom proteins identified were the hemolytic lethal factors stonustoxin and verrucotoxin.

Conclusion:

It is believed that the application of deep sequencing to the analysis of snake venom transcriptomes may represent invaluable insight on their biotechnological potential focusing on candidate molecules.

Genome mapping of a LYST mutation in corn snakes indicates that vertebrate chromatophore vesicles are lysosome-related organelles

Reptiles exhibit a spectacular diversity of skin colors and patterns brought about by the interactions among three chromatophore types: black melanophores with melanin-packed melanosomes, red and yellow xanthophores with pteridine- and/or carotenoid-containing vesicles, and iridophores filled with light-reflecting platelets generating structural colors. Whereas the melanosome, the only color-producing endosome in mammals and birds, has been documented as a lysosome-related organelle, the maturation paths of xanthosomes and iridosomes are unknown. Here, we first use 10x Genomics linked-reads and optical mapping to assemble and annotate a nearly chromosome-quality genome of the corn snake Pantherophis guttatus The assembly is 1.71 Gb long, with an N50 of 16.8 Mb and L50 of 24. Second, we perform mapping-by-sequencing analyses and identify a 3.9-Mb genomic interval where the lavender variant resides. The lavender color morph in corn snakes is characterized by gray, rather than red, blotches on a pink, instead of orange, background. Third, our sequencing analyses reveal a single nucleotide polymorphism introducing a premature stop codon in the lysosomal trafficking regulator gene (LYST) that shortens the corresponding protein by 603 amino acids and removes evolutionary-conserved domains. Fourth, we use light and transmission electron microscopy comparative analyses of wild type versus lavender corn snakes and show that the color-producing endosomes of all chromatophores are substantially affected in the LYST mutant. Our work provides evidence characterizing xanthosomes in xanthophores and iridosomes in iridophores as lysosome-related organelles.

Deimination Protein Profiles in Alligator mississippiensis Reveal Plasma and Extracellular Vesicle-Specific Signatures Relating to Immunity, Metabolic Function, and Gene Regulation

Alligators are crocodilians and among few species that endured the Cretaceous-Paleogene extinction event. With long life spans, low metabolic rates, unusual immunological characteristics, including strong antibacterial and antiviral ability, and cancer resistance, crocodilians may hold information for molecular pathways underlying such physiological traits. Peptidylarginine deiminases (PADs) are a group of calcium-activated enzymes that cause posttranslational protein deimination/citrullination in a range of target proteins contributing to protein moonlighting functions in health and disease. PADs are phylogenetically conserved and are also a key regulator of extracellular vesicle (EV) release, a critical part of cellular communication. As little is known about PAD-mediated mechanisms in reptile immunology, this study was aimed at profiling EVs and protein deimination in Alligator mississippiensis. Alligator plasma EVs were found to be polydispersed in a 50-400-nm size range. Key immune, metabolic, and gene regulatory proteins were identified to be posttranslationally deiminated in plasma and plasma EVs, with some overlapping hits, while some were unique to either plasma or plasma EVs. In whole plasma, 112 target proteins were identified to be deiminated, while 77 proteins were found as deiminated protein hits in plasma EVs, whereof 31 were specific for EVs only, including proteins specific for gene regulatory functions (e.g., histones). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed KEGG pathways specific to deiminated proteins in whole plasma related to adipocytokine signaling, while KEGG pathways of deiminated proteins specific to EVs included ribosome, biosynthesis of amino acids, and glycolysis/gluconeogenesis pathways as well as core histones. This highlights roles for EV-mediated export of deiminated protein cargo with roles in metabolism and gene regulation, also related to cancer. The identification of posttranslational deimination and EV-mediated communication in alligator plasma revealed here contributes to current understanding of protein moonlighting functions and EV-mediated communication in these ancient reptiles, providing novel insight into their unusual immune systems and physiological traits. In addition, our findings may shed light on pathways underlying cancer resistance, antibacterial and antiviral resistance, with translatable value to human pathologies.

Deciphering Ancestral Sex Chromosome Turnovers Based on Analysis of Male Mutation Bias

The age of sex chromosomes is commonly obtained by comparing the substitution rates of XY gametologs. Coupled with phylogenetic reconstructions, one can refine the origin of a sex chromosome system relative to specific speciation events. However, these approaches are insufficient to determine the presence and duration of ancestral sex chromosome systems that were lost in some species. In this study, we worked with genomic and transcriptomic data from mammals and squamates and analyzed the effect of male mutation bias on X-linked sequences in these groups. We searched for signatures indicating whether monotremes shared the same sex chromosomes with placental mammals or whether pleurodonts and acrodonts had a common ancestral sex chromosome system. Our analyses indicate that platypus did not share the XY chromosomes with placental mammals, in agreement with previous work. In contrast, analyses of agamids showed that this lineage maintained the pleurodont XY chromosomes for several million years. We performed multiple simulations using different strengths of male mutation bias to confirm the results. Overall, our work shows that variations in substitution rates due to male mutation bias could be applied to uncover signatures of ancestral sex chromosome systems.

Omics Technologies for Profiling Toxin Diversity and Evolution in Snake Venom: Impacts on the Discovery of Therapeutic and Diagnostic Agents

Snake venoms are primarily composed of proteins and peptides, and these toxins have developed high selectivity to their biological targets. This makes venoms interesting for exploration into protein evolution and structure-function relationships. A single venom protein superfamily can exhibit a variety of pharmacological effects; these variations in activity originate from differences in functional sites, domains, posttranslational modifications, and the formations of toxin complexes. In this review, we discuss examples of how the major venom protein superfamilies have diversified, as well as how newer technologies in the omics fields, such as genomics, transcriptomics, and proteomics, can be used to characterize both known and unknown toxins.Because toxins are bioactive molecules with a rich diversity of activities, they can be useful as therapeutic and diagnostic agents, and successful examples of toxin applications in these areas are also reviewed. With the current rapid pace of technology, snake venom research and its applications will only continue to expand.

Evolution: Fangtastic Venoms Underpin Parasitic Mimicry

Venomous teeth are rare in fishes, which typically utilise spines for defence. A new study reveals the evolutionary origins of fangs and venom in the Nemophini blennies and shows that, in contrast to snakes and lizards, the fangs pre-date the venom.

The Harderian gland transcriptomes of Caraiba andreae, Cubophis cantherigerus and Tretanorhinus variabilis, three colubroid snakes from Cuba

The Harderian gland is a cephalic structure, widely distributed among vertebrates. In snakes, the Harderian gland is anatomically connected to the vomeronasal organ via the nasolacrimal duct, and in some species can be larger than the eyes. The function of the Harderian gland remains elusive, but it has been proposed to play a role in the production of saliva, pheromones, thermoregulatory lipids and growth factors, among others. Here, we have profiled the transcriptomes of the Harderian glands of three non-front-fanged colubroid snakes from Cuba: Caraiba andreae (Cuban Lesser Racer); Cubophis cantherigerus (Cuban Racer); and Tretanorhinus variabilis (Caribbean Water Snake), using Illumina HiSeq2000 100 bp paired-end. In addition to ribosomal and non-characterized proteins, the most abundant transcripts encode putative transport/binding, lipocalin/lipocalin-like, and bactericidal/permeability-increasing-like proteins. Transcripts coding for putative canonical toxins described in venomous snakes were also identified. This transcriptional profile suggests a more complex function than previously recognized for this enigmatic organ.

Biochemical characterization of venom from Pseudoboa neuwiedii (Neuwied's false boa; Xenodontinae; Pseudoboini)

In this work, we examined the proteolytic and phospholipase A₂ (PLA₂) activities of venom from the opisthoglyphous colubrid Pseudoboa neuwiedii. Proteolytic activity (3 and 10 μg of venom) was comparable to that of Bothrops neuwiedii venom but less than Bothrops atrox. This activity was inhibited by EDTA and 1,10-phenanthroline but only slightly affected (≤30% inhibition) by PMSF and AEBSF, indicating it was mediated by snake venom metalloproteinases (SVMPs). The pH and temperature optima for proteolytic activity were 8.0 and 37 °C, respectively. The venom had no esterase activity, whereas PLA₂ activity was similar to B. atrox, greater than B. neuwiedii but less than B. jararacussu. SDS-PAGE revealed venom proteins >100 kDa, 45-70 kDa, 21-24 kDa and ~15 kDa, and mass spectrometry of protein bands revealed SVMPs, cysteine-rich secretory proteins (CRISPs) and PLA₂, but no serine proteinases. In gelatin zymography, the most active bands occurred at 65-68 kDa (seen with 0.05-0.25 μg of venom). Caseinolytic activity occurred at 50-66 kDa and was generally weaker than gelatinolytic activity. RP-HPLC of venom yielded 15 peaks, five of which showed gelatinolytic activity; peak 7 was the most active and apparently contained a P-III class SVMP. The venom showed α-fibrinogenase activity, without affecting the β and γ chains; this activity was inhibited by EDTA and 1,10-phenanthroline. The venom did not clot rat citrated plasma but reduced the rate and extent of coagulation after plasma recalcification. In conclusion, P. neuwiedii venom is highly proteolytic and could potentially affect coagulation in vivo by degrading fibrinogen via SVMPs.

Noise and biases in genomic data may underlie radically different hypotheses for the position of Iguania within Squamata

Squamate reptiles are a major component of vertebrate biodiversity whose crown-clade traces its origin to a narrow window of time in the Mesozoic during which the main subclades diverged in rapid succession. Deciphering phylogenetic relationships among these lineages has proven challenging given the conflicting signals provided by genomic and phenomic data. Most notably, the placement of Iguania has routinely differed between data sources, with morphological evidence supporting a sister relationship to the remaining squamates (Scleroglossa hypothesis) and molecular data favoring a highly nested position alongside snakes and anguimorphs (Toxicofera hypothesis). We provide novel insights by generating an expanded morphological dataset and exploring the presence of phylogenetic signal, noise, and biases in molecular data. Our analyses confirm the presence of strong conflicting signals for the position of Iguania between morphological and molecular datasets. However, we also find that molecular data behave highly erratically when inferring the deepest branches of the squamate tree, a consequence of limited phylogenetic signal to resolve this ancient radiation with confidence. This, in turn, seems to result from a rate of evolution that is too high for historical signals to survive to the present. Finally, we detect significant systematic biases, with iguanians and snakes sharing faster rates of molecular evolution and a similarly biased nucleotide composition. A combination of scant phylogenetic signal, high levels of noise, and the presence of systematic biases could result in the misplacement of Iguania. We regard this explanation to be at least as plausible as the complex scenario of convergence and reversals required for morphological data to be misleading. We further evaluate and discuss the utility of morphological data to resolve ancient radiations, as well as its impact in combined-evidence phylogenomic analyses, with results relevant for the assessment of evidence and conflict across the Tree of Life.

Snake Venom Extracellular vesicles (SVEVs) reveal wide molecular and functional proteome diversity

Proteins constitute almost 95% of snake venom's dry weight and are produced and released by venom glands in a solubilized form during a snake bite. These proteins are responsible for inducing several pharmacological effects aiming to immobilize and initiate the pre-digestion of the prey. This study shows that proteins can be secreted and confined in snake venom extracellular vesicles (SVEVs) presenting a size distribution between 50 nm and 500 nm. SVEVs isolated from lyophilized venoms collected from four different species of snakes (Agkistrodon contortrix contortrix, Crotalus atrox, Crotalus viridis and Crotalus cerberus oreganus) were analyzed by mass spectrometry-based proteomic, which allowed the identification of proteins belonging to eight main functional protein classes such as SVMPs, serine proteinases, PLA₂, LAAO, 5'nucleotidase, C-type lectin, CRISP and Disintegrin. Biochemical assays indicated that SVEVs are functionally active, showing high metalloproteinase and fibrinogenolytic activity besides being cytotoxic against HUVEC cells. Overall, this study comprehensively depicts the protein composition of SVEVs for the first time. In addition, the molecular function of some of the described proteins suggests a central role for SVEVs in the cytotoxicity of the snake venom and sheds new light in the envenomation process.

Venomix: a simple bioinformatic pipeline for identifying and characterizing toxin gene candidates from transcriptomic data

The advent of next-generation sequencing has resulted in transcriptome-based approaches to investigate functionally significant biological components in a variety of non-model organism. This has resulted in the area of “venomics”: a rapidly growing field using combined transcriptomic and proteomic datasets to characterize toxin diversity in a variety of venomous taxa. Ultimately, the transcriptomic portion of these analyses follows very similar pathways after transcriptome assembly often including candidate toxin identification using BLAST, expression level screening, protein sequence alignment, gene tree reconstruction, and characterization of potential toxin function. Here we describe the Python package Venomix, which streamlines these processes using common bioinformatic tools along with ToxProt, a publicly available annotated database comprised of characterized venom proteins. In this study, we use the Venomix pipeline to characterize candidate venom diversity in four phylogenetically distinct organisms, a cone snail (Conidae; Conus sponsalis), a snake (Viperidae; Echis coloratus), an ant (Formicidae; Tetramorium bicarinatum), and a scorpion (Scorpionidae; Urodacus yaschenkoi). Data on these organisms were sampled from public databases, with each original analysis using different approaches for transcriptome assembly, toxin identification, or gene expression quantification. Venomix recovered numerically more candidate toxin transcripts for three of the four transcriptomes than the original analyses and identified new toxin candidates. In summary, we show that the Venomix package is a useful tool to identify and characterize the diversity of toxin-like transcripts derived from transcriptomic datasets. Venomix is available at: https://bitbucket.org/JasonMacrander/Venomix/.

MipLAAO, a new L-amino acid oxidase from the redtail coral snake Micrurus mipartitus

L-amino acid oxidases (LAAOs) are ubiquitous enzymes in nature. Bioactivities described for these enzymes include apoptosis induction, edema formation, induction or inhibition of platelet aggregation, as well as antiviral, antiparasite, and antibacterial actions. With over 80 species, Micrurus snakes are the representatives of the Elapidae family in the New World. Although LAAOs in Micrurus venoms have been predicted by venom gland transcriptomic studies and detected in proteomic studies, no enzymes of this kind have been previously purified from their venoms. Earlier proteomic studies revealed that the venom of M. mipartitus from Colombia contains ∼4% of LAAO. This enzyme, here named MipLAAO, was isolated and biochemically and functionally characterized. The enzyme is found in monomeric form, with an isotope-averaged molecular mass of 59,100.6 Da, as determined by MALDI-TOF. Its oxidase activity shows substrate preference for hydrophobic amino acids, being optimal at pH 8.0. By nucleotide sequencing of venom gland cDNA of mRNA transcripts obtained from a single snake, six isoforms of MipLAAO with minor variations among them were retrieved. The deduced sequences present a mature chain of 483 amino acids, with a predicted pI of 8.9, and theoretical masses between 55,010.9 and 55,121.0 Da. The difference with experimentally observed mass is likely due to glycosylation, in agreement with the finding of three putative N-glycosylation sites in its amino acid sequence. A phylogenetic analysis of MmipLAAO placed this new enzyme within the clade of homologous proteins from elapid snakes, characterized by the conserved Serine at position 223, in contrast to LAAOs from viperids. MmipLAAO showed a potent bactericidal effect on S. aureus (MIC: 2 µg/mL), but not on E. coli. The former activity could be of interest to future studies assessing its potential as antimicrobial agent.

Enter the Dragon: The Dynamic and Multifunctional Evolution of Anguimorpha Lizard Venoms

While snake venoms have been the subject of intense study, comparatively little work has been done on lizard venoms. In this study, we have examined the structural and functional diversification of anguimorph lizard venoms and associated toxins, and related these results to dentition and predatory ecology. Venom composition was shown to be highly variable across the 20 species of Heloderma, Lanthanotus, and Varanus included in our study. While kallikrein enzymes were ubiquitous, they were also a particularly multifunctional toxin type, with differential activities on enzyme substrates and also ability to degrade alpha or beta chains of fibrinogen that reflects structural variability. Examination of other toxin types also revealed similar variability in their presence and activity levels. The high level of venom chemistry variation in varanid lizards compared to that of helodermatid lizards suggests that venom may be subject to different selection pressures in these two families. These results not only contribute to our understanding of venom evolution but also reveal anguimorph lizard venoms to be rich sources of novel bioactive molecules with potential as drug design and development lead compounds.

Coralsnake Venomics: Analyses of Venom Gland Transcriptomes and Proteomes of Six Brazilian Taxa

Venom gland transcriptomes and proteomes of six Micrurus taxa (M. corallinus, M. lemniscatus carvalhoi, M. lemniscatus lemniscatus, M. paraensis, M. spixii spixii, and M. surinamensis) were investigated, providing the most comprehensive, quantitative data on Micrurus venom composition to date, and more than tripling the number of Micrurus venom protein sequences previously available. The six venomes differ dramatically. All are dominated by 2-6 toxin classes that account for 91-99% of the toxin transcripts. The M. s. spixii venome is compositionally the simplest. In it, three-finger toxins (3FTxs) and phospholipases A₂ (PLA₂s) comprise >99% of the toxin transcripts, which include only four additional toxin families at levels ≥0.1%. Micrurus l. lemniscatus venom is the most complex, with at least 17 toxin families. However, in each venome, multiple structural subclasses of 3FTXs and PLA₂s are present. These almost certainly differ in pharmacology as well. All venoms also contain phospholipase B and vascular endothelial growth factors. Minor components (0.1-2.0%) are found in all venoms except that of M. s. spixii. Other toxin families are present in all six venoms at trace levels (<0.005%). Minor and trace venom components differ in each venom. Numerous novel toxin chemistries include 3FTxs with previously unknown 8- and 10-cysteine arrangements, resulting in new 3D structures and target specificities. 9-cysteine toxins raise the possibility of covalent, homodimeric 3FTxs or heterodimeric toxins with unknown pharmacologies. Probable muscarinic sequences may be reptile-specific homologs that promote hypotension via vascular mAChRs. The first complete sequences are presented for 3FTxs putatively responsible for liberating glutamate from rat brain synaptosomes. Micrurus C-type lectin-like proteins may have 6-9 cysteine residues and may be monomers, or homo- or heterodimers of unknown pharmacology. Novel KSPIs, 3× longer than any seen previously, appear to have arisen in three species by gene duplication and fusion. Four species have transcripts homologous to the nociceptive toxin, (MitTx) α-subunit, but all six species had homologs to the β-subunit. The first non-neurotoxic, non-catalytic elapid phospholipase A₂s are reported. All are probably myonecrotic. Phylogenetic analysis indicates that the six taxa diverged 15-35 million years ago and that they split from their last common ancestor with Old World elapines nearly 55 million years ago. Given their early diversification, many cryptic micrurine taxa are anticipated.

Identification and comparative analysis of the epidermal differentiation complex in snakes

The epidermis of snakes efficiently protects against dehydration and mechanical stress. However, only few proteins of the epidermal barrier to the environment have so far been identified in snakes. Here, we determined the organization of the Epidermal Differentiation Complex (EDC), a cluster of genes encoding protein constituents of cornified epidermal structures, in snakes and compared it to the EDCs of other squamates and non-squamate reptiles. The EDC of snakes displays shared synteny with that of the green anole lizard, including the presence of a cluster of corneous beta-protein (CBP)/beta-keratin genes. We found that a unique CBP comprising 4 putative beta-sheets and multiple cysteine-rich EDC proteins are conserved in all snakes and other squamates investigated. Comparative genomics of squamates suggests that the evolution of snakes was associated with a gene duplication generating two isoforms of the S100 fused-type protein, scaffoldin, the origin of distinct snake-specific EDC genes, and the loss of other genes that were present in the EDC of the last common ancestor of snakes and lizards. Taken together, our results provide new insights into the evolution of the skin in squamates and a basis for the characterization of the molecular composition of the epidermis in snakes.

Snake Genome Sequencing: Results and Future Prospects

Snake genome sequencing is in its infancy—very much behind the progress made in sequencing the genomes of humans, model organisms and pathogens relevant to biomedical research, and agricultural species. We provide here an overview of some of the snake genome projects in progress, and discuss the biological findings, with special emphasis on toxinology, from the small number of draft snake genomes already published. We discuss the future of snake genomics, pointing out that new sequencing technologies will help overcome the problem of repetitive sequences in assembling snake genomes. Genome sequences are also likely to be valuable in examining the clustering of toxin genes on the chromosomes, in designing recombinant antivenoms and in studying the epigenetic regulation of toxin gene expression.

Rapid Radiations and the Race to Redundancy: An Investigation of the Evolution of Australian Elapid Snake Venoms

Australia is the stronghold of the front-fanged venomous snake family Elapidae. The Australasian elapid snake radiation, which includes approximately 100 terrestrial species in Australia, as well as Melanesian species and all the world's sea snakes, is less than 12 million years old. The incredible phenotypic and ecological diversity of the clade is matched by considerable diversity in venom composition. The clade's evolutionary youth and dynamic evolution should make it of particular interest to toxinologists, however, the majority of species, which are small, typically inoffensive, and seldom encountered by non-herpetologists, have been almost completely neglected by researchers. The present study investigates the venom composition of 28 species proteomically, revealing several interesting trends in venom composition, and reports, for the first time in elapid snakes, the existence of an ontogenetic shift in the venom composition and activity of brown snakes (Pseudonaja sp.). Trends in venom composition are compared to the snakes' feeding ecology and the paper concludes with an extended discussion of the selection pressures shaping the evolution of snake venom.

Draft genome of the leopard gecko, Eublepharis macularius

Background

Geckos are among the most species-rich reptile groups and the sister clade to all other lizards and snakes. Geckos possess a suite of distinctive characteristics, including adhesive digits, nocturnal activity, hard, calcareous eggshells, and a lack of eyelids. However, one gecko clade, the Eublepharidae, appears to be the exception to most of these ‘rules’ and lacks adhesive toe pads, has eyelids, and lays eggs with soft, leathery eggshells. These differences make eublepharids an important component of any investigation into the underlying genomic innovations contributing to the distinctive phenotypes in ‘typical’ geckos.

Findings

We report high-depth genome sequencing, assembly, and annotation for a male leopard gecko, Eublepharis macularius (Eublepharidae). Illumina sequence data were generated from seven insert libraries (ranging from 170 to 20 kb), representing a raw sequencing depth of 136X from 303 Gb of data, reduced to 84X and 187 Gb after filtering. The assembled genome of 2.02 Gb was close to the 2.23 Gb estimated by k-mer analysis. Scaffold and contig N50 sizes of 664 and 20 kb, respectively, were comparable to the previously published Gekko japonicus genome. Repetitive elements accounted for 42 % of the genome. Gene annotation yielded 24,755 protein-coding genes, of which 93 % were functionally annotated. CEGMA and BUSCO assessment showed that our assembly captured 91 % (225 of 248) of the core eukaryotic genes, and 76 % of vertebrate universal single-copy orthologs.

Conclusions

Assembly of the leopard gecko genome provides a valuable resource for future comparative genomic studies of geckos and other squamate reptiles.

Electronic supplementary material

The online version of this article (doi:10.1186/s13742-016-0151-4) contains supplementary material, which is available to authorized users.

The Snake with the Scorpion's Sting: Novel Three-Finger Toxin Sodium Channel Activators from the Venom of the Long-Glanded Blue Coral Snake (Calliophis bivirgatus)

Millions of years of evolution have fine-tuned the ability of venom peptides to rapidly incapacitate both prey and potential predators. Toxicofera reptiles are characterized by serous-secreting mandibular or maxillary glands with heightened levels of protein expression. These glands are the core anatomical components of the toxicoferan venom system, which exists in myriad points along an evolutionary continuum. Neofunctionalisation of toxins is facilitated by positive selection at functional hotspots on the ancestral protein and venom proteins have undergone dynamic diversification in helodermatid and varanid lizards as well as advanced snakes. A spectacular point on the venom system continuum is the long-glanded blue coral snake (Calliophis bivirgatus), a specialist feeder that preys on fast moving, venomous snakes which have both a high likelihood of prey escape but also represent significant danger to the predator itself. The maxillary venom glands of C. bivirgatus extend one quarter of the snake's body length and nestle within the rib cavity. Despite the snake's notoriety its venom has remained largely unstudied. Here we show that the venom uniquely produces spastic paralysis, in contrast to the flaccid paralysis typically produced by neurotoxic snake venoms. The toxin responsible, which we have called calliotoxin (δ-elapitoxin-Cb1a), is a three-finger toxin (3FTx). Calliotoxin shifts the voltage-dependence of NaV1.4 activation to more hyperpolarised potentials, inhibits inactivation, and produces large ramp currents, consistent with its profound effects on contractile force in an isolated skeletal muscle preparation. Voltage-gated sodium channels (NaV) are a particularly attractive pharmacological target as they are involved in almost all physiological processes including action potential generation and conduction. Accordingly, venom peptides that interfere with NaV function provide a key defensive and predatory advantage to a range of invertebrate venomous species including cone snails, scorpions, spiders, and anemones. Enhanced activation or delayed inactivation of sodium channels by toxins is associated with the extremely rapid onset of tetanic/excitatory paralysis in envenomed prey animals. A strong selection pressure exists for the evolution of such toxins where there is a high chance of prey escape. However, despite their prevalence in other venomous species, toxins causing delay of sodium channel inhibition have never previously been described in vertebrate venoms. Here we show that NaV modulators, convergent with those of invertebrates, have evolved in the venom of the long-glanded coral snake. Calliotoxin represents a functionally novel class of 3FTx and a structurally novel class of NaV toxins that will provide significant insights into the pharmacology and physiology of NaV. The toxin represents a remarkable case of functional convergence between invertebrate and vertebrate venom systems in response to similar selection pressures. These results underscore the dynamic evolution of the Toxicofera reptile system and reinforces the value of using evolution as a roadmap for biodiscovery.

The Deep Origin and Recent Loss of Venom Toxin Genes in Rattlesnakes

The genetic origin of novel traits is a central but challenging puzzle in evolutionary biology. Among snakes, phospholipase A2 (PLA2)-related toxins have evolved in different lineages to function as potent neurotoxins, myotoxins, or hemotoxins. Here, we traced the genomic origin and evolution of PLA2 toxins by examining PLA2 gene number, organization, and expression in both neurotoxic and non-neurotoxic rattlesnakes. We found that even though most North American rattlesnakes do not produce neurotoxins, the genes of a specialized heterodimeric neurotoxin predate the origin of rattlesnakes and were present in their last common ancestor (∼22 mya). The neurotoxin genes were then deleted independently in the lineages leading to the Western Diamondback (Crotalus atrox) and Eastern Diamondback (C. adamanteus) rattlesnakes (∼6 mya), while a PLA2 myotoxin gene retained in C. atrox was deleted from the neurotoxic Mojave rattlesnake (C. scutulatus; ∼4 mya). The rapid evolution of PLA2 gene number appears to be due to transposon invasion that provided a template for non-allelic homologous recombination.

RNA sequence analyses of r-Moj-DM treated cells: TXNIP is required to induce apoptosis of SK-Mel-28

RNA sequencing of untreated and r-Moj-DM treated SK-Mel-28 cells was performed after 6 h, to begin unraveling the apoptotic pathway induced by r-Moj-DM. Bioinformatic analyses of RNA sequencing data yielded 40 genes that were differentially expressed. Nine genes were upregulated and 31 were downregulated. qRT-PCR was used to validate differential expression of 13 genes with known survival or apoptotic-inducing activities. Expression of BNiP3, IGFBP3, PTPSF, Prune 2, TGF-ß, and TXNIP were compared from cells treated with r-Moj-DN (a strong apoptotic inducer) or r-Moj-DA (a non-apoptotic inducer) for 1 h, 2 h, 4 h, and 6 h after treatment. Our results demonstrate that significant differences in expression are only detected after 4 h of treatment. In addition, expression of TXNIP (an apoptotic inducer) remains elevated at 4 h and 6 h only in r-Moj-DN treated cells. Based on the consistency of elevated TXNIP expression, we further studied TXNIP as a novel target of disintegrin activation. Confocal microscopy of anti-TXNIP stained SK-Mel-28 cells suggests nuclear localization of TXNIP after r-Moj-DM treatment. A stable TXNIP knockdown SK-Mel-28 cell line was produced to test TXNIP' role in the apoptotic induction by r-Moj-DM. High cell viability (74.3% ±9.1) was obtained after r-Moj-DM treatment of TXNIP knocked down SK-Mel-28 cells, compared to 34% ±0.187 for untransduced cells. These results suggest that TXNIP is required early in the apoptotic-inducing pathway resulting from r-Moj-DM binding to the αv integrin subunit.

Intragenome Diversity of Gene Families Encoding Toxin-like Proteins in Venomous Animals

The evolution of venoms is the story of how toxins arise and of the processes that generate and maintain their diversity. For animal venoms these processes include recruitment for expression in the venom gland, neofunctionalization, paralogous expansions, and functional divergence. The systematic study of these processes requires the reliable identification of the venom components involved in antagonistic interactions. High-throughput sequencing has the potential of uncovering the entire set of toxins in a given organism, yet the existence of non-venom toxin paralogs and the misleading effects of partial census of the molecular diversity of toxins make necessary to collect complementary evidence to distinguish true toxins from their non-venom paralogs. Here, we analyzed the whole genomes of two scorpions, one spider and one snake, aiming at the identification of the full repertoires of genes encoding toxin-like proteins. We classified the entire set of protein-coding genes into paralogous groups and monotypic genes, identified genes encoding toxin-like proteins based on known toxin families, and quantified their expression in both venom-glands and pooled tissues. Our results confirm that genes encoding toxin-like proteins are part of multigene families, and that these families arise by recruitment events from non-toxin genes followed by limited expansions of the toxin-like protein coding genes. We also show that failing to account for sequence similarity with non-toxin proteins has a considerable misleading effect that can be greatly reduced by comparative transcriptomics. Our study overall contributes to the understanding of the evolutionary dynamics of proteins involved in antagonistic interactions.

Trends in the Evolution of Snake Toxins Underscored by an Integrative Omics Approach to Profile the Venom of the Colubrid Phalotris mertensi

Only few studies on snake venoms were dedicated to deeply characterize the toxin secretion of animals from the Colubridae family, despite the fact that they represent the majority of snake diversity. As a consequence, some evolutionary trends observed in venom proteins that underpinned the evolutionary histories of snake toxins were based on data from a minor parcel of the clade. Here, we investigated the proteins of the totally unknown venom from Phalotris mertensi (Dipsadinae subfamily), in order to obtain a detailed profile of its toxins and to appreciate evolutionary tendencies occurring in colubrid venoms. By means of integrated omics and functional approaches, including RNAseq, Sanger sequencing, high-resolution proteomics, recombinant protein production, and enzymatic tests, we verified an active toxic secretion containing up to 21 types of proteins. A high content of Kunitz-type proteins and C-type lectins were observed, although several enzymatic components such as metalloproteinases and an L-amino acid oxidase were also present in the venom. Interestingly, an arguable venom component of other species was demonstrated as a true venom protein and named svLIPA (snake venom acid lipase). This finding indicates the importance of checking the actual protein occurrence across species before rejecting genes suggested to code for toxins, which are relevant for the discussion about the early evolution of reptile venoms. Moreover, trends in the evolution of some toxin classes, such as simplification of metalloproteinases and rearrangements of Kunitz and Wap domains, parallel similar phenomena observed in other venomous snake families and provide a broader picture of toxin evolution.

Colubrid Venom Composition: An -Omics Perspective

Snake venoms have been subjected to increasingly sensitive analyses for well over 100 years, but most research has been restricted to front-fanged snakes, which actually represent a relatively small proportion of extant species of advanced snakes. Because rear-fanged snakes are a diverse and distinct radiation of the advanced snakes, understanding venom composition among "colubrids" is critical to understanding the evolution of venom among snakes. Here we review the state of knowledge concerning rear-fanged snake venom composition, emphasizing those toxins for which protein or transcript sequences are available. We have also added new transcriptome-based data on venoms of three species of rear-fanged snakes. Based on this compilation, it is apparent that several components, including cysteine-rich secretory proteins (CRiSPs), C-type lectins (CTLs), CTLs-like proteins and snake venom metalloproteinases (SVMPs), are broadly distributed among "colubrid" venoms, while others, notably three-finger toxins (3FTxs), appear nearly restricted to the Colubridae (sensu stricto). Some putative new toxins, such as snake venom matrix metalloproteinases, are in fact present in several colubrid venoms, while others are only transcribed, at lower levels. This work provides insights into the evolution of these toxin classes, but because only a small number of species have been explored, generalizations are still rather limited. It is likely that new venom protein families await discovery, particularly among those species with highly specialized diets.

Insights into the Evolution of a Snake Venom Multi-Gene Family from the Genomic Organization of Echis ocellatus SVMP Genes

The molecular events underlying the evolution of the Snake Venom Metalloproteinase (SVMP) family from an A Disintegrin And Metalloproteinase (ADAM) ancestor remain poorly understood. Comparative genomics may provide decisive information to reconstruct the evolutionary history of this multi-locus toxin family. Here, we report the genomic organization of Echis ocellatus genes encoding SVMPs from the PII and PI classes. Comparisons between them and between these genes and the genomic structures of Anolis carolinensis ADAM28 and E. ocellatus PIII-SVMP EOC00089 suggest that insertions and deletions of intronic regions played key roles along the evolutionary pathway that shaped the current diversity within the multi-locus SVMP gene family. In particular, our data suggest that emergence of EOC00028-like PI-SVMP from an ancestral PII(e/d)-type SVMP involved splicing site mutations that abolished both the 3′ splice AG acceptor site of intron 12* and the 5′ splice GT donor site of intron 13*, and resulted in the intronization of exon 13* and the consequent destruction of the structural integrity of the PII-SVMP characteristic disintegrin domain.

From Mollusks to Medicine: A Venomics Approach for the Discovery and Characterization of Therapeutics from Terebridae Peptide Toxins

Animal venoms comprise a diversity of peptide toxins that manipulate molecular targets such as ion channels and receptors, making venom peptides attractive candidates for the development of therapeutics to benefit human health. However, identifying bioactive venom peptides remains a significant challenge. In this review we describe our particular venomics strategy for the discovery, characterization, and optimization of Terebridae venom peptides, teretoxins. Our strategy reflects the scientific path from mollusks to medicine in an integrative sequential approach with the following steps: (1) delimitation of venomous Terebridae lineages through taxonomic and phylogenetic analyses; (2) identification and classification of putative teretoxins through omics methodologies, including genomics, transcriptomics, and proteomics; (3) chemical and recombinant synthesis of promising peptide toxins; (4) structural characterization through experimental and computational methods; (5) determination of teretoxin bioactivity and molecular function through biological assays and computational modeling; (6) optimization of peptide toxin affinity and selectivity to molecular target; and (7) development of strategies for effective delivery of venom peptide therapeutics. While our research focuses on terebrids, the venomics approach outlined here can be applied to the discovery and characterization of peptide toxins from any venomous taxa.

Identification of hyaluronidase and phospholipase B in Lachesis muta rhombeata venom

Hyaluronidases contribute to local and systemic damages after envenoming, since they act as spreading factors cleaving the hyaluronan presents in the connective tissues of the victim, facilitating the diffusion of venom components. Although hyaluronidases are ubiquitous in snake venoms, they still have not been detected in transcriptomic analysis of the Lachesis venom gland and neither in the proteome of its venom performed previously. This work purified a hyaluronidase from Lachesis muta rhombeata venom whose molecular mass was estimated by SDS-PAGE to be 60 kDa. The hyaluronidase was more active at pH 6 and 37 °C when salt concentration was kept constant and more active in the presence of 0.15 M monovalent ions when the pH was kept at 6. Venom was fractionated by reversed-phase liquid chromatography (RPLC). Edman sequencing after RPLC failed to detect hyaluronidase, but identified a new serine proteinase isoform. The hyaluronidase was identified by mass spectrometry analysis of the protein bands in SDS-PAGE. Additionally, phospholipase B was identified for the first time in Lachesis genus venom. The discovery of new bioactive molecules might contribute to the design of novel drugs and biotechnology products as well as to development of more effective treatments against the envenoming.

Assessing the utility of the Oxford Nanopore MinION for snake venom gland cDNA sequencing

Portable DNA sequencers such as the Oxford Nanopore MinION device have the potential to be truly disruptive technologies, facilitating new approaches and analyses and, in some cases, taking sequencing out of the lab and into the field. However, the capabilities of these technologies are still being revealed. Here we show that single-molecule cDNA sequencing using the MinION accurately characterises venom toxin-encoding genes in the painted saw-scaled viper, Echis coloratus. We find the raw sequencing error rate to be around 12%, improved to 0–2% with hybrid error correction and 3% with de novo error correction. Our corrected data provides full coding sequences and 5′ and 3′ UTRs for 29 of 33 candidate venom toxins detected, far superior to Illumina data (13/40 complete) and Sanger-based ESTs (15/29). We suggest that, should the current pace of improvement continue, the MinION will become the default approach for cDNA sequencing in a variety of species.

Honeybee venom proteome profile of queens and winter bees as determined by a mass spectrometric approach

Venoms of invertebrates contain an enormous diversity of proteins, peptides, and other classes of substances. Insect venoms are characterized by a large interspecific variation resulting in extended lists of venom compounds. The venom composition of several hymenopterans also shows different intraspecific variation. For instance, venom from different honeybee castes, more specifically queens and workers, shows quantitative and qualitative variation, while the environment, like seasonal changes, also proves to be an important factor. The present study aimed at an in-depth analysis of the intraspecific variation in the honeybee venom proteome. In summer workers, the recent list of venom proteins resulted from merging combinatorial peptide ligand library sample pretreatment and targeted tandem mass spectrometry realized with a Fourier transform ion cyclotron resonance mass spectrometer (FT-ICR MS/MS). Now, the same technique was used to determine the venom proteome of queens and winter bees, enabling us to compare it with that of summer bees. In total, 34 putative venom toxins were found, of which two were never described in honeybee venoms before. Venom from winter workers did not contain toxins that were not present in queens or summer workers, while winter worker venom lacked the allergen Api m 12, also known as vitellogenin. Venom from queen bees, on the other hand, was lacking six of the 34 venom toxins compared to worker bees, while it contained two new venom toxins, in particularly serine proteinase stubble and antithrombin-III. Although people are hardly stung by honeybees during winter or by queen bees, these newly identified toxins should be taken into account in the characterization of a putative allergic response against Apis mellifera stings.

Ecological venomics: How genomics, transcriptomics and proteomics can shed new light on the ecology and evolution of venom

Animal venom is a complex cocktail of bioactive chemicals that traditionally drew interest mostly from biochemists and pharmacologists. However, in recent years the evolutionary and ecological importance of venom is realized as this trait has direct and strong influence on interactions between species. Moreover, venom content can be modulated by environmental factors. Like many other fields of biology, venom research has been revolutionized in recent years by the introduction of systems biology approaches, i.e., genomics, transcriptomics and proteomics. The employment of these methods in venom research is known as 'venomics'. In this review we describe the history and recent advancements of venomics and discuss how they are employed in studying venom in general and in particular in the context of evolutionary ecology. We also discuss the pitfalls and challenges of venomics and what the future may hold for this emerging scientific field.

Snake venoms: A brief treatise on etymology, origins of terminology, and definitions

The ancient perceptions of "venomous" and "poisonous snakes", as well as the Indo-European (IE) etymological origins of the term "venom" specifically associated with snakes are considered. Although several ancient cultures perceived snakes as symbols of fecundity and renewal, concurrent beliefs also associated venomous snakes with undesirable human characteristics or as portending non-propitious events. The respective IE roots of the terms "venom" and "poison", "wen" and "poi" refer to desire or the act of ingesting liquids. The origin of the term, "venom", is associated with polytheistic cults that emphasized attainment of desires sometimes assisted by "love potions", a term later interpolated with the word, "poison". Specific interpretation of the term, venom, has varied since its first probable use in the mid-Thirteenth Century. The definition of snake venom has long been contended, and interpretations have often reflected emphasis on the pharmacological or experimental toxicity of medically relevant snake venoms with less regard for the basic biological bases of these venoms, as well as those from snakes with no known medical significance. Several definitions of "snake venom" and their defining criteria are reviewed, and critical consideration is given to traditional criteria that might facilitate the future establishment of a biologically accurate definition.