The UniProtKB/Swiss-Prot Tox-Prot program: A central hub of integrated venom protein data

Redistribution of Ancestral Functions Underlies the Evolution of Venom Production in Marine Predatory Snails

Venom-secreting glands are highly specialized organs evolved throughout the animal kingdom to synthetize and secrete toxins for predation and defense. Venom is extensively studied for its toxin components and application potential; yet, how animals become venomous remains poorly understood. Venom systems therefore offer a unique opportunity to understand the molecular mechanisms underlying functional innovation. Here, we conducted a multispecies multi-tissue comparative transcriptomics analysis of 12 marine predatory gastropod species, including species with venom glands and species with homologous non-venom-producing glands, to examine how specialized functions evolve through gene expression changes. We found that while the venom gland specialized for the mass production of toxins, its homologous glands retained the ancestral digestive functions. The functional divergence and specialization of the venom gland were achieved through a redistribution of its ancestral digestive functions to other organs, specifically the esophagus. This entailed concerted expression changes and accelerated transcriptome evolution across the entire digestive system. The increase in venom gland secretory capacity was achieved through the modulation of an ancient secretory machinery, particularly genes involved in endoplasmic reticulum stress and unfolded protein response. This study shifts the focus from the well-explored evolution of toxins to the lesser-known evolution of the organ and mechanisms responsible for venom production. As such, it contributes to elucidating the molecular mechanisms underlying organ evolution at a fine evolutionary scale, highlighting the specific events that lead to functional divergence.

Venomous Cargo: Diverse Toxin-Related Proteins Are Associated with Extracellular Vesicles in Parasitoid Wasp Venom

Unusual membrane-bound particles are present in the venom of the parasitoid wasps that parasitize Drosophila melanogaster. These venom particles harbor about 400 proteins and suppress the encapsulation of a wasp egg. Whereas the proteins in the particles of Leptopilina boulardi venom modify host hemocyte properties, those in L. heterotoma kill host hemocytes. The mechanisms underlying this differential effect are not well understood. The proteome of the L. heterotoma venom particles has been described before, but that of L. boulardi has not been similarly examined. Using sequence-based programs, we report the presence of conserved proteins in both proteomes with strong enrichment in the endomembrane and exosomal cell components. Extracellular vesicle markers are present in both proteomes, as are numerous toxins. Both proteomes also contain proteins lacking any annotation. Among these, we identified the proteins with structural similarity to the ADP-ribosyltransferase enzymes involved in bacterial virulence. We propose that invertebrate fluids like parasitoid venom contain functional extracellular vesicles that deliver toxins and virulence factors from a parasite to a host. Furthermore, the presence of such vesicles may not be uncommon in the venom of other animals. An experimental verification of the predicted toxin functions will clarify the cellular mechanisms underlying successful parasitism.

Proteomic Diversity of the Sea Anemone Actinia fragacea: Comparative Analysis of Nematocyst Venom, Mucus, and Tissue-Specific Profiles

Sea anemones (Actiniaria, Cnidaria) are promising targets for biomedical research, as they produce unique bioactive compounds, including toxins and antimicrobial peptides (AMPs). However, the diversity and mechanisms underlying their chemical defenses remain poorly understood. In this study, we investigate the proteomic profiles of the unexplored sea anemone Actinia fragacea by analyzing its venom nematocyst extract, tissues, and mucus secretion. A total of 4011 different proteins were identified, clustered into 3383 protein groups. Among the 83 putative toxins detected, actinoporins, neurotoxins, and phospholipase A2 were uncovered, as well as two novel zinc metalloproteinases with two specific domains (ShK) associated with potassium channel inhibition. Common Gene Ontology (GO) terms were related to immune responses, cell adhesion, protease inhibition, and tissue regeneration. Furthermore, 1406 of the 13,276 distinct peptides identified were predicted as potential AMPs, including a putative Aurelin-like AMP localized within the nematocysts. This discovery highlights and strengthens the evidence for a cnidarian-exclusive Aurelin peptide family. Several other bioactive compounds with distinctive defense functions were also detected, including enzymes, pattern recognition proteins (PRPs), and neuropeptides. This study provides the first proteome map of A. fragacea, offering a critical foundation for exploring novel bioactive compounds and valuable insights into its molecular complexity.

Molecular Display of the Animal Meta-Venome for Discovery of Novel Therapeutic Peptides

Animal venoms, distinguished by their unique structural features and potent bioactivities, represent a vast and relatively untapped reservoir of therapeutic molecules. However, limitations associated with comprehensively constructing and expressing highly complex venom and venom-like molecule libraries have precluded their therapeutic evaluation via high-throughput screening. Here, we developed an innovative computational approach to design a highly diverse library of animal venoms and "metavenoms". We used programmable M13 hyperphage display to preserve critical disulfide-bonded structures for highly parallelized single-round biopanning with quantitation via high-throughput DNA sequencing. Our approach led to the discovery of Kunitz-type domain containing proteins that target the human itch receptor Mas-related G-protein coupled receptor member X4, which plays a crucial role in itch perception. Deep learning-based structural homology mining identified two endogenous human homologs, tissue factor pathway inhibitor (TFPI), and serine peptidase inhibitor, Kunitz type 2 (SPINT2), which exhibit agonist-dependent potentiation of Mas-related G-protein coupled receptor member X4. Highly multiplexed screening of animal venoms and metavenoms is therefore a promising approach to uncover new drug candidates.

Current Technologies in Snake Venom Analysis and Applications

This comprehensive review explores the cutting-edge advancements in snake venom research, focusing on the integration of proteomics, genomics, transcriptomics, and bioinformatics. Highlighting the transformative impact of these technologies, the review delves into the genetic and ecological factors driving venom evolution, the complex molecular composition of venoms, and the regulatory mechanisms underlying toxin production. The application of synthetic biology and multi-omics approaches, collectively known as venomics, has revolutionized the field, providing deeper insights into venom function and its therapeutic potential. Despite significant progress, challenges such as the functional characterization of toxins and the development of cost-effective antivenoms remain. This review also discusses the future directions of venom research, emphasizing the need for interdisciplinary collaborations and new technologies (mRNAs, cryo-electron microscopy for structural determinations of toxin complexes, synthetic biology, and other technologies) to fully harness the biomedical potential of venoms and toxins from snakes and other animals.

Revealing the Diversity of Sequences, Structures, and Targets of Peptides from South China Sea Macrodactyla doreensis Based on Transcriptomics

The South China Sea is rich in sea anemone resources, and the protein and peptide components from sea anemone toxins comprise an important treasure trove for researchers to search for leading compounds. This study conducted a comprehensive transcriptomic analysis of the tentacles and column of Macrodactyla doreensis and explored the distribution and diversity of proteins and peptides in depth using bioinformatics, initially constructing a putative protein and peptide database. In this database, typical peptide families are identified through amino acid sequence analysis, and their 3D structures and potential biological activities are revealed through AlphaFold2 modeling and molecular docking. A total of 4239 transcripts were identified, of which the putative protein accounted for 81.53%. The highest content comprised immunoglobulin and a variety of proteases, mainly distributed in the column and related to biological functions. Importantly, the putative peptide accounted for 18.47%, containing ShK domain and Kunitz-type peptides, mainly distributed in the tentacles and related to offensive predatory behavior. Interestingly, 40 putative peptides belonging to eight typical peptide families were identified, and their structures and targets were predicted. This study reveals the diversity and complexity of Macrodactyla doreensis toxins and predicts their structure and targets based on amino acid sequences, providing a feasible approach for research regarding the discovery of peptides with potentially high activity.

Proteomic Analysis and Biochemical Characterization of the Nematocyst Extract of the Hydrozoan Velella velella

The venom contained within cnidarian nematocysts has a complex composition and holds significant potential for biotechnological applications. In this context, one of the most effective methods for studying nematocyst contents is the proteomic approach, which can detect even trace amounts of compounds while minimizing the need for large-scale animal collection, thus helping to preserve ecosystem integrity. This study aimed to provide a comprehensive proteomic and biochemical characterization of the crude nematocyst extract from the common hydrozoan Velella velella. Despite not being harmful to humans, the analysis of the crude venom extract from V. velella brought to the identification of 783 different proteins, categorized into structural components, enzymes, and potential toxins, revealing a qualitative composition of the venom similar to that of other more toxic cnidarians. Biochemical assays confirmed the presence of various active hydrolytic enzymes within the extract, including proteases, phospholipases, hyaluronidases, DNases, and chitinases. These findings pave the road for future studies involving the pharmacological applications of Velella velella venom components through recombinant production and functional testing.

Discovery, Expression, and In Silico Safety Evaluation of Honey Truffle Sweetener, a Sweet Protein Derived from Mattirolomyces terfezioides and Produced by Heterologous Expression in Komagataella phaffii

Honey truffle sweetener (HTS), a 121 amino acid protein is identified as a high-intensity sweetener found naturally occurring in the Hungarian Sweet Truffle Mattirolomyces terfezioides, an edible mushroom used in regional diets. The protein is intensely sweet, but the truffle is difficult to cultivate; therefore, the protein was systematically characterized, and the gene coding for the protein was expressed in a commonly used host yeast Komagataella phaffii. The heterologously expressed protein maintained the structural characteristics and sweet taste of the truffle. Preliminary safety evaluations for use as a food ingredient were performed on the protein including digestibility and in silico approaches for predicting the allergenicity and toxicity of the protein. HTS is predicted to be nonallergenic, nontoxic, and readily digestible. This protein is readily produced by precision fermentation of the host yeast, making it a potential replacement for both added sugars and small molecule high-intensity sweeteners in food.

Comparative Analysis of Tentacle Extract and Nematocyst Venom: Toxicity, Mechanism, and Potential Intervention in the Giant Jellyfish Nemopilema nomurai

The giant jellyfish Nemopilema nomurai sting can cause local and systemic reactions; however, comparative analysis of the tentacle extract (TE) and nematocyst venom extract (NV), and its toxicity, mechanism, and potential intervention are still limited. This study compared venom from TE and NV for their composition, toxicity, and efficacy in vitro and in vivo used RAW264.7 cells and ICR mice. A total of 239 and 225 toxin proteins were identified in TE and NV by proteomics, respectively. Pathological analysis revealed that TE and NV caused heart and liver damage through apoptosis, necrosis, and inflammation, while TE exhibited higher toxicity ex vivo and in vivo. Biochemical markers indicated TE and NV elevated creatine kinase, lactatedehydrogenase, and aspartate aminotransferase, with the TE group showing a more significant increase. Transcriptomics and Western blotting indicated both venoms increased cytokines expression and MAPK signaling pathways. Additionally, 1 mg/kg PACOCF3 (the phospholipase A2 inhibitor) improved survival from 16.7% to 75% in mice. Our results indicate that different extraction methods impact venom activities, tentacle autolysis preserves toxin proteins and their toxicity, and PACOCF3 is a potential antidote, which establishes a good extraction method of jellyfish venom, expands our understanding of jellyfish toxicity, mechanism, and provides a promising intervention.

Identification of molecular targets and small drug candidates for Huntington's disease via bioinformatics and a network-based screening approach

Huntington's disease (HD) is a gradually severe neurodegenerative ailment characterised by an increase of a specific trinucleotide repeat sequence (cytosine-adenine-guanine, CAG). It is passed down as a dominant characteristic that worsens over time, creating a significant risk. Despite being monogenetic, the underlying mechanisms as well as biomarkers remain poorly understood. Furthermore, early detection of HD is challenging, and the available diagnostic procedures have low precision and accuracy. The research was conducted to provide knowledge of the biomarkers, pathways and therapeutic targets involved in the molecular processes of HD using informatic based analysis and applying network-based systems biology approaches. The gene expression profile datasets GSE97100 and GSE74201 relevant to HD were studied. As a consequence, 46 differentially expressed genes (DEGs) were identified. 10 hub genes (TPM1, EIF2S3, CCN2, ACTN1, ACTG2, CCN1, CSRP1, EIF1AX, BEX2 and TCEAL5) were further differentiated in the protein-protein interaction (PPI) network. These hub genes were typically down-regulated. Additionally, DEGs-transcription factors (TFs) connections (e.g. GATA2, YY1 and FOXC1), DEG-microRNA (miRNA) interactions (e.g. hsa-miR-124-3p and has-miR-26b-5p) were also comprehensively forecast. Additionally, related gene ontology concepts (e.g. sequence-specific DNA binding and TF activity) connected to DEGs in HD were identified using gene set enrichment analysis (GSEA). Finally, in silico drug design was employed to find candidate drugs for the treatment HD, and while the possible modest therapeutic compounds (e.g. cortistatin A, 13,16-Epoxy-25-hydroxy-17-cheilanthen-19,25-olide, Hecogenin) against HD were expected. Consequently, the results from this study may give researchers useful resources for the experimental validation of Huntington's diagnosis and therapeutic approaches.

Molecular Display of the Animal Meta-Venome for Discovery of Novel Therapeutic Peptides

Animal venoms, distinguished by their unique structural features and potent bioactivities, represent a vast and relatively untapped reservoir of therapeutic molecules. However, limitations associated with extracting or expressing large numbers of individual venoms and venom-like molecules have precluded their therapeutic evaluation via high throughput screening. Here, we developed an innovative computational approach to design a highly diverse library of animal venoms and "metavenoms". We employed programmable M13 hyperphage display to preserve critical disulfide-bonded structures for highly parallelized single-round biopanning with quantitation via high-throughput DNA sequencing. Our approach led to the discovery of Kunitz type domain containing proteins that target the human itch receptor Mas-related G protein-coupled receptor X4 (MRGPRX4), which plays a crucial role in itch perception. Deep learning-based structural homology mining identified two endogenous human homologs, tissue factor pathway inhibitor (TFPI) and serine peptidase inhibitor, Kunitz type 2 (SPINT2), which exhibit agonist-dependent potentiation of MRGPRX4. Highly multiplexed screening of animal venoms and metavenoms is therefore a promising approach to uncover new drug candidates.

Identification and Evolutionary Analysis of the Widely Distributed CAP Superfamily in Spider Venom

Venom plays a crucial role in the defense and predation of venomous animals. Spiders (Araneae) are among the most successful predators and have a fascinating venom composition. Their venom mainly contains disulfide-rich peptides and large proteins. Here, we analyzed spider venom protein families, utilizing transcriptomic and genomic data, and highlighted their similarities and differences. We show that spiders have specific combinations of toxins for better predation and defense, typically comprising a core toxin expressed alongside several auxiliary toxins. Among them, the CAP superfamily is widely distributed and highly expressed in web-building Araneoidea spiders. Our analysis of evolutionary relationships revealed four subfamilies (subA-subD) of the CAP superfamily that differ in structure and potential functions. CAP proteins are composed of a conserved CAP domain and diverse C-terminal domains. CAP subC shares similar domains with the snake ion channel regulator svCRISP proteins, while CAP subD possesses a sequence similar to that of insect venom allergen 5 (Ag5). Furthermore, we show that gene duplication and selective expression lead to increased expression of CAP subD, making it a core member of the CAP superfamily. This study sheds light on the functional diversity of CAP subfamilies and their evolutionary history, which has important implications for fully understanding the composition of spider venom proteins and the core toxin components of web-building spiders.

Proteotransciptomics of the Most Popular Host Sea Anemone Entacmaea quadricolor Reveals Not All Toxin Genes Expressed by Tentacles Are Recruited into Its Venom Arsenal

While the unique symbiotic relationship between anemonefishes and sea anemones is iconic, it is still not fully understood how anemonefishes can withstand and thrive within the venomous environment of their host sea anemone. In this study, we used a proteotranscriptomics approach to elucidate the proteinaceous toxin repertoire from the most common host sea anemone, Entacmaea quadricolor. Although 1251 different toxin or toxin-like RNA transcripts were expressed in E. quadricolor tentacles (0.05% of gene clusters, 1.8% of expression) and 5375 proteins were detected in milked venom, only 4% of proteins detected in venom were putative toxins (230), and they only represent on average 14% of the normalised protein expression in the milked venom samples. Thus, most proteins in milked venom do not appear to have a toxin function. This work raises the perils of defining a dominant venom phenotype based on transcriptomics data alone in sea anemones, as we found that the dominant venom phenotype differs between the transcriptome and proteome abundance data. E. quadricolor venom contains a mixture of toxin-like proteins of unknown and known function. A newly identified toxin protein family, Z3, rich in conserved cysteines of unknown function, was the most abundant at the RNA transcript and protein levels. The venom was also rich in toxins from the Protease S1, Kunitz-type and PLA2 toxin protein families and contains toxins from eight venom categories. Exploring the intricate venom toxin components in other host sea anemones will be crucial for improving our understanding of how anemonefish adapt to the venomous environment.

Evolutionary Analysis of Cnidaria Small Cysteine-Rich Proteins (SCRiPs), an Enigmatic Neurotoxin Family from Stony Corals and Sea Anemones (Anthozoa: Hexacorallia)

Cnidarians (corals, sea anemones, and jellyfish) produce toxins that play central roles in key ecological processes, including predation, defense, and competition, being the oldest extant venomous animal lineage. Cnidaria small cysteine-rich proteins (SCRiPs) were the first family of neurotoxins detected in stony corals, one of the ocean's most crucial foundation species. Yet, their molecular evolution remains poorly understood. Moreover, the lack of a clear classification system has hindered the establishment of an accurate and phylogenetically informed nomenclature. In this study, we extensively surveyed 117 genomes and 103 transcriptomes of cnidarians to identify orthologous SCRiP gene sequences. We annotated a total of 168 novel putative SCRiPs from over 36 species of stony corals and 12 species of sea anemones. Phylogenetic reconstruction identified four distinct SCRiP subfamilies, according to strict discrimination criteria based on well-supported monophyly with a high percentage of nucleotide and amino acids' identity. Although there is a high prevalence of purifying selection for most SCRiP subfamilies, with few positively selected sites detected, a subset of Acroporidae sequences is influenced by diversifying positive selection, suggesting potential neofunctionalizations related to the fine-tuning of toxin potency. We propose a new nomenclature classification system relying on the phylogenetic distribution and evolution of SCRiPs across Anthozoa, which will further assist future proteomic and functional research efforts.

Venomics Reveals the Venom Complexity of Sea Anemone Heteractis magnifica

The venoms of various sea anemones are rich in diverse toxins, which usually play a dual role in capturing prey and deterring predators. However, the complex components of such venoms have not been well known yet. Here, venomics of integrating transcriptomic and proteomic technologies was applied for the first time to identify putative protein and peptide toxins from different tissues of the representative sea anemone, Heteractis magnifica. The transcriptomic analysis of H. magnifica identified 728 putative toxin sequences, including 442 and 381 from the tentacles and the column, respectively, and they were assigned to 68 gene superfamilies. The proteomic analysis confirmed 101 protein and peptide toxins in the venom, including 91 in the tentacles and 39 in the column. The integrated venomics also confirmed that some toxins such as the ShK-like peptides and defensins are co-expressed in both the tentacles and the column. Meanwhile, a homology analysis was conducted to predict the three-dimensional structures and potential activity of seven representative toxins. Altogether, this venomics study revealed the venom complexity of H. magnifica, which will help deepen our understanding of cnidarian toxins, thereby supporting the in-depth development of valuable marine drugs.

ToxCodAn-Genome: an automated pipeline for toxin-gene annotation in genome assembly of venomous lineages

BACKGROUND

The rapid development of sequencing technologies resulted in a wide expansion of genomics studies using venomous lineages. This facilitated research focusing on understanding the evolution of adaptive traits and the search for novel compounds that can be applied in agriculture and medicine. However, the toxin annotation of genomes is a laborious and time-consuming task, and no consensus pipeline is currently available. No computational tool currently exists to address the challenges specific to toxin annotation and to ensure the reproducibility of the process.

RESULTS

Here, we present ToxCodAn-Genome, the first software designed to perform automated toxin annotation in genomes of venomous lineages. This pipeline was designed to retrieve the full-length coding sequences of toxins and to allow the detection of novel truncated paralogs and pseudogenes. We tested ToxCodAn-Genome using 12 genomes of venomous lineages and achieved high performance on recovering their current toxin annotations. This tool can be easily customized to allow improvements in the final toxin annotation set and can be expanded to virtually any venomous lineage. ToxCodAn-Genome is fast, allowing it to run on any personal computer, but it can also be executed in multicore mode, taking advantage of large high-performance servers. In addition, we provide a guide to direct future research in the venomics field to ensure a confident toxin annotation in the genome being studied. As a case study, we sequenced and annotated the toxin repertoire of Bothrops alternatus, which may facilitate future evolutionary and biomedical studies using vipers as models.

CONCLUSIONS

ToxCodAn-Genome is suitable to perform toxin annotation in the genome of venomous species and may help to improve the reproducibility of further studies. ToxCodAn-Genome and the guide are freely available at https://github.com/pedronachtigall/ToxCodAn-Genome.

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.

Metatranscriptome analysis reveals the putative venom toxin repertoire of the biofouling hydroid Ectopleura larynx

Cnidarians thriving in biofouling communities on aquaculture net pens represent a significant health risk for farmed finfish due to their stinging cells. The toxins coming into contact with the fish, during net cleaning, can adversely affect their behavior, welfare, and survival, with a particularly serious health risk for the gills, causing direct tissue damage such as formation of thrombi and increasing risks of secondary infections. The hydroid Ectopleura larynx is one of the most common fouling organisms in Northern Europe. However, despite its significant economic, environmental, and operational impact on finfish aquaculture, biological information on this species is scarce and its venom composition has never been investigated. In this study, we generated a whole transcriptome of E. larynx, and identified its putative expressed venom toxin proteins (predicted toxin proteins, not functionally characterized) based on in silico transcriptome annotation mining and protein sequence analysis. The results uncovered a broad and diverse repertoire of putative toxin proteins for this hydroid species. Its toxic arsenal appears to include a wide and complex selection of toxin proteins, covering a large panel of potential biological functions that play important roles in envenomation. The putative toxins identified in this species, such as neurotoxins, GTPase toxins, metalloprotease toxins, ion channel impairing toxins, hemorrhagic toxins, serine protease toxins, phospholipase toxins, pore-forming toxins, and multifunction toxins may cause various major deleterious effects in prey, predators, and competitors. These results provide valuable new insights into the venom composition of cnidarians, and venomous marine organisms in general, and offer new opportunities for further research into novel and valuable bioactive molecules for medicine, agronomics and biotechnology.

Pharmacoinformatic Investigation of Silymarin as a Potential Inhibitor against Nemopilema nomurai Jellyfish Metalloproteinase Toxin-like Protein

Jellyfish stings pose a major threat to swimmers and fishermen worldwide. These creatures have explosive cells containing one large secretory organelle called a nematocyst in their tentacles, which contains venom used to immobilize prey. Nemopilema nomurai, a venomous jellyfish belonging to the phylum Cnidaria, produces venom (NnV) comprising various toxins known for their lethal effects on many organisms. Of these toxins, metalloproteinases (which belong to the toxic protease family) play a significant role in local symptoms such as dermatitis and anaphylaxis, as well as systemic reactions such as blood coagulation, disseminated intravascular coagulation, tissue injury, and hemorrhage. Hence, a potential metalloproteinase inhibitor (MPI) could be a promising candidate for reducing the effects of venom toxicity. For this study, we retrieved the Nemopilema nomurai venom metalloproteinase sequence (NnV-MPs) from transcriptome data and modeled its three-dimensional structure using AlphaFold2 in a Google Colab notebook. We employed a pharmacoinformatics approach to screen 39 flavonoids and identify the most potent inhibitor against NnV-MP. Previous studies have demonstrated the efficacy of flavonoids against other animal venoms. Based on our analysis, Silymarin emerged as the top inhibitor through ADMET, docking, and molecular dynamics analyses. In silico simulations provide detailed information on the toxin and ligand binding affinity. Our results demonstrate that Silymarin's strong inhibitory effect on NnV-MP is driven by hydrophobic affinity and optimal hydrogen bonding. These findings suggest that Silymarin could serve as an effective inhibitor of NnV-MP, potentially reducing the toxicity associated with jellyfish envenomation.

A bivalent remipede toxin promotes calcium release via ryanodine receptor activation

Multivalent ligands of ion channels have proven to be both very rare and highly valuable in yielding unique insights into channel structure and pharmacology. Here, we describe a bivalent peptide from the venom of Xibalbanus tulumensis, a troglobitic arthropod from the enigmatic class Remipedia, that causes persistent calcium release by activation of ion channels involved in muscle contraction. The high-resolution solution structure of φ-Xibalbin3-Xt3a reveals a tandem repeat arrangement of inhibitor-cysteine knot (ICK) domains previously only found in spider venoms. The individual repeats of Xt3a share sequence similarity with a family of scorpion toxins that target ryanodine receptors (RyR). Single-channel electrophysiology and quantification of released Ca²⁺ stores within skinned muscle fibers confirm Xt3a as a bivalent RyR modulator. Our results reveal convergent evolution of RyR targeting toxins in remipede and scorpion venoms, while the tandem-ICK repeat architecture is an evolutionary innovation that is convergent with toxins from spider venoms.

ScrepYard: An online resource for disulfide-stabilized tandem repeat peptides

Receptor avidity through multivalency is a highly sought-after property of ligands. While readily available in nature in the form of bivalent antibodies, this property remains challenging to engineer in synthetic molecules. The discovery of several bivalent venom peptides containing two homologous and independently folded domains (in a tandem repeat arrangement) has provided a unique opportunity to better understand the underpinning design of multivalency in multimeric biomolecules, as well as how naturally occurring multivalent ligands can be identified. In previous work, we classified these molecules as a larger class termed secreted cysteine-rich repeat-proteins (SCREPs). Here, we present an online resource; ScrepYard, designed to assist researchers in identification of SCREP sequences of interest and to aid in characterizing this emerging class of biomolecules. Analysis of sequences within the ScrepYard reveals that two-domain tandem repeats constitute the most abundant SCREP domain architecture, while the interdomain "linker" regions connecting the functional domains are found to be abundant in amino acids with short or polar sidechains and contain an unusually high abundance of proline residues. Finally, we demonstrate the utility of ScrepYard as a virtual screening tool for discovery of putatively multivalent peptides, by using it as a resource to identify a previously uncharacterized serine protease inhibitor and confirm its predicted activity using an enzyme assay.

Function-based classification of hazardous biological sequences: Demonstration of a new paradigm for biohazard assessments

Bioengineering applies analytical and engineering principles to identify functional biological building blocks for biotechnology applications. While these building blocks are leveraged to improve the human condition, the lack of simplistic, machine-readable definition of biohazards at the function level is creating a gap for biosafety practices. More specifically, traditional safety practices focus on the biohazards of known pathogens at the organism-level and may not accurately consider novel biodesigns with engineered functionalities at the genetic component-level. This gap is motivating the need for a paradigm shift from organism-centric procedures to function-centric biohazard identification and classification practices. To address this challenge, we present a novel methodology for classifying biohazards at the individual sequence level, which we then compiled to distinguish the biohazardous property of pathogenicity at the whole genome level. Our methodology is rooted in compilation of hazardous functions, defined as a set of sequences and associated metadata that describe coarse-level functions associated with pathogens (e.g., adherence, immune subversion). We demonstrate that the resulting database can be used to develop hazardous “fingerprints” based on the functional metadata categories. We verified that these hazardous functions are found at higher levels in pathogens compared to non-pathogens, and hierarchical clustering of the fingerprints can distinguish between these two groups. The methodology presented here defines the hazardous functions associated with bioengineering functional building blocks at the sequence level, which provide a foundational framework for classifying biological hazards at the organism level, thus leading to the improvement and standardization of current biosecurity and biosafety practices.

Reliability of the In Silico Prediction Approach to In Vitro Evaluation of Bacterial Toxicity

Several pathogens that spread through the air are highly contagious, and related infectious diseases are more easily transmitted through airborne transmission under indoor conditions, as observed during the COVID-19 pandemic. Indoor air contaminated by microorganisms, including viruses, bacteria, and fungi, or by derived pathogenic substances, can endanger human health. Thus, identifying and analyzing the potential pathogens residing in the air are crucial to preventing disease and maintaining indoor air quality. Here, we applied deep learning technology to analyze and predict the toxicity of bacteria in indoor air. We trained the ProtBert model on toxic bacterial and virulence factor proteins and applied them to predict the potential toxicity of some bacterial species by analyzing their protein sequences. The results reflect the results of the in vitro analysis of their toxicity in human cells. The in silico-based simulation and the obtained results demonstrated that it is plausible to find possible toxic sequences in unknown protein sequences.

A Novel Potent Crystalline Chitin Decomposer: Chitin Deacetylase from Acinetobacter schindleri MCDA01

Chitosan is a functional ingredient that is widely used in food chemistry as an emulsifier, flocculant, antioxidant, or preservative. Chitin deacetylases (CDAs) can catalyze the hydrolysis of acetyl groups, making them useful in the clean production of chitosan. However, the high inactivity of crystalline chitin catalyzed by CDAs has been regarded as the technical bottleneck of crystalline chitin deacetylation. Here, we mined the AsCDA gene from the genome of Acinetobacter schindleri MCDA01 and identified a member of the uraD_N-term-dom superfamily, which was a novel chitin deacetylase with the highest deacetylation activity. The AsCDA gene was expressed in Escherichia coli BL21 by IPTG induction, whose activity to colloidal chitin, α-chitin, and β-chitin reached 478.96 U/mg, 397.07 U/mg, and 133.27 U/mg, respectively. In 12 h, the enzymatic hydrolysis of AsCDA removed 63.05% of the acetyl groups from α-chitin to prepare industrial chitosan with a degree of deacetylation higher than 85%. AsCDA, as a potent chitin decomposer in the production of chitosan, plays a positive role in the upgrading of the chitosan industry and the value-added utilization of chitin biological resources.

Investigation of Peptide Toxin Diversity in Ribbon Worms (Nemertea) Using a Transcriptomic Approach

Nemertea is a phylum of nonsegmented worms (supraphylum: Spiralia), also known as ribbon worms. The members of this phylum contain various toxins, including peptide toxins. Here, we provide a transcriptomic analysis of peptide toxins in 14 nemertean species, including Cephalothrix cf. simula, which was sequenced in the current study. The summarized data show that the number of toxin transcripts in the studied nemerteans varied from 12 to 82. The most represented groups of toxins were enzymes and ion channel inhibitors, which, in total, reached a proportion of 72% in some species, and the least represented were pore-forming toxins and neurotoxins, the total proportion of which did not exceed 18%. The study revealed that nemerteans possess a much greater variety of toxins than previously thought and showed that these animals are a promising object for the investigation of venom diversity and evolution, and in the search for new peptide toxins.

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.

A Neuropeptide Y/F-like Polypeptide Derived from the Transcriptome of Turbinaria peltata Suppresses LPS-Induced Astrocytic Inflammation

Neuropeptides are a group of neuronal signaling molecules that regulate physiological and behavioral processes in animals. Here, we used in silico mining to predict the polypeptide composition of available transcriptomic data of Turbinaria peltata. In total, 118 transcripts encoding putative peptide precursors were discovered. One neuropeptide Y/F-like peptide, named TpNPY, was identified and selected for in silico structural, in silico binding, and pharmacological studies. In our study, the anti-inflammation effect of TpNPY was evaluated using an LPS-stimulated C8-D1A astrocyte cell model. Our results demonstrated that TpNPY, at 0.75-3 μM, inhibited LPS-induced NO production and reduced the expression of iNOS in a dose-dependent manner. Furthermore, TpNPY reduced the secretion of proinflammatory cytokines. Additionally, treatment with TpNPY reduced LPS-mediated elevation of ROS production and the intracellular calcium concentration. Further investigation revealed that TpNPY downregulated the IKK/IκB/NF-κB signaling pathway and inhibited expression of the NLRP3 inflammasome. Through molecular docking and using an NPY receptor antagonist, TpNPY was shown to have the ability to interact with the NPY Y1 receptor. On the basis of these findings, we concluded that TpNPY might prevent LPS-induced injury in astrocytes through activation of the NPY-Y1R.

Computer-Aided Analysis of West Sub-Saharan Africa Snakes Venom towards the Design of Epitope-Based Poly-Specific Antivenoms

Snakebite envenomation is a neglected tropical disease that causes over 100,000 deaths each year. The only effective treatment consists of antivenoms derived from animal sera, but these have been deemed with highly variable potency and are usually inaccessible and too costly for victims. The production of antivenoms by venom-independent techniques, such as the immunization with multi-epitope constructs, could circumvent those drawbacks. Herein, we present a knowledge-based pipeline to prioritize potential epitopes of therapeutic relevance from toxins of medically important snakes in West Sub-Saharan Africa. It is mainly based on sequence conservation and protein structural features. The ultimately selected 41 epitopes originate from 11 out of 16 snake species considered of highest medical importance in the region and 3 out of 10 of those considered as secondary medical importance. Echis ocellatus, responsible for the highest casualties in the area, would be covered by 12 different epitopes. Remarkably, this pipeline is versatile and customizable for the analysis of snake venom sequences from any other region of the world.

Venom system variation and the division of labor in the colonial hydrozoan Hydractinia symbiolongicarpus

Cnidarians (jellyfish, hydroids, sea anemones, and corals) possess a unique method for venom production, maintenance, and deployment through a decentralized system composed of different types of venom-filled stinging structures called nematocysts. In many species, nematocyst types are distributed heterogeneously across functionally distinct tissues. This has led to a prediction that different nematocyst types contain specific venom components. The colonial hydrozoan, Hydractinia symbiolongicarpus, is an ideal system to study the functional distribution of nematocyst types and their venoms, given that they display a division of labor through functionally distinct polyps within the colony. Here, we characterized the composition and distribution of nematocysts (cnidome) in the different polyp types and show that the feeding polyp (gastrozooid) has a distinct cnidome compared to the reproductive (gonozooid) and predatory polyp (dactylozooid). We generated a nematocyst-specific reporter line to track nematocyst development (nematogenesis) in H. symbiolongicarpus, and were able to confirm that nematogenesis primarily occurs in the mid-region of the gastrozooid and throughout stolons (tubes of epithelia that connect the polyps in the colony). This reporter line enabled us to isolate a nematocyst-specific lineage of cells for de novo transcriptome assembly, annotate venom-like genes (VLGs) and determine differential expression (DE) across polyp types. We show that a majority of VLGs are upregulated in gastrozooids, consistent with it being the primary site of active nematogenesis. However, despite gastrozooids producing more nematocysts, we found a number of VLGs significantly upregulated in dactylozooids, suggesting that these VLGs may be important for prey-capture. Our transgenic Hydractinia reporter line provides an opportunity to explore the complex interplay between venom composition, nematocyst diversity, and ecological partitioning in a colonial hydrozoan that displays a division of labor.

•

Functionally specific polyp types in Hydractinia symbiolongicarpus have distinct cnidomes.

•

We present a nematocyst-targeted transgenic line for H. symbiolongicarpus, showcasing active areas of nematogenesis.

•

105 venom-like genes (VLGs) were annotated from an assembled nematocyst-enriched transcriptome.

•

Several VLGs were significantly upregulated in feeding polyps, consistent with being a site of active nematogenesis.

•

Differential expression analysis suggests that different polyp types express distinct combinations of VLGs.

Venom profile of the European carpenter bee Xylocopa violacea: Evolutionary and applied considerations on its toxin components

Modern venomics is increasing its focus on hymenopterans such as honeybees, bumblebees, parasitoid wasps, ants and true wasps. However solitary bees remain understudied in comparison and the few available venom studies focus on short melittin-like sequences and antimicrobial peptides. Herein we describe the first comprehensive venom profile of a solitary bee, the violet carpenter bee Xylocopa violacea, by using proteo-transcriptomics. We reveal a diverse and complex venom profile with 43 different protein families identified from dissected venom gland extracts of which 32 are also detected in the defensively injected venom. Melittin and apamin are the most highly secreted components, followed by Phospholipase A2, Icarapin, Secapin and three novel components. Other components, including eight novel protein families, are rather lowly expressed. We further identify multiple forms of apamin-like peptides. The melittin-like sequences of solitary bees separate into two clades, one comprised most sequences from solitary bees including xylopin (the variant in Xylocopa), while sequences from Lasioglossa appear closer related to melittin-like peptides from Bombus (Bombolittins). Our study suggests that more proteo-transcriptomic data from other solitary bees should be complemented with corresponding genome data to fully understand the evolution and complexity of bee venom proteins, and is of a particular need to disentangle the ambiguous phylogenetic relations of short peptides.

Transcriptomic Insights into the Diversity and Evolution of Myxozoa (Cnidaria, Endocnidozoa) Toxin-like Proteins

Myxozoa is a speciose group of endoparasitic cnidarians that can cause severe ecological and economic effects. Their cnidarian affinity is affirmed by genetic relatedness and the presence of nematocysts, historically called "polar capsules". Previous studies have revealed the presence of toxin-like proteins in myxozoans; however, the diversity and evolution of venom in Myxozoa are not fully understood. Here, we performed a comparative analysis using the newly sequenced transcriptomes of five Myxobolidae species as well as some public datasets. Toxin mining revealed that myxozoans have lost most of their toxin families, while most species retained Kunitz, M12B, and CRISP, which may play a role in endoparasitism. The venom composition of Endocnidozoa (Myxozoa + Polypodium) differs from that of free-living cnidarians and may be influenced by ecological and environmental factors. Phylogenetic analyses showed that toxin families of myxozoans and free-living cnidarians were clustered into different clades. Selection analyses showed that purifying selection was the dominant evolutionary pressure in toxins, while they were still influenced by episodic adaptive selection. This suggests that the potency or specificity of a particular toxin or species might increase. Overall, our findings provide a more comprehensive framework for understanding the diversity and evolution of Myxozoa venoms.

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.

Functional, immunological characterization, and anticancer activity of BaMtx: A new Lys49- PLA2 homologue isolated from the venom of Peruvian Bothrops atrox snake (Serpentes: Viperidae)

Bothorps atrox is responsible for most of the ophidism cases in Perú. As part of the envenoming, myotoxicity is one of the most recurrent and destructive effects. In this study, a myotoxin, named BaMtx, was purified from B. atrox venom to elucidate its biological, immunological, and molecular characteristics. BaMtx was purified using CM-Sephadex-C-25 ion-exchange resin and SDS-PAGE analysis showed a unique protein band of 13 kDa or 24 kDa under reducing or non-reducing conditions, respectively. cDNA sequence codified a 122-aa mature protein with high homology with other Lys49-PLA₂s; modeled structure showed a N-terminal helix, a β-wing region, and a C-terminal random coil. This protein has a poor phospholipase A2 enzymatic activity. BaMtx has myotoxic (DMM = 12.30 ± 0.95 μg) and edema-forming (DEM = 26.00 ± 1.15 μg) activities. Rabbit immunization with purified enzyme produced anti-BaMtx antibodies that reduced 50.28 ± 10.15% of myotoxic activity and showed significant cross-reactivity against B. brazili and B pictus venoms. On the other hand, BaMtx exhibits mild anti-proliferative and anti-migratory effects on breast cancer cells, affecting the ROS and NADH levels, which may reduce mitochondrial respiration. These results contribute to the understanding of B. atrox Lys49-PLA2 effects and establish the anticancer potential de BaMtx.

Toxinology in the proteomics era: a review on arachnid venom proteomics

The word venomics was coined to acknowledge the studies that use omics to investigate venom proteins and peptides. Venomics has evolved considerably over the last 20 years. The first works on scorpion or spider venomics were published in the early 2000's. Such studies relied on peptide mass fingerprinting (PMF) to characterize venom complexity. After the introduction of new mass spectrometers with higher resolution, sensitivity and mass accuracy, and the next-generation nucleotide sequencing, the complexity of data reported in research on scorpion and spider venomics increased exponentially, which allowed more comprehensive studies. In the present review article, we covered key publications on scorpion venomics and spider venomics, presenting historical grounds and implemented technologies over the last years. The literature presented in this review was selected after searching the PubMed database using the terms "(scorpion venom) AND (proteome)" for scorpion venomics, and "(spider venom) AND (proteome)" for publications on spider venomics. We presented the key aspects related to proteomics in the covered papers including, but not restricted to, the employed proteomic strategy (i.e., PMF, two-dimensional gel electrophoresis, shotgun/bottom-up and/or top-down/peptidome), and the type of mass spectrometer used. Some conclusions can be drawn from the present study. For example, the scorpion genus Tityus is the most studied concerning venomics, followed by Centruroides; whereas for spiders the studied genera were found more equally distributed. Another interesting conclusion is the lack of high throughput studies on post-translational modifications (PTMs) of scorpion and spider proteins. In our opinion, PTMs should be more studied as they can modulate the activity of scorpion and spider toxins.

Quantitative Insights into the Contribution of Nematocysts to the Adaptive Success of Cnidarians Based on Proteomic Analysis

Nematocysts are secretory organelles in cnidarians that play important roles in predation, defense, locomotion, and host invasion. However, the extent to which nematocysts contribute to adaptation and the mechanisms underlying nematocyst evolution are unclear. Here, we investigated the role of the nematocyst in cnidarian evolution based on eight nematocyst proteomes and 110 cnidarian transcriptomes/genomes. We detected extensive species-specific adaptive mutations in nematocyst proteins (NEMs) and evidence for decentralized evolution, in which most evolutionary events involved non-core NEMs, reflecting the rapid diversification of NEMs in cnidarians. Moreover, there was a 33-55 million year macroevolutionary lag between nematocyst evolution and the main phases of cnidarian diversification, suggesting that the nematocyst can act as a driving force in evolution. Quantitative analysis revealed an excess of adaptive changes in NEMs and enrichment for positively selected conserved NEMs. Together, these findings suggest that nematocysts may be key to the adaptive success of cnidarians and provide a reference for quantitative analyses of the roles of phenotypic novelties in adaptation.

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.

Stingray Venom Proteins: Mechanisms of Action Revealed Using a Novel Network Pharmacology Approach

Animal venoms offer a valuable source of potent new drug leads, but their mechanisms of action are largely unknown. We therefore developed a novel network pharmacology approach based on multi-omics functional data integration to predict how stingray venom disrupts the physiological systems of target animals. We integrated 10 million transcripts from five stingray venom transcriptomes and 848,640 records from three high-content venom bioactivity datasets into a large functional data network. The network featured 216 signaling pathways, 29 of which were shared and targeted by 70 transcripts and 70 bioactivity hits. The network revealed clusters for single envenomation outcomes, such as pain, cardiotoxicity and hemorrhage. We carried out a detailed analysis of the pain cluster representing a primary envenomation symptom, revealing bibrotoxin and cholecystotoxin-like transcripts encoding pain-inducing candidate proteins in stingray venom. The cluster also suggested that such pain-inducing toxins primarily activate the inositol-3-phosphate receptor cascade, inducing intracellular calcium release. We also found strong evidence for synergistic activity among these candidates, with nerve growth factors cooperating with the most abundant translationally-controlled tumor proteins to activate pain signaling pathways. Our network pharmacology approach, here applied to stingray venom, can be used as a template for drug discovery in neglected venomous species.

Peptide Inhibitors of Kv1.5: An Option for the Treatment of Atrial Fibrillation

The human voltage gated potassium channel Kv1.5 that conducts the IKur current is a key determinant of the atrial action potential. Its mutations have been linked to hereditary forms of atrial fibrillation (AF), and the channel is an attractive target for the management of AF. The development of IKur blockers to treat AF resulted in small molecule Kv1.5 inhibitors. The selectivity of the blocker for the target channel plays an important role in the potential therapeutic application of the drug candidate: the higher the selectivity, the lower the risk of side effects. In this respect, small molecule inhibitors of Kv1.5 are compromised due to their limited selectivity. A wide range of peptide toxins from venomous animals are targeting ion channels, including mammalian channels. These peptides usually have a much larger interacting surface with the ion channel compared to small molecule inhibitors and thus, generally confer higher selectivity to the peptide blockers. We found two peptides in the literature, which inhibited IKur: Ts6 and Osu1. Their affinity and selectivity for Kv1.5 can be improved by rational drug design in which their amino acid sequences could be modified in a targeted way guided by in silico docking experiments.

Deep-Sea Anemones Are Prospective Source of New Antimicrobial and Cytotoxic Compounds

The peculiarities of the survival and adaptation of deep-sea organisms raise interest in the study of their metabolites as promising drugs. In this work, the hemolytic, cytotoxic, antimicrobial, and enzyme-inhibitory activities of tentacle extracts from five species of sea anemones (Cnidaria, orders Actiniaria and Corallimorpharia) collected near the Kuril and Commander Islands of the Far East of Russia were evaluated for the first time. The extracts of Liponema brevicorne and Actinostola callosa demonstrated maximal hemolytic activity, while high cytotoxic activity against murine splenocytes and Ehrlich carcinoma cells was found in the extract of Actinostola faeculenta. The extracts of Corallimorphus cf. pilatus demonstrated the greatest activity against Ehrlich carcinoma cells but were not toxic to mouse spleen cells. Sea anemones C. cf. pilatus and Stomphia coccinea are promising sources of antimicrobial and antifungal compounds, being active against Gram-positive bacteria Bacillus subtilis, Staphylococcus aureus, and yeast Candida albicans. Moreover, all sea anemones contain α-galactosidase inhibitors. Peptide mass fingerprinting of L. brevicorne and C. cf. pilatus extracts provided a wide range of peptides, predominantly with molecular masses of 4000–5900 Da, which may belong to a known or new structural class of toxins. The obtained data allow concluding that deep-sea anemones are a promising source of compounds for drug discovery.

Categorizing sequences of concern by function to better assess mechanisms of microbial pathogenesis

To identify sequences with a role in microbial pathogenesis, we assessed the adequacy of their annotation by existing controlled vocabularies and sequence databases. Our goal was to regularize descriptions of microbial pathogenesis for improved integration with bioinformatic applications. Here, we review the challenges of annotating sequences for pathogenic activity. We relate the categorization of more than 2,750 sequences of pathogenic microbes through a controlled vocabulary called Functions of Sequences of Concern (FunSoCs). These allow for an ease of description by both humans and machines. We provide a subset of 220 fully annotated sequences in the supplemental material as examples. The use of this compact (∼30 terms), controlled vocabulary has potential benefits for research in microbial genomics, public health, biosecurity, biosurveillance, and the characterization of new and emerging pathogens.

Current Approaches in Supersecondary Structures Investigation

Proteins expressed during the cell cycle determine cell function, topology, and responses to environmental influences. The development and improvement of experimental methods in the field of structural biology provide valuable information about the structure and functions of individual proteins. This work is devoted to the study of supersecondary structures of proteins and determination of their structural motifs, description of experimental methods for their detection, databases, and repositories for storage, as well as methods of molecular dynamics research. The interest in the study of supersecondary structures in proteins is due to their autonomous stability outside the protein globule, which makes it possible to study folding processes, conformational changes in protein isoforms, and aberrant proteins with high productivity.

Proteo-Transcriptomic Characterization of Sirex nitobei (Hymenoptera: Siricidae) Venom

The wood-boring woodwasp Sirex nitobei is a native pest in Asia, infecting and weakening the host trees in numerous ecological and commercial coniferous forest plantations. In China, hosts of S. nitobei are diverse, so the pest has spread to several provinces of China, resulting in considerable economic and ecological damage. During female oviposition, S. nitobei venom along with arthrospores of the symbiotic fungus Amylostereum areolatum or A. chaetica is injected into host trees, and the combination of these two biological factors causes the death of xylem host trees. The presence of venom alone causes only the yellowing and wilting of needles. In this study, we constructed the venom gland transcriptome of S. nitobei for the first time and a total of 15,036 unigenes were acquired. From the unigenes, 11,560 ORFs were identified and 537 encoding protein sequences with signal peptides at the N-terminus. Then, we used the venomics approach to characterize the venom composition of female S. nitobei and predicted 1095 proteins by liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis. We focused on seven proteins that were both highly expressed in the venom gland transcriptome and predicted in the crude venom proteome. These seven proteins are laccase-2, laccase-3, a protein belonging to the Kazal family, chitooligosaccharidolytic β-N-acetylglucosaminidase, beta-galactosidase, icarapin-like protein, and waprin-Thr1-like protein. Using quantitative real-time PCR (qRT-PCR), we also proved that the genes related to these seven proteins are specifically expressed in the venom glands. Finally, we revealed the functional role of S. nitobei venom in the physiological response of host trees. It can not only promote the colonization of symbiotic fungus but contribute to the development of eggs and larvae. This study provides a deeper understanding of the molecular mechanism of the woodwasp–pine interaction.

ToxDL: deep learning using primary structure and domain embeddings for assessing protein toxicity

MOTIVATION

Genetically engineering food crops involves introducing proteins from other species into crop plant species or modifying already existing proteins with gene editing techniques. In addition, newly synthesized proteins can be used as therapeutic protein drugs against diseases. For both research and safety regulation purposes, being able to assess the potential toxicity of newly introduced/synthesized proteins is of high importance.

RESULTS

In this study, we present ToxDL, a deep learning-based approach for in silico prediction of protein toxicity from sequence alone. ToxDL consists of (i) a module encompassing a convolutional neural network that has been designed to handle variable-length input sequences, (ii) a domain2vec module for generating protein domain embeddings and (iii) an output module that classifies proteins as toxic or non-toxic, using the outputs of the two aforementioned modules. Independent test results obtained for animal proteins and cross-species transferability results obtained for bacteria proteins indicate that ToxDL outperforms traditional homology-based approaches and state-of-the-art machine-learning techniques. Furthermore, through visualizations based on saliency maps, we are able to verify that the proposed network learns known toxic motifs. Moreover, the saliency maps allow for directed in silico modification of a sequence, thus making it possible to alter its predicted protein toxicity.

AVAILABILITY AND IMPLEMENTATION

ToxDL is freely available at http://www.csbio.sjtu.edu.cn/bioinf/ToxDL/. The source code can be found at https://github.com/xypan1232/ToxDL.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Quantitative trait locus analysis of parasitoid counteradaptation to symbiont-conferred resistance

Insect hosts and parasitoids are engaged in an intense struggle of antagonistic coevolution. Infection with heritable bacterial endosymbionts can substantially increase the resistance of aphids to parasitoid wasps, which exerts selection on parasitoids to overcome this symbiont-conferred protection (counteradaptation). Experimental evolution in the laboratory has produced counteradapted populations of the parasitoid wasp Lysiphlebus fabarum. These populations can parasitize black bean aphids (Aphis fabae) protected by the bacterial endosymbiont Hamiltonella defensa, which confers high resistance against L. fabarum. We used two experimentally evolved parasitoid populations to study the genetic architecture of the counteradaptation to symbiont-conferred resistance by QTL analysis. With simple crossing experiments, we showed that the counteradaptation is a recessive trait depending on the maternal genotype. Based on these results, we designed a customized crossing scheme to genotype a mapping population phenotyped for the ability to parasitize Hamiltonella-protected aphids. Using 1835 SNP markers obtained by ddRAD sequencing, we constructed a high-density linkage map consisting of six linkage groups (LGs) with an overall length of 828.3 cM and an average marker spacing of 0.45 cM. We identified a single QTL associated with the counteradaptation to Hamiltonella in L. fabarum on linkage group 2. Out of 120 genes located in this QTL, several genes encoding putative venoms may represent candidates for counteradaptation, as parasitoid wasps inject venoms into their hosts during oviposition.

Two reference-quality sea snake genomes reveal their divergent evolution of adaptive traits and venom systems

True sea snakes (Hydrophiini) are among the last and most successful clades of vertebrates that show secondary marine adaptation, exhibiting diverse phenotypic traits and lethal venom systems. To better understand their evolution, we generated the first chromosome-level genomes of two representative Hydrophiini snakes, Hydrophis cyanocinctus and H. curtus. Through comparative genomics we identified a great expansion of the underwater olfaction-related V2R gene family, consisting of more than 1,000 copies in both snakes. A series of chromosome rearrangements and genomic structural variations were recognized, including large inversions longer than 30 megabase (Mb) on sex chromosomes which potentially affect key functional genes associated with differentiated phenotypes between the two species. By integrating multiomics we found a significant loss of the major weapon for elapid predation, three-finger toxin genes, which displayed a dosage effect in H. curtus. These genetic changes may imply mechanisms that drove the divergent evolution of adaptive traits including prey preferences between the two closely related snakes. Our reference-quality sea snake genomes also enrich the repositories for addressing important issues on the evolution of marine tetrapods, and provide a resource for discovering marine-derived biological products.

TISIGNER.com: web services for improving recombinant protein production

Experiments that are planned using accurate prediction algorithms will mitigate failures in recombinant protein production. We have developed TISIGNER (https://tisigner.com) with the aim of addressing technical challenges to recombinant protein production. We offer three web services, TIsigner (Translation Initiation coding region designer), SoDoPE (Soluble Domain for Protein Expression) and Razor, which are specialised in synonymous optimisation of recombinant protein expression, solubility and signal peptide analysis, respectively. Importantly, TIsigner, SoDoPE and Razor are linked, which allows users to switch between the tools when optimising genes of interest.

A Chromosome-Level Genome Assembly of the Mandarin Fish (Siniperca chuatsi)

The mandarin fish, Siniperca chuatsi, is an economically important perciform species with widespread aquaculture practices in China. Its special feeding habit, acceptance of only live prey fishes, contributes to its delicious meat. However, little is currently known about related genetic mechanisms. Here, we performed whole-genome sequencing and assembled a 758.78 Mb genome assembly of the mandarin fish, with the scaffold and contig N50 values reaching 2.64 Mb and 46.11 Kb, respectively. Approximately 92.8% of the scaffolds were ordered onto 24 chromosomes (Chrs) with the assistance of a previously established genetic linkage map. The chromosome-level genome contained 19,904 protein-coding genes, of which 19,059 (95.75%) genes were functionally annotated. The special feeding behavior of mandarin fish could be attributable to the interaction of a variety of sense organs (such as vision, smell, and endocrine organs). Through comparative genomics analysis, some interesting results were found. For example, olfactory receptor (OR) genes (especially the beta and delta types) underwent a significant expansion, and endocrinology/vision related npy, spexin, and opsin genes presented various functional mutations. These may contribute to the special feeding habit of the mandarin fish by strengthening the olfactory and visual systems. Meanwhile, previously identified sex-related genes and quantitative trait locis (QTLs) were localized on the Chr14 and Chr17, respectively. 155 toxin proteins were predicted from mandarin fish genome. In summary, the high-quality genome assembly of the mandarin fish provides novel insights into the feeding habit of live prey and offers a valuable genetic resource for the quality improvement of this freshwater fish.

Recruitment of toxin-like proteins with ancestral venom function supports endoparasitic lifestyles of Myxozoa

Cnidarians are the oldest lineage of venomous animals and use nematocysts to discharge toxins. Whether venom toxins have been recruited to support parasitic lifestyles in the Endocnidozoa (Myxozoa + Polypodium) is, however, unknown. To examine this issue we variously employed transcriptomic, proteomic, associated molecular phylogenies, and localisation studies on representative primitive and derived myxozoans (Malacosporea and Myxosporea, respectively), Polypodium hydriforme, and the free-living staurozoan Calvadosia cruxmelitensis. Our transcriptomics and proteomics analyses provide evidence for expression and translation of venom toxin homologs in myxozoans. Phylogenetic placement of Kunitz type serine protease inhibitors and phospholipase A2 enzymes reveals modification of toxins inherited from ancestral free-living cnidarian toxins, and that venom diversity is reduced in myxozoans concordant with their reduced genome sizes. Various phylogenetic analyses of the Kunitz-type toxin family in Endocnidozoa suggested lineage-specific gene duplications, which offers a possible mechanism for enhancing toxin diversification. Toxin localisation in the malacosporean Buddenbrockia plumatellae substantiates toxin translation and thus illustrates a repurposing of toxin function for endoparasite development and interactions with hosts, rather than for prey capture or defence. Whether myxozoan venom candidates are expressed in transmission stages (e.g. in nematocysts or secretory vesicles) requires further investigation.

Phylogenetic and Selection Analysis of an Expanded Family of Putatively Pore-Forming Jellyfish Toxins (Cnidaria: Medusozoa)

Many jellyfish species are known to cause a painful sting, but box jellyfish (class Cubozoa) are a well-known danger to humans due to exceptionally potent venoms. Cubozoan toxicity has been attributed to the presence and abundance of cnidarian-specific pore-forming toxins called jellyfish toxins (JFTs), which are highly hemolytic and cardiotoxic. However, JFTs have also been found in other cnidarians outside of Cubozoa, and no comprehensive analysis of their phylogenetic distribution has been conducted to date. Here, we present a thorough annotation of JFTs from 147 cnidarian transcriptomes and document 111 novel putative JFTs from over 20 species within Medusozoa. Phylogenetic analyses show that JFTs form two distinct clades, which we call JFT-1 and JFT-2. JFT-1 includes all known potent cubozoan toxins, as well as hydrozoan and scyphozoan representatives, some of which were derived from medically relevant species. JFT-2 contains primarily uncharacterized JFTs. Although our analyses detected broad purifying selection across JFTs, we found that a subset of cubozoan JFT-1 sequences are influenced by gene-wide episodic positive selection compared with homologous toxins from other taxonomic groups. This suggests that duplication followed by neofunctionalization or subfunctionalization as a potential mechanism for the highly potent venom in cubozoans. Additionally, published RNA-seq data from several medusozoan species indicate that JFTs are differentially expressed, spatially and temporally, between functionally distinct tissues. Overall, our findings suggest a complex evolutionary history of JFTs involving duplication and selection that may have led to functional diversification, including variability in toxin potency and specificity.