Evolution of CRISPs associated with toxicoferan-reptilian venom and mammalian reproduction

Dietary Breadth Predicts Toxin Expression Complexity in the Venoms of North American Gartersnakes

Selection on heritable phenotypic variation has played a prominent role in shaping the remarkable adaptations found across the Tree of Life. Complex ecological traits, such as snake venoms, are thought to be the products of selection because they directly link to fitness and survival. Snake venom increases the efficiency of prey capture and processing and is thus likely under intense selection. While many studies of snake venom have investigated the relationship between venom and diet, they have primarily focused on medically relevant front-fanged snakes. However, recent work has suggested that many non-front fanged snakes also rely on venom for subduing prey, despite having reduced toxicity toward humans. Here, we set out to uncover variation in toxin-producing genes, along with the ecological and evolutionary pressures impacting snake venom characteristics in the North American gartersnakes (Squamata: Natricidae: Thamnophis), a model group of non-front-fanged snakes. We annotated and analyzed Duvernoy's venom gland transcriptomes from 16 species representing all the major lineages within Thamnophis. We then generated measures of complexity of both toxins and dietary breadth. We found strong correlations between the complexity of toxin gene expression and phylogenetic diversity of diet, but no relationship between the complexity of the genetic makeup of the transcriptomes (allelic or sequence variation) and diet complexities. We also found phylogenetic signal associated with venom complexity, suggesting some influence of ancestry on venom characteristics. We suggest that, in non-front-fanged snakes, expression of toxins rather than sequence complexity is under strong selection by dietary diversity. These findings contradict similar studies from front-fanged snakes where increased transcriptomic complexity varies positively with dietary diversity, exposing a potential novel relationship between a complex phenotype-toxin expression-and its selective pressures-diet.

Non-transgenic, PAMAM co-delivery DNA of interactive proteins NbCRVP and NbCalB endows Nicotiana benthamiana with a stronger antiviral effect to RNA viruses

BACKGROUND

Viral diseases continue to pose a major threat to the world's commercial crops. The in-depth exploration and efficient utilization of resistance proteins have become crucial strategies for their control. However, current delivery methods for introducing foreign DNA suffer from host range limitations, low transformation efficiencies, tissue damage, or unavoidable DNA integration into the host genome. The nanocarriers provides a convenient channel for the DNA delivery and functional utilization of disease-resistant proteins.

RESULTS

In this research, we identified a cysteine-rich venom protein (NbCRVP) in Nicotiana benthamiana for the first time. Virus-induced gene silencing and transient overexpression clarified that NbCRVP could inhibit the infection of tobacco mosaic virus, potato virus Y, and cucumber mosaic virus, making it a broad-spectrum antiviral protein. Yeast two-hybrid assay, co-immunoprecipitation, and bimolecular fluorescence complementation revealed that calcium-dependent lipid-binding (CaLB domain) family protein (NbCalB) interacted with NbCRVP to assist NbCRVP playing a stronger antiviral effect. Here, we demonstrated for the first time the efficient co-delivery of DNA expressing NbCRVP and NbCalB into plants using poly(amidoamine) (PAMAM) nanocarriers, achieving stronger broad-spectrum antiviral effects.

CONCLUSIONS

Our work presents a tool for species-independent transfer of two interacting protein DNA into plant cells in a specific ratio for enhanced antiviral effect without transgenic integration, which further demonstrated new strategies for nanocarrier-mediated DNA delivery of disease-resistant proteins.

Varying Modes of Selection Among Toxin Families in the Venoms of the Giant Desert Hairy Scorpions (Hadrurus)

Venoms are primarily believed to evolve under strong diversifying selection resulting from persistent coevolution between predator and prey. Recent research has challenged this hypothesis, proposing that venoms from younger venomous lineages (e.g., snakes and cone snails) are governed predominantly by diversifying selection, while venoms from older venomous lineages (e.g., centipedes, scorpions, and spiders) are under stronger purifying selection. However, most research in older lineages has tested selection at more diverse phylogenetic scales. Although these tests are important for evaluating broad macroevolutionary trends underlying venom evolution, they are less equipped to detect species-level evolutionary trends, which likely have large impacts on venom variation seen at more diverse phylogenetic scales. To test for selection among closely related species from an older venomous lineage, we generated high-throughput venom-gland transcriptomes and venom proteomes for four populations of Giant Desert Hairy Scorpions (Hadrurus), including three Hadrurus arizonensis populations and one Hadrurus spadix population. We detected significant episodic and pervasive diversifying selection across a highly abundant toxin family that likely has a major role in venom function ([Formula: see text]KTxs), providing a contrast to the stronger purifying selection identified from other studies on scorpion venoms. Conversely, we detected weak episodic diversifying and/or stronger purifying selection in four toxin families (non-disulfide bridged peptides, phospholipase A2s, scorpine-like antimicrobial peptides, and serine proteases), most of which were less abundant and likely have ancillary functional roles. Finally, although we detected several major toxin families at disproportionate transcriptomic and/or proteomic abundances, we did not identify significant sex-based variation in Hadrurus venoms.

Venom comparisons of endemic and micro-endemic speckled rattlesnakes Crotalus mitchellii, C. polisi and C. thalassoporus from Baja California Peninsula

A high diversity of rattlesnake species can be found in the Baja California peninsula and the island of the Gulf of California, nevertheless, their venom has been poorly evaluated. The aim of this work was to present the first characterization of endemic Crotalus mitchellii, micro endemic C. polisi and C. thalassoporus venoms. All samples provoke human plasma coagulation showing doses in the rank of 2.3-41.0 μg and also produce rapid hydrolysis of the alpha chain of bovine fibrinogen while the beta chain is attacked at larger incubation periods by C. polisi and especially by C. thalassoporus. Phospholipase activity ranging from 23.2 to 173.8 U/mg. The venoms of C. thalassoporus and C. polisi show very high hemorrhagic activity (from 0.03 to 0.31 μg). A total of 130 toxin-related proteins were identified and classified into ten families. Crotalus mitchellii venom was characterized by high abundance of crotoxin-like and other phospholipase proteins (34.5%) and serine proteinases (29.8%). Crotalus polisi showed a similar proportion of metalloproteinases (34%) and serine proteinases (22.8%) components with important contribution of C-type lectins (14.3%) and CRiSP (14.0%) proteins. Venom of C. thalassoporus is dominated by metalloproteases that amount to more than 66% of total toxin proteins. These results provide a foundation for comprehending the biological, ecological and evolutionary significance of venom composition of speckled rattlesnake from the Baja California peninsula.

The deep-rooted origin of disulfide-rich spider venom toxins

Spider venoms are a complex concoction of enzymes, polyamines, inorganic salts, and disulfide-rich peptides (DRPs). Although DRPs are widely distributed and abundant, their bevolutionary origin has remained elusive. This knowledge gap stems from the extensive molecular divergence of DRPs and a lack of sequence and structural data from diverse lineages. By evaluating DRPs under a comprehensive phylogenetic, structural and evolutionary framework, we have not only identified 78 novel spider toxin superfamilies but also provided the first evidence for their common origin. We trace the origin of these toxin superfamilies to a primordial knot - which we name 'Adi Shakti', after the creator of the Universe according to Hindu mythology - 375 MYA in the common ancestor of Araneomorphae and Mygalomorphae. As the lineages under evaluation constitute nearly 60% of extant spiders, our findings provide fascinating insights into the early evolution and diversification of the spider venom arsenal. Reliance on a single molecular toxin scaffold by nearly all spiders is in complete contrast to most other venomous animals that have recruited into their venoms diverse toxins with independent origins. By comparatively evaluating the molecular evolutionary histories of araneomorph and mygalomorph spider venom toxins, we highlight their contrasting evolutionary diversification rates. Our results also suggest that venom deployment (e.g. prey capture or self-defense) influences evolutionary diversification of DRP toxin superfamilies.

Omics and Remote Homology Integration to Decipher Protein Functionality

In the recent years, several "omics" technologies based on specific biomolecules (from DNA, RNA, proteins, or metabolites) have won growing importance in the scientific field. Despite each omics possess their own laboratorial protocols, they share a background of bioinformatic tools for data integration and analysis. A recent subset of bioinformatic tools, based on available templates or remote homology protocols, allow computational fast and high-accuracy prediction of protein structures. The quickly predict of actually unsolved protein structures, together with late omics findings allow a boost of scientific advances in multiple fields such as cancer, longevity, immunity, mitochondrial function, toxicology, drug design, biosensors, and recombinant protein engineering. In this chapter, we assessed methodological approaches for the integration of omics and remote homology inferences to decipher protein functionality, opening the door to the next era of biological knowledge.

The Fast and the Furriest: Investigating the Rate of Selection on Mammalian Toxins

The evolution of venom and the selection pressures that act on toxins have been increasingly researched within toxinology in the last two decades, in part due to the exceptionally high rates of diversifying selection observed in animal toxins. In 2015, Sungar and Moran proposed the 'two-speed' model of toxin evolution linking evolutionary age of a group to the rates of selection acting on toxins but due to a lack of data, mammals were not included as less than 30 species of venomous mammal have been recorded, represented by elusive species which produce small amounts of venom. Due to advances in genomics and transcriptomics, the availability of toxin sequences from venomous mammals has been increasing. Using branch- and site-specific selection models, we present the rates of both episodic and pervasive selection acting upon venomous mammal toxins as a group for the first time. We identified seven toxin groups present within venomous mammals, representing Chiroptera, Eulipotyphla and Monotremata: KLK1, Plasminogen Activator, Desmallipins, PACAP, CRiSP, Kunitz Domain One and Kunitz Domain Two. All but one group (KLK1) was identified by our results to be evolving under both episodic and pervasive diversifying selection with four toxin groups having sites that were implicated in the fitness of the animal by TreeSAAP (Selection on Amino Acid Properties). Our results suggest that venomous mammal ecology, behaviour or genomic evolution are the main drivers of selection, although evolutionary age may still be a factor. Our conclusion from these results indicates that mammalian toxins are following the two-speed model of selection, evolving predominately under diversifying selection, fitting in with other younger venomous taxa like snakes and cone snails-with high amounts of accumulating mutations, leading to more novel adaptions in their toxins.

Contextual Constraints: Dynamic Evolution of Snake Venom Phospholipase A2

Venom is a dynamic trait that has contributed to the success of numerous organismal lineages. Predominantly composed of proteins, these complex cocktails are deployed for predation and/or self-defence. Many non-toxic physiological proteins have been convergently and recurrently recruited by venomous animals into their toxin arsenal. Phospholipase A₂ (PLA₂) is one such protein and features in the venoms of many organisms across the animal kingdom, including snakes of the families Elapidae and Viperidae. Understanding the evolutionary history of this superfamily would therefore provide insight into the origin and diversification of venom toxins and the evolution of novelty more broadly. The literature is replete with studies that have identified diversifying selection as the sole influence on PLA₂ evolution. However, these studies have largely neglected the structural/functional constraints on PLA₂s, and the ecology and evolutionary histories of the diverse snake lineages that produce them. By considering these crucial factors and employing evolutionary analyses integrated with a schema for the classification of PLA₂s, we uncovered lineage-specific differences in selection regimes. Thus, our work provides novel insights into the evolution of this major snake venom toxin superfamily and underscores the importance of considering the influence of evolutionary and ecological contexts on molecular evolution.

Qualitative Analysis of Proteins in Two Snake Venoms, Gloydius Blomhoffii and Agkistrodon Acutus

Objectives

Snake venom is a complex mixture of various pharmacologically active substances, such as small proteins, peptides, and organic and mineral components. This paper aims to identify and analyse the proteins in common venomous snakes, such as Gloydius blomhoffii (G. blomhoffii) and Agkistrodon acutus (A. acutus), in Korea.

Methods

We used mass spectrometry, electrophoresis, N-terminal sequencing and in-gel digestion to analyse the proteins in these two snake venoms.

Results

We identified eight proteins in G. blomhoffii venom and four proteins in A. acutus venom. The proteins detected in G. blomhoffii and A. acutus venoms were phospholipase A2, snake venom metalloproteinase and cysteine-rich secretory protein. Snake C-type lectin (snaclec) was unique to A. acutus venom.

Conclusion

These data will contribute to the current knowledge of proteins present in the venoms of viper snakes and provide useful information for investigating their therapeutic potential.

Chromosomal-level genome of a sheet-web spider provides insight into the composition and evolution of venom

Spiders are the most abundant venomous predators in the world. Previous research related to spider venom has mostly relied on transcriptomes and proteomes, with only a few high-quality genomes available. This is far from consistent with the species diversity of spiders. In this study, we constructed a high-quality chromosome-level genome assembly of Hylyphantes graminicola, which contained 13 chromosomes, with a genome length of 931.68 Mb and scaffold N50 of 77.07 Mb. Integrating genome, transcriptome, and proteome profiling, we identified a total of 59 coding genes among nine toxin gene families. Among them, Group 7 allergen (ALL7) protein was reported in spider venom for the first time. Its coding genes had a predicted signal peptide and maintained high expression levels in the venom, suggesting that ALL7 plays an important role in venom and maybe is a type of newly discovered venom toxin in the spider. By implementing comparative genomics, we found a similar gene number of main toxin gene families in spiders and the scorpion genome with conservative evolutionary rates, indicating that these toxin genes could be an ancient (~400 million years) and a conserved "basic toolkit" for spiders and scorpions to perform primary defense functions. Obtaining high-quality chromosome-level genomes from spiders not only facilitates venom research and toxin resource application, but also can improve comparative genomic analysis in other important traits, like the evolution of silk or behavior.

Evaluation of Signaling Pathways Profiling in Human Dermal Endothelial Cells Treated by Snake Venom Cysteine-Rich Secretory Proteins (svCRiSPs) from North American Snakes Using Reverse Phase Protein Array (RPPA)

Cysteine-Rich Secretory Proteins (CRiSPs) are typically found in many snake venoms; however, the role that these toxins play in the pathophysiology of snakebites is still unclear. Herein, we compared the effects of snake venom CRiSPs (svCRiSPs) from the most medically important species of North American snakes on endothelial cell permeability and vascular permeability. We used reverse phase protein array (RPPA) to identify key signaling molecules on human dermal lymphatic (HDLECs) and blood (HDBECs) endothelial cells treated with svCRiSPs. The results showed that Css-CRiSP isolated from Crotalus scutulatus scutulatus and App-CRiSP from Agkistrodon piscivorus piscivorus are the most potent causes of increase vascular and endothelial permeability in comparison with other svCRiSPs used in this study. We examined the protein expression levels and their activated phosphorylation states in HDLECs and HDBECs induced by App-CRiSP and Css-CRiSP using RPPA. Interestingly, both App-CRiSP and Css-CRiSP induced caveolin-1 expression in HDBECs. We also found that stimulating HDBECs with Css-CRiSP and App-CRiSP significantly induced the phosphorylation of mTOR and Src, respectively. In HDLECs, Css-CRiSP significantly downregulated the expression of N-Cadherin and phospholipase C-gamma, while App-CRiSP significantly enhanced Akt and JNK phosphorylation. These results suggest that the increased endothelial permeability in HDLECs and HDBECs by Css-CRiSP and App-CRiSP may occur through different pathways.

Characterization of different oligomeric forms of CRISP2 in the perinuclear theca versus the fibrous tail structures of boar spermatozoa

Mammalian sperm carry a variety of highly condensed insoluble protein structures such as the perinuclear theca, the fibrous sheath and the outer dense fibers, which are essential to sperm function. We studied the role of cysteine rich secretory protein 2 (CRISP2); a known inducer of non-pathological protein amyloids, in pig sperm with a variety of techniques. CRISP2, which is synthesized during spermatogenesis, was localized by confocal immunofluorescent imaging in the tail and in the post-acrosomal region of the sperm head. High resolution localization by immunogold labeling electron microscopy (EM) of ultrathin cryosections revealed that CRISP2 was present in the perinuclear theca and neck region of the sperm head, as well as in the outer dense fibers and the fibrous sheath of the sperm tail. Interestingly, we found that under native, non-reducing conditions CRISP2 formed oligomers both in the tail and the head but with different molecular weights and different biochemical properties. The tail oligomers were insensitive to reducing conditions but nearly complete dissociated into monomers under 8 M urea treatment, while the head 250 kDa CRISP2 positive oligomer completely dissociated into CRISP2 monomers under reducing conditions. The head specific dissociation of CRISP2 oligomer is likely a result of the reduction of various sulfhydryl groups in the cysteine rich domain of this protein. The sperm head CRISP2 shared typical solubilization characteristics with other perinuclear theca proteins as was shown with sequential detergent and salt treatments. Thus, CRISP2 is likely to participate in the formation of functional protein complexes in both the sperm tail and sperm head, but with differing oligomeric organization and biochemical properties. Future studies will be devoted to the understand the role of CRISP2 in sperm protein complexes formation and how this contributes to the fertilization processes.

Proteo-Trancriptomic Analyses Reveal a Large Expansion of Metalloprotease-Like Proteins in Atypical Venom Vesicles of the Wasp Meteorus pulchricornis (Braconidae)

Meteorus pulchricornis (Ichneumonoidea, Braconidae) is an endoparasitoid wasp of lepidopteran caterpillars. Its parasitic success relies on vesicles (named M. pulchricornis Virus-Like Particles or MpVLPs) that are synthesized in the venom gland and injected into the parasitoid host along with the venom during oviposition. In order to define the content and understand the biogenesis of these atypical vesicles, we performed a transcriptome analysis of the venom gland and a proteomic analysis of the venom and purified MpVLPs. About half of the MpVLPs and soluble venom proteins identified were unknown and no similarity with any known viral sequence was found. However, MpVLPs contained a large number of proteins labelled as metalloproteinases while the most abundant protein family in the soluble venom was that of proteins containing the Domain of Unknown Function DUF-4803. The high number of these proteins identified suggests that a large expansion of these two protein families occurred in M. pulchricornis. Therefore, although the exact mechanism of MpVLPs formation remains to be elucidated, these vesicles appear to be “metalloproteinase bombs” that may have several physiological roles in the host including modifying the functions of its immune cells. The role of DUF4803 proteins, also present in the venom of other braconids, remains to be clarified.

Old World Vipers-A Review about Snake Venom Proteomics of Viperinae and Their Variations

Fine-tuned by millions of years of evolution, snake venoms have frightened but also fascinated humanity and nowadays they constitute potential resources for drug development, therapeutics and antivenoms. The continuous progress of mass spectrometry techniques and latest advances in proteomics workflows enabled toxinologists to decipher venoms by modern omics technologies, so-called ‘venomics’. A tremendous upsurge reporting on snake venom proteomes could be observed. Within this review we focus on the highly venomous and widely distributed subfamily of Viperinae (Serpentes: Viperidae). A detailed public literature database search was performed (2003–2020) and we extensively reviewed all compositional venom studies of the so-called Old-World Vipers. In total, 54 studies resulted in 89 venom proteomes. The Viperinae venoms are dominated by four major, four secondary, six minor and several rare toxin families and peptides, respectively. The multitude of different venomics approaches complicates the comparison of venom composition datasets and therefore we differentiated between non-quantitative and three groups of quantitative workflows. The resulting direct comparisons within these groups show remarkable differences on the intra- and interspecies level across genera with a focus on regional differences. In summary, the present compilation is the first comprehensive up-to-date database on Viperinae venom proteomes and differentiating between analytical methods and workflows.

An ancient, conserved gene regulatory network led to the rise of oral venom systems

Although oral venom systems are ecologically important characters, how they originated is still unclear. In this study, we show that oral venom systems likely originated from a gene regulatory network conserved across amniotes. This network, which we term the “metavenom network,” comprises over 3,000 housekeeping genes coexpressed with venom and play a role in protein folding and modification. Comparative transcriptomics revealed that the network is conserved between venom glands of snakes and salivary glands of mammals. This suggests that while these tissues have evolved different functions, they share a common regulatory core, that persisted since their common ancestor. We propose several evolutionary mechanisms that can utilize this common regulatory core to give rise to venomous animals from their nonvenomous ancestors.

Oral venom systems evolved multiple times in numerous vertebrates enabling the exploitation of unique predatory niches. Yet how and when they evolved remains poorly understood. Up to now, most research on venom evolution has focused strictly on the toxins. However, using toxins present in modern day animals to trace the origin of the venom system is difficult, since they tend to evolve rapidly, show complex patterns of expression, and were incorporated into the venom arsenal relatively recently. Here we focus on gene regulatory networks associated with the production of toxins in snakes, rather than the toxins themselves. We found that overall venom gland gene expression was surprisingly well conserved when compared to salivary glands of other amniotes. We characterized the “metavenom network,” a network of ∼3,000 nonsecreted housekeeping genes that are strongly coexpressed with the toxins, and are primarily involved in protein folding and modification. Conserved across amniotes, this network was coopted for venom evolution by exaptation of existing members and the recruitment of new toxin genes. For instance, starting from this common molecular foundation, Heloderma lizards, shrews, and solenodon, evolved venoms in parallel by overexpression of kallikreins, which were common in ancestral saliva and induce vasodilation when injected, causing circulatory shock. Derived venoms, such as those of snakes, incorporated novel toxins, though still rely on hypotension for prey immobilization. These similarities suggest repeated cooption of shared molecular machinery for the evolution of oral venom in mammals and reptiles, blurring the line between truly venomous animals and their ancestors.

Integrative multiomics analysis of Premolis semirufa caterpillar venom in the search for molecules leading to a joint disease

The joint disease called pararamosis is an occupational disease caused by accidental contact with bristles of the caterpillar Premolis semirufa. The chronic inflammatory process narrows the joint space and causes alterations in bone structure and cartilage degeneration, leading to joint stiffness. Aiming to determine the bristle components that could be responsible for this peculiar envenomation, in this work we have examined the toxin composition of the caterpillar bristles extract and compared it with the differentially expressed genes (DEGs) in synovial biopsies of patients affected with rheumatoid arthritis (RA) and osteoarthritis (OA). Among the proteins identified, 129 presented an average of 63% homology with human proteins and shared important conserved domains. Among the human homologous proteins, we identified seven DEGs upregulated in synovial biopsies from RA or OA patients using meta-analysis. This approach allowed us to suggest possible toxins from the pararama bristles that could be responsible for starting the joint disease observed in pararamosis. Moreover, the study of pararamosis, in turn, may lead to the discovery of specific pharmacological targets related to the early stages of articular diseases.

The functions of CAP superfamily proteins in mammalian fertility and disease

BACKGROUND

Members of the cysteine-rich secretory proteins (CRISPS), antigen 5 (Ag5) and pathogenesis-related 1 (Pr-1) (CAP) superfamily of proteins are found across the bacterial, fungal, plant and animal kingdoms. Although many CAP superfamily proteins remain poorly characterized, over the past decade evidence has accumulated, which provides insights into the functional roles of these proteins in various processes, including fertilization, immune defence and subversion, pathogen virulence, venom toxicology and cancer biology.

OBJECTIVE AND RATIONALE

The aim of this article is to summarize the current state of knowledge on CAP superfamily proteins in mammalian fertility, organismal homeostasis and disease pathogenesis.

SEARCH METHODS

The scientific literature search was undertaken via PubMed database on all articles published prior to November 2019. Search terms were based on following keywords: 'CAP superfamily', 'CRISP', 'Cysteine-rich secretory proteins', 'Antigen 5', 'Pathogenesis-related 1', 'male fertility', 'CAP and CTL domain containing', 'CRISPLD1', 'CRISPLD2', 'bacterial SCP', 'ion channel regulator', 'CatSper', 'PI15', 'PI16', 'CLEC', 'PRY proteins', 'ASP proteins', 'spermatogenesis', 'epididymal maturation', 'capacitation' and 'snake CRISP'. In addition to that, reference lists of primary and review article were reviewed for additional relevant publications.

OUTCOMES

In this review, we discuss the breadth of knowledge on CAP superfamily proteins with regards to their protein structure, biological functions and emerging significance in reproduction, health and disease. We discuss the evolution of CAP superfamily proteins from their otherwise unembellished prokaryotic predecessors into the multi-domain and neofunctionalized members found in eukaryotic organisms today. At least in part because of the rapid evolution of these proteins, many inconsistencies in nomenclature exist within the literature. As such, and in part through the use of a maximum likelihood phylogenetic analysis of the vertebrate CRISP subfamily, we have attempted to clarify this confusion, thus allowing for a comparison of orthologous protein function between species. This framework also allows the prediction of functional relevance between species based on sequence and structural conservation.

WIDER IMPLICATIONS

This review generates a picture of critical roles for CAP proteins in ion channel regulation, sterol and lipid binding and protease inhibition, and as ligands involved in the induction of multiple cellular processes.

Study of the venom proteome of Vipera ammodytes ammodytes (Linnaeus, 1758): A qualitative overview, biochemical and biological profiling

Vipera ammodytes (Va), is the European venomous snake of the greatest medical importance. We analyzed whole venom proteome of the subspecies V. ammodytes ammodytes (Vaa) from Serbia for the first time using the shotgun proteomics approach and identified 99 proteins belonging to four enzymatic families: serine protease (SVSPs), L-amino acid oxidase (LAAOs), metalloproteinases (SVMPs), group II phospholipase (PLA2s), and five nonenzymatic families: cysteine-rich secretory proteins (CRISPs), C-type lectins (snaclecs), growth factors -nerve (NGFs) and vascular endothelium (VEGFs), and Kunitz-type protease inhibitors (SPIs). Considerable enzymatic activity of LAAO, SVSPs, and SVMPs and a high acidic PLA2 activity was measured implying potential of Vaa to produce haemotoxic, myotoxic, neuro and cardiotoxic effects. Moreover, significant antimicrobial activity of Vaa venom against Gram-negative (Klebsiella pneumoniae, Pseudomonas aeruginosa) and Gram-positive bacteria (Staphylococcus aureus) was found. The crude venom shows considerable potential cytotoxic activity on the C6 and HL60 and a moderate level of potency on B16 cell lines. HeLa cells showed the same sensitivity, while DU 145 and PC-3 are less sensitive than as normal cell line. Our data demonstrated a high complexity of Vaa and considerable enzymatic, antibacterial and cytotoxic activity, implying a great medical potential of Vaa venom as a promising source for new antibacterial and cytostatic agents.

Causes and Consequences of Snake Venom Variation

Snake venoms are mixtures of toxins that vary extensively between and within snake species. This variability has serious consequences for the management of the world’s 1.8 million annual snakebite victims. Advances in ‘omic’ technologies have empowered toxinologists to comprehensively characterize snake venom compositions, unravel the molecular mechanisms that underpin venom variation, and elucidate the ensuing functional consequences. In this review, we describe how such mechanistic processes have resulted in suites of toxin isoforms that cause diverse pathologies in human snakebite victims and we detail how variation in venom composition can result in treatment failure. Finally, we outline current therapeutic approaches designed to circumvent venom variation and deliver next-generation treatments for the world’s most lethal neglected tropical disease.

Evolutionary analysis of genes coding for Cysteine-RIch Secretory Proteins (CRISPs) in mammals

BACKGROUND

Cysteine-RIch Secretory Proteins (CRISP) are expressed in the reproductive tract of mammalian males and are involved in fertilization and related processes. Due to their important role in sperm performance and sperm-egg interaction, these genes are likely to be exposed to strong selective pressures, including postcopulatory sexual selection and/or male-female coevolution. We here perform a comparative evolutionary analysis of Crisp genes in mammals. Currently, the nomenclature of CRISP genes is confusing, as a consequence of discrepancies between assignments of orthologs, particularly due to numbering of CRISP genes. This may generate problems when performing comparative evolutionary analyses of mammalian clades and species. To avoid such problems, we first carried out a study of possible orthologous relationships and putative origins of the known CRISP gene sequences. Furthermore, and with the aim to facilitate analyses, we here propose a different nomenclature for CRISP genes (EVAC1-4, "EVolutionarily-analyzed CRISP") to be used in an evolutionary context.

RESULTS

We found differing selective pressures among Crisp genes. CRISP1/4 (EVAC1) and CRISP2 (EVAC2) orthologs are found across eutherian mammals and seem to be conserved in general, but show signs of positive selection in primate CRISP1/4 (EVAC1). Rodent Crisp1 (Evac3a) seems to evolve under a comparatively more relaxed constraint with positive selection on codon sites. Finally, murine Crisp3 (Evac4), which appears to be specific to the genus Mus, shows signs of possible positive selection. We further provide evidence for sexual selection on the sequence of one of these genes (Crisp1/4) that, unlike others, is thought to be exclusively expressed in male reproductive tissues.

CONCLUSIONS

We found differing selective pressures among CRISP genes and sexual selection as a contributing factor in CRISP1/4 gene sequence evolution. Our evolutionary analysis of this unique set of genes contributes to a better understanding of Crisp function in particular and the influence of sexual selection on reproductive mechanisms in general.

Egg perivitelline fluid proteome of a freshwater snail: Insight into the transition from aquatic to terrestrial egg deposition

RATIONALE

Proteins from the egg perivitelline fluid (PVF) are assumed to play critical roles in embryonic development, but for many groups of animals their identities remain unknown. Identifying egg PVF proteins is a critical step towards understanding their functions including their roles in evolutionary transition in habitats.

METHODS

We applied proteomic and transcriptomic analysis to investigate the PVF proteome of the eggs of Pomacea diffusa, an aerial ovipositing freshwater snail in the family Ampullariidae. The PVF proteins were separated with the sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) method, and proteomic analysis was conducted using an LTQ Velos ion trap mass spectrometer coupled with liquid chromatography. Comparison of PVF proteomes and evolution analyses was performed between P. diffusa and other ampullariids.

RESULTS

In total, 32 egg PVF proteins were identified from P. diffusa. They were categorized as PV1-like subunits, immune-responsive proteins, protein degradation, signaling and binding, transcription and translation, metabolism, oxidation-reduction and proteins with unknown function. Interestingly, the proteome includes a calcium-binding protein important in forming the hard eggshell that enabled the terrestrial transition. However, it does not include PV2, a neurotoxic protein that was assumed to be present in all Pomacea species.

CONCLUSIONS

The PVF proteome data from P. diffusa can help us better understand the roles that reproductive proteins played during the transition from underwater to terrestrial egg deposition. Moreover, they could be useful in comparative studies of the terrestrialization in several groups of animals that occurred independently during their evolution.

Modahl C, Maenaka K, Aoki-Shioi N. Cysteine-Rich Secretory Proteins (CRISPs) From Venomous Snakes: An Overview of the Functional Diversity in A Large and Underappreciated Superfamily

The CAP protein superfamily (Cysteine-rich secretory proteins (CRISPs), Antigen 5 (Ag5), and Pathogenesis-related 1 (PR-1) proteins) is widely distributed, but for toxinologists, snake venom CRISPs are the most familiar members. Although CRISPs are found in the majority of venoms, very few of these proteins have been functionally characterized, but those that have been exhibit diverse activities. Snake venom CRISPs (svCRISPs) inhibit ion channels and the growth of new blood vessels (angiogenesis). They also increase vascular permeability and promote inflammatory responses (leukocyte and neutrophil infiltration). Interestingly, CRISPs in lamprey buccal gland secretions also manifest some of these activities, suggesting an evolutionarily conserved function. As we strive to better understand the functions that CRISPs serve in venoms, it is worth considering the broad range of CRISP physiological activities throughout the animal kingdom. In this review, we summarize those activities, known crystal structures and sequence alignments, and we discuss predicted functional sites. CRISPs may not be lethal or major components of venoms, but given their almost ubiquitous occurrence in venoms and the accelerated evolution of svCRISP genes, these venom proteins are likely to have functions worth investigating.

Evolution of digestive enzyme genes associated with dietary diversity of crabs

Crabs feed on a wide range of items and display diverse feeding strategies. The primary objective of this study was to investigate 10 digestive enzyme genes in representative crabs to provide insights into the genetic basis of feeding habits among crab functional groups. Crabs were classified into three groups based on their feeding habits: herbivores (HV), omnivores (OV), and carnivores (CV). To test whether crabs' feeding adaptations matched adaptive evolution of digestive enzyme genes, we examined the 10 digestive enzyme genes of 12 crab species based on hepatopancreas transcriptome data. Each of the digestive enzyme genes was compared to orthologous sequences using both nucleotide- (i.e., PAML and Datamonkey) and protein-level (i.e., TreeSAAP) approaches. Positive selection genes were detected in HV crabs (AMYA, APN, and MGAM) and CV crabs (APN, CPB, PNLIP, RISC, TRY, and XPD). Additionally, a series of positive selection sites were localized in important functional regions of these digestive enzyme genes. This is the first study to characterize the molecular basis of crabs' digestive enzyme genes based on functional feeding group. Our data suggest that HV crabs have evolved an enhanced digestion capacity for carbohydrates, and CV crabs have acquired digestion capacity for proteins and lipids.

Venomics of the asp viper Vipera aspis aspis from France

The asp viper Vipera aspis aspis is a venomous snake found in France, and despite its medical importance, the complete toxin repertoire produced is unknown. Here, we used a venomics approach to decipher the composition of its venom. Transcriptomic analysis revealed 80 venom-annotated sequences grouped into 16 gene families. Among the most represented toxins were snake venom metalloproteases (23%), phospholipases A2 (15%), serine proteases (13%), snake venom metalloprotease inhibitors (13%) and C-type lectins (12%). LC-MS of venoms revealed similar profiles regardless of the method of extraction (milking vs defensive bite). Proteomic analysis validated 57 venom-annotated transcriptomic sequences (>70%), including one for each of the 16 families, but also identified 7 sequences not initially annotated as venom proteins, including a serine protease, a disintegrin, a glutaminyl-peptide cyclotransferase, a proactivator polypeptide-like and 3 aminopeptidases. Interestingly, phospholipases A2 were the dominant proteins in the venom, among which included an ammodytoxin B-like sequence, which may explain the reported neurotoxicity following some asp viper envenomations. In total, 87 sequences were retrieved from the Vipera aspis aspis transcriptome and proteome, constituting a valuable resource that will help in understanding the toxinological basis of clinical signs of envenoming and for the mining of useful pharmacological compounds. BIOLOGICAL SIGNIFICANCE: The asp viper (Vipera aspis aspis) causes several hundred envenomations annually in France, including unusual cases with neurological signs, resulting in one death per year on average. Here, we performed a proteotranscriptomic analysis of V. a. aspis venom in order to provide a better understanding of its venom composition. We found that, as in other Vipera species, phospholipase A2 dominates in the venom, and the presence of a sequence related to ammodytoxin B may explain the reported neurotoxicity following some asp viper envenomations. Thus, this study will help in informing the toxinological basis of clinical signs of envenoming.

Activity of two key toxin groups in Australian elapid venoms show a strong correlation to phylogeny but not to diet

Background

The relative influence of diet and phylogeny on snake venom activity is a poorly understood aspect of snake venom evolution. We measured the activity of two enzyme toxin groups – phospholipase A₂ (PLA₂), and L-amino acid oxidase (LAAO) – in the venom of 39 species of Australian elapids (40% of terrestrial species diversity) and used linear parsimony and BayesTraits to investigate any correlation between enzyme activity and phylogeny or diet.

Results

PLA₂ activity ranged from 0 to 481 nmol/min/mg of venom, and LAAO activity ranged from 0 to 351 nmol/min/mg. Phylogenetic comparative methods, implemented in BayesTraits showed that enzyme activity was strongly correlated with phylogeny, more so for LAAO activity. For example, LAAO activity was absent in both the Vermicella and Pseudonaja/Oxyuranus clade, supporting previously proposed relationships among these disparate taxa. There was no association between broad dietary categories and either enzyme activity. There was strong evidence for faster initial rates of change over evolutionary time for LAAO (delta parameter mean 0.2), but no such pattern in PLA₂ (delta parameter mean 0.64). There were some exceptions to the phylogenetic patterns of enzyme activity: different PLA₂ activity in the ecologically similar sister-species Denisonia devisi and D. maculata; large inter-specific differences in PLA₂ activity in Hoplocephalus and Austrelaps.

Conclusions

We have shown that phylogeny is a stronger influence on venom enzyme activity than diet for two of the four major enzyme families present in snake venoms. PLA₂ and LAAO activities had contrasting evolutionary dynamics with the higher delta value for PLA₂ Some species/individuals lacked activity in one protein family suggesting that the loss of single protein family may not incur a significant fitness cost.

The Integrity of α-β-α Sandwich Conformation Is Essential for a Novel Adjuvant TFPR1 to Maintain Its Adjuvanticity

TFPR1 is a novel peptide vaccine adjuvant we recently discovered. To define the structural basis and optimize its application as an adjuvant, we designed three different truncated fragments that have removed dominant B epitopes on TFPR1, and evaluated their capacity to activate bone marrow-derived dendritic cells and their adjuvanticity. Results demonstrated that the integrity of an α-β-α sandwich conformation is essential for TFPR1 to maintain its immunologic activity and adjuvanticity. We obtained a functional truncated fragment TFPR-ta ranging from 40-168 aa of triflin that has similar adjuvanticity as TFPR1 but with 2-log fold lower immunogenicity. These results demonstrated a novel approach to evaluate and improve the activity of protein-based vaccine adjuvant.

Beyond the 'big four': Venom profiling of the medically important yet neglected Indian snakes reveals disturbing antivenom deficiencies

Background

Snakebite in India causes the highest annual rates of death (46,000) and disability (140,000) than any other country. Antivenom is the mainstay treatment of snakebite, whose manufacturing protocols, in essence, have remained unchanged for over a century. In India, a polyvalent antivenom is produced for the treatment of envenomations from the so called ‘big four’ snakes: the spectacled cobra (Naja naja), common krait (Bungarus caeruleus), Russell’s viper (Daboia russelii), and saw-scaled viper (Echis carinatus). In addition to the ‘big four’, India is abode to many other species of venomous snakes that have the potential to inflict severe clinical or, even, lethal envenomations in their human bite victims. Unfortunately, specific antivenoms are not produced against these species and, instead, the ‘big four’ antivenom is routinely used for the treatment.

Methods

We characterized the venom compositions, biochemical and pharmacological activities and toxicity profiles (mouse model) of the major neglected yet medically important Indian snakes (E. c. sochureki, B. sindanus, B. fasciatus, and two populations of N. kaouthia) and their closest ‘big four’ congeners. By performing WHO recommended in vitro and in vivo preclinical assays, we evaluated the efficiencies of the commercially marketed Indian antivenoms in recognizing venoms and neutralizing envenomations by these neglected species.

Findings

As a consequence of dissimilar ecologies and diet, the medically important snakes investigated exhibited dramatic inter- and intraspecific differences in their venom profiles. Currently marketed antivenoms were found to exhibit poor dose efficacy and venom recognition potential against the ‘neglected many’. Premium Serums antivenom failed to neutralise bites from many of the neglected species and one of the ‘big four’ snakes (North Indian population of B. caeruleus).

Conclusions

This study unravels disturbing deficiencies in dose efficacy and neutralisation capabilities of the currently marketed Indian antivenoms, and emphasises the pressing need to develop region-specific snakebite therapy for the ‘neglected many’.

Snakebite is a ‘neglected tropical disease’ that majorly affects the rural populations in developing countries. India bears the brunt of snakebites with over 46,000 deaths and 140,000 disabilities, annually. A significant number of these bites are attributed to the widely distributed ‘big four’ snakes, namely spectacled cobra (Naja naja), common krait (Bungarus caeruleus), Russell’s viper (Daboia russelii), and saw-scaled viper (Echis carinatus). The commercial antivenoms marketed in India are only manufactured against these four species, while neglecting many other medically relevant snakes with restricted geographic distribution. Snakebite pathology is dependent on the venom composition of the population/species, which can, in turn, vary intra- and inter-specifically. Though this variation severely limits the cross-population/species antivenom efficacy, envenomations by the neglected snakes in India are treated with the ‘big four’ antivenom. Therefore, to unravel the underlying venom variability, we investigated venom proteomic, biochemical/pharmacological and toxicity profiles of the major neglected Indian snakes and their ‘big four’ relatives. To assess the effectiveness of the ‘big four’ antivenom in treating bites from these neglected snakes, we performed preclinical experiments, which revealed alarming inadequacies of the commercial antivenoms. Our findings accentuate the compelling necessity for the innovation of highly efficacious next-generation snakebite therapy in India.

TFPR1 acts as an immune regulator and an efficient adjuvant for proteins and peptides by activating immune cells, primarily through TLR2

Triflin, a non-toxic protein found in the venom of the Habu snake, belongs to the CRISP (cysteine-rich secretory protein) family, which comprises two domains: a C-terminal cysteine-rich domain (CRD) and an N-terminal pathogenesis-related-1 (PR-1) domain. The function of the highly structurally conserved PR-1 domain is unknown. Here, we successfully expressed the PR-1 domain of triflin (hereafter called TFPR1) in E. coli. Animal experiments showed that TFPR1 augmented Th1-biased antibody- and cell-mediated immune responses in mice immunized with two protein antigens (OVA and HBsAg) or a peptide antigen (HIV-1 pep). A flow cytometry-based binding assay and in vitro stimulation with TFPR1 showed that it triggered Th1-biased proinflammatory and immunoregulatory cytokine secretion primarily by binding to B cells and macrophages within the mouse splenocyte population. Quantitative RT-PCR, antibody blocking assays using a specific anti-mTLR2 antibody, and stimulatory experiments in vitro using splenocytes from TLR2-KO mice demonstrated that TFPR1 activated murine immune cells, primarily by stimulating toll-like receptor 2 (TLR2). These results suggest that TFPR1 acts as a novel immune modulator and potent adjuvant primarily by activating TLR2. Thus, the PR-1-based core domain might play a role in immune regulation.

Mapping the chemical and sequence space of the ShKT superfamily

The ShKT superfamily is widely distributed throughout nature and encompasses a wide range of documented functions and processes, from modulation of potassium channels to involvement in morphogenesis pathways. Cysteine-rich secretory proteins (CRISPs) contain a cysteine-rich domain (CRD) at the C-terminus that is similar in structure to the ShK fold. Despite the structural similarity of the CRD and ShK-like domains, we know little of the sequence-function relationships in these families. Here, for the first time, we examine the evolution of the biophysical properties of sequences within the ShKT superfamily in relation to function, with a focus on the ShK-like superfamily. ShKT data were sourced from published sequences in the protein family database, in addition to new ShK-like sequences from the Australian speckled anemone (Oulactis sp.). Our analysis clearly delineates the ShK-like family from the CRDs of CRISP proteins. The four CRISP subclusters separate out into the main phyla of Mammalia, Insecta and Reptilia. The ShK-like family is in turn composed of seven subclusters, the largest of which contains members from across the eukaryotes, with a continuum of intermediate properties. Smaller sub-clusters contain specialised members such as nematode ShK-like sequences. Several of these ShKT sub-clusters contain no functionally characterised sequences. This chemical space analysis should be useful as a guide to select sequences for functional studies and to gain insight into the evolution of these highly divergent sequences with an ancient conserved fold.

Epididymal cysteine-rich secretory proteins are required for epididymal sperm maturation and optimal sperm function

STUDY QUESTION

What is the role of epididymal cysteine-rich secretory proteins (CRISPs) in male fertility?

SUMMARY ANSWER

While epididymal CRISPs are not absolutely required for male fertility, they are required for optimal sperm function.

WHAT IS KNOWN ALREADY

CRISPs are members of the CRISP, Antigen 5 and Pathogenesis related protein 1 (CAP) superfamily and are characterized by the presence of an N-terminal CAP domain and a C-terminal CRISP domain. CRISPs are highly enriched in the male reproductive tract of mammals, including in the epididymis. Within humans there is one epididymal CRISP, CRISP1, whereas in mice there are two, CRISP1 and CRISP4.

STUDY DESIGN, SIZE, DURATION

In order to define the role of CRISPs within the epididymis, Crisp1 and Crisp4 knockout mouse lines were produced then interbred to produce Crisp1 and 4 double knockout (DKO) mice, wherein the expression of all epididymal CRISPs was ablated. Individual and DKO models were then assessed, relative to their own strain-specific wild type littermates for fertility, and sperm output and functional competence at young (10-12 weeks of age) and older ages (22-24 weeks). Crisp1 and 4 DKO and control mice were also compared for their ability to bind to the zona pellucida and achieve fertilization.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Knockout mouse production was achieved using modified embryonic stem cells and standard methods. The knockout of individual genes was confirmed at a mRNA (quantitative PCR) and protein (immunochemistry) level. Fertility was assessed using breeding experiments and a histological assessment of testes and epididymal tissue. Sperm functional competence was assessed using a computer assisted sperm analyser, induction of the acrosome reaction using progesterone followed by staining for acrosome contents, using immunochemical and western blotting to assess the ability of sperm to manifest tyrosine phosphorylation under capacitating conditions and using sperm-zona pellucida binding assays and IVF methods. A minimum of three biological replicates were used per assay and per genotype.

MAIN RESULTS AND THE ROLE OF CHANCE

While epididymal CRISPs are not absolutely required for male fertility, their production results in enhanced sperm function and, depending on context, CRISP1 and CRISP4 act redundantly or autonomously. Specifically, CRISP1 is the most important CRISP in the establishment of normally motile sperm, whereas CRISP4 acts to enhance capacitation-associated tyrosine phosphorylation, and CRISP1 and CRISP4 act together to establish normal acrosome function. Both are required to achieve optimal sperm-egg interaction. The presence of immune infiltrates into the epididymis of older, but not younger, DKO animals also suggests epididymal CRISPs function to produce an immune privileged environment for maturing sperm within the epididymis.

LIMITATIONS REASONS FOR CAUTION

Caution should be displayed in the translation of mouse-derived data into the human wherein the histology of the epididymis is someone what different. The mice used in the study were housed in a specific pathogen-free environment and were thus not exposed to the full range of environmental challenges experienced by wild mice or humans. As such, the role of CRISPs in the maintenance of an immune privileged environment, for example, may be understated.

WIDER IMPLICATIONS OF THE FINDINGS

The combined deletion of Crisp1 and Crisp4 in mice is equivalent to the removal of all CRISP expression in humans. As such, these data suggest that mammalian CRISPs, including that in humans, function to enhance sperm function and thus male fertility. These data also suggest that in the presence of an environmental challenge, CRISPs help to maintain an immune privileged environment and thus, protect against immune-mediated male infertility.

LARGE SCALE DATA

Not applicable.

STUDY FUNDING AND COMPETING INTEREST(S)

This study was funded by the National Health and Medical Research Council, the Victorian Cancer Agency and a scholarship from the Chinese Scholarship Council. The authors have no conflicts of interest to declare.

Proteomic Characterization of Two Medically Important Malaysian Snake Venoms, Calloselasma rhodostoma (Malayan Pit Viper) and Ophiophagus hannah (King Cobra)

Calloselasma rhodostoma (CR) and Ophiophagus hannah (OH) are two medically important snakes found in Malaysia. While some studies have described the biological properties of these venoms, feeding and environmental conditions also influence the concentration and distribution of snake venom toxins, resulting in variations in venom composition. Therefore, a combined proteomic approach using shotgun and gel filtration chromatography, analyzed by tandem mass spectrometry, was used to examine the composition of venoms from these Malaysian snakes. The analysis revealed 114 proteins (15 toxin families) and 176 proteins (20 toxin families) in Malaysian Calloselasma rhodostoma and Ophiophagus hannah species, respectively. Flavin monoamine oxidase, phospholipase A₂, phosphodiesterase, snake venom metalloproteinase, and serine protease toxin families were identified in both venoms. Aminopeptidase, glutaminyl-peptide cyclotransferase along with ankyrin repeats were identified for the first time in CR venom, and insulin, c-type lectins/snaclecs, hepatocyte growth factor, and macrophage colony-stimulating factor together with tumor necrosis factor were identified in OH venom for the first time. Our combined proteomic approach has identified a comprehensive arsenal of toxins in CR and OH venoms. These data may be utilized for improved antivenom production, understanding pathological effects of envenoming, and the discovery of biologically active peptides with medical and/or biotechnological value.

Evolutionary Genetics of Hypoxia and Cold Tolerance in Mammals

Low oxygen and fluctuant ambient temperature pose serious challenges to mammalian survival. Physiological adaptations in mammals to hypoxia and low temperatures have been intensively investigated, yet their underlying molecular mechanisms need further exploration. Independent invasions of high-altitude plateaus, subterranean burrows and marine environments by different mammals provide opportunities to conduct such analyses. Here, we focused on six genes in the hypoxia inducible factor (HIF) pathway and two non-shivering thermogenesis (NST)-related genes [PPAR co-activator 1 (PGC-1) and uncoupling protein 1 (UCP1)] in representative species of pikas and other mammals to understand whether these loci were targeted by natural selection during independent invasions to conditions characterized by hypoxia and temperature fluctuations by high-altitude, subterranean and marine mammals. Our analyses revealed pervasive positive selection signals in the HIF pathway genes of mammals occupying high-altitude, subterranean and aquatic ecosystems; however, the mechanisms underlying their independent adaptations to hypoxic environments varied by taxa, since different genes were positively selected in each taxon and expression levels of individual genes varied among species. Additionally, parallel amino acid substitutions were also detected in hypoxia-tolerant mammals, indicating that convergent evolution may play a role in their independent adaptations to hypoxic environments. However, divergent evolutionary histories of NST-related genes were noted, since significant evidence of positive selection was observed in PGC-1 and UCP1 in high-altitude species and subterranean rodents; however, UCP1 may have already lost its function in diving cetaceans, which may be related to the thick blubber layer of adipose and connective tissue in these mammals.

Snake venom components in medicine: From the symbolic rod of Asclepius to tangible medical research and application

Both mythologically and logically, snakes have always fascinated man. Snakes have attracted both awe and fear not only because of the elegant movement of their limbless bodies, but also because of the potency of their deadly venoms. Practically, in 2017, the world health organization (WHO) listed snake envenomation as a high priority neglected disease, as snakes inflict up to 2.7 million poisonous bites, around 100.000 casualties, and about three times as many invalidities on man. The venoms of poisonous snakes are a cocktail of potent compounds which specifically and avidly target numerous essential molecules with high efficacy. The individual effects of all venom toxins integrate into lethal dysfunctions of almost any organ system. It is this efficacy and specificity of each venom component, which after analysis of its structure and activity may serve as a potential lead structure for chemical imitation. Such toxin mimetics may help in influencing a specific body function pharmaceutically for the sake of man's health. In this review article, we will give some examples of snake venom components which have spurred the development of novel pharmaceutical compounds. Moreover, we will provide examples where such snake toxin-derived mimetics are in clinical use, trials, or consideration for further pharmaceutical exploitation, especially in the fields of hemostasis, thrombosis, coagulation, and metastasis. Thus, it becomes clear why a snake captured its symbolic place at the Asclepius rod with good reason still nowadays.

Venomics: A Mini-Review

Venomics is the integration of proteomic, genomic and transcriptomic approaches to study venoms. Advances in these approaches have enabled increasingly more comprehensive analyses of venoms to be carried out, overcoming to some extent the limitations imposed by the complexity of the venoms and the small quantities that are often available. Advances in bioinformatics and high-throughput functional assay screening approaches have also had a significant impact on venomics. A combination of all these techniques is critical for enhancing our knowledge on the complexity of venoms and their potential therapeutic and agricultural applications. Here we highlight recent advances in these fields and their impact on venom analyses.

Evolutionary changes of Hox genes and relevant regulatory factors provide novel insights into mammalian morphological modifications

The diversity of body plans of mammals accelerates the innovation of lifestyles and the extensive adaptation to different habitats, including terrestrial, aerial and aquatic habitats. However, the genetic basis of those phenotypic modifications, which have occurred during mammalian evolution, remains poorly explored. In the present study, we synthetically surveyed the evolutionary pattern of Hox clusters that played a powerful role in the morphogenesis along the head–tail axis of animal embryos and the main regulatory factors (Mll, Bmi1 and E2f6) that control the expression of Hox genes. A deflected density of repetitive elements and lineage‐specific radical mutations of Mll have been determined in marine mammals with morphological changes, suggesting that evolutionary changes may alter Hox gene expression in these lineages, leading to the morphological modification of these lineages. Although no positive selection was detected at certain ancestor nodes of lineages, the increased ω values of Hox genes implied the relaxation of functional constraints of these genes during the mammalian evolutionary process. More importantly, 49 positively‐selected sites were identified in mammalian lineages with phenotypic modifications, indicating adaptive evolution acting on Hox genes and regulatory factors. In addition, 3 parallel amino acid substitutions in some Hox genes were examined in marine mammals, which might be responsible for their streamlined body.

What killed Karl Patterson Schmidt? Combined venom gland transcriptomic, venomic and antivenomic analysis of the South African green tree snake (the boomslang), Dispholidus typus

Background

Non-front-fanged colubroid snakes comprise about two-thirds of extant ophidian species. The medical significance of the majority of these snakes is unknown, but at least five species have caused life-threatening or fatal human envenomings. However, the venoms of only a small number of species have been explored.

Methods

A combined venomic and venom gland transcriptomic approach was employed to characterise of venom of Dispholidus typus (boomslang), the snake that caused the tragic death of Professor Karl Patterson Schmidt. The ability of CroFab™ antivenom to immunocapture boomslang venom proteins was investigated using antivenomics.

Results

Transcriptomic-assisted proteomic analysis identified venom proteins belonging to seven protein families: three-finger toxin (3FTx); phospholipase A₂ (PLA₂); cysteine-rich secretory proteins (CRISP); snake venom (SV) serine proteinase (SP); C-type lectin-like (CTL); SV metalloproteinases (SVMPs); and disintegrin-like/cysteine-rich (DC) proteolytic fragments. CroFab™ antivenom efficiently immunodepleted some boomslang SVMPs.

Conclusions

The present work is the first to address the overall proteomic profile of D. typus venom. This study allowed us to correlate the toxin composition with the toxic activities of the venom. The antivenomic analysis suggested that the antivenom available at the time of the unfortunate accident could have exhibited at least some immunoreactivity against the boomslang SVMPs responsible for the disseminated intravascular coagulation syndrome that caused K.P. Schmidt's fatal outcome.

General significance

This study may stimulate further research on other non-front-fanged colubroid snake venoms capable of causing life-threatening envenomings to humans, which in turn should contribute to prevent fatal human accidents, such as that unfortunately suffered by K.P. Schmidt.

•

The venom proteome of Dispholidus typus (boomslang) is reported.

•

Transcriptomic-assisted proteomic analysis identified venom proteins belonging to seven protein families.

•

Boomslang venom proteome is dominated (75%) by snake venom PIII-metalloproteinases (PIII-SVMPs).

•

CroFab™ antivenom efficiently immunodepleted some boomslang PIII-SVMPs.

Understanding Biological Roles of Venoms Among the Caenophidia: The Importance of Rear-Fanged Snakes

Snake venoms represent an adaptive trophic response to the challenges confronting a limbless predator for overcoming combative prey, and this chemical means of subduing prey shows several dominant phenotypes. Many front-fanged snakes, particularly vipers, feed on various vertebrate and invertebrate prey species, and some of their venom components (e.g., metalloproteinases, cobratoxin) appear to have been selected for "broad-brush" incapacitation of different prey taxa. Using proteomic and genomic techniques, the compositional diversity of front-fanged snakes is becoming well characterized; however, this is not the case for most rear-fanged colubroid snakes. Because these species consume a high diversity of prey, and because venoms are primarily a trophic adaptation, important clues for understanding specific selective pressures favoring venom component composition will be found among rear-fanged snake venoms. Rear-fanged snakes typically (but not always) produce venoms with lower complexity than front-fanged snakes, and there are even fewer dominant (and, arguably, biologically most relevant) venom protein families. We have demonstrated taxon-specific toxic effects, where lizards and birds show high susceptibility while mammals are largely unaffected, for both Old World and New World rear-fanged snakes, strongly indicating a causal link between toxin evolution and prey preference. New data are presented on myotoxin a, showing that the extremely rapid paralysis induced by this rattlesnake toxin is specific for rodents, and that myotoxin a is ineffectual against lizards. Relatively few rear-fanged snake venoms have been characterized, and basic natural history data are largely lacking, but directed sampling of specialized species indicates that novel compounds are likely among these specialists, particularly among those species feeding on invertebrate prey such as scorpions and centipedes. Because many of the more than 2200 species of colubroid snakes are rear-fanged, and many possess a Duvernoy's venom gland, understanding the nature of their venoms is foundational to understanding venom evolution in advanced snakes.

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.

Canopy Venom: Proteomic Comparison among New World Arboreal Pit-Viper Venoms

Central and South American pitvipers, belonging to the genera Bothrops and Bothriechis, have independently evolved arboreal tendencies. Little is known regarding the composition and activity of their venoms. In order to close this knowledge gap, venom proteomics and toxin activity of species of Bothriechis, and Bothrops (including Bothriopsis) were investigated through established analytical methods. A combination of proteomics and bioactivity techniques was used to demonstrate a similar diversification of venom composition between large and small species within Bothriechis and Bothriopsis. Increasing our understanding of the evolution of complex venom cocktails may facilitate future biodiscoveries.

Venomics of Vipera berus berus to explain differences in pathology elicited by Vipera ammodytes ammodytes envenomation: Therapeutic implications

Vipera berus berus (Vbb) is the most widely distributed and Vipera ammodytes ammodytes (Vaa) the most venomous viper in Europe. In particular areas of the Old continent their toxic bites constitute a considerable public health problem. To make the current envenomation therapy more effective we have analysed the proteome of Vbb venom and compared it with that of Vaa. We found the proteome of Vbb to be much less complex and to contain smaller levels of particularly snaclecs and sPLA2s. Snaclecs are probably responsible for thrombocytopenia. The neurotoxic sPLA2s, ammodytoxins, are responsible for the most specific feature of the Vaa venom poisoning - induction of signs of neurotoxicity in patients. These molecules were not found in Vbb venom. Both venoms induce haemorrhage and coagulopathy in man. As Vaa and Vbb venoms possess homologous P-III snake venom metalloproteinases, the main haemorrhagic factors, the severity of the haemorrhage is dictated by concentration and specific activity of these molecules. The much greater anticoagulant effect of Vaa venom than that of Vbb venom lies in its higher extrinsic pathway coagulation factor-proteolysing activity and content of ammodytoxins which block the prothrombinase complex formation.

BIOLOGICAL SIGNIFICANCE

Envenomations by venomous snakes constitute a considerable public health problem worldwide, and also in Europe. In the submitted work we analysed the venom proteome of Vipera berus berus (Vbb), the most widely distributed venomous snake in Europe and compared it with the venom proteome of the most venomous viper in Europe, Vipera ammodytes ammodytes (Vaa). We have offered a possible explanation, at the molecular level, for the differences in clinical pictures inflicted by the Vbb and Vaa venoms. We have provided an explanation for the effectiveness of treatment of Vbb envenomation by Vaa antiserum and explained why full protection of Vaa venom poisoning by Vbb antiserum should not be always expected, especially not in cases of severe poisoning. The latter makes a strong case for Vaa antiserum production as we are faced with its shortage due to ceasing of production of two most frequently used products.

Identification and characterization of toxins in the venom gland of the Chinese bird spider, Haplopelma hainanum, by transcriptomic analysis

Tarantula venoms provide a model system for studying toxin selectivity, structure-activity relationships and molecular evolution of peptide toxins. Previous studies have identified a large number of peptide toxins in the venom of the Chinese bird spider Haplopelma hainanum, generally regarded as a highly venomous spider. However, the lack of available RNA-seq transcriptomic and genomic data is an obstacle to understanding its venom at the molecular level. In this study, we investigated the venom gland transcriptome of H. hainanum by RNA-seq, in the absence of an available genomic sequence. We identified 201 potential toxins among 57 181 de novo assembled transcripts, including knottins, Kunitz-type toxins, enzymes and other proteins. We systematically identified most of the knottins and Kunitz-type toxins, some of which showed strongly biased expression in the venom gland, including members of the huwentoxin-1, huwentoxin-2 and magi-1 families. We also discovered several novel potential toxins. These data demonstrate the high molecular and structural diversity in the venom toxins of H. hainanum. This study offers a useful strategy for exploring the complex components of spider venoms.

Venom gland transcriptome analyses of two freshwater stingrays (Myliobatiformes: Potamotrygonidae) from Brazil

Stingrays commonly cause human envenoming related accidents in populations of the sea, near rivers and lakes. Transcriptomic profiles have been used to elucidate components of animal venom, since they are capable of providing molecular information on the biology of the animal and could have biomedical applications. In this study, we elucidated the transcriptomic profile of the venom glands from two different freshwater stingray species that are endemic to the Paraná-Paraguay basin in Brazil, Potamotrygon amandae and Potamotrygon falkneri. Using RNA-Seq, we identified species-specific transcripts and overlapping proteins in the venom gland of both species. Among the transcripts related with envenoming, high abundance of hyaluronidases was observed in both species. In addition, we built three-dimensional homology models based on several venom transcripts identified. Our study represents a significant improvement in the information about the venoms employed by these two species and their molecular characteristics. Moreover, the information generated by our group helps in a better understanding of the biology of freshwater cartilaginous fishes and offers clues for the development of clinical treatments for stingray envenoming in Brazil and around the world. Finally, our results might have biomedical implications in developing treatments for complex diseases.