Defensive Toxins Research Articles

The effects of diatom polyunsaturated aldehydes on embryonic and larval zebrafish (Danio rerio).

Marine diatoms are pervasive in many planktonic and benthic environments and represent an important food source for a wide range of species. Some diatoms produce polyunsaturated aldehydes (PUAs) as defensive toxins. PUA exposure is known to reduce the fecundity of invertebrate grazers like copepods and echinoderm larvae, but little is known about the effects of PUAs on vertebrates. Many fish species are likely to come into close contact with diatoms. Many may deposit eggs on diatom-coated substrates, consume diatoms as larvae, and/or feed heavily on zooplankters that may be gut-loaded with diatoms. The purpose of this study was to test whether dissolved diatom PUAs affect the early life stages of a model fish, Danio rerio (zebrafish). To test this, zebrafish embryos and larvae were exposed to proportionally increasing mixtures of the dissolved diatom PUAs 2E,4E-decadienal, 2E,4E-octadienal and 2E,4E-heptadienal. Under PUA exposure, three metrics of fitness were assessed: embryo heart rate, larval size at hatch, and pre-feeding mortality rate. Zebrafish embryos exposed at 24 h post fertilization (hpf) experienced decreased average heart rate after 2 days of PUA exposure. Embryos 24 hpf exposed to PUA mixtures for 6 days showed a reduction in size in comparison to embryos from controls. Embryos exposed to PUAs from 2 hpf until death showed lower survivorship compared to larvae in controls. The results of this study suggest that larval fish that are sympatric with PUA producing diatoms during their embryonic and larval stages may be susceptible to detrimental effects from PUA exposure.

Characterizing the evolution of defense in a tripartite marine symbiosis using adaptive dynamics

Abstract The evolution and maintenance of symbiotic systems remains a fascinating puzzle. While the coevolutionary dynamics of bipartite (host–symbiont) systems are well-studied, the dynamics of more complex systems have only recently garnered attention with increasing technological advances. We model a tripartite system inspired by the marine symbiotic relationship between the alga Bryopsis sp., its intracellular defensive bacterial symbiont “Candidatus Endobryopsis kahalalidifaciens,” which produces a toxin that protects the alga against fish herbivores, and the sea-slug Elysia rufescens (Zan et al., 2019), which is not deterred by the toxin. We disentangle the role of selection on different actors within this system by investigating evolutionary scenarios where defense evolves as (i) a host-controlled trait that reduces algal reproductive ability; (ii) a symbiont-controlled trait that impacts symbiont transmission; and (iii) a trait jointly controlled by both host and symbiont. Optimal investment in defensive toxins varies based on the characteristics of the host, symbiont, and sea slug; and evolutionary trajectories are modulated by trade-off shape, i.e., a strongly decelerating trade-off between defense and symbiont transmission can drive symbiont diversification via evolutionary branching. Increasing slug herbivory reduces host investment in defense to favor reproduction, while symbiont investment in defense first declines and then increases as host density declines to the degree that horizontal symbiont transmission is no longer beneficial. Increasing vertical transmission selects for reduced defense by the host when it evolves as a jointly controlled trait, as a result of investment by the symbiont. Our theoretical exploration of the evolution of defensive symbiosis in scenarios involving interactions with multiple herbivores provides a first window into the origin and maintenance of the Bryopsis sp. system, and adds another piece to the puzzle of the evolution of symbiotic systems.

Predator-induced defense decreases growth rate and photoprotective capacity in a nitrogen-limited dinoflagellate, Alexandrium minutum

Some dinoflagellates produce toxic secondary metabolites that correlate with increased resistance to grazers. The allocation costs of toxin production have been repeatedly addressed, but with conflicting results. Few studies have considered the potential costs of this defense to the photosystem, even though defense toxins (e.g., karlotoxins and brevetoxins) are closely linked to the photoprotective process. Here, we used chemical cues from copepods to induce paralytic shellfish toxin (PST) production in resource-limited Alexandrium minutum and quantitatively determined the growth rate and potential trade-offs with the photosystem process. The results show that grazer-induced, more toxic A. minutum had larger cell volume, lower cell division rate, and lower pigment content under nitrogen-limited conditions than control cells. In addition, predator cues led to a lower relative abundance of photoprotective xanthophylls and a reduced de-epoxidation efficiency of the xanthophyll cycle under high light conditions, reducing the ability of the cells to resist photodamage. Decreased photoprotective capacity may reflect an overlooked defense cost of toxin production.

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.

Deep divergences among inconspicuously colored clades of Epipedobates poison frogs

Poison frogs (Dendrobatidae) are famous for their aposematic species, having a combination of diverse color patterns and defensive skin toxins, yet most species in this family are inconspicuously colored and considered non-aposematic. Epipedobates is among the youngest genus-level clades of Dendrobatidae that includes both aposematic and inconspicuous species. Using Sanger-sequenced mitochondrial and nuclear markers, we demonstrate deep genetic divergences among inconspicuous species of Epipedobates but relatively shallow genetic divergences among conspicuous species. Our phylogenetic analysis includes broad geographic sampling of the inconspicuous lineages typically identified as E. boulengeri and E. espinosai, which reveals two putative new species, one in west-central Colombia (E. sp. 1) and the other in north-central Ecuador (E. aff. espinosai). We conclude that E. darwinwallacei is a junior subjective synonym of E. espinosai. We also clarify the geographic distributions of inconspicuous Epipedobates species including the widespread E. boulengeri. We provide a qualitative assessment of the phenotypic diversity in each nominal species, with a focus on the color and pattern of inconspicuous species. We conclude that Epipedobates contains eight known valid species, six of which are inconspicuous. A relaxed molecular clock analysis suggests that the most recent common ancestor of Epipedobates is ∼11.1 million years old, which nearly doubles previous estimates. Last, genetic information points to a center of species diversity in the Chocó at the southwestern border of Colombia with Ecuador. A Spanish translation of this text is available in the supplementary materials.

The first chromosome-level genome assembly and transcriptome sequencing provide insights into cantharidin production of the blister beetles.

Blister beetles (Coleoptera: Meloidae) produce a natural defensive toxin cantharidin (CTD), which has been used for various cancer treatments and other diseases. Currently, the lack of chromosome-level reference genomes in Meloidae limits further understanding of the mechanism of CTD biosynthesis and environmental adaptation. In this study, the chromosome-level genome assembly of Mylabris phalerata was generated based on PacBio and Hi-C sequencing. This reference genome was about 136.68 Mb in size with contig N50 of 9.17 Mb and composed of 12 chromosomes. In comparison to six other Coleoptera insects, M. phalerata exhibited multiple expanded gene families enriched in juvenile hormone (JH) biosynthetic process pathway, farnesol dehydrogenase activity, and cytochrome P450, which may be related to CTD biosynthesis. Consistently, the transcriptomic analysis suggested the "terpenoid backbone biosynthesis" pathway and "the juvenile hormone" as putative core pathways of CTD biosynthesis and presented eight up-regulated differential expression genes in male adults as candidate genes. It is possible that the restricted feeding niche and lifestyle of M. phalerata were the cause of the gene family's contraction of odorant binding proteins. The ABC transporters (ABCs) related to exporting bound toxins out of the cell and the resistance to the self-secreted toxins (e.g. CTD) were also contracted, possibly due to other self-protection strategies in M. phalerata. A foundation of understanding CTD biosynthesis and environmental adaptation of blister beetles will be established by our reference genome and discoveries.

The pederin-producing bacteria density dynamics in Paederus fuscipes at different developmental stages.

Pederin, a defensive toxin in Paederus fuscipes, is produced by an uncultured Gram-negative symbiont, which establishes a stable symbiotic relationship with a female host before completion of metamorphosis. However, the transmission process of pederin-producing bacteria (PPB) in P. fuscipes at different life stages remains unknown. Herein, the PPB population dynamics and transcriptome atlas for P. fuscipes development (egg, first-instar larva, second-instar larva, pupa, and newly emerged female and male) were characterised. We found that a microbial layer containing PPB covered the eggshell, which could be sterilised by smearing the eggshell with streptomycin. Maternal secretions over the eggshell are likely the main PPB acquisition route for P. fuscipes offspring. The PPB density in eggs was significantly higher than that in other life stages (p < 0.05), which demonstrated that the beetle mothers gave more PPB than the larvae acquired. Physiological changes (hatching and eclosion) led to a decreased PPB density in P. fuscipes. Pattern recognition receptors related to Gram-negative bacteria recognition were identified from P. fuscipes transcriptomes across various life stages, which might be used to screen genes involved in PPB regulation. These results will help advance future efforts to determine the molecular mechanisms of PPB colonisation of P. fuscipes.

Investigation of sex expression profiles and the cantharidin biosynthesis genes in two blister beetles.

Cantharidin (CTD) is a well-established defensive toxin synthesized by blister beetles, displaying both therapeutic potential and toxicity. Among these beetles, Hycleus cichorii and Hycleus phaleratus are the two most commercially significant species due to their capacity to produce CTD in males. In this investigation, we conducted a gene expression profiling analysis of male and female individuals of these two species, utilizing the Illumina Hiseq4000 platform. We identified 7,983 expressed genes, including 2,823 differentially expressed genes (DEGs) shared by both male and female blister beetles. Nineteen genes related to CTD biosynthesis in the terpenoid backbone biosynthesis pathway were identified, including hydroxymethylglutaryl-CoA reductase (HMGR; EC:1.1.1.34), which demonstrated a significant correlation with CTD content. Furthermore, hydroxymethylglutaryl-CoA synthase (HMGS; EC:2.3.3.10) and isopentenyl-diphosphate Delta-isomerase (IDI; EC:5.3.3.2) were also found to be significantly up-regulated in males. Comparative analysis revealed that NADP+-dependent farnesol dehydrogenase (FOHSDR; EC:1.1.1.216) and farnesyl diphosphate synthase (FDPS; EC:2.5.1.1) had the highest copy number in these beetles, significantly higher than the copy number of the other four non-Meloidae insects. The analysis of the protein-protein interaction network of genes related to CTD biosynthesis revealed that the acetyl-CoA C-acetyltransferase (ACAT; EC:2.3.1.9) gene was the central gene, exhibiting greater expression in male blister beetles than in females. This study offers novel insights into the mechanisms of CTD biosynthesis in blister beetles and enhances our comprehensions of the association between particular genes and CTD content.

A nutrition-defence trade-off drives diet choice in a toxic plant generalist.

Plant toxicity shapes the dietary choices of herbivores. Especially when herbivores sequester plant toxins, they may experience a trade-off between gaining protection from natural enemies and avoiding toxicity. The availability of toxins for sequestration may additionally trade off with the nutritional quality of a potential food source for sequestering herbivores. We hypothesized that diet mixing might allow a sequestering herbivore to balance nutrition and defence (via sequestration of plant toxins). Accordingly, here we address diet mixing and sequestration of large milkweed bugs (Oncopeltus fasciatus) when they have differential access to toxins (cardenolides) in their diet. In the absence of toxins from a preferred food (milkweed seeds), large milkweed bugs fed on nutritionally adequate non-toxic seeds, but supplemented their diet by feeding on nutritionally poor, but cardenolide-rich milkweed leaf and stem tissues. This dietary shift corresponded to reduced insect growth but facilitated sequestration of defensive toxins. Plant production of cardenolides was also substantially induced by bug feeding on leaf and stem tissues, perhaps benefitting this cardenolide-resistant herbivore. Thus, sequestration appears to drive diet mixing in this toxic plant generalist, even at the cost of feeding on nutritionally poor plant tissue.

The plant toxin 4-methylsulfinylbutyl isothiocyanate decreases herbivore performance and modulates cellular and humoral immunity.

Insect herbivores frequently encounter plant defense molecules, but the physiological and ecological consequences for their immune systems are not fully understood. The majority of studies attempting to relate levels of plant defensive chemistry to herbivore immune responses have used natural population or species-level variation in plant defensive chemistry. Yet, this potentially confounds the effects of plant defense chemistry with other potential plant trait differences that may affect the expression of herbivore immunity. We used an artificial diet containing known quantities of a plant toxin (4-methylsulfinylbutyl isothiocyanate; 4MSOB-ITC or ITC, a breakdown product of the glucosinolate glucoraphanin upon herbivory) to explicitly explore the effects of a plant toxin on the cellular and humoral immune responses of the generalist herbivore Trichoplusia ni (Lepidoptera: Noctuidae) that frequently feeds on glucosinolate-containing plants. Caterpillars feeding on diets with high concentrations of ITC experienced reduced survivorship and growth rates. High concentrations of ITC suppressed the appearance of several types of hemocytes and melanization activity, which are critical defenses against parasitic Hymenoptera and microbial pathogens. In terms of T. ni humoral immunity, only the antimicrobial peptide (AMP) genes lebocin and gallerimycin were significantly upregulated in caterpillars fed on diets containing high levels of ITC relative to caterpillars that were provided with ITC-free diet. Surprisingly, challenging caterpillars with a non-pathogenic strain of Escherichia coli resulted in the upregulation of the AMP gene cecropin. Feeding on high concentrations of plant toxins hindered caterpillar development, decreased cellular immunity, but conferred mixed effects on humoral immunity. Our findings provide novel insights into the effects of herbivore diet composition on insect performance demonstrating the role of specific plant defense toxins that shape herbivore immunity and trophic interactions.

Coevolutionary escalation led to differentially adapted paralogs of an insect's Na,K-ATPase optimizing resistance to host plant toxins.

Cardiac glycosides are chemical defence toxins known to fatally inhibit the Na,K-ATPase (NKA) throughout the animal kingdom. Several animals, however, have evolved target-site insensitivity through substitutions in the otherwise highly conserved cardiac glycoside binding pocket of the NKA. The large milkweed bug, Oncopeltus fasciatus, shares a long evolutionary history with cardiac glycoside containing plants that led to intricate adaptations. Most strikingly, several duplications of the bugs' NKA1α gene provided the opportunity for differential resistance-conferring substitutions and subsequent sub-functionalization of the enzymes. Here, we analysed cardiac glycoside resistance and ion pumping activity of nine functional NKA α/β-combinations of O. fasciatus expressed in cell culture. We tested the enzymes with two structurally distinct cardiac glycosides, calotropin, a host plant compound, and ouabain, a standard cardiac glycoside. The identity and number of known resistance-conferring substitutions in the cardiac glycoside binding site significantly impacted activity and toxin resistance in the three α-subunits. The β-subunits also influenced the enzymes' characteristics, yet to a lesser extent. Enzymes containing the more ancient αC-subunit were inhibited by both compounds but much more strongly by the host plant toxin calotropin than by ouabain. The sensitivity to calotropin was diminished in enzymes containing the more derived αB and αA, which were only marginally inhibited by both cardiac glycosides. This trend culminated in αAβ1 having higher resistance against calotropin than against ouabain. These results support the coevolutionary escalation of plant defences and herbivore tolerance mechanisms. The possession of multiple paralogs additionally mitigates pleiotropic effects by compromising between ion pumping activity and resistance.

Intrinsic Factors Associated with Dietary Toxin Quantity and Concentration in the Nuchal Glands of a Natricine Snake Rhabdophis Tigrinus.

The snake Rhabdophis tigrinus sequesters cardiotonic steroids, bufadienolides (BDs), from ingested toads and stores them in the nuchal glands as defensive toxins. It has previously been shown that there are individual differences in the total quantity of BDs stored in the nuchal glands of adult R. tigrinus and that BD quantities and profiles of R. tigrinus exhibit geographic variation. However, no previous study has examined the total quantity of BDs as a percentage of body mass (relative BD quantity) and the concentration of BDs in the nuchal gland fluid (BD gland concentration). In addition, intrinsic factors that are associated with relative BD quantity and BD concentration have not been examined within a single population. We collected 158 adult snakes from an area of central Japan from May to October and analyzed their BD quantities by UV analysis. We assessed individual differences in BD quantity, relative BD quantity and BD gland concentration. We found that 1) in approximately 60% of the 158 individuals, the BD gland concentration was greater than 50%; 2) body length and body condition are positively correlated with relative BD quantity and BD gland concentration; 3) even in a single population, individual differences of BD quantity are large, and are greater in females than in males; and 4) relative BD quantity and BD gland concentration of females during the gestation season are lower than those during the non-gestation season.

Does the Glucocorticoid Stress Response Make Toads More Toxic? An Experimental Study on the Regulation of Bufadienolide Toxin Synthesis.

Chemical defense is a crucial component of fitness in many organisms, yet the physiological regulation of defensive toxin synthesis is poorly understood, especially in vertebrates. Bufadienolides, the main defensive compounds of toads, are toxic to many predators and other natural enemies, and their synthesis can be upregulated by stressors, including predation risk, high conspecific density, and pollutants. Thus, higher toxin content may be the consequence of a general endocrine stress response in toads. Therefore, we hypothesized that bufadienolide synthesis may be stimulated by elevated levels of corticosterone (CORT), the main glucocorticoid hormone of amphibians, or by upstream regulators that stimulate CORT production. To test these alternatives, we treated common toad tadpoles with exogenous CORT (exoCORT) or metyrapone (MTP, a CORT-synthesis inhibitor that stimulates upstream regulators of CORT by negative feedback) in the presence or absence of predation cues for 2 or 6 days, and subsequently measured their CORT release rates and bufadienolide content. We found that CORT release rates were elevated by exoCORT, and to a lesser extent also by MTP, regardless of treatment length. Bufadienolide content was significantly decreased by treatment with exoCORT for 6 days but was unaffected by exposure to exoCORT for 2 days or to MTP for either 6 or 2 days. The presence or absence of predation cues affected neither CORT release rate nor bufadienolide content. Our results suggest that changes in bufadienolide synthesis in response to environmental challenges are not driven by CORT but may rather be regulated by upstream hormones of the stress response.

Distribution characteristics of traditional Chinese medicine syndromes in 4 367 adult influenza patients: a Meta-analysis

To systematically evaluate the distribution characteristics of traditional Chinese medicine (TCM) syndromes in adult influenza patients and to provide a basis for the TCM syndrome differentiation of influenza. The CNKI, CBM, Wanfang, VIP, PubMed, Embase, Cochrane Library databases were searched to collect cross-sectional studies on the distribution pattern of TCM syndromes in adult patients with influenza. The risk of bias assessment tool for cross-sectional studies developed by the Joanna Briggs Institute (JBI) evidence-based health care center was used to evaluate the literature quality, and the Stata 15.1 software was used to conduct a Meta-analysis of the pooled effect sizes of the included studies. A total of 11 studies with 4 367 influenza patients were included. Quality assessment results of JBI showed that the risk bias was higher in the sample size calculation, and the description of sampling modalities and response rate was unclear. There were 17 influenza syndromes after specification, and a single group rate Meta-analysis was performed of the syndromes with ≥ 50 incident cases showed that there were 9 syndromes with an incidence ≥ 10% and statistical significance, the top 5 syndromes were syndrome of wind and heat invading the defense [n = 1 583, RATE = 34.3%, 95% confidence interval (95%CI) was 22.2%-46.3%], syndrome of exterior cold and interior heat (n = 1 122, RATE = 36.1%, 95%CI was 21.2%-51.1%), syndrome of wind-cold fettering the exterior (n = 860, RATE = 19.4%, 95%CI was 10.7%-28.0%), syndrome of heat and toxin in the lung (n = 217, RATE = 17.1%, 95%CI was 9.1%-25.0%), and syndrome of disease involving both defense phase and qi phase (n = 184, RATE = 38.8%, 95%CI was 14.2%-63.5%). The results of the subgroup analysis in different geographical regions showed that the frequency of distribution of syndrome of wind and heat invading the defense and heat and toxin in the lung was higher in the South (RATE: 36.5%, 18.6%) than in the North (RATE: 30.9%, 15.4%), and the frequency of distribution of syndrome of wind-cold fettering the exterior and exterior cold and interior heat in the North (RATE: 23.8%, 40.1%) was higher than that in the South (RATE: 15.7%, 32.3%). There are 9 common TCM syndromes of influenza, including wind and heat invading the defense syndrome, exterior cold and interior heat syndrome, wind-cold fettering the exterior syndrome, heat and toxin in the lung syndrome, disease involving both defense phase and qi phase syndrome, wind and heat complicated by dampnessinvading the surface syndrome, wind and cold complicated by dampnessinvading the surface syndrome, defense phase syndrome and dampness and heatinvading the surface syndrome, which can provide a reference for the TCM syndrome differentiation and treatment of influenza.

The Detoxification Enzymatic Responses of Plutella xylostella (Lepidoptera: Plutellidae) to Cantharidin

Plutella xylostella (L.) (Lepidoptera: Plutellidae) is one of the most destructive pests of Brassicaceae vegetables. Cantharidin is an insect-derived defensive toxin, which has been reported to have toxicity to a variety of pests and especially lepidopteran pests. Although the toxicity of cantharidin on P. xylostella has been demonstrated, there is little information available on the specific detoxification response of P. xylostella against cantharidin. This study investigates the enzymatic response (including serine/threonine phosphatases [PSPs], carboxylesterases [CarEs], glutathione-S-transferases [GSTs], and cytochrome P450 monooxygenases [P450]) in P. xylostella to the sublethal and low lethal concentrations of cantharidin (LC10 and LC25). Results showed that the inhibitory activity of PSPs was increased and then decreased in vivo, while PSPs activity could be almost completely inhibited in vitro. Interestingly, the activities of detoxification enzymes (GST, CarE, and P450) in P. xylostella displayed a trend of decreasing and then increasing after exposure to the two concentrations of cantharidin. Notably, the increase in P450 enzyme activity was the most significant. The increasing trend of detoxification enzyme activity was congruent with the recovery trend of PSPs activity. This study contributes to our understanding of the detoxification mechanism of cantharidin in P. xylostella and helps in the further development of biogenic agents.

Surveillance of venomous arthropods on a Nigerian Public University campus

Background: Venomous arthropods are arthropod species that use toxins for defense and offense. The objective of this study was to assess the biodiversity and abundance of venomous arthropods on the premises of Nnamdi Azikiwe University (NAU), Awka, Nigeria. Materials and Methods: Pitfall traps, sticky, and bait traps, beating and knockdown, and active search and handpicking with forceps were the methods used to collect various arthropods on the university campus. The data collected during this study were analyzed using the statistical package for social sciences version 23.0. Results: A total of 1,070 venomous arthropods were collected from the study areas, and they belong to four classes, six orders, fourteen families, and twenty-two species of the phylum Arthropoda. Conclusion: This study demonstrates that venomous arthropods are widespread on the NAU, Awka campus. Regular fumigation of the university premises is important in preventing human injuries, and infectious diseases carried by some arthropods.

Physiological constraints dictate toxin spatial heterogeneity in snake venom glands

BackgroundVenoms are ecological innovations that have evolved numerous times, on each occasion accompanied by the co-evolution of specialised morphological and behavioural characters for venom production and delivery. The close evolutionary interdependence between these characters is exemplified by animals that control the composition of their secreted venom. This ability depends in part on the production of different toxins in different locations of the venom gland, which was recently documented in venomous snakes. Here, we test the hypothesis that the distinct spatial distributions of toxins in snake venom glands are an adaptation that enables the secretion of venoms with distinct ecological functions.ResultsWe show that the main defensive and predatory peptide toxins are produced in distinct regions of the venom glands of the black-necked spitting cobra (Naja nigricollis), but these distributions likely reflect developmental effects. Indeed, we detected no significant differences in venom collected via defensive ‘spitting’ or predatory ‘biting’ events from the same specimens representing multiple lineages of spitting cobra. We also found the same spatial distribution of toxins in a non-spitting cobra and show that heterogeneous toxin distribution is a feature shared with a viper with primarily predatory venom.ConclusionsOur findings suggest that heterogeneous distributions of toxins are not an adaptation to controlling venom composition in snakes. Instead, it likely reflects physiological constraints on toxin production by the venom glands, opening avenues for future research on the mechanisms of functional differentiation of populations of protein-secreting cells within adaptive contexts.

Evolution, Expression Patterns, and Distribution of Novel Ribbon Worm Predatory and Defensive Toxins.

Ribbon worms are active predators that use an eversible proboscis to inject venom into their prey and defend themselves with toxic epidermal secretions. Previous work on nemertean venom has largely focused on just a few species and has not investigated the different predatory and defensive secretions in detail. Consequently, our understanding of the composition and evolution of ribbon worm venoms is still very limited. Here, we present a comparative study of nemertean venom combining RNA-seq differential gene expression analyses of venom-producing tissues, tandem mass spectrometry-based proteomics of toxic secretions, and mass spectrometry imaging of proboscis sections, to shed light onto the composition and evolution of predatory and defensive toxic secretions in Antarctonemertes valida. Our analyses reveal a wide diversity of putative defensive and predatory toxins with tissue-specific gene expression patterns and restricted distributions to the mucus and proboscis proteomes respectively, suggesting that ribbon worms produce distinct toxin cocktails for predation and defense. Our results also highlight the presence of numerous lineage-specific toxins, indicating that venom evolution is highly divergent across nemerteans, producing toxin cocktails that might be finely tuned to subdue different prey. Our data also suggest that the hoplonemertean proboscis is a highly specialized predatory organ that seems to be involved in a variety of biological functions besides predation, including secretion and sensory perception. Overall, our results advance our knowledge into the diversity and evolution of nemertean venoms and highlight the importance of combining different types of data to characterize toxin composition in understudied venomous organisms.

Transcriptomic Analysis Reveals Potential Candidate Pathways and Genes Involved in Toxin Biosynthesis in True Toads.

Synthesized chemical defenses have broadly evolved across countless taxa and are important in shaping evolutionary and ecological interactions within ecosystems. However, the underlying genomic mechanisms by which these organisms synthesize and utilize their toxins are relatively unknown. Herein, we use comparative transcriptomics to uncover potential toxin synthesizing genes and pathways, as well as interspecific patterns of toxin synthesizing genes across 10 species of North American true toads (Bufonidae). Upon assembly and annotation of the 10 transcriptomes, we explored patterns of relative gene expression and possible protein-protein interactions across the species to determine what genes and/or pathways may be responsible for toxin synthesis. We also tested our transcriptome dataset for signatures of positive selection to reveal how selection may be acting upon potential toxin producing genes. We assembled high-quality transcriptomes of the bufonid parotoid gland, a tissue not often investigated in other bufonid-related RNAseq studies. We found several genes involved in metabolic and biosynthetic pathways (e.g., steroid biosynthesis, terpenoid backbone biosynthesis, isoquinoline biosynthesis, and glucosinolate biosynthesis) that were functionally enriched and/or relatively expressed across the 10 focal species that may be involved in the synthesis of alkaloid and steroid toxins, as well as other small metabolic compounds that cause distastefulness in bufonids. We hope that our study lays a foundation for future studies to explore the genomic underpinnings and specific pathways of toxin synthesis in toads, as well as at the macroevolutionary scale across numerous taxa that produce their own defensive toxins.

Choosy cannibals: Targeted consumption of conspecific hatchlings by larval cane toads is triggered by species-specific defensive toxins.

In many species, cannibalism is uncommon and involves nonselective consumption of conspecifics as well as heterospecifics. However, within their invasive Australian range, cane toad larvae (Rhinella marina) specifically target and voraciously consume the eggs and hatchlings of conspecifics, often extirpating entire clutches. In contrast, toad larvae rarely consume the eggs and hatchlings of native frogs. Here, we use laboratory studies to demonstrate that this selective consumption is triggered by species‐specific chemical cues: maternally‐invested bufadienolide toxins that otherwise defend cane toad eggs and hatchlings against predators. We find that these cues stimulate feeding behaviors in toad tadpoles, such that the addition of bufadienolide toxins to the water column increases predation on eggs, not only of conspecifics, but also of native anuran species that are otherwise usually ignored. In contrast, we find that cannibalism rates on conspecific hatchlings are high and unaffected by the addition of bufadienolide cues. The maternally‐invested toxins present in conspecific eggs may therefore be more easily detected post‐hatching, at which point tadpole feeding behaviors are induced whether or not additional toxin cues are present. As bufadienolide cues have previously been found to attract toad tadpoles to vulnerable hatchlings, our present findings demonstrate that the same toxin cues that attract cannibalistic tadpoles also induce them to feed, thereby facilitating cannibalism through multiple behavioral effects. Because native fauna do not produce bufadienolide toxins, the species specificity of these chemical cues in the Australian landscape may have facilitated the evolution of targeted (species‐specific) cannibalism in invasive cane toad populations. Thus, these bufadienolide toxins confer cost (increased vulnerability to cannibalism in early life‐stages) as well as benefit (reduced vulnerability to predation by other taxa).