Web of venom: exploration of big data resources in animal toxin research.

Research on animal venoms and their components spans multiple disciplines, including biology, biochemistry, bioinformatics, pharmacology, medicine, and more. Manipulating and analyzing the diverse array of data required for venom research can be challenging, and relevant tools and resources are often dispersed across different online platforms, making them less accessible to nonexperts. In this article, we address the multifaceted needs of the scientific community involved in venom and toxin-related research by identifying and discussing web resources, databases, and tools commonly used in this field. We have compiled these resources into a comprehensive table available on the VenomZone website (https://venomzone.expasy.org/10897). Furthermore, we highlight the challenges currently faced by researchers in accessing and using these resources and emphasize the importance of community-driven interdisciplinary approaches. We conclude by underscoring the significance of enhancing standards, promoting interoperability, and encouraging data and method sharing within the venom research community.

Exploring Top-Down Mass Spectrometric Approaches To Probe Forest Cobra (Naja melanoleuca) Venom Proteoforms.

Snake venoms are comprised of bioactive proteins and peptides that facilitate severe snakebite envenomation symptoms. A comprehensive understanding of venom compositions and the subtle heterogeneity therein is important. While bottom-up proteomics has been the well-established approach to catalogue venom compositions, top-down proteomics has emerged as a complementary strategy to characterize venom heterogeneity at the intact protein level. However, top-down proteomics has not been as widely implemented in the snake venom field as bottom-up proteomics, with various emerging top-down methods yet to be developed for venom systems. Here, we have explored three main top-down mass spectrometry methodologies in a proof-of-concept study to characterize selected three-finger toxin and phospholipase A2 proteoforms from the forest cobra (Naja melanoleuca) venom. We demonstrated the utility of a data-independent acquisition mode "MSE" for untargeted fragmentation on a chromatographic time scale and its improvement in protein sequence coverage compared to conventional targeted tandem mass spectrometry analysis. We also showed that protein identification can be further improved using a hybrid fragmentation approach, combining electron-capture dissociation and collision-induced dissociation. Lastly, we reported the promising application of multifunctional cyclic ion mobility separation and post-ion mobility fragmentation on snake venom proteins for the first time.

Where the 'ruber' meets the road: Using the genome of the Red Diamond Rattlesnake to unravel the evolutionary processes driving venom evolution.

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

Melittin and phospholipase A2: Promising anti-cancer candidates from bee venom.

As the research on cancer-related treatment deepens, integrating traditional therapies with emerging interventions reveals new therapeutic possibilities. Melittin and phospholipase A2, the primary anti-cancer components of bee venom, are currently gaining increasing attention. This article reviews the various formulations of melittin in cancer therapy and its potential applications in clinical treatments. The reviewed formulations include melittin analogs, hydrogels, adenoviruses, fusion toxins, fusion peptides/proteins, conjugates, liposomes, and nanoparticles. The article also explored the collaborative therapeutic effects of melittin with natural products, synthetic drugs, radiotherapy, and gene expression regulatory strategies. Phospholipase A2 plays a key role in bee venom anti-cancer strategy due to its unique biological activity. Using an extensive literature review and the latest scientific results, this paper explores the current state and challenges of this field, with the aim to provide new perspectives that guide future research and potential clinical applications. This will further promote the application of bee venom in cancer therapy.

Obscure properties of a traditional herb Pittosporum neelgherrense used to treat snakebite envenoming against Daboia russelli venoms.

In tropical nations, snakebite envenomation is a significant public health issue with negative human and social effects. This is due to three factors: 1) more species of the most hazardous snakes are present; 2) emergency medical assistance is not readily available; and 3) inadequate health care. The problems caused by snakebite envenomation have been partially resolved by immunotherapy. An extensive collection of medicinal herbs is recognized to have antivenomous properties in traditional medicine. However, very few species have undergone scientific investigation, and even fewer have had their active components separated and structurally and functionally defined. In this work, the anti-venom potential of hot and cold aqueous extracts from Pittosporum neelgherrense is evaluated using an in-vitro model. The experimental results showed that 4H-pyran-4-one, 2,3-dihydro-3,5-dihydroxy-6-methyl-(11.20), 1-Undecanol (16.38), Lauryl acetate (18.25), and Cyclotridecane (19.14) were phytochemical substances whose chemical structures were recognized by GCMS. The Direct and Indirect hemorrhagic activity was found to be completely neutralized by P. neelgherrense extract (44.61% hot plant extract & 55.38% cold plant extract) and the zone (2.4 mm), respectively. The neutralization of venoms was indicated by the zone (0.5-0.9 cm) of hydrolysis production of proteolytic activity. Additionally, the results of the gelatine liquefaction study demonstrated that clot formation was not triggered by venom at low concentrations (50:50) but was instead brought on by higher concentrations. The present study suggested that the neutralization of venom by hot water extracts of P. neelgherrense is a potentially therapeutic application.

Insights into the Evolutionary Dynamics: Characterization of Disintegrin and Metalloproteinase Proteins in the Venom Gland Transcriptome of the Hemiscorpius lepturus Scorpion.

BACKGROUND: The Disintegrin and Metalloproteinase (ADAM) family, also known as the metalloproteinase/disintegrin/cysteine-rich (MDC) proteins, includes both secreted and transmembrane molecules involved in critical biological processes, such as cell migration, adhesion, and signaling. This study aimed to investigate the evolutionary relationships and structural characteristics of disintegrin and metalloproteinase proteins identified in the venom gland transcriptome of the scorpion Hemiscorpius lepturus. METHODS: Using bioinformatics tools, we analyzed the open reading frame, conserved motifs, and primary, secondary, and tertiary structures of these proteins. Five proteins, named HLDis- Met1, HLDisMet2, HLDisMet3, HLDisMet4, and HLDisMet5, were identified. Their predicted 3- D structures were within normal ranges (Z-score between -4 to -9). RESULTS: Phylogenetic analysis revealed that HLDisMet1 shares similarities with proteins from various spider species (Nephila pilipes, Argiope bruennichi, Araneus ventricosus, and Trichonephila inaurata madagascariensis), HLDisMet2 with the scorpion Centruroides sculpturatus, HLDis- Met4 with the scorpion Tityus serrulatus, and HLDisMet5 with several snake species (Python bivittatus, Vipera anatolica senliki, Protobothrops mucrosquamatus, and Naja naja). CONCLUSION: These findings highlight the significant similarities between HLDisMet proteins and those found in other venomous species, suggesting a complex and diverse evolutionary pathway for venom components. The cross-species conservation observed may indicate a convergent evolutionary strategy, where different species independently develop similar venom components to adapt to similar ecological niches or prey types. This study highlights the evolutionary significance of venom diversification and its potential applications in understanding venom biology across different species.

Health-related quality of life in severe hypersensitivity reactions: focus on severe allergic asthma and hymenoptera venom anaphylaxis-a cross-sectional study.

Background: Growing evidence reveals the important role of clinical psychological factors in chronic-immune diseases. The aim of this study was to investigate Health-Related Quality of Life (HR-QoL), depression, anxiety, and alexithymia in patients with severe hypersensitivity reactions such as Severe Allergic Asthma (SAA) and Hymenoptera Venom Anaphylaxis (HVA). Methods: The Short-Form Health Survey-36 (SF-36), the Beck Depression Inventory Questionnaire (BDI-II), the Hamilton Anxiety Rating Scale (HAM-A) and the Toronto Alexithymia Scale (TAS-20) were used to assess HR-QoL and clinical psychological features of patients with SAA and HVA. Results: Overall, 78 patients were recruited. Patients with SAA (n = 35) reported lower scores for physical functioning [65 (58-75) vs. 90 (85-95); p = <0.001], role limitations due to physical health [25 (0-50) vs. 62 (50-75); p = 0.004], bodily pain [47.5 (41.1-61.3) vs. 55.5 (55-96); p = 0.001], general health [40 (30-60) vs. 70 (50-80); p = 0.0003] and social functioning [50 (37.5-62.5) vs. 62.5 (54.9-75); p = 0.007] while higher scores for depressive symptoms [14 (11-15.4) vs. (9.5 (6-15.4); p = 0.05)] compared to HVA patients (n = 43). All the dimensions of SF-36 were negatively correlated with anxiety (r from -0.26 to -0.66; p all < 0.01) and depressive symptoms (r from -0.44 to -0.73; p all < 0.001). Alexithymia was negatively correlated with vitality (r = -0.28; p = 0.02) and mental health (r = -027; p = 0.03). Additionally, patients with alexithymia (38% of participants) showed higher levels of depressive symptoms [9.5 (10-19) vs. 14 (6-13.9); p = 0.005] and anxiety levels [31 (27.9-35) vs. 24 (16-33.9); p = 0.02]; they also showed less vitality [40 (39.9-50) vs. 55 (50-60) p = 0.01], social functioning [50 (37.5-62.5) vs. 62.5 (50 vs. 75); p = 0.01] and mental health [48 (44-60) vs. 68 (56-76); p = 0.004]. Conclusion: Clinical psychological features due to severe hypersensitive reactions may contribute to the patient's perceived HR-QoL. Focused clinical psychological interventions should be promoted to improve the clinical management of such conditions.

Optimizing Nav1.7-Targeted Analgesics: Revealing Off-Target Effects of Spider Venom-Derived Peptide Toxins and Engineering Strategies for Improvement.

The inhibition of Nav1.7 is a promising strategy for the development of analgesic treatments. Spider venom-derived peptide toxins are recognized as significant sources of Nav1.7 inhibitors. However, their development has been impeded by limited selectivity. In this study, eight peptide toxins from three distinct spider venom Nav channel families demonstrated robust inhibition of hNav1.7, rKv4.2, and rKv4.3 (rKv4.2/4.3) currents, exhibiting a similar mode of action. The analysis of structure and function relationship revealed a significant overlap in the pharmacophore responsible for inhibiting hNav1.7 and rKv4.2 by HNTX-III, although Lys25 seems to play a more pivotal role in the inhibition of rKv4.2/4.3. Pharmacophore-guided rational design is employed for the development of an mGpTx1 analogue, mGpTx1-SA, which retains its inhibition of hNav1.7 while significantly reducing its inhibition of rKv4.2/4.3 and eliminating cardiotoxicity. Moreover, mGpTx1-SA demonstrates potent analgesic effects in both inflammatory and neuropathic pain models, accompanied by an improved in vivo safety profile. The results suggest that off-target inhibition of rKv4.2/4.3 by specific spider peptide toxins targeting hNav1.7 may arise from a conserved binding motif. This insight promises to facilitate the design of hNav1.7-specific analgesics, aimed at minimizing rKv4.2/4.3 inhibition and associated toxicity, thereby enhancing their suitability for therapeutic applications.

Potential mitochondrial ROS-mediated damage induced by chitosan nanoparticles bee venom-loaded on cancer cell lines.

Recently, numerous studies have confirmed the importance of chitosan nanoparticles (CNP) as a viable drug delivery carrier for increasing the efficacy of anticancer drugs in cancer treatment. It is a macromolecule and natural biopolymer compound, more stable and safer in use than metal nanoparticles. Bee venom (BV), a form of defense venom, has been shown to have anti-tumor, neuroprotective, anti-inflammatory, analgesic, and anti-infectivity properties. Moreover, the regulation of cell death has been linked to reactive oxygen species (ROS)-mediated cell apoptosis, which induces mitochondrial damage and ER stress through oxidative stress events. Therefore, this study aimed to illustrate the ROS-mediated effect on the cancer cells treatment with CNP-loaded BV (CNP-BV) and explained the adverse effects of ROS generation on Mitochondria and ER. We have found that the targeted CNP-BV were high in cytotoxicity against MCF-7 (IC50 437.2 µg/mL) and HepG2 (IC50 109.5 µg/mL) through the induction of massive generation of ROS, which in turn results in activating the mitochondrial cascade and ER stress. These results highlighted the role of ROS generation in inducing apoptosis in cancer cells.

The state of parasitoid wasp genomics.

Parasitoid wasps represent a group of parasitic insects with high species diversity that have played a pivotal role in biological control and evolutionary studies. Over the past 20 years, developments in genomics have greatly enhanced our understanding of the biology of these species. Technological leaps in sequencing have facilitated the improvement of genome quality and quantity, leading to the availability of hundreds of parasitoid wasp genomes. Here, we summarize recent progress in parasitoid wasp genomics, focusing on the evolution of genome size (GS) and the genomic basis of several key traits. We also discuss the contributions of genomics in studying venom evolution and endogenization of viruses. Finally, we advocate for increased sequencing and functional research to better understand parasitoid biology and enhance biological control.

Taipan envenoming south of the border.

We present a case of severe taipan envenoming in northern New South Wales in a 68-year-old man. He developed severe neurotoxicity requiring intubation and ventilation, venom-induced consumption coagulopathy, myotoxicity and thrombotic microangiopathy with acute kidney injury requiring dialysis. He was administered brown and tiger snake antivenom consistent with guidelines and snake occurrence in the region. Taipan venom was detected in serum (72 ng/ml) following concern about the severity of neurotoxicity, clinical toxicology consultation and a concurrent report of a taipan in the area. Based on this it would be prudent to stock and consider treating with polyvalent antivenom in north-eastern New South Wales and south-eastern Queensland.

[Research Progress in Detection of Bee Venom Allergens].

Hymenopteran insect stings are a risk factor that cannot be ignored for the people allergic to hymenopteran venoms.In China,the current diagnostic tools cannot provide accurate information to identify sensitized insects,thus affecting clinical diagnosis and treatment.Honeybee is a common hymenopteran insect.Due to its wide distribution,large number,and complex venom composition,researchers have carried out recombination schemes for the main allergens of honeybee venom,laying a theoretical foundation for the detection of allergens.The development of diagnostic technologies for allergen components can accurately detect bee venom allergens,providing a new set of clinical diagnosis and treatment schemes for the population allergic to bee venom.

Engineering of pH-dependent antigen binding properties for toxin-targeting IgG1 antibodies using light-chain shuffling.

Immunoglobulin G (IgG) antibodies that bind their cognate antigen in a pH-dependent manner (acid-switched antibodies) can release their bound antigen for degradation in the acidic environment of endosomes, while the IgGs are rescued by the neonatal Fc receptor (FcRn). Thus, such IgGs can neutralize multiple antigens over time and therefore be used at lower doses than their non-pH-responsive counterparts. Here, we show that light-chain shuffling combined with phage display technology can be used to discover IgG1 antibodies with increased pH-dependent antigen binding properties, using the snake venom toxins, myotoxin II and α-cobratoxin, as examples. We reveal differences in how the selected IgG1s engage their antigens and human FcRn and show how these differences translate into distinct cellular handling properties related to their pH-dependent antigen binding phenotypes and Fc-engineering for improved FcRn binding. Our study showcases the complexity of engineering pH-dependent antigen binding IgG1s and demonstrates the effects on cellular antibody-antigen recycling.

Novel GPIb-independent platelet aggregation induced by botrocetin: implications for diagnosis and antithrombotic therapy.

BACKGROUND: Snake venom botrocetin facilitates von Willebrand factor (VWF) binding to platelet GPIbα and has been widely used for the diagnosis of von Willebrand disease and GPIb-related disorders. Botrocetin is also commonly employed for the development/characterization of antithrombotics targeting the GPIb-VWF axis. OBJECTIVES: To explore the alternative receptor(s)/mechanisms that participate in botrocetin-induced platelet aggregation. METHODS: The effects of botrocetin on platelet aggregation were examined using platelets from wild-type, VWF- and fibrinogen-deficient, GPIbα-deficient, IL4Rα/GPIbα-transgenic, ITGA2B and ITGB3-deficient mice, and Bernard-Soulier syndrome and healthy human samples. Platelet-fibrinogen and platelet-VWF interaction were measured using flow cytometry. GPIbα-VWF binding was evaluated utilizing enzyme-linked immunosorbent assay. Botrocetin-αIIbß3 and botrocetin-GPIbα interactions were measured using enzyme-linked immunosorbent assay and fluorescence anisotropy assays. Heparinized whole blood from healthy donors was examined for thrombus formation and growth in a perfusion chamber. RESULTS: Botrocetin could induce aggregation of platelets from a Bernard-Soulier syndrome patient and GPIbα-deficient mice as well as platelets lacking the N-terminal extracellular domain of GPIbα. Botrocetin could interact with αIIbß3 and facilitated αIIbß3-VWF interaction independent of GPIb. Botrocetin competitively bound to the ligand-binding domain of activated rather than resting αIIbß3. Although botrocetin-induced platelet aggregation requires VWF, strikingly, in the absence of VWF, botrocetin blocked fibrinogen and other ligand binding to αIIbß3 and inhibited platelet aggregation and thrombus formation. Consistently, recombinant botrocetin defective in VWF binding inhibited αIIbß3- and GPIb-mediated platelet aggregation, spreading, and thrombus formation. CONCLUSION: Our study provides insights into avoiding the misdiagnosis of GPIb-related disorders and developing botrocetin mutants as potential new antithrombotics that may simultaneously target both αIIbß3 and GPIbα.

Modulation of response to braconid wasp venom by adipokinetic hormone in Drosophila melanogaster.

The minute wasp Habrobracon hebetor venom (HH venom) is a potent cocktail of toxins that paralyzes the victim's muscles and suppresses humoral and cellular immunity. This study examined the effect of HH venom on specific biochemical, physiological, and ultrastructural characteristics of the thoracic and nervous (CNS) tissues of Drosophila melanogaster under in vitro conditions. Venom treatment modulated the activities of superoxide dismutase (SOD) and catalase (CAT), endogenous Drome-AKH level, and affected the relative viability of the cells. Additionally, it reduced the expression of genes related to the immune system in the CNS, including Keap1, Relish, Nox, Eiger, Gadd45, and Domeless, as well as in the thoracic muscles, except for Nox. Besides, venom treatment led to deteriorative changes in the ultrastructure of muscle cells, particularly affecting the mitochondria. When venom and Drosophila melanogaster-adipokinetic hormone (Drome-AKH) were applied together, the effects of the venom alone were often modulated. The harmful effect of the venom on SOD activity was relatively reduced and the activity returned to a level similar to that of the control. In the CNS, the simultaneous application of venom and hormones abolished the suppression of previously reported immune-related genes (except for Gadd45), whereas in the muscles, this was only true for Eiger. Additionally, Drome-AKH restored cell structure to a level comparable to that of the control and lessened the harmful effects of HH venom on muscle mitochondria. These findings suggest a general body response of D. melanogaster to HH venom and a partial defensive role of Drome-AKH in this process.

The geographical distribution of scorpions, implication of venom toxins, envenomation, and potential therapeutics in Southern and Northern Africa.

Scorpions are predatory arachnids whose venomous sting primarily affects people in tropical and subtropical regions. Most scorpion stings can only cause localized pain without severe envenomation. Less than one-third of the stings cause systemic envenoming and possibly lead to death. About 350,000 scorpion stings in Northern Africa are recorded yearly, resulting in about 810 deaths. In Eastern/Southern Africa, there are about 79,000 stings recorded yearly, resulting in 245 deaths. Farmers and those living in poverty-stricken areas are among the most vulnerable to getting stung by scorpions. However, compared to adults, children are at greater risk of severe envenomation. Scorpion venom is made up of complex mixtures dominated by peptides and proteins that confer its potency and toxicity. These venom toxins have intra- and interspecies variations associated with the scorpion's habitat, sex, diet, and age. These variations alter the activity of antivenoms used to treat scorpion sting envenomation. Thus, the study of the proteome composition of medically important scorpion venoms needs to be scaled up along their geographical distribution and contributions to envenomation in Southern and Northern Africa. This will help the production of safer, more effective, and broad-spectrum antivenoms within these regions. Here, we review the clinical implications of scorpion sting envenomation in Southern and Northern Africa. We further highlight the compositions of scorpion venoms and tools used in scorpion venomics. We discuss current antivenoms used against scorpion sting envenomation and suggestions for future production of better antivenoms or alternatives. Finally, we discuss the therapeutic properties of scorpion venom.

Reported snakebite mortality and state compensation payments in Madhya Pradesh, India, from 2020 to 2022.

BACKGROUND: India experiences the highest snakebite burden globally, with 58 000 predicted deaths annually. The central Indian state of Madhya Pradesh is thought to have a substantial snakebite burden and provides compensation to families who can demonstrate by postmortem and hospital treatment reports that their relatives have died due to snakebite. This study represents the first report on the frequency of distribution of compensation for snakebite deaths in Madhya Pradesh. METHODS: Statewide snakebite death compensation data from 2020-2021 and 2021-2022, provided by the Madhya Pradesh health authorities, were analysed alongside interviews with 15 families that described the events that ultimately led to their compensation claims. RESULTS: Compensation was paid to a total of 5728 families, with a total value equating to 22 912 Lakhs (approximately US${\$}$27.94 million). Families described commonly recognised snakebite risk factors and behaviours in the events that resulted in their relatives' deaths. CONCLUSIONS: The snakebite burden in Madhya Pradesh is significant, both in terms of mortality and economic expenditure of the state. Sustained investment in preventative interventions, as well as monitoring of the rate of compensation payouts due to snakebite death as a measure of intervention effectiveness, should be considered to substantially reduce snakebite incidence and mortality.

Dual functionality of pathogenesis-related proteins: defensive role in plants versus immunosuppressive role in pathogens.

Plants respond to pathogen exposure by activating the expression of a group of defense-related proteins known as Pathogenesis-Related (PR) proteins, initially discovered in the 1970s. These PR proteins are categorized into 17 distinct families, denoted as PR1-PR17. Predominantly secreted, most of these proteins execute their defensive roles within the apoplastic space. Several PR proteins possess well-defined enzymatic functions, such as ß-glucanase (PR2), chitinases (PR3, 4, 8, 11), proteinase (PR7), or RNase (PR10). Enhanced resistance against pathogens is observed upon PR protein overexpression, while their downregulation renders plants more susceptible to pathogen infections. Many of these proteins exhibit antimicrobial activity in vitro, and due to their compact size, some are classified as antimicrobial peptides. Recent research has unveiled that phytopathogens, including nematodes, fungi, and phytophthora, employ analogous proteins to bolster their virulence and suppress plant immunity. This raises a fundamental question: how can these conserved proteins act as antimicrobial agents when produced by the host plant but simultaneously suppress plant immunity when generated by the pathogen? In this hypothesis, we investigate PR proteins produced by pathogens, which we term "PR-like proteins," and explore potential mechanisms by which this class of virulence factors operate. Preliminary data suggests that these proteins may form complexes with the host's own PR proteins, thereby interfering with their defense-related functions. This analysis sheds light on the intriguing interplay between plant and pathogen-derived PR-like proteins, providing fresh insights into the intricate mechanisms governing plant-pathogen interactions.

Spider and scorpion knottins targeting voltage-gated sodium ion channels in pain signaling.

In sensory neurons that transmit pain signals, whether acute or chronic, voltage-gated sodium channels (VGSCs) are crucial for regulating excitability. NaV1.1, NaV1.3, NaV1.6, NaV1.7, NaV1.8, and NaV1.9 have been demonstrated and defined their functional roles in pain signaling based on their biophysical properties and distinct patterns of expression in each subtype of sensory neurons. Scorpions and spiders are traditional Chinese medicinal materials, belonging to the arachnid class. Most of the studied species of them have evolved venom peptides that exhibit a wide variety of knottins specifically targeting VGSCs with subtype selectivity and conformational specificity. This review provides an overview on the exquisite knottins from scorpion and spider venoms targeting pain-related NaV channels, describing the sequences and the structural features as well as molecular determinants that influence their selectivity on special subtype and at particular conformation, with an aim for the development of novel research tools on NaV channels and analgesics with minimal adverse effects.

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.