Next Article in Journal
Dissecting Toxicity: The Venom Gland Transcriptome and the Venom Proteome of the Highly Venomous Scorpion Centruroides limpidus (Karsch, 1879)
Previous Article in Journal
Toxicity and Toxin Composition of the Greater Blue-Ringed Octopus Hapalochlaena lunulata from Ishigaki Island, Okinawa Prefecture, Japan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Proteomic Analysis of Human Blister Fluids Following Envenomation by Three Snake Species in India: Differential Markers for Venom Mechanisms of Action

1
Department of Microbiology, Immunology and Cancer Biology, University of Virginia, School of Medicine, Charlottesville, VA 22908-734, USA
2
Department of Cell Biology, Brazilian Center of Protein Research, University of Brasilia, Brasilia/DF 70297-400, Brazil
3
Little Flower Hospital, Angamaly 683572, India
4
Sree Naryana Institute of Medical Science, Kerala 683594, India
5
Instituto Clodomiro Picado, School of Microbiology, University of Costa Rica, San José 11501-2060, Costa Rica
*
Author to whom correspondence should be addressed.
Submission received: 11 April 2019 / Revised: 25 April 2019 / Accepted: 27 April 2019 / Published: 30 April 2019
(This article belongs to the Section Animal Venoms)

Abstract

:
Skin blistering as a result of snakebite envenomation is characteristic of some bites, however little is known regarding the mechanism of blister formation or the composition of the blister fluid. In order to investigate if blister fluid proteomes from humans suffering snakebite envenomation could provide insights on the pathophysiology of these skin alterations, blister fluid was collected from six patients upon presentation at a clinic in India bitten by three species of snakes, Daboia russelii (3), Hypnale hypnale (2), or Naja naja (1). Standard clinical data were recorded throughout the treatment. Approximately 805 proteins were identified in blister fluids using proteomic analyses. Informatics analyses of the proteomes identified the top biological response categories as: platelet degranulation, innate immune response, receptor-mediated endocytosis, complement activation, and blood coagulation. Hierarchical clustering did not show a clear segregation of patients’ proteomes being associated with the species of snake involved, suggesting that either the proteomic profiles described reflect a general response to venom-induced tissue damage or more patient data sets will be required to observe significant differences. Finally, it is of interest that venom proteins were also identified in the blister fluids suggesting that this fluid may serve as a reservoir of venom biologically active proteins/toxins, and as such, may indicate the clinical value of removing blister fluid to attenuate further tissue damage.
Key Contribution: The proteomic analysis of blister fluid collected from people suffering snakebite envenomation is a valuable tool to understand the pathophysiology of envenomation; to identify venom components; and to provide insights into new therapeutic options.

1. Introduction

Snakebite envenomations represent a significant public health problem observed on a global basis, but particularly in sub-Saharan Africa, Asia, Latin America, and parts of Oceania, and have been included by the World Health Organization (WHO) in the list of neglected tropical diseases [1]. India is the country that bears the highest number of cases and it has been estimated that 50,000 people die each year in this country due to snakebites [2].
Among the many species of venomous snakes that inhabit India, four have been considered as causing the heaviest toll: the Indian cobra (Naja naja), the common krait (Bungarus caeruleus), the Russell’s viper (Daboia russelii), and the saw-scaled viper (Echis carinatus); these species comprise the group known as “the big four” [3]. In addition, it has been observed that, in certain regions, the hump-nosed pit viper species (Hypnale hypnale) inflicts a high number of bites and is capable of causing serious accidents, underscoring the need to further study the epidemiology of snakebites in India [4]. Although polyvalent snake antivenoms (SAV) are produced by several manufacturers in India against venoms of the “big four” species, envenomations continue to cause fatalities and pathophysiological sequelae in thousands of people a year. Thus, there is a need to further expand our knowledge on the various aspects of snakebite envenomation in this country, including a better understanding of the pathophysiology and mechanisms of action of venom-induced effects. A rapid, personalized and highly sensitive diagnosis of envenomation and identification of the venom is one of the essential needs for effective injury management. Thus, efforts have been directed at the identification of markers of tissue damage, especially at the injury site, since they may provide valuable information about the mechanisms of disease and the condition and prognosis of the patient, all of which may lead to a more effective treatment [5,6,7]
One of the characteristic manifestations of the local tissue damage inflicted by venoms of snakes from the family Viperidae, and by some species of the family Elapidae, is the formation of blisters, resulting in the accumulation of a proteinaceous fluid as a consequence of the collection of wound exudate [1]. Thereby, the study of wound exudate and blister fluids as a potential source of biomarkers have emerged as a new approach to study the effects of venoms and may assist in the diagnosis and treatment of these patients [5,8]. The proteomic analysis of fluids collected from injured tissue, and from tissue undergoing healing processes, has allowed the identification of markers that differentiate between healing and non-healing wounds, with mediators characteristic of tissue formation or mediators characteristic for a persistent inflammatory and tissue damaging response, respectively [7]. Moreover, immunochemical analysis of blister fluid from an envenomated patient has been used to identify the presence of venom components in the fluid [9].
Since viperid snake venoms induce the formation of blisters and an inflammatory exudate in the damaged tissues, a number of experimental studies have addressed the analysis of the proteomes of exudates collected from mice injected with venoms or isolated toxins. A comparative proteomic analysis of wound exudates caused by BaP1, a snake venom metalloproteinase (SVMP), and Mtx-I, a myotoxic phospholipase A2 (PLA2), showed that the SVMP causes degradation of nonfibrillar collagens, whereas PLA2-mediated tissue damage results in the presence of fibrillar collagen type I fragments; apolipoproteins A-I, A-IV, and E; and fibronectin in the exudates [5]. Moreover, the use of inhibitors demonstrated that myotoxic PLA2s and hemorrhagic SVMPs are involved in blister formation and skin damage induced by Bothrops asper venom [10]. Furthermore, the presence of matrix metalloproteinases (MMPs) in these exudates was observed, hence suggesting the combined action of endogenous and exogenous proteinases in these pathological events [11,12].
In a recent study, the exudates generated by SVMPs of the PI, PII, and PIII classes revealed specific proteins associated with the mechanism of action of the SVMPs [13]. Exudates from tissue injected with PI and PII SVMPs contained a higher abundance of keratins as compared to samples from mice injected with a PIII SVMP. In agreement, a PI SVMP has been shown to induce blistering in the skin of mice [14]. Hence the study of the composition of exudates and blister fluid has proven to be an information-rich approach to assess pathological processes taking place in tissues injected with snake venoms and toxins.
Despite these advances at the experimental level, to our knowledge, there have been no attempts to investigate the composition of blister fluid or exudates in humans resulting from snakebite envenomations. The proteome analysis of fluids collected from affected tissue is therefore an experimental approach with great potential for biomarker identification by allowing the assessment of tissue alterations associated with the pathology of envenomation. In this study, utilizing mass spectrometry, we investigated the composition of blister fluid collected from patients envenomated by three species of Indian snakes, i.e., Daboia russelii and Hypnale sp. (family Viperidae), or Naja naja (family Elapidae).

2. Results and Discussion

2.1. Clinical Features of Envenomation

The basic demographic and clinical information of the six patients included in this study are shown in Table 1. Three of them were bitten by D. russelii, two by Hypnale sp, and one by Naja naja. Ages ranged from 2 to 67 years old. The time lapse between the bite and the arrival to the hospital ranged between 2 and 8.5 h. Envenomation by D. russelii was characterized by local manifestations (swelling, blisters), and one of them presented evidence of systemic capillary leakage syndrome (parotid swelling, conjunctival chemosis) and thrombocytopenia. Three patients developed unclottable blood, as evidenced by prolongation of the 20 min whole blood clotting test, and acute kidney injury (AKI) was present in two of them, as evidenced by elevated serum creatinine concentration (Table 1).
The two patients suffering bites by Hypnale sp. presented local manifestations (local pain, swelling and discoloration in the skin), and one of them indicated abdominal pain. One patient showed a prolonged 20 min whole blood clotting test, thrombocytopenia, and AKI, with elevated serum creatinine concentration (Table 1). The patient suffering from a cobra (N. naja) bite presented with local alterations (swelling and bluish discolouration) and manifestations of systemic neurotoxicity (ptosis and external ophtalmoplegia).
Patients suffering envenomations by D. russelii and N. naja received 20 vials (200 mL total) of polyspecific antivenom, whereas no antivenom was administered to the patients bitten by Hypnale sp. since this antivenom did not include Hypnale sp. venom in the immunizing mixture. The patient bitten by N. naja was treated with polyspecific antivenom (200 mL) and neostigmine, and those who developed acute kidney injury underwent hemodialysis (Table 1). All patients survived.

2.2. Proteomic Analysis of Blisters

Blister fluid obtained from these patients were first analyzed using SDS-PAGE gel electrophoresis. For some samples, there appeared to be a correlation between the gel profiles and the species of snakes causing the bites. For others, however, there were differences even when comparing blisters obtained from patients bitten by the same species. This was observed in the case of the two patients bitten by Hypnale sp. (Figure 1). After electrophoresis, bands digestion, and LC-MS assay, the proteomic data were obtained and various analyses were performed. Mass spectrometry results identified a total of 805 different proteins in blisters collected from the patients (Supplementary Table S1) Many of the proteins detected in this analysis have also been described in laboratory studies collecting exudate fluid from the tissue of mice injected with the venom of the Central American pit viper Bothrops asper or with purified hemorrhagic SVMPs or myotoxic PLA2s isolated from this and other viperid venoms [5,13,15]. The individual protein identification was considered significant with a minimum of two peptides detected and protein identification probability above 95%.
The most abundant proteins based on their quantitative value are presented in Table 2 and the 15 most abundant proteins in each exudate are highlighted in gray. Informatics analysis of the proteomes identified the top biological response categories as: platelet degranulation, negative regulation of endopeptidase activity, innate immune response, proteolysis, receptor-mediated endocytosis, complement activation, extracellular matrix organization, and blood coagulation (Figure 2). Such events illustrate the local and systemic events triggered in the course of envenomation in terms of change in metabolism, secretion, and synthesis of diverse components as a response to the action of snake venom in the tissues.
Proteins were classified in groups according to functional characteristics, using Gene Ontology within the DAVID Gene Functional Classification Tool (http://david.abcc.ncifcrf.gov). From the pathophysiological standpoint, it was of interest to analyze proteins of the coagulation system (Table 3) and from the extracellular matrix (ECM) (Table 4), as coagulopathies and ECM degradation are characteristic of viperid snakebite envenomation. Regarding proteins of the coagulation system, fibrinogen and plasminogen were the most abundant, with differences noted between patients. Blister fluid from patient 3 (bitten by D. russelii), who showed the most severe clinical manifestations of envenomation (Table 1), presented a higher abundance for fibrinogen alpha-chain and plasminogen. Interestingly, the sample from the patient bitten by the cobra N. naja, also contained fibrinogen, although cobra venoms do not cause significant effects on coagulation. This finding is probably a consequence of plasma extravasation following cobra venom-induced local tissue damage.
The analysis of extracellular matrix (ECM) proteins provided interesting findings (Table 4). Samples from patients 3 (D. russelii) and 4 (Hypnale sp.) presented fragments of heparan sulphate proteoglycan core protein and type IV collagen, two key components of basement membranes. Hydrolysis of these proteins has been associated, in experimental studies in mice, with the action of hemorrhagic SVMPs and with the pathogenesis of microvascular damage leading to extravasation [13,18]. The presence of these proteins in blister fluid is indicative of basement membrane degradation, possibly associated with local hemorrhage and blistering. Again, in the case of patient 3 bitten by D. russelii, these findings correlated with the severe manifestations of envenoming. In this patient, tenascin-X and subunit alpha-5 of laminin were also detected in the blister fluid, further corroborating the ECM damage induced by D. russelii venom. Fragments of tenascin-X were also found in exudates collected from mice injected with hemorrhagic SVMPs [5,13,18]. In addition, blister fluid from patients 4 (Hypnale sp.) and 1 (Naja naja) contained fragments of type VI collagen. Experimental studies in mice injected with SVMPs have shown similar findings in exudates collected from damaged muscle tissue [13,18]. Type VI collagen is known to connect the basement membrane with fibrillar collagens, hence integrating the basement membrane with the surrounding matrix. Hydrolysis of this collagen may also have implications for the pathogenesis of local tissue damage. Our findings underscore the value of the proteomic analysis of blister fluid to assess ECM degradation in snakebite envenoming.

2.3. Protein Clustering

Hierarchical clustering created from the proteomic profile did not show a clear segregation on the basis of the species of snake causing envenomation (Figure 3). This may have a number of possible explanations including: (a) The set of proteins present in blisters is mostly a consequence of tissue damage and the inflammatory response to tissue pathology, regardless of the venom responsible for the damage. The three species involved in this study are known to induce local tissue damage. (b) The set of proteins found is highly dependent on the severity of the damage and not on the species causing the bite. (c) The number of patient blisters analyzed by a particular snake was insufficient to power a significant observation of differences in response. Lastly, (d) the proteomics of blister fluid likely is also time dependent and thus these samples collected at different times might mask differences due to the species of snake. As indicated in Table 1, there was a wide variation in the severity of envenomation within the patients studied. In agreement, patient 3, who suffered a severe envenomation from D. russelii, presented the highest abundance of some clotting proteins and of ECM proteins, in contrast with the other two patients affected by this species. Whatever the case, it is necessary to study the proteomics of blister fluid of a larger number of patients showing a spectrum of severity and at different times after the bite in order to have a comprehensive view of the changes in the proteomic profile of blister fluid in different circumstances of envenomation.

2.4. Identification of Proteins from the Snake Venoms in Blisters

Venom proteins were identified in the blister fluids using proteomics and Western blotting. Using proteomics, venom proteins of various families were identified in two patients bitten by D. russelii and N. naja (Table 5). In the case of patient 3 (D. russelii), this is in agreement with the clinical manifestations that showed a high severity envenomation (Table 1). On the other hand, detection of venom proteins using immunoblotting performed with Indian polyspecific antivenom was able to identify snake venom proteins in several samples. The three Russell’s viper samples exhibited proteins imaged by Western blot in the high molecular mass range, plus some between 40 and 75 kDa and another of around 25 kDa. The sample from the N. naja case, in turn, presented not only a greater intensity but also a higher number of bands (Figure 4). Blister fluid from the patient bitten by N. naja was therefore the sample with a higher abundance of venom detected using both mass spectrometry and Western blot analyses. The abundance of venom proteins present in the blister fluid in the case of the patient bitten by N. naja may be related to the fact that cobras usually inject their venoms subcutaneously, whereas viperid venoms generally inject venom deeper in the tissues. The lower amount of proteins identified by proteomics in some samples as compared to Western blot may be due to the scarcity of venom sequence databases for mass spectrometry analysis, which has been a challenge for all venom proteomic studies. On the other hand, the low reactivity of some samples in the Western blot might be due to variations in the immunoreactivity of antivenom antibodies against some venom components. Immunochemical detection of venom of the Taiwanese snake Protobothrops mucrosquamatus was described in the blister fluid collected from a patient [9].

3. Conclusions

Previous studies carried out in mice have described the potential of proteomics analysis of wound exudate fluid to provide a “window” to explore the pathological and inflammatory events occurring in tissues injected with snake venoms. The present investigation expanded this type of analysis to the fluid of blisters present in human patients proximal to the anatomical site of venom injection in bites inflicted by three different snake species present in India. The results corroborate the value of proteomics analysis for detecting a variety of tissue plasma and inflammatory proteins of different origin as a consequence of the tissue-damaging action of these venoms and the subsequent tissue response through inflammation and repair. Although some differences were observed between samples from patients bitten by different species, overall there was no clear segregation on the basis of the offending species via hierarchical clustering analysis. Therefore, it would be of value in the future to expand these analyses to a larger group of patients in order to assess whether specific tissue markers may characterize different types of envenomation, and whether different proteomic patterns occur depending on the severity and time-course of envenomation. Notably, venom was detected in the blister fluid, as observed by both mass spectrometric and immunochemical analyses. Overall, our findings underscore the use of proteomics analysis of exudate and blister fluids of patients suffering snakebite envenomations for further advancing the understanding of the pathophysiology of this neglected tropical disease. Furthermore, the presence of venom in blister fluid suggests that it may represent a reservoir of venom for further diffusion into the tissues, hence supporting the concept that aseptic removal of blister fluid may reduce the extent and time-dependence of venom-induced tissue damage.

4. Materials and Methods

4.1. Clinical Evaluation and Treatment

Upon admission to hospital, routine clinical and laboratory evaluations were carried out for all patients. The hospital’s routine treatment protocols were followed, including, when appropriate, the administration of polyspecific Indian antivenom, and various interventions that depended on the clinical evolution of each case. Identification of the offending snake was done via direct examination of the specimen, brought in by the patient or his/her relatives.

4.2. Sample Collection

Fluid was aspirated from the blister under sterile conditions using a syringe. Two to five milliliters of volume was typically collected from a single blister.

4.3. Ethics Statement

The Ethical Committee of the Little Flower Hospital and Research Centre, Algamaly, India, authorized the project (“Blister fluid of snake bite victims analysis,” IRB session held on May 6, 2012). An informed consent was obtained from participating patients.

4.4. Proteomic Analysis of the Exudates

4.4.1. SDS-PAGE

Lyophilized blister samples were kept frozen (−18 °C) until time of use. Protein quantification was performed using a micro BCA protein assay kit (Thermo Scientific, Waltham, WA, USA),. Twenty micrograms of protein were further resuspended in Laemmli buffer, applied to a 12% precast electrophoresis gel (Bio-Rad, Hercules, CA, USA), separated, and stained with 0.1% Coomassie Brilliant Blue in 40% methanol, 10% acetic acid for about 1.5 h, then distained with 40% methanol, 10% acetic acid, then in 5% acetic acid.

4.4.2. Protein Extraction and Digestion

Gel lanes were cut in ten equal-sized slices, which were destained for 3 h, and the proteins reduced with 10 mM dithiothreitol (DTT) and alkylated with 50 mM iodoacetamide at room temperature. Gel slices were washed with 100mM ammonium bicarbonate, dehydrated with acetonitrile, and dried in a speed vac. A solution of Promega modified trypsin (20 ng/μL) in 50mM ammonium bicarbonate was used to rehydratate the samples for 30 min on ice. Excess trypsin solution was removed and the digestion was carried on for an additional 18 h at 37 °C. Tryptic peptides were twice extracted from gel slices with 30 μL of a 50% acetonitrile/5% formic acid solution. The extracts were dried to a volume of 15 μL for mass spectrometric analysis.

4.4.3. Proteomic Analysis

A Thermo Electron Orbitrap Velos ETD mass spectrometer system was used to perform the LC/MS/MS experiments. Analytical columns were made by packing 0.5 cm of irregular C18 Beads (YMC Gel ODS-A, 12 nm, I-10-25 μm) followed by 7.5 cm Jupiter 10 μm C18 packing material (Phenomenex, Torrance, CA, USA) into 360 × 75 μm fused silica (Polymicro Technologies, Phoenix, AZ, USA) behind a bottleneck. Aliquots of 7 μL were loaded directly onto these columns and then eluted into the mass spectrometer at 0.5 μL/min using a 1h gradient consisting of acetonitrile/0.1M acetic acid (2–90% acetonitrile). The instrument was set to a full MS (m/z 300–1600) resolution of 60,000 and programmed to acquire a cycle of one mass spectrum followed by collision-induced dissociation (CID) MS/MS performed in the ion trap on the twenty most abundant ions in a data-dependent mode. Dynamic exclusion was enabled with an exclusion list of 400 masses, duration of 60 seconds, and repeat count of 1. The electrospray voltage was set to 2.4 kV, and the capillary temperature was 265 °C. Peak lists were generated from the raw data against the Uniprot Human and NR database from July 2014 using the Sequest search algorithm in Proteome Discoverer 1.4.1. Analysis of the spectra generated was performed using carbamidomethylation on cysteine as a fixed modification, oxidation of methionine as a variable modification, 10 ppm parent tolerance, and 1 Da fragment tolerance. All hits were required to be fully tryptic.
For analysis of the results and validation of peptide and protein identifications, data obtained were exported to Scaffold (version 4.3.2, Proteome Software Inc., Portland, OR, USA). Protein identifications were filtered using Xcorr cutoff values dependent on charge state (+1 >1.8, +2 >2.2, +3 >2.5, and +4 >3.5). Relative quantization of proteins was accomplished by grouping all data from the 10 gel slices for a particular sample in Scaffold and then displaying the quantitative value. This number gives an average total of non-grouped spectral counts for a protein divided by the total non-grouping spectral counts for the 10 mass spectral runs from the gels slices from each lane (http://www.proteomesoftware.com/) and allows a relative quantitative comparison between specific protein from different samples.
Hierarchical clustering analysis was performed using ClustVis (http://biit.cs.ut.ee/clustvis/) to determine any clusters of samples or proteins. To create heat maps, the list of proteins obtained using proteomic analysis was log transformed and both rows and columns were clustered using Euclidean distance and average linkage. DAVID (https://david.ncifcrf.gov/) was applied to the list of the most abundant proteins in order to identify gene ontology terms and KEGG pathways or over-represented biological processes.

4.5. Western Blot

Detection of venom proteins in human wound exudate was done using Western blotting. Exudate samples were separated under non-reducing conditions on 12% SDS-PAGE gels and were transferred to a nitrocellulose membrane. Immunodetection was performed with polyvalent VINS polyspecific antivenom (VINS Bioproducts Limited, India, batch No. 01AS13100). The reaction was detected using an anti-horse peroxidase antibody and a chemiluminiscent substrate. This antivenom is an F(ab’)2 preparation, produced from the plasma of horses immunized with venoms of the Indian species Naja naja, Bungarus caeruleus, Daboia russelii, and Echis carinatus.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6651/11/5/246/s1, Table S1: Total proteins identified by proteomic analysis. Proteins were identified from the raw data using the Sequest search algorithm Proteome Discovery 1.4.1 against a non-specific database and the functional annotation provided using DAVID Bioinformatics Resources 6.7.

Author Contributions

Conceptualization, J.K.J., J.W.F.; methodology, J.K.A.M., J.K.J., J.M., T.E., A.R., J.M.G., J.W.F.; data analysis, J.K.A.M., J.K.J., J.M., T.E., A.R., J.M.G., J.W.F.; data curation, J.K.A.M., T.E., A.R., J.M.G., J.W.F.; writing—original draft preparation, J.K.A.M., A.R., J.M.G., J.W.F.; writing—review and editing, J.K.A.M., J.K.J., J.M., T.E., A.R., J.M.G., J.W.F.; funding acquisition, A.R., J.W.F.

Funding

This research was funded by Vicerrectoría de Insvestigación (grant number 741-B6-125) and the University of Virginia, School of Medicine.

Acknowledgments

The authors thank the Biomolecular Analysis Facility of the University of Virginia, USA, for the proteomic analyses.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Gutiérrez, J.M.; Calvete, J.J.; Habib, A.G.; Harrison, R.A.; Williams, D.J.; Warrell, D.A. Snakebite envenoming. Nat. Rev. Dis. Prim. 2017, 3, 17063. [Google Scholar] [CrossRef] [PubMed]
  2. Mohapatra, B.; Warrell, D.A.; Suraweera, W.; Bhatia, P.; Dhingra, N.; Jotkar, R.M.; Rodriguez, P.S.; Mishra, K.; Whitaker, R.; Jha, P. Snakebite mortality in India: A nationally representative mortality survey. PLoS Negl. Trop. Dis. 2011, 5, 1–8. [Google Scholar] [CrossRef] [PubMed]
  3. Warrell, D.A.; Gutiérrez, J.M.; Calvete, J.J.; Williams, D. New approaches & technologies of venomics to meet the challenge of human envenoming by snakebites in India. Indian J. Med. Res. 2013, 138, 38–59. [Google Scholar] [PubMed]
  4. Simpson, I.D.; Norris, R.L. Snakes of medical importance in India: Is the concept of the “Big 4” still relevant and useful? Wilderness Environ. Med. 2007, 18, 2–9. [Google Scholar] [CrossRef] [PubMed]
  5. Escalante, T.; Rucavado, A.; Pinto, A.F.M.; Terra, R.M.S.; Gutiérrez, J.M.; Fox, J.W. Wound exudate as a proteomic window to reveal different mechanisms of tissue damage by snake venom toxins. J. Proteome Res. 2009, 8, 5120–5131. [Google Scholar] [CrossRef] [PubMed]
  6. Kool, J.; Reubsaet, L.; Wesseldijk, F.; Maravilha, R.T.; Pinkse, M.W.; D’Santos, C.S.; Van Hilten, J.J.; Zijlstra, F.J.; Heck, A.J.R. Suction blister fluid as potential body fluid for biomarker proteins. Proteomics 2007, 7, 3638–3650. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Eming, S.A.; Koch, M.; Krieger, A.; Brachvogel, B.; Kreft, S.; Bruckner-Tuderman, L.; Krieg, T.; Shannon, J.D.; Fox, J.W. Differential Proteomic Analysis Distinguishes Tissue Repair Biomarker Signatures in Wound Exudates Obtained from Normal Healing and Chronic Wounds. J. Proteome Res. 2010, 9, 4758–4766. [Google Scholar] [CrossRef] [PubMed]
  8. Herrera, C.; Escalante, T.; Rucavado, A.; Fox, J.W.; Gutiérrez, J.M. Metalloproteinases in disease: Identification of biomarkers of tissue damage through proteomics. Expert Rev. Proteomics 2018, 15, 967–982. [Google Scholar] [CrossRef] [PubMed]
  9. Lin, C.-C.; Wang, P.-J.; Liu, C.-C. Venom concentrations in blisters and hemorrhagic bullae in a patient bitten by a Taiwan habu (Protobothrops mucrosquamatus). Rev. Soc. Bras. Med. Trop. 2019, 52, e20180160. [Google Scholar] [CrossRef] [PubMed]
  10. Rucavado, A.; Escalante, T.; Shannon, J.; Gutiérrez, J.M.; Fox, J.W. Proteomics of Wound Exudate in Snake Venom-Induced Pathology: Search for Biomarkers To Assess Tissue Damage and Therapeutic Success. J. Proteome Res. 2011, 10, 1987–2005. [Google Scholar] [CrossRef] [PubMed]
  11. Rucavado, A.; Núñez, J.; Gutiérrez, J.M. Blister formation and skin damage induced by BaP1, a haemorrhagic metalloproteinase from the venom of the snake Bothrops asper. Int. J. Exp. Pathol. 1998, 79, 245–254. [Google Scholar] [PubMed]
  12. Herrera, C.; Macêdo, J.K.A.; Feoli, A.; Escalante, T.; Rucavado, A.; Gutiérrez, J.M.; Fox, J.W. Muscle Tissue Damage Induced by the Venom of Bothrops asper: Identification of Early and Late Pathological Events through Proteomic Analysis. PLoS Negl. Trop. Dis. 2016, 10, e0004599. [Google Scholar] [CrossRef] [PubMed]
  13. Herrera, C.; Escalante, T.; Voisin, M.B.; Rucavado, A.; Morazán, D.; Macêdo, J.K.A.; Calvete, J.J.; Sanz, L.; Nourshargh, S.; Gutiérrez, J.M.; Fox, J.W. Tissue Localization and Extracellular Matrix Degradation by PI, PII and PIII Snake Venom Metalloproteinases: Clues on the Mechanisms of Venom-Induced Hemorrhage. PLoS Negl. Trop. Dis. 2015, 9, 1–20. [Google Scholar] [CrossRef] [PubMed]
  14. Jiménez, N.; Escalante, T.; Gutiérrez, J.M.; Rucavado, A. Skin Pathology Induced by Snake Venom Metalloproteinase: Acute Damage, Revascularization, and Re-epithelization in a Mouse Ear Model. J. Invest. Dermatol. 2008, 128, 2421–2428. [Google Scholar] [CrossRef] [PubMed]
  15. Rucavado, A.; Escalante, T.; Shannon, J.D.; Ayala-Castro, C.N.; Villalta, M.; Gutiérrez, J.M.; Fox, J.W. Efficacy of IgG and F(ab′) 2 antivenoms to neutralize snake venom-induced local tissue damage as assessed by the proteomic analysis of wound exudate. J. Proteome Res. 2012, 11, 292–305. [Google Scholar] [CrossRef] [PubMed]
  16. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Bioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37, 1–13. [Google Scholar] [CrossRef] [PubMed]
  17. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef] [PubMed]
  18. Escalante, T.; Ortiz, N.; Rucavado, A.; Sanchez, E.F.; Richardson, M.; Fox, J.W.; Gutiérrez, J.M. Role of collagens and perlecan in microvascular stability: Exploring the mechanism of capillary vessel damage by snake venom metalloproteinases. PLoS ONE 2011, 6, e28017. [Google Scholar] [CrossRef] [PubMed]
  19. Metsalu, T.; Vilo, J. ClustVis: A web tool for visualizing clustering of multivariate data using Principal Component Analysis and heatmap. Nucleic Acids Res. 2015, 43, W566–W570. [Google Scholar] [CrossRef] [PubMed]
Figure 1. SDS-PAGE profile of samples of blister fluid collected from patients. Samples were run on a 12% gel and stained with Coomassie brilliant blue. Samples from patients bitten by N. naja (lane 1); Hypnale sp. (lanes 2 and 4); D. russelii (lanes 3, 5, and 6). The left lane corresponds to molecular mass markers. The red lines indicate the locations in which the bands were cut for analysis using LC-MS/MS.
Figure 1. SDS-PAGE profile of samples of blister fluid collected from patients. Samples were run on a 12% gel and stained with Coomassie brilliant blue. Samples from patients bitten by N. naja (lane 1); Hypnale sp. (lanes 2 and 4); D. russelii (lanes 3, 5, and 6). The left lane corresponds to molecular mass markers. The red lines indicate the locations in which the bands were cut for analysis using LC-MS/MS.
Toxins 11 00246 g001
Figure 2. Top biological responses, identified via DAVID bionformatic analysis, of the proteomes of blisters’ samples collected from the patients [16,17]. The number of proteins detected corresponding to each process is indicated.
Figure 2. Top biological responses, identified via DAVID bionformatic analysis, of the proteomes of blisters’ samples collected from the patients [16,17]. The number of proteins detected corresponding to each process is indicated.
Toxins 11 00246 g002
Figure 3. Hierarchical clusters created from the proteomic profile of blisters’ samples collected from the patients. Analysis was performed using the Clutvis software [19]. The colors represent a scale on the level of expression that goes from red, representing high abundance proteins, to blue, representing low abundant proteins. Pit viper (Pv) corresponds to Hypnale sp., Rv corresponds to Daboia russelii, and C corresponds to Naja naja.
Figure 3. Hierarchical clusters created from the proteomic profile of blisters’ samples collected from the patients. Analysis was performed using the Clutvis software [19]. The colors represent a scale on the level of expression that goes from red, representing high abundance proteins, to blue, representing low abundant proteins. Pit viper (Pv) corresponds to Hypnale sp., Rv corresponds to Daboia russelii, and C corresponds to Naja naja.
Toxins 11 00246 g003
Figure 4. Detection of venom proteins in human blisters using Western blotting. Blister samples were separated under non-reducing conditions on 12% SDS-PAGE gels and were transferred to a nitrocellulose membrane. Immunodetection was performed with polyspecific VINS antivenom. The reaction was detected using an anti-horse peroxidase antibody and a chemiluminescent substrate. Lanes 1 to 6 correspond to samples of blister fluid (see Table 1). Nn: N. naja venom, Dr: D. russelii venom run in parallel to the blister samples.
Figure 4. Detection of venom proteins in human blisters using Western blotting. Blister samples were separated under non-reducing conditions on 12% SDS-PAGE gels and were transferred to a nitrocellulose membrane. Immunodetection was performed with polyspecific VINS antivenom. The reaction was detected using an anti-horse peroxidase antibody and a chemiluminescent substrate. Lanes 1 to 6 correspond to samples of blister fluid (see Table 1). Nn: N. naja venom, Dr: D. russelii venom run in parallel to the blister samples.
Toxins 11 00246 g004
Table 1. Clinical and laboratory parameters of the six patients whose blister fluid was analyzed in this work.
Table 1. Clinical and laboratory parameters of the six patients whose blister fluid was analyzed in this work.
Snake SpeciesNaja najaHypnale sp.Daboia russeliiHypnale sp.Daboia russeliiDaboia russelii
Patient number123456
SexMaleMaleMaleFemaleMaleFemale
Age53 years old62 years old22 years old2 years old23 years old67 years old
Bite siteLeft thumbLeft big toeLeft legLeft index fingerLeft ankleRight leg
Time between the bite and arrival in the hospital2 h4 h 20 min8 h4 h 15 min8 h 35 min3 h 20 min
SymptomsPtosis and external ophthalmoplegia developed 1 h after arrival to the hospital. Local pain, swelling, and bluish discoloration of left thumb. No nausea, vomiting nor abdominal pain.Local pain, swelling, and bluish discoloration of left big toe. Nausea, vomiting, and abdominal pain present.Disorientation, left leg swelling, and blisters. Under the influence of alcohol at the time of bite. Bilateral parotid swelling and conjunctival chemosis. Nausea, vomiting, and painLocal pain, swelling, and discoloration and blister on the left index finger. No nausea, vomiting, nor abdominal pain.Local pain and swelling. Proteinuria. Nausea, vomiting, and abdominal pain present.Local pain and swelling Bite mark present. Nausea and vomiting present, and no abdominal pain present.
Platelet count177,000/µL70,000/µL20,000/µL300,000/µL160,000/µL120,000/µL
Hemoglobin12.8 g/dL10 g/dL15.7 g/dL11.7 g/dL17.9 g/dL14.1 g/dL
Serum creatinine0.86 mg/dL3.47 mg/dL7.94 mg/dL0.41 mg/dL2.40 mg/dL0.80 mg/dL
Serum K+4.1 mEq/L3 mEq/L4.8 mEq/Ll4 mEq/L6.7 mEq/L
Serum Na+ 126 mEq/L133 mEq/L130 mEq/L131 mEq/L
Whole blood clotting timeNot determined>20 min>20 min12 min>20 min D dimer 4 µg/mL (positive) (normal range <0.5 µg/mL)>20 min
Treatment20 vials antivenom, with neostigmine given 1 h after antivenom; ophthalmoplegia and ptosis disappeared.No antivenom was administered. Patient had acute kidney injury, underwent hemodialysis.20 vials of antivenom. Patient had acute kidney injury, underwent hemodialysis.Patient treated symptomatically. No antivenom was administered.20 vials of antivenom. Patient had acute kidney injury, underwent hemodialysis.20 vials of antivenom. Early adverse reaction to antivenom that was managed with adrenaline.
ObservationsMultiple tourniquets appliedNo tourniquet was appliedNo tourniquet was appliedTourniquet appliedTraditional medicine received prior to reaching the hospital. Tourniquet applied.No tourniquet was applied
Table 2. The most abundant proteins identified in each blister’s samples via proteomic analysis. The top 15 proteins by sample are highlighted in gray. Proteins were identified from the raw data using the Sequest search algorithm Proteome Discovery 1.4.1 against a nonspecific database and the functional annotation provided by DAVID Bioinformatics Resources 6.7.
Table 2. The most abundant proteins identified in each blister’s samples via proteomic analysis. The top 15 proteins by sample are highlighted in gray. Proteins were identified from the raw data using the Sequest search algorithm Proteome Discovery 1.4.1 against a nonspecific database and the functional annotation provided by DAVID Bioinformatics Resources 6.7.
Identified ProteinsAccessionD.russelliHypnale sp.N.naja
356241
Cluster of serum albuminP02768189218371620114225251891
Cluster of alpha-1-antitrypsinP010098513937625136102
Cluster of serotransferrinP02787134926746443366
Cluster of hemoglobin subunit deltaP02042034166414486105
Cluster of keratin, type II cytoskeletal 1P042642726465364155371
Cluster of immunoglobulin kappa constantP018342534528264359306
Cluster of complement C3P01024035527118220148
Cluster of hemoglobin subunit alphaP6990543221891844655
Cluster of alpha-1-antichymotrypsinP0101134952205215858
Cluster of alpha-2-macroglobulinP01023117718224345126
Cluster of immunoglobulin gamma-1 heavy chainP0DOX52926219081301249
Cluster of trypsinP00761140100441068891
Cluster of immunoglobulin lambda constant 2P0DOY21215112429193149
Cluster of carbonic anhydrase 1P0091501082305422249
Cluster of keratin, type I cytoskeletal 10P1364503924243104205
Cluster of transthyretinP027669677444110349
Cluster of haptoglobinQ625584150186178631
Cluster of alpha-1-acid glycoprotein 1P0276320359271536545
Cluster of antithrombin-IIIP010083429252096731
Cluster of isoform 3 of plasma protease C1 inhibitorP05155-33730101194546
Cluster of complement C4-BP0C0L5235853214522
Keratin, type I cytoskeletal 9P3552703533835285
Cluster of hemopexinP02790123973138065
Cluster of ubiquitin-40S ribosomal protein S27aP682020811139816
Cluster of peroxiredoxin-2P321192650817101
Cluster of Immunoglobulin heavy constant muP018710424537625
Cluster of fibronectinP0275115972857852
Protein S100-A9P067024555696514
Cluster of fibrinogen gamma chainP0267929245525553
Cluster of fibrinogen beta chainP02675534900023
Immunoglobulin heavy constant alpha 1P0187612707332474
Cluster of fibrinogen alpha chainP02671114224821218
Cluster of isoform 3 of vitamin D-binding proteinP02774-30324851357
Immunoglobulin heavy constant gamma 2P0185901311013076136
Cluster of superoxide dismutaseP004412711158014
Protein S100-A8P05109323946126
Cluster of ceruloplasminP004501040410242
Cluster of catalaseP04040011732080
Cluster of gelsolinP06396037306428
Immunoglobulin constant gamma 4P01861880503170140
Cluster of plasminogenP007472341877124
Cluster of apolipoprotein A-IG3QY9834127110
Protein AMBPP027602068359312
Zinc-alpha-2-glycoproteinP253115426179259
Cluster of immunoglobulin kappa variable 3-20P016190232233021
Cluster of carbonic anhydrase 2P00918010273010
Alpha-1B-glycoproteinP0421732321253
Cluster of immunoglobulin heavy variable 3-43DP0DP040241214940
ATP-binding cassette sub-family F member 1Q8NE7143562211
Cluster of cartilage oligomeric matrix proteinP4974722308634
Table 3. Proteins in blister fluid in relation to coagulation. Proteins were identified from the raw data using the Sequest search algorithm Proteome Discovery 1.4.1 against the Swiss-Prot database.
Table 3. Proteins in blister fluid in relation to coagulation. Proteins were identified from the raw data using the Sequest search algorithm Proteome Discovery 1.4.1 against the Swiss-Prot database.
Coagulation ProteinsAccessionMWD.russeliiHypnale sp.N.naja
356241
Cluster of fibrinogen alpha chainP0267195 kDa114224821218
Cluster of fibrinogen beta chainP0267556 kDa534900023
Cluster of fibrinogen gamma chainP0267952 kDa29245525553
Cluster of plasminogenP0074791 kDa2341877124
Cluster of prothrombinP0073470 kDa400000
Coagulation factor XIIP0074868 kDa800000
Coagulation factor XIII B chainP0516076 kDa000100
Factor XIIa inhibitorP5044852 kDa000100
Heparin cofactor 2P0554657 kDa042010
Isoform LMW of kininogen-1P01042-248 kDa6600000
Kininogen-1O0867773 kDa100000
von Willebrand factorP04275309 kDa100700
Table 4. Proteins in blister fluid related to the extracellular matrix. Proteins were identified from the raw data using the Sequest search algorithm Proteome Discovery 1.4.1 against the human basic database.
Table 4. Proteins in blister fluid related to the extracellular matrix. Proteins were identified from the raw data using the Sequest search algorithm Proteome Discovery 1.4.1 against the human basic database.
Extracellular Matrix ProteinsAccessionMWD.russeliiHypnale sp.N.naja
356241
Basement membrane-specific heparan sulfate proteoglycan core proteinP98160469 kDa66002205
Cluster of aggrecan core proteinP16112261 kDa100000
Cluster of cartilage oligomeric matrix proteinP4974783 kDa2200000
Cluster of laminin subunit alpha-3Q16787367 kDa400000
Cluster of laminin subunit gamma-2Q13753131 kDa800000
Cluster of nidogen-2Q14112151 kDa700000
Cluster of thrombospondin-1P07996129 kDa021201
Collagen alpha-1(I) chainP02452139 kDa100600
Collagen alpha-1(I) chain (fragments)C0HJN388 kDa000100
Collagen alpha-1(III) chainP0425894 kDa000500
Collagen alpha-1(XIV) chainQ05707194 kDa2100000
Collagen alpha-1(XVI) chainQ07092158 kDa000200
Collagen alpha-1(XXVIII) chainQ2UY09117 kDa100000
Collagen alpha-2(IV) chainP08572168 kDa34001701
Collagen alpha-3(VI) chainP12111344 kDa30023023
Isoform 5 of tenascin-XP22105-4459 kDa5600000
Isoform C of proteoglycan 4Q92954-3141 kDa800000
Laminin subunit alpha-4P97927202 kDa800100
Laminin subunit alpha-5O15230400 kDa3200100
LumicanP5188438 kDa030000
Matrix metalloproteinase-9P1478078 kDa030000
Nidogen-1P14543136 kDa300610
TenascinP24821241 kDa000100
Thrombospondin-4P49744108 kDa000200
Table 5. Proteins from snake venoms identified in blister fluids using proteomic analysis. Proteins were identified from the raw data using the Sequest search algorithm Proteome Discovery 1.4.1 against a nonspecific database.
Table 5. Proteins from snake venoms identified in blister fluids using proteomic analysis. Proteins were identified from the raw data using the Sequest search algorithm Proteome Discovery 1.4.1 against a nonspecific database.
Venom ProteinsAccessionMWD.russeliiHypnale sp.N.naja
356241
Acidic phospholipase A2 inhibitor chain HPD-1A4VBF015 kDa300000
Alpha-cobratoxinP013918 kDa0000016
Beta-fibrinogenaseE0Y41928 kDa100000
Cluster of acidic phospholipase A2P1544513 kDa0100075
Cluster of acidic phospholipase A2 Drk-a1A8CG8615 kDa2700000
Cluster of acidic phospholipase A2 Tbo-E6Q2YHJ316 kDa0002046
Cluster of cytotoxin 2aQ9PST49 kDa0000046
Cluster of factor V activator RVV-V alphaP1896426 kDa1400000
Cluster of Kunitz-type serine protease inhibitor C3A8Y7N69 kDa1000003
Cysteine-rich venom protein kaouthin-1P8480527 kDa0000010
Hemorrhagic metalloproteinase-disintegrin-like kaouthiaginP8294244 kDa000003
Kunitz-type serine protease inhibitor 2Q2ES4910 kDa4800000
Kunitz-type serine protease inhibitor 2P009907 kDa1400000
Kunitz-type serine protease inhibitor C6A8Y7N910 kDa400000
Kunitz-type serine protease inhibitor DrKIn-IH6VC0510 kDa800006
Long neurotoxin 1P256688 kDa0000019
Muscarinic toxin-like protein 2P824637 kDa000002
Muscarinic toxin-like protein 3P824648 kDa000006
Neutral phospholipase A2 RVV-PFIIc’ (fragment)P0DKX12 kDa1000000
Phospholipase A2 1 (fragment)P865292 kDa300000
Serine protease VLSP-3E0Y42028 kDa800000
Snake venom vascular endothelial growth factor toxin VR-1’P0DL4213 kDa1000000
Thrombin-like enzyme 1A7LAC629 kDa300000
Venom nerve growth factorP3089413 kDa2300000
Venom nerve growth factor 2Q5YF8927 kDa0000011
Zinc metalloproteinase-disintegrin-like daborhagin-KB8K1W070 kDa100000

Share and Cite

MDPI and ACS Style

Macêdo, J.K.A.; Joseph, J.K.; Menon, J.; Escalante, T.; Rucavado, A.; Gutiérrez, J.M.; Fox, J.W. Proteomic Analysis of Human Blister Fluids Following Envenomation by Three Snake Species in India: Differential Markers for Venom Mechanisms of Action. Toxins 2019, 11, 246. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins11050246

AMA Style

Macêdo JKA, Joseph JK, Menon J, Escalante T, Rucavado A, Gutiérrez JM, Fox JW. Proteomic Analysis of Human Blister Fluids Following Envenomation by Three Snake Species in India: Differential Markers for Venom Mechanisms of Action. Toxins. 2019; 11(5):246. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins11050246

Chicago/Turabian Style

Macêdo, Jéssica K. A., Joseph K. Joseph, Jaideep Menon, Teresa Escalante, Alexandra Rucavado, José María Gutiérrez, and Jay W. Fox. 2019. "Proteomic Analysis of Human Blister Fluids Following Envenomation by Three Snake Species in India: Differential Markers for Venom Mechanisms of Action" Toxins 11, no. 5: 246. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins11050246

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop