What's in a mass?

This short essay pretends to make the reader reflect on the concept of biological mass and on the added value that the determination of this molecular property of a protein brings to the interpretation of evolutionary and translational snake venomics research. Starting from the premise that the amino acid sequence is the most distinctive primary molecular characteristics of any protein, the thesis underlying the first part of this essay is that the isotopic distribution of a protein's molecular mass serves to unambiguously differentiate it from any other of an organism's proteome. In the second part of the essay, we discuss examples of collaborative projects among our laboratories, where mass profiling of snake venom PLA2 across conspecific populations played a key role revealing dispersal routes that determined the current phylogeographic pattern of the species.

Danger in the Canopy. Comparative Proteomics and Bioactivities of the Venoms of the South American Palm Pit Viper Bothrops bilineatus Subspecies bilineatus and smaragdinus and Antivenomics of B. b. bilineatus (Rondônia) Venom against the Brazilian Pentabothropic Antivenom.

We report a structural and functional proteomics characterization of venoms of the two subspecies (Bothrops bilineatusbilineatus and B. b. smaragdinus) of the South American palm pit viper from the Brazilian state of Rondônia and B. b. smaragdinus from Perú. These poorly known arboreal and mostly nocturnal generalist predators are widely distributed in lowland rainforests throughout the entire Amazon region, where they represent an important cause of snakebites. The three B. bilineatus spp. venom samples exhibit overall conserved proteomic profiles comprising components belonging to 11 venom protein classes, with PIII (34-40% of the total venom proteins) and PI (8-18%) SVMPs and their endogenous tripeptide inhibitors (SVMPi, 8-10%); bradykinin-potentiating-like peptides (BBPs, 10.7-15%); snake venom serine proteinases (SVSP, 5.5-14%); C-type lectin-like proteins (CTL, 3-10%); phospholipases A2 (PLA2, 2.8-7.6%); cysteine-rich secretory proteins (CRISP, 0.9-2.8%); l-amino acid oxidases (LAO, 0.9-5%) representing the major components of their common venom proteomes. Comparative analysis of the venom proteomes of the two geographic variants of B. b. smaragdinus with that of B. b. bilineatus revealed that the two Brazilian taxa share identical molecules between themselves but not with Peruvian B. b. smaragdinus, suggesting hybridization between the geographically close, possibly sympatric, Porto Velho (RO, BR) B. b. smaragdinus and B. b. bilineatus parental populations. However, limited sampling does not allow determining the frequency of this event. The toxin arsenal of the South American palm pit vipers may account for the in vitro recorded collagenolytic, caseinolytic, PLA2, l-amino acid oxidase, thrombin-like and factor X-activating activities, and the clinical features of South American palm pit viper envenomings, i.e., local and progressively ascending pain, shock and loss of consciousness, spontaneous bleeding, and profound coagulopathy. The remarkable cross-reactivity of the Brazilian pentabothropic SAB antivenom toward the heterologous B. b. bilineatus venom suggests that the paraspecific antigenic determinants should have been already present in the venom of the last common ancestor of the Bothrops ″jararaca″ and ″taeniatus″ clades, about 8.5 Mya in the mid-late Miocene epoch of the Cenozoic era. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the data set identifiers PXD020043, PXD020026, and PXD020013.

The molecular basis of venom resistance in a rattlesnake-squirrel predator-prey system.

Understanding how interspecific interactions mould the molecular basis of adaptations in coevolving species is a long-sought goal of evolutionary biology. Venom in predators and venom resistance proteins in prey are coevolving molecular phenotypes, and while venoms are highly complex mixtures it is unclear if prey respond with equally complex resistance traits. Here, we use a novel molecular methodology based on protein affinity columns to capture and identify candidate blood serum resistance proteins ("venom interactive proteins" [VIPs]) in California Ground Squirrels (Otospermophilus beecheyi) that interact with venom proteins from their main predator, Northern Pacific Rattlesnakes (Crotalus o. oreganus). This assay showed that serum-based resistance is both population- and species-specific, with serum proteins from ground squirrels showing higher binding affinities for venom proteins of local snakes compared to allopatric individuals. Venom protein specificity assays identified numerous and diverse candidate prey resistance VIPs but also potential targets of venom in prey tissues. Many specific VIPs bind to multiple snake venom proteins and, conversely, single venom proteins bind multiple VIPs, demonstrating that a portion of the squirrel blood serum "resistome" involves broad-based inhibition of nonself proteins and suggests that resistance involves a toxin scavenging mechanism. Analyses of rates of evolution of VIP protein homologues in related mammals show that most of these proteins evolve under purifying selection possibly due to molecular constraints that limit the evolutionary responses of prey to rapidly evolving snake venom proteins. Our method represents a general approach to identify specific proteins involved in co-evolutionary interactions between species at the molecular level.

When one phenotype is not enough: divergent evolutionary trajectories govern venom variation in a widespread rattlesnake species.

Understanding the origin and maintenance of phenotypic variation, particularly across a continuous spatial distribution, represents a key challenge in evolutionary biology. For this, animal venoms represent ideal study systems: they are complex, variable, yet easily quantifiable molecular phenotypes with a clear function. Rattlesnakes display tremendous variation in their venom composition, mostly through strongly dichotomous venom strategies, which may even coexist within a single species. Here, through dense, widespread population-level sampling of the Mojave rattlesnake, Crotalus scutulatus, we show that genomic structural variation at multiple loci underlies extreme geographical variation in venom composition, which is maintained despite extensive gene flow. Unexpectedly, neither diet composition nor neutral population structure explain venom variation. Instead, venom divergence is strongly correlated with environmental conditions. Individual toxin genes correlate with distinct environmental factors, suggesting that different selective pressures can act on individual loci independently of their co-expression patterns or genomic proximity. Our results challenge common assumptions about diet composition as the key selective driver of snake venom evolution and emphasize how the interplay between genomic architecture and local-scale spatial heterogeneity in selective pressures may facilitate the retention of adaptive functional polymorphisms across a continuous space.

Integrated Venomics and Venom Gland Transcriptome Analysis of Juvenile and Adult Mexican Rattlesnakes Crotalus simus, C. tzabcan, and C. culminatus Revealed miRNA-modulated Ontogenetic Shifts.

Adult rattlesnakes within genus Crotalus express one of two distinct venom phenotypes, type I (hemorrhagic) and type II (neurotoxic). In Costa Rican Central American rattlesnake, ontogenetic changes in the concentration of miRNAs modulate venom type II to type I transition. Venomics and venom gland transcriptome analyses showed that adult C. simus and C. tzabcan expressed intermediate patterns between type II and type I venoms, whereas C. culminatus had a canonical type I venom. Neonate/juvenile and adult Mexican rattlesnakes showed notable inter- and intraspecific variability in the number, type, abundance and ontogenetic shifts of the transcriptional and translational venom gland activities. These results support a role for miRNAs in the ontogenetic venom compositional changes in the three congeneric Mexican rattlesnakes. It is worth noting the finding of dual-action miRNAs, which silence the translation of neurotoxic heterodimeric PLA2 crotoxin and acidic PLA2 mRNAs while simultaneously up-regulating SVMP-targeting mRNAs. Dual transcriptional regulation potentially explains the existence of mutually exclusive crotoxin-rich (type-II) and SVMP-rich (type-I) venom phenotypic dichotomy among rattlesnakes. Our results support the hypothesis that alterations of the distribution of miRNAs, modulating the translational activity of venom gland toxin-encoding mRNAs in response to an external cue, may contribute to the mechanism generating adaptive venom variability.

What killed Karl Patterson Schmidt? Combined venom gland transcriptomic, venomic and antivenomic analysis of the South African green tree snake (the boomslang), Dispholidus typus.

BACKGROUND: Non-front-fanged colubroid snakes comprise about two-thirds of extant ophidian species. The medical significance of the majority of these snakes is unknown, but at least five species have caused life-threatening or fatal human envenomings. However, the venoms of only a small number of species have been explored. METHODS: A combined venomic and venom gland transcriptomic approach was employed to characterise of venom of Dispholidus typus (boomslang), the snake that caused the tragic death of Professor Karl Patterson Schmidt. The ability of CroFab™ antivenom to immunocapture boomslang venom proteins was investigated using antivenomics. RESULTS: Transcriptomic-assisted proteomic analysis identified venom proteins belonging to seven protein families: three-finger toxin (3FTx); phospholipase A2 (PLA2); cysteine-rich secretory proteins (CRISP); snake venom (SV) serine proteinase (SP); C-type lectin-like (CTL); SV metalloproteinases (SVMPs); and disintegrin-like/cysteine-rich (DC) proteolytic fragments. CroFab™ antivenom efficiently immunodepleted some boomslang SVMPs. CONCLUSIONS: The present work is the first to address the overall proteomic profile of D. typus venom. This study allowed us to correlate the toxin composition with the toxic activities of the venom. The antivenomic analysis suggested that the antivenom available at the time of the unfortunate accident could have exhibited at least some immunoreactivity against the boomslang SVMPs responsible for the disseminated intravascular coagulation syndrome that caused K.P. Schmidt's fatal outcome. GENERAL SIGNIFICANCE: This study may stimulate further research on other non-front-fanged colubroid snake venoms capable of causing life-threatening envenomings to humans, which in turn should contribute to prevent fatal human accidents, such as that unfortunately suffered by K.P. Schmidt.

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.

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.

Constructing comprehensive venom proteome reference maps for integrative venomics.

OBJECTIVE: Understanding the molecular basis of complex adaptive traits, such as snake venom, demands qualitative and quantitative comparisons of the temporal and spatial patterns of venom variation. Here, we assessed the proof-of-concept that locus-resolved reference venom proteome maps can be achieved through efficient pre-MS venom proteome decomplexation, peptide-centric MS/MS analysis and species-specific database searching. METHODS: Venom proteome components were fractionated and quantified by RP-HPLC, SDS-PAGE and 2DE prior to LC-MS/MS matching against a species-specific transcriptomic dataset. RESULTS: Combination of RP-HPLC/SDS-PAGE and 2DE followed by LC-MS/MS showed the existence of ∼178-180 venom protein species generated from ∼48 unique transcripts. CONCLUSIONS: Our results underscore that if sufficient pre-MS and MS efforts are applied, comprehensive venom maps can be achieved. And - equally important - dissociating the venom decomplexing steps from the protein identification process represents the key to achieving a quantitative and locus-resolved insight of the venom proteome.

Integrated "omics" profiling indicates that miRNAs are modulators of the ontogenetic venom composition shift in the Central American rattlesnake, Crotalus simus simus.

BACKGROUND: Understanding the processes that drive the evolution of snake venom is a topic of great research interest in molecular and evolutionary toxinology. Recent studies suggest that ontogenetic changes in venom composition are genetically controlled rather than environmentally induced. However, the molecular mechanisms underlying these changes remain elusive. Here we have explored the basis and level of regulation of the ontogenetic shift in the venom composition of the Central American rattlesnake, Crotalus s. simus using a combined proteomics and transcriptomics approach. RESULTS: Proteomic analysis showed that the ontogenetic shift in the venom composition of C. s. simus is essentially characterized by a gradual reduction in the expression of serine proteinases and PLA2 molecules, particularly crotoxin, a ß-neurotoxic heterodimeric PLA2, concominantly with an increment of PI and PIII metalloproteinases at age 9-18 months. Comparison of the transcriptional activity of the venom glands of neonate and adult C. s. simus specimens indicated that their transcriptomes exhibit indistinguisable toxin family profiles, suggesting that the elusive mechanism by which shared transcriptomes generate divergent venom phenotypes may operate post-transcriptionally. Specifically, miRNAs with frequency count of 1000 or greater exhibited an uneven distribution between the newborn and adult datasets. Of note, 590 copies of a miRNA targeting crotoxin B-subunit was exclusively found in the transcriptome of the adult snake, whereas 1185 copies of a miRNA complementary to a PIII-SVMP mRNA was uniquely present in the newborn dataset. These results support the view that age-dependent changes in the concentration of miRNA modulating the transition from a crotoxin-rich to a SVMP-rich venom from birth through adulthood can potentially explain what is observed in the proteomic analysis of the ontogenetic changes in the venom composition of C. s. simus. CONCLUSIONS: Existing snake venom toxins are the result of early recruitment events in the Toxicofera clade of reptiles by which ordinary genes were duplicated, and the new genes selectively expressed in the venom gland and amplified to multigene families with extensive neofunctionalization throughout the approximately 112-125 million years of ophidian evolution. Our findings support the view that understanding the phenotypic diversity of snake venoms requires a deep knowledge of the mechanisms regulating the transcriptional and translational activity of the venom gland. Our results suggest a functional role for miRNAs. The impact of specific miRNAs in the modulation of venom composition, and the integration of the mechanisms responsible for the generation of these miRNAs in the evolutionary landscape of the snake's venom gland, are further challenges for future research.

Snake venomics of Crotalus tigris: the minimalist toxin arsenal of the deadliest Nearctic rattlesnake venom. Evolutionary Clues for generating a pan-specific antivenom against crotalid type II venoms [corrected].

We report the proteomic and antivenomic characterization of Crotalus tigris venom. This venom exhibits the highest lethality for mice among rattlesnakes and the simplest toxin proteome reported to date. The venom proteome of C. tigris comprises 7-8 gene products from 6 toxin families; the presynaptic ß-neurotoxic heterodimeric PLA(2), Mojave toxin, and two serine proteinases comprise, respectively, 66 and 27% of the C. tigris toxin arsenal, whereas a VEGF-like protein, a CRISP molecule, a medium-sized disintegrin, and 1-2 PIII-SVMPs each represent 0.1-5% of the total venom proteome. This toxin profile really explains the systemic neuro- and myotoxic effects observed in envenomated animals. In addition, we found that venom lethality of C. tigris and other North American rattlesnake type II venoms correlates with the concentration of Mojave toxin A-subunit, supporting the view that the neurotoxic venom phenotype of crotalid type II venoms may be described as a single-allele adaptation. Our data suggest that the evolutionary trend toward neurotoxicity, which has been also reported for the South American rattlesnakes, may have resulted by pedomorphism. The ability of an experimental antivenom to effectively immunodeplete proteins from the type II venoms of C. tigris, Crotalus horridus , Crotalus oreganus helleri, Crotalus scutulatus scutulatus, and Sistrurus catenatus catenatus indicated the feasibility of generating a pan-American anti-Crotalus type II antivenom, suggested by the identification of shared evolutionary trends among South and North American Crotalus species.

Interpopulational variation and ontogenetic shift in the venom composition of Lataste's viper (Vipera latastei, Boscá 1878) from northern Portugal.

Lataste's viper (Vipera latastei) is a venomous European viper endemic to the Iberian Peninsula, recognised as medically important by the World Health Organization. To date, no comprehensive characterisation of this species' venom has been reported. Here, we analysed the venoms of juvenile and adult specimens of V. latastei from two environmentally different populations from northern Portugal. Using bottom-up venomics, we produced six venom proteomes (three per population) from vipers belonging to both age classes (i.e., two juveniles and four adults), and RP-HPLC profiles of 54 venoms collected from wild specimens. Venoms from juveniles and adults differed in their chromatographic profiles and relative abundances of their toxins, suggesting the occurrence of ontogenetic changes in venom composition. Specifically, snake venom metalloproteinase (SVMP) was the most abundant toxin family in juvenile venoms, while snake venom serine proteinases (SVSPs), phospholipases A2 (PLA2s), and C-type lectin-like (CTLs) proteins were the main toxins comprising adult venoms. The RP-HPLC venom profiles were found to vary significantly between the two sampled localities, indicating geographic variability. Furthermore, the presence/absence of certain peaks in the venom chromatographic profiles appeared to be significantly correlated also to factors like body size and sex of the vipers. Our findings show that V. latastei venom is a variable phenotype. The intraspecific differences we detected in its composition likely mirror changes in the feeding ecology of this species, taking place during different life stages and under different environmental pressures. SIGNIFICANCE: Lataste's viper (Vipera latastei) is a medically important viper endemic to the Iberian Peninsula, inhabiting different habitats and undergoing a marked ontogenetic dietary shift. In the current study, we report the first proteomic analysis of V. latastei venom from two environmentally different localities in northern Portugal. Our bottom-up venomic analyses show that snake venom serine proteinases (SVSPs), phospholipases A2 (PLA2s), and C-type lectin-like (CTLs) proteins are the major components of adult V. latastei venom. The comparative analysis of young and adult venoms suggests the occurrence of ontogenetic shift in toxin abundances, with snake venom metalloproteinases (SVMPs) being the predominant toxins in juvenile venoms. Moreover, geographic venom variation between the two studied populations is also detected, with our statistical analyses suggesting that factors like body size and sex of the vipers are possibly at play in its determination. Our work represents the first assessment of the composition of V. latastei venom, and the first step towards a better understanding of the drivers behind its variability.

Monoclonal-Based Antivenomics Reveals Conserved Neutralizing Epitopes in Type I PLA2 Molecules from Coral Snakes.

For over a century, polyclonal antibodies have been used to treat snakebite envenoming and are still considered by the WHO as the only scientifically validated treatment for snakebites. Nevertheless, moderate innovations have been introduced to this immunotherapy. New strategies and approaches to understanding how antibodies recognize and neutralize snake toxins represent a challenge for next-generation antivenoms. The neurotoxic activity of Micrurus venom is mainly due to two distinct protein families, three-finger toxins (3FTx) and phospholipases A2 (PLA2). Structural conservation among protein family members may represent an opportunity to generate neutralizing monoclonal antibodies (mAbs) against family-conserved epitopes. In this work, we sought to produce a set of monoclonal antibodies against the most toxic components of M. altirostris venom. To this end, the crude venom was fractionated, and its major toxic proteins were identified and used to generate a panel of five mAbs. The specificity of these mAbs was characterized by ELISA and antivenomics approaches. Two of the generated mAbs recognized PLA2 epitopes. They inhibited PLA2 catalytic activity and showed paraspecific neutralization against the myotoxicity from the lethal effect of Micrurus and Naja venoms' PLA2s. Epitope conservation among venom PLA2 molecules suggests the possibility of generating pan-PLA2 neutralizing antibodies.

Venomic and antivenomic analyses of the Central American coral snake, Micrurus nigrocinctus (Elapidae).

The proteome of the venom of Micrurus nigrocinctus (Central American coral snake) was analyzed by a "venomics" approach. Nearly 50 venom peaks were resolved by RP-HPLC, revealing a complex protein composition. Comparative analyses of venoms from individual specimens revealed that such complexity is an intrinsic feature of this species, rather than the sum of variable individual patterns of simpler composition. Proteins related to eight distinct families were identified by MS/MS de novo peptide sequencing or N-terminal sequencing: phospholipase A(2) (PLA(2)), three-finger toxin (3FTx), l-amino acid oxidase, C-type lectin/lectin-like, metalloproteinase, serine proteinase, ohanin, and nucleotidase. PLA(2)s and 3FTxs are predominant, representing 48 and 38% of the venom proteins, respectively. Within 3FTxs, several isoforms of short-chain α-neurotoxins as well as muscarinic-like toxins and proteins with similarity to long-chain κ-2 bungarotoxin were identified. PLA(2)s are also highly diverse, and a toxicity screening showed that they mainly exert myotoxicity, although some are lethal and may contribute to the known presynaptic neurotoxicity of this venom. An antivenomic characterization of a therapeutic monospecific M. nigrocinctus equine antivenom revealed differences in immunorecognition of venom proteins that correlate with their molecular mass, with the weakest recognition observed toward 3FTxs.

Snake venomics of African spitting cobras: toxin composition and assessment of congeneric cross-reactivity of the pan-African EchiTAb-Plus-ICP antivenom by antivenomics and neutralization approaches.

Venomic analysis of the venoms of Naja nigricollis, N. katiensis, N. nubiae, N. mossambica, and N. pallida revealed similar compositional trends. The high content of cytotoxins and PLA(2)s may account for the extensive tissue necrosis characteristic of the envenomings by these species. The high abundance of a type I α-neurotoxin in N. nubiae may be responsible for the high lethal toxicity of this venom (in rodents). The ability of EchiTAb-Plus-ICP antivenom to immunodeplete and neutralize the venoms of African spitting cobras was assessed by antivenomics and neutralization tests. It partially immunodepleted 3FTx and PLA(2)s and completely immunodepleted SVMPs and CRISPs in all venoms. The antivenom neutralized the dermonecrotic and PLA(2) activities of all African Naja venoms, whereas lethality was eliminated in the venoms of N. nigricollis, N. mossambica, and N. pallida but not in those of N. nubiae and N. katiensis. The lack of neutralization of lethality of N. nubiae venom may be of medical relevance only in relatively populous areas of the Saharan region. The impaired activity of EchiTAb-Plus-ICP against N. katiensis may not represent a major concern. This species is sympatric with N. nigricollis in many regions of Africa, although very few bites have been attributed to it.

NMR structure and dynamics of recombinant wild type and mutated jerdostatin, a selective inhibitor of integrin α1ß1.

NMR analysis of four recombinant jerdostatin molecules was assessed to define the structural basis of two naturally occurring gain-of-function events: C-terminal dipeptide processing and mutation of the active residue K21 to arginine. Removal of the highly mobile and a bulky C-terminal dipeptide produced pronounced chemical shift changes in the sequentially unconnected but spatially nearby α(1)ß(1) inhibitory loop. Analysis of chemical shift divergence and (15)N backbone relaxation dynamics indicated differences in motions in the picosecond to nanosecond time scale, and the higher T(2) rate of S25, S26, and H27 of rJerK21 point to a slowdown in the microsecond to millisecond motions of these residues when compared with rJerR21. The evidence presented in this article converges on the hypothesis that dynamic differences between the α(1)ß(1) recognition loops of rJerR21 and rJerK21 may influence the thermodynamics of their receptor recognition and binding. A decrease in the µs-ms time scale may impair the binding affinity by reducing the rate of possible conformations that the rJerK21 can adopt in this time scale.

Profiling the venom gland transcriptomes of Costa Rican snakes by 454 pyrosequencing.

BACKGROUND: A long term research goal of venomics, of applied importance for improving current antivenom therapy, but also for drug discovery, is to understand the pharmacological potential of venoms. Individually or combined, proteomic and transcriptomic studies have demonstrated their feasibility to explore in depth the molecular diversity of venoms. In the absence of genome sequence, transcriptomes represent also valuable searchable databases for proteomic projects. RESULTS: The venom gland transcriptomes of 8 Costa Rican taxa from 5 genera (Crotalus, Bothrops, Atropoides, Cerrophidion, and Bothriechis) of pitvipers were investigated using high-throughput 454 pyrosequencing. 100,394 out of 330,010 masked reads produced significant hits in the available databases. 5.165,220 nucleotides (8.27%) were masked by RepeatMasker, the vast majority of which corresponding to class I (retroelements) and class II (DNA transposons) mobile elements. BLAST hits included 79,991 matches to entries of the taxonomic suborder Serpentes, of which 62,433 displayed similarity to documented venom proteins. Strong discrepancies between the transcriptome-computed and the proteome-gathered toxin compositions were obvious at first sight. Although the reasons underlaying this discrepancy are elusive, since no clear trend within or between species is apparent, the data indicate that individual mRNA species may be translationally controlled in a species-dependent manner. The minimum number of genes from each toxin family transcribed into the venom gland transcriptome of each species was calculated from multiple alignments of reads matched to a full-length reference sequence of each toxin family. Reads encoding ORF regions of Kazal-type inhibitor-like proteins were uniquely found in Bothriechis schlegelii and B. lateralis transcriptomes, suggesting a genus-specific recruitment event during the early-Middle Miocene. A transcriptome-based cladogram supports the large divergence between A. mexicanus and A. picadoi, and a closer kinship between A. mexicanus and C. godmani. CONCLUSIONS: Our comparative next-generation sequencing (NGS) analysis reveals taxon-specific trends governing the formulation of the venom arsenal. Knowledge of the venom proteome provides hints on the translation efficiency of toxin-coding transcripts, contributing thereby to a more accurate interpretation of the transcriptome. The application of NGS to the analysis of snake venom transcriptomes, may represent the tool for opening the door to systems venomics.

The earless monitor lizard Lanthanotus borneensis - A venomous animal?

Based on its mandibular gland secretion, the earless monitor lizard, Lanthanotus borneensis, has been considered a venomous animal like other members of the Toxicofera group, including Heloderma. In the present study, the gland structure and teeth of L. borneensis were examined by micro-tomography (µCT) and scanning electron microscopy (SEM), respectively, and proteomic analysis of the gland extract was performed. The mandibular gland consists of six compartments with separate ducts. The pleurodont teeth of the lower and upper jaw are not grooved but possess a sharp ridge on the anterior surface. Proteomic analysis of the gland extract confirmed previous studies that kallikrein enzymes are the major biologically active components. In view of the lizard's biology, its mandibular gland secretion is obviously not needed for prey capture or defence. It seems not justified the labelling of L. borneensis as a venomous animal. However, definitively answering this question requires toxinological studies on natural prey.

Venomics of the poorly studied hognosed pitvipers Porthidium arcosae and Porthidium volcanicum.

We report the first proteomics analyses of the venoms of two poorly studied snakes, the Manabi hognosed pitviper Porthidium arcosae endemic to the western coastal province of Manabí (Ecuador), and the Costa Rican hognosed pitviper P. volcanicum with distribution restricted to South Pacific Costa Rica and western Panamá. These venom proteomes share a conserved compositional pattern reported in four other congeneric species within the clade of South American Porthidium species, P. nasutum, P. lansbergii, P. ophryomegas, and P. porrasi. The paraspecific immunorecognition profile of antivenoms produced in Costa Rica (ICP polyvalent), Perú (Instituto Nacional de Salud) and Brazil (soro antibotrópico pentavalente, SAB, from Instituto Butantan) against the venom of P. arcosae was investigated through a third-generation antivenomics approach. The maximal venom-binding capacities of the investigated antivenoms were 97.1 mg, 21.8 mg, and 25.7 mg of P. arcosae venom proteins per gram of SAB, ICP, and INS-PERU antibody molecules, respectively, which translate into 28.4 mg, 13.1 mg, and 15.2 mg of total venom proteins bound per vial of SAB, ICP, and INS-PERU AV. The antivenomics results suggest that 21.8%, 7.8% and 6.1% of the SAB, ICP, and INS-PERU antibody molecules recognized P. arcosae venom toxins. The SAB antivenom neutralized P. arcosae venom's lethality in mice with an ED50 of 31.3 mgV/g SAB AV. This preclinical neutralization paraspecificity points to Brazilian SAB as a promising candidate for the treatment of envenomings by Ecuadorian P. arcosae. BIOLOGICAL SIGNIFICANCE: Assessing the preclinical efficacy profile of antivenoms against homologous and heterologous medically relevant snake venoms represents an important goal towards defining the biogeographic range of their clinical utility. This is particularly relevant in regions, such as Mesoamerica, where a small number of pharmaceutical companies produce antivenoms against the venoms of a small number of species of maximum medical relevance among the local rich herpetofauna, leaving a wide range of snakes of secondary medical relevance, but also causing life-threatening human envenomings without nominal clinical coverage. This work is part of a larger project aiming at mapping the immunological characteristics of antivenoms generated in Latin American countries towards venoms of such poorly studied snakes of the local and neighboring countries' herpetofauna. Here we report the proteomics characterization of the Manabi hognosed pitviper Porthidium arcosae endemic to the western coastal province of Manabí (Ecuador), and the Costa Rican hognosed pitviper P. volcanicum with distribution restricted to southwestern Costa Rica, the antivenomics assessment of three bothropoid commercial antivenoms produced in Costa Rica, Perú, and Brazil against the venom components of P. arcosae, and the in vivo capacity of the Brazilian soro antibotrópico pentavalente (SAB) from Instituto Butantan to neutralize the murine lethality of P. arcosae venom. The preclinical paraspecific ED50 of 31.3 mg of P. arcosae venom per gram of antivenom points to Brazilian SAB as a promising candidate for the treatment of envenomings by the Manabi hognosed pitviper P. arcosae.