Ohanin, a novel protein from king cobra venom, induces hypolocomotion and hyperalgesia in mice

A retrospective analysis of epidemiology, clinical features of envenomation, and in-patient management of snakebites in a model secondary hospital of Assam, North-east India

Assam, a Northeastern State of India, is inhabited by several venomous snake species causing substantial morbidity and mortality. The data on the epidemiology of snakebites and their management is underreported in this region. Hence, a secondary health-based retrospective study was carried out at Demow Model Hospital, Sivasagar, Assam, to evaluate the clinical and epidemiological profile of snakebite cases reported in this rural hospital and their management. Snakebites occurring between April 2018 to August 2022 were reviewed based on socio-demographic details of the patient, clinical symptoms, and treatment using a standard questionnaire. Out of the 1011 registered snakebite cases, 139 patients (13.7%) counted for venomous bites, among which 92 patients (66.19%) accounted for viper bites (green pit viper and Salazar's pit viper), and 30 patients (21.5%) were bitten by elapid snakes (Indian monocled Cobra, banded krait, and greater/lesser black krait). A maximum number of snakebite cases (80.5%) were reported from the interior rural villages and documented from July to September (51.3%). Elapid snake envenomed patients, except one, were successfully treated with commercial antivenom, neostigmine, and glycopyrrolate. Because commercial polyvalent antivenom against "Big Four" venomous snakes of India showed poor neutralization of pit-vipers envenomation; therefore, pit-viper bite patients were treated with repurposed drugs magnesium sulfate and glycerin compression dressing. Adverse serum reactions were reported only in 3 (11.1%) cases. The preventive measures and facilities adopted at the Demow Model Hospital significantly reduce snakebite death and morbidity; therefore, they can be s practised across various states in India as a prototype.

Venom comparisons of endemic and micro-endemic speckled rattlesnakes Crotalus mitchellii, C. polisi and C. thalassoporus from Baja California Peninsula

A high diversity of rattlesnake species can be found in the Baja California peninsula and the island of the Gulf of California, nevertheless, their venom has been poorly evaluated. The aim of this work was to present the first characterization of endemic Crotalus mitchellii, micro endemic C. polisi and C. thalassoporus venoms. All samples provoke human plasma coagulation showing doses in the rank of 2.3-41.0 μg and also produce rapid hydrolysis of the alpha chain of bovine fibrinogen while the beta chain is attacked at larger incubation periods by C. polisi and especially by C. thalassoporus. Phospholipase activity ranging from 23.2 to 173.8 U/mg. The venoms of C. thalassoporus and C. polisi show very high hemorrhagic activity (from 0.03 to 0.31 μg). A total of 130 toxin-related proteins were identified and classified into ten families. Crotalus mitchellii venom was characterized by high abundance of crotoxin-like and other phospholipase proteins (34.5%) and serine proteinases (29.8%). Crotalus polisi showed a similar proportion of metalloproteinases (34%) and serine proteinases (22.8%) components with important contribution of C-type lectins (14.3%) and CRiSP (14.0%) proteins. Venom of C. thalassoporus is dominated by metalloproteases that amount to more than 66% of total toxin proteins. These results provide a foundation for comprehending the biological, ecological and evolutionary significance of venom composition of speckled rattlesnake from the Baja California peninsula.

Venomics and antivenomics of Indian spectacled cobra (Naja naja) from the Western Ghats

Venom proteome profiling of Naja naja from the Western Ghats region in Kerala was achieved through SDS-PAGE and RP-HPLC followed by Q-TOF LC-MS/MS analysis, incorporating PEAKS and Novor assisted de novo sequencing methodologies. A total of 115 proteins distributed across 17 different enzymatic and non-enzymatic venom protein families were identified through conventional and 39 peptides through homology-driven proteomics approaches. Fourteen peptides derived through de novo complements the Mascot data indicating the importance of homology-driven approaches in improving protein sequence information. Among the protein families identified, glutathione peroxidase and endonuclease were reported for the first time in the Indian cobra venom. Immunological cross-reactivity assessed using Indian polyvalent antivenoms suggested that VINS showed better EC₅₀ (2.48 µg/mL) value than that of PSAV (6.04 µg/mL) and Virchow (6.03 µg/mL) antivenoms. Western blotting experiments indicated that all the antivenoms elicited poor binding specificities, especially towards low molecular mass proteins. Second-generation antivenomics studies revealed that VINS antivenom was less efficient to detect many low molecular mass proteins such as three-finger toxins and Kunitz-type serine protease Inhibitors. Taken together, the present study enabled a large-scale characterization of the venom proteome of Naja naja from the Western Ghats and emphasized the need for developing more efficient antivenoms.

Snake Venom Proteomics of Samar Cobra (Naja samarensis) from the Southern Philippines: Short Alpha-Neurotoxins as the Dominant Lethal Component Weakly Cross-Neutralized by the Philippine Cobra Antivenom

The Samar Cobra, Naja samarensis, is endemic to the southern Philippines and is a WHO-listed Category 1 venomous snake species of medical importance. Envenomation caused by N. samarensis results in neurotoxicity, while there is no species-specific antivenom available for its treatment. The composition and neutralization of N. samarensis venom remain largely unknown to date. This study thus aimed to investigate the venom proteome of N. samarensis for a comprehensive profiling of the venom composition, and to examine the immunorecognition as well as neutralization of its toxins by a hetero-specific antivenom. Applying C₁₈ reverse-phase high-performance liquid chromatography (RP-HPLC) and tandem mass spectrometry (LC-MS/MS), three-finger toxins (3FTx) were shown to dominate the venom proteome by 90.48% of total venom proteins. Other proteins in the venom comprised snake venom metalloproteinases, phospholipases A_2, cysteine-rich secretory proteins, venom nerve growth factors, L-amino acid oxidases and vespryn, which were present at much lower abundances. Among all, short-chain alpha-neurotoxins (SαNTX) were the most highly expressed toxin within 3FTx family, constituting 65.87% of the total venom proteins. The SαNTX is the sole neurotoxic component of the venom and has an intravenous median lethal dose (LD₅₀) of 0.18 μg/g in mice. The high abundance and low LD₅₀ support the potent lethal activity of N. samarensis venom. The hetero-specific antivenom, Philippine Cobra Antivenom (PCAV, raised against Naja philippinensis) were immunoreactive toward the venom and its protein fractions, including the principal SαNTX. In efficacy study, PCAV was able to cross-neutralize the lethality of SαNTX albeit the effect was weak with a low potency of 0.20 mg/ml (defined as the amount of toxin completely neutralized per milliliter of the antivenom). With a volume of 5 ml, each vial of PCAV may cross-neutralize approximately 1 mg of the toxin in vivo. The findings support the potential para-specific use of PCAV in treating envenomation caused by N. samarensis while underscoring the need to improve the potency of its neutralization activity, especially against the highly lethal alpha-neurotoxins.

Major histocompatibility complex genes and locus organization in the Komodo dragon (Varanus komodoensis)

We performed a meta-analysis of the newly assembled Komodo dragon (Varanus komodoensis) genome to characterize the major histocompatibility complex (MHC) of the species. The MHC gene clusters of the Komodo dragon are gene dense, complex, and contain counterparts of many genes of the human MHC. Our analysis identified 20 contigs encompassing ~ 6.9 Mbp of sequence with 223 annotated genes of which many are predicted orthologs to the genes of the human MHC. These MHC contigs range in size from 13.2 kb to 21.5 Mbp, contain an average of one gene per 30 kb, and are thought to occur on at least two chromosomes. Eight contigs, each > 100 kb, could be aligned to the human MHC based on gene content, and these represent gene clusters found in each of the recognized mammalian MHC subregions. The MHC of the Komodo dragon shares organizational features of other non-mammalian taxa. Multiple class Iα and class IIβ genes are indicated, with linkage between classical class I and immunoproteasome genes and between framework class I genes and genes associated with the mammalian class III subregion. These findings are supported in both Komodo genome assemblies and provide new insight into the MHC organization of these unique squamate reptiles.

King Cobra and snakebite envenomation: on the natural history, human-snake relationship and medical importance of Ophiophagus hannah

King Cobra (Ophiophagus hannah) has a significant place in many cultures, and is a medically important venomous snake in the world. Envenomation by this snake is highly lethal, manifested mainly by neurotoxicity and local tissue damage. King Cobra may be part of a larger species complex, and is widely distributed across Southeast Asia, southern China, northern and eastern regions as well as the Western Ghats of India, indicating potential geographical variation in venom composition. There is, however, only one species-specific King Cobra antivenom available worldwide that is produced in Thailand, using venom from the snake of Thai origin. Issues relating to the management of King Cobra envenomation (e.g., variation in the composition and toxicity of the venom, limited availability and efficacy of antivenom), and challenges faced in the research of venom (in particular proteomics), are rarely addressed. This article reviews the natural history and sociocultural importance of King Cobra, cases of snakebite envenomation caused by this species, current practice of management (preclinical and clinical), and major toxinological studies of the venom with a focus on venom proteomics, toxicity and neutralization. Unfortunately, epidemiological data of King Cobra bite is scarce, and venom proteomes reported in various studies revealed marked discrepancies in details. Challenges, such as inconsistency in snake venom sampling, varying methodology of proteomic analysis, lack of mechanistic and antivenomic studies, and controversy surrounding antivenom use in treating King Cobra envenomation are herein discussed. Future directions are proposed, including the effort to establish a standard, comprehensive Pan-Asian proteomic database of King Cobra venom, from which the venom variation can be determined. Research should be undertaken to characterize the toxin antigenicity, and to develop an antivenom with improved efficacy and wider geographical utility. The endeavors are aligned with the WHO´s roadmap that aims to reduce the disease burden of snakebite by 50% before 2030.

Proteomic Investigations of Two Pakistani Naja Snake Venoms Species Unravel the Venom Complexity, Posttranslational Modifications, and Presence of Extracellular Vesicles

Latest advancement of omics technologies allows in-depth characterization of venom compositions. In the present work we present a proteomic study of two snake venoms of the genus Naja i.e., Naja naja (black cobra) and Naja oxiana (brown cobra) of Pakistani origin. The present study has shown that these snake venoms consist of a highly diversified proteome. Furthermore, the data also revealed variation among closely related species. High throughput mass spectrometric analysis of the venom proteome allowed to identify for the N. naja venom 34 protein families and for the N. oxiana 24 protein families. The comparative evaluation of the two venoms showed that N. naja consists of a more complex venom proteome than N. oxiana venom. Analysis also showed N-terminal acetylation (N-ace) of a few proteins in both venoms. To the best of our knowledge, this is the first study revealing this posttranslational modification in snake venom. N-ace can shed light on the mechanism of regulation of venom proteins inside the venom gland. Furthermore, our data showed the presence of other body proteins, e.g., ankyrin repeats, leucine repeats, zinc finger, cobra serum albumin, transferrin, insulin, deoxyribonuclease-2-alpha, and other regulatory proteins in these venoms. Interestingly, our data identified Ras-GTpase type of proteins, which indicate the presence of extracellular vesicles in the venom. The data can support the production of distinct and specific anti-venoms and also allow a better understanding of the envenomation and mechanism of distribution of toxins. Data are available via ProteomeXchange with identifier PXD018726.

Quantitative proteomics of Naja annulifera (sub-Saharan snouted cobra) venom and neutralization activities of two antivenoms in Africa

Envenomation by Naja annulifera (snouted cobra), a non-spitting African cobra, can result in local tissue damage and fatal paralysis but a species-specific antivenom treatment is currently lacking. In this study, we investigated the quantitative proteome of N. annulifera venom, incorporating HPLC and LC-MS/MS to elucidate the venom toxicity. The immunoreactivities and in vivo neutralization activities of two hetero-specific antivenom products (Premium Serums Pan Africa polyvalent antivenom, PANAF and VINS African polyvalent antivenom, VAPAV) against the venom were subsequently examined. N. annulifera venom comprises 10 toxin families, with three-finger toxin (3FTx) being the most abundantly expressed (~78%). Within 3FTx, cytotoxin is the most dominant form and made up three-quarter of the venom bulk (~74%), whereas alpha-neurotoxins constitute <4% of the total venom proteins. Phospholipase A2 was undetected in the venom proteome, consistent with the unusual absence of PLA2 from the venoms of cobras in the Uraeus subgenus. In ELISA, PANAF and VAPAV showed comparable immunoreactivity toward the protein antigens of N. annulifera venom. These antivenoms, despite being raised against hetero-specific venoms, were capable of cross-neutralizing the lethal effect of N. annulifera venom in mice, with PANAF being marginally more potent.

The Indian cobra reference genome and transcriptome enables comprehensive identification of venom toxins

Snakebite envenoming is a serious and neglected tropical disease that kills ~100,000 people annually. High-quality, genome-enabled comprehensive characterization of toxin genes will facilitate development of effective humanized recombinant antivenom. We report a de novo near-chromosomal genome assembly of Naja naja, the Indian cobra, a highly venomous, medically important snake. Our assembly has a scaffold N50 of 223.35 Mb, with 19 scaffolds containing 95% of the genome. Of the 23,248 predicted protein-coding genes, 12,346 venom-gland-expressed genes constitute the ‘venom-ome’ and this included 139 genes from 33 toxin families. Among the 139 toxin genes were 19 ‘venom-ome-specific toxins’ (VSTs) that showed venom-gland-specific expression, and these probably encode the minimal core venom effector proteins. Synthetic venom reconstituted through recombinant VST expression will aid in the rapid development of safe and effective synthetic antivenom. Additionally, our genome could serve as a reference for snake genomes, support evolutionary studies and enable venom-driven drug discovery.

Analysis of a near-chromosomal genome assembly and transcriptome profiling of the Indian cobra identifies genes expressed in the venom glands. These data should help develop a new antivenom.

Proteomic insights into short neurotoxin-driven, highly neurotoxic venom of Philippine cobra (Naja philippinensis) and toxicity correlation of cobra envenomation in Asia

The Philippine cobra, Naja philippinensis, is a WHO Category 1 venomous snake of medical importance responsible for fatal envenomation in the northern Philippines. To elucidate the venom proteome and pathophysiology of envenomation, N. philippinensis venom proteins were decomplexed with reverse-phase high-performance liquid chromatography, and protein fractions were subsequently digested with trypsin, followed by nano-liquid chromatography-tandem mass spectrometry analysis and data mining. Three-finger toxins (3FTX, 66.64% of total venom proteins) and phospholipases A₂ (PLA₂, 22.88%) constitute the main bulk of venom proteome. Other proteins are present at low abundances (<4% each); these include metalloproteinase, serine protease, cobra venom factor, cysteine-rich secretory protein, vespryn, phosphodiesterase, 5' nucleotidase and nerve growth factor. In the three-finger toxin family, the alpha-neurotoxins comprise solely short neurotoxins (SNTX, 44.55%), supporting that SNTX is the principal toxin responsible for neuromuscular paralysis and lethality reported in clinical envenomation. Cytotoxins (CTX) are the second most abundant 3FTX proteins in the venom (21.31%). The presence of CTX correlates with the venom cytotoxic effect, which is more prominent in murine cells than in human cells. From the practical standpoint, SNTX-driven neuromuscular paralysis is significant in N. philippinensis envenomation. Antivenom production and treatment should be tailored accordingly to ensure effective neutralization of SNTX. BIOLOGICAL SIGNIFICANCE: The venom proteome of Naja philippinensis, the Philippine cobra, is unravelled for the first time. Approximately half the protein bulk of the venom is made up of short neurotoxins (44.55% of the total venom proteins). As the only alpha-neurotoxins present in the venom, short neurotoxins are the causative toxins of the post-synaptic blockade and fast-onset neuromuscular paralysis in N. philippinensis envenomation. A substantial amount of cytotoxins (21.31%) was also detected in N. philippinensis venom, supporting that the venom can be cytotoxic although the effect is much weaker in human cells compared to murine cells. The finding is consistent with the low incidence of local tissue necrosis in N. philippinensis envenomation, although this does not negate the need for monitoring and care of bite wound in the patients.

Exploring the Diversity and Novelty of Toxin Genes in Naja sumatrana, the Equatorial Spitting Cobra from Malaysia through De Novo Venom-Gland Transcriptomics

The equatorial spitting cobra, Naja sumatrana, is a distinct species of medically important venomous snakes, listed as WHO Category 1 in Southeast Asia. The diversity of its venom genes has not been comprehensively examined, although a few toxin sequences annotated to Naja sputatrix were reported previously through cloning studies. To investigate this species venom genes’ diversity, de novo venom-gland transcriptomics of N. sumatrana from West Malaysia was conducted using next-generation sequencing technology. Genes encoding toxins represented only 60 of the 55,396 transcripts, but were highly expressed, contributing to 79.22% of total gene expression (by total FPKM) in the venom-glands. The toxin transcripts belong to 21 families, and 29 transcripts were further identified as full-length. Three-finger toxins (3FTx) composed of long, short, and non-conventional groups, constituted the majority of toxin transcripts (91.11% of total toxin FPKM), followed by phospholipase A₂ (PLA₂, 7.42%)—which are putatively pro-inflammatory and cytotoxic. The remaining transcripts in the 19 families were expressed at extremely low levels. Presumably, these toxins were associated with ancillary functions. Our findings unveil the diverse toxin genes unique to N. sumatrana, and provide insights into the pathophysiology of N. sumatrana envenoming.

Proteomic Characterization of Two Medically Important Malaysian Snake Venoms, Calloselasma rhodostoma (Malayan Pit Viper) and Ophiophagus hannah (King Cobra)

Calloselasma rhodostoma (CR) and Ophiophagus hannah (OH) are two medically important snakes found in Malaysia. While some studies have described the biological properties of these venoms, feeding and environmental conditions also influence the concentration and distribution of snake venom toxins, resulting in variations in venom composition. Therefore, a combined proteomic approach using shotgun and gel filtration chromatography, analyzed by tandem mass spectrometry, was used to examine the composition of venoms from these Malaysian snakes. The analysis revealed 114 proteins (15 toxin families) and 176 proteins (20 toxin families) in Malaysian Calloselasma rhodostoma and Ophiophagus hannah species, respectively. Flavin monoamine oxidase, phospholipase A₂, phosphodiesterase, snake venom metalloproteinase, and serine protease toxin families were identified in both venoms. Aminopeptidase, glutaminyl-peptide cyclotransferase along with ankyrin repeats were identified for the first time in CR venom, and insulin, c-type lectins/snaclecs, hepatocyte growth factor, and macrophage colony-stimulating factor together with tumor necrosis factor were identified in OH venom for the first time. Our combined proteomic approach has identified a comprehensive arsenal of toxins in CR and OH venoms. These data may be utilized for improved antivenom production, understanding pathological effects of envenoming, and the discovery of biologically active peptides with medical and/or biotechnological value.

In-Depth Venome of the Brazilian Rattlesnake Crotalus durissus terrificus: An Integrative Approach Combining Its Venom Gland Transcriptome and Venom Proteome

Snake venoms are complex mixtures mainly composed of proteins and small peptides. Crotoxin is one of the most studied components from Crotalus venoms, but many other components are less known due to their low abundance. The venome of Crotalus durissus terrificus, the most lethal Brazilian snake, was investigated by combining its venom gland transcriptome and proteome to create a holistic database of venom compounds unraveling novel toxins. We constructed a cDNA library from C. d. terrificus venom gland using the Illumina platform and investigated its venom proteome through high resolution liquid chromotography-tandem mass spectrometry. After integrating data from both data sets, more than 30 venom components classes were identified by the transcriptomic analysis and 15 of them were detected in the venom proteome. However, few of them (PLA₂, SVMP, SVSP, and VEGF) were relatively abundant. Furthermore, only seven expressed transcripts contributed to ∼82% and ∼73% of the abundance in the transcriptome and proteome, respectively. Additionally, novel venom proteins are reported, and we highlight the importance of using different databases to perform the data integration and discuss the structure of the venom components-related transcripts identified. Concluding, this research paves the way for novel investigations and discovery of future pharmacological agents or targets in the antivenom therapy.

Proteomic Deep Mining the Venom of the Red-Headed Krait, Bungarus flaviceps

The use of -omics technologies allows for the characterization of snake venom composition at a fast rate and at high levels of detail. In the present study, we investigated the protein content of Red-headed Krait (Bungarus flaviceps) venom. This analysis revealed a high diversity of snake venom protein families, as evidenced by high-throughput mass spectrometric analysis. We found all six venom protein families previously reported in a transcriptome study of the venom gland of B. flaviceps, including phospholipases A₂ (PLA₂s), Kunitz-type serine proteinase inhibitors (KSPIs), three-finger toxins (3FTxs), cysteine-rich secretory proteins (CRISPs), snaclecs, and natriuretic peptides. A combined approach of automated database searches and de novo sequencing of tandem mass spectra, followed by sequence similarity searches, revealed the presence of 12 additional toxin families. De novo sequencing alone was able to identify 58 additional peptides, and this approach contributed significantly to the comprehensive description of the venom. Abundant protein families comprise 3FTxs (22.3%), KSPIs (19%), acetylcholinesterases (12.6%), PLA₂s (11.9%), venom endothelial growth factors (VEGFs, 8.4%), nucleotidases (4.3%), and C-type lectin-like proteins (snaclecs, 3.3%); an additional 11 toxin families are present at significantly lower concentrations, including complement depleting factors, a family not previously detected in Bungarus venoms. The utility of a multifaceted approach toward unraveling the proteome of snake venoms, employed here, allowed detection of even minor venom components. This more in-depth knowledge of the composition of B. flaviceps venom facilitates a better understanding of snake venom molecular evolution, in turn contributing to more effective treatment of krait bites.

Marine Leech Anticoagulant Diversity and Evolution

Leeches (Annelida: Hirudinea) possess powerful salivary anticoagulants and, accordingly, are frequently employed in modern, authoritative medicine. Members of the almost exclusively marine family Piscicolidae account for 20% of leech species diversity, and they feed on host groups (e.g., sharks) not encountered by their freshwater and terrestrial counterparts. Moreover, some species of Ozobranchidae feed on endangered marine turtles and have been implicated as potential vectors for the tumor-associated turtle herpesvirus. In spite of their ecological importance and unique host associations, there is a distinct paucity of data regarding the salivary transcriptomes of either of these families. Using next-generation sequencing, we profiled transcribed, putative anticoagulants and other salivary bioactive compounds that have previously been linked to blood feeding from 7 piscicolid species (3 elasmobranch feeders; 4 non-cartilaginous fish feeders) and 1 ozobranchid species (2 samples). In total, 149 putative anticoagulants and bioactive loci were discovered in varying constellations throughout the different samples. The putative anticoagulants showed a broad spectrum of described antagonistic pathways, such as inhibition of factor Xa and platelet aggregation, which likely have similar bioactive roles in marine fish and turtles. A transcript with homology to ohanin, originally isolated from king cobras, was found in Cystobranchus vividus but is otherwise unknown from leeches. Estimation of selection pressures for the putative anticoagulants recovered evidence for both positive and purifying selection along several isolated branches in the gene trees, and positive selection was also estimated for a few select codons in a variety of marine species. Similarly, phylogenetic analyses of the amino acid sequences for several anticoagulants indicated divergent evolution.

Coralsnake Venomics: Analyses of Venom Gland Transcriptomes and Proteomes of Six Brazilian Taxa

Venom gland transcriptomes and proteomes of six Micrurus taxa (M. corallinus, M. lemniscatus carvalhoi, M. lemniscatus lemniscatus, M. paraensis, M. spixii spixii, and M. surinamensis) were investigated, providing the most comprehensive, quantitative data on Micrurus venom composition to date, and more than tripling the number of Micrurus venom protein sequences previously available. The six venomes differ dramatically. All are dominated by 2-6 toxin classes that account for 91-99% of the toxin transcripts. The M. s. spixii venome is compositionally the simplest. In it, three-finger toxins (3FTxs) and phospholipases A₂ (PLA₂s) comprise >99% of the toxin transcripts, which include only four additional toxin families at levels ≥0.1%. Micrurus l. lemniscatus venom is the most complex, with at least 17 toxin families. However, in each venome, multiple structural subclasses of 3FTXs and PLA₂s are present. These almost certainly differ in pharmacology as well. All venoms also contain phospholipase B and vascular endothelial growth factors. Minor components (0.1-2.0%) are found in all venoms except that of M. s. spixii. Other toxin families are present in all six venoms at trace levels (<0.005%). Minor and trace venom components differ in each venom. Numerous novel toxin chemistries include 3FTxs with previously unknown 8- and 10-cysteine arrangements, resulting in new 3D structures and target specificities. 9-cysteine toxins raise the possibility of covalent, homodimeric 3FTxs or heterodimeric toxins with unknown pharmacologies. Probable muscarinic sequences may be reptile-specific homologs that promote hypotension via vascular mAChRs. The first complete sequences are presented for 3FTxs putatively responsible for liberating glutamate from rat brain synaptosomes. Micrurus C-type lectin-like proteins may have 6-9 cysteine residues and may be monomers, or homo- or heterodimers of unknown pharmacology. Novel KSPIs, 3× longer than any seen previously, appear to have arisen in three species by gene duplication and fusion. Four species have transcripts homologous to the nociceptive toxin, (MitTx) α-subunit, but all six species had homologs to the β-subunit. The first non-neurotoxic, non-catalytic elapid phospholipase A₂s are reported. All are probably myonecrotic. Phylogenetic analysis indicates that the six taxa diverged 15-35 million years ago and that they split from their last common ancestor with Old World elapines nearly 55 million years ago. Given their early diversification, many cryptic micrurine taxa are anticipated.

Comparison of proteomic profiles of the venoms of two of the 'Big Four' snakes of India, the Indian cobra (Naja naja) and the common krait (Bungarus caeruleus), and analyses of their toxins

Snake venoms are mixtures of biologically-active proteins and peptides, and several studies have described the characteristics of some of these toxins. However, complete proteomic profiling of the venoms of many snake species has not yet been done. The Indian cobra (Naja naja) and common krait (Bungarus caeruleus) are elapid snake species that are among the 'Big Four' responsible for the majority of human snake envenomation cases in India. As understanding the composition and complexity of venoms is necessary for successful treatment of envenomation in humans, we utilized three different proteomic profiling approaches to characterize these venoms: i) one-dimensional SDS-PAGE coupled with in-gel tryptic digestion and electrospray tandem mass spectrometry (ESI-LC-MS/MS) of individual protein bands; ii) in-solution tryptic digestion of crude venoms coupled with ESI-LC-MS/MS; and iii) separation by gel-filtration chromatography coupled with tryptic digestion and ESI-LC-MS/MS of separated fractions. From the generated data, 81 and 46 different proteins were identified from N. naja and B. caeruleus venoms, respectively, belonging to fifteen different protein families. Venoms from both species were found to contain a variety of phospholipases A₂ and three-finger toxins, whereas relatively higher numbers of snake venom metalloproteinases were found in N. naja compared to B. caeruleus venom. The analyses also identified less represented venom proteins including L-amino acid oxidases, cysteine-rich secretory proteins, 5'-nucleotidases and venom nerve growth factors. Further, Kunitz-type serine protease inhibitors, cobra venom factors, phosphodiesterases, vespryns and aminopeptidases were identified in the N. naja venom, while acetylcholinesterases and hyaluronidases were found in the B. caeruleus venom. We further analyzed protein coverage (Lys/Arg rich and poor regions as well as potential glycosylation sites) using in-house software. These studies expand our understanding of the proteomes of the venoms of these two medically-important species.

Comparative venom gland transcriptomics of Naja kaouthia (monocled cobra) from Malaysia and Thailand: elucidating geographical venom variation and insights into sequence novelty

Background

The monocled cobra (Naja kaouthia) is a medically important venomous snake in Southeast Asia. Its venom has been shown to vary geographically in relation to venom composition and neurotoxic activity, indicating vast diversity of the toxin genes within the species. To investigate the polygenic trait of the venom and its locale-specific variation, we profiled and compared the venom gland transcriptomes of N. kaouthia from Malaysia (NK-M) and Thailand (NK-T) applying next-generation sequencing (NGS) technology.

Methods

The transcriptomes were sequenced on the Illumina HiSeq platform, assembled and followed by transcript clustering and annotations for gene expression and function. Pairwise or multiple sequence alignments were conducted on the toxin genes expressed. Substitution rates were studied for the major toxins co-expressed in NK-M and NK-T.

Results and discussion

The toxin transcripts showed high redundancy (41–82% of the total mRNA expression) and comprised 23 gene families expressed in NK-M and NK-T, respectively (22 gene families were co-expressed). Among the venom genes, three-finger toxins (3FTxs) predominated in the expression, with multiple sequences noted. Comparative analysis and selection study revealed that 3FTxs are genetically conserved between the geographical specimens whilst demonstrating distinct differential expression patterns, implying gene up-regulation for selected principal toxins, or alternatively, enhanced transcript degradation or lack of transcription of certain traits. One of the striking features that elucidates the inter-geographical venom variation is the up-regulation of α-neurotoxins (constitutes ∼80.0% of toxin’s fragments per kilobase of exon model per million mapped reads (FPKM)), particularly the long-chain α-elapitoxin-Nk2a (48.3%) in NK-T but only 1.7% was noted in NK-M. Instead, short neurotoxin isoforms were up-regulated in NK-M (46.4%). Another distinct transcriptional pattern observed is the exclusively and abundantly expressed cytotoxin CTX-3 in NK-T. The findings suggested correlation with the geographical variation in proteome and toxicity of the venom, and support the call for optimising antivenom production and use in the region. Besides, the current study uncovered full and partial sequences of numerous toxin genes from N. kaouthia which have not been reported hitherto; these include N. kaouthia-specific l-amino acid oxidase (LAAO), snake venom serine protease (SVSP), cystatin, acetylcholinesterase (AChE), hyaluronidase (HYA), waprin, phospholipase B (PLB), aminopeptidase (AP), neprilysin, etc. Taken together, the findings further enrich the snake toxin database and provide deeper insights into the genetic diversity of cobra venom toxins.

Peptide Fraction pOh2 Exerts Antiadipogenic Activity through Inhibition of C/EBP-α and PPAR-γ Expression in 3T3-L1 Adipocytes

Many studies have comprehensively examined the venom of Ophiophagus hannah snake. Its venom comprises different compounds exhibiting a wide range of pharmacological activities. In this investigation, four peptide fractions (PFs), ranging from 3 kDa to 10 kDa, isolated from the Vietnamese snake venom of O. hannah were separated by HPLC and investigated for their inhibitory activity on adipogenesis in 3T3-L1 adipocytes. The most effective PF was then further purified, generating two peptides, pOh1 and pOh2. Upon investigation of these two peptides on 3T3-L1 adipocytes, it was revealed that, at 10 μg/mL, pOh2 was able to inhibit the lipid accumulation in 3T3-L1 adipocytes by up to 56%, without affecting cell viability. Furthermore, the pOh2 downregulated the gene expression of important transcription factors C/EBP-α and PPAR-γ. In addition, aP2 and GPDH adipocyte-specific markers were also significantly reduced compared to untreated differentiated cells. Taken together, pOh2 inhibited the expression of key transcription factors C/EBP-α and PPAR-γ and their target genes, aP2 and GPDH, thereby blocking the adipocyte differentiation. In conclusion, this novel class of peptide might have potential for in vivo antiobesity effects.

Proteomic analysis to unravel the complex venom proteome of eastern India Naja naja: Correlation of venom composition with its biochemical and pharmacological properties

The complex venom proteome of the eastern India (EI) spectacled cobra (Naja naja) was analyzed using tandem mass spectrometry of cation-exchange venom fractions. About 75% of EI N. naja venom proteins were <18kDa and cationic at physiological pH of blood. SDS-PAGE (non-reduced) analysis indicated that in the native state venom proteins either interacted with each-other or self-aggregated resulting in the formation of higher molecular mass complexes. Proteomic analysis revealed that 43 enzymatic and non-enzymatic proteins in EI N. naja venom with a percent composition of about 28.4% and 71.6% respectively were distributed over 15 venom protein families. The three finger toxins (63.8%) and phospholipase A₂s (11.4%) were the most abundant families of non-enzymatic and enzymatic proteins, respectively. nanoLC-ESI-MS/MS analysis demonstrated the occurrence of acetylcholinesterase, phosphodiesterase, cholinesterase and snake venom serine proteases in N. naja venom previously not detected by proteomic analysis. ATPase, ADPase, hyaluronidase, TAME, and BAEE-esterase activities were detected by biochemical analysis; however, due to a limitation in the protein database depository they were not identified in EI N. naja venom by proteomic analysis. The proteome composition of EI N. naja venom was well correlated with its in vitro and in vivo pharmacological properties in experimental animals and envenomed human.

BIOLOGICAL SIGNIFICANCE

Proteomic analysis reveals the complex and diverse protein profile of EI N. naja venom which collectively contributes to the severe pathophysiological manifestation upon cobra envenomation. The study has also aided in comprehending the compositional variation in venom proteins of N. naja within the Indian sub-continent. In addition, this study has also identified several enzymes in EI N. naja venom which were previously uncharacterized by proteomic analysis of Naja venom.

Mapping Proteoforms and Protein Complexes From King Cobra Venom Using Both Denaturing and Native Top-down Proteomics

Characterizing whole proteins by top-down proteomics avoids a step of inference encountered in the dominant bottom-up methodology when peptides are assembled computationally into proteins for identification. The direct interrogation of whole proteins and protein complexes from the venom of Ophiophagus hannah (king cobra) provides a sharply clarified view of toxin sequence variation, transit peptide cleavage sites and post-translational modifications (PTMs) likely critical for venom lethality. A tube-gel format for electrophoresis (called GELFrEE) and solution isoelectric focusing were used for protein fractionation prior to LC-MS/MS analysis resulting in 131 protein identifications (18 more than bottom-up) and a total of 184 proteoforms characterized from 14 protein toxin families. Operating both GELFrEE and mass spectrometry to preserve non-covalent interactions generated detailed information about two of the largest venom glycoprotein complexes: the homodimeric l-amino acid oxidase (∼130 kDa) and the multichain toxin cobra venom factor (∼147 kDa). The l-amino acid oxidase complex exhibited two clusters of multiproteoform complexes corresponding to the presence of 5 or 6 N-glycans moieties, each consistent with a distribution of N-acetyl hexosamines. Employing top-down proteomics in both native and denaturing modes provides unprecedented characterization of venom proteoforms and their complexes. A precise molecular inventory of venom proteins will propel the study of snake toxin variation and the targeted development of new antivenoms or other biotherapeutics.

Functional characterizations of venom phenotypes in the eastern diamondback rattlesnake (Crotalus adamanteus) and evidence for expression-driven divergence in toxic activities among populations

Phenotypes frequently vary across and within species. The connection between specific phenotypic effects and function, however, is less understood despite being essential to our understanding of the adaptive process. Snake venoms are ideal for identifying functionally important phenotypic variation because venom variation is common, and venoms can be functionally characterized through simple assays and toxicity measurements. Previous work with the eastern diamondback rattlesnake (Crotalus adamanteus) used multivariate statistical approaches to identify six unique venom phenotypes. We functionally characterized hemolytic, gelatinase, fibrinogenolytic, and coagulant activity for all six phenotypes, as well as one additional venom, to determine if the statistically significant differences in toxin expression levels previously documented corresponded to differences in venom activity. In general, statistical differences in toxin expression predicted the identified functional differences, or lack thereof, in toxic activity, demonstrating that the statistical approach used to characterize C. adamanteus venoms was a fair representation of biologically meaningful differences. Minor differences in activity not accounted for by the statistical model may be the result of amino-acid differences and/or post-translational modifications, but overall we were able to link variation in protein expression levels to variation in function as predicted by multivariate statistical approaches.

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

Stingrays commonly cause human envenoming related accidents in populations of the sea, near rivers and lakes. Transcriptomic profiles have been used to elucidate components of animal venom, since they are capable of providing molecular information on the biology of the animal and could have biomedical applications. In this study, we elucidated the transcriptomic profile of the venom glands from two different freshwater stingray species that are endemic to the Paraná-Paraguay basin in Brazil, Potamotrygon amandae and Potamotrygon falkneri. Using RNA-Seq, we identified species-specific transcripts and overlapping proteins in the venom gland of both species. Among the transcripts related with envenoming, high abundance of hyaluronidases was observed in both species. In addition, we built three-dimensional homology models based on several venom transcripts identified. Our study represents a significant improvement in the information about the venoms employed by these two species and their molecular characteristics. Moreover, the information generated by our group helps in a better understanding of the biology of freshwater cartilaginous fishes and offers clues for the development of clinical treatments for stingray envenoming in Brazil and around the world. Finally, our results might have biomedical implications in developing treatments for complex diseases.

Unveiling the elusive and exotic: Venomics of the Malayan blue coral snake (Calliophis bivirgata flaviceps)

The venom proteome of the Malayan blue coral snake, Calliophis bivirgata flaviceps from west Malaysia was investigated by 1D-SDS-PAGE and shotgun-LCMS/MS. A total of 23 proteins belonging to 11 protein families were detected from the venom proteome. For the toxin proteins, the venom consists mainly of phospholipase A2 (41.1%), cytotoxin (22.6%), SVMPs (18.7%) and vespryns (14.6%). However, in contrast to the venoms of New World coral snakes and most elapids, there was no post-synaptic α-neurotoxin detected. The proteome also revealed a relatively high level of phosphodiesterase (1.3%), which may be associated with the reported high level of adenosine in the venom. Also detected were 5'-nucleotidase (0.3%), hyaluronidase (0.1%) and cysteine-type endopeptide inhibitor (0.6%). Enzymatic studies confirmed the presence of phospholipase A2, phosphodiesterase, 5'-nucleotidase and acetylcholinesterase activities but not l-amino acid oxidase activity. The venom exhibited moderate cytotoxic activity against CRL-2648 fibroblast cell lines (IC50=62.14±0.87 μg/mL) and myotoxicity in mice, presumably due to the action of its cytotoxin or its synergistic action with phospholipase A2. Interestingly, the venom lethality could be cross-neutralized by a neurotoxic bivalent antivenom from Taiwan. Together, the findings provide insights into the composition and functions of the venom of this exotic oriental elapid snake.

BIOLOGICAL SIGNIFICANCE

While venoms of the New World coral snake have been extensively studied, literature pertaining to the Old World or Asiatic coral snake venoms remains lacking. This could be partly due to the inaccessibility to the venom of this rare species and infrequent cases of envenomation reported. This study identified and profiled the venom proteome of the Malayan blue coral snake (C. b. flaviceps) through SDS-PAGE and a high-resolution nano-LCMS/MS method, detailing the types and abundance of proteins found in the venom. The biological and toxic activities of the venom were also investigated, offering functional correlation to the venom proteome studied. Of note, the venom contains a unique toxin profile predominated with phospholipase A2 and cytotoxin with no detectable post-synaptic neurotoxin. The venom is moderately lethal to mice and the fatal effect could be cross-neutralized by a heterologous elapid bivalent antivenom from Taiwan. The findings enrich snake toxin databases and provide insights into the composition and pathogenesis of the venom of this exotic species.

Venom-gland transcriptome and venom proteome of the Malaysian king cobra (Ophiophagus hannah)

Background

The king cobra (Ophiophagus hannah) is widely distributed throughout many parts of Asia. This study aims to investigate the complexity of Malaysian Ophiophagus hannah (MOh) venom for a better understanding of king cobra venom variation and its envenoming pathophysiology. The venom gland transcriptome was investigated using the Illumina HiSeq™ platform, while the venom proteome was profiled by 1D-SDS-PAGE-nano-ESI-LCMS/MS.

Results

Transcriptomic results reveal high redundancy of toxin transcripts (3357.36 FPKM/transcript) despite small cluster numbers, implying gene duplication and diversification within restricted protein families. Among the 23 toxin families identified, three-finger toxins (3FTxs) and snake-venom metalloproteases (SVMPs) have the most diverse isoforms. These 2 toxin families are also the most abundantly transcribed, followed in descending order by phospholipases A₂ (PLA₂s), cysteine-rich secretory proteins (CRISPs), Kunitz-type inhibitors (KUNs), and L-amino acid oxidases (LAAOs). Seventeen toxin families exhibited low mRNA expression, including hyaluronidase, DPP-IV and 5’-nucleotidase that were not previously reported in the venom-gland transcriptome of a Balinese O. hannah. On the other hand, the MOh proteome includes 3FTxs, the most abundantly expressed proteins in the venom (43 % toxin sbundance). Within this toxin family, there are 6 long-chain, 5 short-chain and 2 non-conventional 3FTx. Neurotoxins comprise the major 3FTxs in the MOh venom, consistent with rapid neuromuscular paralysis reported in systemic envenoming. The presence of toxic enzymes such as LAAOs, SVMPs and PLA₂ would explain tissue inflammation and necrotising destruction in local envenoming. Dissimilarities in the subtypes and sequences between the neurotoxins of MOh and Naja kaouthia (monocled cobra) are in agreement with the poor cross-neutralization activity of N. kaouthia antivenom used against MOh venom. Besides, the presence of cobra venom factor, nerve growth factors, phosphodiesterase, 5’-nucleotidase, and DPP-IV in the venom proteome suggests its probable hypotensive action in subduing prey.

Conclusion

This study reports the diversity and abundance of toxins in the venom of the Malaysian king cobra (MOh). The results correlate with the pathophysiological actions of MOh venom, and dispute the use of Naja cobra antivenoms to treat MOh envenomation. The findings also provide a deeper insight into venom variations due to geography, which is crucial for the development of a useful pan-regional antivenom.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1828-2) contains supplementary material, which is available to authorized users.

Venomics, lethality and neutralization of Naja kaouthia (monocled cobra) venoms from three different geographical regions of Southeast Asia

Previous studies showed that venoms of the monocled cobra, Naja kaouthia from Thailand and Malaysia are substantially different in their median lethal doses. The intraspecific venom variations of N. kaouthia, however, have not been fully elucidated. Here we investigated the venom proteomes of N. kaouthia from Malaysia (NK-M), Thailand (NK-T) and Vietnam (NK-V) through reverse-phase HPLC, SDS-PAGE and tandem mass spectrometry. The venom proteins comprise 13 toxin families, with three-finger toxins being the most abundant (63-77%) and the most varied (11-18 isoforms) among the three populations. NK-T has the highest content of neurotoxins (50%, predominantly long neurotoxins), followed by NK-V (29%, predominantly weak neurotoxins and some short neurotoxins), while NK-M has the least (18%, some weak neurotoxins but less short and long neurotoxins). On the other hand, cytotoxins constitute the main bulk of toxins in NK-M and NK-V venoms (up to 45% each), but less in NK-T venom (27%). The three venoms show different lethal potencies that generally reflect the proteomic findings. Despite the proteomic variations, the use of Thai monovalent and Neuro polyvalent antivenoms for N. kaouthia envenomation in the three regions is appropriate as the different venoms were neutralized by the antivenoms albeit at different degrees of effectiveness.

BIOLOGICAL SIGNIFICANCE

Biogeographical variations were observed in the venom proteome of monocled cobra (Naja kaouthia) from Malaysia, Thailand and Vietnam. The Thai N. kaouthia venom is particularly rich in long neurotoxins, while the Malaysian and Vietnamese specimens were predominated with cytotoxins. The differentially expressed toxin profile accounts for the discrepancy in the lethal dose of the venom from different populations. Commercially available Thai antivenoms (monovalent and polyvalent) were able to neutralize the three venoms at different effective doses, hence supporting their uses in the three regions. While dose adjustment according to geographical region seems possible, changes to standard recommended dosage should only be made if further study validates that the monocled cobras within a population do not exhibit remarkable inter-individual venom variation.

Proteomic characterization and comparison of Malaysian Bungarus candidus and Bungarus fasciatus venoms

Kraits (Bungarus spp.) are highly venomous elapids that are only found in Asia. In the current study, 103 and 86 different proteins were identified from Bungarus candidus and Bungarus fasciatus venoms, respectively. These proteins were classified into 18 different venom protein families. Both venoms were found to contain a high percentage of three finger toxins, phospholipase A2 enzymes and Kunitz-type inhibitors. Smaller number of high molecular weight enzymes such as L-amino acid oxidase, hyaluronidases, and acetylcholinesterase were also detected in the venoms. We also detected some unique proteins that were not known to be present in these venoms. The presence of a natriuretic peptide, vespryn, and serine protease families was detected in B. candidus venom. We also detected the presence of subunit A and B of β-bungarotoxin and α-bungarotoxin which had not been previously found in B. fasciatus venom. Understanding the proteome composition of Malaysian krait species will provide useful information on unique toxins and proteins which are present in the venoms. This knowledge will assist in the management of krait envenoming. In addition, these proteins may have potential use as research tools or as drug-design templates.

BIOLOGICAL SIGNIFICANCE

This study has revealed the proteome composition of Malaysian B. candidus and B. fasciatus venoms, two medically important snake species in Asia. Information on the venom proteome of these species will provide useful information for krait bite management and aid in antivenom selection. Venom proteome profiles of these venoms showed that there are significant differences in the venom protein family compositions. Detection of proteins and peptides that have not been documented in these species such as natriuretic peptides, vespryn and serine proteases provides new knowledge on the composition of these venoms. The roles of these new proteins and peptides in krait envenoming are still unknown. Discovery of these proteins and peptides may also be useful for future research tool and therapeutic development.

Ophiophagus hannah venom: proteome, components bound by Naja kaouthia antivenin and neutralization by N. kaouthia neurotoxin-specific human ScFv

Venomous snakebites are an important health problem in tropical and subtropical countries. King cobra (Ophiophagushannah) is the largest venomous snake found in South and Southeast Asia. In this study, the O. hannah venom proteome and the venom components cross-reactive to N. kaouthia monospecific antivenin were studied. O. hannah venom consisted of 14 different protein families, including three finger toxins, phospholipases, cysteine-rich secretory proteins, cobra venom factor, muscarinic toxin, L-amino acid oxidase, hypothetical proteins, low cysteine protein, phosphodiesterase, proteases, vespryn toxin, Kunitz, growth factor activators and others (coagulation factor, endonuclease, 5’-nucleotidase). N. kaouthia antivenin recognized several functionally different O. hannah venom proteins and mediated paratherapeutic efficacy by rescuing the O. hannah envenomed mice from lethality. An engineered human ScFv specific to N. kaouthia long neurotoxin (NkLN-HuScFv) cross-neutralized the O. hannah venom and extricated the O. hannah envenomed mice from death in a dose escalation manner. Homology modeling and molecular docking revealed that NkLN-HuScFv interacted with residues in loops 2 and 3 of the neurotoxins of both snake species, which are important for neuronal acetylcholine receptor binding. The data of this study are useful for snakebite treatment when and where the polyspecific antivenin is not available. Because the supply of horse-derived antivenin is limited and the preparation may cause some adverse effects in recipients, a cocktail of recombinant human ScFvs for various toxic venom components shared by different venomous snakes, exemplified by the in vitro produced NkLN-HuScFv in this study, should contribute to a possible future route for an improved alternative to the antivenins.

Proteomic characterization of venom of the medically important Southeast Asian Naja sumatrana (Equatorial spitting cobra)

The proteome of Naja sumatrana (Equatorial spitting cobra) venom was investigated by shotgun analysis and a combination of ion-exchange chromatography and reverse phase HPLC. Shotgun analysis revealed the presence of 39 proteins in the venom while the chromatographic approach identified 37 venom proteins. The results indicated that, like other Asiatic cobra venoms, N. sumatrana contains large number of three finger toxins and phospholipases A2, which together constitute 92.1% by weight of venom protein. However, only eight of the toxins can be considered as major venom toxins. These include two phospholipases A2, three neurotoxins (two long neurotoxins and a short neurotoxin) and three cardiotoxins. The eight major toxins have relative abundance of 1.6-27.2% venom proteins and together account for 89.8% (by weight) of total venom protein. Other venom proteins identified include Zn-metalloproteinase-disintegrin, Thaicobrin, CRISP, natriuretic peptide, complement depleting factors, cobra venom factors, venom nerve growth factor and cobra serum albumin. The proteome of N. sumatrana venom is similar to proteome of other Asiatic cobra venoms but differs from that of African spitting cobra venom. Our results confirm that the main toxic action of N. sumatrana venom is neurotoxic but the large amount of cardiotoxins and phospholipases A2 are likely to contribute significantly to the overall pathophysiological action of the venom. The differences in toxin distribution between N. sumatrana venom and African spitting cobra venoms suggest possible differences in the pathophysiological actions of N. sumatrana venom and the African spitting cobra venoms, and explain why antivenom raised against Asiatic cobra venom is not effective against African spitting cobra venoms.

Intracellular antibody immunity and the cytosolic Fc receptor TRIM21

Until recently, it was thought that antibody effector mechanisms were mediated purely by Fc receptors expressed on professional cells, following capture of immune complexes in the extracellular space. Recently a new Fc receptor, TRIM21, was discovered that is expressed by cells of all histogenetic lineages and which mediates immune responses intracellularly. This new receptor possesses many unique structural and functional properties. TRIM21 binds both IgG and IgM, interacts primarily with the CH3 rather than CH2 domain and engages two heavy chains simultaneously. This latter property allows TRIM21 to bind antibodies with a higher affinity than any other Fc receptor. TRIM21 is cytosolic, has both effector and signalling functions and is exquisitely conserved in mammals. The discovery of this missing part of humoral immunity has important implications for where and how antibodies work.

Characterization of inflamin, the first member of a new family of snake venom proteins that induces inflammation

Unlike other sea snakes, Aipysurus eydouxii feeds exclusively on fish eggs. This unusual feeding habit prompted us to search for unique transcripts in their venom glands. In the present study we expressed a novel cysteine-rich secretory protein containing 94 amino acid residues that was identified in its cDNA library. As it induced inflammation and writhing in animals, this protein was named inflamin. It induced two waves of prostanoid production. The first wave peaked at 10 min and 6-oxo PGF1α (prostaglandin F1α) (6-keto PGF1α) was the major product. The second wave, specifically of 6-oxo PGF1α and PGE2 (prostanglandin E2), started after 2 h. In RAW 264.7 cells, COX-1 (cyclo-oxygenase-1) activity showed a transient increase at 10 min and is responsible for the first wave, but its expression was unaffected. COX-2 was induced after 3 h and is responsible for the second wave. Using specific inhibitors, we showed that cPLA2 (calcium-dependent phospholipase A2), and not sPLA2 (secretory phospholipase A2), iPLA2 (calcium-independent phospholipase A2) or DAG (diacylglycerol) lipase, plays a key role in arachidonate release. The cPLA2 activity showed a transient increase of 62% at 10 min; this increase was due to its phosphorylation and not due to an increase in its expression. Thus inflamin, the first member of a new family of snake venom proteins, leads to an increase in the cPLA2 and COX-1 activity resulting in inflammation and pain.