Kinetic analysis of effects of temperature and time on the regulation of venom expression in Bungarus multicinctus

Intraspecific venom variation in the medically important puff adder (Bitis arietans): Comparative venom gland transcriptomics, in vitro venom activity and immunological recognition by antivenom

BACKGROUND

Variation in snake venoms is well documented, both between and within species, with intraspecific venom variation often correlated with geographically distinct populations. The puff adder, Bitis arietans, is widely distributed across sub-Saharan Africa and into the Arabian Peninsula where it is considered a leading cause of the ~310,000 annual snakebites across the region, with its venom capable of causing substantial morbidity and mortality. Despite its medical importance and wide geographic distribution, there is little known about venom variation between different B. arietans populations and the potential implications of this variation on antivenom efficacy.

METHODOLOGY

We applied a range of analyses, including venom gland transcriptomics, in vitro enzymatic assays and reverse phase chromatography to comparatively analyse B. arietans venoms originating from Nigeria, Tanzania, and South Africa. Immunological assays and in vitro enzymatic neutralisation assays were then applied to investigate the impact of venom variation on the potential efficacy of three antivenom products; SAIMR Polyvalent, EchiTAb-Plus and Fav-Afrique.

FINDINGS

Through the first comparison of venom gland transcriptomes of B. arietans from three geographically distinct regions (Nigeria, Tanzania, and South Africa), we identified substantial variation in toxin expression. Findings of venom variation were further supported by chromatographic venom profiling, and the application of enzymatic assays to quantify the activity of three pathologically relevant toxin families. However, the use of western blotting, ELISA, and in vitro enzymatic inhibition assays revealed that variation within B. arietans venom does not appear to substantially impact upon the efficacy of three African polyvalent antivenoms.

CONCLUSIONS

The large distribution and medical importance of B. arietans makes this species ideal for understanding venom variation and the impact this has on therapeutic efficacy. The findings in this study highlight the likelihood for considerable venom toxin variation across the range of B. arietans, but that this may not dramatically impact upon the utility of treatment available in the region.

Polyploidization of Indotyphlops braminus: evidence from isoform-sequencing

BACKGROUND

Indotyphlops braminus, the only known triploid parthenogenetic snake, is a compelling species for revealing the mechanism of polyploid emergence in vertebrates.

METHODS

In this study, we applied PacBio isoform sequencing technology to generate the first full-length transcriptome of I. braminus, aiming to improve the understanding of the molecular characteristics of this species.

RESULTS

A total of 51,849 nonredundant full-length transcript assemblies (with an N50 length of 2980 bp) from I. braminus were generated and fully annotated using various gene function databases. Our analysis provides preliminary evidence supporting a recent genome duplication event in I. braminus. Phylogenetic analysis indicated that the divergence of I. braminus subgenomes occurred approximately 11.5 ~ 15 million years ago (Mya). The full-length transcript resource generated as part of this research will facilitate transcriptome analysis and genomic evolution studies in the future.

Distinct regulatory networks control toxin gene expression in elapid and viperid snakes

BACKGROUND

Venom systems are ideal models to study genetic regulatory mechanisms that underpin evolutionary novelty. Snake venom glands are thought to share a common origin, but there are major distinctions between venom toxins from the medically significant snake families Elapidae and Viperidae, and toxin gene regulatory investigations in elapid snakes have been limited. Here, we used high-throughput RNA-sequencing to profile gene expression and microRNAs between active (milked) and resting (unmilked) venom glands in an elapid (Eastern Brown Snake, Pseudonaja textilis), in addition to comparative genomics, to identify cis- and trans-acting regulation of venom production in an elapid in comparison to viperids (Crotalus viridis and C. tigris).

RESULTS

Although there is conservation in high-level mechanistic pathways regulating venom production (unfolded protein response, Notch signaling and cholesterol homeostasis), there are differences in the regulation of histone methylation enzymes, transcription factors, and microRNAs in venom glands from these two snake families. Histone methyltransferases and transcription factor (TF) specificity protein 1 (Sp1) were highly upregulated in the milked elapid venom gland in comparison to the viperids, whereas nuclear factor I (NFI) TFs were upregulated after viperid venom milking. Sp1 and NFI cis-regulatory elements were common to toxin gene promoter regions, but many unique elements were also present between elapid and viperid toxins. The presence of Sp1 binding sites across multiple elapid toxin gene promoter regions that have been experimentally determined to regulate expression, in addition to upregulation of Sp1 after venom milking, suggests this transcription factor is involved in elapid toxin expression. microRNA profiles were distinctive between milked and unmilked venom glands for both snake families, and microRNAs were predicted to target a diversity of toxin transcripts in the elapid P. textilis venom gland, but only snake venom metalloproteinase transcripts in the viperid C. viridis venom gland. These results suggest differences in toxin gene posttranscriptional regulation between the elapid P. textilis and viperid C. viridis.

CONCLUSIONS

Our comparative transcriptomic and genomic analyses between toxin genes and isoforms in elapid and viperid snakes suggests independent toxin regulation between these two snake families, demonstrating multiple different regulatory mechanisms underpin a venomous phenotype.

Diversity and Evolutionary Analysis of Venom Insulin Derived from Cone Snails

Cone snails possess a diverse array of novel peptide toxins, which selectively target ion channels and receptors in the nervous and cardiovascular systems. These numerous novel peptide toxins are a valuable resource for future marine drug development. In this review, we compared and analyzed the sequence diversity, three-dimensional structural variations, and evolutionary aspects of venom insulin derived from different cone snail species. The comparative analysis reveals that there are significant variations in the sequences and three-dimensional structures of venom insulins from cone snails with different feeding habits. Notably, the venom insulin of some piscivorous cone snails exhibits a greater similarity to humans and zebrafish insulins. It is important to emphasize that these venom insulins play a crucial role in the predatory strategies of these cone snails. Furthermore, a phylogenetic tree was constructed to trace the lineage of venom insulin sequences, shedding light on the evolutionary interconnections among cone snails with diverse diets.

Comparing morphological and secretory aspects of cephalic glands among the New World coral snakes brings novel insights on their biological roles

Oral and other cephalic glands have been surveyed by several studies with distinct purposes. Despite the wide diversity and medical relevance of the New World coral snakes, studies focusing on understanding the biological roles of the glands within this group are still scarce. Specifically, the venom glands of some coral snakes were previously investigated but all other cephalic glands remain uncharacterized. In this sense, performing morphological and molecular analysis of these glands may help better understand their biological role. Here, we studied the morphology of the venom, infralabial, rictal, and harderian glands of thirteen species of Micrurus and Micruroides euryxanthus. We also performed a molecular characterization of these glands from selected species of Micrurus using transcriptomic and proteomic approaches. We described substantial morphological variation in the cephalic glands of New World coral snakes and structural evidence for protein-secreting cells in the inferior rictal glands. Our molecular analysis revealed that the venom glands, as expected, are majorly devoted to toxin production, however, the infralabial and inferior rictal glands also expressed some toxin genes at low to medium levels, despite the marked morphological differences. On the other hand, the harderian glands were dominated by the expression of lipocalins, but do not produce toxins. Our integrative analysis, including the prediction of biological processes and pathways, helped decipher some important traits of cephalic glands and better understand their biology.

The structural and functional divergence of a neglected three-finger toxin subfamily in lethal elapids

Bungarus multicinctus is a widely distributed and medically important elapid snake that produces lethal neurotoxic venom. To study and enhance existing antivenom, we explore the complete repertoire of its toxin genes based on de novo chromosome-level assembly and multi-tissue transcriptome data. Comparative genomic analyses suggest that the three-finger toxin family (3FTX) may evolve through the neofunctionalization of flanking LY6E. A long-neglected 3FTX subfamily (i.e., MKA-3FTX) is also investigated. Only one MKA-3FTX gene, which evolves a different protein conformation, is under positive selection and actively transcribed in the venom gland, functioning as a major toxin effector together with MKT-3FTX subfamily homologs. Furthermore, this lethal snake may acquire self-resistance to its β-bungarotoxin via amino acid replacements on fast-evolving KCNA2. This study provides valuable resources for further evolutionary and structure-function studies of snake toxins, which are fundamental for the development of effective antivenoms and drug candidates.

The Enzymatic Core of Scorpion Venoms

Enzymes are an integral part of animal venoms. Unlike snakes, in which enzymes play a primary role in envenomation, in scorpions, their function appears to be ancillary in most species. Due to this, studies on the diversity of scorpion venom components have focused primarily on the peptides responsible for envenomation (toxins) and a few others (e.g., antimicrobials), while enzymes have been overlooked. In this work, a comprehensive study on enzyme diversity in scorpion venoms was performed by transcriptomic and proteomic techniques. Enzymes of 63 different EC types were found, belonging to 330 orthogroups. Of them, 24 ECs conform the scorpion venom enzymatic core, since they were determined to be present in all the studied scorpion species. Transferases and lyases are reported for the first time. Novel enzymes, which can play different roles in the venom, including direct toxicity, as venom spreading factors, activators of venom components, venom preservatives, or in prey pre-digestion, were described and annotated. The expression profile for transcripts coding for venom enzymes was analyzed, and shown to be similar among the studied species, while being significantly different from their expression pattern outside the telson.

A Combined Transcriptomics and Proteomics Approach Reveals the Differences in the Predatory and Defensive Venoms of the Molluscivorous Cone Snail Cylinder ammiralis (Caenogastropoda: Conidae)

Venoms are complex mixtures of proteins that have evolved repeatedly in the animal kingdom. Cone snail venoms represent one of the best studied venom systems. In nature, this venom can be dynamically adjusted depending on its final purpose, whether to deter predators or hunt prey. Here, the transcriptome of the venom gland and the proteomes of the predation-evoked and defensive venoms of the molluscivorous cone snail Cylinder ammiralis were catalogued. A total of 242 venom-related transcripts were annotated. The conotoxin superfamilies presenting more different peptides were O1, O2, T, and M, which also showed high expression levels (except T). The three precursors of the J superfamily were also highly expressed. The predation-evoked and defensive venoms showed a markedly distinct profile. A total of 217 different peptides were identified, with half of them being unique to one venom. A total of 59 peptides ascribed to 23 different protein families were found to be exclusive to the predatory venom, including the cono-insulin, which was, for the first time, identified in an injected venom. A total of 43 peptides from 20 protein families were exclusive to the defensive venom. Finally, comparisons of the relative abundance (in terms of number of peptides) of the different conotoxin precursor superfamilies showed that most of them present similar abundance regardless of the diet.

Serum metabolomics of Bama miniature pigs bitten by Bungarus multicinctus

Bungarus multicinctus is one of the top ten venomous snakes in China, and its bite causes acute and severe diseases, but its pathophysiology remains poorly elucidated. Thus, an animal model of Bungarus multicinctus bite was established by intramuscular injection of 30μg/kg of Bungarus multicinctus venom, and then the serum metabolites were subsequently screened, identified and validated by ultra-performance liquid chromatography-quadrupole-time of flight-mass spectrometry (UPLC-Q-TOF-MS) methods to explore the potential biomakers and possible metabolic pathways. Untargeted metabolomics analysis showed that 36 and 38 endogenous metabolites levels changed in ESI+ and ESI-, respectively, KEGG pathway analysis showed that 5 metabolic pathways, including mineral absorption, central carbon metabolism in cancer, protein digestion and absorption, aminoacyl-tRNA biosynthesis and ABC transporters might be closely related to Bungarus multicinctus bite. Targeted metabolomics analysis showed that there were significant differences in serum D-proline, L-leucine and L-glutamine after Bungarus multicinctus bite (P < 0.05). In addition, receiver operating characteristic (ROC) analysis showed that the diagnostic efficiency of L-Glutamine was superior to other potential biomarkers and the AUC value was 0.944. Moreover, we found evidence for differences in the pathophysiology of glutamine between Bungarus multicinctus bite group and normal group, specifically with the content of glutamine synthetase (GS) and glutaminase (GLS). Taken together, the current study has successfully established an animal model of Bungarus multicinctus bite, and further identified the links between the metabolic perturbations and the pathophysiology and the potential diagnostic biomakers of Bungarus multicinctus bite, which provided valuable insights for studying the mechanism of Bungarus multicinctus bite.

Impact of Temperature on the Immune Interaction between a Parasitoid Wasp and Drosophila Host Species

Temperature is particularly important for ectotherms, including endoparasitoid wasps that develop inside another ectotherm host. In this study, we tested the impact of three temperatures (20 °C, 25 °C and 30 °C) on the host-parasitoid immune interaction using two Drosophila host species (Drosophila melanogaster and D. yakuba) and two parasitoid lines of Leptopilina boulardi. Drosophila's immune defense against parasitoids consists of the formation of a melanized capsule surrounding the parasitoid egg. To counteract this response, Leptopilina parasitoids rely on the injection of venom during oviposition. Here, we tested the effect of temperature on parasitic success and host encapsulation capacity in response to a parasitoid egg or other foreign body. Increased temperature either promoted or did not affect the parasitic success, depending on the parasitoid-host pairs considered. The mechanisms behind the higher success seemed to vary depending on whether the temperature primarily affected the host immune response or also affected the parasitoid counter-immune response. Next, we tested the effect of parasitoid rearing temperature on its success and venom composition. Venom composition varied strongly with temperature for both parasitoid lines, partially consistent with a change in their parasitic success. Overall, temperature may have a significant impact on the host-parasitoid immune interaction.

The influence of ecological factors on cnidarian venoms

Venom research is often focussed on medical relevance, novel compounds and venom evolution, whilst studying the relationship between a venom and its environment – venom ecology - has been conducted to a lesser extent. Given the projected environmental changes envisioned to occur with global warming, it is pertinent now more than ever, to highlight this topic. Here we review literature examining the influence of ecological factors such as environmental temperature, salinity, ontogeny, geographic location and diet on cnidarian venoms. This review provides an exclusive focus on the cnidarian phylum and encompasses all available published, peer-reviewed literature to our knowledge regarding the ecological factors influencing venom. We find a startling lack of research into the effects of both environmental and biological factors on venoms, with very few to no studies available per category. Importantly, research does exist that suggest these ecological processes may influence other marine or terrestrial venoms, thus we recommend future research is needed to explore this concept in cnidarians.

•

Cnidarian toxins are significantly affected by environment and biology, yet literature on the subject is scarce.

•

Temperature, ontogeny, geographic location and diet can influence cnidarian venoms.

•

Salinity can influence other marine toxins, but literature is lacking for cnidarians.

•

More venom ecology research is needed in medically important species, if medical treatments are to advance.

A Wolf in Another Wolf's Clothing: Post-Genomic Regulation Dictates Venom Profiles of Medically-Important Cryptic Kraits in India

The Common Krait (Bungarus caeruleus) shares a distribution range with many other ‘phenotypically-similar’ kraits across the Indian subcontinent. Despite several reports of fatal envenomings by other Bungarus species, commercial Indian antivenoms are only manufactured against B. caeruleus. It is, therefore, imperative to understand the distribution of genetically distinct lineages of kraits, the compositional differences in their venoms, and the consequent impact of venom variation on the (pre)clinical effectiveness of antivenom therapy. To address this knowledge gap, we conducted phylogenetic and comparative venomics investigations of kraits in Southern and Western India. Phylogenetic reconstructions using mitochondrial markers revealed a new species of krait, Romulus’ krait (Bungarus romulusi sp. nov.), in Southern India. Additionally, we found that kraits with 17 mid-body dorsal scale rows in Western India do not represent a subspecies of the Sind Krait (B. sindanus walli) as previously believed, but are genetically very similar to B. sindanus in Pakistan. Furthermore, venom proteomics and comparative transcriptomics revealed completely contrasting venom profiles. While the venom gland transcriptomes of all three species were highly similar, venom proteomes and toxicity profiles differed significantly, suggesting the prominent role of post-genomic regulatory mechanisms in shaping the venoms of these cryptic kraits. In vitro venom recognition and in vivo neutralisation experiments revealed a strong negative impact of venom variability on the preclinical performance of commercial antivenoms. While the venom of B. caeruleus was neutralised as per the manufacturer’s claim, performance against the venoms of B. sindanus and B. romulusi was poor, highlighting the need for regionally-effective antivenoms in India.