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.

Reduced Toxicity of Centruroides vittatus (Say, 1821) May Result from Lowered Sodium β Toxin Gene Expression and Toxin Protein Production

Body tissue and venom glands from an eastern population of the scorpion Centruroides vittatus (Say, 1821) were homogenized and molecular constituents removed to characterize putative sodium β toxin gene diversity, RT-qPCR, transcriptomic, and proteomic variation. We cloned sodium β toxins from genomic DNA, conducted RT-qPCR experiments with seven sodium β toxin variants, performed venom gland tissue RNA-seq, and isolated venom proteins for mass spectrophotometry. We identified >70 putative novel sodium β toxin genes, 111 toxin gene transcripts, 24 different toxin proteins, and quantified sodium β toxin gene expression variation among individuals and between sexes. Our analyses contribute to the growing evidence that venom toxicity among scorpion taxa and their populations may be associated with toxin gene diversity, specific toxin transcripts variation, and subsequent protein production. Here, slight transcript variation among toxin gene variants may contribute to the major toxin protein variation in individual scorpion venom composition.

Expression of Brown and Southern Black Widow Spider (Araneae: Theridiidae) Latrotoxins Is Tissue- and Life Stage-Specific for α-Latroinsectotoxins and δ-Latroinsectotoxins and Is Ubiquitous for α-Latrotoxins

Widow spiders are widely known for their potent venom toxins that make them among the few spiders of medical concern. The latrotoxins are the most well-studied widow toxins and include both the vertebrate-specific latrotoxins and the insect-specific latroinsectotoxins (LITs). Previous studies have shown that toxins are not limited to expression in the venom glands of adult spiders; however, gaps exist in latrotoxin screening across all life stages for brown widows, Latrodectus geometricus and southern black widows, Latrodectus mactans. In this study, we screened male and female venom gland, cephalothorax, and abdomen tissues, spiderling cephalothorax and abdomen tissues, and eggs of both L. geometricus and L. mactans, for the presence of three latrotoxins: α-latrotoxin (α-LTX), and α- and δ-latroinsectotoxins (α/δ-LITs). Widows were locally collected. Extracted RNA was used to prepare cDNA that was analyzed by PCR for the presence or absence of latrotoxin expression. Results show that expression profiles between the two species are very similar but not identical. Expression of α-LTX was found in all life stages in all tissues examined for both species. For both species, no LIT expression was detected in eggs and variable patterns of α-LIT expression were detected in spiderlings and adults. Notably, δ-LIT could only be detected in females for both species. Our results show that latrotoxin expression profiles differ within and between widow species. Data on their expression distribution provide further insight into the specific latrotoxins that contribute to toxicity profiles for each life stage in each species and their specific role in widow biology.

Hierarchical chromatin features reveal the toxin production in Bungarus multicinctus

Background Bungarus multicinctus, from which a classical Chinese medicine is produced, is known as the most venomous land snake in the world, but the chromatin organization and transcription factor activity during venom replenishment progress have not been explored yet. This study aimed to determine the roles of chromatin structure in toxin activity via bioinformatics and experimental validation. Methods Chromosome conformation capture (Hi-C) analysis was used to examine interactions among chromosomes and identify different scales of chromatin during envenomation in B. multicinctus. Correlations between epigenetic modifications and chromatin structure were verified through ChIP-seq analysis. RNA-seq was used to validate the influence of variation in chromatin structure and gene expression levels on venom production and regulation. Results Our results suggested that intra-chromosomal interactions are more intense than inter-chromosomal interactions among the control group, 3-day group of venom glands and muscles. Through this, we found that compartmental transition was correlated with chromatin interactions. Interestingly, the up-regulated genes in more compartmental switch regions reflect the function of toxin activity. Topologically associated domain (TAD) boundaries enriched with histone modifications are associated with different distributions of genes and the expression levels. Toxin-coding genes in the same loop are highly expressed, implying that the importance of epigenetic regulation during envenomination. On a smaller scale, the epigenetic markers affect transcriptional regulation by controlling the recruitment/inhibition of transcription initiation complexes. Conclusions Chromatin structure and epigenetic modifications could play a vital status role in the mechanisms of venom regulation in B. multicinctus.

Pectoral-fin glands and delivery apparatus in the catfish genus Brachyrhamdia Myers, 1927 (Siluriformes: Heptapteridae)

The Siluriformes, popularly known as catfishes, are probably the vertebrate group with the highest diversity of venomous animals, even though only approximately a hundred venomous catfishes are reported to date. Venomous catfishes might present a delivery system apparatus, formed by an unbranched ray at the leading edge of pectoral and dorsal fins (spine), which can be stiffened and pungent, while venom glands can be present at the surface of such spines and/or the axillary region. This work investigated the presence, morphology and distribution of glands and pectoral-fin delivery apparatus in the heptapterid Brachyrhamdia genus. Pectoral-fin spine external morphology was compared across all valid species in the genus, histological sections of the pectoral-fin spine and axillary regions of B. heteropleura and B. marthae were produced, and dissections of the pectoral girdle region of the mentioned species were analyzed. The histological sections confirmed the presence of pectoral-fin glands at the surface of the pectoral-fin spine of Brachyrhamdia species, and cellular morphology indicates these glands are probably venomous. Also, we found a piriform gland at the axillary region, whose cell morphology is like the reported for other catfishes. However, we cannot currently confirm or deny axillary gland participation in the venom delivery apparatus. This work constitutes the first report of venom glands in Brachyrhamdia, and the first description of Heptapteridae axillary glands.

Origin and Characterization of Extracellular Vesicles Present in the Spider Venom of Ornithoctonus hainana

Extracellular vesicles (EVs), including exosomes and microvesicles, are membranous vesicles released from nearly all cellular types. They contain various bioactive molecules, and their molecular composition varies depending on their cellular origin. As research into venomous animals has progressed, EVs have been discovered in the venom of snakes and parasitic wasps. Although vesicle secretion in spider venom glands has been observed, these secretory vesicles’ origin and biological properties are unknown. In this study, the origin of the EVs from Ornithoctonus hainana venom was observed using transmission electron microscopy (TEM). The Ornithoctonus hainana venom extracellular vesicles (HN-EVs) were isolated and purified by density gradient centrifugation. HN-EVs possess classic membranous vesicles with a size distribution ranging from 50 to 150 nm and express the arthropod EV marker Tsp29Fb. The LC-MS/MS analysis identified a total of 150 proteins, which were divided into three groups according to their potential function: conservative vesicle transport-related proteins, virulence-related proteins, and other proteins of unknown function. Functionally, HN-EVs have hyaluronidase activity and inhibit the proliferation of human umbilical vein endothelial cells (HUVECs) by affecting the cytoskeleton and cell cycle. Overall, this study investigates the biological characteristics of HN-EVs for the first time and sheds new light on the envenomation process of spider venom.

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.

Genome of the pincer wasp Gonatopus flavifemur reveals unique venom evolution and a dual adaptation to parasitism and predation

Background Hymenoptera comprise extremely diverse insect species with extensive variation in their life histories. The Dryinidae, a family of solitary wasps of Hymenoptera, have evolved innovations that allow them to hunt using venom and a pair of chelae developed from the fore legs that can grasp prey. Dryinidae larvae are also parasitoids of Auchenorrhyncha, a group including common pests such as planthoppers and leafhoppers. Both of these traits make them effective and valuable for pest control, but little is yet known about the genetic basis of its dual adaptation to parasitism and predation. Results We sequenced and assembled a high-quality genome of the dryinid wasp Gonatopus flavifemur, which at 636.5 Mb is larger than most hymenopterans. The expansion of transposable elements, especially DNA transposons, is a major contributor to the genome size enlargement. Our genome-wide screens reveal a number of positively selected genes and rapidly evolving proteins involved in energy production and motor activity, which may contribute to the predatory adaptation of dryinid wasp. We further show that three female-biased, reproductive-associated yellow genes, in response to the prey feeding behavior, are significantly elevated in adult females, which may facilitate the egg production. Venom is a powerful weapon for dryinid wasp during parasitism and predation. We therefore analyze the transcriptomes of venom glands and describe specific expansions in venom Idgf-like genes and neprilysin-like genes. Furthermore, we find the LWS2-opsin gene is exclusively expressed in male G. flavifemur, which may contribute to partner searching and mating. Conclusions Our results provide new insights into the genome evolution, predatory adaptation, venom evolution, and sex-biased genes in G. flavifemur, and present genomic resources for future in-depth comparative analyses of hymenopterans that may benefit pest control.

Convergent evolution of venom gland transcriptomes across Metazoa

Animals have repeatedly evolved specialized organs and anatomical structures to produce and deliver a cocktail 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 the first 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 turns, 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 the first step towards an understanding of the molecular mechanisms underlying the repeated evolution of one of the most successful adaptive traits in the animal kingdom.