Cloning and expression of calglandulin, a new EF-hand protein from the venom glands of Bothrops insularis snake in E. coli

EP4-induced mitochondrial localization and cell migration mediated by CALML6 in human oral squamous cell carcinoma

Lymph node metastasis, primarily caused by the migration of oral squamous cell carcinoma (OSCC) cells, stands as a crucial prognostic marker. We have previously demonstrated that EP4, a subtype of the prostaglandin E2 (PGE2) receptor, orchestrates OSCC cell migration via Ca2+ signaling. The exact mechanisms by which EP4 influences cell migration through Ca2+ signaling, however, is unclear. Our study aims to clarify how EP4 controls OSCC cell migration through this pathway. We find that activating EP4 with an agonist (ONO-AE1-473) increased intracellular Ca2+ levels and the migration of human oral cancer cells (HSC-3), but not human gingival fibroblasts (HGnF). Further RNA sequencing linked EP4 to calmodulin-like protein 6 (CALML6), whose role remains undefined in OSCC. Through protein-protein interaction network analysis, a strong connection is identified between CALML6 and calcium/calmodulin-dependent protein kinase kinase 2 (CaMKK2), with EP4 activation also boosting mitochondrial function. Overexpressing EP4 in HSC-3 cells increases experimental lung metastasis in mice, whereas inhibiting CaMKK2 with STO-609 markedly lowers these metastases. This positions CaMKK2 as a potential new target for treating OSCC metastasis. Our findings highlight CALML6 as a pivotal regulator in EP4-driven mitochondrial respiration, affecting cell migration and metastasis via the CaMKK2 pathway.

Safety evaluation of Neurospora crassa mycoprotein for use as a novel meat alternative and enhancer

Cultivation of filamentous fungi to produce sustainable, nutrient rich meat replacements has recently attracted significant commercial and research interest. Here, we report evidence for the safety and nutritional value of Neurospora crassa mycoprotein, a whole mycelium food ingredient produced by fermentation and minimal downstream processing. N. crassa has a long history of human use in fermented foods and in molecular biology research. A survey of studies that used N. crassa in animal feed revealed no adverse effects to the health of the animals. Furthermore, a review of the literature found no reports of confirmed allergenicity or toxicity in humans involving N. crassa. Genomic toxigenicity analysis and in vitro testing did not identify any toxins in N. crassa mycoprotein. Two independent genomic allergenicity studies did not identify proteins that would be considered a particular risk for allergenic potential. Furthermore, nutritional analysis demonstrated that N. crassa mycoprotein is a good source of complete protein and is rich in fiber, potassium, and iron. Taken together, the presented data and the history of human use without evidence of human or animal harm indicate that foods containing N. crassa can generally be regarded as safe.

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.

An Integrated Proteomic and Transcriptomic Analysis Reveals the Venom Complexity of the Bullet Ant Paraponera clavata

A critical hurdle in ant venom proteomic investigations is the lack of databases to comprehensively and specifically identify the sequence and function of venom proteins and peptides. To resolve this, we used venom gland transcriptomics to generate a sequence database that was used to assign the tandem mass spectrometry (MS) fragmentation spectra of venom peptides and proteins to specific transcripts. This was performed alongside a shotgun liquid chromatography-mass spectrometry (LC-MS/MS) analysis of the venom to confirm that these assigned transcripts were expressed as proteins. Through the combined transcriptomic and proteomic investigation of Paraponera clavata venom, we identified four times the number of proteins previously identified using 2D-PAGE alone. In addition to this, by mining the transcriptomic data, we identified several novel peptide sequences for future pharmacological investigations, some of which conform with inhibitor cysteine knot motifs. These types of peptides have the potential to be developed into pharmaceutical or bioinsecticide peptides.

Identification of unusual peptides with new Cys frameworks in the venom of the cold-water sea anemone Cnidopus japonicus

Sea anemones (Actiniaria) are intensely popular objects of study in venomics. Order Actiniaria includes more than 1,000 species, thus presenting almost unlimited opportunities for the discovery of novel biologically active molecules. The venoms of cold-water sea anemones are studied far less than the venoms of tropical sea anemones. In this work, we analysed the molecular venom composition of the cold-water sea anemone Cnidopus japonicus. Two sets of NGS data from two species revealed molecules belonging to a variety of structural classes, including neurotoxins, toxin-like molecules, linear polypeptides (Cys-free), enzymes, and cytolytics. High-throughput proteomic analyses identified 27 compounds that were present in the venoms. Some of the toxin-like polypeptides exhibited novel Cys frameworks. To characterise their function in the venom, we heterologously expressed 3 polypeptides with unusual Cys frameworks (designated CjTL7, CjTL8, and AnmTx Cj 1c-1) in E. coli. Toxicity tests revealed that the CjTL8 polypeptide displays strong crustacean-specific toxicity, while AnmTx Cj 1c-1 is toxic to both crustaceans and insects. Thus, an improved NGS data analysis algorithm assisted in the identification of toxins with unusual Cys frameworks showing no homology according to BLAST. Our study shows the advantage of combining omics analysis with functional tests for active polypeptide discovery.

Proteomics and antivenomics of Papuan black snake (Pseudechis papuanus) venom with analysis of its toxicological profile and the preclinical efficacy of Australian antivenoms

The Papuan black snake (Pseudechis papuanus Serpentes: Elapidae) is endemic to Papua New Guinea, Indonesian Papua and Australia's Torres Strait Islands. We have investigated the biological activity and proteomic composition of its venom. The P. papuanus venom proteome is dominated by a variety (n≥18) of PLA₂s, which together account for ~90% of the venom proteins, and a set of low relative abundance proteins, including a short-neurotoxic 3FTx (3.1%), 3-4 PIII-SVMPs (2.8%), 3 cysteine-rich secretory proteins (CRISP; 2.3%) 1-3 l-amino acid oxidase (LAAO) molecules (1.6%). Probing of a P. papuanus cDNA library with specific primers resulted in the elucidation of the full-length nucleotide sequences of six new toxins, including vespryn and NGF not found in the venom proteome, and a calglandulin protein involved in toxin expression with the venom glands. Intravenous injection of P. papuanus venom in mice induced lethality, intravascular haemolysis, pulmonary congestion and oedema, and anticoagulation after intravenous injection, and these effects are mainly due to the action of PLA₂s. This study also evaluated the in vivo preclinical efficacy of Australian black snake and polyvalent Seqirus antivenoms. These antivenoms were effective in neutralising the lethal, PLA₂ and anticoagulant activities of P. papuanus venom in mice. On the other hand, all of the Seqirus antivenoms tested using an antivenomic approach exhibited strong immunorecognition of all the venom components. These preclinical results suggest that Australian Seqirus¹ antivenoms may provide paraspecific protection against P. papuanus venom in humans.

SIGNIFICANCE PARAGRAPH

The toxicological profile and proteomic composition of the venom of the Papuan black snake, Pseudechis papuanus, a large diurnal snake endemic to the southern coast of New Guinea and a handful of close offshore islands, were investigated. Intravenous injection of P. papuanus venom in mice induced intravascular hemolysis, pulmonary congestion and edema, anticoagulation, and death. These activities could be assigned to the set of PLA₂ molecules, which dominate the P. papuanus venom proteome. This study also showed that Australian Seqirus black snake or polyvalent antivenoms were effective in neutralising the lethal, PLA₂ and anticoagulant activities of the venom. These preclinical results support the continued recommendation of these Seqirus antivenoms in the clinical management of P. papuanus envenoming in Australia, Papua New Guinea or Indonesian Papua Province.

Combined snake venomics and venom gland transcriptomic analysis of Bothropoides pauloensis

Unraveling the repertoire of venom toxins of Bothropoides pauloensis was assessed by snake venomics and venom gland transcriptomic surveys. Both approaches yielded converging overall figures, pointing to metalloproteinases (~37%), PLA(2)s (26-32%), and vasoactive (bradykinin-potentiating) peptides (12-17%) as the major toxin classes. The high occurrence of SVMPs, PLA(2) molecules, vasoactive peptides, along with serine proteinases, explains the local and systemic effects observed in envenomations by B. pauloensis. Minor (<3%) C-type lectin, serine proteinase, L-amino acid oxidase, nerve growth factor, and CRISP molecules were also identified in the transcriptome and the proteome. Low abundance (0.3%) EST singletons coding for vascular endothelial growth factor (svVEGF), ohanin, hyaluronidase, and 5' nucleotidase were found only in the venom gland cDNA library. At the molecular level, the transcriptomic and proteomic datasets display low compositional concordance. In particular, although there is good agreement between transcriptome and proteome in the identity of BPPs, PLA(2) molecules and L-amino acid oxidase, both datasets strongly depart in their C-type lectin and SVMP complements. These data support the view that venom composition is influenced by transcriptional and translational mechanisms and emphasize the value of combining proteomic and transcriptomic approaches to acquire a more complete understanding of the toxinological profile and natural history of the snake venom.

A transcriptomic analysis of gene expression in the venom gland of the snake Bothrops alternatus (urutu)

Background

The genus Bothrops is widespread throughout Central and South America and is the principal cause of snakebite in these regions. Transcriptomic and proteomic studies have examined the venom composition of several species in this genus, but many others remain to be studied. In this work, we used a transcriptomic approach to examine the venom gland genes of Bothrops alternatus, a clinically important species found in southeastern and southern Brazil, Uruguay, northern Argentina and eastern Paraguay.

Results

A cDNA library of 5,350 expressed sequence tags (ESTs) was produced and assembled into 838 contigs and 4512 singletons. BLAST searches of relevant databases showed 30% hits and 70% no-hits, with toxin-related transcripts accounting for 23% and 78% of the total transcripts and hits, respectively. Gene ontology analysis identified non-toxin genes related to general metabolism, transcription and translation, processing and sorting, (polypeptide) degradation, structural functions and cell regulation. The major groups of toxin transcripts identified were metalloproteinases (81%), bradykinin-potentiating peptides/C-type natriuretic peptides (8.8%), phospholipases A₂ (5.6%), serine proteinases (1.9%) and C-type lectins (1.5%). Metalloproteinases were almost exclusively type PIII proteins, with few type PII and no type PI proteins. Phospholipases A₂ were essentially acidic; no basic PLA₂ were detected. Minor toxin transcripts were related to L-amino acid oxidase, cysteine-rich secretory proteins, dipeptidylpeptidase IV, hyaluronidase, three-finger toxins and ohanin. Two non-toxic proteins, thioredoxin and double-specificity phosphatase Dusp6, showed high sequence identity to similar proteins from other snakes. In addition to the above features, single-nucleotide polymorphisms, microsatellites, transposable elements and inverted repeats that could contribute to toxin diversity were observed.

Conclusions

Bothrops alternatus venom gland contains the major toxin classes described for other Bothrops venoms based on trancriptomic and proteomic studies. The predominance of type PIII metalloproteinases agrees with the well-known hemorrhagic activity of this venom, whereas the lower content of serine proteases and C-type lectins could contribute to less marked coagulopathy following envenoming by this species. The lack of basic PLA₂ agrees with the lower myotoxicity of this venom compared to other Bothrops species with these toxins. Together, these results contribute to our understanding of the physiopathology of envenoming by this species.

Venomics of the spider Ornithoctonus huwena based on transcriptomic versus proteomic analysis

The spider Ornithoctonus huwena is a venomous spider found in southern China. Its venom is a complex mixture of numerous biologically active components. In this study, 41 novel unique transcripts encoding cellular proteins or other possible venom components were generated from the previously constructed cDNA library. These proteins were also annotated by KOG (eukaryotic orthologous group) and GO (gene ontology) terms. A novel cellular transcript contig encoding an EF-hand protein (named HWEFHP1) was found, which might be involved in the secretion of toxins in the venom glands. In order to have an overview of the molecular diversity of the O. huwena venom, the datasets of all the transcripts, peptides and proteins known so far were analyzed. A comparison of the data obtained through a proteomic versus a transcriptomic approach, revealed that only 15 putative cystine knot toxins (CKTs) were identified by both approaches, 29 transcripts coding for CKTs were found in the transcriptome but not as translated peptides in the venom proteome. However, no cellular protein with identical molecular weight was identified by both approaches. Our data may contribute to a deeper understanding of the biology and ecology of O. huwena and the relationship between structure and function of individual toxins.

Transcriptome analysis of the venom glands of the Chinese wolf spider Lycosa singoriensis

The wolf spider Lycosa singoriensis is a hunting spider with a widespread distribution in northwest China. The venom gland of spiders, which is a very specialized secretory tissue, can secrete abundant and complex toxin components. To extensively examine the transcripts expressed in the venom glands of L. singoriensis, we generated 833 expressed sequence tags (ESTs) from a directional cDNA library. Toxin-like sequences account for 69.1% of these ESTs, 17.3% are similar to cellular transcripts and 13.6% have no significant similarity to any known sequences. Here, we identified 223 novel toxin-like sequences, which can be classified into six different superfamilies; that means a novel potential source of ligands for varied ion channels was discovered. With the aid of Gene Ontology terms and homology to eukaryotic orthologous groups, the annotation of cellular transcripts revealed some cellular processes important for the toxin secretion of venom glands including protein synthesis, protein folding, tuned post-translational processing and trafficking, etc.

Recombinant expression and affinity purification of snake venom gland parvalbumin in Escherichia coli

Parvalbumins (PV) are small, acidic, water soluble and calcium-binding proteins generally present in muscular and nervous tissues. In the present study, we identified and characterized a cDNA clone encoding PV, named AplPV, from a snake (Agkistrodon piscivorus leucostoma) venom gland cDNA library. AplPV belongs to EF-hand proteins with six alpha-helices constituting three EF-hand domains. The deduced amino acid sequence of AplPV is 91% and 68% identical to the previously characterized PVs of Boa constrictor and Cyprinus carpio, respectively. The full-length cDNA was subcloned into the expression vector pGEX and transformed into Escherichia coli (E.coli) to produce recombinant protein. The bacterially expressed GST-AplPV fusion protein was highly expressed, and effectively purified by Glutathione-Sepharose affinity chromatography. A high concentration of thrombin protease specifically cleaved and removed the GST tag from fusion protein, and further purified by Benzamidine column for removal of thrombin protease. As a result, the 12 kDa AplPV recombinant protein alone was purified. To investigate the tissue-specific biological occurrence of AplPV, a polyclonal antibody (anti-AplPV-antibody) was raised against GST-AplPV fusion protein in rabbit. Western blot analysis revealed that immunoreactive bands were exhibited in both recombinant protein and samples of venom glands, but not in any crude venom. This specific occurrence indicates a specialized function of AplPV in snake venom glands.

Alpha1-adrenoceptors trigger the snake venom production cycle in secretory cells by activating phosphatidylinositol 4,5-bisphosphate hydrolysis and ERK signaling pathway

Loss of venom from the venom gland after biting or manual extraction leads to morphological changes in venom secreting cells and the start of a cycle of production of new venom. We have previously shown that stimulation of both alpha- and beta-adrenoceptors in the secretory cells of the venom gland is essential for the onset of the venom production cycle in Bothrops jararaca. We investigated the signaling pathway by which the alpha-adrenoceptor initiates the venom production cycle. Our results show that the alpha(1)-adrenoceptor subtype is present in venom gland of the snake. In quiescent cells, stimulation of alpha(1)-adrenoceptor with phenylephrine increased the total inositol phosphate concentration, and this effect was blocked by the phospholipase C inhibitor U73122. Phenylephrine mobilized Ca(2+) from thapsigargin-sensitive stores and increased protein kinase C activity. In addition, alpha(1)-adrenoceptor stimulation increased the activity of ERK 1/2, partially via protein kinase C. Using RT-PCR approach we obtained a partial sequence of a snake alpha(1)-adrenoceptor (260 bp) with higher identity with alpha(1D) and alpha(1B)-adrenoceptors from different species. These results suggest that alpha(1)-adrenoceptors in the venom secreting cells are probably coupled to a G(q) protein and trigger the venom production cycle by activating the phosphatidylinositol 4,5-bisphosphate and ERK signaling pathway.

Expression pattern of three-finger toxin and phospholipase A2 genes in the venom glands of two sea snakes, Lapemis curtus and Acalyptophis peronii: comparison of evolution of these toxins in land snakes, sea kraits and sea snakes

Background

Snake venom composition varies widely both among closely related species and within the same species, based on ecological variables. In terrestrial snakes, such variation has been proposed to be due to snakes' diet. Land snakes target various prey species including insects (arthropods), lizards (reptiles), frogs and toads (amphibians), birds (aves), and rodents (mammals), whereas sea snakes target a single vertebrate class (fishes) and often specialize on specific types of fish. It is therefore interesting to examine the evolution of toxins in sea snake venoms compared to that of land snakes.

Results

Here we describe the expression of toxin genes in the venom glands of two sea snakes, Lapemis curtus (Spine-bellied Sea Snake) and Acalyptophis peronii (Horned Sea Snake), two members of a large adaptive radiation which occupy very different ecological niches. We constructed cDNA libraries from their venom glands and sequenced 214 and 192 clones, respectively. Our data show that despite their explosive evolutionary radiation, there is very little variability in the three-finger toxin (3FTx) as well as the phospholipase A₂ (PLA₂) enzymes, the two main constituents of Lapemis curtus and Acalyptophis peronii venom. To understand the evolutionary trends among land snakes, sea snakes and sea kraits, pairwise genetic distances (intraspecific and interspecific) of 3FTx and PLA₂ sequences were calculated. Results show that these proteins appear to be highly conserved in sea snakes in contrast to land snakes or sea kraits, despite their extremely divergent and adaptive ecological radiation.

Conclusion

Based on these results, we suggest that streamlining in habitat and diet in sea snakes has possibly kept their toxin genes conserved, suggesting the idea that prey composition and diet breadth may contribute to the diversity and evolution of venom components.

Cellular autoreactivity against heat shock protein 60 in renal transplant patients: peripheral and graft-infiltrating responses

Autoreactivity to heat shock protein 60 (Hsp60) has been implicated in the pathogenesis and regulation of chronic inflammation, especially in autoimmune diseases. In transplantation, there is a lack of information regarding the cytokine profile and specificity of cells that recognize self-Hsp60 as well as the kinetics of autoreactivity following transplantation. We studied the cellular reactivity of peripheral and graft-infiltrating lymphocytes against Hsp60 in renal transplant patients. Cytokine production induced by this protein in peripheral blood mononuclear cells indicated a predominance of interleukin (IL)-10 during the late post-transplantation period, mainly in response to intermediate and C-terminal peptides. Patients with chronic rejection presented reactivity to Hsp60 with a higher IL-10/interferon (IFN)-gamma ratio compared to long-term clinically stable patients. Graft-infiltrating T cell lines, cocultured with antigen-presenting cells, preferentially produced IL-10 after Hsp60 stimulation. These results suggest that, besides its proinflammatory activity, autoreactivity to Hsp60 in transplantation may also have a regulatory role.

Transcriptome analysis of Deinagkistrodon acutus venomous gland focusing on cellular structure and functional aspects using expressed sequence tags

Background

The snake venom gland is a specialized organ, which synthesizes and secretes the complex and abundant toxin proteins. Though gene expression in the snake venom gland has been extensively studied, the focus has been on the components of the venom. As far as the molecular mechanism of toxin secretion and metabolism is concerned, we still knew a little. Therefore, a fundamental question being arisen is what genes are expressed in the snake venom glands besides many toxin components?

Results

To examine extensively the transcripts expressed in the venom gland of Deinagkistrodon acutus and unveil the potential of its products on cellular structure and functional aspects, we generated 8696 expressed sequence tags (ESTs) from a non-normalized cDNA library. All ESTs were clustered into 3416 clusters, of which 40.16% of total ESTs belong to recognized toxin-coding sequences; 39.85% are similar to cellular transcripts; and 20.00% have no significant similarity to any known sequences. By analyzing cellular functional transcripts, we found high expression of some venom related genes and gland-specific genes, such as calglandulin EF-hand protein gene and protein disulfide isomerase gene. The transcripts of creatine kinase and NADH dehydrogenase were also identified at high level. Moreover, abundant cellular structural proteins similar to mammalian muscle tissues were also identified. The phylogenetic analysis of two snake venom toxin families of group III metalloproteinase and serine protease in suborder Colubroidea showed an early single recruitment event in the viperids evolutionary process.

Conclusion

Gene cataloguing and profiling of the venom gland of Deinagkistrodon acutus is an essential requisite to provide molecular reagents for functional genomic studies needed for elucidating mechanisms of action of toxins and surveying physiological events taking place in the very specialized secretory tissue. So this study provides a first global view of the genetic programs for the venom gland of Deinagkistrodon acutus described so far and an insight into molecular mechanism of toxin secreting.

All sequences data reported in this paper have been submitted into the public database [GenBank: DV556511-DV565206].

Cloning, characterization, and structural analysis of a C-type lectin from Bothrops insularis (BiL) venom

Lectins are carbohydrate-binding molecules that mediate a variety of biological processes. In this work, we identify and characterize a lectin from Bothrops insularis venom, with respect to its biochemical properties and theoretical structure. Initially, from a venom gland cDNA library, we cloned and sequenced a cDNA encoding a protein with high identity to snake venom lectins. A lectin molecule was purified to homogeneity from the venom by affinity column and gel filtration. This protein named BiL displayed hemagglutinating activity that was inhibited by galactose, lactose, and EDTA. Mass spectrometry analysis and sodium dodecyl sulfate-polyacrylamide gel electrophoresis showed that BiL is a disulfide-linked dimeric protein consisting of monomers with 16,206 m/z. The amino acid sequence, deduced from its cDNA sequence, was confirmed by Edman sequencing and by peptide mass fingerprint analysis. BiL shows similarity to other C-type lectin family members. Modeling studies provide insights into BiL dimeric structure and its structural determinants for carbohydrate and calcium binding.