Snake venom dipeptidyl peptidase IV: taxonomic distribution and quantitative variation

Proteomic Identification and Quantification of Snake Venom Biomarkers in Venom and Plasma Extracellular Vesicles

The global exploration of snakebites requires the use of quantitative omics approaches to characterize snake venom as it enters into the systemic circulation. These omics approaches give insights into the venom proteome, but a further exploration is warranted to analyze the venom-reactome for the identification of snake venom biomarkers. The recent discovery of extracellular vesicles (EVs), and their critical cellular functions, has presented them as intriguing sources for biomarker discovery and disease diagnosis. Herein, we purified EV’s from the snake venom (svEVs) of Crotalus atrox and C. oreganus helleri, and from plasma of BALB/c mice injected with venom from each snake using EVtrap in conjunction with quantitative mass spectrometry for the proteomic identification and quantification of svEVs and plasma biomarkers. Snake venom EVs from C. atrox and C. o. helleri were highly enriched in 5′ nucleosidase, L-amino acid oxidase, and metalloproteinases. In mouse plasma EVs, a bioinformatic analysis for revealed upregulated responses involved with cytochrome P450, lipid metabolism, acute phase inflammation immune, and heat shock responses, while downregulated proteins were associated with mitochondrial electron transport, NADH, TCA, cortical cytoskeleton, reticulum stress, and oxidative reduction. Altogether, this analysis will provide direct evidence for svEVs composition and observation of the physiological changes of an envenomated organism.

An integrated transcriptomic and proteomic approach to identify the main Torymus sinensis venom components

During oviposition, ectoparasitoid wasps not only inject their eggs but also a complex mixture of proteins and peptides (venom) in order to regulate the host physiology to benefit their progeny. Although several endoparasitoid venom proteins have been identified, little is known about the components of ectoparasitoid venom. To characterize the protein composition of Torymus sinensis Kamijo (Hymenoptera: Torymidae) venom, we used an integrated transcriptomic and proteomic approach and identified 143 venom proteins. Moreover, focusing on venom gland transcriptome, we selected additional 52 transcripts encoding putative venom proteins. As in other parasitoid venoms, hydrolases, including proteases, phosphatases, esterases, and nucleases, constitute the most abundant families in T. sinensis venom, followed by protease inhibitors. These proteins are potentially involved in the complex parasitic syndrome, with different effects on the immune system, physiological processes and development of the host, and contribute to provide nutrients to the parasitoid progeny. Although additional in vivo studies are needed, initial findings offer important information about venom factors and their putative host effects, which are essential to ensure the success of parasitism.

A Wolf in Another Wolf's Clothing: Post-Genomic Regulation Dictates Venom Profiles of Medically-Important Cryptic Kraits in India

The Common Krait (Bungarus caeruleus) shares a distribution range with many other ‘phenotypically-similar’ kraits across the Indian subcontinent. Despite several reports of fatal envenomings by other Bungarus species, commercial Indian antivenoms are only manufactured against B. caeruleus. It is, therefore, imperative to understand the distribution of genetically distinct lineages of kraits, the compositional differences in their venoms, and the consequent impact of venom variation on the (pre)clinical effectiveness of antivenom therapy. To address this knowledge gap, we conducted phylogenetic and comparative venomics investigations of kraits in Southern and Western India. Phylogenetic reconstructions using mitochondrial markers revealed a new species of krait, Romulus’ krait (Bungarus romulusi sp. nov.), in Southern India. Additionally, we found that kraits with 17 mid-body dorsal scale rows in Western India do not represent a subspecies of the Sind Krait (B. sindanus walli) as previously believed, but are genetically very similar to B. sindanus in Pakistan. Furthermore, venom proteomics and comparative transcriptomics revealed completely contrasting venom profiles. While the venom gland transcriptomes of all three species were highly similar, venom proteomes and toxicity profiles differed significantly, suggesting the prominent role of post-genomic regulatory mechanisms in shaping the venoms of these cryptic kraits. In vitro venom recognition and in vivo neutralisation experiments revealed a strong negative impact of venom variability on the preclinical performance of commercial antivenoms. While the venom of B. caeruleus was neutralised as per the manufacturer’s claim, performance against the venoms of B. sindanus and B. romulusi was poor, highlighting the need for regionally-effective antivenoms in India.

Venom composition of adult Western Diamondback Rattlesnakes (Crotalus atrox) maintained under controlled diet and environmental conditions shows only minor changes

Many species of snakes produce venom as a chemical means of procuring potentially fractious prey. Studies have increasingly focused on venom compositional variation between and within individual snakes of the same species/subspecies, with significant differences often being observed. This variation in composition has been attributed to differences in age, season, diet, and environment, suggesting that these factors could help explain the inter- and intra-specific variation found in some snake venoms, perhaps via some type of feedback mechanism(s). To address several of these possible sources of variation, this study utilized wild-caught Western Diamondback Rattlesnakes (Crotalus atrox) from Cochise Co., AZ. Sixteen adult C. atrox were maintained in the lab on a diet of NSA mice for eight months to determine whether venom composition changed in captivity under a static diet in a stable environment. Reducing 1-D SDS-PAGE, fibrinogen degradation assays, reversed-phase HPLC, and MALDI-TOF mass spectrometry revealed only minor differences over time within individuals. Venom L-amino acid oxidase (LAAO) and phosphodiesterase activities significantly increased over the course of captivity, with no changes occurring in azocasein metalloproteinase, kallikrein-like serine proteinase (KLSP), or thrombin-like serine proteinase (TLSP) activities. Snake total length was positively correlated with TLSP activity and negatively correlated with LAAO and KLSP activity. There was typically a much higher degree of variation between individuals than within individuals for all analyses performed and measurements collected. Because the overall "fingerprint" of each snake's venom remained more/less constant, it is concluded that biologically significant changes in venom composition did not occur within individual C. atrox as a function of captivity/diet. However, this study does indicate that differences in activity levels do occur in minor venom enzyme components, but the differences observed are likely to be of minimal significance to the production of antivenom or to subsequent treatment of human envenomations.

In-Depth Venome of the Brazilian Rattlesnake Crotalus durissus terrificus: An Integrative Approach Combining Its Venom Gland Transcriptome and Venom Proteome

Snake venoms are complex mixtures mainly composed of proteins and small peptides. Crotoxin is one of the most studied components from Crotalus venoms, but many other components are less known due to their low abundance. The venome of Crotalus durissus terrificus, the most lethal Brazilian snake, was investigated by combining its venom gland transcriptome and proteome to create a holistic database of venom compounds unraveling novel toxins. We constructed a cDNA library from C. d. terrificus venom gland using the Illumina platform and investigated its venom proteome through high resolution liquid chromotography-tandem mass spectrometry. After integrating data from both data sets, more than 30 venom components classes were identified by the transcriptomic analysis and 15 of them were detected in the venom proteome. However, few of them (PLA₂, SVMP, SVSP, and VEGF) were relatively abundant. Furthermore, only seven expressed transcripts contributed to ∼82% and ∼73% of the abundance in the transcriptome and proteome, respectively. Additionally, novel venom proteins are reported, and we highlight the importance of using different databases to perform the data integration and discuss the structure of the venom components-related transcripts identified. Concluding, this research paves the way for novel investigations and discovery of future pharmacological agents or targets in the antivenom therapy.

The high molecular weight dipeptidyl peptidase IV Pol d 3 is a major allergen of Polistes dominula venom

Hymenoptera venom allergy can cause severe anaphylaxis in untreated patients. Polistes dominula is an important elicitor of venom allergy in Southern Europe as well as in the United States. Due to its increased spreading to more moderate climate zones, Polistes venom allergy is likely to gain importance also in these areas. So far, only few allergens of Polistes dominula venom were identified as basis for component-resolved diagnostics. Therefore, this study aimed to broaden the available panel of important Polistes venom allergens. The 100 kDa allergen Pol d 3 was identified by mass spectrometry and found to be a dipeptidyl peptidase IV. Recombinantly produced Pol d 3 exhibited sIgE-reactivity with approximately 66% of Polistes venom-sensitized patients. Moreover, its clinical relevance was supported by the potent activation of basophils from allergic patients. Cross-reactivity with the dipeptidyl peptidases IV from honeybee and yellow jacket venom suggests the presence of exclusive as well as conserved IgE epitopes. The obtained data suggest a pivotal role of Pol d 3 as sensitizing component of Polistes venom, thus supporting its status as a major allergen of clinical relevance. Therefore, Pol d 3 might become a key element for proper diagnosis of Polistes venom allergy.

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

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

Comparative venom gland transcriptomics of Naja kaouthia (monocled cobra) from Malaysia and Thailand: elucidating geographical venom variation and insights into sequence novelty

Background

The monocled cobra (Naja kaouthia) is a medically important venomous snake in Southeast Asia. Its venom has been shown to vary geographically in relation to venom composition and neurotoxic activity, indicating vast diversity of the toxin genes within the species. To investigate the polygenic trait of the venom and its locale-specific variation, we profiled and compared the venom gland transcriptomes of N. kaouthia from Malaysia (NK-M) and Thailand (NK-T) applying next-generation sequencing (NGS) technology.

Methods

The transcriptomes were sequenced on the Illumina HiSeq platform, assembled and followed by transcript clustering and annotations for gene expression and function. Pairwise or multiple sequence alignments were conducted on the toxin genes expressed. Substitution rates were studied for the major toxins co-expressed in NK-M and NK-T.

Results and discussion

The toxin transcripts showed high redundancy (41–82% of the total mRNA expression) and comprised 23 gene families expressed in NK-M and NK-T, respectively (22 gene families were co-expressed). Among the venom genes, three-finger toxins (3FTxs) predominated in the expression, with multiple sequences noted. Comparative analysis and selection study revealed that 3FTxs are genetically conserved between the geographical specimens whilst demonstrating distinct differential expression patterns, implying gene up-regulation for selected principal toxins, or alternatively, enhanced transcript degradation or lack of transcription of certain traits. One of the striking features that elucidates the inter-geographical venom variation is the up-regulation of α-neurotoxins (constitutes ∼80.0% of toxin’s fragments per kilobase of exon model per million mapped reads (FPKM)), particularly the long-chain α-elapitoxin-Nk2a (48.3%) in NK-T but only 1.7% was noted in NK-M. Instead, short neurotoxin isoforms were up-regulated in NK-M (46.4%). Another distinct transcriptional pattern observed is the exclusively and abundantly expressed cytotoxin CTX-3 in NK-T. The findings suggested correlation with the geographical variation in proteome and toxicity of the venom, and support the call for optimising antivenom production and use in the region. Besides, the current study uncovered full and partial sequences of numerous toxin genes from N. kaouthia which have not been reported hitherto; these include N. kaouthia-specific l-amino acid oxidase (LAAO), snake venom serine protease (SVSP), cystatin, acetylcholinesterase (AChE), hyaluronidase (HYA), waprin, phospholipase B (PLB), aminopeptidase (AP), neprilysin, etc. Taken together, the findings further enrich the snake toxin database and provide deeper insights into the genetic diversity of cobra venom toxins.

Global Transcriptome Analysis of the Tentacle of the Jellyfish Cyanea capillata Using Deep Sequencing and Expressed Sequence Tags: Insight into the Toxin- and Degenerative Disease-Related Transcripts

BACKGROUND

Jellyfish contain diverse toxins and other bioactive components. However, large-scale identification of novel toxins and bioactive components from jellyfish has been hampered by the low efficiency of traditional isolation and purification methods.

RESULTS

We performed de novo transcriptome sequencing of the tentacle tissue of the jellyfish Cyanea capillata. A total of 51,304,108 reads were obtained and assembled into 50,536 unigenes. Of these, 21,357 unigenes had homologues in public databases, but the remaining unigenes had no significant matches due to the limited sequence information available and species-specific novel sequences. Functional annotation of the unigenes also revealed general gene expression profile characteristics in the tentacle of C. capillata. A primary goal of this study was to identify putative toxin transcripts. As expected, we screened many transcripts encoding proteins similar to several well-known toxin families including phospholipases, metalloproteases, serine proteases and serine protease inhibitors. In addition, some transcripts also resembled molecules with potential toxic activities, including cnidarian CfTX-like toxins with hemolytic activity, plancitoxin-1, venom toxin-like peptide-6, histamine-releasing factor, neprilysin, dipeptidyl peptidase 4, vascular endothelial growth factor A, angiotensin-converting enzyme-like and endothelin-converting enzyme 1-like proteins. Most of these molecules have not been previously reported in jellyfish. Interestingly, we also characterized a number of transcripts with similarities to proteins relevant to several degenerative diseases, including Huntington's, Alzheimer's and Parkinson's diseases. This is the first description of degenerative disease-associated genes in jellyfish.

CONCLUSION

We obtained a well-categorized and annotated transcriptome of C. capillata tentacle that will be an important and valuable resource for further understanding of jellyfish at the molecular level and information on the underlying molecular mechanisms of jellyfish stinging. The findings of this study may also be used in comparative studies of gene expression profiling among different jellyfish species.

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

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

Quantitative high-throughput profiling of snake venom gland transcriptomes and proteomes (Ovophis okinavensis and Protobothrops flavoviridis)

Background

Advances in DNA sequencing and proteomics have facilitated quantitative comparisons of snake venom composition. Most studies have employed one approach or the other. Here, both Illumina cDNA sequencing and LC/MS were used to compare the transcriptomes and proteomes of two pit vipers, Protobothrops flavoviridis and Ovophis okinavensis, which differ greatly in their biology.

Results

Sequencing of venom gland cDNA produced 104,830 transcripts. The Protobothrops transcriptome contained transcripts for 103 venom-related proteins, while the Ovophis transcriptome contained 95. In both, transcript abundances spanned six orders of magnitude. Mass spectrometry identified peptides from 100% of transcripts that occurred at higher than contaminant (e.g. human keratin) levels, including a number of proteins never before sequenced from snakes. These transcriptomes reveal fundamentally different envenomation strategies. Adult Protobothrops venom promotes hemorrhage, hypotension, incoagulable blood, and prey digestion, consistent with mammalian predation. Ovophis venom composition is less readily interpreted, owing to insufficient pharmacological data for venom serine and metalloproteases, which comprise more than 97.3% of Ovophis transcripts, but only 38.0% of Protobothrops transcripts. Ovophis venom apparently represents a hybrid strategy optimized for frogs and small mammals.

Conclusions

This study illustrates the power of cDNA sequencing combined with MS profiling. The former quantifies transcript composition, allowing detection of novel proteins, but cannot indicate which proteins are actually secreted, as does MS. We show, for the first time, that transcript and peptide abundances are correlated. This means that MS can be used for quantitative, non-invasive venom profiling, which will be beneficial for studies of endangered species.

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.

Dipeptidyl peptidase IV activity in the blood of the American alligator (Alligator mississippiensis)

Dipeptidyl peptidase IV (DPP4) enzyme activity, associated with cell surface and soluble CD26, was measured in alligator plasma and whole blood. DPP4 activity was higher in whole blood than in the plasma, presumably due to the inclusion of both membrane-bound CD26 on the surface of T-cells and sCD26 in the plasma. The plasma DPP4 activity was measured within 5 min after the addition of the substrate, and maximal accumulation of product was reached at 120 min. Alligator whole blood DPP4 activity peaked at 120 min. In addition, the DPP4 activity in both plasma and whole blood was temperature-dependent, with moderate activities at lower temperatures (5-20 degrees C) and peak activity at 40 degrees C. Alligator plasma DPP4 displayed classical Michaelis-Menten kinetics, with a V(max) value of 267+/-9 nmol/min, and a K(m) of 21.1+/-1.9 microM, while these values were 286+/-14 nmol/min and 57.1+/-3.7 microM, respectively, for whole blood. The proteolytic activity was inhibited in a concentration-dependent manner by a specific DPP4 inhibitor, diprotin A, which indicates that the activity is probably due to the presence of DPP4. The significance of the presence of this activity in alligator blood is still unknown.