Molecular evolution of toxin genes in Elapidae snakes

Distinct regulatory networks control toxin gene expression in elapid and viperid snakes

BACKGROUND

Venom systems are ideal models to study genetic regulatory mechanisms that underpin evolutionary novelty. Snake venom glands are thought to share a common origin, but there are major distinctions between venom toxins from the medically significant snake families Elapidae and Viperidae, and toxin gene regulatory investigations in elapid snakes have been limited. Here, we used high-throughput RNA-sequencing to profile gene expression and microRNAs between active (milked) and resting (unmilked) venom glands in an elapid (Eastern Brown Snake, Pseudonaja textilis), in addition to comparative genomics, to identify cis- and trans-acting regulation of venom production in an elapid in comparison to viperids (Crotalus viridis and C. tigris).

RESULTS

Although there is conservation in high-level mechanistic pathways regulating venom production (unfolded protein response, Notch signaling and cholesterol homeostasis), there are differences in the regulation of histone methylation enzymes, transcription factors, and microRNAs in venom glands from these two snake families. Histone methyltransferases and transcription factor (TF) specificity protein 1 (Sp1) were highly upregulated in the milked elapid venom gland in comparison to the viperids, whereas nuclear factor I (NFI) TFs were upregulated after viperid venom milking. Sp1 and NFI cis-regulatory elements were common to toxin gene promoter regions, but many unique elements were also present between elapid and viperid toxins. The presence of Sp1 binding sites across multiple elapid toxin gene promoter regions that have been experimentally determined to regulate expression, in addition to upregulation of Sp1 after venom milking, suggests this transcription factor is involved in elapid toxin expression. microRNA profiles were distinctive between milked and unmilked venom glands for both snake families, and microRNAs were predicted to target a diversity of toxin transcripts in the elapid P. textilis venom gland, but only snake venom metalloproteinase transcripts in the viperid C. viridis venom gland. These results suggest differences in toxin gene posttranscriptional regulation between the elapid P. textilis and viperid C. viridis.

CONCLUSIONS

Our comparative transcriptomic and genomic analyses between toxin genes and isoforms in elapid and viperid snakes suggests independent toxin regulation between these two snake families, demonstrating multiple different regulatory mechanisms underpin a venomous phenotype.

Toxinology provides multidirectional and multidimensional opportunities: A personal perspective

In nature, toxins have evolved as weapons to capture and subdue the prey or to counter predators or competitors. When they are inadvertently injected into humans, they cause symptoms ranging from mild discomfort to debilitation and death. Toxinology is the science of studying venoms and toxins that are produced by a wide variety of organisms. In the past, the structure, function and mechanisms of most abundant and/or most toxic components were characterized to understand and to develop strategies to neutralize their toxicity. With recent technical advances, we are able to evaluate and determine the toxin profiles using transcriptomes of venom glands and proteomes of tiny amounts of venom. Enormous amounts of data from these studies have opened tremendous opportunities in many directions of basic and applied research. The lower costs for profiling venoms will further fuel the expansion of toxin database, which in turn will provide greater exciting and bright opportunities in toxin research.

[Venoms and medical research]

Animal venoms are complex chemical cocktails, comprising a wide range of biologically active reticulated peptides that target with high selectivity and efficacy a variety of enzymes, membrane receptors, ion channels...Venoms can therefore be seen as large natural libraries of biologically active molecules that are continuously selected and highly refined by the evolution process, up to the point where every molecule is endowed with pharmacological properties that are highly valuable in the context of human use and drug development. Therefore, venom exploration constitutes a prerequisite to drug discovery. However, mass spectrometry and transcriptomics via NGS (Next Generation Sequencing) studies have shown the presence of up to 1000 peptides in the venom of single species of cone snails and spiders. Therefore the global animal venom resource can be seen as a collection of more than 50 to 100 000 000 peptides and proteins of which only ~5000 are known. That extraordinary "Eldorado" of bio-optimized compounds justifies the development of more global and cutting-edge strategies and technologies to explore this resource more efficiently than actually. De novo developed approaches and recently obtained results will be described.

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.

Molding the business end of neurotoxins by diversifying evolution

A diverse range of organisms utilize neurotoxins that target specific ion channels and modulate their activity. Typically, toxins are clustered into several multigene families, providing an organism with the upper hand in the never-ending predator-prey arms race. Several gene families, including those encoding certain neurotoxins, have been subject to diversifying selection forces, resulting in rapid gene evolution. Here we sought a spatial pattern in the distribution of both diversifying and purifying selection forces common to neurotoxin gene families. Utilizing the mechanistic empirical combination model, we analyzed various toxin families from different phyla affecting various receptors and relying on diverse modes of action. Through this approach, we were able to detect clear correlations between the pharmacological surface of a toxin and rapidly evolving domains, rich in positively selected residues. On the other hand, patches of negatively selected residues were restricted to the nontoxic face of the molecule and most likely help in stabilizing the tertiary structure of the toxin. We thus propose a mutual evolutionary strategy of venomous animals in which adaptive molecular evolution is directed toward the toxin active surface. Furthermore, we propose that the binding domains of unstudied toxins could be readily predicted using evolutionary considerations.

Venom gland transcriptomes of two elapid snakes (Bungarus multicinctus and Naja atra) and evolution of toxin genes

BACKGROUND

Kraits (genus Bungarus) and cobras (genus Naja) are two representative toxic genera of elapids in the old world. Although they are closely related genera and both of their venoms are very toxic, the compositions of their venoms are very different. To unveil their detailed venoms and their evolutionary patterns, we constructed venom gland cDNA libraries and genomic bacterial artificial chromosome (BAC) libraries for Bungarus multicinctus and Naja atra, respectively. We sequenced about 1500 cDNA clones for each of the venom cDNA libraries and screened BAC libraries of the two snakes by blot analysis using four kinds of toxin probes; i.e., three-finger toxin (3FTx), phospholipase A2 (PLA2), kunitz-type protease inhibitor (Kunitz), and natriuretic peptide (NP).

RESULTS

In total, 1092 valid expressed sequences tags (ESTs) for B. multicinctus and 1166 ESTs for N. atra were generated. About 70% of these ESTs can be annotated as snake toxin transcripts. 3FTx (64.5%) and β bungarotoxin (25.1%) comprise the main toxin classes in B. multicinctus, while 3FTx (95.8%) is the dominant toxin in N. atra. We also observed several less abundant venom families in B. multicinctus and N. atra, such as PLA2, C-type lectins, and Kunitz. Peculiarly a cluster of NP precursors with tandem NPs was detected in B. multicinctus. A total of 71 positive toxin BAC clones in B. multicinctus and N. atra were identified using four kinds of toxin probes (3FTx, PLA2, Kunitz, and NP), among which 39 3FTx-positive BACs were sequenced to reveal gene structures of 3FTx toxin genes.

CONCLUSIONS

Based on the toxin ESTs and 3FTx gene sequences, the major components of B. multicinctus venom transcriptome are neurotoxins, including long chain alpha neurotoxins (α-ntx) and the recently originated β bungarotoxin, whereas the N. atra venom transcriptome mainly contains 3FTxs with cytotoxicity and neurotoxicity (short chain α-ntx). The data also revealed that tandem duplications contributed the most to the expansion of toxin multigene families. Analysis of nonsynonymous to synonymous nucleotide substitution rate ratios (dN/dS) indicates that not only multigene toxin families but also other less abundant toxins might have been under rapid diversifying evolution.

The Toxicogenomic Multiverse: Convergent Recruitment of Proteins Into Animal Venoms

Throughout evolution, numerous proteins have been convergently recruited into the venoms of various animals, including centipedes, cephalopods, cone snails, fish, insects (several independent venom systems), platypus, scorpions, shrews, spiders, toxicoferan reptiles (lizards and snakes), and sea anemones. The protein scaffolds utilized convergently have included AVIT/colipase/prokineticin, CAP, chitinase, cystatin, defensins, hyaluronidase, Kunitz, lectin, lipocalin, natriuretic peptide, peptidase S1, phospholipase A2, sphingomyelinase D, and SPRY. Many of these same venom protein types have also been convergently recruited for use in the hematophagous gland secretions of invertebrates (e.g., fleas, leeches, kissing bugs, mosquitoes, and ticks) and vertebrates (e.g., vampire bats). Here, we discuss a number of overarching structural, functional, and evolutionary generalities of the protein families from which these toxins have been frequently recruited and propose a revised and expanded working definition for venom. Given the large number of striking similarities between the protein compositions of conventional venoms and hematophagous secretions, we argue that the latter should also fall under the same definition.

QSAR for RNases and theoretic-experimental study of molecular diversity on peptide mass fingerprints of a new Leishmania infantum protein

The toxicity and low success of current treatments for Leishmaniosis determines the search of new peptide drugs and/or molecular targets in Leishmania pathogen species (L. infantum and L. major). For example, Ribonucleases (RNases) are enzymes relevant to several biologic processes; then, theoretical and experimental study of the molecular diversity of Peptide Mass Fingerprints (PMFs) of RNases is useful for drug design. This study introduces a methodology that combines QSAR models, 2D-Electrophoresis (2D-E), MALDI-TOF Mass Spectroscopy (MS), BLAST alignment, and Molecular Dynamics (MD) to explore PMFs of RNases. We illustrate this approach by investigating for the first time the PMFs of a new protein of L. infantum. Here we report and compare new versus old predictive models for RNases based on Topological Indices (TIs) of Markov Pseudo-Folding Lattices. These group of indices called Pseudo-folding Lattice 2D-TIs include: Spectral moments pi ( k )(x,y), Mean Electrostatic potentials xi ( k )(x,y), and Entropy measures theta ( k )(x,y). The accuracy of the models (training/cross-validation) was as follows: xi ( k )(x,y)-model (96.0%/91.7%)>pi ( k )(x,y)-model (84.7/83.3) > theta ( k )(x,y)-model (66.0/66.7). We also carried out a 2D-E analysis of biological samples of L. infantum promastigotes focusing on a 2D-E gel spot of one unknown protein with M<20, 100 and pI <7. MASCOT search identified 20 proteins with Mowse score >30, but not one >52 (threshold value), the higher value of 42 was for a probable DNA-directed RNA polymerase. However, we determined experimentally the sequence of more than 140 peptides. We used QSAR models to predict RNase scores for these peptides and BLAST alignment to confirm some results. We also calculated 3D-folding TIs based on MD experiments and compared 2D versus 3D-TIs on molecular phylogenetic analysis of the molecular diversity of these peptides. This combined strategy may be of interest in drug development or target identification.

Assembling an arsenal, the scorpion way

Background

For survival, scorpions depend on a wide array of short neurotoxic polypeptides. The venoms of scorpions from the most studied group, the Buthida, are a rich source of small, 23–78 amino acid-long peptides, well packed by either three or four disulfide bridges that affect ion channel function in excitable and non-excitable cells.

Results

In this work, by constructing a toxin transcripts data set from the venom gland of the scorpion Buthus occitanus israelis, we were able to follow the evolutionary path leading to mature toxin diversification and suggest a mechanism for leader peptide hyper-conservation. Toxins from each family were more closely related to one another than to toxins from other species, implying that fixation of duplicated genes followed speciation, suggesting early gene conversion events. Upon fixation, the mature toxin-coding domain was subjected to diversifying selection resulting in a significantly higher substitution rate that can be explained solely by diversifying selection. In contrast to the mature peptide, the leader peptide sequence was hyper-conserved and characterized by an atypical sub-neutral synonymous substitution rate. We interpret this as resulting from purifying selection acting on both the peptide and, as reported here for the first time, the DNA sequence, to create a toxin family-specific codon bias.

Conclusion

We thus propose that scorpion toxin genes were shaped by selective forces acting at three levels, namely (1) diversifying the mature toxin, (2) conserving the leader peptide amino acid sequence and intriguingly, (3) conserving the leader DNA sequences.

Accelerated exchange of exon segments in Viperid three-finger toxin genes (Sistrurus catenatus edwardsii; Desert Massasauga)

Background

Snake venoms consist primarily of proteins and peptides showing a myriad of potent biological activities which have been shaped by both adaptive and neutral selective forces. Venom proteins are encoded by multigene families that have evolved through a process of gene duplication followed by accelerated evolution in the protein coding region.

Results

Here we report five gene structures of three-finger toxins from a viperid snake, Sistrurus catenatus edwardsii. These toxin genes are structured similarly to elapid and hydrophiid three-finger toxin genes, with two introns and three exons. Both introns and exons show distinct patterns of segmentation, and the insertion/deletion of segments may define their evolutionary history. The segments in introns, when present, are highly similar to their corresponding segments in other members of the gene family. In contrast, some segments in the exons show high similarity, while others are often distinctly different among corresponding regions of the isoforms.

Conclusion

Ordered, conserved exon structure strongly suggests that segments in corresponding regions in exons have been exchanged with distinctly different ones during the evolution of these genes. Such a "switching" of segments in exons may result in drastically altering the molecular surface topology and charge, and hence the molecular targets of these three-finger toxins. Thus the phenomenon of accelerated segment switch in exons to alter targeting (ASSET) may play an important role in the evolution of three-finger toxins, resulting in a family of toxins with a highly conserved structural fold but widely varying biological activities.

Diversity of Toxic Components from the Venom of the Evolutionarily Distinct Black Whip Snake, Demansia vestigiata

Included among the more than 300 species of elapid snakes worldwide is the Australian genus Demansia, or whip snakes. Despite evidence to suggest adverse clinical outcomes from envenomation by these snakes, together with confusion on their true phylogenetic relationship to other Australian elapids, not a single toxin sequence has previously been reported from the venom of a Demansia species. We describe here a combined proteomic and transcriptomic approach characterizing the venom from the black whip snake, Demansia vestigiata. A total of 13 distinct toxin families were identified, including homologues of all of the major toxic components previously reported from the venom of other Australian elapids, such as factor X-like prothrombin activators, neurotoxins, phospholipases, cysteine rich secretory proteins, textilinin-like molecules, nerve growth factors, l-amino acid oxidases, vespryns, 5' nucleotidases, metalloproteinases, and C-type lectins as well as a novel dipeptidyl peptidase family. Phylogenetic analysis of these sequences revealed an early evolutionary split of the black whip snake from all other characterized Australian snakes, with a low degree of sequence identity between D. vestigiata and the other snakes, across all toxin families. The results of this study have important implications not only for the further characterization of venom from whip snakes, but also for our understanding of the evolutionary relationship of Australian snake species.