The rise of genomics in snake venom research: recent advances and future perspectives

Proteogenomics-guided functional venomics resolves the toxin arsenal and activity of Deinagkistrodon acutus venom

Snakebite primarily impacts rural communities of Africa, Asia, and Latin America. The sharp-nosed viper (Deinagkistrodon acutus) is among the snakes of highest medical importance in Asia. Despite various studies on its venom using modern venomics techniques, a comprehensive understanding of composition and function of this species' venom remains lacking. We combined proteogenomics with extensive bioactivity profiling to present the first genome-level catalogue of D. acutus venom proteins and their exochemistry. Our analysis identified an unusually simple venom containing 45 components from 20 distinct protein families. Relative toxin abundances indicate that C-type lectin and C-type lectin-related protein (CTL), snake venom metalloproteinase (svMP), snake venom serine protease (svSP), and phospholipase A2 (PLA2) constitute 90 % of the venom. Bioassays targeting key aspects of viperid envenomation showed considerable concentration-dependent cytotoxicity, particularly in kidney and lung cells, and potent protease and PLA2 activity. Factor Xa and thrombin activities were minor, and no plasmin activity was observed. Effects on haemolysis, intracellular calcium (Ca2+) release, and nitric oxide (NO) synthesis were negligible. Our analysis provides the first holistic genome-based overview of the toxin arsenal of D. acutus, predicting the molecular and functional basis of its life-threatening effects, and opens novel avenues for treating envenomation by this highly dangerous snake.

Where the "ruber" Meets the Road: Using the Genome of the Red Diamond Rattlesnake to Unravel the Evolutionary Processes Driving Venom Evolution

Understanding the proximate and ultimate causes of phenotypic variation is fundamental in evolutionary research, as such variation provides the substrate for selection to act upon. Although trait variation can arise due to selection, the importance of neutral processes is sometimes understudied. We presented the first reference-quality genome of the Red Diamond Rattlesnake (Crotalus ruber) and used range-wide 'omic data to estimate the degree to which neutral and adaptive evolutionary processes shaped venom evolution. We characterized population structure and found substantial genetic differentiation across two populations, each with distinct demographic histories. We identified significant differentiation in venom expression across age classes with substantially reduced but discernible differentiation across populations. We then used conditional redundancy analysis to test whether venom expression variation was best predicted by neutral divergence patterns or geographically variable (a)biotic factors. Snake size was the most significant predictor of venom variation, with environment, prey availability, and neutral sequence variation also identified as significant factors, though to a lesser degree. By directly including neutrality in the model, our results confidently highlight the predominant, yet not singular, role of life history in shaping venom evolution.

Quantification of snake venom proteomes by mass spectrometry-considerations and perspectives

The advent of soft ionization mass spectrometry-based proteomics in the 1990s led to the development of a new dimension in biology that conceptually allows for the integral analysis of whole proteomes. This transition from a reductionist to a global-integrative approach is conditioned to the capability of proteomic platforms to generate and analyze complete qualitative and quantitative proteomics data. Paradoxically, the underlying analytical technique, molecular mass spectrometry, is inherently nonquantitative. The turn of the century witnessed the development of analytical strategies to endow proteomics with the ability to quantify proteomes of model organisms in the sense of "an organism for which comprehensive molecular (genomic and/or transcriptomic) resources are available." This essay presents an overview of the strategies and the lights and shadows of the most popular quantification methods highlighting the common misuse of label-free approaches developed for model species' when applied to quantify the individual components of proteomes of nonmodel species (In this essay we use the term "non-model" organisms for species lacking comprehensive molecular (genomic and/or transcriptomic) resources, a circumstance that, as we detail in this review-essay, conditions the quantification of their proteomes.). We also point out the opportunity of combining elemental and molecular mass spectrometry systems into a hybrid instrumental configuration for the parallel identification and absolute quantification of venom proteomes. The successful application of this novel mass spectrometry configuration in snake venomics represents a proof-of-concept for a broader and more routine application of hybrid elemental/molecular mass spectrometry setups in other areas of the proteomics field, such as phosphoproteomics, metallomics, and in general in any biological process where a heteroatom (i.e., any atom other than C, H, O, N) forms integral part of its mechanism.

V-ToCs (Venom Toxin Clustering): A tool for the investigation of sequence and structure similarities in snake venom toxins

Recently, there has been a major push toward the development of next-generation treatments against snakebite envenoming. However, unlike current antivenoms that rely on animal-derived polyclonal antibodies, most of these novel approaches are reliant on an in-depth understanding of the over 2000 known snake venom toxins. Indeed, by identifying similarities (i.e., conserved epitopes) across these different toxins, it is possible to design cross-reactive treatments, such as broadly-neutralising antibodies, that target these similarities. Therefore, in this project, we built an automated pipeline that generates sequence and structural distance matrices and homology trees across all available snake venom toxin sequences and structures. To facilitate analysis, we also developed a user-friendly and high-throughput visualisation tool, coined "Venom TOxin CluStering" (V-ToCs). This tool allows researchers to easily investigate sequence and structure patterns in snake venom toxins for a wide array of purposes, such as elucidating toxin evolution, and will also hopefully help guide the discovery and development of increasingly broadly-neutralising antivenoms in the near future.

Dabsylated Bradykinin Is Cleaved by Snake Venom Proteases from Echis ocellatus

The vasoactive peptide bradykinin (BK) is an important member of the renin-angiotensin system. Its discovery is tightly interwoven with snake venom research, because it was first detected in plasma following the addition of viper venom. While the fact that venoms liberate BK from a serum globulin fraction is well described, its destruction by the venom has largely gone unnoticed. Here, BK was found to be cleaved by snake venom metalloproteinases in the venom of Echis ocellatus, one of the deadliest snakes, which degraded its dabsylated form (DBK) in a few minutes after Pro7 (RPPGFSP↓FR). This is a common cleavage site for several mammalian proteases such as ACE, but is not typical for matrix metalloproteinases. Residual protease activity < 5% after addition of EDTA indicated that DBK is also cleaved by serine proteases to a minor extent. Mass spectrometry-based protein analysis provided spectral proof for several peptides of zinc metalloproteinase-disintegrin-like Eoc1, disintegrin EO4A, and three serine proteases in the venom.

Hidden in the sand: Phylogenomics unravel an unexpected evolutionary history for the desert-adapted vipers of the genus Cerastes

The desert vipers of the genus Cerastes are a small clade of medically important venomous snakes within the family Viperidae. According to published morphological and molecular studies, the group is comprised by four species: two morphologically similar and phylogenetically sister taxa, the African horned viper (Cerastes cerastes) and the Arabian horned viper (Cerastes gasperettii); a more distantly related species, the Saharan sand viper (Cerastes vipera), and the enigmatic Böhme's sand viper (Cerastes boehmei), only known from a single specimen in captivity allegedly captured in Central Tunisia. In this study, we sequenced one mitochondrial marker (COI) as well as genome-wide data (ddRAD sequencing) from 28 and 41 samples, respectively, covering the entire distribution range of the genus to explore the population genomics, phylogenomic relationships and introgression patterns within the genus Cerastes. Additionally, and to provide insights into the mode of diversification of the group, we carried out niche overlap analyses considering climatic and habitat variables. Both nuclear phylogenomic reconstructions and population structure analyses have unveiled an unexpected evolutionary history for the genus Cerastes, which sharply contradicts the morphological similarities and previously published mitochondrial approaches. Cerastes cerastes and C. vipera are recovered as sister taxa whilst C. gasperettii is a sister taxon to the clade formed by these two species. We found a relatively high niche overlap (OI > 0.7) in both climatic and habitat variables between C. cerastes and C. vipera, contradicting a potential scenario of sympatric speciation. These results are in line with the introgression found between the northwestern African populations of C. cerastes and C. vipera. Finally, our genomic data confirms the existence of a lineage of C. cerastes in Arabia. All these results highlight the importance of genome-wide data over few genetic markers to study the evolutionary history of species.

The genome assembly and annotation of the white-lipped tree pit viper Trimeresurus albolabris

Trimeresurus albolabris, also known as the white-lipped pit viper or white-lipped tree viper, is a highly venomous snake distributed across Southeast Asia and the cause of many snakebite cases. In this study, we report the first whole genome assembly of T. albolabris obtained with next-generation sequencing from a specimen collected in Mengzi, Yunnan, China. After genome sequencing and assembly, the genome of this male T. albolabris individual was 1.51 Gb in length and included 38.42% repeat-element content. Using this genome, 21,695 genes were identified, and 99.17% of genes could be annotated using gene functional databases. Our genome assembly and annotation process was validated using a phylogenetic tree, which included six species and focused on single-copy genes of nuclear genomes. This research will contribute to future studies on Trimeresurus biology and the genetic basis of snake venom.

Comparative Evolutionary Genomics in Insects

Genome sequencing quality, in terms of both read length and accuracy, is constantly improving. By combining long-read sequencing technologies with various scaffolding techniques, chromosome-level genome assemblies are now achievable at an affordable price for non-model organisms. Insects represent an exciting taxon for studying the genomic underpinnings of evolutionary innovations, due to ancient origins, immense species-richness, and broad phenotypic diversity. Here we summarize some of the most important methods for carrying out a comparative genomics study on insects. We describe available tools and offer concrete tips on all stages of such an endeavor from DNA extraction through genome sequencing, annotation, and several evolutionary analyses. Along the way we describe important insect-specific aspects, such as DNA extraction difficulties or gene families that are particularly difficult to annotate, and offer solutions. We describe results from several examples of comparative genomics analyses on insects to illustrate the fascinating questions that can now be addressed in this new age of genomics research.

Genome assembly and annotation of the Brown-Spotted Pit viper Protobothrops mucrosquamatus

The Brown-Spotted Pit viper (Protobothrops mucrosquamatus), also known as the Chinese habu, is a widespread and highly venomous snake distributed from Northeastern India to Eastern China. Genomics research can contribute to our understanding of venom components and natural selection in vipers. Here, we collected, sequenced and assembled the genome of a male P. mucrosquamatus individual from China. We generated a highly continuous reference genome, with a length of 1.53 Gb and 41.18% of repeat elements content. Using this genome, we identified 24,799 genes, 97.97% of which could be annotated. We verified the validity of our genome assembly and annotation process by generating a phylogenetic tree based on the nuclear genome single-copy genes of six other reptile species. The results of our research will contribute to future studies on Protobothrops biology and the genetic basis of snake venom.

Slowly Making Sense: A Review of the Two-Step Venom System within Slow (Nycticebus spp.) and Pygmy Lorises (Xanthonycticebus spp.)

Since the early 2000s, studies of the evolution of venom within animals have rapidly expanded, offering new revelations on the origins and development of venom within various species. The venomous mammals represent excellent opportunities to study venom evolution due to the varying functional usages, the unusual distribution of venom across unrelated mammals and the diverse variety of delivery systems. A group of mammals that excellently represents a combination of these traits are the slow (Nycticebus spp.) and pygmy lorises (Xanthonycticebus spp.) of south-east Asia, which possess the only confirmed two-step venom system. These taxa also present one of the most intriguing mixes of toxic symptoms (cytotoxicity and immunotoxicity) and functional usages (intraspecific competition and ectoparasitic defence) seen in extant animals. We still lack many pieces of the puzzle in understanding how this venom system works, why it evolved what is involved in the venom system and what triggers the toxic components to work. Here, we review available data building upon a decade of research on this topic, focusing especially on why and how this venom system may have evolved. We discuss that research now suggests that venom in slow lorises has a sophisticated set of multiple uses in both intraspecific competition and the potential to disrupt the immune system of targets; we suggest that an exudate diet reveals several toxic plants consumed by slow and pygmy lorises that could be sequestered into their venom and which may help heal venomous bite wounds; we provide the most up-to-date visual model of the brachial gland exudate secretion protein (BGEsp); and we discuss research on a complement component 1r (C1R) protein in saliva that may solve the mystery of what activates the toxicity of slow and pygmy loris venom. We conclude that the slow and pygmy lorises possess amongst the most complex venom system in extant animals, and while we have still a lot more to understand about their venom system, we are close to a breakthrough, particularly with current technological advances.

The genome assembly and annotation of the many-banded krait, Bungarus multicinctus

Snakes are a vital component of wildlife resources and are widely distributed across the globe. The many-banded krait Bungarus multicinctus is a highly venomous snake found across Southern Asia and central and southern China. Snakes are an ancient reptile group, and their genomes can provide important clues for understanding the evolutionary history of reptiles. Additionally, genomic resources play a crucial role in comprehending the evolution of all species. However, snake genomic resources are still scarce. Here, we present a highly contiguous genome of B. multicinctus with a size of 1.51 Gb. The genome contains a repeat content of 40.15%, with a total length exceeding 620 Mb. Additionally, we annotated a total of 24,869 functional genes. This research is of great significance for comprehending the evolution of B. multicinctus and provides genomic information on the genes involved in venom gland functions.

Unlocking the potential of snake venom-based molecules against the malaria, Chagas disease, and leishmaniasis triad

Malaria, leishmaniasis and Chagas disease are vector-borne protozoal infections with a disproportionately high impact on the most fragile societies in the world, and despite malaria-focused research gained momentum in the past two decades, both trypanosomiases and leishmaniases remain neglected tropical diseases. Affordable effective drugs remain the mainstay of tackling this burden, but toxicicty, inneficiency against later stage disease, and drug resistance issues are serious shortcomings. One strategy to overcome these hurdles is to get new therapeutics or inspiration in nature. Indeed, snake venoms have been recognized as valuable sources of biomacromolecules, like peptides and proteins, with antiprotozoal activity. This review highlights major snake venom components active against at least one of the three aforementioned diseases, which include phospholipases A2, metalloproteases, L-amino acid oxidases, lectins, and oligopeptides. The relevance of this repertoire of biomacromolecules and the bottlenecks in their clinical translation are discussed considering approaches that should increase the success rate in this arduous task. Overall, this review underlines how venom-derived biomacromolecules could lead to pioneering antiprotozoal treatments and how the drug landscape for neglected diseases may be revolutionized by a closer look at venoms. Further investigations on poorly studied venoms is needed and could add new therapeutics to the pipeline.

A single-chain variable fragment selected against a conformational epitope of a recombinantly produced snake toxin using phage display

Phage display technology is a powerful tool for selecting monoclonal antibodies against a diverse set of antigens. Within toxinology, however, it remains challenging to generate monoclonal antibodies against many animal toxins, as they are difficult to obtain from venom. Recombinant toxins have been proposed as a solution to overcome this challenge, but so far, few have been used as antigens to generate neutralizing antibodies. Here, we describe the recombinant expression of α-cobratoxin in E. coli and its successful application as an antigen in a phage display selection campaign. From this campaign, an scFv (single chain variable fragment) was isolated with similar binding affinity to a control scFv generated against the native toxin. The selected scFv recognizes a structural epitope, enabling it to inhibit the interaction between the acetylcholine receptor and the native toxin in vitro. This approach represents the first entirely in vitro antibody selection strategy for generating neutralizing monoclonal antibodies against a snake toxin.

High-throughput proteomics and in vitro functional characterization of the 26 medically most important elapids and vipers from sub-Saharan Africa

Venomous snakes are important parts of the ecosystem, and their behavior and evolution have been shaped by their surrounding environments over the eons. This is reflected in their venoms, which are typically highly adapted for their biological niche, including their diet and defense mechanisms for deterring predators. Sub-Saharan Africa is rich in venomous snake species, of which many are dangerous to humans due to the high toxicity of their venoms and their ability to effectively deliver large amounts of venom into their victims via their bite. In this study, the venoms of 26 of sub-Saharan Africa's medically most relevant elapid and viper species were subjected to parallelized toxicovenomics analysis. The analysis included venom proteomics and in vitro functional characterization of whole venom toxicities, enabling a robust comparison of venom profiles between species. The data presented here corroborate previous studies and provide biochemical details for the clinical manifestations observed in envenomings by the 26 snake species. Moreover, two new venom proteomes (Naja anchietae and Echis leucogaster) are presented here for the first time. Combined, the presented data can help shine light on snake venom evolutionary trends and possibly be used to further improve or develop novel antivenoms.

Polymerase chain reaction-based snake origin tracing in commercial venom crystals by targeting the mitochondrial D-loop

Snake venom is a valuable raw material for numerous therapeutic formulations because of its life-saving pharmacological potential. However, due to their high price, fake "snake venoms" have captured a significant portion of the global market, and there is currently no reliable reported DNA-based method available for quickly distinguishing between fakes and originals. Therefore, in this study, a set of newly designed snake-specific universal primers targeting mitochondrial D-loop fragments were employed to detect snake origins in commercial venom crystals by only simplex polymerase chain reaction analysis. Under the optimal thermal cycling conditions, only the 145-149 bp snake-specific mitochondrial D-loop fragments from pure and mixed backgrounds were amplified by the newly designed primers. Specificity was achieved by confirming no DNA amplification occurred in the DNA admixture of ten different chordates, and universality by individual DNA amplification of nine different snakes. The primers that efficiently amplified the minimum mitochondrial DNA contained in a total of 10^-2 ng in a 10.0 μl reaction were also successfully able to detect the snake origin in commercial cobra venom crystals. These findings suggest that the newly designed primers can be used to differentiate the original and fake commercial snake venom crystals in order to achieve the highest standards of snake venom-based medications through amplifying the snake-specific mitochondrial D-loop fragments.