The phylogeny, phylogeography, and diversification history of the westernmost Asian cobra (Serpentes: Elapidae: Naja oxiana) in the Trans-Caspian region

Biogeography of the Iranian snakes

The events of the Cenozoic era such as mountain formation caused Iran to become one of the most amazing biodiversity hotspots in the world today. This pioneering study on Iranian snake biogeography integrates historical and ecological analyses. A phylogeographic review traces speciation and dispersal, while cluster analysis with a new snake checklist assesses faunistic similarities within Iran and its surroundings. Jaccard and Sorenson indices generate similarity dendrograms, Indicator Species Analysis pinpoints regional key species, and Endemism index calculates regional endemism rates, enriching our knowledge of Iran's species diversity. Phylogeographic analyses identify four biogeographical corridors for snake ingress into Iran: the Arabian region through southwestern Iran, the Western Asian mountainous transition zone via northwestern Iran, the Turanian region into northeastern Iran, and the Indus River Valley into southeastern and eastern Iran. Dendrogram analysis divides snake fauna into three groups. The first group associates western Zagros and Khuzestan fauna with the Sahara and Arabian regions. The second group links Kopet Dagh and Turkmen Steppe fauna with the Turanian region, and Central Plateau and Baluchistan fauna with the Iranian region. The third group connects northwest highlands, Alborz and Zagros mountains, and Caspian Sea coasts with the Western Asian Mountain transition zone. The study validates broad biogeographic patterns via ecoregional associations and indicator species analysis, providing finer resolution. Species like Platyceps najadum in Caspian Hyrcanian mixed forests exemplify ecoregional alignment, while Zagros and Alborz mountains exhibit unique faunal indicators, indicating species-level divergence. Shared indicators among widespread ecoregions reflect habitat continuity; exclusive indicators emphasize regional distinctiveness. Despite endemic species prevalence, they seldom act as significant indicators due to various factors. Our research confirms the Zagros Mountains, Khuzestan Plain, Alborz Mountains, and Persian Gulf coasts as snake diversity hotspots, marked by higher species richness compared to other Iranian regions.

The Toxin Diversity, Cytotoxicity, and Enzymatic Activity of Cape Cobra (Naja nivea) Venom

"True" cobras (genus Naja) are among the venomous snakes most frequently involved in snakebite accidents in Africa and Asia. The Cape cobra (Naja nivea) is one of the African cobras of highest medical importance, but much remains to be learned about its venom. Here, we used a shotgun proteomics approach to better understand the qualitative composition of N. nivea venom and tested its cytotoxicity and protease activity as well as its effect on intracellular Ca2+ release and NO synthesis. We identified 156 venom components representing 17 protein families, with the dominant ones being three-finger toxins, mostly of the short-chain type. Two-thirds of the three-finger toxin entries identified were assigned as cytotoxins, while the remainder were categorized as neurotoxins, including short-chain, long-chain, and ancestral three-finger toxins. We also identified snake venom metalloproteinases and members of CRISP, l-amino acid oxidase, and other families. Protease activity and its effect on intracellular Ca2+ release and NO synthesis were low. Phospholipase A2 activity was surprisingly high, despite this toxin family being marginally recovered in the analyzed venom. Cytotoxicity was relevant only at higher venom concentrations, with macrophage and neuroblastoma cell lines showing the lowest viability. These results are in line with the predominantly neurotoxic envenomation symptoms caused by Cape cobra bites. The present overview of the qualitatively complex and functionally intriguing venom of N. nivea may provide insights into the pathobiochemistry of this species' venom.

Museomics help resolving the phylogeny of snowfinches (Aves, Passeridae, Montifringilla and allies)

Historical specimens from museum collections provide a valuable source of material also from remote areas or regions of conflict that are not easily accessible to scientists today. With this study, we are providing a taxon-complete phylogeny of snowfinches using historical DNA from whole skins of an endemic species from Afghanistan, the Afghan snowfinch, Pyrgilauda theresae. To resolve the strong conflict between previous phylogenetic hypotheses, we generated novel mitogenome sequences for selected taxa and genome-wide SNP data using ddRAD sequencing for all extant snowfinch species endemic to the Qinghai-Tibet Plateau (QTP) and for an extended intraspecific sampling of the sole Central and Western Palearctic snowfinch species (Montifringilla nivalis). Our phylogenetic reconstructions unanimously refuted the previously suggested paraphyly of genus Pyrgilauda. Misplacement of one species-level taxon (Onychostruthus tazcanowskii) in previous snowfinch phylogenies was undoubtedly inferred from chimeric mitogenomes that included heterospecific sequence information. Furthermore, comparison of novel and previously generated sequence data showed that the presumed sister-group relationship between M. nivalis and the QTP endemic M. henrici was suggested based on flawed taxonomy. Our phylogenetic reconstructions based on genome-wide SNP data and on mitogenomes were largely congruent and supported reciprocal monophyly of genera Montifringilla and Pyrgilauda with monotypic Onychostruthus being sister to the latter. The Afghan endemic P. theresae likely originated from a rather ancient Pliocene out-of-Tibet dispersal probably from a common ancestor with P. ruficollis. Our extended trans-Palearctic sampling for the white-winged snowfinch, M. nivalis, confirmed strong lineage divergence between an Asian and a European clade dated to 1.5 - 2.7 million years ago (mya). Genome-wide SNP data suggested subtle divergence among European samples from the Alps and from the Cantabrian mountains.

Kukri snakes Oligodon Fitzinger, 1826 of the Western Palearctic with the resurrection of Contia transcaspica Nikolsky, 1902 (Reptilia, Squamata, Colubridae)

The kukri snakes of the genus Oligodon Fitzinger, 1826 reach the westernmost limits of their distribution in Middle and Southwest Asia (Afghanistan, Iran, and Turkmenistan), and the Palearctic portions of Pakistan. In this article, we review the systematics and distribution of the two species native to this region, Oligodon arnensis (Shaw, 1802) and Oligodon taeniolatus (Jerdon, 1853) based on an integrative approach combining morphological, molecular, and species distribution modeling (SDM) data. Phylogenetic analyses recover O. taeniolatus populations from Iran and Turkmenistan in a clade with the O. arnensis species complex, rendering the former species paraphyletic relative to O. taeniolatus sensu stricto on the Indian subcontinent. To correct this, we resurrect the name Contia transcaspica Nikolsky, 1902 from the synonymy of O. taeniolatus and assign it to populations in Middle-Southwest Asia. So far, Oligodon transcaspicus comb. et stat. nov. is known only from the Köpet-Dag Mountain Range of northeast Iran and southern Turkmenistan, but SDM mapping suggests it may have a wider range. Genetic samples of O. "arnensis" from northern Pakistan are nested in a clade sister to the recently described Oligodon churahensis Mirza, Bhardwaj & Patel, 2021, and are phylogenetically separate from O. arnensis sensu stricto in south India and Sri Lanka. Based on morphological similarity, the Afghanistan and Pakistan populations are assigned to Oligodon russelius (Daudin, 1803) and we synonymize O. churahensis with this species. Our investigation leads us to remove O. taeniolatus from the snake fauna of Afghanistan, Iran, and Turkmenistan, with the consequence that only Oligodon transcaspicus comb. et stat. nov. and O. russelius are present in these countries. Additional studies are needed to resolve the taxonomy of the O. taeniolatus and O. arnensis species complexes on the Indian subcontinent, and an updated key for both groups is provided.

Medically important snakes and snakebite envenoming in Iran

Snakebite is a common health condition in Iran with a diverse snake fauna, especially in tropical southern and mountainous western areas of the country with plethora of snake species. The list of medically important snakes, circumstances and effects of their bite, and necessary medical care require critical appraisal and should be updated regularly. This study aims to review and map the distributions of medically important snake species of Iran, re-evaluate their taxonomy, review their venomics, describe the clinical effects of envenoming, and discuss medical management and treatment, including the use of antivenom. Nearly 350 published articles and 26 textbooks with information on venomous and mildly venomous snake species and snakebites of Iran, were reviewed, many in Persian (Farsi) language, making them relatively inaccessible to an international readership. This has resulted in a revised updated list of Iran's medically important snake species, with taxonomic revisions of some, compilation of their morphological features, remapping of their geographical distributions, and description of species-specific clinical effects of envenoming. Moreover, the antivenom manufactured in Iran is discussed, together with treatment protocols that have been developed for the hospital management of envenomed patients.

Multilevel Comparison of Indian Naja Venoms and Their Cross-Reactivity with Indian Polyvalent Antivenoms

Snake envenoming is caused by many biological species, rather than a single infectious agent, each with a multiplicity of toxins in their venom. Hence, developing effective treatments is challenging, especially in biodiverse and biogeographically complex countries such as India. The present study represents the first genus-wide proteomics analysis of venom composition across Naja species (N. naja, N. oxiana, and N. kaouthia) found in mainland India. Venom proteomes were consistent between individuals from the same localities in terms of the toxin families present, but not in the relative abundance of those in the venom. There appears to be more compositional variation among N. naja from different locations than among N. kaouthia. Immunoblotting and in vitro neutralization assays indicated cross-reactivity with Indian polyvalent antivenom, in which antibodies raised against N. naja are present. However, we observed ineffective neutralization of PLA2 activities of N. naja venoms from locations distant from the source of immunizing venoms. Antivenom immunoprofiling by antivenomics revealed differential antigenicity of venoms from N. kaouthia and N. oxiana, and poor reactivity towards 3FTxs and PLA2s. Moreover, there was considerable variation between antivenoms from different manufacturers. These data indicate that improvements to antivenom manufacturing in India are highly desirable.

Description of a New Cobra (Naja Laurenti, 1768; Squamata, Elapidae) from China with Designation of a Neotype for Naja atra

Taxonomic frameworks for medically important species such as cobras (genus Naja Laurenti, 1768; Squamata, Elapidae) are essential for the medical treatment of snake bites and accurate antivenin development. In this paper, we described the former N. kaouthia populations recorded from China as a new species and designated a neotype for N. atra-based morphological and mitochondrial phylogenetic analysis. The new species N. fuxisp. nov. was morphologically diagnosed from N. kaouthia by (1) regular single narrow crossband present on the middle and posterior parts of the dorsum (3-15, 7.9 ± 2.7, n = 32) and the dorsal surface of the tail (1-6, 4.2 ± 1.1, n = 32) of both adults and juveniles, buff-colored with dark fringes on both edges, vs. South Asian populations (n = 39) and Southeast Asian populations (n = 35) without cross bands, with irregular cross bands or multiple light-colored crossbands pairs, or densely woven lines; (2) small scales between the posterior chin shields, usually three (40%) or two (37%), rarely four (13%), or one (10%) (n = 30) vs. mostly one (81%) and rarely two (19%) (n = 28); (3) ventrals 179-205 (195.4 ± 6.7, n = 33) vs. South Asian populations 179-199 (188.7 ± 5.9, n = 12); Southeast Asian populations 168-186 (177.8 ± 4.9, n = 18). Phylogenetically, the new species forms an independent sister clade to the clade including N. atra, N. kaouthia, N. oxiana and N. sagittifera. Furthermore, the subspecies N. naja polyocellata should be resurrected and recognized as a full species, N. polyocellatacomb. nov., and the subspecies N. sumatrana miolepis should be resurrected.

Biogeographic Inferences on the Evolutionary History of the King Cobra (Ophiophagus hannah, Cantor 1836) Species Complex

King cobra (Ophiophagus hannah) is a snake widely distributed through southeastern tropical Asia, but in two separate subpopulations: one located in the Western Ghats (western Indian Peninsula) and the other much more extensive, ranging between the southern slopes of the Himalayas, Assam, Indochina to southeastern China. Similarly, it also appears in numerous tropical archipelagos such as Indonesia, the Philippines, and the Andaman Islands, but surprisingly it is absent from other large islands like Sri Lanka and Taiwan. In this study, we evaluated how climate could be shaping the distribution of this snake and estimated the future distribution of the species utilizing ecological niche modelling. To evaluate the effect of paleoclimatic conditions on the genetic structure of this species we performed Bayesian phylogenetic analysis under a molecular clock using mitochondrial DNA. Our analyses indicated that the current distribution of O. hannah is strongly influenced by the availability of humid climate conditions. King cobras have a long evolutionary history reflected in the appearance of four main mitochondrial lineages before the Pliocene (the Western Ghats, southeastern mainland Asia, Luzon, and Indonesia), congruently with paleoclimatic models that indicated the availability of suitable conditions for this species in these refugia during the glacial cycles. Climate history could explain the absence of O. hannah in Sri Lanka and Taiwan due to the absence of suitable climatic corridors when these islands were connected to the mainland (20,000 years ago). Future projections (2050‒2070) did not suggest significant range shifts in the region, even considering the worst global warming scenarios.

Snake Venomics: Fundamentals, Recent Updates, and a Look to the Next Decade

Venomic research, powered by techniques adapted from proteomics, transcriptomics, and genomics, seeks to unravel the diversity and complexity of venom through which knowledge can be applied in the treatment of envenoming, biodiscovery, and conservation. Snake venom proteomics is most extensively studied, but the methods varied widely, creating a massive amount of information which complicates data comparison and interpretation. Advancement in mass spectrometry technology, accompanied by growing databases and sophisticated bioinformatic tools, has overcome earlier limitations of protein identification. The progress, however, remains challenged by limited accessibility to samples, non-standardized quantitative methods, and biased interpretation of -omic data. Next-generation sequencing (NGS) technologies enable high-throughput venom-gland transcriptomics and genomics, complementing venom proteomics by providing deeper insights into the structural diversity, differential expression, regulation and functional interaction of the toxin genes. Venomic tissue sampling is, however, difficult due to strict regulations on wildlife use and transfer of biological materials in some countries. Limited resources for techniques and funding are among other pertinent issues that impede the progress of venomics, particularly in less developed regions and for neglected species. Genuine collaboration between international researchers, due recognition of regional experts by global organizations (e.g., WHO), and improved distribution of research support, should be embraced.

Venomics of the Enigmatic Andaman Cobra (Naja sagittifera) and the Preclinical Failure of Indian Antivenoms in Andaman and Nicobar Islands

The Andaman and Nicobar Islands are an abode to a diversity of flora and fauna, including the many endemic species of snakes, such as the elusive Andaman cobra (Naja sagittifera). However, the ecology and evolution of venomous snakes inhabiting these islands have remained entirely uninvestigated. This study aims to bridge this knowledge gap by investigating the evolutionary history of N. sagittifera and its venom proteomic, biochemical and toxicity profile. Phylogenetic reconstructions confirmed the close relationship between N. sagittifera and the Southeast Asian monocellate cobra (N. kaouthia). Overlooking this evolutionary history, a polyvalent antivenom manufactured using the venom of the spectacled cobra (N. naja) from mainland India is used for treating N. sagittifera envenomations. Comparative evaluation of venoms of these congeners revealed significant differences in their composition, functions and potencies. Given the close phylogenetic relatedness between N. sagittifera and N. kaouthia, we further assessed the cross-neutralising efficacy of Thai monovalent N. kaouthia antivenom against N. sagittifera venoms. Our findings revealed the inadequate preclinical performance of the Indian polyvalent and Thai monovalent antivenoms in neutralising N. sagittifera venoms. Moreover, the poor efficacy of the polyvalent antivenom against N. naja venom from southern India further revealed the critical need to manufacture region-specific Indian antivenoms.

Climate change and the increase of human population will threaten conservation of Asian cobras

Asian cobras (genus Naja) are venomous snakes distributed from the Middle East to Southeast Asia. Because cobras often live near humans settlements, they are responsible for a large part of snakebite incidents and as such pose a challenge for public health systems. In the light of growing human populations, correctly mapping the present and future ranges of Asian cobras is therefore important for both biological conservation and public health management. Here, we mapped the potential climatic niches of ten Asian cobra species for both the present and the future, with the aim to quantify changes in climate and human population densities relative to their current and future ranges. Our analyses reveal that cobras that are adapted to dry climates and inhabit islands have narrow climatic niches, while those of mainland species with larger geographic ranges are much wider. We also found a higher degree of fragmentation of future cobra distributions; within the next 50 years, Asian cobras will lose an average of around 60% of their current suitable climatic range. In the near future, Naja mandalayensis, N. sputatrix, N. samarensis, and N. philippinensis are likely to have no accessible suitable climate space left. Besides, a further increase of human populations in this region may also exponentially accelerate the effects of anthropogenic impacts. Solutions for conservation may involve awareness and appropriate use of law to overcome the rate of habitat degradation and the increase of animal trade of Asian cobras, while promoting investment on health systems to avoid snakebite fatalities.