Next Article in Journal
Action of Botulinum Neurotoxin E Type in Experimental Epilepsies
Previous Article in Journal
Biotics (Pre-, Pro-, Post-) and Uremic Toxicity: Implications, Mechanisms, and Possible Therapies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Malaysian and Chinese King Cobra Venom Cytotoxicity in Melanoma and Neonatal Foreskin Fibroblasts Is Mediated by Age and Geography

by
Bianca Op den Brouw
1,†,
Manuel A. Fernandez-Rojo
2,3,†,
Tom Charlton
4,
Bryan G. Fry
1,*,‡ and
Maria P. Ikonomopoulou
5,6,*,‡
1
Venom Evolution Lab, School of Biological Sciences, The University of Queensland, Brisbane, QLD 4072, Australia
2
Hepatic Regenerative Medicine Group, Madrid Institute for Advanced Studies in Food, E28049 Madrid, Spain
3
Diamantina Institute, The University of Queensland, Brisbane, QLD 4072, Australia
4
Department of Natural Sciences, Manchester Metropolitan University, All Saints Building, Manchester M15 6BH, UK
5
Translational Venomics Group, Madrid Institute for Advanced Studies in Food, E28049 Madrid, Spain
6
Institute for Molecular Bioscience, The University of Queensland, Brisbane, QLD 4072, Australia
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
These authors contributed equally to this work.
Toxins 2023, 15(9), 549; https://doi.org/10.3390/toxins15090549
Submission received: 4 July 2023 / Revised: 24 August 2023 / Accepted: 31 August 2023 / Published: 4 September 2023
(This article belongs to the Section Animal Venoms)

Abstract

:
Snake venoms constitute a complex, rapidly evolving trait, whose composition varies between and within populations depending on geographical location, age and preys (diets). These factors have determined the adaptive evolution for predatory success and link venom heterogeneity with prey specificity. Moreover, understanding the evolutionary drivers of animal venoms has streamlined the biodiscovery of venom-derived compounds as drug candidates in biomedicine and biotechnology. The king cobra (Ophiophagus hannah; Cantor, 1836) is distributed in diverse habitats, forming independent populations, which confer differing scale markings, including between hatchlings and adults. Furthermore, king cobra venoms possess unique cytotoxic properties that are used as a defensive trait, but their toxins may also have utility as promising anticancer-agent candidates. However, the impact of geographical distribution and age on these potential venom applications has been typically neglected. In this study, we hypothesised that ontogenetic venom variation accompanies the morphological distinction between hatchlings and adults. We used non-transformed neonatal foreskin (NFF) fibroblasts to examine and compare the variability of venom cytotoxicity between adult captive breeding pairs from Malaysian and Chinese lineages, along with that of their progeny upon hatching. In parallel, we assessed the anticancer potential of these venoms in human-melanoma-patient-derived cells (MM96L). We found that in a geographical distribution and gender-independent manner, venoms from hatchlings were significantly less cytotoxic than those from adults (NFF; ~Log EC50: 0.5–0.6 vs. 0.2–0.35 mg/mL). This is consistent with neonates occupying a semifossorial habitat, while adults inhabit more above-ground habitats and are therefore more conspicuous to potential predators. We also observed that Malaysian venoms exhibited a slightly higher cytotoxicity than those from the Chinese cobra cohorts (NFF; Log EC50: 0.1–0.3 vs. 0.3–0.4 mg/mL), which is consistent with Malaysian king cobras being more strongly aposematically marked. These variations are therefore suggestive of differential anti-predator strategies associated with the occupation of distinct niches. However, all cobra venoms were similarly cytotoxic in both melanoma cells and fibroblasts, limiting their potential medical applications in their native forms.
Key Contribution: The venom of hatchlings is significantly less potent than that from adults King cobras; The Malaysian cobra venom exhibits higher cytotoxicity than that from the Chinese cobra venom.

1. Introduction

Venom is generally defined as a functional trait: a complex secretion produced by a specialised gland tissue that is delivered into the body of another animal via a wound for use in antagonistic interactions [1]. By this biological definition, the majority of snake species are venomous. However, while the concept of “a venomous snake” usually invokes medical notions of significant danger to humans, the toxic secretions of most snakes are mild, and only a small proportion poses a serious threat to humans. Furthermore, the innumerable bioactive molecules (i.e., toxins) comprising these secretions possess pharmacological properties, such as analgesic, metabolic, immunomodulating, anti-thrombotic or anticancer ones, that can be exploited as therapeutics and developed into life-saving medicines [2,3].
The scale of mining venoms for drug candidates is immense, as this mixture often contains hundreds of toxins and isoforms that undergo accelerated rates of evolution [4]. This biochemical complexity underpins the striking diversity in toxin expression found across species, populations, and even within individuals throughout their lifetime [5]. This variability deepens the reservoir of pharmacological potential while highlighting new medical applications [6].
Studying the natural history and evolution of venomous species provides key insights regarding the biological sources reflecting toxin diversity. Snake venoms function at the ecological and biochemical interface between predator and prey, aiming mainly to facilitate predation [7]. This has led investigations of drivers in snake venom variation to focus on prey preferences [6,8,9,10]. However, this approach overlooks the ecological complexity that underlies the evolution of a phenotype. Numerous biotic and abiotic factors, such as behaviour, morphology, prey diversity, population size, biogeography, and climate, variably interact to produce a “fingerprint” of macro- and microevolutionary processes within a population. Consequently, it shapes the biochemical architecture of both predator and prey. Accordingly, inter-populational differences in any of these components can be indicators of cryptic diversity [11,12,13,14,15,16].
The king cobra (Ophiophagus hannah; Cantor, 1836) belongs to the Elapidae family and is widely distributed from India to southern China and Southeast Asia in numerous geographically isolated mainland and island populations. These populations form at least four genetically distinct lineages, which are likely separate species [17]. They occupy a diversity of ecotypes across their distribution, including rainforest, alpine forest, plantations, paddy fields, and mangrove swamps, across tropical, subtropical, and temperate climes [18,19,20]. Morphological variation in coloration and banding is evident between lineages, as well as between juveniles and adults (Figure 1) and between sexes for some localities [21].
The king cobra is not a “true” cobra (Naja sp.), despite being named as such. This misnomer likely derives from the cobra-like hooding behaviour that it adopts when feeling threatened, its immense length of 5–6 metres [18], and its tendency to feed primarily on other snakes—including highly venomous species [18]. Adults yield tremendous quantities of venom dominated by neurotoxic three-finger-toxins (3FTx), along with lesser quantities of toxins such as phospholipases (PLA2), snake venom metalloproteases (SVMP), and L-amino acid oxidases (LAAO), the latter of which are evolutionarily distant from those of other snake LAAOs and have shown promise as anticancer agents [22,23,24,25,26]. Venom variation between individuals within Thailand has been documented, as well as between pooled venoms of adults from Thai, Malaysian, Chinese, and Indonesian populations, with significant variation in antivenom effectiveness [27,28,29]. Notably, the Malaysian population has been shown to be the most cytotoxic [12] and least neurotoxic [30]. Distinctions between the venoms of juveniles and adults or between genders are still lacking for any geographical population.
The king cobra complex represents an ecologically and toxicologically interesting group that harbours a significant degree of venom diversity and pharmacological potential. We hypothesised that ontogenetic venom variation accompanies morphological alterations evident between adult and hatchling snakes as a reflection of differences in niche occupation. In this study, we examined the cytotoxicity of juvenile and adult king cobras from Chinese and Malaysian lineages using neonatal foreskin (NFF) cells and discussed these results in the context of king cobra ecology. We additionally assessed their anticancer potential in human melanoma (MM96L) cells.

2. Results

To examine cytotoxicity in non-transformed cells, we treated neonatal foreskin fibroblasts (NFFs) with venoms for 24 h at 2.5 µg/mL. It was shown that gender was not a significant factor determining the diversity of cytotoxicity between Chinese and Malaysian king cobras (Figure 2). Indeed, venoms from adult male and female snakes demonstrated high levels of cytotoxicity, with the Chinese venoms reducing proportions of viable cells down to 15 ± 3.0% (male) and 19 ± 6.4% (female) (unpaired t-test: t(4) = 0.8599, p = 0.438353) and those of Malaysia doing so down to 9 ± 3.5% (male) and 8 ± 3.1% (unpaired t-test: t(4) = 0.3689, p = 0.730875) (Figure 2).
However, venom-induced cytotoxicity was mediated by age and geographical distribution in cobras. Unlike adults—females and males—venom from hatchlings exhibited less impact on NFF cell viability (Chinese King Cobra, 91 ± 8.2%; Malaysian King cobra, 65 ± 29.3% of viability (two-way ANOVA: F3,16 = 40, p < 0.0001)). This difference was already apparent when we assessed the venom of a juvenile (2 year old) Malaysian snake. The latter reduced cell viability to 14 ± 1.0%, a similar cytotoxicity to that observed in adult snake venoms (Tukey’s, p = 0.0002). Moreover, geographical origin was found to be a significant co-factor in cytotoxicity but a weaker source of venom diversity. Indeed, Malaysian venoms exhibited a more profound cytotoxicity than their Chinese counterparts (two-way ANOVA: F1,16 = 6.892, p = 0.0184).
Cell imaging of NFF cells with light microscopy during the first four hours of incubation with cobra venoms showed that at high concentrations (50 μg/mL), their cytotoxicity was obvious. Indeed, there was evidence of cell stress, such as changes in cellular morphology that exhibited a rounder shape and lower confluency, which translated to decreased cell-to-cell contact and lesser cell adhesion after one hour of incubation. These symptoms were more apparent after two hours in adult venoms from both species (Figure 3). Interestingly, although cell imaging showed that venoms from hatchlings exhibited cytotoxicity at high concentrations, cells were visually less affected at 4 h of treatment in comparison to adult snake venoms. A limited effect in the viability of non-transformed NFF cells is an essential asset of anticancer drug candidates. In our study, this was evident in hatchlings’ venoms. Therefore, venom samples were assessed further by increasing the incubation period to 48 h in both healthy (NFF) and melanoma (MM96L) cells at a range of venom concentrations.
At concentrations of ≥5 μg/mL, all venoms, including those from hatchlings, diminished the viability of melanoma cells to fewer than 25%, showcasing their potency (Figure 4B). Of interest, at this incubation period, the age-associated variations previously observed at 24 h of incubation were greatly reduced or negated in non-transformed fibroblasts. Indeed, all venoms equally decreased the proportion of viable NFF cells to ≤30% at 48 h (Figure 4A), but their cytotoxicity tended to be higher in melanoma cells, indicating a preference for targeting tumorigenic cells.
At high venom concentrations and at an increased incubation period, the geographic-associated variations in cytotoxicity in NFF cells observed at 2.5 μg/mL were also greatly lessened. However, the EC50s calculated from the dose-response data indicated that the Malaysian venoms were, broadly speaking, slightly more potent than their Chinese counterparts in both cell lines (Table 1). Although venoms generally diminished the viability in melanoma cells to a slightly greater extent than in fibroblasts, their EC50s implied a similarly high potency.

3. Discussion

Venoms from the Malaysian king cobras generally demonstrated a greater cytotoxicity than their Chinese counterpart snakes, which is consistent both with the Malaysian population being more strongly marked with a defensive colouration and with previous studies showing a higher level of cytotoxicity and consequent lower levels of neurotoxins present in the venoms [12,30]. Of interest, the venom of hatchlings from both localities possessed significantly weaker cytotoxic properties than that of adults. This is in accordance with their relative conspicuousness to predators, with juveniles being semi-fossorial, while adults are more active in the open and therefore more likely to encounter potential predators. In addition, most venoms demonstrated a slightly greater, but not significant, potency in melanoma when compared with non-transformed fibroblasts. As crude venoms were used, the observed cell cytotoxicity could be the result of a potent toxin or the synergistic effect of two or multiple toxins. However, the observed broad-spectrum potency in healthy and cancerous cells suggests that any anticancer therapeutic potential of these venoms should be formulated with direct delivery approaches into tumors. Indeed, a limitation of investigating the venom lies in the quantity needed to extensively characterise it. In addition, the complexity of the crude venom’s components prompts one to identify the toxin reflecting the biological activity. However, an even more promising strategy is the chemical synthesis of the toxin of interest. Chemically synthesised compounds could be further optimised for their efficacy and safety and could therefore set the foundations for the next generation of anticancer drug candidates. Promising examples of such toxins targeting various tumor cells are reported in the literature, including peptides inspired from Lys49 phospholipase A2 snake venoms [31,32,33].
While king cobra venoms do possess cytotoxic 3FTxs (CTX), LAAOs are primarily responsible for their cytotoxicity [34,35]. King cobra LAAOs are exceptionally potent and are known to possess a number of unique drug assets, including heat stability and retention of activity over extended periods of time [24,36]. Accordingly, the potency of LAAOs in hatchlings’ venoms at 48 h could account for the observed cell death.
Age is a common factor of venom variation among snakes and is often linked to parallel shifts in prey and, potentially, diet [37]. Unfortunately, little knowledge exists on the ecology of hatchling king cobras, including their diet composition. The literature suggests that cytotoxins in hooding Asian and African elapids are defensive, rather than predatory, venom components [11,12]. Previously, we reported that cytotoxicity evolved primarily as a defensive innovation and that it has co-evolved twice alongside hooding behavior and independently of spitting [12]. Newly hatched king cobras are markedly more distinct in patterning and coloration than adults and are thought to be largely arboreal during the first year of their life, transitioning to a semi-arboreal lifestyle between 1–2 years of age, and becoming mostly terrestrial as adults, though they may also hunt arboreally [38]. The temporal increase in venom cytotoxicity observed in this study generally coincides with these shifts in niche occupation. With the exception of instances of cannibalism, large adult king cobras have few natural predators, though subadults may be vulnerable to predation by terrestrial carnivores [39], with species such as mongooses (Herpestidae) and hog badgers (Mustelidae) showcasing exceptional skills in overpowering venomous snakes and resistance to snake venom neurotoxins [40]. Consequently, pain-inducing cytotoxins form a crucial line of defence against such mammalian predators, in this instance during the transitional period from the occupation of arboreal to terrestrial environments. Juveniles may be easily overcome by predators of sufficient size, including carnivorous birds, monitor lizards, and even other snakes [38]. Cytotoxic venom may therefore be of limited defensive use to hatchlings, who likely rely instead on predator avoidance through their arboreality and the crypsis conferred by their distinct banding in these light-mottled arboreal habitats.
The extent of dietary differences among king cobra populations is understudied, though divergence in prey preference (e.g., from primarily snakes to primarily monitor lizards) has been observed in some localities [41], along with significant differences in taxon-specific venom lethality among populations [29]. Further studies into ontogenetic variations in the ecology and venom difference of these snakes are needed to test this hypothesis, though assessing ontogenetic and geographic differences in neurotoxicity alongside cytotoxicity may offer some broad insights into the role of predatory versus defensive venom components.
While gender was not found to be a factor of venom variation in this study, no firm conclusions can be drawn on this variable, as only one adult of each gender was tested for each location, and venom is known to vary significantly among individuals [29]. Such individual variation is perhaps unsurprising, given the discrepancies between habitat preference and home range found between individuals of a given locality [20,41]. However, sexual dimorphism and sex-linked differences in behaviour are evident within king cobras, and therefore broad gender-linked trends in venom composition may also occur. Increasing the sample size is necessary to adequately test this hypothesis.
The king cobra complex contains venom variation that is worthy of concerted research attention. Although the complex is yet to be taxonomically revised, the venom variation observed in this study reinforces the need to clearly define origin and age in snake samples in future work. Research on the ecology and venom composition of each lineage of these snakes is essential for a further assessment of their pharmacological potential and for the conservation of these biologically important snakes.

4. Conclusions

We observed that age was a contributor of cytotoxicity in Malaysian and Chinese king cobra. We also reported that the venom of hatchlings was significantly less potent than that of adults, based on the cytotoxicity results in fibroblasts. This is in alignment with their ecology, where neonates are born as arboreal before transitioning to a semi-arboreal lifestyle and finally transforming into terrestrial adults. The temporal increase in venom cytotoxicity coincides with these shifts in niche occupation, which reflect their defensive behaviour as adults. Therefore, this work sheds light on the association between ecology and cytotoxicity in king cobras. It also generates new knowledge on the impact of geographical distribution, probably due to diet preferences, on venom cytotoxicity. From a therapeutic discovery perspective, the observation that venom from hatchlings is less cytotoxic than that from adults in fibroblasts, but still effective at targeting melanoma cells, might suggest that the former could be more suitable for further evaluation and development as a drug candidate.

5. Materials and Methods

5.1. Venoms

Venoms were obtained from captive specimens of O. hannah from a private collection in the UK. All venom study protocols in this work were performed with the University of Queensland Biosafety Approval #IBC134BSBS2015 and University of Queensland Animal Ethics Approval 2021/AE000075. Pairs of adult snakes from Malaysian and Chinese lineages were bred in captivity. The venom from the adults and that of their young was extracted within approximately eight weeks after hatchlings emerged. Additionally, venom was extracted from a 2 year old captive-bred Malaysian female king cobra hatched from a previous breeding by the same adult pair of Malaysian adults. Venoms were immediately flash-frozen in liquid nitrogen, lyophilised, and stored at −80 °C. Prior to testing, the venoms of sibling hatchlings were reconstituted in deionised water and pooled.

5.2. Cytotoxicity Analyses

We used human melanoma cells (MM96L) and human neonatal foreskin fibroblasts (NFFs), supplied by QIMR Berghofer Medical Research Institute. Cells were maintained in RPMI medium supplemented with 1% penicillin–streptomycin and foetal calf serum (FCS), as previously described [6,42,43]. Cell viability was evaluated using colourimetric MTT (Thiazolyl Blue Tetrazolium Bromide; Sigma-Aldrich, St Louis, MO, USA) assays [6,42,43]. The following venom concentrations of 1 µg/mL, 2.5 µg/mL, 5 µg/mL, 10 µg/mL, or 50 µg/mL were added to cells (n = 4) and incubated for 24 h or 48 h. We used 10% of dimethyl sulfoxide (DMSO) as a positive control (100% toxicity). The absorbance was read at 570 nm. We performed three independent experiments with a minimum of three replicates per treatment. The readings from all treatments were normalised against their untreated control cells, which were considered as 100% viability and then subtracted from wells only containing media.

5.3. Light Microscopy

Cells were plated in 24-well plates (8000/well) and cultured for 0 to 4 h in the presence or absence of Chinese or Malaysian (babies and adults of both genders) cobra venoms at a concentration of 50 µg/mL. We conducted three independent experiments of three replicates each. Then, images of the cells were taken at the mentioned times using a Leica DMIL LED epifluorescence microscope and using the LAS V413 software version number LASV4.13.0. Panel figures were set using ADOBE photoshop CS3.

5.4. Statistics

Data were graphed and analysed using GraphPad Prism version 9 for iOS (GraphPad Software, San Diego, CA, USA). For the MTT assays, venom cytotoxicity data were normalised to the vehicle-treated cells (100% cell viability) to calculate the percentage of viable cells following venom treatment. Data were then tested for normality using Shapiro–Wilk tests (significance set to p < 0.05). To assess differences between genders, the percentages of viable cells were compared using unpaired t-tests (males vs. females), since data were normally distributed. As no difference was evident, cell viability data from adult snakes were combined for subsequent statistical analyses. Two-way ANOVAs with Tukey’s post hoc comparisons were then run to evaluate geography and age-linked venom variation. To calculate venom EC50s, dose-response data were transformed (venom concentrations log-transformed) and normalised (smallest mean 0%, largest mean 100%), and then curves were subjected to non-linear regression using the [log(agonist) vs. normalised response − variable slope] prism function.

Author Contributions

Conceptualization, B.O.d.B., B.G.F. and M.P.I.; Methodology, M.A.F.-R. and M.P.I.; Resources, M.A.F.-R., T.C. and M.P.I. Writing—Original Draft Preparation, B.O.d.B.; Writing—Review & Editing, B.O.d.B., M.A.F.-R., T.C., B.G.F. and M.P.I.; Funding Acquisition, M.P.I. All authors have read and agreed to the published version of the manuscript.

Funding

M.A. Fernandez-Rojo was supported by the Talento Program from the Madrid Government of Spain (Grants No. T1-BIO-1854 and 2020-5A_BIO-19724). M.P. Ikonomopoulou was supported by the AMAROUT Marie Curie program (Nº 291803-AMAROUT II), the TALENTO Program (2018-T1/BIO-11262; 2022-5A/BIO-24228) and Youth Funding (PEJ-2020-AI/BIO-17904) by the Regional Madrid Government of Spain, as well as the National Plan (PID2021-126691OB-I00) by the Spanish Government.

Institutional Review Board Statement

This work was performed with the University of Queensland Biosafety Approval #IBC134BSBS2015 and University of Queensland Animal Ethics Approval 2021/AE000075.

Informed Consent Statement

Not applicable

Data Availability Statement

The data presented in this study are available upon request to the corresponding authors.

Acknowledgments

The cell lines were provided by Glen Boyle from QIMR Berghofer.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jackson, T.N.; Fry, B.G. A tricky trait: Applying the fruits of the “function debate” in the philosophy of biology to the “venom debate” in the science of toxinology. Toxins 2016, 8, 263. [Google Scholar] [CrossRef] [PubMed]
  2. von Reumont, B.M.; Anderluh, G.; Antunes, A.; Ayvazyan, N.; Beis, D.; Caliskan, F.; Crnković, A.; Damm, M.; Dutertre, S.; Ellgaard, L.; et al. Modern venomics-Current insights, novel methods, and future perspectives in biological and applied animal venom research. Gigascience 2022, 11, giac048. [Google Scholar]
  3. Smith, J.J.; Lau, C.H.Y.; Herzig, V.; Ikonomopoulou, M.P.; King, G.F.; Rash, L.D. Therapeutic applications of spider-venom peptides. In Venoms to Drugs: Venom as a Source for the Development of Human Therapeutics; King, G.F., Ed.; Royal Society of Chemistry: London, UK, 2014; Chapter 8. [Google Scholar]
  4. Shibata, H.; Chijiwa, T.; Oda-Ueda, N.; Nakamura, H.; Yamaguchi, K.; Hattori, S.; Matsubara, K.; Matsuda, Y.; Yamashita, A.; Isomoto, A. The habu genome reveals accelerated evolution of venom protein genes. Sci. Rep. 2018, 8, 11300. [Google Scholar] [CrossRef] [PubMed]
  5. Casewell, N.R.; Wüster, W.; Vonk, F.J.; Harrison, R.A.; Fry, B.G. Complex cocktails: The evolutionary novelty of venoms. Trends Ecol. Evol. 2013, 28, 219–229. [Google Scholar] [CrossRef] [PubMed]
  6. Op den Brouw, B.; Coimbra, F.C.P.; Bourke, L.A.; Huynh, T.M.; Vlecken, D.H.W.; Ghezellou, P.; Visser, J.C.; Dobson, J.S.; Fernandez-Rojo, M.A.; Ikonomopoulou, M.P.; et al. Extensive Variation in the Activities of Pseudocerastes and Eristicophis Viper Venoms Suggests Divergent Envenoming Strategies Are Used for Prey Capture. Toxins 2021, 13, 112. [Google Scholar] [CrossRef]
  7. Fry, B.G.; Sunagar, K.; Casewell, N.R.; Kochva, E.; Roelants, K.; Scheib, H.; Wüster, W.; Vidal, N.; Young, B.; Burbrink, F.; et al. The origin and evolution of the Toxicofera reptile venom system. In Venomous Reptiles and Their Toxins: Evolution, Pathophysiology and Biodiscovery; Fry, B.G., Ed.; Oxford University Press: New York, NY, USA, 2015; Chapter 1; pp. 1–31. [Google Scholar]
  8. Barlow, A.; Pook, C.E.; Harrison, R.A.; Wüster, W. Coevolution of diet and prey-specific venom activity supports the role of selection in snake venom evolution. Proc. R. Soc. B Biol. Sci. 2009, 276, 2443–2449. [Google Scholar] [CrossRef]
  9. Bernardoni, J.L.; Sousa, L.F.; Wermelinger, L.S.; Lopes, A.S.; Prezoto, B.C.; Serrano, S.M.T.; Zingali, R.B.; Moura-da-Silva, A.M. Functional Variability of Snake Venom Metalloproteinases: Adaptive Advantages in Targeting Different Prey and Implications for Human Envenomation. PLoS ONE 2014, 9, e109651. [Google Scholar] [CrossRef]
  10. Healy, K.; Carbone, C.; Jackson, A.L. Snake venom potency and yield are associated with prey-evolution, predator metabolism and habitat structure. Ecol. Lett. 2019, 22, 527–537. [Google Scholar] [CrossRef]
  11. Kazandjian, T.D.; Petras, D.; Robinson, S.D.; Thiel, J.V.; Greene, H.W.; Arbuckle, K.; Barlow, A.; Carter, D.A.; Wouters, R.M.; Whiteley, G.; et al. Convergent evolution of pain-inducing defensive venom components in spitting cobras. Science 2021, 371, 386–390. [Google Scholar] [CrossRef]
  12. Panagides, N.; Jackson, T.N.W.; Ikonomopoulou, M.P.; Arbuckle, K.; Pretzler, R.; Yang, D.C.; Ali, S.A.; Koludarov, I.; Dobson, J.; Sanker, B.; et al. How the Cobra Got Its Flesh-Eating Venom: Cytotoxicity as a Defensive Innovation and Its Co-Evolution with Hooding, Aposematic Marking, and Spitting. Toxins 2017, 9, 103. [Google Scholar] [CrossRef]
  13. Strickland, J.L.; Smith, C.F.; Mason, A.J.; Schield, D.R.; Borja, M.; Castañeda-Gaytán, G.; Spencer, C.L.; Smith, L.L.; Trápaga, A.; Bouzid, N.M. Evidence for divergent patterns of local selection driving venom variation in Mojave Rattlesnakes (Crotalus scutulatus). Sci. Rep. 2018, 8, 17622. [Google Scholar] [CrossRef]
  14. Zancolli, G.; Calvete, J.J.; Cardwell, M.D.; Greene, H.W.; Hayes, W.K.; Hegarty, M.J.; Herrmann, H.-W.; Holycross, A.T.; Lannutti, D.I.; Mulley, J.F. When one phenotype is not enough: Divergent evolutionary trajectories govern venom variation in a widespread rattlesnake species. Proc. R. Soc. B 2019, 286, 2018–2735. [Google Scholar] [CrossRef] [PubMed]
  15. Aird, S.D.; Arora, J.; Barua, A.; Qiu, L.; Terada, K.; Mikheyev, A.S. Population genomic analysis of a pitviper reveals microevolutionary forces underlying venom chemistry. Genome Biol. Evol. 2017, 9, 2640–2649. [Google Scholar] [CrossRef] [PubMed]
  16. Schendel, V.; Rash, L.D.; Jenner, R.A.; Undheim, E.A. The diversity of venom: The importance of behavior and venom system morphology in understanding its ecology and evolution. Toxins 2019, 11, 666. [Google Scholar] [CrossRef] [PubMed]
  17. Gowri Shankar, P.; Swamy, P.; Williams, R.C.; Ganesh, S.R.; Moss, M.; Höglund, J.; Das, I.; Sahoo, G.; Vijayakumar, S.P.; Shanker, K.; et al. King or royal family? Testing for species boundaries in the King Cobra, Ophiophagus hannah (Cantor, 1836), using morphology and multilocus DNA analyses. Mol. Phylogenet. Evol. 2021, 165, 107300. [Google Scholar] [CrossRef]
  18. Lim, K.K.; Leong, T.M.; Lim, F.L. The king cobra, Ophiophagus hannah (cantor) in singapore (Reptilia: Squamata: Elapidae). Nat. Singap. 2011, 4, 143–156. [Google Scholar]
  19. Bashir, T.; Poudyal, K.; Bhattacharya, T.; Sathyakumar, S.; Subba, J. Sighting of King Cobra Ophiophagus hannah in Sikkim, India: A new altitude record for the northeast. J. Threat. Taxa 2010, 2, 990–991. [Google Scholar] [CrossRef]
  20. Rao, C.; Talukdar, G.; Choudhury, B.; Shankar, P.G.; Whitaker, R.; Goode, M. Habitat use of king cobra (Ophiophagus hannah) in a heterogeneous landscape matrix in the tropical forests of the Western Ghats, India. Hamadryad 2013, 36, 69–79. [Google Scholar]
  21. Shankar, P.G.; Ganesh, S.; Whitaker, R.; Prashanth, P. King Cobra Ophiophagus hannah (Cantor, 1836) encounters in human-modified rainforests of the Western Ghats, India. Hamadryad 2013, 36, 62–68. [Google Scholar]
  22. Jin, Y.; Lee, W.-H.; Zeng, L.; Zhang, Y. Molecular characterization of l-amino acid oxidase from king cobra venom. Toxicon 2007, 50, 479–489. [Google Scholar] [CrossRef]
  23. Petras, D.; Heiss, P.; Süssmuth, R.D.; Calvete, J.J. Venom Proteomics of Indonesian King Cobra, Ophiophagus hannah: Integrating Top-Down and Bottom-Up Approaches. J. Proteome Res. 2015, 14, 2539–2556. [Google Scholar] [CrossRef] [PubMed]
  24. Lee, M.L.; Fung, S.Y.; Chung, I.; Pailoor, J.; Cheah, S.H.; Tan, N.H. King cobra (Ophiophagus hannah) venom L-amino acid oxidase induces apoptosis in PC-3 cells and suppresses PC-3 solid tumor growth in a tumor xenograft mouse model. Int. J. Med. Sci. 2014, 11, 593–601. [Google Scholar] [CrossRef]
  25. Tan, C.H.; Tan, K.Y.; Fung, S.Y.; Tan, N.H. Venom-gland transcriptome and venom proteome of the Malaysian king cobra (Ophiophagus hannah). BMC Genom. 2015, 16, 687. [Google Scholar] [CrossRef]
  26. Melani, R.D.; Skinner, O.S.; Fornelli, L.; Domont, G.B.; Compton, P.D.; Kelleher, N.L. Mapping Proteoforms and Protein Complexes From King Cobra Venom Using Both Denaturing and Native Top-down Proteomics. Mol. Cell. Proteom. 2016, 15, 2423–2434. [Google Scholar] [CrossRef] [PubMed]
  27. Wongtay, P.; Sangtanoo, P.; Sangvanich, P.; Karnchanatat, A. Variation in the Protein Composition and Biological Activity of King Cobra (Ophiophagus hannah) Venoms. Protein J. 2019, 38, 565–575. [Google Scholar] [CrossRef]
  28. Tan, K.Y.; Ng, T.S.; Bourges, A.; Ismail, A.K.; Maharani, T.; Khomvilai, S.; Sitprija, V.; Tan, N.H.; Tan, C.H. Geographical variations in king cobra (Ophiophagus hannah) venom from Thailand, Malaysia, Indonesia and China: On venom lethality, antivenom immunoreactivity and in vivo neutralization. Acta Trop. 2020, 203, 105311. [Google Scholar] [CrossRef]
  29. Chang, H.-C.; Tsai, T.-S.; Tsai, I.-H. Functional proteomic approach to discover geographic variations of king cobra venoms from Southeast Asia and China. J. Proteom. 2013, 89, 141–153. [Google Scholar] [CrossRef]
  30. Chandrasekara, U.; Harris, R.J.; Fry, B.G. The Target Selects the Toxin: Specific Amino Acids in Snake-Prey Nicotinic Acetylcholine Receptors That Are Selectively Bound by King Cobra Venoms. Toxins 2022, 14, 528. [Google Scholar] [CrossRef]
  31. Araya, C.; Lomonte, B. Antitumor effects of cationic synthetic peptides derived from Lys49 phospholipase A2 homologues of snake venoms. Cell Biol. Int. 2007, 31, 263–268. [Google Scholar] [CrossRef]
  32. Peña-Carrillo, M.S.; Pinos-Tamayo, E.A.; Mendes, B.; Domínguez-Borbor, C.; Proaño-Bolaños, C.; Miguel, D.C.; Almeida, J.R. Dissection of phospholipases A2 reveals multifaceted peptides targeting cancer cells, Leishmania and bacteria. Bioorganic Chem. 2021, 114, 105041. [Google Scholar] [CrossRef]
  33. Almeida, J.R.; Mendes, B.; Lancellotti, M.; Franchi, G.C., Jr.; Passos, Ó.; Ramos, M.J.; Fernandes, P.A.; Alves, C.; Vale, N.; Gomes, P.; et al. Lessons from a Single Amino Acid Substitution: Anticancer and Antibacterial Properties of Two Phospholipase A2-Derived Peptides. Curr. Issues Mol. Biol. 2022, 44, 46–62. [Google Scholar] [CrossRef]
  34. Ahn, M.Y.; Lee, B.M.; Kim, Y.S. Characterization and cytotoxicity of l-amino acid oxidase from the venom of king cobra (Ophiophagus hannah). Int. J. Biochem. Cell Biol. 1997, 29, 911–919. [Google Scholar] [CrossRef]
  35. Ahn, M.Y.; Lee, B.M.; Kim, Y.S. Cytotoxicity and L-amino acid oxidase activity of animal venoms. Arch. Pharm. Res. 1997, 20, 13–16. [Google Scholar] [CrossRef] [PubMed]
  36. Tan, N.H.; Saifuddin, M.N. Isolation and characterization of an unusual form of L-amino acid oxidase from King cobra (Ophiophagus hannah) venom. Biochem. Int. 1989, 19, 937–944. [Google Scholar] [PubMed]
  37. Andrade, D.V.; Abe, A.S. Relationship of venom ontogeny and diet in Bothrops. Herpetologica 1999, 55, 200–204. [Google Scholar]
  38. Charlton, T. King Cobra Natural History and Captive Management; Natural History Publications Borneo: Kota Kinabalu, Malaysia, 2018. [Google Scholar]
  39. Marshall, B.; Strine, C.; Jones, M.; Theodorou, A.; Amber, E.; Waengsothorn, S.; Suwanwaree, P.; Goode, M. Hits Close to Home: Repeated Persecution of King Cobras (Ophiophagus hannah) in Northeastern Thailand. Trop. Conserv. Sci. 2018, 11, 1940082918818401. [Google Scholar] [CrossRef]
  40. Khan, M.A.; Dashevsky, D.; Kerkkamp, H.; Kordiš, D.; de Bakker, M.A.G.; Wouters, R.; van Thiel, J.; op den Brouw, B.; Vonk, F.J.; Kini, R.M.; et al. Widespread Evolution of Molecular Resistance to Snake Venom α-Neurotoxins in Vertebrates. Toxins 2020, 12, 638. [Google Scholar] [CrossRef]
  41. Marshall, B.M.; Strine, C.T.; Jones, M.D.; Artchawakom, T.; Silva, I.; Suwanwaree, P.; Goode, M. Space fit for a king: Spatial ecology of king cobras (Ophiophagus hannah) in Sakaerat Biosphere Reserve, Northeastern Thailand. Amphib. Reptil. 2019, 40, 163–178. [Google Scholar] [CrossRef]
  42. op den Brouw, B.; Ghezellou, P.; Casewell, N.R.; Ali, S.A.; Fathinia, B.; Fry, B.G.; Bos, M.H.A.; Ikonomopoulou, M.P. Pharmacological Characterisation of Pseudocerastes and Eristicophis Viper Venoms Reveal Anticancer (Melanoma) Properties and a Potentially Novel Mode of Fibrinogenolysis. Int. J. Mol. Sci. 2021, 22, 6896. [Google Scholar] [CrossRef]
  43. Goldenberg, J.; Cipriani, V.; Jackson, T.N.; Arbuckle, K.; Debono, J.; Dashevsky, D.; Panagides, N.; Ikonomopoulou, M.P.; Koludarov, I.; Li, B.; et al. Proteomic and functional variation within black snake venoms (Elapidae: Pseudechis). Comp. Biochem. Physiol. C Toxicol. Pharmacol. 2018, 205, 53–61. [Google Scholar] [CrossRef]
Figure 1. Representative image of adult (left) and juvenile (right) Malaysian king cobra provided by Tom Charlton.
Figure 1. Representative image of adult (left) and juvenile (right) Malaysian king cobra provided by Tom Charlton.
Toxins 15 00549 g001
Figure 2. Remaining proportion (%) of viable neonatal foreskin (NFF) cells after 24 h of incubation with 2.5 µg/mL of king cobra (Ophiophagus hannah sp.) venoms from Malaysian (dark blue) and Chinese (yellow) lineages as well as untreated cells (grey)quantified using MTT assays. Data are presented as mean ± SD and are the result of at least three replicates, * p  <  0.05, ** p < 0.0002 & *** p  <  0.0001. Data are normalised to the vehicle (negative) control (100% viability); 10% DMSO was used as a positive control.
Figure 2. Remaining proportion (%) of viable neonatal foreskin (NFF) cells after 24 h of incubation with 2.5 µg/mL of king cobra (Ophiophagus hannah sp.) venoms from Malaysian (dark blue) and Chinese (yellow) lineages as well as untreated cells (grey)quantified using MTT assays. Data are presented as mean ± SD and are the result of at least three replicates, * p  <  0.05, ** p < 0.0002 & *** p  <  0.0001. Data are normalised to the vehicle (negative) control (100% viability); 10% DMSO was used as a positive control.
Toxins 15 00549 g002
Figure 3. Light-microscopy representative images of neonatal foreskin fibroblast (NFF) cells at 63× magnification, incubated at 0–4 h with 50 µg/mL of king cobra (Ophiophagus hannah sp.) venoms from Chinese and Malaysian lineages. Scale bar = 20 µm.
Figure 3. Light-microscopy representative images of neonatal foreskin fibroblast (NFF) cells at 63× magnification, incubated at 0–4 h with 50 µg/mL of king cobra (Ophiophagus hannah sp.) venoms from Chinese and Malaysian lineages. Scale bar = 20 µm.
Toxins 15 00549 g003
Figure 4. The proportion (%) of viable neonatal foreskin fibroblasts (NFFs) (Column (A)) and melanoma (MM96L) cells (Column (B)) after 48 h of incubation with 1–50 μg/mL of king cobra (Ophiophagus hannah sp.) venoms from Malaysian (dark blue) and Chinese (yellow) lineages, quantified using MTT assays. Data are presented as mean ± SD and are the result of at least three independent experiments, performed in triplicates and normalised to vehicle controls (100% viability). All treatments were compared to untreated (control, grey) cells. **** p  <  0.0001.
Figure 4. The proportion (%) of viable neonatal foreskin fibroblasts (NFFs) (Column (A)) and melanoma (MM96L) cells (Column (B)) after 48 h of incubation with 1–50 μg/mL of king cobra (Ophiophagus hannah sp.) venoms from Malaysian (dark blue) and Chinese (yellow) lineages, quantified using MTT assays. Data are presented as mean ± SD and are the result of at least three independent experiments, performed in triplicates and normalised to vehicle controls (100% viability). All treatments were compared to untreated (control, grey) cells. **** p  <  0.0001.
Toxins 15 00549 g004
Table 1. Venom cytotoxicity (EC50s) calculated from the remaining proportion (%) of viable neonatal foreskin fibroblasts (NFFs) and melanoma (MM96L) cells at 48 h of incubation with 1–50 μg/mL of king cobra (Ophiophagus hannah sp.) venoms derived from Malaysian and Chinese lineages and quantified using MTT assays. Data are the result of best-fit values, following non-linear regression of at least three independent experiments conducted in three replicates.
Table 1. Venom cytotoxicity (EC50s) calculated from the remaining proportion (%) of viable neonatal foreskin fibroblasts (NFFs) and melanoma (MM96L) cells at 48 h of incubation with 1–50 μg/mL of king cobra (Ophiophagus hannah sp.) venoms derived from Malaysian and Chinese lineages and quantified using MTT assays. Data are the result of best-fit values, following non-linear regression of at least three independent experiments conducted in three replicates.
Log EC50 (μg/mL)
ChineseMalaysian
MM96LNFFMM96LNFF
Adult (M)0.3 ± 0.060.4 ± 0.050.2 ± 0.070.3 ± 0.07
Adult (F)0.2 ± 0.080.3 ± 0.080.2 ± 0.070.1 ± 0.09
2 years (F)N/AN/A0.3 ± 0.050.2 ± 0.10
Hatchlings0.3 ± 0.080.5 ± 0.040.5 ± 0.050.6 ± 4.60 *
* ambiguous fit.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Op den Brouw, B.; Fernandez-Rojo, M.A.; Charlton, T.; Fry, B.G.; Ikonomopoulou, M.P. Malaysian and Chinese King Cobra Venom Cytotoxicity in Melanoma and Neonatal Foreskin Fibroblasts Is Mediated by Age and Geography. Toxins 2023, 15, 549. https://doi.org/10.3390/toxins15090549

AMA Style

Op den Brouw B, Fernandez-Rojo MA, Charlton T, Fry BG, Ikonomopoulou MP. Malaysian and Chinese King Cobra Venom Cytotoxicity in Melanoma and Neonatal Foreskin Fibroblasts Is Mediated by Age and Geography. Toxins. 2023; 15(9):549. https://doi.org/10.3390/toxins15090549

Chicago/Turabian Style

Op den Brouw, Bianca, Manuel A. Fernandez-Rojo, Tom Charlton, Bryan G. Fry, and Maria P. Ikonomopoulou. 2023. "Malaysian and Chinese King Cobra Venom Cytotoxicity in Melanoma and Neonatal Foreskin Fibroblasts Is Mediated by Age and Geography" Toxins 15, no. 9: 549. https://doi.org/10.3390/toxins15090549

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