Molecular evolution of vertebrate neurotrophins: co-option of the highly conserved nerve growth factor gene into the advanced snake venom arsenalf

Assessing Target Specificity of the Small Molecule Inhibitor MARIMASTAT to Snake Venom Toxins: A Novel Application of Thermal Proteome Profiling

New treatments that circumvent the pitfalls of traditional antivenom therapies are critical to address the problem of snakebite globally. Numerous snake venom toxin inhibitors have shown promising cross-species neutralization of medically significant venom toxins in vivo and in vitro. The development of high-throughput approaches for the screening of such inhibitors could accelerate their identification, testing, and implementation and thus holds exciting potential for improving the treatments and outcomes of snakebite envenomation worldwide. Energetics-based proteomic approaches, including thermal proteome profiling and proteome integral solubility alteration (PISA) assays, represent "deep proteomics" methods for high throughput, proteome-wide identification of drug targets and ligands. In the following study, we apply thermal proteome profiling and PISA methods to characterize the interactions between venom toxin proteoforms in Crotalus atrox (Western Diamondback Rattlesnake) and the snake venom metalloprotease (SVMP) inhibitor marimastat. We investigate its venom proteome-wide effects and characterize its interactions with specific SVMP proteoforms, as well as its potential targeting of non-SVMP venom toxin families. We also compare the performance of PISA thermal window and soluble supernatant with insoluble precipitate using two inhibitor concentrations, providing the first demonstration of the utility of a sensitive high-throughput PISA-based approach to assess the direct targets of small molecule inhibitors for snake venom.

Investigating In Situ Expression of Neurotrophic Factors and Partner Proteins in Irreversible Pulpitis

INTRODUCTION

In situ assessments of neurotrophic factors and their associated molecular partners have not been explored to date, particularly in humans. The present investigation aimed to explore the expressional dysregulation of neurotrophic factors (nerve growth factor [NGF], brain derived neurotrophic factor [BDNF], and NT4/5), their receptors (TrkA and TrkB), and their modulators (USP36 and Nedd4-2) directly in irreversibly inflamed human pulp tissues.

METHODS

Forty samples each of healthy and irreversibly inflamed pulp were extirpated for the study. Immunohistochemical examinations were carried out for the anatomic changes and expression of neurotrophic factors and partner proteins. Expression was digitally quantified using the IHC profiler module of ImageJ and deduced as optical density. Statistical analyses were carried out by GraphPad Prism.

RESULTS

Decrease in nuclear and vessel diameters was observed in irreversibly inflamed pulp tissues. NGF and BDNF were found to be significantly upregulated in symptomatic irreversible pulpitis (SIP), whereas no significant difference was observed in the expression of TrkA and TrkB. Expression of Nedd4-2, USP36, and TrkA was found positively correlated with the NGF in healthy pulp tissues. However, in SIP, positive correlation was only observed between the expression of USP36 and NGF. Among the ligands, BDNF expression was found positively correlated with NGF in healthy pulp but not with NT4/5. In the case of SIP, no correlation was observed between any neurotrophic factors.

CONCLUSIONS

Upregulation of NGF, BDNF, USP36 and Nedd4-2 in SIP indicates dysregulation in the molecular events underlying the disease biology and could be exploited as potential markers for the disease diagnosis.

Evaluation of Neutralization Potential of Naja naja and Daboiarusselii Snake Venom by Root Extract of Cyanthillium cinereum

Aim

One of the main reasons for the death due to snake bites is the non-availability of antivenoms in the regions where they are needed. The use of medicinal plants and plant-based natural products as an alternative to antivenom will become a milestone in snake bite envenomation. The present study investigates the in vitro antivenom properties of Cyanthillium cinereum root extracts.

Materials and methods

The C. cinereum root's aqueous extract was prepared by the Soxhlet extraction method, and phytochemical screening was performed to detect the presence of various bioactive compounds. Thin-layer chromatography (TLC) and gas chromatography-mass spectrometry (GC-MS) analysis were performed for the detection and identification of phytochemical constituents. In this study, an in vitro model is used to assess the antivenom capability of aqueous extract. Venom toxicity and neutralization assays were as follows: An in vitro pharmacological evaluation was performed by direct hemolysis assay, indirect hemolytic assay, proteolytic activity, neutralization of procoagulant activity, and gelatin liquefaction method.

Results

Qualitative analysis of phytochemicals by the standard method showed the presence of various phytochemical constituents. Also, GC-MS analysis showed the presence of three major compounds that possess antivenom activity from the obtained 60 bioactive compounds, and their chemical structures were also determined. Venom protein profiling was performed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) analysis. The plant extract was able to neutralize the Naja naja (N. naja) and Daboia russelii (D. russelii) venom induced hemolysis and it was reduced below 50 and 40%, respectively and the extract was also able to reduce the hemolytic halo produced by venoms. Procoagulant activity and gelatin liquefaction assay showed that venom-induced clotting was neutralized by increasing the root extract concentration sufficiently.

Conclusion

The aqueous extract of the root of C. cinereum showed potent in vitro venom-neutralizing activity, and it can be used as a formidable therapeutic agent against N. naja and D. russelii envenomation.

How to cite this article

Suji S, Dinesh MD, Keerthi KU, Anagha KP, Arya J, Anju KV. Evaluation of Neutralization Potential of Naja naja and Daboia russelii Snake Venom by Root Extract of Cyanthillium cinereum. Indian J Crit Care Med 2023;27(11):821-829.

A current perspective on snake venom composition and constituent protein families

Snake venoms are heterogeneous mixtures of proteins and peptides used for prey subjugation. With modern proteomics there has been a rapid expansion in our knowledge of snake venom composition, resulting in the venom proteomes of 30% of vipers and 17% of elapids being characterised. From the reasonably complete proteomic coverage of front-fanged snake venom composition (179 species-68 species of elapids and 111 species of vipers), the venoms of vipers and elapids contained 42 different protein families, although 18 were only reported in < 5% of snake species. Based on the mean abundance and occurrence of the 42 protein families, they can be classified into 4 dominant, 6 secondary, 14 minor, and 18 rare protein families. The dominant, secondary and minor categories account for 96% on average of a snake's venom composition. The four dominant protein families are: phospholipase A₂ (PLA₂), snake venom metalloprotease (SVMP), three-finger toxins (3FTx), and snake venom serine protease (SVSP). The six secondary protein families are: L-amino acid oxidase (LAAO), cysteine-rich secretory protein (CRiSP), C-type lectins (CTL), disintegrins (DIS), kunitz peptides (KUN), and natriuretic peptides (NP). Venom variation occurs at all taxonomic levels, including within populations. The reasons for venom variation are complex, as variation is not always associated with geographical variation in diet. The four dominant protein families appear to be the most important toxin families in human envenomation, being responsible for coagulopathy, neurotoxicity, myotoxicity and cytotoxicity. Proteomic techniques can be used to investigate the toxicological profile of a snake venom and hence identify key protein families for antivenom immunorecognition.

The genes from the pseudoautosomal region 1 (PAR1) of the mammalian sex chromosomes: Synteny, phylogeny and selection

Sex chromosomes recombine restrictly in their homologous area, the pseudoautosomal region (PAR), represented by PAR1 and PAR2, which behave like an autosome in both pairing and recombination. The PAR1, common to most of the eutherian mammals, is located at the terminus of the sex chromosomes short arm and exhibit recombination rates ~20 times higher than the autosomes. Here, we assessed the interspecific evolutionary genomic dynamics of 15 genes of the PAR1 across 41 mammalian genera (representing six orders). The strong negative selection detected in most of the assessed groups reinforces the presence of evolutionary constraints, imposed by the important function of the PAR1 genes. Indeed, mutations in these genes are associated with various diseases in humans, including stature problems (Klinefelter Syndrome), leukemia and mental diseases. Yet, a few genes exhibiting positive selection (ω-value >1) were depicted in Rodentia (ASMT and ZBED1) and Primates (CRLF2 and CSF2RA). Rodents have the smallest described PAR1, while that of simian primates/humans underwent a 3 to 5 fold size reduction. The assessment of the PAR1 genes synteny revealed differences among the mammalian species, especially in the Rodentia order where chromosomic translocations from the sex chromosomes to the autosomes were observed. Such syntenic changes may be an evidence of the rapid evolution in rodents, as previous referred in other papers, also depicted by their increased branch lengths in the phylogenetic analyses. Concluding, we suggest that genome migration is an important factor influencing the evolution of mammals and may result in changes of the selective pressures operating on the genome.

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.

A Meta-Analysis of the Protein Components in Rattlesnake Venom

The specificity and potency of venom components give them a unique advantage in developing various pharmaceutical drugs. Though venom is a cocktail of proteins, rarely are the synergy and association between various venom components studied. Understanding the relationship between various components of venom is critical in medical research. Using meta-analysis, we observed underlying patterns and associations in the appearance of the toxin families. For Crotalus, Dis has the most associations with the following toxins: PDE; BPP; CRL; CRiSP; LAAO; SVMP P-I and LAAO; SVMP P-III and LAAO. In Sistrurus venom, CTL and NGF have the most associations. These associations can predict the presence of proteins in novel venom and understand synergies between venom components for enhanced bioactivity. Using this approach, the need to revisit the classification of proteins as major components or minor components is highlighted. The revised classification of venom components is based on ubiquity, bioactivity, the number of associations, and synergies. The revised classification can be expected to trigger increased research on venom components, such as NGF, which have high biomedical significance. Using hierarchical clustering, we observed that the genera's venom compositions were similar, based on functional characteristics rather than phylogenetic relationships.

Putative Antimicrobial Peptides of the Posterior Salivary Glands from the Cephalopod Octopus vulgaris Revealed by Exploring a Composite Protein Database

Cephalopods, successful predators, can use a mixture of substances to subdue their prey, becoming interesting sources of bioactive compounds. In addition to neurotoxins and enzymes, the presence of antimicrobial compounds has been reported. Recently, the transcriptome and the whole proteome of the Octopus vulgaris salivary apparatus were released, but the role of some compounds—e.g., histones, antimicrobial peptides (AMPs), and toxins—remains unclear. Herein, we profiled the proteome of the posterior salivary glands (PSGs) of O. vulgaris using two sample preparation protocols combined with a shotgun-proteomics approach. Protein identification was performed against a composite database comprising data from the UniProtKB, all transcriptomes available from the cephalopods’ PSGs, and a comprehensive non-redundant AMPs database. Out of the 10,075 proteins clustered in 1868 protein groups, 90 clusters corresponded to venom protein toxin families. Additionally, we detected putative AMPs clustered with histones previously found as abundant proteins in the saliva of O. vulgaris. Some of these histones, such as H2A and H2B, are involved in systemic inflammatory responses and their antimicrobial effects have been demonstrated. These results not only confirm the production of enzymes and toxins by the O. vulgaris PSGs but also suggest their involvement in the first line of defense against microbes.

The functional equivalence of two variants of the suboptimal choice task: choice proportion and response latency as measures of value

In the suboptimal-choice task, birds systematically choose the leaner but informative option (suboptimal) over the richer but non-informative option (optimal). The task has two variations. In the standard task, the optimal option includes two terminal link stimuli. In the original task, it includes a single terminal link stimulus. Two models, the temporal information account (Cunningham and Shahan, J Exp Psychol Anim Learn Cogn 44:1-22, 2018) and the ∆-∑ hypothesis (González et al., J Exp Anal Behav 113:591-608, 2020), presuppose that these procedures are equivalent, but no formal comparison is available. Here we test whether or not these procedures are functionally equivalent. One group of pigeons was trained with the standard procedure, another group with the original procedure, and a third group was trained with a hybrid of the other two (i.e., the two options were the optimal links of the standard and original procedures). Our findings indicate that the number of terminal link stimuli in the optimal option is inconsequential vis-à-vis choice. Moreover, our findings also indicate that latencies to respond are a sensitive metric of value and choice. As predicted by the Sequential Choice Model, we were able to predict simultaneous choices from the latencies of sequential choices and observed a substantial shortening of latencies during simultaneous choices.

Comparative venom proteomics of banded krait (Bungarus fasciatus) from five geographical locales: Correlation of venom lethality, immunoreactivity and antivenom neutralization

The banded krait, Bungarus fasciatus is a medically important venomous snake in Asia. The wide distribution of this species in Southeast Asia and southern China indicates potential geographical variation of the venom which may impact the clinical management of snakebite envenomation. This study investigated the intraspecific venom variation of B. fasciatus from five geographical locales through a venom decomplexing proteomic approach, followed by toxinological and immunological studies. The venom proteomes composed of a total of 9 toxin families, comprising 22 to 31 proteoforms at varying abundances. The predominant proteins were phospholipase A₂ (including beta-bungarotoxin), Kunitz-type serine protease inhibitor (KSPI) and three-finger toxins (3FTx), which are toxins that cause neurotoxicity and lethality. The venom lethality varied with geographical origins of the snake, with intravenous median lethal doses (LD₅₀) ranging from 0.45-2.55 µg/g in mice. The Thai Bungarus fasciatus monovalent antivenom (BFMAV) demonstrated a dose-dependent increasing immunological binding activity toward all venoms; however, its in vivo neutralization efficacy varied vastly with normalized potency values ranging from 3 to 28 mg/g, presumably due to the compositional differences of dominant proteins in the different venoms. The findings support that antivenom use should be optimized in different geographical areas. The development of a pan-regional antivenom may be a more sustainable solution for the treatment of snakebite envenomation.

Avian Binocularity and Adaptation to Nocturnal Environments: Genomic Insights from a Highly Derived Visual Phenotype

Typical avian eyes are phenotypically engineered for photopic vision (daylight). In contrast, the highly derived eyes of the barn owl (Tyto alba) are adapted for scotopic vision (dim light). The dramatic modifications distinguishing barn owl eyes from other birds include: 1) shifts in frontal orientation to improve binocularity, 2) rod-dominated retina, and 3) enlarged corneas and lenses. Some of these features parallel mammalian eye patterns, which are hypothesized to have initially evolved in nocturnal environments. Here, we used an integrative approach combining phylogenomics and functional phenotypes of 211 eye-development genes across 48 avian genomes representing most avian orders, including the stem lineage of the scotopic-adapted barn owl. Overall, we identified 25 eye-development genes that coevolved under intensified or relaxed selection in the retina, lens, cornea, and optic nerves of the barn owl. The agtpbp1 gene, which is associated with the survival of photoreceptor populations, was pseudogenized in the barn owl genome. Our results further revealed that barn owl retinal genes responsible for the maintenance, proliferation, and differentiation of photoreceptors experienced an evolutionary relaxation. Signatures of relaxed selection were also observed in the lens and cornea morphology-associated genes, suggesting that adaptive evolution in these structures was essentially structural. Four eye-development genes (ephb1, phactr4, prph2, and rs1) evolved in positive association with the orbit convergence in birds and under relaxed selection in the barn owl lineage, likely contributing to an increased reliance on binocular vision in the barn owl. Moreover, we found evidence of coevolutionary interactions among genes that are expressed in the retina, lens, and optic nerve, suggesting synergetic adaptive events. Our study disentangles the genomic changes governing the binocularity and low-light perception adaptations of barn owls to nocturnal environments while revealing the molecular mechanisms contributing to the shift from the typical avian photopic vision to the more-novel scotopic-adapted eye.

Neofunctionalization of the UCP1 mediated the non-shivering thermogenesis in the evolution of small-sized placental mammals

The acquisition of UCP1-mediated non-shivering thermogenesis (NST) was an important event during the evolution of mammals. Here, we assessed the thermogenic neofunctionalization that occurred in the mammalian UCP1, by performing detailed comparative evolutionary genomics analyses (including phylogenetic and selection analyses) of the UCP family members across all major vertebrate classes. Heterogeneously distributed positive selection signatures were found in several UCPs, being preferably located in the mitochondrial matrix domains. Additionally, comparisons with non-mammalian orthologs showed increased evolutionary rates of the mammalian UCP1, not observable in the phylogenetically related UCP2 and UCP3 paralogs. Also, parallel signatures of episodic positive selection (ω > 1) were found in the ancestral branches of both Glires (rodents and lagomorphs) and Afroinsectivores (afrosoricids and macroscelids), underlining the importance of the UCP1 thermogenic activity in these mammalian groups. Finally, we hypothesize that the independent positive selection events that occurred in these two lineages resulted in two UCP1-mediated NST approaches, namely the cold acute response in the Glires and the reproduction success enhancement in the Afroinsectivores.

The Urgent Need to Develop Novel Strategies for the Diagnosis and Treatment of Snakebites

Snakebite envenoming (SBE) is a priority neglected tropical disease, which kills in excess of 100,000 people per year. Additionally, many millions of survivors also suffer through disabilities and long-term health consequences. The only treatment for SBE, antivenom, has a number of major associated problems, not least, adverse reactions and limited availability. This emphasises the necessity for urgent improvements to the management of this disease. Administration of antivenom is too frequently based on symptomatology, which results in wasting crucial time. The majority of SBE-affected regions rely on broad-spectrum polyvalent antivenoms that have a low content of case-specific efficacious immunoglobulins. Research into small molecular therapeutics such as varespladib/methyl-varespladib (PLA2 inhibitors) and batimastat/marimastat (metalloprotease inhibitors) suggest that such adjunctive treatments could be hugely beneficial to victims. Progress into toxin-specific monoclonal antibodies as well as alternative binding scaffolds such as aptamers hold much promise for future treatment strategies. SBE is not implicit during snakebite, due to venom metering. Thus, the delay between bite and symptom presentation is critical and when symptoms appear it may often already be too late to effectively treat SBE. The development of reliable diagnostical tools could therefore initiate a paradigm shift in the treatment of SBE. While the complete eradication of SBE is an impossibility, mitigation is in the pipeline, with new treatments and diagnostics rapidly emerging. Here we critically review the urgent necessity for the development of diagnostic tools and improved therapeutics to mitigate the deaths and disabilities caused by SBE.

Comprehensive Study of the Proteome and Transcriptome of the Venom of the Most Venomous European Viper: Discovery of a New Subclass of Ancestral Snake Venom Metalloproteinase Precursor-Derived Proteins

The nose-horned viper, its nominotypical subspecies Vipera ammodytes ammodytes (Vaa), in particular, is, medically, one of the most relevant snakes in Europe. The local and systemic clinical manifestations of poisoning by the venom of this snake are the result of the pathophysiological effects inflicted by enzymatic and nonenzymatic venom components acting, most prominently, on the blood, cardiovascular, and nerve systems. This venom is a very complex mixture of pharmacologically active proteins and peptides. To help improve the current antivenom therapy toward higher specificity and efficiency and to assist drug discovery, we have constructed, by combining transcriptomic and proteomic analyses, the most comprehensive library yet of the Vaa venom proteins and peptides. Sequence analysis of the venom gland cDNA library has revealed the presence of messages encoding 12 types of polypeptide precursors. The most abundant are those for metalloproteinase inhibitors (MPis), bradykinin-potentiating peptides (BPPs), and natriuretic peptides (NPs) (all three on a single precursor), snake C-type lectin-like proteins (snaclecs), serine proteases (SVSPs), P-II and P-III metalloproteinases (SVMPs), secreted phospholipases A₂ (sPLA₂s), and disintegrins (Dis). These constitute >88% of the venom transcriptome. At the protein level, 57 venom proteins belonging to 16 different protein families have been identified and, with SVSPs, sPLA₂s, snaclecs, and SVMPs, comprise ∼80% of all venom proteins. Peptides detected in the venom include NPs, BPPs, and inhibitors of SVSPs and SVMPs. Of particular interest, a transcript coding for a protein similar to P-III SVMPs but lacking the MP domain was also found at the protein level in the venom. The existence of such proteins, also supported by finding similar venom gland transcripts in related snake species, has been demonstrated for the first time, justifying the proposal of a new P-IIIe subclass of ancestral SVMP precursor-derived proteins.

Upregulation of interleukin-8 and activin A induces osteoclastogenesis in ameloblastoma

Ameloblastoma is a common odontogenic benign tumor located in the jaws and is characterized by severe local bone destruction. The current study aimed to investigate the effect of interactions between tumor cells and bone marrow stromal cells (BMSCs) on osteoclast formation in ameloblastoma. The impact of ameloblastoma/BMSC interactions on cytokine production, gene expression and osteoclastogenesis was examined using an immortalized ameloblastoma cell line that the authors' previously established. The results demonstrated that interactions between ameloblastoma cells and BMSCs increased interleukin (IL)‑8 and activin A secretion by BMSCs. IL‑8 expression in BMSCs was modulated by tumor‑derived tumor necrosis factor‑α and IL‑8 contributed to osteoclast formation not only directly but also by stimulating receptor activator of NF‑κB ligand (RANKL) expression in BMSCs. Activin A secretion in BMSCs was stimulated by ameloblastoma cells via cell‑to‑cell‑mediated activation of c‑Jun N‑terminal kinase activation, acting as a cofactor of RANKL to induce osteoclast formation and function. The present study highlights the critical role of communication between BMSCs and ameloblastoma cells in bone resorption in ameloblastoma.

Transcriptomic Characterization of the South American Freshwater Stingray Potamotrygon motoro Venom Apparatus

Venomous animals are found through a wide taxonomic range including cartilaginous fish such as the freshwater stingray Potamotrygon motoro occurring in South America, which can injure people and cause venom-related symptoms. Ensuring the efficacy of drug development to treat stingray injuries can be assisted by the knowledge of the venom composition. Here we performed a detailed transcriptomic characterization of the venom gland of the South American freshwater stingray Potamotrygon motoro. The transcripts retrieved showed 418 hits to venom components (comparably to 426 and 396 hits in other two Potamotrygon species), with high expression levels of hyaluronidase, cystatin and calglandulin along with hits uniquely found in P. motoro such as DELTA-alicitoxin-Pse1b, Augerpeptide hhe53 and PI-actitoxin-Aeq3a. We also identified undescribed molecules with extremely high expression values with sequence similarity to the SE-cephalotoxin and Rapunzel genes. Comparative analyses showed that despite being closely related, there may be significant variation among the venoms of freshwater stingrays, highlighting the importance of considering elicit care in handling different envenomation cases. Since hyaluronidase represents a major component of fish venom, we have performed phylogenetic and selective pressure analyses of this gene/protein across all fish with the available information. Results indicated an independent recruitment of the hyaluronidase into the stingray venom relative to that of venomous bony fish. The hyaluronidase residues were found to be mostly under negative selection, but 18 sites showed evidence of diversifying positive selection (P < 0.05). Our data provides new insight into stingray venom variation, composition, and selective pressure in hyaluronidase.

Dynamic Interplay of Innate and Adaptive Immunity During Sterile Retinal Inflammation: Insights From the Transcriptome

The pathogenesis of many retinal degenerations, such as age-related macular degeneration (AMD), is punctuated by an ill-defined network of sterile inflammatory responses. The delineation of innate and adaptive immune milieu among the broad leukocyte infiltrate, and the gene networks, which construct these responses, are poorly described in the eye. Using photo-oxidative damage in a rodent model of subretinal inflammation, we employed a novel RNA-sequencing framework to map the global gene network signature of retinal leukocytes. This revealed a previously uncharted interplay of adaptive immunity during subretinal inflammation, including prolonged enrichment of myeloid and lymphocyte migration, antigen presentation, and the alternative arm of the complement cascade involving Factor B. We demonstrate Factor B-deficient mice are protected against macrophage infiltration and subretinal inflammation. Suppressing the drivers of retinal leukocyte proliferation, or their capacity to elicit complement responses, may help preserve retinal structure and function during sterile inflammation in diseases such as AMD.

Anti-Inflammatory and Immune Regulatory Actions of Naja naja atra Venom

Naja naja atra venom (NNAV) is composed of various proteins, peptides, and enzymes with different biological and pharmacological functions. A number of previous studies have reported that NNAV exerts potent analgesic effects on various animal models of pain. The clinical studies using whole venom or active components have confirmed that NNAV is an effective and safe medicine for treatment of chronic pain. Furthermore, recent studies have demonstrated that NNAV has anti-inflammatory and immune regulatory actions in vitro and in vivo. In this review article, we summarize recent studies of NNAV and its components on inflammation and immunity. The main new findings in NNAV research show that it may enhance innate and humoral immune responses while suppressing T lymphocytes-mediated cellular immunity, thus suggesting that NNAV and its active components may have therapeutic values in the treatment of inflammatory and autoimmune diseases.

Snake Genome Sequencing: Results and Future Prospects

Snake genome sequencing is in its infancy—very much behind the progress made in sequencing the genomes of humans, model organisms and pathogens relevant to biomedical research, and agricultural species. We provide here an overview of some of the snake genome projects in progress, and discuss the biological findings, with special emphasis on toxinology, from the small number of draft snake genomes already published. We discuss the future of snake genomics, pointing out that new sequencing technologies will help overcome the problem of repetitive sequences in assembling snake genomes. Genome sequences are also likely to be valuable in examining the clustering of toxin genes on the chromosomes, in designing recombinant antivenoms and in studying the epigenetic regulation of toxin gene expression.

Clusters of Multiple Mutations: Incidence and Molecular Mechanisms

It has been long understood that mutation distribution is not completely random across genomic space and in time. Indeed, recent surprising discoveries identified multiple simultaneous mutations occurring in tiny regions within chromosomes while the rest of the genome remains relatively mutation-free. Mechanistic elucidation of these phenomena, called mutation showers, mutation clusters, or kataegis, in parallel with findings of abundant clustered mutagenesis in cancer genomes, is ongoing. So far, the combination of factors most important for clustered mutagenesis is the induction of DNA lesions within unusually long and persistent single-strand DNA intermediates. In addition to being a fascinating phenomenon, clustered mutagenesis also became an indispensable tool for identifying a previously unrecognized major source of mutation in cancer, APOBEC cytidine deaminases. Future research on clustered mutagenesis may shed light onto important mechanistic details of genome maintenance, with potentially profound implications for human health.

Snake venomics of Micrurus alleni and Micrurus mosquitensis from the Caribbean region of Costa Rica reveals two divergent compositional patterns in New World elapids

Protein composition, toxicity, and neutralization of the venoms of Micrurus alleni and Micrurus mosquitensis, two sympatric monadal coral snakes found in humid environments of the Caribbean region of Costa Rica, were studied. Proteomic profiling revealed that these venoms display highly divergent compositions: the former dominated by three-finger toxins (3FTx) and the latter by phospholipases A2 (PLA2). Protein family abundances correlated with enzymatic and toxic characteristics of the venoms. Selective inhibition experiments showed that PLA2s play only a marginal role in the lethal effect of M. alleni venom, but have a major role in M. mosquitensis venom. Proteomic data gathered from other Micrurus species evidenced that the two divergent venom phenotypes are recurrent, and may constitute a general trend across New World elapids. Further, M. mosquitensis, but not M. alleni, venom contains PLA2-like/Kunitz-type inhibitor complex(es) that resemble the ASIC1a/2-activating MitTx heterodimeric toxin isolated from Micrurus tener venom. The evolutionary origin and adaptive relevance of the puzzling phenotypic variability of Micrurus venoms remain to be understood. An antivenom against the PLA2-predominant Micrurus nigrocinctus venom strongly cross-recognized and neutralized M. mosquitensis venom, but only weakly M. alleni venom.

Restriction and recruitment-gene duplication and the origin and evolution of snake venom toxins

Snake venom has been hypothesized to have originated and diversified through a process that involves duplication of genes encoding body proteins with subsequent recruitment of the copy to the venom gland, where natural selection acts to develop or increase toxicity. However, gene duplication is known to be a rare event in vertebrate genomes, and the recruitment of duplicated genes to a novel expression domain (neofunctionalization) is an even rarer process that requires the evolution of novel combinations of transcription factor binding sites in upstream regulatory regions. Therefore, although this hypothesis concerning the evolution of snake venom is very unlikely and should be regarded with caution, it is nonetheless often assumed to be established fact, hindering research into the true origins of snake venom toxins. To critically evaluate this hypothesis, we have generated transcriptomic data for body tissues and salivary and venom glands from five species of venomous and nonvenomous reptiles. Our comparative transcriptomic analysis of these data reveals that snake venom does not evolve through the hypothesized process of duplication and recruitment of genes encoding body proteins. Indeed, our results show that many proposed venom toxins are in fact expressed in a wide variety of body tissues, including the salivary gland of nonvenomous reptiles and that these genes have therefore been restricted to the venom gland following duplication, not recruited. Thus, snake venom evolves through the duplication and subfunctionalization of genes encoding existing salivary proteins. These results highlight the danger of the elegant and intuitive “just-so story” in evolutionary biology.

Evolution of an ancient venom: recognition of a novel family of cnidarian toxins and the common evolutionary origin of sodium and potassium neurotoxins in sea anemone

Despite Cnidaria (sea anemones, corals, jellyfish, and hydroids) being the oldest venomous animal lineage, structure-function relationships, phyletic distributions, and the molecular evolutionary regimes of toxins encoded by these intriguing animals are poorly understood. Hence, we have comprehensively elucidated the phylogenetic and molecular evolutionary histories of pharmacologically characterized cnidarian toxin families, including peptide neurotoxins (voltage-gated Na(+) and K(+) channel-targeting toxins: NaTxs and KTxs, respectively), pore-forming toxins (actinoporins, aerolysin-related toxins, and jellyfish toxins), and the newly discovered small cysteine-rich peptides (SCRiPs). We show that despite long evolutionary histories, most cnidarian toxins remain conserved under the strong influence of negative selection-a finding that is in striking contrast to the rapid evolution of toxin families in evolutionarily younger lineages, such as cone snails and advanced snakes. In contrast to the previous suggestions that implicated SCRiPs in the biomineralization process in corals, we demonstrate that they are potent neurotoxins that are likely involved in the envenoming function, and thus represent the first family of neurotoxins from corals. We also demonstrate the common evolutionary origin of type III KTxs and NaTxs in sea anemones. We show that type III KTxs have evolved from NaTxs under the regime of positive selection, and likely represent a unique evolutionary innovation of the Actinioidea lineage. We report a correlation between the accumulation of episodically adaptive sites and the emergence of novel pharmacological activities in this rapidly evolving neurotoxic clade.

Testing the Toxicofera: comparative transcriptomics casts doubt on the single, early evolution of the reptile venom system

The identification of apparently conserved gene complements in the venom and salivary glands of a diverse set of reptiles led to the development of the Toxicofera hypothesis - the single, early evolution of the venom system in reptiles. However, this hypothesis is based largely on relatively small scale EST-based studies of only venom or salivary glands and toxic effects have been assigned to only some putative Toxicoferan toxins in some species. We set out to examine the distribution of these proposed venom toxin transcripts in order to investigate to what extent conservation of gene complements may reflect a bias in previous sampling efforts. Our quantitative transcriptomic analyses of venom and salivary glands and other body tissues in five species of reptile, together with the use of available RNA-Seq datasets for additional species, shows that the majority of genes used to support the establishment and expansion of the Toxicofera are in fact expressed in multiple body tissues and most likely represent general maintenance or "housekeeping" genes. The apparent conservation of gene complements across the Toxicofera therefore reflects an artefact of incomplete tissue sampling. We therefore conclude that venom has evolved multiple times in reptiles.

IMPACT_S: integrated multiprogram platform to analyze and combine tests of selection

Among the major goals of research in evolutionary biology are the identification of genes targeted by natural selection and understanding how various regimes of evolution affect the fitness of an organism. In particular, adaptive evolution enables organisms to adapt to changing ecological factors such as diet, temperature, habitat, predatory pressures and prey abundance. An integrative approach is crucial for the identification of non-synonymous mutations that introduce radical changes in protein biochemistry and thus in turn influence the structure and function of proteins. Performing such analyses manually is often a time-consuming process, due to the large number of statistical files generated from multiple approaches, especially when assessing numerous taxa and/or large datasets. We present IMPACT_S, an easy-to-use Graphical User Interface (GUI) software, which rapidly and effectively integrates, filters and combines results from three widely used programs for assessing the influence of selection: Codeml (PAML package), Datamonkey and TreeSAAP. It enables the identification and tabulation of sites detected by these programs as evolving under the influence of positive, neutral and/or negative selection in protein-coding genes. IMPACT_S further facilitates the automatic mapping of these sites onto the three-dimensional structures of proteins. Other useful tools incorporated in IMPACT_S include Jmol, Archaeopteryx, Gnuplot, PhyML, a built-in Swiss-Model interface and a PDB downloader. The relevance and functionality of IMPACT_S is shown through a case study on the toxicoferan-reptilian Cysteine-rich Secretory Proteins (CRiSPs). IMPACT_S is a platform-independent software released under GPLv3 license, freely available online from http://impact-s.sourceforge.net.

Extreme venom variation in Middle Eastern vipers: a proteomics comparison of Eristicophis macmahonii, Pseudocerastes fieldi and Pseudocerastes persicus

Venoms of the viperid sister genera Eristicophis and Pseudocerastes are poorly studied despite their anecdotal reputation for producing severe or even lethal envenomations. This is due in part to the remote and politically unstable regions that they occupy. All species contained are sit and wait ambush feeders. Thus, this study examined their venoms through proteomics techniques in order to establish if this feeding ecology, and putatively low levels of gene flow, have resulted in significant variations in venom profile. The techniques indeed revealed extreme venom variation. This has immediate implications as only one antivenom is made (using the venom of Pseudocerastes persicus) yet the proteomic variation suggests that it would be of only limited use for the other species, even the sister species Pseudocerastes fieldi. The high degree of variation however also points toward these species being rich resources for novel compounds which may have use as lead molecules in drug design and development.

BIOLOGICAL SIGNIFICANCE

These results show extreme venom variation between these closely related snakes. These results have direct implications for the treatment of the envenomed patient.

Diversification of a single ancestral gene into a successful toxin superfamily in highly venomous Australian funnel-web spiders

BACKGROUND

Spiders have evolved pharmacologically complex venoms that serve to rapidly subdue prey and deter predators. The major toxic factors in most spider venoms are small, disulfide-rich peptides. While there is abundant evidence that snake venoms evolved by recruitment of genes encoding normal body proteins followed by extensive gene duplication accompanied by explosive structural and functional diversification, the evolutionary trajectory of spider-venom peptides is less clear.

RESULTS

Here we present evidence of a spider-toxin superfamily encoding a high degree of sequence and functional diversity that has evolved via accelerated duplication and diversification of a single ancestral gene. The peptides within this toxin superfamily are translated as prepropeptides that are posttranslationally processed to yield the mature toxin. The N-terminal signal sequence, as well as the protease recognition site at the junction of the propeptide and mature toxin are conserved, whereas the remainder of the propeptide and mature toxin sequences are variable. All toxin transcripts within this superfamily exhibit a striking cysteine codon bias. We show that different pharmacological classes of toxins within this peptide superfamily evolved under different evolutionary selection pressures.

CONCLUSIONS

Overall, this study reinforces the hypothesis that spiders use a combinatorial peptide library strategy to evolve a complex cocktail of peptide toxins that target neuronal receptors and ion channels in prey and predators. We show that the ω-hexatoxins that target insect voltage-gated calcium channels evolved under the influence of positive Darwinian selection in an episodic fashion, whereas the κ-hexatoxins that target insect calcium-activated potassium channels appear to be under negative selection. A majority of the diversifying sites in the ω-hexatoxins are concentrated on the molecular surface of the toxins, thereby facilitating neofunctionalisation leading to new toxin pharmacology.

Venom down under: dynamic evolution of Australian elapid snake toxins

Despite the unparalleled diversity of venomous snakes in Australia, research has concentrated on a handful of medically significant species and even of these very few toxins have been fully sequenced. In this study, venom gland transcriptomes were sequenced from eleven species of small Australian elapid snakes, from eleven genera, spanning a broad phylogenetic range. The particularly large number of sequences obtained for three-finger toxin (3FTx) peptides allowed for robust reconstructions of their dynamic molecular evolutionary histories. We demonstrated that each species preferentially favoured different types of α-neurotoxic 3FTx, probably as a result of differing feeding ecologies. The three forms of α-neurotoxin [Type I (also known as (aka): short-chain), Type II (aka: long-chain) and Type III] not only adopted differential rates of evolution, but have also conserved a diversity of residues, presumably to potentiate prey-specific toxicity. Despite these differences, the different α-neurotoxin types were shown to accumulate mutations in similar regions of the protein, largely in the loops and structurally unimportant regions, highlighting the significant role of focal mutagenesis. We theorize that this phenomenon not only affects toxin potency or specificity, but also generates necessary variation for preventing/delaying prey animals from acquiring venom-resistance. This study also recovered the first full-length sequences for multimeric phospholipase A₂ (PLA₂) ‘taipoxin/paradoxin’ subunits from non-Oxyuranus species, confirming the early recruitment of this extremely potent neurotoxin complex to the venom arsenal of Australian elapid snakes. We also recovered the first natriuretic peptides from an elapid that lack the derived C-terminal tail and resemble the plesiotypic form (ancestral character state) found in viper venoms. This provides supporting evidence for a single early recruitment of natriuretic peptides into snake venoms. Novel forms of kunitz and waprin peptides were recovered, including dual domain kunitz-kunitz precursors and the first kunitz-waprin hybrid precursors from elapid snakes. The novel sequences recovered in this study reveal that the huge diversity of unstudied venomous Australian snakes are of considerable interest not only for the investigation of venom and whole organism evolution but also represent an untapped bioresource in the search for novel compounds for use in drug design and development.

Evolution stings: the origin and diversification of scorpion toxin peptide scaffolds

The episodic nature of natural selection and the accumulation of extreme sequence divergence in venom-encoding genes over long periods of evolutionary time can obscure the signature of positive Darwinian selection. Recognition of the true biocomplexity is further hampered by the limited taxon selection, with easy to obtain or medically important species typically being the subject of intense venom research, relative to the actual taxonomical diversity in nature. This holds true for scorpions, which are one of the most ancient terrestrial venomous animal lineages. The family Buthidae that includes all the medically significant species has been intensely investigated around the globe, while almost completely ignoring the remaining non-buthid families. Australian scorpion lineages, for instance, have been completely neglected, with only a single scorpion species (Urodacus yaschenkoi) having its venom transcriptome sequenced. Hence, the lack of venom composition and toxin sequence information from an entire continent’s worth of scorpions has impeded our understanding of the molecular evolution of scorpion venom. The molecular origin, phylogenetic relationships and evolutionary histories of most scorpion toxin scaffolds remain enigmatic. In this study, we have sequenced venom gland transcriptomes of a wide taxonomical diversity of scorpions from Australia, including buthid and non-buthid representatives. Using state-of-art molecular evolutionary analyses, we show that a majority of CSα/β toxin scaffolds have experienced episodic influence of positive selection, while most non-CSα/β linear toxins evolve under the extreme influence of negative selection. For the first time, we have unraveled the molecular origin of the major scorpion toxin scaffolds, such as scorpion venom single von Willebrand factor C-domain peptides (SV-SVC), inhibitor cystine knot (ICK), disulphide-directed beta-hairpin (DDH), bradykinin potentiating peptides (BPP), linear non-disulphide bridged peptides and antimicrobial peptides (AMP). We have thus demonstrated that even neglected lineages of scorpions are a rich pool of novel biochemical components, which have evolved over millions of years to target specific ion channels in prey animals, and as a result, possess tremendous implications in therapeutics.

Three-fingered RAVERs: Rapid Accumulation of Variations in Exposed Residues of snake venom toxins

Three-finger toxins (3FTx) represent one of the most abundantly secreted and potently toxic components of colubrid (Colubridae), elapid (Elapidae) and psammophid (Psammophiinae subfamily of the Lamprophidae) snake venom arsenal. Despite their conserved structural similarity, they perform a diversity of biological functions. Although they are theorised to undergo adaptive evolution, the underlying diversification mechanisms remain elusive. Here, we report the molecular evolution of different 3FTx functional forms and show that positively selected point mutations have driven the rapid evolution and diversification of 3FTx. These diversification events not only correlate with the evolution of advanced venom delivery systems (VDS) in Caenophidia, but in particular the explosive diversification of the clade subsequent to the evolution of a high pressure, hollow-fanged VDS in elapids, highlighting the significant role of these toxins in the evolution of advanced snakes. We show that Type I, II and III α-neurotoxins have evolved with extreme rapidity under the influence of positive selection. We also show that novel Oxyuranus/Pseudonaja Type II forms lacking the apotypic loop-2 stabilising cysteine doublet characteristic of Type II forms are not phylogenetically basal in relation to other Type IIs as previously thought, but are the result of secondary loss of these apotypic cysteines on at least three separate occasions. Not all 3FTxs have evolved rapidly: κ-neurotoxins, which form non-covalently associated heterodimers, have experienced a relatively weaker influence of diversifying selection; while cytotoxic 3FTx, with their functional sites, dispersed over 40% of the molecular surface, have been extremely constrained by negative selection. We show that the a previous theory of 3FTx molecular evolution (termed ASSET) is evolutionarily implausible and cannot account for the considerable variation observed in very short segments of 3FTx. Instead, we propose a theory of Rapid Accumulation of Variations in Exposed Residues (RAVER) to illustrate the significance of point mutations, guided by focal mutagenesis and positive selection in the evolution and diversification of 3FTx.