Bitis gabonica (Gaboon viper) snake venom gland: toward a catalog for the full-length transcripts (cDNA) and proteins

High-throughput proteomics and in vitro functional characterization of the 26 medically most important elapids and vipers from sub-Saharan Africa

Venomous snakes are important parts of the ecosystem, and their behavior and evolution have been shaped by their surrounding environments over the eons. This is reflected in their venoms, which are typically highly adapted for their biological niche, including their diet and defense mechanisms for deterring predators. Sub-Saharan Africa is rich in venomous snake species, of which many are dangerous to humans due to the high toxicity of their venoms and their ability to effectively deliver large amounts of venom into their victims via their bite. In this study, the venoms of 26 of sub-Saharan Africa's medically most relevant elapid and viper species were subjected to parallelized toxicovenomics analysis. The analysis included venom proteomics and in vitro functional characterization of whole venom toxicities, enabling a robust comparison of venom profiles between species. The data presented here corroborate previous studies and provide biochemical details for the clinical manifestations observed in envenomings by the 26 snake species. Moreover, two new venom proteomes (Naja anchietae and Echis leucogaster) are presented here for the first time. Combined, the presented data can help shine light on snake venom evolutionary trends and possibly be used to further improve or develop novel antivenoms.

Bibliometric Analysis of Literature in Snake Venom-Related Research Worldwide (1933–2022)

Simple Summary

Around the world, snake envenomation poses a serious health risk. Proteins with pharmacological effects are present in snake venom. Recent studies elaborate snake venom and its potential application, including as a cancer drug and antibacterial substances. Our study aimed to analyze the global profile of the literature in snake venom research from documents indexed in the Scopus database between 1933 and 2022. In total, 2999 documents were published with Brazil showing the highest productivity. Antivenom, proteomics, and transcriptomics are emerging as hot topics on a global scale. The present study offers a distinctive overview of snake venom research conducted worldwide.

Abstract

Snake envenomation is a severe economic and health concern affecting countries worldwide. Snake venom carries a wide variety of small peptides and proteins with various immunological and pharmacological properties. A few key research areas related to snake venom, including its applications in treating cancer and eradicating antibiotic-resistant bacteria, have been gaining significant attention in recent years. The goal of the current study was to analyze the global profile of literature in snake venom research. This study presents a bibliometric review of snake venom-related research documents indexed in the Scopus database between 1933 and 2022. The overall number of documents published on a global scale was 2999, with an average annual production of 34 documents. Brazil produced the highest number of documents (n = 729), followed by the United States (n = 548), Australia (n = 240), and Costa Rica (n = 235). Since 1963, the number of publications has been steadily increasing globally. At a worldwide level, antivenom, proteomics, and transcriptomics are growing hot issues for research in this field. The current research provides a unique overview of snake venom research at global level from 1933 through 2022, and it may be beneficial in guiding future research.

Synthesis and Pharmacological Characterization of Visabron, a Backbone Cyclic Peptide Dual Antagonist of α4β1 (VLA-4)/α9β1 Integrin for Therapy of Multiple Sclerosis

Integrins α4β1/ α9β1 are important in the pathogenesis and progression of inflammatory and autoimmune diseases by their roles in leukocyte activation and trafficking. Natalizumab, a monoclonal antibody selectively targeting α4β1 integrin and blocking leukocyte trafficking to the central nervous system, is an immunotherapy for multiple sclerosis (MS). However, due to its adverse effects associated with chronic treatment, alternative strategies using small peptide mimetic inhibitors are being sought. In the present study, we synthesized and characterized visabron c (4-4), a backbone cyclic octapeptide based on the sequence TMLD, a non-RGD unique α4β1 integrin recognition sequence motif derived from visabres, a proteinous disintegrin from the viper venom. Visabron c (4-4) was selected from a minilibrary with conformational diversity based on its potency and selectivity in functional adhesion cellular assays. Visabron c (4-4)'s serum stability, pharmacokinetics, and therapeutic effects following ip injection were assessed in an experimental autoimmune encephalomyelitis (EAE) animal model. Furthermore, visabron c (4-4)'s lack of toxic effects in mice was verified by blood analysis, tissue pathology, immunogenicity, and "off-target" effects, indicating its significant tolerability and lack of immunogenicity. Visabron c (4-4) can be delivered systemically. The in vitro and in vivo data justify visabron c (4-4) as a safe alternative peptidomimetic lead compound/drug to monoclonal anti-α4 integrin antibodies, steroids, and other immunosuppressant drugs. Moreover, visabron c (4-4) design may pave the way for developing new therapies for a variety of other inflammatory and/or autoimmune diseases.

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

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

Snake venom vascular endothelial growth factors (svVEGFs): Unravelling their molecular structure, functions, and research potential

Vascular endothelial growth factor (VEGF) is a key regulator of angiogenesis, a physiological process characterized by the formation of new vessels from a preexisting endothelium. VEGF has also been implicated in pathologic states, such as neoplasias, intraocular neovascular disorders, among other conditions. VEGFs are distributed in seven different families: VEGF-A, B, C, D, and PIGF (placental growth factor), which are identified in mammals; VEGF-E, which are encountered in viruses; and VEGF-F or svVEGF (snake venom VEGF) described in snake venoms. This is the pioneer review of svVEGF family, exploring its distribution among the snake venoms, molecular structure, main functions, and potential applications.

Selected Emerging Infectious Diseases of Squamata: An Update

This article details emerging infectious diseases that have devastating impacts on captive and wild squamates. Treatment advances have been attempted for Cryptosporidium infections in squamates. Gram-positive bacteria, Devriesea agamarum and Austwickia chelonae, are contributing to severe disease in captive and now in wild reptiles, some critically endangered. Nannizziposis, Paranannizziopsis, and Ophidiomyces continue to cause fatal disease as primary pathogens in wild and captive populations of squamates and sphenodontids. Nidovirus, bornavirus, paramyxovirus, sunshine virus, and arenavirus have emerged to be significant causes of neurorespiratory disease in snakes. Controlled studies evaluating environmental stability, disinfection, transmission control, and treatment are lacking.

Update on immunopathology of bornavirus infections in humans and animals

Knowledge on bornaviruses has expanded tremendously during the last decade through detection of novel bornaviruses and endogenous bornavirus-like elements in many eukaryote genomes, as well as by confirmation of insectivores as reservoir species for classical Borna disease virus 1 (BoDV-1). The most intriguing finding was the demonstration of the zoonotic potential of lethal human bornavirus infections caused by a novel bornavirus of different squirrel species (variegated squirrel 1 bornavirus, VSBV-1) and by BoDV-1 known as the causative agent for the classical Borna disease in horses and sheep. Whereas a T cell-mediated immunopathology has already been confirmed as key disease mechanism for infection with BoDV-1 by experimental studies in rodents, the underlying pathomechanisms remain less clear for human bornavirus infections, infection with other bornaviruses or infection of reservoir species. Thus, an overview of current knowledge on the pathogenesis of bornavirus infections focusing on immunopathology is given.

A serine protease-associated lectin in the cytolytic system of blowfly (Chrysomya megacephala) larvae: Evidence and characterization

Cytolytic activity against invading microorganisms is one of the innate forms of immunity in invertebrates. A serine protease-associated sialic acid-specific cytolytic lectin was purified using glutaraldehyde-fixed ox erythrocytes from the larval extract of blowfly (Chrysomya megacephala). The purified lectin lysed vertebrate erythrocytes with effective haemolysis of ox red blood cells (RBCs) in an isotonic medium. The degree of haemolytic (HL) activity of the purified cytolytic lectin depended on its concentration, pH, temperature, and calcium ions. It was sensitive to ethylenediaminetetraacetic acid. The native molecular mass of the C-type lectin was 260 ± 26 kDa, comprising four different polypeptide subunits of 75 kDa (pI ~8), 69 kDa (pI ~7.0), 61 kDa (pI ~5.3), and 55 kDa (pI ~4.6). The association between the C-type lectin and serine protease was confirmed by MALDI-TOF-MS analysis that revealed its homology in the same spectral peak as well as the proteases and phenylmethylsulphonyl fluoride inhibition of HL activity. Haemolysis inhibition by N-acetylneuraminic acid and other sugars revealed the properties of the lectin. The purified lectin distorted the integrity of ox RBCs and Paenalcaligenes hermetiae. This in vitro study documents the presence of a cytolytic system in blowfly (C. megacephala) larvae for the clearance of invading microbial pathogens in their feeding niche.

Investigation of the estuarine stonefish (Synanceia horrida) venom composition

The Estuarine stonefish (Synanceia horrida) is recognised as one of the most venomous fish species in the world but the overall venom composition has yet to be investigated using in-depth transcriptomic and proteomic methods. To date, known venom components are restricted to a hyaluronidase and a large, pore-forming toxin known as Stonustoxin (SNTX). Transcriptomic sequencing of the venom gland resulted in over 170,000 contigs with only 0.4% that were homologous to putative venom proteins. Integration of the transcriptomic data with proteomic data from the S. horrida venom confirmed the hyaluronidase and SNTX to be present, together with several other protein families including major contributions from C-type lectins. Other protein families observed included peroxiredoxin and several minor protein families such as Golgi-associated plant pathogenesis related proteins, tissue pathway factor inhibitors, and Kazal-type serine protease inhibitors that, although not putative venom proteins, may contribute to the venom's adverse effects. BIOLOGICAL SIGNIFICANCE: Proteomic analysis of milked Synanceia horrida venom, paired with transcriptomic analysis of the venom gland tissue revealed for the first time the composition of one of the world's most dangerous fish venoms. The results demonstrate that the venom is relatively less complex compared to other well-studied venomous animals with a number of unique proteins not previously found in animal venoms.

Characterization of the honeybee venom proteins C1q-like protein and PVF1 and their allergenic potential

Honeybee (Apis mellifera) venom (HBV) represents an ideal model to study the role of particular venom components in allergic reactions in sensitized individuals as well as in the eusociality of Hymenoptera species. The aim of this study was to further characterize the HBV components C1q-like protein (C1q) and PDGF/VEGF-like factor 1 (PVF1). C1q and PVF1 were produced as recombinant proteins in insect cells. Their allergenic properties were examined by determining the level of specific IgE antibodies in the sera of HBV-allergic patients (n = 26) as well as by their capacity to activate patients' basophils (n = 11). Moreover, the transcript heterogeneity of PVF1 was analyzed. It could be demonstrated that at least three PVF1 variants are present in the venom gland, which all result from alternative splicing of one transcript. Additionally, recombinant C1q and PVF1 from Spodoptera frugiperda insect cells exhibited specific IgE reactivity with approximately 38.5% of sera of HBV-allergic patients. Interestingly, both proteins were unable to activate basophils of the patients, questioning their role in the context of clinically relevant sensitization. Recombinant C1q and PVF1 can build the basis for a deeper understanding of the molecular mechanisms of Hymenoptera venoms. Moreover, the conflicting results between IgE sensitization and lack of basophil activation, might in the future contribute to the identification of factors that determine the allergenic potential of proteins.

Accelerated evolution of toxin genes: Exonization and intronization in snake venom disintegrin/metalloprotease genes

Toxin genes in animals undergo accelerated evolution compared to non-toxin genes to be effective and competitive in prey capture, as well as to enhance their predator defense. Several mechanisms have been proposed to explain this unusual phenomenon. These include (a) frequent mutations in exons compared to introns and nonsynonymous substitutions in exons; (b) high frequency of point mutations are due to the presence of more unstable triplets in exons compared to introns; (c) Accelerated Segment Switch in Exons to alter Targeting (ASSET); (d) Rapid Accumulation of Variations in Exposed Residues (RAVERs); (e) alteration in intron-exon boundary; (f) deletion of exon; and (g) loss/gain of domains through recombination. By systematic analyses of snake venom disintegrin/metalloprotease genes, I describe a new mechanism in the evolution of these genes through exonization and intronization. In the evolution of RTS/KTS disintegrins, a new exon (10a) is formed in intron 10 of the disintegrin/metalloprotease gene. Unlike more than 90% new exons that are from repetitive elements in introns, exon 10a originated from a non-repetitive element. To incorporate exon 10a, part of the exon 11 is intronized to retain the open reading frame. This is the first case of simultaneous exonization and intronization within a single gene. This new mechanism alters the function of toxins through drastic changes to the molecular surface via insertion of new exons and deletion of exons.

Divergent bornaviruses from Australian carpet pythons with neurological disease date the origin of extant Bornaviridae prior to the end-Cretaceous extinction

Tissue samples from Australian carpet pythons (Morelia spilota) with neurological disease were screened for viruses using next-generation sequencing. Coding complete genomes of two bornaviruses were identified with the gene order 3'-N-X-P-G-M-L, representing a transposition of the G and M genes compared to other bornaviruses and most mononegaviruses. Use of these viruses to search available vertebrate genomes enabled recognition of further endogenous bornavirus-like elements (EBLs) in diverse placental mammals, including humans. Codivergence patterns and shared integration sites revealed an ancestral laurasiatherian EBLG integration (77 million years ago [MYA]) and a previously identified afrotherian EBLG integration (83 MYA). The novel python bornaviruses clustered more closely with these EBLs than with other exogenous bornaviruses, suggesting that these viruses diverged from previously known bornaviruses prior to the end-Cretaceous (K-Pg) extinction, 66 MYA. It is possible that EBLs protected mammals from ancient bornaviral disease, providing a selective advantage in the recovery from the K-Pg extinction. A degenerate PCR primer set was developed to detect a highly conserved region of the bornaviral polymerase gene. It was used to detect 15 more genetically distinct bornaviruses from Australian pythons that represent a group that is likely to contain a number of novel species.

Computational Studies of Snake Venom Toxins

Most snake venom toxins are proteins, and participate to envenomation through a diverse array of bioactivities, such as bleeding, inflammation, and pain, cytotoxic, cardiotoxic or neurotoxic effects. The venom of a single snake species contains hundreds of toxins, and the venoms of the 725 species of venomous snakes represent a large pool of potentially bioactive proteins. Despite considerable discovery efforts, most of the snake venom toxins are still uncharacterized. Modern bioinformatics tools have been recently developed to mine snake venoms, helping focus experimental research on the most potentially interesting toxins. Some computational techniques predict toxin molecular targets, and the binding mode to these targets. This review gives an overview of current knowledge on the ~2200 sequences, and more than 400 three-dimensional structures of snake toxins deposited in public repositories, as well as of molecular modeling studies of the interaction between these toxins and their molecular targets. We also describe how modern bioinformatics have been used to study the snake venom protein phospholipase A2, the small basic myotoxin Crotamine, and the three-finger peptide Mambalgin.

Comparison of proteomic profiles of the venoms of two of the 'Big Four' snakes of India, the Indian cobra (Naja naja) and the common krait (Bungarus caeruleus), and analyses of their toxins

Snake venoms are mixtures of biologically-active proteins and peptides, and several studies have described the characteristics of some of these toxins. However, complete proteomic profiling of the venoms of many snake species has not yet been done. The Indian cobra (Naja naja) and common krait (Bungarus caeruleus) are elapid snake species that are among the 'Big Four' responsible for the majority of human snake envenomation cases in India. As understanding the composition and complexity of venoms is necessary for successful treatment of envenomation in humans, we utilized three different proteomic profiling approaches to characterize these venoms: i) one-dimensional SDS-PAGE coupled with in-gel tryptic digestion and electrospray tandem mass spectrometry (ESI-LC-MS/MS) of individual protein bands; ii) in-solution tryptic digestion of crude venoms coupled with ESI-LC-MS/MS; and iii) separation by gel-filtration chromatography coupled with tryptic digestion and ESI-LC-MS/MS of separated fractions. From the generated data, 81 and 46 different proteins were identified from N. naja and B. caeruleus venoms, respectively, belonging to fifteen different protein families. Venoms from both species were found to contain a variety of phospholipases A₂ and three-finger toxins, whereas relatively higher numbers of snake venom metalloproteinases were found in N. naja compared to B. caeruleus venom. The analyses also identified less represented venom proteins including L-amino acid oxidases, cysteine-rich secretory proteins, 5'-nucleotidases and venom nerve growth factors. Further, Kunitz-type serine protease inhibitors, cobra venom factors, phosphodiesterases, vespryns and aminopeptidases were identified in the N. naja venom, while acetylcholinesterases and hyaluronidases were found in the B. caeruleus venom. We further analyzed protein coverage (Lys/Arg rich and poor regions as well as potential glycosylation sites) using in-house software. These studies expand our understanding of the proteomes of the venoms of these two medically-important species.

Full-Length Venom Protein cDNA Sequences from Venom-Derived mRNA: Exploring Compositional Variation and Adaptive Multigene Evolution

Envenomation of humans by snakes is a complex and continuously evolving medical emergency, and treatment is made that much more difficult by the diverse biochemical composition of many venoms. Venomous snakes and their venoms also provide models for the study of molecular evolutionary processes leading to adaptation and genotype-phenotype relationships. To compare venom complexity and protein sequences, venom gland transcriptomes are assembled, which usually requires the sacrifice of snakes for tissue. However, toxin transcripts are also present in venoms, offering the possibility of obtaining cDNA sequences directly from venom. This study provides evidence that unknown full-length venom protein transcripts can be obtained from the venoms of multiple species from all major venomous snake families. These unknown venom protein cDNAs are obtained by the use of primers designed from conserved signal peptide sequences within each venom protein superfamily. This technique was used to assemble a partial venom gland transcriptome for the Middle American Rattlesnake (Crotalus simus tzabcan) by amplifying sequences for phospholipases A₂, serine proteases, C-lectins, and metalloproteinases from within venom. Phospholipase A₂ sequences were also recovered from the venoms of several rattlesnakes and an elapid snake (Pseudechis porphyriacus), and three-finger toxin sequences were recovered from multiple rear-fanged snake species, demonstrating that the three major clades of advanced snakes (Elapidae, Viperidae, Colubridae) have stable mRNA present in their venoms. These cDNA sequences from venom were then used to explore potential activities derived from protein sequence similarities and evolutionary histories within these large multigene superfamilies. Venom-derived sequences can also be used to aid in characterizing venoms that lack proteomic profiles and identify sequence characteristics indicating specific envenomation profiles. This approach, requiring only venom, provides access to cDNA sequences in the absence of living specimens, even from commercial venom sources, to evaluate important regional differences in venom composition and to study snake venom protein evolution.

This work demonstrates that full-length venom protein messenger RNAs are present in secreted venoms and can be used to acquire full-length protein sequences of toxins from both front-fanged (Elapidae, Viperidae) and rear-fanged (Colubridae) snake venoms, eliminating the need to use venom glands. Full-length transcripts were obtained from venom samples that were fresh, newly lyophilized, old, field desiccated or commercially prepared, representing a significant advance over previous attempts which produced only partial sequence transcripts. Transcripts for all major venom protein families (metalloproteinases, serine proteases, C-type lectins, phospholipases A₂ and three-finger toxins) responsible for clinically significant snakebite symptoms were obtained from venoms. These sequences aid in the identification and characterization of venom proteome profiles, allowing for the identification of peptide sequences, specific isoforms, and novel venom proteins. The application of this technique will help to provide venom protein sequences for many snake species, including understudied rear-fanged snakes. Venom protein transcripts offer important insights into potential snakebite envenomation profiles and the molecular evolution of venom protein multigene families. By requiring only venom to obtain venom protein cDNAs, the approach detailed here will provide access to cDNA-based protein sequences from commercial and other venom sources, facilitating study of snake venom protein composition and evolution.

Molecular Characterization of Three Novel Phospholipase A₂ Proteins from the Venom of Atheris chlorechis, Atheris nitschei and Atheris squamigera

Secretory phospholipase A₂ (sPLA₂) is known as a major component of snake venoms and displays higher-order catalytic hydrolysis functions as well as a wide range of pathological effects. Atheris is not a notoriously dangerous genus of snakes although there are some reports of fatal cases after envenomation due to the effects of coagulation disturbances and hemorrhaging. Molecular characterization of Atheris venom enzymes is incomplete and there are only a few reports in the literature. Here, we report, for the first time, the cloning and characterization of three novel cDNAs encoding phospholipase A₂ precursors (one each) from the venoms of the Western bush viper (Atheris chlorechis), the Great Lakes bush viper (Atheris nitschei) and the Variable bush viper (Atheris squamigera), using a “shotgun cloning” strategy. Open-reading frames of respective cloned cDNAs contained putative 16 residue signal peptides and mature proteins composed of 121 to 123 amino acid residues. Alignment of mature protein sequences revealed high degrees of structural conservation and identity with Group II venom PLA₂ proteins from other taxa within the Viperidae. Reverse-phase High Performance Liquid Chromatography (HPLC) profiles of these three snake venoms were obtained separately and chromatographic fractions were assessed for phospholipase activity using an egg yolk suspension assay. The molecular masses of mature proteins were all identified as approximately 14 kDa. Mass spectrometric analyses of the fractionated oligopeptides arising from tryptic digestion of intact venom proteins, was performed for further structural characterization.

Gene expression profiling of the venom gland from the Venezuelan mapanare (Bothrops colombiensis) using expressed sequence tags (ESTs)

BACKGROUND

Bothrops colombiensis is a highly dangerous pit viper and responsible for over 70% of snakebites in Venezuela. Although the composition in B. colombiensis venom has been identified using a proteome analysis, the venom gland transcriptome is currently lacking.

RESULTS

We constructed a cDNA library from the venom gland of B. colombiensis, and a set of 729 high quality expressed sequence tags (ESTs) was identified. A total number of 344 ESTs (47.2% of total ESTs) was related to toxins. The most abundant toxin transcripts were metalloproteinases (37.5%), phospholipases A2s (PLA2, 29.7%), and serine proteinases (11.9%). Minor toxin transcripts were linked to waprins (5.5%), C-type lectins (4.1%), ATPases (2.9%), cysteine-rich secretory proteins (CRISP, 2.3%), snake venom vascular endothelium growth factors (svVEGF, 2.3%), L-amino acid oxidases (2%), and other putative toxins (1.7%). While 160 ESTs (22% of total ESTs) coded for translation proteins, regulatory proteins, ribosomal proteins, elongation factors, release factors, metabolic proteins, and immune response proteins. Other proteins detected in the transcriptome (87 ESTs, 11.9% of total ESTs) were undescribed proteins with unknown functions. The remaining 138 (18.9%) cDNAs had no match with known GenBank accessions.

CONCLUSION

This study represents the analysis of transcript expressions and provides a physical resource of unique genes for further study of gene function and the development of novel molecules for medical applications.

Antibodies against Venom of the Snake Deinagkistrodon acutus

Snake venom protein from Deinagkistrodon acutus (DA protein), one of the major venomous species in Taiwan, causes hemorrhagic symptoms that can lead to death. Although horse-derived antivenin is a major treatment, relatively strong and detrimental side effects are seen occasionally. In our study, yolk immunoglobulin (IgY) was purified from eggs, and DA protein was recognized using Western blotting and an enzyme-linked immunosorbent assay (ELISA), similar to therapeutic horse antivenin. The ELISA also indicated that specific IgY antibodies were elicited after the fifth booster, plateaued, and lasted for at least 3 months. To generate monoclonal single-chain variable fragment (scFv) antibodies, we used phage display technology to construct two libraries with short or long linkers, containing 6.24 × 10(8) and 5.28 × 10(8) transformants, respectively. After four rounds of biopanning, the eluted phage titer increased, and the phage-based ELISA indicated that the specific clones were enriched. Nucleotide sequences of 30 individual clones expressing scFv were analyzed and classified into four groups that all specifically recognized the DA venom protein. Furthermore, based on mass spectrometry, the scFv-bound protein was deduced to be snake venom metalloproteinase proteins. Most importantly, both IgY and mixed scFv inhibited the lethal effect in mice injected with the minimum lethal dosage of the DA protein. We suggest that together, these antibodies could be applied to the development of diagnostic agents or treatments for snakebite envenomation in the future.

New proline-rich oligopeptides from the venom of African adders: Insights into the hypotensive effect of the venoms

BACKGROUND

The snakes from the Bitis genus are some of the most medically important venomous snakes in sub Saharan Africa, however little is known about the composition and effects of these snake venom peptides. Considering that the victims with Bitis genus snakes have exacerbate hypotension and cardiovascular disorders, we investigated here the presence of angiotensin-converting enzyme modulators on four different species of venoms.

METHODS

The peptide fractions from Bitis gabonica gabonica, Bitis nasicornis, Bitis gabonica rhinoceros and Bitis arietans which showed inhibitory activity on angiotensin-converting enzyme were subjected to mass spectrometry analysis. Eight proline-rich peptides were synthetized and their potencies were evaluated in vitro and in vivo.

RESULTS

The MS analysis resulted in over 150 sequences, out of which 32 are new proline-rich oligopeptides, and eight were selected for syntheses. For some peptides, inhibition assays showed inhibitory potentials of cleavage of angiotensin I ten times greater when compared to bradykinin. In vivo tests showed that all peptides decreased mean arterial pressure, followed by tachycardia in 6 out of 8 of the tests.

CONCLUSION

We describe here some new and already known proline-rich peptides, also known as bradykinin-potentiating peptides. Four synthetic peptides indicated a preferential inhibition of angiotensin-converting enzyme C-domain. In vivo studies show that the proline-rich oligopeptides are hypotensive molecules.

GENERAL SIGNIFICANCE

Although proline-rich oligopeptides are known molecules, we present here 32 new sequences that are inhibitors of the angiotensin-converting enzyme and consistent with the symptoms of the victims of Bitis spp, who display severe hypotension.

VTBuilder: a tool for the assembly of multi isoform transcriptomes

Background

Within many research areas, such as transcriptomics, the millions of short DNA fragments (reads) produced by current sequencing platforms need to be assembled into transcript sequences before they can be utilized. Despite recent advances in assembly software, creating such transcripts from read data harboring isoform variation remains challenging. This is because current approaches fail to identify all variants present or they create chimeric transcripts within which relationships between co-evolving sites and other evolutionary factors are disrupted. We present VTBuilder, a tool for constructing non-chimeric transcripts from read data that has been sequenced from sources containing isoform complexity.

Results

We validated VTBuilder using reads simulated from 54 Sanger sequenced transcripts (SSTs) expressed in the venom gland of the saw scaled viper, Echis ocellatus. The SSTs were selected to represent genes from major co-expressed toxin groups known to harbor isoform variants. From the simulated reads, VTBuilder constructed 55 transcripts, 50 of which had a greater than 99% sequence similarity to 48 of the SSTs. In contrast, using the popular assembler tool Trinity (r2013-02-25), only 14 transcripts were constructed with a similar level of sequence identity to just 11 SSTs. Furthermore VTBuilder produced transcripts with a similar length distribution to the SSTs while those produced by Trinity were considerably shorter. To demonstrate that our approach can be scaled to real world data we assembled the venom gland transcriptome of the African puff adder Bitis arietans using paired-end reads sequenced on Illumina’s MiSeq platform. VTBuilder constructed 1481 transcripts from 5 million reads and, following annotation, all major toxin genes were recovered demonstrating reconstruction of complex underlying sequence and isoform diversity.

Conclusion

Unlike other approaches, VTBuilder strives to maintain the relationships between co-evolving sites within the constructed transcripts, and thus increases transcript utility for a wide range of research areas ranging from transcriptomics to phylogenetics and including the monitoring of drug resistant parasite populations. Additionally, improving the quality of transcripts assembled from read data will have an impact on future studies that query these data. VTBuilder has been implemented in java and is available, under the GPL GPU V0.3 license, from http:// http://www.lstmed.ac.uk/vtbuilder.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0389-8) contains supplementary material, which is available to authorized users.

Exploring the hidden honeybee (Apis mellifera) venom proteome by integrating a combinatorial peptide ligand library approach with FTMS

At present, 30 compounds have been described in the venom of the honeybee, and 16 of them were confirmed by mass spectrometry. Previous studies typically combined 2-D PAGE with MALDI-TOF/TOF MS, a technology which now appears to lack sensitivity to detect additional venom compounds. Here, we report an in-depth study of the honeybee venom proteome using a combinatorial peptide ligand library sample pretreatment to enrich for minor components followed by shotgun LC-FT-ICR MS analysis. This strategy revealed an unexpectedly rich venom composition: in total 102 proteins and peptides were found, with 83 of them never described in bee venom samples before. Based on their predicted function and subcellular location, the proteins could be divided into two groups. A group of 33 putative toxins is proposed to contribute to venom activity by exerting toxic functions or by playing a role in social immunity. The other group, considered as venom trace molecules, appears to be secreted for their functions in the extracellular space, or is unintentionally secreted by the venom gland cells due to insufficient protein recycling or co-secretion with other compounds. In conclusion, our approach allowed to explore the hidden honeybee venom proteome and extended the list of potential venom allergens.

BIOLOGICAL SIGNIFICANCE

This study dug deeper into the complex honeybee venom proteome than ever before by applying a highly performing sample pretreatment and mass spectrometric technology. We present putative biological functions for all identified compounds, largely extending our knowledge of the venom toxicity. In addition, this study offers a long list of potential new venom allergens.

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.

Genome Sequence of a Bornavirus Recovered from an African Garter Snake (Elapsoidea loveridgei)

Bornaviruses are known to infect mammals and birds, and they have been associated with disease in both groups of animals. Here, we report the genome sequence of a bornavirus identified in a wild-caught Loveridge’s garter snake (Elapsoidea loveridgei).

Snake venom L-amino acid oxidases: trends in pharmacology and biochemistry

L-amino acid oxidases are enzymes found in several organisms, including venoms of snakes, where they contribute to the toxicity of ophidian envenomation. Their toxicity is primarily due to enzymatic activity, but other mechanisms have been proposed recently which require further investigation. L-amino acid oxidases exert biological and pharmacological effects, including actions on platelet aggregation and the induction of apoptosis, hemorrhage, and cytotoxicity. These proteins present a high biotechnological potential for the development of antimicrobial, antitumor, and antiprotozoan agents. This review provides an overview of the biochemical properties and pharmacological effects of snake venom L-amino acid oxidases, their structure/activity relationship, and supposed mechanisms of action described so far.

Combined proteomic and transcriptomic interrogation of the venom gland of Conus geographus uncovers novel components and functional compartmentalization

Cone snails are highly successful marine predators that use complex venoms to capture prey. At any given time, hundreds of toxins (conotoxins) are synthesized in the secretory epithelial cells of the venom gland, a long and convoluted organ that can measure 4 times the length of the snail's body. In recent years a number of studies have begun to unveil the transcriptomic, proteomic and peptidomic complexity of the venom and venom glands of a number of cone snail species. By using a combination of DIGE, bottom-up proteomics and next-generation transcriptome sequencing the present study identifies proteins involved in envenomation and conotoxin maturation, significantly extending the repertoire of known (poly)peptides expressed in the venom gland of these remarkable animals. We interrogate the molecular and proteomic composition of different sections of the venom glands of 3 specimens of the fish hunter Conus geographus and demonstrate regional variations in gene expression and protein abundance. DIGE analysis identified 1204 gel spots of which 157 showed significant regional differences in abundance as determined by biological variation analysis. Proteomic interrogation identified 342 unique proteins including those that exhibited greatest fold change. The majority of these proteins also exhibited significant changes in their mRNA expression levels validating the reliability of the experimental approach. Transcriptome sequencing further revealed a yet unknown genetic diversity of several venom gland components. Interestingly, abundant proteins that potentially form part of the injected venom mixture, such as echotoxins, phospholipase A2 and con-ikots-ikots, classified into distinct expression clusters with expression peaking in different parts of the gland. Our findings significantly enhance the known repertoire of venom gland polypeptides and provide molecular and biochemical evidence for the compartmentalization of this organ into distinct functional entities.

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 genesis of an exceptionally lethal venom in the timber rattlesnake (Crotalus horridus) revealed through comparative venom-gland transcriptomics

Background

Snake venoms generally show sequence and quantitative variation within and between species, but some rattlesnakes have undergone exceptionally rapid, dramatic shifts in the composition, lethality, and pharmacological effects of their venoms. Such shifts have occurred within species, most notably in Mojave (Crotalus scutulatus), South American (C. durissus), and timber (C. horridus) rattlesnakes, resulting in some populations with extremely potent, neurotoxic venoms without the hemorrhagic effects typical of rattlesnake bites.

Results

To better understand the evolutionary changes that resulted in the potent venom of a population of C. horridus from northern Florida, we sequenced the venom-gland transcriptome of an animal from this population for comparison with the previously described transcriptome of the eastern diamondback rattlesnake (C. adamanteus), a congener with a more typical rattlesnake venom. Relative to the toxin transcription of C. adamanteus, which consisted primarily of snake-venom metalloproteinases, C-type lectins, snake-venom serine proteinases, and myotoxin-A, the toxin transcription of C. horridus was far simpler in composition and consisted almost entirely of snake-venom serine proteinases, phospholipases A₂, and bradykinin-potentiating and C-type natriuretic peptides. Crotalus horridus lacked significant expression of the hemorrhagic snake-venom metalloproteinases and C-type lectins. Evolution of shared toxin families involved differential expansion and loss of toxin clades within each species and pronounced differences in the highly expressed toxin paralogs. Toxin genes showed significantly higher rates of nonsynonymous substitution than nontoxin genes. The expression patterns of nontoxin genes were conserved between species, despite the vast differences in toxin expression.

Conclusions

Our results represent the first complete, sequence-based comparison between the venoms of closely related snake species and reveal in unprecedented detail the rapid evolution of snake venoms. We found that the difference in venom properties resulted from major changes in expression levels of toxin gene families, differential gene-family expansion and loss, changes in which paralogs within gene families were expressed at high levels, and higher nonsynonymous substitution rates in the toxin genes relative to nontoxins. These massive alterations in the genetics of the venom phenotype emphasize the evolutionary lability and flexibility of this ecologically critical trait.

Evolutionarily conserved interaction between the phosphoproteins and X proteins of bornaviruses from different vertebrate species

Bornavirus, a non-segmented, negative-strand RNA viruses, is currently classified into several genetically distinct genotypes, such as Borna disease virus (BDV) and avian bornaviruses (ABVs). Recent studies revealed that bornavirus genotypes show unique sequence variability in the putative 5′ untranslated region (5′ UTR) of X/P mRNA, a bicistronic mRNA for the X protein and phosphoprotein (P). In this study, to understand the evolutionary relationship among the bornavirus genotypes, we investigated the functional interaction between the X and P proteins of four bornavirus genotypes, BDV, ABV genotype 4 and 5 and reptile bornavirus (RBV), the putative 5′ UTRs of which exhibit variation in the length. Immunofluorescence and immunoprecipitation analyses using mammalian and avian cell lines revealed that the X proteins of bornaviruses conserve the ability to facilitate the export of P from the nucleus to the cytoplasm via interaction with P. Furthermore, we showed that inter-genotypic interactions may occur between X and P among the genotypes, except for X of RBV. In addition, a BDV minireplicon assay demonstrated that the X and P proteins of ABVs, but not RBV, can affect the polymerase activity of BDV. This study demonstrates that bornaviruses may have conserved the fundamental function of a regulatory protein during their evolution, whereas RBV has evolved distinctly from the other bornavirus genotypes.

Biochemical and biological characterization of two serine proteinases from Colombian Crotalus durissus cumanensis snake venom

Two clotting serine proteinases, named Cdc SI and Cdc SII, were isolated and characterized for the first time from Colombian Crotalus durissus cumanensis snake venom. The enzymes were purified using two chromatographic steps: molecular exclusion on Sephacryl S-200 and RP-HPLC on C8 Column. The molecular masses of the proteins, determined by MALDI-TOF mass spectrometry, were 28,561.4 and 28,799.2 Da for Cdc SI and Cdc SII, respectively. The aim of the present study was to evaluate enzymatic, coagulant and toxic properties of the two enzymes. The serine proteinases hydrolyzed specific chromogenic substrate (BaPNA) and exhibited a Michaelis-Menten behavior. Cdc SI had V(max) of 0.038 ± 0.003 nmol/min and K(M) of 0.034 ± 0.017 mM, while Cdc SII displayed values of V(max) of 0.267 ± 0.011 nmol/min and K(M) of 0.145 ± 0.023 mM. N-terminal sequences were VIGGDEXNIN and VIGGDICNINEHNFLVALYE for Cdc SI and Cdc SII, respectively. Molecular masses, N-terminal sequences, inhibition assays, and enzymatic profile suggest that Cdc SI and Cdc SII belong to the family of snake venom thrombin-like enzymes. These serine proteinases differed in their clotting activity on human plasma, showing a minimum coagulant dose of 25 μg and 0.571 μg for Cdc SI and Cdc SII, respectively. Enzymes also showed coagulant activity on bovine fibrinogen and degraded chain α of this protein. Toxins lack hemorrhagic and myotoxic activities, but are capable to induce defibrin(ogen)ation, moderate edema, and an increase in vascular permeability. These serine proteinases may contribute indirectly to the local hemorrhage induced by metalloproteinases, by causing blood clotting disturbances, and might also contribute to cardiovascular alterations characteristic of patients envenomed by C. d. cumanensis in Colombia.

Sequence and phylogenetic analysis of viper venom serine proteases

Snakebites are a major neglected tropical disease responsible for as many as 95000 deaths every year worldwide. Viper venom serine proteases disrupt haemostasis of prey and victims by affecting various stages of the blood coagulation system. A better understanding of their sequence, structure, function and phylogenetic relationships will improve the knowledge on the pathological conditions and aid in the development of novel therapeutics for treating snakebites. A large dataset for all available viper venom serine proteases was developed and analysed to study various features of these enzymes. Despite the large number of venom serine protease sequences available, only a small proportion of these have been functionally characterised. Although, they share some of the common features such as a C-terminal extension, GWG motif and disulphide linkages, they vary widely between each other in features such as isoelectric points, potential N-glycosylation sites and functional characteristics. Some of the serine proteases contain substitutions for one or more of the critical residues in catalytic triad or primary specificity pockets. Phylogenetic analysis clustered all the sequences in three major groups. The sequences with substitutions in catalytic triad or specificity pocket clustered together in separate groups. Our study provides the most complete information on viper venom serine proteases to date and improves the current knowledge on the sequence, structure, function and phylogenetic relationships of these enzymes. This collective analysis of venom serine proteases will help in understanding the complexity of envenomation and potential therapeutic avenues.

Rhinocetin, a venom-derived integrin-specific antagonist inhibits collagen-induced platelet and endothelial cell functions

Snaclecs are small non-enzymatic proteins present in viper venoms reported to modulate hemostasis of victims through effects on platelets, vascular endothelial, and smooth muscle cells. In this study, we have isolated and functionally characterized a snaclec that we named "rhinocetin" from the venom of West African gaboon viper, Bitis gabonica rhinoceros. Rhinocetin was shown to comprise α and β chains with the molecular masses of 13.5 and 13 kDa, respectively. Sequence and immunoblot analysis of rhinocetin confirmed this to be a novel snaclec. Rhinocetin inhibited collagen-stimulated activation of human platelets in a dose-dependent manner but displayed no inhibitory effects on glycoprotein VI (collagen receptor) selective agonist, CRP-XL-, ADP-, or thrombin-induced platelet activation. Rhinocetin antagonized the binding of monoclonal antibodies against the α2 subunit of integrin α2β1 to platelets and coimmunoprecipitation analysis confirmed integrin α2β1 as a target for this venom protein. Rhinocetin inhibited a range of collagen-induced platelet functions such as fibrinogen binding, calcium mobilization, granule secretion, aggregation, and thrombus formation. It also inhibited integrin α2β1-dependent functions of human endothelial cells. Together, our data suggest rhinocetin to be a modulator of integrin α2β1 function and thus may provide valuable insights into the role of this integrin in physiological and pathophysiological scenarios, including hemostasis, thrombosis, and envenomation.

Venom peptide analysis of Vipera ammodytes meridionalis (Viperinae) and Bothrops jararacussu (Crotalinae) demonstrates subfamily-specificity of the peptidome in the family Viperidae

Snake venom peptidomes are valuable sources of pharmacologically active compounds. We analyzed the peptidic fractions (peptides with molecular masses < 10,000 Da) of venoms of Vipera ammodytes meridionalis (Viperinae), the most toxic snake in Europe, and Bothrops jararacussu (Crotalinae), an extremely poisonous snake of South America. Liquid chromatography/mass spectrometry (LC/MS), direct infusion electrospray mass spectrometry (ESI-MS) and matrix-assisted desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS) were applied to characterize the peptides of both snake venoms. 32 bradykinin-potentiating peptides (BPPs) were identified in the Crotalinae venom and their sequences determined. 3 metalloproteinase inhibitors, 10 BPPs and a Kunitz-type inhibitor were observed in the Viperinae venom peptidome. Variability in the C-terminus of homologous BPPs was observed, which can influence the pharmacological effects. The data obtained so far show a subfamily specificity of the venom peptidome in the Viperidae family: BPPs are the major peptide component of the Crotalinae venom peptidome lacking Kunitz-type inhibitors (with one exception) while the Viperinae venom, in addition to BPPs, can contain peptides of the bovine pancreatic trypsin inhibitor family. We found indications for a post-translational phosphorylation of serine residues in Bothrops jararacussu venom BPP (S[combining low line]QGLPPGPPIP), which could be a regulatory mechanism in their interactions with ACE, and might influence the hypotensive effect. Homology between venom BPPs from Viperidae snakes and venom natriuretic peptide precursors from Elapidae snakes suggests a structural similarity between the respective peptides from the peptidomes of both snake families. The results demonstrate that the venoms of both snakes are rich sources of peptides influencing important physiological systems such as blood pressure regulation and hemostasis. The data can be used for pharmacological and medical applications.

Evolutionary analysis of novel serine proteases in the venom gland transcriptome of Bitis gabonica rhinoceros

BACKGROUND

Serine proteases are major components of viper venom and target various stages of the blood coagulation system in victims and prey. A better understanding of the diversity of serine proteases and other enzymes present in snake venom will help to understand how the complexity of snake venom has evolved and will aid the development of novel therapeutics for treating snake bites.

METHODOLOGY AND PRINCIPAL FINDINGS

Four serine protease-encoding genes from the venom gland transcriptome of Bitis gabonica rhinoceros were amplified and sequenced. Mass spectrometry suggests the four enzymes corresponding to these genes are present in the venom of B. g. rhinoceros. Two of the enzymes, rhinocerases 2 and 3 have substitutions to two of the serine protease catalytic triad residues and are thus unlikely to be catalytically active, though they may have evolved other toxic functions. The other two enzymes, rhinocerases 4 and 5, have classical serine protease catalytic triad residues and thus are likely to be catalytically active, however they have glycine rather than the more typical aspartic acid at the base of the primary specificity pocket (position 189). Based on a detailed analysis of these sequences we suggest that alternative splicing together with individual amino acid mutations may have been involved in their evolution. Changes within amino acid segments which were previously proposed to undergo accelerated change in venom serine proteases have also been observed.

CONCLUSIONS AND SIGNIFICANCE

Our study provides further insight into the diversity of serine protease isoforms present within snake venom and discusses their possible functions and how they may have evolved. These multiple serine protease isoforms with different substrate specificities may enhance the envenomation effects and help the snake to adapt to new habitats and diets. Our findings have potential for helping the future development of improved therapeutics for snake bites.

A high-throughput venom-gland transcriptome for the Eastern Diamondback Rattlesnake (Crotalus adamanteus) and evidence for pervasive positive selection across toxin classes

Despite causing considerable human mortality and morbidity, animal toxins represent a valuable source of pharmacologically active macromolecules, a unique system for studying molecular adaptation, and a powerful framework for examining structure-function relationships in proteins. Snake venoms are particularly useful in the latter regard as they consist primarily of a moderate number of proteins and peptides that have been found to belong to just a handful of protein families. As these proteins and peptides are produced in dedicated glands, transcriptome sequencing has proven to be an effective approach to identifying the expressed toxin genes. We generated a venom-gland transcriptome for the Eastern Diamondback Rattlesnake (Crotalus adamanteus) using Roche 454 sequencing technology. In the current work, we focus on transcripts encoding toxins. We identified 40 unique toxin transcripts, 30 of which have full-length coding sequences, and 10 have only partial coding sequences. These toxins account for 24% of the total sequencing reads. We found toxins from 11 previously described families of snake-venom toxins and have discovered two putative, previously undescribed toxin classes. The most diverse and highly expressed toxin classes in the C. adamanteus venom-gland transcriptome are the serine proteinases, metalloproteinases, and C-type lectins. The serine proteinases are the most abundant class, accounting for 35% of the toxin sequencing reads. Metalloproteinases are the most diverse; 11 different forms have been identified. Using our sequences and those available in public databases, we detected positive selection in seven of the eight toxin families for which sufficient sequences were available for the analysis. We find that the vast majority of the genes that contribute directly to this vertebrate trait show evidence for a role for positive selection in their evolutionary history.

Cloning and characterisation of novel cystatins from elapid snake venom glands

Snake venoms contain a complex mixture of polypeptides that modulate prey homeostatic mechanisms through highly specific and targeted interactions. In this study we have identified and characterised cystatin-like cysteine-protease inhibitors from elapid snake venoms for the first time. Novel cystatin sequences were cloned from 12 of 13 elapid snake venom glands and the protein was detected, albeit at very low levels, in a total of 22 venoms. One highly conserved isoform, which displayed close sequence identity with family 2 cystatins, was detected in each elapid snake. Crude Austrelaps superbus (Australian lowland copperhead) snake venom inhibited papain, and a recombinant form of A. superbus cystatin inhibited cathepsin L ≅ papain > cathepsin B, with no inhibition observed for calpain or legumain. While snake venom cystatins have truncated N-termini, sequence alignment and structural modelling suggested that the evolutionarily conserved Gly-11 of family 2 cystatins, essential for cysteine protease inhibition, is conserved in snake venom cystatins as Gly-3. This was confirmed by mutagenesis at the Gly-3 site, which increased the dissociation constant for papain by 10(4)-fold. These data demonstrate that elapid snake venom cystatins are novel members of the type 2 family. The widespread, low level expression of type 2 cystatins in snake venom, as well as the presence of only one highly conserved isoform in each species, imply essential housekeeping or regulatory roles for these proteins.

A transcriptomic analysis of gene expression in the venom gland of the snake Bothrops alternatus (urutu)

Background

The genus Bothrops is widespread throughout Central and South America and is the principal cause of snakebite in these regions. Transcriptomic and proteomic studies have examined the venom composition of several species in this genus, but many others remain to be studied. In this work, we used a transcriptomic approach to examine the venom gland genes of Bothrops alternatus, a clinically important species found in southeastern and southern Brazil, Uruguay, northern Argentina and eastern Paraguay.

Results

A cDNA library of 5,350 expressed sequence tags (ESTs) was produced and assembled into 838 contigs and 4512 singletons. BLAST searches of relevant databases showed 30% hits and 70% no-hits, with toxin-related transcripts accounting for 23% and 78% of the total transcripts and hits, respectively. Gene ontology analysis identified non-toxin genes related to general metabolism, transcription and translation, processing and sorting, (polypeptide) degradation, structural functions and cell regulation. The major groups of toxin transcripts identified were metalloproteinases (81%), bradykinin-potentiating peptides/C-type natriuretic peptides (8.8%), phospholipases A₂ (5.6%), serine proteinases (1.9%) and C-type lectins (1.5%). Metalloproteinases were almost exclusively type PIII proteins, with few type PII and no type PI proteins. Phospholipases A₂ were essentially acidic; no basic PLA₂ were detected. Minor toxin transcripts were related to L-amino acid oxidase, cysteine-rich secretory proteins, dipeptidylpeptidase IV, hyaluronidase, three-finger toxins and ohanin. Two non-toxic proteins, thioredoxin and double-specificity phosphatase Dusp6, showed high sequence identity to similar proteins from other snakes. In addition to the above features, single-nucleotide polymorphisms, microsatellites, transposable elements and inverted repeats that could contribute to toxin diversity were observed.

Conclusions

Bothrops alternatus venom gland contains the major toxin classes described for other Bothrops venoms based on trancriptomic and proteomic studies. The predominance of type PIII metalloproteinases agrees with the well-known hemorrhagic activity of this venom, whereas the lower content of serine proteases and C-type lectins could contribute to less marked coagulopathy following envenoming by this species. The lack of basic PLA₂ agrees with the lower myotoxicity of this venom compared to other Bothrops species with these toxins. Together, these results contribute to our understanding of the physiopathology of envenoming by this species.

Serine protease variants encoded by Echis ocellatus venom gland cDNA: cloning and sequencing analysis

Envenoming by Echis saw-scaled viper is the leading cause of death and morbidity in Africa due to snake bite. Despite its medical importance, there have been few investigations into the toxin composition of the venom of this viper. Here, we report the cloning of cDNA sequences encoding four groups or isoforms of the haemostasis-disruptive Serine protease proteins (SPs) from the venom glands of Echis ocellatus. All these SP sequences encoded the cysteine residues scaffold that form the 6-disulphide bonds responsible for the characteristic tertiary structure of venom serine proteases. All the Echis ocellatus EoSP groups showed varying degrees of sequence similarity to published viper venom SPs. However, these groups also showed marked intercluster sequence conservation across them which were significantly different from that of previously published viper SPs. Because viper venom SPs exhibit a high degree of sequence similarity and yet exert profoundly different effects on the mammalian haemostatic system, no attempt was made to assign functionality to the new Echis ocellatus EoSPs on the basis of sequence alone. The extraordinary level of interspecific and intergeneric sequence conservation exhibited by the Echis ocellatus EoSPs and analogous serine proteases from other viper species leads us to speculate that antibodies to representative molecules should neutralise (that we will exploit, by epidermal DNA immunization) the biological function of this important group of venom toxins in vipers that are distributed throughout Africa, the Middle East, and the Indian subcontinent.

Purification and functional characterisation of rhiminopeptidase A, a novel aminopeptidase from the venom of Bitis gabonica rhinoceros

Background

Snake bite is a major neglected public health issue within poor communities living in the rural areas of several countries throughout the world. An estimated 2.5 million people are bitten by snakes each year and the cost and lack of efficacy of current anti-venom therapy, together with the lack of detailed knowledge about toxic components of venom and their modes of action, and the unavailability of treatments in rural areas mean that annually there are around 125,000 deaths worldwide. In order to develop cheaper and more effective therapeutics, the toxic components of snake venom and their modes of action need to be clearly understood. One particularly poorly understood component of snake venom is aminopeptidases. These are exo-metalloproteases, which, in mammals, are involved in important physiological functions such as the maintenance of blood pressure and brain function. Although aminopeptidase activities have been reported in some snake venoms, no detailed analysis of any individual snake venom aminopeptidases has been performed so far. As is the case for mammals, snake venom aminopeptidases may also play important roles in altering the physiological functions of victims during envenomation. In order to further understand this important group of snake venom enzymes we have isolated, functionally characterised and analysed the sequence-structure relationships of an aminopeptidase from the venom of the large, highly venomous West African gaboon viper, Bitis gabonica rhinoceros.

Methodology and Principal Findings

The venom of B. g. rhinoceros was fractionated by size exclusion chromatography and fractions with aminopeptidase activities were isolated. Fractions with aminopeptidase activities showed a pure protein with a molecular weight of 150 kDa on SDS-PAGE. In the absence of calcium, this purified protein had broad aminopeptidase activities against acidic, basic and neutral amino acids but in the presence of calcium, it had only acidic aminopeptidase activity (APA). Together with the functional data, mass spectrometry analysis of the purified protein confirmed this as an aminopeptidase A and thus this has been named as rhiminopeptidase A. The complete gene sequence of rhiminopeptidase A was obtained by sequencing the PCR amplified aminopeptidase A gene from the venom gland cDNA of B. g. rhinoceros. The gene codes for a predicted protein of 955 amino acids (110 kDa), which contains the key amino acids necessary for functioning as an aminopeptidase A. A structural model of rhiminopeptidase A shows the structure to consist of 4 domains: an N-terminal saddle-shaped β domain, a mixed α and β catalytic domain, a β-sandwich domain and a C-terminal α helical domain.

Conclusions

This study describes the discovery and characterisation of a novel aminopeptidase A from the venom of B. g. rhinoceros and highlights its potential biological importance. Similar to mammalian aminopeptidases, rhiminopeptidase A might be capable of playing roles in altering the blood pressure and brain function of victims. Furthermore, it could have additional effects on the biological functions of other host proteins by cleaving their N-terminal amino acids. This study points towards the importance of complete analysis of individual components of snake venom in order to develop effective therapies for snake bites.

Snake bite is a major neglected public health issue causing an estimated 125,000 deaths each year, predominantly within poor communities living in rural areas of countries in South East Asia and Africa. Current treatments for snake bites are costly and have limited effectiveness, thus there is a need to develop novel therapeutics. In order to do this the toxic components of snake venom need to be clearly understood. Enzymes called aminopeptidases have been noticed in several snake venoms, but their functions have not been characterised. Related enzymes are also present in mammals, where they are involved in the maintenance of blood pressure and brain function. To further understand this important group of enzymes within snake venom we have purified and analysed the function and structure of an aminopeptidase from the venom of the West African gaboon viper. Our results suggest that this enzyme could also affect the maintenance of blood pressure and brain function in victims of snake bites. Along with other snake venom components, aminopeptidases might be a potential therapeutic target for developing novel treatments for snake bites.

Transcriptomic analysis of the venom gland of the red-headed krait (Bungarus flaviceps) using expressed sequence tags

Background

The Red-headed krait (Bungarus flaviceps, Squamata: Serpentes: Elapidae) is a medically important venomous snake that inhabits South-East Asia. Although the venoms of most species of the snake genus Bungarus have been well characterized, a detailed compositional analysis of B. flaviceps is currently lacking.

Results

Here, we have sequenced 845 expressed sequence tags (ESTs) from the venom gland of a B. flaviceps. Of the transcripts, 74.8% were putative toxins; 20.6% were cellular; and 4.6% were unknown. The main venom protein families identified were three-finger toxins (3FTxs), Kunitz-type serine protease inhibitors (including chain B of β-bungarotoxin), phospholipase A₂ (including chain A of β-bungarotoxin), natriuretic peptide (NP), CRISPs, and C-type lectin.

Conclusion

The 3FTxs were found to be the major component of the venom (39%). We found eight groups of unique 3FTxs and most of them were different from the well-characterized 3FTxs. We found three groups of Kunitz-type serine protease inhibitors (SPIs); one group was comparable to the classical SPIs and the other two groups to chain B of β-bungarotoxins (with or without the extra cysteine) based on sequence identity. The latter group may be functional equivalents of dendrotoxins in Bungarus venoms. The natriuretic peptide (NP) found is the first NP for any Asian elapid, and distantly related to Australian elapid NPs. Our study identifies several unique toxins in B. flaviceps venom, which may help in understanding the evolution of venom toxins and the pathophysiological symptoms induced after envenomation.

Purification and functional characterisation of rhinocerase, a novel serine protease from the venom of Bitis gabonica rhinoceros

Background

Serine proteases are a major component of viper venoms and are thought to disrupt several distinct elements of the blood coagulation system of envenomed victims. A detailed understanding of the functions of these enzymes is important both for acquiring a fuller understanding of the pathology of envenoming and because these venom proteins have shown potential in treating blood coagulation disorders.

Methodology/Principal Findings

In this study a novel, highly abundant serine protease, which we have named rhinocerase, has been isolated and characterised from the venom of Bitis gabonica rhinoceros using liquid phase isoelectric focusing and gel filtration. Like many viper venom serine proteases, this enzyme is glycosylated; the estimated molecular mass of the native enzyme is approximately 36kDa, which reduces to 31kDa after deglycosylation. The partial amino acid sequence shows similarity to other viper venom serine proteases, but is clearly distinct from the sequence of the only other sequenced serine protease from Bitis gabonica. Other viper venom serine proteases have been shown to exert distinct biological effects, and our preliminary functional characterization of rhinocerase suggest it to be multifunctional. It is capable of degrading α and β chains of fibrinogen, dissolving plasma clots and of hydrolysing a kallikrein substrate.

Conclusions/Significance

A novel multifunctional viper venom serine protease has been isolated and characterised. The activities of the enzyme are consistent with the known in vivo effects of Bitis gabonica envenoming, including bleeding disorders, clotting disorders and hypotension. This study will form the basis for future research to understand the mechanisms of serine protease action, and examine the potential for rhinocerase to be used clinically to reduce the risk of human haemostatic disorders such as heart attacks and strokes.

Endogenous non-retroviral RNA virus elements in mammalian genomes

Retroviruses are the only group of viruses known to have left a fossil record, in the form of endogenous proviruses, and approximately 8% of the human genome is made up of these elements. Although many other viruses, including non-retroviral RNA viruses, are known to generate DNA forms of their own genomes during replication, none has been found as DNA in the germline of animals. Bornaviruses, a genus of non-segmented, negative-sense RNA virus, are unique among RNA viruses in that they establish persistent infection in the cell nucleus. Here we show that elements homologous to the nucleoprotein (N) gene of bornavirus exist in the genomes of several mammalian species, including humans, non-human primates, rodents and elephants. These sequences have been designated endogenous Borna-like N (EBLN) elements. Some of the primate EBLNs contain an intact open reading frame (ORF) and are expressed as mRNA. Phylogenetic analyses showed that EBLNs seem to have been generated by different insertional events in each specific animal family. Furthermore, the EBLN of a ground squirrel was formed by a recent integration event, whereas those in primates must have been formed more than 40 million years ago. We also show that the N mRNA of a current mammalian bornavirus, Borna disease virus (BDV), can form EBLN-like elements in the genomes of persistently infected cultured cells. Our results provide the first evidence for endogenization of non-retroviral virus-derived elements in mammalian genomes and give novel insights not only into generation of endogenous elements, but also into a role of bornavirus as a source of genetic novelty in its host.

Venomics of the spider Ornithoctonus huwena based on transcriptomic versus proteomic analysis

The spider Ornithoctonus huwena is a venomous spider found in southern China. Its venom is a complex mixture of numerous biologically active components. In this study, 41 novel unique transcripts encoding cellular proteins or other possible venom components were generated from the previously constructed cDNA library. These proteins were also annotated by KOG (eukaryotic orthologous group) and GO (gene ontology) terms. A novel cellular transcript contig encoding an EF-hand protein (named HWEFHP1) was found, which might be involved in the secretion of toxins in the venom glands. In order to have an overview of the molecular diversity of the O. huwena venom, the datasets of all the transcripts, peptides and proteins known so far were analyzed. A comparison of the data obtained through a proteomic versus a transcriptomic approach, revealed that only 15 putative cystine knot toxins (CKTs) were identified by both approaches, 29 transcripts coding for CKTs were found in the transcriptome but not as translated peptides in the venom proteome. However, no cellular protein with identical molecular weight was identified by both approaches. Our data may contribute to a deeper understanding of the biology and ecology of O. huwena and the relationship between structure and function of individual toxins.