Dramatic expansion of the black widow toxin arsenal uncovered by multi-tissue transcriptomics and venom proteomics.

BackgroundAnimal venoms attract enormous interest given their potential for pharmacological discovery and understanding the evolution of natural chemistries. Next-generation transcriptomics and proteomics provide unparalleled, but underexploited, capabilities for venom characterization. We combined multi-tissue RNA-Seq with mass spectrometry and bioinformatic analyses to determine venom gland specific transcripts and venom proteins from the Western black widow spider (Latrodectus hesperus) and investigated their evolution.ResultsWe estimated expression of 97,217 L. hesperus transcripts in venom glands relative to silk and cephalothorax tissues. We identified 695 venom gland specific transcripts (VSTs), many of which BLAST and GO term analyses indicate may function as toxins or their delivery agents. ~38% of VSTs had BLAST hits, including latrotoxins, inhibitor cystine knot toxins, CRISPs, hyaluronidases, chitinase, and proteases, and 59% of VSTs had predicted protein domains. Latrotoxins are venom toxins that cause massive neurotransmitter release from vertebrate or invertebrate neurons. We discovered ≥ 20 divergent latrotoxin paralogs expressed in L. hesperus venom glands, significantly increasing this biomedically important family. Mass spectrometry of L. hesperus venom identified 49 proteins from VSTs, 24 of which BLAST to toxins. Phylogenetic analyses showed venom gland specific gene family expansions and shifts in tissue expression.ConclusionsQuantitative expression analyses comparing multiple tissues are necessary to identify venom gland specific transcripts. We present a black widow venom specific exome that uncovers a trove of diverse toxins and associated proteins, suggesting a dynamic evolutionary history. This justifies a reevaluation of the functional activities of black widow venom in light of its emerging complexity

Comparative venom-gland transcriptomics and venom proteomics of four Sidewinder Rattlesnake (\u3ci\u3eCrotalus cerastes\u3c/i\u3e) lineages reveal little differential expression despite individual variation

Changes in gene expression can rapidly influence adaptive traits in the early stages of lineage diversification. Venom is an adaptive trait comprised of numerous toxins used for prey capture and defense. Snake venoms can vary widely between conspecific populations, but the influence of lineage diversification on such compositional differences are unknown. To explore venom differentiation in the early stages of lineage diversification, we used RNA-seq and mass spectrometry to characterize Sidewinder Rattlesnake (Crotalus cerastes) venom. We generated the first venom-gland transcriptomes and complementary venom proteomes for eight individuals collected across the United States and tested for expression differences across life history traits and between subspecific, mitochondrial, and phylotranscriptomic hypotheses. Sidewinder venom was comprised primarily of hemorrhagic toxins, with few cases of differential expression attributable to life history or lineage hypotheses. However, phylotranscriptomic lineage comparisons more than doubled instances of significant expression differences compared to all other factors. Nevertheless, only 6.4% of toxins were differentially expressed overall, suggesting that shallow divergence has not led to major changes in Sidewinder venom composition. Our results demonstrate the need for consensus venom-gland transcriptomes based on multiple individuals and highlight the potential for discrepancies in differential expression between different phylogenetic hypotheses

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

A proteomics and transcriptomics investigation of the venom from the Barychelid spider Trittame loki (brush-foot trapdoor)

Although known for their potent venom and ability to prey upon both invertebrate and vertebrate species, the Barychelidae spider family has been entirely neglected by toxinologists. In striking contrast, the sister family Theraphosidae (commonly known as tarantulas), which last shared a most recent common ancestor with Barychelidae over 200 million years ago, has received much attention, accounting for 25% of all the described spider toxins while representing only 2% of all spider species. In this study, we evaluated for the first time the venom arsenal of a barychelid spider, Trittame loki, using transcriptomic, proteomic, and bioinformatic methods. The venom was revealed to be dominated by extremely diverse inhibitor cystine knot (ICK)/knottin peptides, accounting for 42 of the 46 full-length toxin precursors recovered in the transcriptomic sequencing. In addition to documenting differential rates of evolution adopted by different ICK/knottin toxin lineages, we discovered homologues with completely novel cysteine skeletal architecture. Moreover, acetylcholinesterase and neprilysin were revealed for the first time as part of the spider-venom arsenal and CAP (CRiSP/Allergen/PR-1) were identified for the first time in mygalomorph spider venoms. These results not only highlight the extent of venom diversification in this neglected ancient spider lineage, but also reinforce the idea that unique venomous lineages are rich pools of novel biomolecules that may have significant applied uses as therapeutics and/or insecticides

Chaperone-mediated native folding of a β-scorpion toxin in the periplasm of E.coli

Background: Animal neurotoxin peptides are valuable probes for investigating ion channel structure/function relationships and represent lead compounds for novel therapeutics and insecticides. However, misfolding and aggregation are common outcomes when toxins containing multiple disulfides are expressed in bacteria. Methods: The ß-scorpion peptide toxin Bj-xtrIT from Hottentotta judaica and four chaperone enzymes (DsbA, DsbC, SurA and FkpA) were co-secreted into the oxidizing environment of the E.coli periplasm. Expressed Bj-xtrIT was purified and analyzed by HPLC and FPLC chromatography. Its thermostability was assessed using synchrotron radiation circular dichroism spectroscopy and its crystal structure was determined. Results: Western blot analysis showed that robust expression was only achieved when cells co-expressed the chaperones. The purified samples were homogenous and monodisperse and the protein was thermostable. The crystal structure of the recombinant toxin confirmed that it adopts the native disulfide connectivity and fold. Conclusions: The chaperones enabled correct folding of the four-disulfide-bridged Bj-xtrIT toxin. There was no apparent sub-population of misfolded Bj-xtrIT, which attests to the effectiveness of this expression method. General Significance: We report the first example of a disulfide-linked scorpion toxin natively folded during bacterial expression. This method eliminates downstream processing steps such as oxidative refolding or cleavage of a fusion-carrier and therefore enables efficient production of insecticidal Bj-xtrIT. Periplasmic chaperone activity may produce native folding of other extensively disulfide-reticulated proteins including animal neurotoxins. This work is therefore relevant to venomics and studies of a wide range of channels and receptors

Trait differentiation and modular toxin expression in palm-pitvipers

Background Modularity is the tendency for systems to organize into semi-independent units and can be a key to the evolution and diversification of complex biological systems. Snake venoms are highly variable modular systems that exhibit extreme diversification even across very short time scales. One well-studied venom phenotype dichotomy is a trade-off between neurotoxicity versus hemotoxicity that occurs through the high expression of a heterodimeric neurotoxic phospholipase A2 (PLA2) or snake venom metalloproteinases (SVMPs). We tested whether the variation in these venom phenotypes could occur via variation in regulatory sub-modules through comparative venom gland transcriptomics of representative Black-Speckled Palm-Pitvipers (Bothriechis nigroviridis) and Talamancan Palm-Pitvipers (B. nubestris). Results We assembled 1517 coding sequences, including 43 toxins for B. nigroviridis and 1787 coding sequences including 42 toxins for B. nubestris. The venom gland transcriptomes were extremely divergent between these two species with one B. nigroviridis exhibiting a primarily neurotoxic pattern of expression, both B. nubestris expressing primarily hemorrhagic toxins, and a second B. nigroviridis exhibiting a mixed expression phenotype. Weighted gene coexpression analyses identified six submodules of transcript expression variation, one of which was highly associated with SVMPs and a second which contained both subunits of the neurotoxic PLA2 complex. The sub-module association of these toxins suggest common regulatory pathways underlie the variation in their expression and is consistent with known patterns of inheritance of similar haplotypes in other species. We also find evidence that module associated toxin families show fewer gene duplications and transcript losses between species, but module association did not appear to affect sequence diversification. Conclusion Sub-modular regulation of expression likely contributes to the diversification of venom phenotypes within and among species and underscores the role of modularity in facilitating rapid evolution of complex traits

Transcriptome and venom proteome of the box jellyfish Chironex fleckeri

Background: The box jellyfish, Chironex fleckeri, is the largest and most dangerous cubozoan jellyfish to humans. It produces potent and rapid-acting venom and its sting causes severe localized and systemic effects that are potentially life-threatening. In this study, a combined transcriptomic and proteomic approach was used to identify C. fleckeri proteins that elicit toxic effects in envenoming. Results: More than 40,000,000 Illumina reads were used to de novo assemble ~34,000 contiguous cDNA sequences and ~20,000 proteins were predicted based on homology searches, protein motifs, gene ontology and biological pathway mapping. More than 170 potential toxin proteins were identified from the transcriptome on the basis of homology to known toxins in publicly available sequence databases. MS/MS analysis of C. fleckeri venom identified over 250 proteins, including a subset of the toxins predicted from analysis of the transcriptome. Potential toxins identified using MS/MS included metalloproteinases, an alpha-macroglobulin domain containing protein, two CRISP proteins and a turripeptide-like protease inhibitor. Nine novel examples of a taxonomically restricted family of potent cnidarian pore-forming toxins were also identified. Members of this toxin family are potently haemolytic and cause pain, inflammation, dermonecrosis, cardiovascular collapse and death in experimental animals, suggesting that these toxins are responsible for many of the symptoms of C. fleckeri envenomation. Conclusions: This study provides the first overview of a box jellyfish transcriptome which, coupled with venom proteomics data, enhances our current understanding of box jellyfish venom composition and the molecular structure and function of cnidarian toxins. The generated data represent a useful resource to guide future comparative studies, novel protein/peptide discovery and the development of more effective treatments for jellyfish stings in humans. (Length: 300)

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

Artículo 10 páginas, 3 figuras 1 tablaUnderstanding the origin and maintenance of phenotypic variation, particularly across a continuous spatial distribution, represents a key challenge in evolutionary biology. For this, animal venoms represent ideal study systems: they are complex, variable, yet easily quantifiable molecular phenotypes with a clear function. Rattlesnakes display tremendous variation in their venom composition, mostly through strongly dichotomous venom strategies, which may even coexist within a single species. Here, through dense, widespread population-level sampling of the Mojave rattlesnake, Crotalus scutulatus, we show that genomic structural variation at multiple loci underlies extreme geographical variation in venom composition, which is maintained despite extensive gene flow. Unexpectedly, neither diet composition nor neutral population structure explain venom variation. Instead, venom divergence is strongly correlated with environmental conditions. Individual toxin genes correlate with distinct environmental factors, suggesting that different selective pressures can act on individual loci independently of their co-expression patterns or genomic proximity. Our results challenge common assumptions about diet composition as the key selective driver of snake venom evolution and emphasize how the interplay between genomic architecture and local-scale spatial heterogeneity in selective pressures may facilitate the retention of adaptive functional polymorphisms across a continuous space.Funding: Leverhulme Trust Grant RPG 2013-315 to WW, Santander Early Career Research Scholarship to GZ, Ministerio de Economía y Competitividad Grant BFU2013-42833-P to JJC.Peer reviewe

Unusual accelerated rate of deletions and insertions in toxin genes in the venom glands of the pygmy copperhead (Austrelaps labialis) from kangaroo island

<p>Abstract</p> <p>Background</p> <p>Toxin profiling helps in cataloguing the toxin present in the venom as well as in searching for novel toxins. The former helps in understanding potential pharmacological profile of the venom and evolution of toxins, while the latter contributes to understanding of novel mechanisms of toxicity and provide new research tools or prototypes of therapeutic agents.</p> <p>Results</p> <p>The pygmy copperhead (<it>Austrelaps labialis</it>) is one of the less studied species. In this present study, an attempt has been made to describe the toxin profile of <it>A. labialis </it>from Kangaroo Island using the cDNA library of its venom glands. We sequenced 658 clones which represent the common families of toxin genes present in snake venom. They include (a) putative long-chain and short-chain neurotoxins, (b) phospholipase A₂, (c) Kunitz-type protease inhibitor, (d) CRISPs, (e) C-type lectins and (f) Metalloproteases. In addition, we have also identified a novel protein with two Kunitz-type domains in tandem similar to bikunin.</p> <p>Conclusion</p> <p>Interestingly, the cDNA library reveals that most of the toxin families (17 out of 43 toxin genes; ~40%) have truncated transcripts due to insertion or deletion of nucleotides. These truncated products might not be functionally active proteins. However, cellular trancripts from the same venom glands are not affected. This unusual higher rate of deletion and insertion of nucleotide in toxin genes may be responsible for the lower toxicity of <it>A. labialis </it>venom of Kangroo Island and have significant effect on evolution of toxin genes.</p