From genome to "venome": molecular origin and evolution of the snake venom proteome inferred from phylogenetic analysis of toxin sequences and related body proteins

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

The identification of apparently conserved gene complements in the venom and salivary glands of a diverse set of reptiles led to the development of the Toxicofera hypothesis - the single, early evolution of the venom system in reptiles. However, this hypothesis is based largely on relatively small scale EST-based studies of only venom or salivary glands and toxic effects have been assigned to only some putative Toxicoferan toxins in some species. We set out to examine the distribution of these proposed venom toxin transcripts in order to investigate to what extent conservation of gene complements may reflect a bias in previous sampling efforts. Our quantitative transcriptomic analyses of venom and salivary glands and other body tissues in five species of reptile, together with the use of available RNA-Seq datasets for additional species, shows that the majority of genes used to support the establishment and expansion of the Toxicofera are in fact expressed in multiple body tissues and most likely represent general maintenance or "housekeeping" genes. The apparent conservation of gene complements across the Toxicofera therefore reflects an artefact of incomplete tissue sampling. We therefore conclude that venom has evolved multiple times in reptiles.

Venom-related transcripts from Bothrops jararaca tissues provide novel molecular insights into the production and evolution of snake venom

Attempts to reconstruct the evolutionary history of snake toxins in the context of their co-option to the venom gland rarely account for nonvenom snake genes that are paralogous to toxins, and which therefore represent important connectors to ancestral genes. In order to reevaluate this process, we conducted a comparative transcriptomic survey on body tissues from a venomous snake. A nonredundant set of 33,000 unigenes (assembled transcripts of reference genes) was independently assembled from six organs of the medically important viperid snake Bothrops jararaca, providing a reference list of 82 full-length toxins from the venom gland and specific products from other tissues, such as pancreatic digestive enzymes. Unigenes were then screened for nontoxin transcripts paralogous to toxins revealing 1) low level coexpression of approximately 20% of toxin genes (e.g., bradykinin-potentiating peptide, C-type lectin, snake venom metalloproteinase, snake venom nerve growth factor) in body tissues, 2) the identity of the closest paralogs to toxin genes in eight classes of toxins, 3) the location and level of paralog expression, indicating that, in general, co-expression occurs in a higher number of tissues and at lower levels than observed for toxin genes, and 4) strong evidence of a toxin gene reverting back to selective expression in a body tissue. In addition, our differential gene expression analyses identify specific cellular processes that make the venom gland a highly specialized secretory tissue. Our results demonstrate that the evolution and production of venom in snakes is a complex process that can only be understood in the context of comparative data from other snake tissues, including the identification of genes paralogous to venom toxins.

RNA-seq and high-definition mass spectrometry reveal the complex and divergent venoms of two rear-fanged colubrid snakes

Background

Largely because of their direct, negative impacts on human health, the venoms of front-fanged snakes of the families Viperidae and Elapidae have been extensively characterized proteomically, transcriptomically, and pharmacologically. However, relatively little is known about the molecular complexity and evolution of the venoms of rear-fanged colubrid snakes, which are, with a few notable exceptions, regarded as harmless to humans. Many of these snakes have venoms with major effects on their preferred prey, and their venoms are probably as critical to their survival as those of front-fanged elapids and viperids.

Results

We sequenced the venom-gland transcriptomes from a specimen of Hypsiglena (Desert Night Snake; family Colubridae, subfamily Dipsadinae) and of Boiga irregularis (Brown Treesnake; family Colubridae, subfamily Colubrinae) and verified the transcriptomic results proteomically by means of high-definition mass spectrometry. We identified nearly 3,000 nontoxin genes for each species. For B. irregularis, we found 108 putative toxin transcripts in 46 clusters with <1% nucleotide divergence, and for Hypsiglena we identified 79 toxin sequences that were grouped into 33 clusters. Comparisons of the venoms revealed divergent venom types, with Hypsiglena possessing a viper-like venom dominated by metalloproteinases, and B. irregularis having a more elapid-like venom, consisting primarily of three-finger toxins.

Conclusions

Despite the difficulty of procuring venom from rear-fanged species, we were able to complete all analyses from a single specimen of each species without pooling venom samples or glands, demonstrating the power of high-definition transcriptomic and proteomic approaches. We found a high level of divergence in the venom types of two colubrids. These two venoms reflected the hemorrhagic/neurotoxic venom dichotomy that broadly characterizes the difference in venom strategies between elapids and viperids.

Evolution of endothelin receptors in vertebrates

Endothelin receptors are G protein coupled receptors (GPCRs) of the β-group of rhodopsin receptors that bind to endothelin ligands, which are 21 amino acid long peptides derived from longer prepro-endothelin precursors. The most basal Ednr-like GPCR is found outside vertebrates in the cephalochordate amphioxus, but endothelin ligands are only present among vertebrates, including the lineages of jawless vertebrates (lampreys and hagfishes), cartilaginous vertebrates (sharks, rays, and chimaeras), and bony vertebrates (ray-finned fishes and lobe-finned vertebrates including tetrapods). A bona fide endothelin system is thus a vertebrate-specific innovation with important roles for regulating the cardiovascular system, renal and pulmonary processes, as well as for the development of the vertebrate-specific neural crest cell population and its derivatives. Expectedly, dysregulation of endothelin receptors and the endothelin system leads to a multitude of human diseases. Despite the importance of different types of endothelin receptors for vertebrate development and physiology, current knowledge on endothelin ligand-receptor interactions, on the expression of endothelin receptors and their ligands, and on the functional roles of the endothelin system for embryonic development and in adult vertebrates is very much biased towards amniote vertebrates. Recent analyses from a variety of vertebrate lineages, however, have shown that the endothelin system in lineages such as teleost fish and lampreys is more diverse and is divergent from the mammalian endothelin system. This diversity is mainly based on differential evolution of numerous endothelin system components among vertebrate lineages generated by two rounds of whole genome duplication (three in teleosts) during vertebrate evolution. Here we review current understanding of the evolutionary history of the endothelin receptor family in vertebrates supplemented with surveys on the endothelin receptor gene complement of newly available genome assemblies from phylogenetically informative taxa. Our assessment further highlights the diversity of the vertebrate endothelin system and calls for detailed functional and pharmacological analyses of the endothelin system beyond tetrapods.

NMR structure of bitistatin – a missing piece in the evolutionary pathway of snake venom disintegrins

Extant disintegrins, as found in the venoms of Viperidae and Crotalidae snakes (vipers and rattlesnakes, represent a family of polypeptides that block the function of β1 and β3 integrin receptors, both potently and with a high degree of selectivity. This toxin family owes its origin to the neofunctionalization of the extracellular region of an ADAM (a disintegrin and metalloprotease) molecule recruited into the snake venom gland proteome in the Jurassic. The evolutionary structural diversification of the disintegrin scaffold, from the ancestral long disintegrins to the more recently evolved medium-sized, dimeric and short disintegrins, involved the stepwise loss of pairs of class-specific disulfide linkages and the processing of the N-terminal region. NMR and crystal structures of medium-sized, dimeric and short disintegrins have been solved. However, the structure of a long disintegrin remained unknown. The present study reports the NMR solution structures of two disulfide bond conformers of the long disintegrin bitistatin from the African puff adder Bitis arietans. The findings provide insight into how a structural domain of the extracellular region of an ADAM molecule, recruited into and selectively expressed in the snake venom gland proteome as a PIII metalloprotease in the Jurassic, has subsequently been tranformed into a family of integrin receptor antagonists.

Crystallization and preliminary X-ray diffraction studies of an L-amino-acid oxidase from Lachesis muta venom

Snake-venom proteins form multi-component defence systems by the recruitment and rapid evolution of nonvenomous proteins and hence serve as model systems to understand the structural modifications that result in toxicity. L-Amino-acid oxidases (LAAOs) are encountered in a number of snake venoms and have been implicated in the inhibition of platelet aggregation, cytotoxicity, haemolysis, apoptosis and haemorrhage. An L-amino-acid oxidase from Lachesis muta venom has been purified and crystallized. The crystals belonged to space group P2₁, with unit-cell parameters a=66.05, b=79.41, c=100.52 Å, β=96.55°. The asymmetric unit contained two molecules and the structure has been determined and partially refined at 3.0 Å resolution.

Identification and phylogeny of Arabian snakes: Comparison of venom chromatographic profiles versus 16S rRNA gene sequences

Identification of snake species is important for various reasons including the emergency treatment of snake bite victims. We present a simple method for identification of six snake species using the gel filtration chromatographic profiles of their venoms. The venoms of Echis coloratus, Echis pyramidum, Cerastes gasperettii, Bitis arietans, Naja arabica, and Walterinnesia aegyptia were milked, lyophilized, diluted and centrifuged to separate the mucus from the venom. The clear supernatants were filtered and chromatographed on fast protein liquid chromatography (FPLC). We obtained the 16S rRNA gene sequences of the above species and performed phylogenetic analysis using the neighbor-joining method. The chromatograms of venoms from different snake species showed peculiar patterns based on the number and location of peaks. The dendrograms generated from similarity matrix based on the presence/absence of particular chromatographic peaks clearly differentiated Elapids from Viperids. Molecular cladistics using 16S rRNA gene sequences resulted in jumping clades while separating the members of these two families. These findings suggest that chromatographic profiles of snake venoms may provide a simple and reproducible chemical fingerprinting method for quick identification of snake species. However, the validation of this methodology requires further studies on large number of specimens from within and across species.

Venom gland transcriptomics for identifying, cataloging, and characterizing venom proteins in snakes

Snake venoms are cocktails of protein toxins that play important roles in capture and digestion of prey. Significant qualitative and quantitative variation in snake venom composition has been observed among and within species. Understanding these variations in protein components is instrumental in interpreting clinical symptoms during human envenomation and in searching for novel venom proteins with potential therapeutic applications. In the last decade, transcriptomic analyses of venom glands have helped in understanding the composition of various snake venoms in great detail. Here we review transcriptomic analysis as a powerful tool for understanding venom profile, variation and evolution.

Expression of venom gene homologs in diverse python tissues suggests a new model for the evolution of snake venom

Snake venom gene evolution has been studied intensively over the past several decades, yet most previous studies have lacked the context of complete snake genomes and the full context of gene expression across diverse snake tissues. We took a novel approach to studying snake venom evolution by leveraging the complete genome of the Burmese python, including information from tissue-specific patterns of gene expression. We identified the orthologs of snake venom genes in the python genome, and conducted detailed analysis of gene expression of these venom homologs to identify patterns that differ between snake venom gene families and all other genes. We found that venom gene homologs in the python are expressed in many different tissues outside of oral glands, which illustrates the pitfalls of using transcriptomic data alone to define "venom toxins." We hypothesize that the python may represent an ancestral state prior to major venom development, which is supported by our finding that the expansion of venom gene families is largely restricted to highly venomous caenophidian snakes. Therefore, the python provides insight into biases in which genes were recruited for snake venom systems. Python venom homologs are generally expressed at lower levels, have higher variance among tissues, and are expressed in fewer organs compared with all other python genes. We propose a model for the evolution of snake venoms in which venom genes are recruited preferentially from genes with particular expression profile characteristics, which facilitate a nearly neutral transition toward specialized venom system expression.

Extreme venom variation in Middle Eastern vipers: a proteomics comparison of Eristicophis macmahonii, Pseudocerastes fieldi and Pseudocerastes persicus

Venoms of the viperid sister genera Eristicophis and Pseudocerastes are poorly studied despite their anecdotal reputation for producing severe or even lethal envenomations. This is due in part to the remote and politically unstable regions that they occupy. All species contained are sit and wait ambush feeders. Thus, this study examined their venoms through proteomics techniques in order to establish if this feeding ecology, and putatively low levels of gene flow, have resulted in significant variations in venom profile. The techniques indeed revealed extreme venom variation. This has immediate implications as only one antivenom is made (using the venom of Pseudocerastes persicus) yet the proteomic variation suggests that it would be of only limited use for the other species, even the sister species Pseudocerastes fieldi. The high degree of variation however also points toward these species being rich resources for novel compounds which may have use as lead molecules in drug design and development.

BIOLOGICAL SIGNIFICANCE

These results show extreme venom variation between these closely related snakes. These results have direct implications for the treatment of the envenomed patient.

Gene structure, regulatory control, and evolution of black widow venom latrotoxins

Black widow venom contains α-latrotoxin, infamous for causing intense pain. Combining 33 kb of Latrodectus hesperus genomic DNA with RNA-Seq, we characterized the α-latrotoxin gene and discovered a paralog, 4.5 kb downstream. Both paralogs exhibit venom gland specific transcription, and may be regulated post-transcriptionally via musashi-like proteins. A 4 kb intron interrupts the α-latrotoxin coding sequence, while a 10 kb intron in the 3' UTR of the paralog may cause non-sense-mediated decay. Phylogenetic analysis confirms these divergent latrotoxins diversified through recent tandem gene duplications. Thus, latrotoxin genes have more complex structures, regulatory controls, and sequence diversity than previously proposed.

Proteomic analyses of courtship pheromones in the redback salamander, Plethodon cinereus

The evolutionary success of plethodontid salamanders for ~100 MY is due partly to the use of courtship pheromones that regulate female receptivity. In ~90 % of plethodontid species, males deliver pheromones by "scratching" a female's dorsum, where pheromones diffuse transdermally into the bloodstream. However, in a single clade, representing ~10 % of Plethodon spp., males apply pheromones to the female's nares for olfactory delivery. Molecular studies have identified three major pheromone families: Plethodontid Receptivity Factor (PRF), Plethodontid Modulating Factor (PMF), and Sodefrin Precursor-like Factor (SPF). SPF and PMF genes are relatively ancient and found in all plethodontid species; however, PRF is found exclusively in the genus Plethodon - which includes species with transdermal, olfactory, and intermediate delivery behaviors. While previous proteomic analyses suggested PRF and PMF are dominant in slapping species and SPF is dominant in non-Plethodon scratching species, it was unclear how protein expression of different pheromone components may vary across delivery modes within Plethodon. Therefore, the aim of this study was to proteomically characterize the pheromones of a key scratching species in this evolutionary transition, Plethodon cinereus. Using mass spectrometry-based techniques, our data support the functional replacement of SPF by PRF in Plethodon spp. and an increase in PMF gene duplication events in both lineage-dependent and delivery-dependent manners. Novel glycosylation was observed on P. cinereus PRFs, which may modulate the metabolism and/or mechanism of action for PRF in scratching species. Cumulatively, these molecular data suggest that the replacement of pheromone components (e.g., SPF by PRF) preceded the evolutionary transition of the functional complex from transdermal to olfactory delivery.

Restriction and recruitment-gene duplication and the origin and evolution of snake venom toxins

Snake venom has been hypothesized to have originated and diversified through a process that involves duplication of genes encoding body proteins with subsequent recruitment of the copy to the venom gland, where natural selection acts to develop or increase toxicity. However, gene duplication is known to be a rare event in vertebrate genomes, and the recruitment of duplicated genes to a novel expression domain (neofunctionalization) is an even rarer process that requires the evolution of novel combinations of transcription factor binding sites in upstream regulatory regions. Therefore, although this hypothesis concerning the evolution of snake venom is very unlikely and should be regarded with caution, it is nonetheless often assumed to be established fact, hindering research into the true origins of snake venom toxins. To critically evaluate this hypothesis, we have generated transcriptomic data for body tissues and salivary and venom glands from five species of venomous and nonvenomous reptiles. Our comparative transcriptomic analysis of these data reveals that snake venom does not evolve through the hypothesized process of duplication and recruitment of genes encoding body proteins. Indeed, our results show that many proposed venom toxins are in fact expressed in a wide variety of body tissues, including the salivary gland of nonvenomous reptiles and that these genes have therefore been restricted to the venom gland following duplication, not recruited. Thus, snake venom evolves through the duplication and subfunctionalization of genes encoding existing salivary proteins. These results highlight the danger of the elegant and intuitive "just-so story" in evolutionary biology.

Scorpions: a presentation

Scorpions, at least the species of the family Buthidæ whose venoms are better known, appear as animals that have evolved very little over time. The composition of their venoms is relatively simple as most toxins have a common structural motif that is found in other venoms from primitive species. Moreover, all the scorpion venom toxins principally act on membrane ionic channels of excitable cells. The results of recent works lead to the conclusion that in scorpions there is a close relationship between venomous function and innate immune function both remarkably efficient.

Molecular cloning of a hyaluronidase from Bothrops pauloensis venom gland

Background

Hyaluronate is one of the major components of extracellular matrix from vertebrates whose breakdown is catalyzed by the enzyme hyaluronidase. These enzymes are widely described in snake venoms, in which they facilitate the spreading of the main toxins in the victim’s body during the envenoming. Snake venoms also present some variants (hyaluronidases-like substances) that are probably originated by alternative splicing, even though their relevance in envenomation is still under investigation. Hyaluronidases-like proteins have not yet been purified from any snake venom, but the cDNA that encodes these toxins was already identified in snake venom glands by transcriptomic analysis. Herein, we report the cloning and in silico analysis of the first hyaluronidase-like proteins from a Brazilian snake venom.

Methods

The cDNA sequence of hyaluronidase was cloned from the transcriptome of Bothrops pauloensis venom glands. This sequence was submitted to multiple alignment with other related sequences by ClustalW. A phylogenetic analysis was performed using MEGA 4 software by the neighbor joining (NJ) method.

Results

The cDNA from Bothrops pauloensis venom gland that corresponds to hyaluronidase comprises 1175 bp and codifies a protein containing 194 amino acid residues. The sequence, denominated BpHyase, was identified as hyaluronidase-like since it shows high sequence identities (above 83%) with other described snake venom hyaluronidase-like sequences. Hyaluronidases-like proteins are thought to be products of alternative splicing implicated in deletions of central amino acids, including the catalytic residues. Structure-based sequence alignment of BpHyase to human hyaluronidase hHyal-1 demonstrates a loss of some key secondary structures. The phylogenetic analysis indicates an independent evolution of BpHyal when compared to other hyaluronidases. However, these toxins might share a common ancestor, thus suggesting a broad hyaluronidase-like distribution among venomous snakes.

Conclusions

This work is the first report of a cDNA sequence of hyaluronidase from Brazilian snake venoms. Moreover, the in silico analysis of its deduced amino acid sequence opens new perspectives about the biological function of hyaluronidases-like proteins and may direct further studies comprising their isolation and/or recombinant production, as well as their structural and functional characterization.

Medically important differences in snake venom composition are dictated by distinct postgenomic mechanisms

Variation in venom composition is a ubiquitous phenomenon in snakes and occurs both interspecifically and intraspecifically. Venom variation can have severe outcomes for snakebite victims by rendering the specific antibodies found in antivenoms ineffective against heterologous toxins found in different venoms. The rapid evolutionary expansion of different toxin-encoding gene families in different snake lineages is widely perceived as the main cause of venom variation. However, this view is simplistic and disregards the understudied influence that processes acting on gene transcription and translation may have on the production of the venom proteome. Here, we assess the venom composition of six related viperid snakes and compare interspecific changes in the number of toxin genes, their transcription in the venom gland, and their translation into proteins secreted in venom. Our results reveal that multiple levels of regulation are responsible for generating variation in venom composition between related snake species. We demonstrate that differential levels of toxin transcription, translation, and their posttranslational modification have a substantial impact upon the resulting venom protein mixture. Notably, these processes act to varying extents on different toxin paralogs found in different snakes and are therefore likely to be as important as ancestral gene duplication events for generating compositionally distinct venom proteomes. Our results suggest that these processes may also contribute to altering the toxicity of snake venoms, and we demonstrate how this variability can undermine the treatment of a neglected tropical disease, snakebite.

Structural insights into the evolution of a sexy protein: novel topology and restricted backbone flexibility in a hypervariable pheromone from the red-legged salamander, Plethodon shermani

In response to pervasive sexual selection, protein sex pheromones often display rapid mutation and accelerated evolution of corresponding gene sequences. For proteins, the general dogma is that structure is maintained even as sequence or function may rapidly change. This phenomenon is well exemplified by the three-finger protein (TFP) superfamily: a diverse class of vertebrate proteins co-opted for many biological functions - such as components of snake venoms, regulators of the complement system, and coordinators of amphibian limb regeneration. All of the >200 structurally characterized TFPs adopt the namesake "three-finger" topology. In male red-legged salamanders, the TFP pheromone Plethodontid Modulating Factor (PMF) is a hypervariable protein such that, through extensive gene duplication and pervasive sexual selection, individual male salamanders express more than 30 unique isoforms. However, it remained unclear how this accelerated evolution affected the protein structure of PMF. Using LC/MS-MS and multidimensional NMR, we report the 3D structure of the most abundant PMF isoform, PMF-G. The high resolution structural ensemble revealed a highly modified TFP structure, including a unique disulfide bonding pattern and loss of secondary structure, that define a novel protein topology with greater backbone flexibility in the third peptide finger. Sequence comparison, models of molecular evolution, and homology modeling together support that this flexible third finger is the most rapidly evolving segment of PMF. Combined with PMF sequence hypervariability, this structural flexibility may enhance the plasticity of PMF as a chemical signal by permitting potentially thousands of structural conformers. We propose that the flexible third finger plays a critical role in PMF:receptor interactions. As female receptors co-evolve, this flexibility may allow PMF to still bind its receptor(s) without the immediate need for complementary mutations. Consequently, this unique adaptation may establish new paradigms for how receptor:ligand pairs co-evolve, in particular with respect to sexual conflict.

Clawing through evolution: toxin diversification and convergence in the ancient lineage Chilopoda (centipedes)

Despite the staggering diversity of venomous animals, there seems to be remarkable convergence in regard to the types of proteins used as toxin scaffolds. However, our understanding of this fascinating area of evolution has been hampered by the narrow taxonomical range studied, with entire groups of venomous animals remaining almost completely unstudied. One such group is centipedes, class Chilopoda, which emerged about 440 Ma and may represent the oldest terrestrial venomous lineage next to scorpions. Here, we provide the first comprehensive insight into the chilopod "venome" and its evolution, which has revealed novel and convergent toxin recruitments as well as entirely new toxin families among both high- and low molecular weight venom components. The ancient evolutionary history of centipedes is also apparent from the differences between the Scolopendromorpha and Scutigeromorpha venoms, which diverged over 430 Ma, and appear to employ substantially different venom strategies. The presence of a wide range of novel proteins and peptides in centipede venoms highlights these animals as a rich source of novel bioactive molecules. Understanding the evolutionary processes behind these ancient venom systems will not only broaden our understanding of which traits make proteins and peptides amenable to neofunctionalization but it may also aid in directing bioprospecting efforts.

Tracing monotreme venom evolution in the genomics era

The monotremes (platypuses and echidnas) represent one of only four extant venomous mammalian lineages. Until recently, monotreme venom was poorly understood. However, the availability of the platypus genome and increasingly sophisticated genomic tools has allowed us to characterize platypus toxins, and provides a means of reconstructing the evolutionary history of monotreme venom. Here we review the physiology of platypus and echidna crural (venom) systems as well as pharmacological and genomic studies of monotreme toxins. Further, we synthesize current ideas about the evolution of the venom system, which in the platypus is likely to have been retained from a venomous ancestor, whilst being lost in the echidnas. We also outline several research directions and outstanding questions that would be productive to address in future research. An improved characterization of mammalian venoms will not only yield new toxins with potential therapeutic uses, but will also aid in our understanding of the way that this unusual trait evolves.

Next-generation snake venomics: protein-locus resolution through venom proteome decomplexation

Venom research has been continuously enhanced by technological advances. High-throughput technologies are changing the classical paradigm of hypothesis-driven research to technology-driven approaches. However, the thesis advocated in this paper is that full proteome coverage at locus-specific resolution requires integrating the best of both worlds into a protocol that includes decomplexation of the venom proteome prior to liquid chromatography-tandem mass spectrometry matching against a species-specific transcriptome. This approach offers the possibility of proof-checking the species-specific contig database using proteomics data. Immunoaffinity chromatography constitutes the basis of an antivenomics workflow designed to quantify the extent of cross-reactivity of antivenoms against homologous and heterologous venom toxins. In the author's view, snake venomics and antivenomics form part of a biology-driven conceptual framework to unveil the genesis and natural history of venoms, and their within- and between-species toxicological and immunological divergences and similarities. Understanding evolutionary trends across venoms represents the Rosetta Stone for generating broad-ranging polyspecific antivenoms.

Diversification of a single ancestral gene into a successful toxin superfamily in highly venomous Australian funnel-web spiders

BACKGROUND

Spiders have evolved pharmacologically complex venoms that serve to rapidly subdue prey and deter predators. The major toxic factors in most spider venoms are small, disulfide-rich peptides. While there is abundant evidence that snake venoms evolved by recruitment of genes encoding normal body proteins followed by extensive gene duplication accompanied by explosive structural and functional diversification, the evolutionary trajectory of spider-venom peptides is less clear.

RESULTS

Here we present evidence of a spider-toxin superfamily encoding a high degree of sequence and functional diversity that has evolved via accelerated duplication and diversification of a single ancestral gene. The peptides within this toxin superfamily are translated as prepropeptides that are posttranslationally processed to yield the mature toxin. The N-terminal signal sequence, as well as the protease recognition site at the junction of the propeptide and mature toxin are conserved, whereas the remainder of the propeptide and mature toxin sequences are variable. All toxin transcripts within this superfamily exhibit a striking cysteine codon bias. We show that different pharmacological classes of toxins within this peptide superfamily evolved under different evolutionary selection pressures.

CONCLUSIONS

Overall, this study reinforces the hypothesis that spiders use a combinatorial peptide library strategy to evolve a complex cocktail of peptide toxins that target neuronal receptors and ion channels in prey and predators. We show that the ω-hexatoxins that target insect voltage-gated calcium channels evolved under the influence of positive Darwinian selection in an episodic fashion, whereas the κ-hexatoxins that target insect calcium-activated potassium channels appear to be under negative selection. A majority of the diversifying sites in the ω-hexatoxins are concentrated on the molecular surface of the toxins, thereby facilitating neofunctionalisation leading to new toxin pharmacology.

Molecular models of the Mojave rattlesnake (Crotalus scutulatus scutulatus) venom metalloproteinases reveal a structural basis for differences in hemorrhagic activities

Rattlesnake venom can differ in composition and in metalloproteinase-associated activities. The molecular basis for this intra-species variation in Crotalus scutulatus scutulatus (Mojave rattlesnake) remains an enigma. To understand the molecular basis for intra-species variation of metalloproteinase-associated activities, we modeled the three-dimensional structures of four metalloproteinases based on the amino acid sequence of four variations of the proteinase domain of the C. s. scutulatus metalloproteinase gene (GP1, GP2, GP3, and GP4). For comparative purposes, we modeled the atrolysin metalloproteinases of C. atrox as well. All molecular models shared the same topology. While the atrolysin metalloproteinase molecular models contained highly conserved substrate binding sites, the Mojave rattlesnake metalloproteinases showed higher structural divergence when superimposed onto each other. The highest structural divergence among the four C. s. scutulatus molecular models was located at the northern cleft wall and the S'1-pocket of the substrate binding site, molecular regions that modulate substrate selectivity. Molecular dynamics and field potential maps for each C. s. scutulatus metalloproteinase model demonstrated that the non-hemorrhagic metalloproteinases (GP2 and GP3) contain highly basic molecular and field potential surfaces while the hemorrhagic metalloproteinases GP1 and atrolysin C showed extensive acidic field potential maps and shallow but less dynamic active site pockets. Hence, differences in the spatial arrangement of the northern cleft wall, the S'1-pocket, and the physico-chemical environment surrounding the catalytic site contribute to differences in metalloproteinase activities in the Mojave rattlesnake. Our results provide a structural basis for variation of metalloproteinase-associated activities in the rattlesnake venom of the Mojave rattlesnake.

State of the art in the studies on crotamine, a cell penetrating peptide from South American rattlesnake

Animal venoms comprise a naturally selected cocktail of bioactive peptides/proteins and other molecules, each of which playing a defined role thanks to the highly specific interactions with diverse molecular targets found in the prey. Research focused on isolation, structural, and functional characterizations of novel natural biologics (bioactive peptides/proteins from natural sources) has a long way to go through from the basic science to clinical applications. Herein, we overview the structural and functional characteristics of the myoneurotoxin crotamine, firstly isolated from the South American rattlesnake venom. Crotamine is the first venom peptide classified as a natural cell penetrating and antimicrobial peptide (CPP and AMP) with a more pronounced antifungal activity. In contrast to other known natural CPPs and AMPs, crotamine demonstrates a wide spectrum of biological activities with potential biotechnological and therapeutic values. More recent studies have demonstrated the selective in vitro anticancer activity of crotamine. In vivo, using a murine melanoma model, it was shown that crotamine delays tumor implantation, inhibits tumor cells proliferation, and also increases the survival of mice engrafted with subcutaneous melanoma. The structural and functional properties and also the possible biotechnological applications of minimized molecules derived from crotamine are also discussed.

Experimental conversion of a defensin into a neurotoxin: implications for origin of toxic function

Scorpion K(+) channel toxins and insect defensins share a conserved three-dimensional structure and related biological activities (defense against competitors or invasive microbes by disrupting their membrane functions), which provides an ideal system to study how functional evolution occurs in a conserved structural scaffold. Using an experimental approach, we show that the deletion of a small loop of a parasitoid venom defensin possessing the "scorpion toxin signature" (STS) can remove steric hindrance of peptide-channel interactions and result in a neurotoxin selectively inhibiting K(+) channels with high affinities. This insect defensin-derived toxin adopts a hallmark scorpion toxin fold with a common cysteine-stabilized α-helical and β-sheet motif, as determined by nuclear magnetic resonance analysis. Mutations of two key residues located in STS completely diminish or significantly decrease the affinity of the toxin on the channels, demonstrating that this toxin binds to K(+) channels in the same manner as scorpion toxins. Taken together, these results provide new structural and functional evidence supporting the predictability of toxin evolution. The experimental strategy is the first employed to establish an evolutionary relationship of two distantly related protein families.

Venom down under: dynamic evolution of Australian elapid snake toxins

Despite the unparalleled diversity of venomous snakes in Australia, research has concentrated on a handful of medically significant species and even of these very few toxins have been fully sequenced. In this study, venom gland transcriptomes were sequenced from eleven species of small Australian elapid snakes, from eleven genera, spanning a broad phylogenetic range. The particularly large number of sequences obtained for three-finger toxin (3FTx) peptides allowed for robust reconstructions of their dynamic molecular evolutionary histories. We demonstrated that each species preferentially favoured different types of α-neurotoxic 3FTx, probably as a result of differing feeding ecologies. The three forms of α-neurotoxin [Type I (also known as (aka): short-chain), Type II (aka: long-chain) and Type III] not only adopted differential rates of evolution, but have also conserved a diversity of residues, presumably to potentiate prey-specific toxicity. Despite these differences, the different α-neurotoxin types were shown to accumulate mutations in similar regions of the protein, largely in the loops and structurally unimportant regions, highlighting the significant role of focal mutagenesis. We theorize that this phenomenon not only affects toxin potency or specificity, but also generates necessary variation for preventing/delaying prey animals from acquiring venom-resistance. This study also recovered the first full-length sequences for multimeric phospholipase A₂ (PLA₂) ‘taipoxin/paradoxin’ subunits from non-Oxyuranus species, confirming the early recruitment of this extremely potent neurotoxin complex to the venom arsenal of Australian elapid snakes. We also recovered the first natriuretic peptides from an elapid that lack the derived C-terminal tail and resemble the plesiotypic form (ancestral character state) found in viper venoms. This provides supporting evidence for a single early recruitment of natriuretic peptides into snake venoms. Novel forms of kunitz and waprin peptides were recovered, including dual domain kunitz-kunitz precursors and the first kunitz-waprin hybrid precursors from elapid snakes. The novel sequences recovered in this study reveal that the huge diversity of unstudied venomous Australian snakes are of considerable interest not only for the investigation of venom and whole organism evolution but also represent an untapped bioresource in the search for novel compounds for use in drug design and development.

The king cobra genome reveals dynamic gene evolution and adaptation in the snake venom system

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.

Molecular evolution of vertebrate neurotrophins: co-option of the highly conserved nerve growth factor gene into the advanced snake venom arsenalf

Neurotrophins are a diverse class of structurally related proteins, essential for neuronal development, survival, plasticity and regeneration. They are characterized by major family members, such as the nerve growth factors (NGF), brain-derived neurotrophic factors (BDNF) and neurotrophin-3 (NT-3), which have been demonstrated here to lack coding sequence variations and follow the regime of negative selection, highlighting their extremely important conserved role in vertebrate homeostasis. However, in stark contrast, venom NGF secreted as part of the chemical arsenal of the venomous advanced snake family Elapidae (and to a lesser extent Viperidae) have characteristics consistent with the typical accelerated molecular evolution of venom components. This includes a rapid rate of diversification under the significant influence of positive-selection, with the majority of positively-selected sites found in the secreted β-polypeptide chain (74%) and on the molecular surface of the protein (92%), while the core structural and functional residues remain highly constrained. Such focal mutagenesis generates active residues on the toxin molecular surface, which are capable of interacting with novel biological targets in prey to induce a myriad of pharmacological effects. We propose that caenophidian NGFs could participate in prey-envenoming by causing a massive release of chemical mediators from mast cells to mount inflammatory reactions and increase vascular permeability, thereby aiding the spread of other toxins and/or by acting as proapoptotic factors. Despite their presence in reptilian venom having been known for over 60 years, this is the first evidence that venom-secreted NGF follows the molecular evolutionary pattern of other venom components, and thus likely participates in prey-envenomation.

Three-fingered RAVERs: Rapid Accumulation of Variations in Exposed Residues of snake venom toxins

Three-finger toxins (3FTx) represent one of the most abundantly secreted and potently toxic components of colubrid (Colubridae), elapid (Elapidae) and psammophid (Psammophiinae subfamily of the Lamprophidae) snake venom arsenal. Despite their conserved structural similarity, they perform a diversity of biological functions. Although they are theorised to undergo adaptive evolution, the underlying diversification mechanisms remain elusive. Here, we report the molecular evolution of different 3FTx functional forms and show that positively selected point mutations have driven the rapid evolution and diversification of 3FTx. These diversification events not only correlate with the evolution of advanced venom delivery systems (VDS) in Caenophidia, but in particular the explosive diversification of the clade subsequent to the evolution of a high pressure, hollow-fanged VDS in elapids, highlighting the significant role of these toxins in the evolution of advanced snakes. We show that Type I, II and III α-neurotoxins have evolved with extreme rapidity under the influence of positive selection. We also show that novel Oxyuranus/Pseudonaja Type II forms lacking the apotypic loop-2 stabilising cysteine doublet characteristic of Type II forms are not phylogenetically basal in relation to other Type IIs as previously thought, but are the result of secondary loss of these apotypic cysteines on at least three separate occasions. Not all 3FTxs have evolved rapidly: κ-neurotoxins, which form non-covalently associated heterodimers, have experienced a relatively weaker influence of diversifying selection; while cytotoxic 3FTx, with their functional sites, dispersed over 40% of the molecular surface, have been extremely constrained by negative selection. We show that the a previous theory of 3FTx molecular evolution (termed ASSET) is evolutionarily implausible and cannot account for the considerable variation observed in very short segments of 3FTx. Instead, we propose a theory of Rapid Accumulation of Variations in Exposed Residues (RAVER) to illustrate the significance of point mutations, guided by focal mutagenesis and positive selection in the evolution and diversification of 3FTx.

Atractaspis aterrima toxins: the first insight into the molecular evolution of venom in side-stabbers

Although snake venoms have been the subject of intense research, primarily because of their tremendous potential as a bioresource for design and development of therapeutic compounds, some specific groups of snakes, such as the genus Atractaspis, have been completely neglected. To date only limited number of toxins, such as sarafotoxins have been well characterized from this lineage. In order to investigate the molecular diversity of venom from Atractaspis aterrima—the slender burrowing asp, we utilized a high-throughput transcriptomic approach completed with an original bioinformatics analysis pipeline. Surprisingly, we found that Sarafotoxins do not constitute the major ingredient of the transcriptomic cocktail; rather a large number of previously well-characterized snake venom-components were identified. Notably, we recovered a large diversity of three-finger toxins (3FTxs), which were found to have evolved under the significant influence of positive selection. From the normalized and non-normalized transcriptome libraries, we were able to evaluate the relative abundance of the different toxin groups, uncover rare transcripts, and gain new insight into the transcriptomic machinery. In addition to previously characterized toxin families, we were able to detect numerous highly-transcribed compounds that possess all the key features of venom-components and may constitute new classes of toxins.

Short toxin-like proteins attack the defense line of innate immunity

ClanTox (classifier of animal toxins) was developed for identifying toxin-like candidates from complete proteomes. Searching mammalian proteomes for short toxin-like proteins (coined TOLIPs) revealed a number of overlooked secreted short proteins with an abundance of cysteines throughout their sequences. We applied bioinformatics and data-mining methods to infer the function of several top predicted candidates. We focused on cysteine-rich peptides that adopt the fold of the three-finger proteins (TFPs). We identified a cluster of duplicated genes that share a structural similarity with elapid neurotoxins, such as α-bungarotoxin. In the murine proteome, there are about 60 such proteins that belong to the Ly6/uPAR family. These proteins are secreted or anchored to the cell membrane. Ly6/uPAR proteins are associated with a rich repertoire of functions, including binding to receptors and adhesion. Ly6/uPAR proteins modulate cell signaling in the context of brain functions and cells of the innate immune system. We postulate that TOLIPs, as modulators of cell signaling, may be associated with pathologies and cellular imbalance. We show that proteins of the Ly6/uPAR family are associated with cancer diagnosis and malfunction of the immune system.

Squeezers and leaf-cutters: differential diversification and degeneration of the venom system in toxicoferan reptiles

Although it has been established that all toxicoferan squamates share a common venomous ancestor, it has remained unclear whether the maxillary and mandibular venom glands are evolving on separate gene expression trajectories or if they remain under shared genetic control. We show that identical transcripts are simultaneously expressed not only in the mandibular and maxillary glands, but also in the enigmatic snake rictal gland. Toxin molecular frameworks recovered in this study were three-finger toxin (3FTx), CRiSP, crotamine (beta-defensin), cobra venom factor, cystatin, epididymal secretory protein, kunitz, L-amino acid oxidase, lectin, renin aspartate protease, veficolin, and vespryn. We also discovered a novel low-molecular weight disulfide bridged peptide class in pythonid snake glands. In the iguanian lizards, the most highly expressed are potentially antimicrobial in nature (crotamine (beta-defensin) and cystatin), with crotamine (beta-defensin) also the most diverse. However, a number of proteins characterized from anguimorph lizards and caenophidian snakes with hemotoxic or neurotoxic activities were recruited in the common toxicoferan ancestor and remain expressed, albeit in low levels, even in the iguanian lizards. In contrast, the henophidian snakes express 3FTx and lectin toxins as the dominant transcripts. Even in the constricting pythonid and boid snakes, where the glands are predominantly mucous-secreting, low-levels of toxin transcripts can be detected. Venom thus appears to play little role in feeding behavior of most iguanian lizards or the powerful constricting snakes, and the low levels of expression argue against a defensive role. However, clearly the incipient or secondarily atrophied venom systems of these taxa may be a source of novel compounds useful in drug design and discovery.

Phylogeny-based comparative analysis of venom proteome variation in a clade of rattlesnakes (Sistrurus sp.)

A long-standing question in evolutionary studies of snake venoms is the extent to which phylogenetic divergence and diet can account for between-species differences in venom composition. Here we apply phylogeny-based comparative methods to address this question. We use data on venom variation generated using proteomic techniques for all members of a small clade of rattlesnakes (Sistrurus sp.) and two outgroups for which phylogenetic and diet information is available. We first complete the characterization of venom variation for all members of this clade with a “venomic” analysis of pooled venoms from two members of this genus, S. milarius streckeri and S. m. milarius. These venoms exhibit the same general classes of proteins as those found in other Sistrurus species but differ in their relative abundances of specific protein families. We then test whether there is significant phylogenetic signal in the relative abundances of major venom proteins across species and if diet (measured as percent mammals and lizards among all prey consumed) covaries with venom composition after phylogenetic divergence is accounted for. We found no evidence for significant phylogenetic signal in venom variation: K values for seven snake venom proteins and two composite venom variables [PC 1 and 2]) were all nonsignificant and lower (mean = 0.11+0.06 sd) than mean K values (>0.35) previously reported for a wide range of morphological, life history, physiological and behavioral traits from other species. Finally, analyses based on Phylogenetic Generalized Least Squares (PGLS) methods reveal that variation in abundance of some venom proteins, most strongly CRISP is significantly related to snake diet. Our results demonstrate that venom variation in these snakes is evolutionarily a highly labile trait even among very closely-related taxa and that natural selection acting through diet variation may play a role in molding the relative abundance of specific venom proteins.

Dracula's children: molecular evolution of vampire bat venom

While vampire bat oral secretions have been the subject of intense research, efforts have concentrated only on two components: DSPA (Desmodus rotundus salivary plasminogen activator) and Draculin. The molecular evolutionary history of DSPA has been elucidated, while conversely draculin has long been known from only a very small fragment and thus even the basic protein class was not even established. Despite the fact that vampire bat venom has a multitude of effects unaccounted by the documented bioactivities of DSPA and draculin, efforts have not been made to establish what other bioactive proteins are secreted by their submaxillary gland. In addition, it has remained unclear whether the anatomically distinct anterior and posterior lobes of the submaxillary gland are evolving on separate gene expression trajectories or if they remain under the shared genetic control. Using a combined proteomic and transcriptomic approach, we show that identical proteins are simultaneously expressed in both lobes. In addition to recovering the known structural classes of DSPA, we recovered a novel DSPA isoform as well as obtained a very large sequence stretch of draculin and thus established that it is a mutated version of the lactotransferrin scaffold. This study reveals a much more complex secretion profile than previously recognised. In addition to obtaining novel versions of scaffolds convergently recruited into other venoms (allergen-like, CRiSP, kallikrein, Kunitz, lysozyme), we also documented novel expression of small peptides related to calcitonin, PACAP, and statherin. Other overexpressed protein types included BPI-fold, lacritin, and secretoglobin. Further, we investigate the molecular evolution of various vampire bat venom-components and highlight the dominant role of positive selection in the evolution of these proteins. Conspicuously many of the proteins identified in the proteome were found to be homologous to proteins with known activities affecting vasodilation and platelet aggregation. We show that vampire bat venom proteins possibly evade host immune response by the mutation of the surface chemistry through focal mutagenesis under the guidance of positive Darwinian selection. These results not only contribute to the body of knowledge regarding haematophagous venoms but also provide a rich resource for novel lead compounds for use in drug design and development.

BIOLOGICAL SIGNIFICANCE

These results have direct implications in understanding the molecular evolutionary history of vampire bat venom. The unusual peptides discovered reinforce the value of studying such neglected taxon for biodiscovery.

Cloning and characterisation of three novel disintegrin precursors from the venoms of three Atheris species: Atheris chlorechis, Atheris nitschei and Atheris squamigera

Snake venom constitutes one of the most complex mixtures of naturally-occurring toxic proteins/polypeptides and a large number of these possess very profound biological activities. Disintegrins, that are commonly found in viper venoms, are low molecular weight proteins that usually contain an -Arg-Gly-Asp- (-RGD-) motif that is known to be involved in cell adhesion ligand recognition, binding specifically to cell surface integrin receptors and also exhibiting platelet anti-aggregation activity. Here, we report for the first time, the successful cloning of three cDNAs encoding disintegrin precursors from lyophilised venom-derived libraries of Atheris chlorechis, Atheris nitschei and Atheris squamigera, respectively. All of these disintegrins belong to the short-coding class and all exhibit high degrees of structural identity, both in their amino acid sequences and in the arrangement of their functional domains. Mass spectrometric analyses of the HPLC-separated/in-gel digested venom proteins was performed to characterise the mature disintegrins as expressed in the venom proteome. Studies on both the structures and conserved sites within these disintegrins are of considerable theoretical interest in the field of biological evolution and in the development of new research tools or novel templates for drug design.

Venom proteomic characterization and relative antivenom neutralization of two medically important Pakistani elapid snakes (Bungarus sindanus and Naja naja)

Intra- and interspecific variation in venom composition has been shown to have a major effect upon the efficacy of antivenoms. Due to the absence of domestically produced antivenoms, Pakistan is wholly reliant upon antivenoms produced in other countries, such as India. However, the efficacy of these antivenoms in neutralising the venoms of Pakistani snakes has not been ascertained. This is symptomatic of the general state of toxicological research in this country, which has a myriad of highly toxic and medically important venomous animals. Thus, there is a dire need for knowledge regarding the fundamental proteomics of these venoms and applied knowledge of the relative efficacy of foreign antivenoms. Here we present the results of our proteomic research on two medically important snakes of Pakistan: Bungarus sindanus and Naja naja. Indian Polyvalent Antivenom (Bharat Serums and Vaccines Ltd), which is currently marketed for use in Pakistan, was completely ineffective against either Pakistani species. In addition to the expected pre- and post-synaptic neurotoxic activity, the venom of the Pakistan population of N. naja was shown to be quite divergent from other populations of this species in being potently myotoxic. These results highlight the importance of studying divergent species and isolated populations, where the same data not only elucidates clinical problems in need of immediate attention, but also uncovers sources for novel toxins with potentially useful activities.

BIOLOGICAL SIGNIFICANCE

Pakistan Bungarus sindanus and Naja naja venoms are differentially complex. Naja naja is potently myotoxic. Neither venom is neutralized by Indian antivenom. These results have direct implications for the treatment of envenomed patients in Pakistan. The unusually myotoxic effects of Naja naja demonstrates the value of studying remote populations for biodiscovery.

Structural and phylogenetic studies with MjTX-I reveal a multi-oligomeric toxin--a novel feature in Lys49-PLA2s protein class

The mortality caused by snakebites is more damaging than many tropical diseases, such as dengue haemorrhagic fever, cholera, leishmaniasis, schistosomiasis and Chagas disease. For this reason, snakebite envenoming adversely affects health services of tropical and subtropical countries and is recognized as a neglected disease by the World Health Organization. One of the main components of snake venoms is the Lys49-phospholipases A₂, which is catalytically inactive but possesses other toxic and pharmacological activities. Preliminary studies with MjTX-I from Bothrops moojeni snake venom revealed intriguing new structural and functional characteristics compared to other bothropic Lys49-PLA₂s. We present in this article a comprehensive study with MjTX-I using several techniques, including crystallography, small angle X-ray scattering, analytical size-exclusion chromatography, dynamic light scattering, myographic studies, bioinformatics and molecular phylogenetic analyses.Based in all these experiments we demonstrated that MjTX-I is probably a unique Lys49-PLA₂, which may adopt different oligomeric forms depending on the physical-chemical environment. Furthermore, we showed that its myotoxic activity is dramatically low compared to other Lys49-PLA₂s, probably due to the novel oligomeric conformations and important mutations in the C-terminal region of the protein. The phylogenetic analysis also showed that this toxin is clearly distinct from other bothropic Lys49-PLA₂s, in conformity with the peculiar oligomeric characteristics of MjTX-I and possible emergence of new functionalities inresponse to environmental changes and adaptation to new preys.

Molecular phylogeny and evolution of the proteins encoded by coleoid (cuttlefish, octopus, and squid) posterior venom glands

In this study, we report for the first time a detailed evaluation of the phylogenetic history and molecular evolution of the major coleoid toxins: CAP, carboxypeptidase, chitinase, metalloprotease GON-domain, hyaluronidase, pacifastin, PLA2, SE-cephalotoxin and serine proteases, with the carboxypeptidase and GON-domain documented for the first time in the coleoid venom arsenal. We show that although a majority of sites in these coleoid venom-encoding genes have evolved under the regime of negative selection, a very small proportion of sites are influenced by the transient selection pressures. Moreover, nearly 70 % of these episodically adapted sites are confined to the molecular surface, highlighting the importance of variation of the toxin surface chemistry. Coleoid venoms were revealed to be as complex as other venoms that have traditionally been the recipient of the bulk of research efforts. The presence of multiple peptide/protein types in coleoids similar to those present in other animal venoms identifies a convergent strategy, revealing new information as to what characteristics make a peptide/protein type amenable for recruitment into chemical arsenals. Coleoid venoms have significant potential not only for understanding fundamental aspects of venom evolution but also as an untapped source of novel toxins for use in drug design and discovery.

Sequence analysis and phylogenetic study of some toxin proteins of snakes and related non-toxin proteins of chordates

Snakes are equipped with their venomic armory to tackle different prey and predators in adverse natural world. The venomic composition of snakes is a mix of biologically active proteins and polypeptides. Among different components snake venom cytotoxins and short neurotoxin are non-enzymatic polypeptide candidates with in the venom. These two components structurally resembled to three-finger protein superfamily specific scaffold. Different non-toxin family members of three-finger protein superfamily are involved in different biological roles. In the present study we analyzed the snake venom cytotoxins, short neurotoxins and related non-toxin proteins of different chordates in terms of amino acid sequence level diversification profile, polarity profile of amino acid sequences, conserved pattern of amino acids and phylogenetic relationship of these toxin and nontoxin protein sequences. Sequence alignment analysis demonstrates the polarity specific molecular enrichment strategy for better system adaptivity. Occurrence of amino acid substitution is high in number in toxin sequences. In non-toxin body proteins there are less amino acid substitutions. With the help of conserved residues these proteins maintain the three-finger protein scaffold. Due to system specific adaptation toxin and non-toxin proteins exhibit a varied type of amino acid residue distribution in sequence stretch. Understanding of Natural invention scheme (recruitment of venom proteins from normal body proteins) may help us to develop futuristic engineered bio-molecules with remedial properties.

The "Vampirome": Transcriptome and proteome analysis of the principal and accessory submaxillary glands of the vampire bat Desmodus rotundus, a vector of human rabies

Vampire bats are notorious for being the sole mammals that strictly feed on fresh blood for their survival. While their saliva has been historically associated with anticoagulants, only one antihemostatic (plasminogen activator) has been molecularly and functionally characterized. Here, RNAs from both principal and accessory submaxillary (submandibular) salivary glands of Desmodus rotundus were extracted, and ~200 million reads were sequenced by Illumina. The principal gland was enriched with plasminogen activators with fibrinolytic properties, members of lipocalin and secretoglobin families, which bind prohemostatic prostaglandins, and endonucleases, which cleave neutrophil-derived procoagulant NETs. Anticoagulant (tissue factor pathway inhibitor, TFPI), vasodilators (PACAP and C-natriuretic peptide), and metalloproteases (ADAMTS-1) were also abundantly expressed. Members of the TSG-6 (anti-inflammatory), antigen 5/CRISP, and CCL28-like (antimicrobial) protein families were also sequenced. Apyrases (which remove platelet agonist ADP), phosphatases (which degrade procoagulant polyphosphates), and sphingomyelinase were found at lower transcriptional levels. Accessory glands were enriched with antimicrobials (lysozyme, defensin, lactotransferrin) and protease inhibitors (TIL-domain, cystatin, Kazal). Mucins, heme-oxygenase, and IgG chains were present in both glands. Proteome analysis by nano LC-MS/MS confirmed that several transcripts are expressed in the glands. The database presented herein is accessible online at http://exon.niaid.nih.gov/transcriptome/D_rotundus/Supplemental-web.xlsx. These results reveal that bat saliva emerges as a novel source of modulators of vascular biology.

BIOLOGICAL SIGNIFICANCE

Vampire bat saliva emerges as a novel source of antihemostatics which modulate several aspects of vascular biology.

Structural characteristics and evolution of the Protobothrops elegans pancreatic phospholipase A2 gene in contrast with those of Protobothrops genus venom phospholipase A2 genes

The nucleotide sequence of the gene encoding Protobothrops elegans (Crotalinae) pancreatic phospholipase A(2) (PLA(2)), abbreviated PePancPLA(2), was determined by means of inverted PCR techniques. Since its deduced amino acid sequence contains a pancreatic loop and shows high similarity to that of Laticauda semifasciata (Elapinae) group IB pancreatic PLA(2), PePancPLA(2) is classified into group IB PLA(2). The nucleotide sequences of the PePancPLA(2) gene, the L. semifasciata group IB pancreatic PLA(2) gene, and the L. semifasciata group IA venom PLA(2) gene are similar to one another but greatly dissimilar to those of Protobothrops genus (Crotalinae) group II venom PLA(2) genes, suggesting that the Elapinae group IB PLA(2) gene and the group IA PLA(2) gene appeared after Elapinae was established, and that the Crotalinae group II venom PLA(2) genes came into existence before Elapinae and Crotalinae diverged. A phylogenetic analysis of their amino acid sequences confirms this.

Differential evolution and neofunctionalization of snake venom metalloprotease domains

Snake venom metalloproteases (SVMP) are composed of five domains: signal peptide, propeptide, metalloprotease, disintegrin, and cysteine-rich. Secreted toxins are typically combinatorial variations of the latter three domains. The SVMP-encoding genes of Psammophis mossambicus venom are unique in containing only the signal and propeptide domains. We show that the Psammophis SVMP propeptide evolves rapidly and is subject to a high degree of positive selection. Unlike Psammophis, some species of Echis express both the typical multidomain and the unusual monodomain (propeptide only) SVMP, with the result that a lower level of variation is exerted upon the latter. We showed that most mutations in the multidomain Echis SVMP occurred in the protease domain responsible for proteolytic and hemorrhagic activities. The cysteine-rich and disintegrin-like domains, which are putatively responsible for making the P-III SVMPs more potent than the P-I and P-II forms, accumulate the remaining variation. Thus, the binding sites on the molecule's surface are evolving rapidly whereas the core remains relatively conserved. Bioassays conducted on two post-translationally cleaved novel proline-rich peptides from the P. mossambicus propeptide domain showed them to have been neofunctionalized for specific inhibition of mammalian a7 neuronal nicotinic acetylcholine receptors. We show that the proline rich postsynaptic specific neurotoxic peptides from Azemiops feae are the result of convergent evolution within the precursor region of the C-type natriuretic peptide instead of the SVMP. The results of this study reinforce the value of studying obscure venoms for biodiscovery of novel investigational ligands.

Complex cocktails: the evolutionary novelty of venoms

Venoms have evolved on numerous occasions throughout the animal kingdom. These 'biochemical weapon systems' typically function to facilitate, or protect the producing animal from, predation. Most venomous animals remain unstudied despite venoms providing model systems for investigating predator-prey interactions, molecular evolution, functional convergence, and novel targets for pharmaceutical discovery. Through advances in 'omic' technologies, venom composition data have recently become available for several venomous lineages, revealing considerable complexity in the processes responsible for generating the genetic and functional diversity observed in many venoms. Here, we review these recent advances and highlight the ecological and evolutionary novelty of venom systems.