Snake Venom Gland Organoids.

Wnt dependency and Lgr5 expression define multiple mammalian epithelial stem cell types. Under defined growth factor conditions, such adult stem cells (ASCs) grow as 3D organoids that recapitulate essential features of the pertinent epithelium. Here, we establish long-term expanding venom gland organoids from several snake species. The newly assembled transcriptome of the Cape coral snake reveals that organoids express high levels of toxin transcripts. Single-cell RNA sequencing of both organoids and primary tissue identifies distinct venom-expressing cell types as well as proliferative cells expressing homologs of known mammalian stem cell markers. A hard-wired regional heterogeneity in the expression of individual venom components is maintained in organoid cultures. Harvested venom peptides reflect crude venom composition and display biological activity. This study extends organoid technology to reptilian tissues and describes an experimentally tractable model system representing the snake venom gland.

Recombinant snake antivenoms get closer to the clinic.

Snakebite envenomings kill ~100 000 victims each year and leave many more with permanent sequelae. Antivenoms have been available for more than 125 years but are in need of innovation. A new study by Khalek et al. highlights broadly neutralizing human monoclonal antibodies (mAbs) that might be used to develop recombinant antivenoms with superior therapeutic benefits.

The genetic regulatory architecture and epigenomic basis for age-related changes in rattlesnake venom.

Developmental phenotypic changes can evolve under selection imposed by age- and size-related ecological differences. Many of these changes occur through programmed alterations to gene expression patterns, but the molecular mechanisms and gene-regulatory networks underlying these adaptive changes remain poorly understood. Many venomous snakes, including the eastern diamondback rattlesnake (Crotalus adamanteus), undergo correlated changes in diet and venom expression as snakes grow larger with age, providing models for identifying mechanisms of timed expression changes that underlie adaptive life history traits. By combining a highly contiguous, chromosome-level genome assembly with measures of expression, chromatin accessibility, and histone modifications, we identified cis-regulatory elements and trans-regulatory factors controlling venom ontogeny in the venom glands of C. adamanteus. Ontogenetic expression changes were significantly correlated with epigenomic changes within genes, immediately adjacent to genes (e.g., promoters), and more distant from genes (e.g., enhancers). We identified 37 candidate transcription factors (TFs), with the vast majority being up-regulated in adults. The ontogenetic change is largely driven by an increase in the expression of TFs associated with growth signaling, transcriptional activation, and circadian rhythm/biological timing systems in adults with corresponding epigenomic changes near the differentially expressed venom genes. However, both expression activation and repression contributed to the composition of both adult and juvenile venoms, demonstrating the complexity and potential evolvability of gene regulation for this trait. Overall, given that age-based trait variation is common across the tree of life, we provide a framework for understanding gene-regulatory-network-driven life-history evolution more broadly.

DeTox: a pipeline for the detection of toxins in venomous organisms.

Venomous organisms have independently evolved the ability to produce toxins 101 times during their evolutionary history, resulting in over 200 000 venomous species. Collectively, these species produce millions of toxins, making them a valuable resource for bioprospecting and understanding the evolutionary mechanisms underlying genetic diversification. RNA-seq is the preferred method for characterizing toxin repertoires, but the analysis of the resulting data remains challenging. While early approaches relied on similarity-based mapping to known toxin databases, recent studies have highlighted the importance of structural features for toxin detection. The few existing pipelines lack an integration between these complementary approaches, and tend to be difficult to run for non-experienced users. To address these issues, we developed DeTox, a comprehensive and user-friendly tool for toxin research. It combines fast execution, parallelization and customization of parameters. DeTox was tested on published transcriptomes from gastropod mollusks, cnidarians and snakes, retrieving most putative toxins from the original articles and identifying additional peptides as potential toxins to be confirmed through manual annotation and eventually proteomic analysis. By integrating a structure-based search with similarity-based approaches, DeTox allows the comprehensive characterization of toxin repertoire in poorly-known taxa. The effect of the taxonomic bias in existing databases is minimized in DeTox, as mirrored in the detection of unique and divergent toxins that would have been overlooked by similarity-based methods. DeTox streamlines toxin annotation, providing a valuable tool for efficient identification of venom components that will enhance venom research in neglected taxa.

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.

Horizontal gene transfer underlies the painful stings of asp caterpillars (Lepidoptera: Megalopygidae).

Larvae of the genus Megalopyge (Lepidoptera: Zygaenoidea: Megalopygidae), known as asp or puss caterpillars, produce defensive venoms that cause severe pain. Here, we present the anatomy, chemistry, and mode of action of the venom systems of caterpillars of two megalopygid species, the Southern flannel moth Megalopyge opercularis and the black-waved flannel moth Megalopyge crispata. We show that megalopygid venom is produced in secretory cells that lie beneath the cuticle and are connected to the venom spines by canals. Megalopygid venoms consist of large aerolysin-like pore-forming toxins, which we have named megalysins, and a small number of peptides. The venom system differs markedly from those of previously studied venomous zygaenoids of the family Limacodidae, suggestive of an independent origin. Megalopygid venom potently activates mammalian sensory neurons via membrane permeabilization and induces sustained spontaneous pain behavior and paw swelling in mice. These bioactivities are ablated by treatment with heat, organic solvents, or proteases, indicating that they are mediated by larger proteins such as the megalysins. We show that the megalysins were recruited as venom toxins in the Megalopygidae following horizontal transfer of genes from bacteria to the ancestors of ditrysian Lepidoptera. Megalopygids have recruited aerolysin-like proteins as venom toxins convergently with centipedes, cnidarians, and fish. This study highlights the role of horizontal gene transfer in venom evolution.

Mast cell-derived prostaglandin D2 limits the subcutaneous absorption of honey bee venom in mice.

Mast cells play pivotal roles in innate host defenses against venom. Activated mast cells release large amounts of prostaglandin D2 (PGD2). However, the role of PGD2 in such host defense remains unclear. We found that c-kit-dependent and c-kit-independent mast cell-specific hematopoietic prostaglandin D synthase (H-pgds) deficiency significantly exacerbated honey bee venom (BV)-induced hypothermia and increased mortality rates in mice. BV absorption via postcapillary venules in the skin was accelerated upon endothelial barrier disruption resulting in increased plasma venom concentrations. These results suggest that mast cell-derived PGD2 may enhance host defense against BV and save lives by inhibiting BV absorption into circulation.

Peptide toxins that target vertebrate voltage-gated sodium channels underly the painful stings of harvester ants.

Harvester ants (genus Pogonomyrmex) are renowned for their stings which cause intense, long-lasting pain, and other neurotoxic symptoms in vertebrates. Here, we show that harvester ant venoms are relatively simple and composed largely of peptide toxins. One class of peptides is primarily responsible for the long-lasting local pain of envenomation via activation of peripheral sensory neurons. These hydrophobic, cysteine-free peptides potently modulate mammalian voltage-gated sodium (NaV) channels, reducing the voltage threshold for activation and inhibiting channel inactivation. These toxins appear to have evolved specifically to deter vertebrates.

Prey Shifts Drive Venom Evolution in Cone Snails.

Venom systems are complex traits that have independently emerged multiple times in diverse plant and animal phyla. Within each venomous lineage there typically exists interspecific variation in venom composition where several factors have been proposed as drivers of variation, including phylogeny and diet. Understanding these factors is of broad biological interest and has implications for the development of antivenom therapies and venom-based drug discovery. Because of their high species richness and the presence of several major evolutionary prey shifts, venomous marine cone snails (genus Conus) provide an ideal system to investigate drivers of interspecific venom variation. Here, by analyzing the venom gland expression profiles of ∼3,000 toxin genes from 42 species of cone snail, we elucidate the role of prey-specific selection pressures in shaping venom variation. By analyzing overall venom composition and individual toxin structures, we demonstrate that the shifts from vermivory to piscivory in Conus are complemented by distinct changes in venom composition independent of phylogeny. In vivo injections of venom from piscivorous cone snails in fish further showed a higher potency compared with venom of nonpiscivores demonstrating a selective advantage. Together, our findings provide compelling evidence for the role of prey shifts in directing the venom composition of cone snails and expand our understanding of the mechanisms of venom variation and diversification.

Sea Anemone Membrane Attack Complex/Perforin Superfamily Demonstrates an Evolutionary Transitional State between Venomous and Developmental Functions.

Gene duplication is a major force driving evolutionary innovation. A classic example is generating new animal toxins via duplication of physiological protein-encoding genes and recruitment into venom. While this process drives the innovation of many animal venoms, reverse recruitment of toxins into nonvenomous cells remains unresolved. Using comparative genomics, we find members of the Membrane Attack Complex and Perforin Family (MAC) have been recruited into venom-injecting cells (cnidocytes), in soft and stony corals and sea anemones, suggesting that the ancestral MAC was a cnidocyte expressed toxin. Further investigation into the model sea anemone Nematostella vectensis reveals that three members have undergone Nematostella-specific duplications leading to their reverse recruitment into endomesodermal cells. Furthermore, simultaneous knockdown of all three endomesodermally expressed MACs leads to mis-development, supporting that these paralogs have nonvenomous function. By resolving the evolutionary history and function of MACs in Nematostella, we provide the first proof for reverse recruitment from venom to organismal development.

Whole transcriptome sequencing reveals the activity of the PLA2 family members in Androctonus crassicauda (Scorpionida: Buthidae) venom gland.

Phospholipase A2 is the most abundant venom gland enzyme, whose activity leads to the activation of the inflammatory response by accumulating lipid mediators. This study aimed to identify, classify, and investigate the properties of venom PLA2 isoforms. Then, the present findings were confirmed by chemically measuring the activity of PLA2. The sequences representing PLA2 annotation were extracted from the Androctonus crassicauda transcriptome dataset using BLAS searches against the local PLA2 database. We found several cDNA sequences of PLA2 classified and named by conducting multiple searches as platelet-activating factor acetylhydrolases, calcium-dependent PLA2s, calcium-independent PLA2s, and secreted PLA2s. The largest and smallest isoforms of these proteins range between approximately 70.34 kDa (iPLA2) and 17.75 kDa (cPLA2). Among sPLA2 isoforms, sPLA2GXIIA and sPLA2G3 with ORF encoding 169 and 299 amino acids are the smallest and largest secreted PLA2, respectively. These results collectively suggested that A. crassicauda venom has PLA2 activity, and the members of this protein family may have important biological roles in lipid metabolism. This study also revealed the interaction between members of PLA2s in the PPI network. The results of this study would greatly help with the classification, evolutionary relationships, and interactions between PLA2 family proteins in the gene network.

Selective block of human Kv1.1 channels and an epilepsy-associated gain-of-function mutation by AETX-K peptide.

Dysfunction of the human voltage-gated K+ channel Kv1.1 has been associated with epilepsy, multiple sclerosis, episodic ataxia, myokymia, and cardiorespiratory dysregulation. We report here that AETX-K, a sea anemone type I (SAK1) peptide toxin we isolated from a phage display library, blocks Kv1.1 with high affinity (Ki ~ 1.6 pM) and notable specificity, inhibiting other Kv channels we tested a million-fold less well. Nuclear magnetic resonance (NMR) was employed both to determine the three-dimensional structure of AETX-K, showing it to employ a classic SAK1 scaffold while exhibiting a unique electrostatic potential surface, and to visualize AETX-K bound to the Kv1.1 pore domain embedded in lipoprotein nanodiscs. Study of Kv1.1 in Xenopus oocytes with AETX-K and point variants using electrophysiology demonstrated the blocking mechanism to employ a toxin-channel configuration we have described before whereby AETX-K Lys23 , two positions away on the toxin interaction surface from the classical blocking residue, enters the pore deeply enough to interact with K+ ions traversing the pathway from the opposite side of the membrane. The mutant channel Kv1.1-L296 F is associated with pharmaco-resistant multifocal epilepsy in infants because it significantly increases K+ currents by facilitating opening and slowing closure of the channels. Consistent with the therapeutic potential of AETX-K for Kv1.1 gain-of-function-associated diseases, AETX-K at 4 pM decreased Kv1.1-L296 F currents to wild-type levels; further, populations of heteromeric channels formed by co-expression Kv1.1 and Kv1.2, as found in many neurons, showed a Ki of ~10 nM even though homomeric Kv1.2 channels were insensitive to the toxin (Ki > 2000 nM).

Revealing molecular determinants governing mambalgin-3 pharmacology at acid-sensing ion channel 1 variants.

Acid-sensing ion channels (ASICs) are trimeric proton-gated cation channels that play a role in neurotransmission and pain sensation. The snake venom-derived peptides, mambalgins, exhibit potent analgesic effects in rodents by inhibiting central ASIC1a and peripheral ASIC1b. Despite their distinct species- and subtype-dependent pharmacology, previous structure-function studies have focussed on the mambalgin interaction with ASIC1a. Currently, the specific channel residues responsible for this pharmacological profile, and the mambalgin pharmacophore at ASIC1b remain unknown. Here we identify non-conserved residues at the ASIC1 subunit interface that drive differences in the mambalgin pharmacology from rat ASIC1a to ASIC1b, some of which likely do not make peptide binding interactions. Additionally, an amino acid variation below the core binding site explains potency differences between rat and human ASIC1. Two regions within the palm domain, which contribute to subtype-dependent effects for mambalgins, play key roles in ASIC gating, consistent with subtype-specific differences in the peptides mechanism. Lastly, there is a shared primary mambalgin pharmacophore for ASIC1a and ASIC1b activity, with certain peripheral peptide residues showing variant-specific significance for potency. Through our broad mutagenesis studies across various species and subtype variants, we gain a more comprehensive understanding of the pharmacophore and the intricate molecular interactions that underlie ligand specificity. These insights pave the way for the development of more potent and targeted peptide analogues required to advance our understating of human ASIC1 function and its role in disease.

Exploring oak processionary caterpillar induced lepidopterism (Part 1): unveiling molecular insights through transcriptomics and proteomics.

Lepidopterism, a skin inflammation condition caused by direct or airborne exposure to irritating hairs (setae) from processionary caterpillars, is becoming a significant public health concern. Recent outbreaks of the oak processionary caterpillar (Thaumetopoea processionea) have caused noteworthy health and economic consequences, with a rising frequency expected in the future, exacerbated by global warming promoting the survival of the caterpillar. Current medical treatments focus on symptom relief due to the lack of an effective therapy. While the source is known, understanding the precise causes of symptoms remain incomplete understood. In this study, we employed an advanced method to extract venom from the setae and identify the venom components through high-quality de novo transcriptomics, venom proteomics, and bioinformatic analysis. A total of 171 venom components were identified, including allergens, odorant binding proteins, small peptides, enzymes, enzyme inhibitors, and chitin biosynthesis products, potentially responsible for inflammatory and allergic reactions. This work presents the first comprehensive proteotranscriptomic database of T. processionea, contributing to understanding the complexity of lepidopterism. Furthermore, these findings hold promise for advancing therapeutic approaches to mitigate the global health impact of T. processionea and related caterpillars.

A novel broad spectrum venom metalloproteinase autoinhibitor in the rattlesnake Crotalus atrox evolved via a shift in paralog function.

The complexity of snake venom composition reflects adaptation to the diversity of prey and may be driven at times by a coevolutionary arms race between snakes and venom-resistant prey. However, many snakes are also resistant to their own venom due to serum-borne inhibitors of venom toxins, which raises the question of how snake autoinhibitors maintain their efficacy as venom proteins evolve. To investigate this potential three-way arms race among venom, prey, and autoinhibitors, we have identified and traced the evolutionary origin of serum inhibitors of snake venom metalloproteinases (SVMPs) in the Western Diamondback rattlesnake Crotalus atrox which possesses the largest known battery of SVMP genes among crotalids examined. We found that C. atrox expresses five members of a Fetuin A-related metalloproteinase inhibitor family but that one family member, FETUA-3, is the major SVMP inhibitor that binds to approximately 20 different C. atrox SVMPs and inhibits activities of all three SVMP classes. We show that the fetua-3 gene arose deep within crotalid evolution before the origin of New World species but, surprisingly, fetua-3 belongs to a different paralog group than previously identified SVMP inhibitors in Asian and South American crotalids. Conversely, the C. atrox FETUA-2 ortholog of previously characterized crotalid SVMP inhibitors shows limited activity against C. atrox SVMPs. These results reveal that there has been a functional evolutionary shift in the major SVMP inhibitor in the C. atrox lineage as the SVMP family expanded and diversified in the Crotalus lineage. This broad-spectrum inhibitor may be of potential therapeutic interest.

Convergent evolution of venom gland transcriptomes across Metazoa.

Animals have repeatedly evolved specialized organs and anatomical structures to produce and deliver a mixture of potent bioactive molecules to subdue prey or predators-venom. This makes it one of the most widespread, convergent functions in the animal kingdom. Whether animals have adopted the same genetic toolkit to evolved venom systems is a fascinating question that still eludes us. Here, we performed a comparative analysis of venom gland transcriptomes from 20 venomous species spanning the main Metazoan lineages to test whether different animals have independently adopted similar molecular mechanisms to perform the same function. We found a strong convergence in gene expression profiles, with venom glands being more similar to each other than to any other tissue from the same species, and their differences closely mirroring the species phylogeny. Although venom glands secrete some of the fastest evolving molecules (toxins), their gene expression does not evolve faster than evolutionarily older tissues. We found 15 venom gland-specific gene modules enriched in endoplasmic reticulum stress and unfolded protein response pathways, indicating that animals have independently adopted stress response mechanisms to cope with mass production of toxins. This, in turn, activates regulatory networks for epithelial development, cell turnover, and maintenance, which seem composed of both convergent and lineage-specific factors, possibly reflecting the different developmental origins of venom glands. This study represents a first step toward an understanding of the molecular mechanisms underlying the repeated evolution of one of the most successful adaptive traits in the animal kingdom.

Host and venom evolution in parasitoid wasps: does independently adapting to the same host shape the evolution of the venom gland transcriptome?

BACKGROUND: Venoms have repeatedly evolved over 100 occasions throughout the animal tree of life, making them excellent systems for exploring convergent evolutionary novelty. Growing evidence supports that venom evolution is predominantly driven by prey or host-related selection pressures, and the expression patterns of venom glands reflect adaptive evolution. However, it remains elusive whether the evolution of expression patterns in venom glands is likewise a convergent evolution driven by their prey/host species. RESULTS: We utilized parasitoid wasps that had independently adapted to Drosophila hosts as models to investigate the convergent evolution of venom gland transcriptomes in 19 hymenopteran species spanning ~ 200 million years of evolution. Comparative transcriptome analysis reveals that the global expression patterns among the venom glands of Drosophila parasitoid wasps do not achieve higher similarity compared to non-Drosophila parasitoid wasps. Further evolutionary analyses of expression patterns at the single gene, orthogroup, and Gene Ontology (GO) term levels indicate that some orthogroups/GO terms show correlation with the Drosophila parasitoid wasps. However, these groups rarely include genes highly expressed in venom glands or putative venom genes in the Drosophila parasitoid wasps. CONCLUSIONS: Our study suggests that convergent evolution may not play a predominant force shaping gene expression levels in the venom gland of the Drosophila parasitoid wasps, offering novel insights into the co-evolution between venom and prey/host.

The venom and telopodal defence systems of the centipede Lithobius forficatus are functionally convergent serial homologues.

BACKGROUND: Evolution of novelty is a central theme in evolutionary biology, yet studying the origins of traits with an apparently discontinuous origin remains a major challenge. Venom systems are a well-suited model for the study of this phenomenon because they capture several aspects of novelty across multiple levels of biological complexity. However, while there is some knowledge on the evolution of individual toxins, not much is known about the evolution of venom systems as a whole. One way of shedding light on the evolution of new traits is to investigate less specialised serial homologues, i.e. repeated traits in an organism that share a developmental origin. This approach can be particularly informative in animals with repetitive body segments, such as centipedes. RESULTS: Here, we investigate morphological and biochemical aspects of the defensive telopodal glandular organs borne on the posterior legs of venomous stone centipedes (Lithobiomorpha), using a multimethod approach, including behavioural observations, comparative morphology, proteomics, comparative transcriptomics and molecular phylogenetics. We show that the anterior venom system and posterior telopodal defence system are functionally convergent serial homologues, where one (telopodal defence) represents a model for the putative early evolutionary state of the other (venom). Venom glands and telopodal glandular organs appear to have evolved from the same type of epidermal gland (four-cell recto-canal type) and while the telopodal defensive secretion shares a great degree of compositional overlap with centipede venoms in general, these similarities arose predominantly through convergent recruitment of distantly related toxin-like components. Both systems are composed of elements predisposed to functional innovation across levels of biological complexity that range from proteins to glands, demonstrating clear parallels between molecular and morphological traits in the properties that facilitate the evolution of novelty. CONCLUSIONS: The evolution of the lithobiomorph telopodal defence system provides indirect empirical support for the plausibility of the hypothesised evolutionary origin of the centipede venom system, which occurred through functional innovation and gradual specialisation of existing epidermal glands. Our results thus exemplify how continuous transformation and functional innovation can drive the apparent discontinuous emergence of novelties on higher levels of biological complexity.

From birth to bite: the evolutionary ecology of India's medically most important snake venoms.

BACKGROUND: Snake venoms can exhibit remarkable inter- and intraspecific variation. While diverse ecological and environmental factors are theorised to explain this variation, only a handful of studies have attempted to unravel their precise roles. This knowledge gap not only impedes our understanding of venom evolution but may also have dire consequences on snakebite treatment. To address this shortcoming, we investigated the evolutionary ecology of venoms of Russell's viper (Daboia russelii) and spectacled cobra (Naja naja), India's two clinically most important snakes responsible for an alarming number of human deaths and disabilities. METHODOLOGY: Several individuals (n = 226) of D. russelii and N. naja belonging to multiple clutches (n = 9) and their mothers were maintained in captivity to source ontogenetic stage-specific venoms. Using various in vitro and in vivo assays, we assessed the significance of prey, ontogeny and sex in driving venom composition, function, and potency. RESULTS: Considerable ontogenetic shifts in venom profiles were observed in D. russelii, with the venoms of newborns being many times as potent as juveniles and adults against mammalian (2.3-2.5 ×) and reptilian (2-10 ×) prey. This is the first documentation of the ontogenetic shift in viperine snakes. In stark contrast, N. naja, which shares a biogeographic distribution similar to D. russelii, deployed identical biochemical cocktails across development. Furthermore, the binding kinetics of cobra venom toxins against synthetic target receptors from various prey and predators shed light on the evolutionary arms race. CONCLUSIONS: Our findings, therefore, provide fascinating insights into the roles of ecology and life history traits in shaping snake venoms.

Venomics and Peptidomics of Palearctic Vipers: A Clade-Wide Analysis of Seven Taxa of the Genera Vipera, Montivipera, Macrovipera, and Daboia across Türkiye.

Snake venom variations are a crucial factor to understand the consequences of snakebite envenoming worldwide, and therefore it is important to know about toxin composition alterations between taxa. Palearctic vipers of the genera Vipera, Montivipera, Macrovipera, and Daboia have high medical impacts across the Old World. One hotspot for their occurrence and diversity is Türkiye, located on the border between continents, but many of their venoms remain still understudied. Here, we present the venom compositions of seven Turkish viper taxa. By complementary mass spectrometry-based bottom-up and top-down workflows, the venom profiles were investigated on proteomics and peptidomics level. This study includes the first venom descriptions of Vipera berus barani, Vipera darevskii, Montivipera bulgardaghica albizona, and Montivipera xanthina, as well as the first snake venomics profiles of Turkish Macrovipera lebetinus obtusa, and Daboia palaestinae, including an in-depth reanalysis of M. bulgardaghica bulgardaghica venom. Additionally, we identified the modular consensus sequence pEXW(PZ)1-2P(EI)/(KV)PPLE for bradykinin-potentiating peptides in viper venoms. For better insights into variations and potential impacts of medical significance, the venoms were compared against other Palearctic viper proteomes, including the first genus-wide Montivipera venom comparison. This will help the risk assessment of snakebite envenoming by these vipers and aid in predicting the venoms' pathophysiology and clinical treatments.