Discovery of toxin-encoding genes from the false viper Macropisthodon rudis, a rear-fanged snake, by transcriptome analysis of venom gland

Nationwide and long-term epidemiological research of snakebite envenomation in Taiwan during 2002-2014 based on the use of snake antivenoms: A study utilizing National Health Insurance Database

INTRODUCTION

In Taiwan, six venomous snake species with medical importance have been found; however, long-term epidemiological data of snakebite envenomation (SBE) is lacking. This study aimed to explore the epidemiology of SBE based on the distribution and use of different antivenoms in different parts of Taiwan to facilitate the development of prevention strategies and resource allocation.

METHODS AND RESULTS

This retrospective study was conducted using the Taiwan National Health Insurance Research Database from 2002 to 2014. A total of 12,542 patients were treated with antivenoms. The directly standardized cumulative incidence was 3.6 cases per 100,000 individuals based on the 2000 World Standard Population. The incidence of SBEs peaked in the summer (35.9%). The relative risk (RR) of male patients versus female patients was 2.5 (p < 0.0001). The RRs of patients aged 18-64 and ≥65 years versus those aged <18 years were 6.0 (p < 0.0001) and 14.3 (p < 0.0001), respectively. Furthermore, the RR of eastern Taiwan versus northern Taiwan was 6.8 (p < 0.0001). The RR of agricultural workers versus laborers was 5.5 (p < 0.0001). Compared with patients envenomed by Trimeresurus stejnegeri stejnegeri or Protobothrops mucrosquamatus, those envenomed by Naja atra or Bungarus multicinctus multicinctus were more likely to occur in central (adjusted odds ratio [aOR] = 2.6, p < 0.0001) or southern (aOR = 3.2, p < 0.0001) Taiwan, but less frequently among agricultural workers (aOR = 0.6, p < 0.0001). The overall case-fatality rate was 0.11%.

CONCLUSIONS

Among Asian countries, Taiwan had low incidence and case-fatality rates of SBE. Risk factors included male gender, old age, summer season, being in eastern Taiwan, and being an agricultural worker. Differences of the epidemiological findings between snake species should be focused on when developing strategies for snakebite prevention.

Omics Technologies for Profiling Toxin Diversity and Evolution in Snake Venom: Impacts on the Discovery of Therapeutic and Diagnostic Agents

Snake venoms are primarily composed of proteins and peptides, and these toxins have developed high selectivity to their biological targets. This makes venoms interesting for exploration into protein evolution and structure-function relationships. A single venom protein superfamily can exhibit a variety of pharmacological effects; these variations in activity originate from differences in functional sites, domains, posttranslational modifications, and the formations of toxin complexes. In this review, we discuss examples of how the major venom protein superfamilies have diversified, as well as how newer technologies in the omics fields, such as genomics, transcriptomics, and proteomics, can be used to characterize both known and unknown toxins.Because toxins are bioactive molecules with a rich diversity of activities, they can be useful as therapeutic and diagnostic agents, and successful examples of toxin applications in these areas are also reviewed. With the current rapid pace of technology, snake venom research and its applications will only continue to expand.

Evaluating the Performance of De Novo Assembly Methods for Venom-Gland Transcriptomics

Venom-gland transcriptomics is a key tool in the study of the evolution, ecology, function, and pharmacology of animal venoms. In particular, gene-expression variation and coding sequences gained through transcriptomics provide key information for explaining functional venom variation over both ecological and evolutionary timescales. The accuracy and usefulness of inferences made through transcriptomics, however, is limited by the accuracy of the transcriptome assembly, which is a bioinformatic problem with several possible solutions. Several methods have been employed to assemble venom-gland transcriptomes, with the Trinity assembler being the most commonly applied among them. Although previous evidence of variation in performance among assembly software exists, particularly regarding recovery of difficult-to-assemble multigene families such as snake venom metalloproteinases, much work to date still employs a single assembly method. We evaluated the performance of several commonly used de novo assembly methods for the recovery of both nontoxin transcripts and complete, high-quality venom-gene transcripts across eleven snake and four scorpion transcriptomes. We varied k-mer sizes used by some assemblers to evaluate the impact of k-mer length on transcript recovery. We showed that the recovery of nontoxin transcripts and toxin transcripts is best accomplished through different assembly software, with SDT at smaller k-mer lengths and Trinity being best for nontoxin recovery and a combination of SeqMan NGen and a seed-and-extend approach implemented in Extender as the best means of recovering a complete set of toxin transcripts. In particular, Extender was the only means tested capable of assembling multiple isoforms of the diverse snake venom metalloproteinase family, while traditional approaches such as Trinity recovered at most one metalloproteinase transcript. Our work demonstrated that traditional metrics of assembly performance are not predictive of performance in the recovery of complete and high quality toxin genes. Instead, effective venom-gland transcriptomic studies should combine and quality-filter the results of several assemblers with varying algorithmic strategies.

Computational Studies of Snake Venom Toxins

Most snake venom toxins are proteins, and participate to envenomation through a diverse array of bioactivities, such as bleeding, inflammation, and pain, cytotoxic, cardiotoxic or neurotoxic effects. The venom of a single snake species contains hundreds of toxins, and the venoms of the 725 species of venomous snakes represent a large pool of potentially bioactive proteins. Despite considerable discovery efforts, most of the snake venom toxins are still uncharacterized. Modern bioinformatics tools have been recently developed to mine snake venoms, helping focus experimental research on the most potentially interesting toxins. Some computational techniques predict toxin molecular targets, and the binding mode to these targets. This review gives an overview of current knowledge on the ~2200 sequences, and more than 400 three-dimensional structures of snake toxins deposited in public repositories, as well as of molecular modeling studies of the interaction between these toxins and their molecular targets. We also describe how modern bioinformatics have been used to study the snake venom protein phospholipase A2, the small basic myotoxin Crotamine, and the three-finger peptide Mambalgin.

Targeting Metastasis with Snake Toxins: Molecular Mechanisms

Metastasis involves the migration of cancer cells from a primary tumor to invade and establish secondary tumors in distant organs, and it is the main cause for cancer-related deaths. Currently, the conventional cytostatic drugs target the proliferation of malignant cells, being ineffective in metastatic disease. This highlights the need to find new anti-metastatic drugs. Toxins isolated from snake venoms are a natural source of potentially useful molecular scaffolds to obtain agents with anti-migratory and anti-invasive effects in cancer cells. While there is greater evidence concerning the mechanisms of cell death induction of several snake toxin classes on cancer cells; only a reduced number of toxin classes have been reported on (i.e., disintegrins/disintegrin-like proteins, C-type lectin-like proteins, C-type lectins, serinproteases, cardiotoxins, snake venom cystatins) as inhibitors of adhesion, migration, and invasion of cancer cells. Here, we discuss the anti-metastatic mechanisms of snake toxins, distinguishing three targets, which involve (1) inhibition of extracellular matrix components-dependent adhesion and migration, (2) inhibition of epithelial-mesenchymal transition, and (3) inhibition of migration by alterations in the actin/cytoskeleton network.

Colubrid Venom Composition: An -Omics Perspective

Snake venoms have been subjected to increasingly sensitive analyses for well over 100 years, but most research has been restricted to front-fanged snakes, which actually represent a relatively small proportion of extant species of advanced snakes. Because rear-fanged snakes are a diverse and distinct radiation of the advanced snakes, understanding venom composition among "colubrids" is critical to understanding the evolution of venom among snakes. Here we review the state of knowledge concerning rear-fanged snake venom composition, emphasizing those toxins for which protein or transcript sequences are available. We have also added new transcriptome-based data on venoms of three species of rear-fanged snakes. Based on this compilation, it is apparent that several components, including cysteine-rich secretory proteins (CRiSPs), C-type lectins (CTLs), CTLs-like proteins and snake venom metalloproteinases (SVMPs), are broadly distributed among "colubrid" venoms, while others, notably three-finger toxins (3FTxs), appear nearly restricted to the Colubridae (sensu stricto). Some putative new toxins, such as snake venom matrix metalloproteinases, are in fact present in several colubrid venoms, while others are only transcribed, at lower levels. This work provides insights into the evolution of these toxin classes, but because only a small number of species have been explored, generalizations are still rather limited. It is likely that new venom protein families await discovery, particularly among those species with highly specialized diets.