Background: Snakes are found on every continent in the world except Antarctica, and on smaller land masses. Being ecologically important, they also cause a large number of bites, leading to millions of deaths. Mitochondrial and nuclear gene sequences are being used to identify, characterize, and infer genetic biodiversity among different snake species. Furthermore, phylogenetics helps in inferring the relationships and evolutionary histories among these species. Black cobra is one of the four most venomous snakes in Pakistan. Four mitochondrial (ND4, Cytochrome b, 12S rRNA, and 16S rRNA) and four nuclear (C-mos, RAG-1, BDNF, and NT3) genes were used to trace diversity and infer the phylogenetic relationship of black cobra in Pakistan.

Results: Almost similar phylogenies were obtained through maximum likelihood and Bayesian inference, showing two species of cobra in Pakistan, namely, black cobra (Naja naja) and brown cobra (Naja oxiana). All Naja species were divided into three clades: black cobra (N. naja) and brown cobra (N. oxiana) cladding with different species of Naja; N. naja (Pakistan) cladding with N. naja from Nepal; and N. oxiana showed close relationship with Naja kaouthia from Thailand and Naja siamensis from Thailand.

Conclusion: It was confirmed genetically that there are two cobra species in Pakistan, i.e., black and brown cobras. This study will help in not only genetic conservation but also developing anti-venom against snake species.