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

<p>Abstract</p> <p>Background</p> <p>The genus <it>Bothrops </it>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 <it>Bothrops alternatus</it>, a clinically important species found in southeastern and southern Brazil, Uruguay, northern Argentina and eastern Paraguay.</p> <p>Results</p> <p>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.</p> <p>Conclusions</p> <p><it>Bothrops alternatus </it>venom gland contains the major toxin classes described for other <it>Bothrops </it>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 <it>Bothrops </it>species with these toxins. Together, these results contribute to our understanding of the physiopathology of envenoming by this species.</p

The Indian cobra reference genome and transcriptome enables comprehensive identification of venom toxins

Snakebite envenoming is a serious and neglected tropical disease that kills ~100,000 people annually. High-quality, genome-enabled comprehensive characterization of toxin genes will facilitate development of effective humanized recombinant antivenom. We report a de novo near-chromosomal genome assembly of Naja naja, the Indian cobra, a highly venomous, medically important snake. Our assembly has a scaffold N50 of 223.35 Mb, with 19 scaffolds containing 95% of the genome. Of the 23,248 predicted protein-coding genes, 12,346 venom-gland-expressed genes constitute the \u27venom-ome\u27 and this included 139 genes from 33 toxin families. Among the 139 toxin genes were 19 \u27venom-ome-specific toxins\u27 (VSTs) that showed venom-gland-specific expression, and these probably encode the minimal core venom effector proteins. Synthetic venom reconstituted through recombinant VST expression will aid in the rapid development of safe and effective synthetic antivenom. Additionally, our genome could serve as a reference for snake genomes, support evolutionary studies and enable venom-driven drug discovery