The venom-gland transcriptome of the eastern coral snake (Micrurus fulvius) reveals high venom complexity in the intragenomic evolution of venoms

Expression of venom gene homologs in diverse python tissues suggests a new model for the evolution of snake venom

Snake venom gene evolution has been studied intensively over the past several decades, yet most previous studies have lacked the context of complete snake genomes and the full context of gene expression across diverse snake tissues. We took a novel approach to studying snake venom evolution by leveraging the complete genome of the Burmese python, including information from tissue-specific patterns of gene expression. We identified the orthologs of snake venom genes in the python genome, and conducted detailed analysis of gene expression of these venom homologs to identify patterns that differ between snake venom gene families and all other genes. We found that venom gene homologs in the python are expressed in many different tissues outside of oral glands, which illustrates the pitfalls of using transcriptomic data alone to define "venom toxins." We hypothesize that the python may represent an ancestral state prior to major venom development, which is supported by our finding that the expansion of venom gene families is largely restricted to highly venomous caenophidian snakes. Therefore, the python provides insight into biases in which genes were recruited for snake venom systems. Python venom homologs are generally expressed at lower levels, have higher variance among tissues, and are expressed in fewer organs compared with all other python genes. We propose a model for the evolution of snake venoms in which venom genes are recruited preferentially from genes with particular expression profile characteristics, which facilitate a nearly neutral transition toward specialized venom system expression.

A Polychaete's powerful punch: venom gland transcriptomics of Glycera reveals a complex cocktail of toxin homologs

Glycerids are marine annelids commonly known as bloodworms. Bloodworms have an eversible proboscis adorned with jaws connected to venom glands. Bloodworms prey on invertebrates, and it is known that the venom glands produce compounds that can induce toxic effects in animals. Yet, none of these putative toxins has been characterized on a molecular basis. Here we present the transcriptomic profiles of the venom glands of three species of bloodworm, Glycera dibranchiata, Glycera fallax and Glycera tridactyla, as well as the body tissue of G. tridactyla. The venom glands express a complex mixture of transcripts coding for putative toxin precursors. These transcripts represent 20 known toxin classes that have been convergently recruited into animal venoms, as well as transcripts potentially coding for Glycera-specific toxins. The toxins represent five functional categories: Pore-forming and membrane-disrupting toxins, neurotoxins, protease inhibitors, other enzymes, and CAP domain toxins. Many of the transcripts coding for putative Glycera toxins belong to classes that have been widely recruited into venoms, but some are homologs of toxins previously only known from the venoms of scorpaeniform fish and monotremes (stonustoxin-like toxin), turrid gastropods (turripeptide-like peptides), and sea anemones (gigantoxin I-like neurotoxin). This complex mixture of toxin homologs suggests that bloodworms employ venom while predating on macroscopic prey, casting doubt on the previously widespread opinion that G. dibranchiata is a detritivore. Our results further show that researchers should be aware that different assembly methods, as well as different methods of homology prediction, can influence the transcriptomic profiling of venom glands.

Toxin diversity revealed by a transcriptomic study of Ornithoctonus huwena

Spider venom comprises a mixture of compounds with diverse biological activities, which are used to capture prey and defend against predators. The peptide components bind a broad range of cellular targets with high affinity and selectivity, and appear to have remarkable structural diversity. Although spider venoms have been intensively investigated over the past few decades, venomic strategies to date have generally focused on high-abundance peptides. In addition, the lack of complete spider genomes or representative cDNA libraries has presented significant limitations for researchers interested in molecular diversity and understanding the genetic mechanisms of toxin evolution. In the present study, second-generation sequencing technologies, combined with proteomic analysis, were applied to determine the diverse peptide toxins in venom of the Chinese bird spider Ornithoctonus huwena. In total, 626 toxin precursor sequences were retrieved from transcriptomic data. All toxin precursors clustered into 16 gene superfamilies, which included six novel superfamilies and six novel cysteine patterns. A surprisingly high number of hypermutations and fragment insertions/deletions were detected, which accounted for the majority of toxin gene sequences with low-level expression. These mutations contribute to the formation of diverse cysteine patterns and highly variable isoforms. Furthermore, intraspecific venom variability, in combination with variable transcripts and peptide processing, contributes to the hypervariability of toxins in venoms, and associated rapid and adaptive evolution of toxins for prey capture and defense.

Next-generation snake venomics: protein-locus resolution through venom proteome decomplexation

Venom research has been continuously enhanced by technological advances. High-throughput technologies are changing the classical paradigm of hypothesis-driven research to technology-driven approaches. However, the thesis advocated in this paper is that full proteome coverage at locus-specific resolution requires integrating the best of both worlds into a protocol that includes decomplexation of the venom proteome prior to liquid chromatography-tandem mass spectrometry matching against a species-specific transcriptome. This approach offers the possibility of proof-checking the species-specific contig database using proteomics data. Immunoaffinity chromatography constitutes the basis of an antivenomics workflow designed to quantify the extent of cross-reactivity of antivenoms against homologous and heterologous venom toxins. In the author's view, snake venomics and antivenomics form part of a biology-driven conceptual framework to unveil the genesis and natural history of venoms, and their within- and between-species toxicological and immunological divergences and similarities. Understanding evolutionary trends across venoms represents the Rosetta Stone for generating broad-ranging polyspecific antivenoms.

Eastern coral snake Micrurus fulvius venom toxicity in mice is mainly determined by neurotoxic phospholipases A2

Here we show for the first time that the venom from an elapid (Micrurus fulvius) contains three finger toxin (3FTxs) peptides with low toxicity but high content of lethal phospholipases A2 (PLA2). The intravenous venom LD50 in mice was 0.3μg/g. Fractionation on a C18 column yielded 22 fractions; in terms of abundance, 58.3% of them were components of 13-14kDa and 24.9% were molecules of 6-7kDa. Two fractions with PLA2 activity represented 33.4% of the whole venom and were the most lethal fractions. Fractions with low molecular mass (<7000Da) partially and reversibly blocked the nicotinic acetylcholine receptor (nAChR), with the exception of one that blocked it completely. The fraction that blocked 100% contained two protein species whose dose-response was determined; the IC50s were 13±1 and 9.5±0.3nM. Despite the apparent effect on nAChR none of the low molecular mass fractions were lethal in mice, at concentrations of 1μg/g. From 2D-PAGE and LC-MS/MS, we identified fourteen species of PLA2, four protein species of C-type lectin, three zinc metalloproteinases, one phosphodiesterase and one 3FTx. The N-terminal amino acid sequence of fractions with biological interest was obtained.

BIOLOGICAL SIGNIFICANCE

In contrast with coral snake venoms from South America, M. fulvius has minor amounts of low molecular mass components, but high content of PLA2, which is responsible for the venom lethality of this species. The results reported here contribute to better understanding of envenomation development and to improve antivenom design and production. These findings break from the paradigm that neurotoxicity caused by Micrurus venoms is mainly attributable to 3FTx neurotoxins and encourage future studies on Micrurus evolution and venom specialization. This article is part of a Special Issue entitled Non-model organisms.

Proteomic and biochemical analyses of short-tailed pit viper (Gloydius brevicaudus) venom: age-related variation and composition-activity correlation

We conducted an in-depth analysis of the proteomic and biochemical profiles of the venom of neonate and adult short-tailed pit vipers (Gloydius brevicaudus). Identified proteins were assigned to a few main toxin families. Disintegrin, phospholipase A2 (PLA2), serine proteinase, cysteine-rich secretory protein, C-type lectin-like protein, l-amino acid oxidase and snake venom metalloproteinase (SVMP) were detected in both venoms, while 5'-nucleotidase was detected only in the adult venom. SVMP was the predominant protein family in both venoms (neonate: 65.7%; adult: 64.4%), followed by PLA2 (neonate: 13.4%; adult: 25.0%). Antivenomic analysis revealed that commercial G. brevicaudus antivenom almost neutralized the chromatographic peaks with medium and high molecular masses in both venoms, but did not completely recognize peaks with low molecular mass. Toxicological and enzymatic activities show remarkable age-related variation in G. brevicaudus venom, probably resulting from variation in venom composition. Our data demonstrate age-related variation across venomics, antivenomics and biochemical profiles of G. brevicaudus venom, and have implications for the management of G. brevicaudus bites, including improving antivenom preparation by combining both venoms.

BIOLOGICAL SIGNIFICANCE

This study investigates the composition and biochemical activity of neonate and adult Gloydius brevicaudus venoms. We found remarkable age-related variation in venom biological activity, likely the result of variation in venom composition. Antivenomics analysis was used to explore difference in neonate and adult G. brevicaudus venoms. Our findings have implications for the diagnosis and clinical management of G. brevicaudus bites, and the design of venom mixtures that will increase the efficacy of commercial antivenom. This article is part of a Special Issue entitled: Proteomics of non-model organisms.