Venom gland transcriptome analyses of two freshwater stingrays (Myliobatiformes: Potamotrygonidae) from Brazil

Bioactive Properties of Venoms Isolated from Whiptail Stingrays and the Search for Molecular Mechanisms and Targets

The venom-containing barb attached to their 'whip-like' tail provides stingrays a defensive mechanism for evading predators such as sharks. From human encounters, dermal stingray envenomation is characterized by intense pain often followed by tissue necrosis occurring over several days to weeks. The bioactive components in stingray venoms (SRVs) and their molecular targets and mechanisms that mediate these complex responses are not well understood. Given the utility of venom-derived proteins from other venomous species for biomedical and pharmaceutical applications, we set out to characterize the bioactivity of SRV extracts from three local species that belong to the Dasyatoidea 'whiptail' superfamily. Multiple cell-based assays were used to quantify and compare the in vitro effects of these SRVs on different cell lines. All three SRVs demonstrated concentration-dependent growth-inhibitory effects on three different human cell lines tested. In contrast, a mouse fibrosarcoma cell line was markedly resistant to all three SRVs, indicating the molecular target(s) for mediating the SRV effects are not expressed on these cells. The multifunctional SRV responses were characterized by an acute disruption of cell adhesion leading to apoptosis. These findings aim to guide future investigations of individual SRV proteins and their molecular targets for potential use in biomedical applications.

Health assessment of rice cultivated and harvested from plasma-irradiated seeds

This study provides the health effects assessment of rice cultivated from plasma-irradiated seeds. The rice (Oryza sativa L.) cultivated from seeds with plasma irradiation showed a growth improvement (slope-ratios of with plasma to without plasma were 1.066, 1.042, and 1.255 for tiller, and earing, and ripening periods, respectively) and an 4% increase in yield. The cultivated rice was used for repeated oral administrations to mice for 4-week period. Distilled water and rice cultivated from seeds without plasma irradiation were also used as control. The weights of the lung, kidney, liver, and spleen, with corresponding average values of 0.22 g, 0.72 g, 2.1 g, and 0.17 g for w/ plasma group and 0.22 g, 0.68 g, 2.16 g, and 0.14 g for w/o plasma group, respectively, showing no effect due to the administration of rice cultivated from plasma-irradiated seeds. Nutritional status, liver function, kidney function, and lipid, neutral fat profiles, and glucose metabolism have no significant difference between with and without plasma groups. These results show no obvious subacute effects were observed on rice grains cultivated and harvested from the mother plant that experienced growth improvement by plasma irradiation. This study provides a new finding that there is no apparent adverse health effect on the grains harvested from the plasma-irradiated seeds.

Detection of NO3- introduced in plasma-irradiated dry lettuce seeds using liquid chromatography-electrospray ionization quantum mass spectrometry (LC-ESI QMS)

Discharge plasma irradiates seeds with reactive oxygen and nitrogen species (RONS). However, RONS introduced in seeds by plasma irradiation have not been successfully detected thus far. This study provides experimental evidence that nitrate ion NO3- is introduced in lettuce seeds as RONS upon irradiation with atmospheric-pressure air dielectric barrier discharge plasma. Plasma irradiation for 5 min promotes seed germination. The components of the plasma-irradiated seeds were examined using electrospray ionization quantum mass spectrometry (ESI QMS), which revealed that the plasma irradiation introduced an ion with a mass of 62 m/z in detectable amounts. This ion was identified as NO3- by liquid chromatography (LC), multiple wavelength detector (MWD), and LC-ESI QMS. A one-dimensional simulation at electron temperature Te = 1 eV, electron density Ne = 1013/m3, and gas temperature Tg = 300 K indicated the introduction of NO3-, involving nitric oxide NO. NO3- is one of the most important ions that trigger signal transduction for germination when introduced in seeds. The scanning electron microscopy (SEM) images revealed that there was no change on the surface of the seeds after plasma irradiation. Plasma irradiation is an effective method of introducing NO3- in seeds in a dry process without causing damage.

Emerging cold plasma treatment and machine learning prospects for seed priming: a step towards sustainable food production

Seeds are vulnerable to physical and biological stresses during the germination process. Seed priming strategies can alleviate such stresses. Seed priming is a technique of treating and drying seeds prior to germination in order to accelerate the metabolic process of germination. Multiple benefits are offered by seed priming techniques, such as reducing fertilizer use, accelerating seed germination, and inducing systemic resistance in plants, which are both cost-effective and eco-friendly. For seed priming, cold plasma (CP)-mediated priming could be an innovative alternative to synthetic chemical treatments. CP priming is an eco-friendly, safe and economical, yet relatively less explored technique towards the development of seed priming. In this review, we discussed in detail the application of CP technology for seed priming to enhance germination, the quality of seeds, and the production of crops in a sustainable manner. Additionally, the combination treatment of CP with nanoparticle (NP) priming is also discussed. The large numbers of parameters need to be monitored and optimized during CP treatment to achieve the desired priming results. Here, we discussed a new perspective of machine learning for modeling plasma treatment parameters in agriculture for the development of synergistic protocols for different types of seed priming.

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.

Biochemical and biological characterization of the Hypanus americanus mucus: A perspective on stingray immunity and toxins

Stingrays skin secretions are largely studied due to the human envenoming medical relevance of the sting puncture that evolves to inflammatory events, including necrosis. Such toxic effects can be correlated to the biochemical composition of the sting mucus, according to the literature. Fish skin plays important biological roles, such as the control of the osmotic pressure gradient, protection against mechanical forces and microorganism infections. The mucus, on the other hand, is a rich and complex fluid, acting on swimming, nutrition and the innate immune system. The elasmobranch's epidermis is a tissue composed mainly by mucus secretory cells, and marine stingrays have already been described to present secretory glands spread throughout the body. Little is known about the biochemical composition of the stingray mucus, but recent studies have corroborated the importance of mucus in the envenomation process. Aiming to assess the mucus composition, a new non-invasive mucus collection method was developed that focused on peptides and proteins, and biological assays were performed to analyze the toxic and immune activities of the Hypanus americanus mucus. Pathophysiological characterization showed the presence of peptidases on the mucus, as well as the induction of edema and leukocyte recruitment in mice. The fractionated mucus improved phagocytosis on macrophages and showed antimicrobial activity against T. rubrumç. neoformans and C. albicans in vitro. The proteomic analyses showed the presence of immune-related proteins like actin, histones, hemoglobin, and ribosomal proteins. This protein pattern is similar to those reported for other fish mucus and stingray venoms. This is the first report depicting the Hypanus stingray mucus composition, highlighting its biochemical composition and importance for the stingray immune system and the possible role on the envenomation process.

Transcriptomic Characterization of the South American Freshwater Stingray Potamotrygon motoro Venom Apparatus

Venomous animals are found through a wide taxonomic range including cartilaginous fish such as the freshwater stingray Potamotrygon motoro occurring in South America, which can injure people and cause venom-related symptoms. Ensuring the efficacy of drug development to treat stingray injuries can be assisted by the knowledge of the venom composition. Here we performed a detailed transcriptomic characterization of the venom gland of the South American freshwater stingray Potamotrygon motoro. The transcripts retrieved showed 418 hits to venom components (comparably to 426 and 396 hits in other two Potamotrygon species), with high expression levels of hyaluronidase, cystatin and calglandulin along with hits uniquely found in P. motoro such as DELTA-alicitoxin-Pse1b, Augerpeptide hhe53 and PI-actitoxin-Aeq3a. We also identified undescribed molecules with extremely high expression values with sequence similarity to the SE-cephalotoxin and Rapunzel genes. Comparative analyses showed that despite being closely related, there may be significant variation among the venoms of freshwater stingrays, highlighting the importance of considering elicit care in handling different envenomation cases. Since hyaluronidase represents a major component of fish venom, we have performed phylogenetic and selective pressure analyses of this gene/protein across all fish with the available information. Results indicated an independent recruitment of the hyaluronidase into the stingray venom relative to that of venomous bony fish. The hyaluronidase residues were found to be mostly under negative selection, but 18 sites showed evidence of diversifying positive selection (P < 0.05). Our data provides new insight into stingray venom variation, composition, and selective pressure in hyaluronidase.

Phylotranscriptomics suggests the jawed vertebrate ancestor could generate diverse helper and regulatory T cell subsets

BACKGROUND

The cartilaginous fishes diverged from other jawed vertebrates ~ 450 million years ago (mya). Despite this key evolutionary position, the only high-quality cartilaginous fish genome available is for the elephant shark (Callorhinchus milii), a chimaera whose ancestors split from the elasmobranch lineage ~ 420 mya. Initial analysis of this resource led to proposals that key components of the cartilaginous fish adaptive immune system, most notably their array of T cell subsets, was primitive compared to mammals. This proposal is at odds with the robust, antigen-specific antibody responses reported in elasmobranchs following immunization. To explore this discrepancy, we generated a multi-tissue transcriptome for small-spotted catshark (Scyliorhinus canicula), a tractable elasmobranch model for functional studies. We searched this, and other newly available sequence datasets, for CD4+ T cell subset-defining genes, aiming to confirm the presence or absence of each subset in cartilaginous fishes.

RESULTS

We generated a new transcriptome based on a normalised, multi-tissue RNA pool, aiming to maximise representation of tissue-specific and lowly expressed genes. We utilized multiple transcriptomic datasets and assembly variants in phylogenetic reconstructions to unambiguously identify several T cell subset-specific molecules in cartilaginous fishes for the first time, including interleukins, interleukin receptors, and key transcription factors. Our results reveal the inability of standard phylogenetic reconstruction approaches to capture the site-specific evolutionary processes of fast-evolving immune genes but show that site-heterogeneous mixture models can adequately do so.

CONCLUSIONS

Our analyses reveal that cartilaginous fishes are capable of producing a range of CD4+ T cell subsets comparable to that of mammals. Further, that the key molecules required for the differentiation and functioning of these subsets existed in the jawed vertebrate ancestor. Additionally, we highlight the importance of considering phylogenetic diversity and, where possible, utilizing multiple datasets for individual species, to accurately infer gene presence or absence at higher taxonomic levels.

A transcriptomic investigation of appetite-regulation and digestive processes in giant grouper Epinephelus lanceolatus during early larval development

The giant grouper Epinephelus lanceolatus is an ecologically vulnerable species with high market demand. However, efforts to improve larval husbandry are hindered by a lack of knowledge surrounding larval developmental physiology. To address this shortfall, a transcriptomic approach was applied to larvae between 1 and 14 days post hatch (dph) to characterise the molecular ontogenesis of genes that influence appetite and digestion. Appetite regulating factors were detected from 1 dph, including neuropeptide Y, nesfatin-1, cocaine and amphetamine regulated transcript, cholecystokinin and pituitary adenylate cyclase activating peptide and the expression level of several genes changed sharply with the onset of exogenous feeding. The level of expression for proteases, chitinases, lipases and amylases typically followed one of two expression patterns, a general increase as development progressed, or an inverted U-shape with maximal expression at c. 6 dph. Similarly, the tendency among both expression patterns was for the level of expression to increase around the time of mouth-opening. There was also evidence to suggest the presence of putative isoforms for several digestion-related genes. We have provided an insight into appetite-regulation and digestive processes in groupers during early larval development and have developed a transcriptomic database that will aid future efforts to rear this species in an aquaculture setting.

Venomix: a simple bioinformatic pipeline for identifying and characterizing toxin gene candidates from transcriptomic data

The advent of next-generation sequencing has resulted in transcriptome-based approaches to investigate functionally significant biological components in a variety of non-model organism. This has resulted in the area of “venomics”: a rapidly growing field using combined transcriptomic and proteomic datasets to characterize toxin diversity in a variety of venomous taxa. Ultimately, the transcriptomic portion of these analyses follows very similar pathways after transcriptome assembly often including candidate toxin identification using BLAST, expression level screening, protein sequence alignment, gene tree reconstruction, and characterization of potential toxin function. Here we describe the Python package Venomix, which streamlines these processes using common bioinformatic tools along with ToxProt, a publicly available annotated database comprised of characterized venom proteins. In this study, we use the Venomix pipeline to characterize candidate venom diversity in four phylogenetically distinct organisms, a cone snail (Conidae; Conus sponsalis), a snake (Viperidae; Echis coloratus), an ant (Formicidae; Tetramorium bicarinatum), and a scorpion (Scorpionidae; Urodacus yaschenkoi). Data on these organisms were sampled from public databases, with each original analysis using different approaches for transcriptome assembly, toxin identification, or gene expression quantification. Venomix recovered numerically more candidate toxin transcripts for three of the four transcriptomes than the original analyses and identified new toxin candidates. In summary, we show that the Venomix package is a useful tool to identify and characterize the diversity of toxin-like transcripts derived from transcriptomic datasets. Venomix is available at: https://bitbucket.org/JasonMacrander/Venomix/.

Evaluating the Performance of De Novo Assembly Methods for Venom-Gland Transcriptomics

Venom-gland transcriptomics is a key tool in the study of the evolution, ecology, function, and pharmacology of animal venoms. In particular, gene-expression variation and coding sequences gained through transcriptomics provide key information for explaining functional venom variation over both ecological and evolutionary timescales. The accuracy and usefulness of inferences made through transcriptomics, however, is limited by the accuracy of the transcriptome assembly, which is a bioinformatic problem with several possible solutions. Several methods have been employed to assemble venom-gland transcriptomes, with the Trinity assembler being the most commonly applied among them. Although previous evidence of variation in performance among assembly software exists, particularly regarding recovery of difficult-to-assemble multigene families such as snake venom metalloproteinases, much work to date still employs a single assembly method. We evaluated the performance of several commonly used de novo assembly methods for the recovery of both nontoxin transcripts and complete, high-quality venom-gene transcripts across eleven snake and four scorpion transcriptomes. We varied k-mer sizes used by some assemblers to evaluate the impact of k-mer length on transcript recovery. We showed that the recovery of nontoxin transcripts and toxin transcripts is best accomplished through different assembly software, with SDT at smaller k-mer lengths and Trinity being best for nontoxin recovery and a combination of SeqMan NGen and a seed-and-extend approach implemented in Extender as the best means of recovering a complete set of toxin transcripts. In particular, Extender was the only means tested capable of assembling multiple isoforms of the diverse snake venom metalloproteinase family, while traditional approaches such as Trinity recovered at most one metalloproteinase transcript. Our work demonstrated that traditional metrics of assembly performance are not predictive of performance in the recovery of complete and high quality toxin genes. Instead, effective venom-gland transcriptomic studies should combine and quality-filter the results of several assemblers with varying algorithmic strategies.

De novo characterization of venom apparatus transcriptome of Pardosa pseudoannulata and analysis of its gene expression in response to Bt protein

Background

Pardosa pseudoannulata is a prevailing spider species, and has been regarded as an important bio-control agent of insect pests in farmland of China. However, the available genomic and transcriptomic databases of P. pseudoannulata and their venom are limited, which severely hampers functional genomic analysis of P. pseudoannulata. Recently high-throughput sequencing technology has been proved to be an efficient tool for profiling the transcriptome of relevant non-target organisms exposed to Bacillus thuringiensis (Bt) protein through food webs.

Results

In this study, the transcriptome of the venom apparatus was analyzed. A total of 113,358 non-redundant unigenes were yielded, among which 34,041 unigenes with complete or various length encoding regions were assigned biological function annotations and annotated with gene ontology and karyotic orthologous group terms. In addition, 3726 unigenes involved in response to stimulus and 720 unigenes associated with immune-response pathways were identified. Furthermore, we investigated transcriptomic changes in the venom apparatus using tag-based DGE technique. A total of 1724 differentially expressed genes (DEGs) were detected, while 75 and 372 DEGs were functionally annotated with KEGG pathways and GO terms, respectively. qPCR analyses were performed to verify the DEGs directly or indirectly related to immune and stress responses, including genes encoding heat shock protein, toll-like receptor, GST and NADH dehydrogenase.

Conclusion

This is the first study conducted to specifically investigate the venom apparatus of P. pseudoannulata in response to Bt protein exposure through tritrophic chain. A substantial fraction of transcript sequences was generated by high-throughput sequencing of the venom apparatus of P. pseudoannulata. Then a comparative transcriptome analysis showing a large number of candidate genes involved in immune response were identified by the tag-based DGE technology. This transcriptome dataset will provide a comprehensive sequence resource for furture molecular genetic research of the venom apparatus of P. pseudoannulata.

Electronic supplementary material

The online version of this article (10.1186/s12896-017-0392-z) contains supplementary material, which is available to authorized users.

Stingray venom activates IL-33 producing cardiomyocytes, but not mast cell, to promote acute neutrophil-mediated injury

One of the hallmarks of acute inflammation is neutrophil infiltration of tissues. We investigated molecular mechanisms implicated in acute neutrophilic inflammation induced by the venom of a freshwater stingray (Potamotrygon cf. henlei) in mice. Ray venom induced early mobilization of neutrophil in the microvasculature of cremaster mice and infiltration of the peritoneal cavity 2 hours after injury, in a dose-response manner. IL-1β, IL-6, TNF-α, and KC were produced. The neutrophilic infiltration did not occur in mice with ST2 receptor and MyD88 adapters neutralized, or in those with PI3K and p38 MAPK signaling blocked. Drastic reduction of neutrophil infiltration to peritoneal cavities was observed in ST2^−/−, TLR2/TLR4^−/−, MyD88^−/−, TRIF^−/− and IL-17A^−/− mice, and a partial reduction was observed in IL-18R^−/− mice. Mast cell Kit W(sh)/W(sh)-, AHR-, NLRP3-, ICE-, IL-1β-, P2RX7-, CD39-, IL-17RA-, and TBX21 KO mice retain the ability to induce neutrophilia in peritoneal cavity after ray venom injection. IL-6 and TNF-α alone were insufficient for promote neutrophilia in the absence of ST2 signaling. Finally, abundant production of IL-33 by cardiomyocytes was observed. These results refine our understanding of the importance of the IL-33/ST2 axis and IL-33-producing cardiomyocytes in the early acute neutrophilia induced by freshwater stingray venoms.

Revisiting venom of the sea anemone Stichodactyla haddoni: Omics techniques reveal the complete toxin arsenal of a well-studied sea anemone genus

More than a century of research on sea anemone venoms has shown that they contain a diversity of biologically active proteins and peptides. However, recent omics studies have revealed that much of the venom proteome remains unexplored. We used, for the first time, a combination of proteomic and transcriptomic techniques to obtain a holistic overview of the venom arsenal of the well-studied sea anemone Stichodactyla haddoni. A purely search-based approach to identify putative toxins in a transcriptome from tentacles regenerating after venom extraction identified 508 unique toxin-like transcripts grouped into 63 families. However, proteomic analysis of venom revealed that 52 of these toxin families are likely false positives. In contrast, the combination of transcriptomic and proteomic data enabled positive identification of 23 families of putative toxins, 12 of which have no homology known proteins or peptides. Our data highlight the importance of using proteomics of milked venom to correctly identify venom proteins/peptides, both known and novel, while minimizing false positive identifications from non-toxin homologues identified in transcriptomes of venom-producing tissues. This work lays the foundation for uncovering the role of individual toxins in sea anemone venom and how they contribute to the envenomation of prey, predators, and competitors.

BIOLOGICAL SIGNIFICANCE

Proteomic analysis of milked venom combined with analysis of a tentacle transcriptome revealed the full extent of the venom arsenal of the sea anemone Stichodactyla haddoni. This combined approach led to the discovery of 12 entirely new families of disulfide-rich peptides and proteins in a genus of anemones that have been studied for over a century.

Wasp venomic: Unravelling the toxins arsenal of Polybia paulista venom and its potential pharmaceutical applications

Polybia paulista (Hymenoptera: Vespidae) is a neotropical social wasp from southeast Brazil. As most social Hymenoptera, venom from P. paulista comprises a complex mixture of bioactive toxins ranging from low molecular weight compounds to peptides and proteins. Several efforts have been made to elucidate the molecular composition of the P. paulista venom. Data derived from proteomic, peptidomic and allergomic analyses has enhanced our understanding of the whole envenoming process caused by the insect sting. The combined use of bioinformatics, -omics- and molecular biology tools have allowed the identification, characterization, in vitro synthesis and recombinant expression of several wasp venom toxins. Some of these P. paulista - derived bioactive compounds have been evaluated for the rational design of antivenoms and the improvement of allergy specific diagnosis and immunotherapy. Molecular characterization of crude venom extract has enabled the description and isolation of novel toxins with potential biotechnological applications. Here, we review the different approaches that have been used to unravel the venom composition of P. paulista. We also describe the main groups of P. paulista - venom toxins currently identified and analyze their potential in the development of component-resolved diagnosis of allergy, and in the rational design of antivenoms and novel bioactive drugs.

Intragenome Diversity of Gene Families Encoding Toxin-like Proteins in Venomous Animals

The evolution of venoms is the story of how toxins arise and of the processes that generate and maintain their diversity. For animal venoms these processes include recruitment for expression in the venom gland, neofunctionalization, paralogous expansions, and functional divergence. The systematic study of these processes requires the reliable identification of the venom components involved in antagonistic interactions. High-throughput sequencing has the potential of uncovering the entire set of toxins in a given organism, yet the existence of non-venom toxin paralogs and the misleading effects of partial census of the molecular diversity of toxins make necessary to collect complementary evidence to distinguish true toxins from their non-venom paralogs. Here, we analyzed the whole genomes of two scorpions, one spider and one snake, aiming at the identification of the full repertoires of genes encoding toxin-like proteins. We classified the entire set of protein-coding genes into paralogous groups and monotypic genes, identified genes encoding toxin-like proteins based on known toxin families, and quantified their expression in both venom-glands and pooled tissues. Our results confirm that genes encoding toxin-like proteins are part of multigene families, and that these families arise by recruitment events from non-toxin genes followed by limited expansions of the toxin-like protein coding genes. We also show that failing to account for sequence similarity with non-toxin proteins has a considerable misleading effect that can be greatly reduced by comparative transcriptomics. Our study overall contributes to the understanding of the evolutionary dynamics of proteins involved in antagonistic interactions.