Phylogenomic insights into the diversity and evolution of Palearctic vipers.

Despite decades of molecular research, phylogenetic relationships in Palearctic vipers (genus Vipera) still essentially rely on a few loci, such as mitochondrial barcoding genes. Here we examined the diversity and evolution of Vipera with ddRAD-seq data from 33 representative species and subspecies. Phylogenomic analyses of âˆ¼ 1.1 Mb recovered nine major clades corresponding to known species/species complexes which are generally consistent with the mitochondrial phylogeny, albeit with a few deep discrepancies that highlight past hybridization events. The most spectacular case is the Italian-endemic V. walser, which is grouped with the alpine genetic diversity of V. berus in the nuclear tree despite carrying a divergent mitogenome related to the Caucasian V. kaznakovi complex. Clustering analyses of SNPs suggest potential admixture between diverged Iberian taxa (V. aspis zinnikeri and V. seoanei), and confirm that the Anatolian V. pontica corresponds to occasional hybrids between V. (ammodytes) meridionalis and V. kaznakovi. Finally, all analyzed lineages of the V. berus complex (including V. walser and V. barani) form vast areas of admixture and may be delimited as subspecies. Our study sets grounds for future taxonomic and phylogeographic surveys on Palearctic vipers, a group of prime interest for toxinological, ecological, biogeographic and conservation research.

Phylogenomics of Psammodynastes and Buhoma (Elapoidea: Serpentes), with the description of a new Asian snake family.

Asian mock vipers of the genus Psammodynastes and African forest snakes of the genus Buhoma are two genera belonging to the snake superfamily Elapoidea. The phylogenetic placements of Psammodynastes and Buhoma within Elapoidea has been extremely unstable which has resulted in their uncertain and debated taxonomy. We used ultraconserved elements and traditional nuclear and mitochondrial markers to infer the phylogenetic relationships of these two genera with other elapoids. Psammodynastes, for which a reference genome has been sequenced, were found, with strong branch support, to be a relatively early diverging split within Elapoidea that is sister to a clade consisting of Elapidae, Micrelapidae and Lamprophiidae. Hence, we allocate Psammodynastes to its own family, Psammodynastidae new family. However, the phylogenetic position of Buhoma could not be resolved with a high degree of confidence. Attempts to identify the possible sources of conflict in the rapid radiation of elapoid snakes suggest that both hybridisation/introgression during the rapid diversification, including possible ghost introgression, as well as incomplete lineage sorting likely have had a confounding role. The usual practice of combining mitochondrial loci with nuclear genomic data appears to mislead phylogeny reconstructions in rapid radiation scenarios, especially in the absence of genome scale data.

What's in a mass?

This short essay pretends to make the reader reflect on the concept of biological mass and on the added value that the determination of this molecular property of a protein brings to the interpretation of evolutionary and translational snake venomics research. Starting from the premise that the amino acid sequence is the most distinctive primary molecular characteristics of any protein, the thesis underlying the first part of this essay is that the isotopic distribution of a protein's molecular mass serves to unambiguously differentiate it from any other of an organism's proteome. In the second part of the essay, we discuss examples of collaborative projects among our laboratories, where mass profiling of snake venom PLA2 across conspecific populations played a key role revealing dispersal routes that determined the current phylogeographic pattern of the species.

Vipers of Major clinical relevance in Europe: Taxonomy, venom composition, toxicology and clinical management of human bites.

Snakebites in Europe are mostly due to bites from Viperidae species of the genus Vipera. This represents a neglected public health hazard with poorly defined incidence, morbidity and mortality. In Europe, fourteen species of "true vipers" (subfamily Viperinae) are present, eleven of which belong to the genus Vipera. Amongst these, the main medically relevant species due to their greater diffusion across Europe and the highest number of registered snakebites are six, namely: Vipera ammodytes, V. aspis, V. berus, V. latastei, V. seoanei and V. ursinii. Generally speaking, viper venom composition is characterised by many different toxin families, like phospholipases A2, snake venom serine proteases, snake venom metalloproteases, cysteine-rich secretory proteins, C-type lectins, disintegrins, haemorrhagic factors and coagulation inhibitors. A suspected snakebite is often associated with severe pain, erythema, oedema and, subsequently, the onset of an ecchymotic area around one or two visible fang marks. In the field, the affected limb should be immobilised and mildly compressed with a bandage, which can then be removed once the patient is being treated in hospital. The clinician should advise the patient to remain calm to reduce blood circulation and, therefore, decrease the spread of the toxins. In the case of pain, an analgesic therapy can be administered, the affected area can be treated with hydrogen peroxide or clean water. However, anti-inflammatory drugs and disinfection with alcohol or alcoholic substances should be avoided. For each patient, clinical chemistry and ECG are always a pre-requisite as well as the evaluation of the tetanus immunisation status and for which immunisation may be provided if needed. The treatment of any clinical complication, due to the envenomation, does not differ from treatments of emergency nature. Antivenom is recommended when signs of systemic envenomation exist or in case of advanced local or systemic progressive symptoms. Recommendations for future work concludes. The aim of this review is to support clinicians for the clinical management of viper envenomation, through taxonomic keys for main species identification, description of venom composition and mode of action of known toxins and provide a standardised clinical protocol and antivenom administration.

An ecological analysis of snakes captured by C.J.P. Ionides in eastern Africa in the mid-1900s.

Historical data can clarify ecological attributes of fauna in sites that have subsequently been altered by anthropogenic activities. We used the 1960s notebooks of pioneering East African naturalist C.J.P. Ionides to extract quantitative information on captures of 484 snakes of five species (rhinoceros vipers Bitis nasicornis, black mambas Dendroaspis polylepis, Jameson's mambas D. jamesoni, water cobras Naja annulata, and eastern forest cobras N. subfulva). High capture rates suggest high abundances of all species. The relative numbers of each species collected changed over the years and differed seasonally, reflecting targeting by Ionides. Sex ratios and age-class distributions differed among species and were affected by factors such as month of collection and time of day. Habitat use was affected by species, sex and body size: for example, arboreality became less common with increasing body size in the rhinoceros viper and black mamba, and males were found in arboreal sites more often than were females. In both D. jamesoni and D. polylepis, adult males and females were recorded together in September-October, suggesting reproductive activity at this time of year. Although fragmentary, the data from Ionides' notebooks provide a unique glimpse into ecological patterns of snakes within an African landscape half a century ago.

Protein structure of the venom in nine species of snake: from bio-compounds to possible healing agents.

Due to its various structures in bio-compounds, snake venom is the indisputable result of evolutionary stages of molecules with an increasingly complex structure, high specificity, and of great importance for medicine because of their potential. The present study proposed an underpinning examination of venom composition from nine species of venomous snakes using a useful and replicable methodology. The objective was the extension of the evaluation of protein fractions in the field up to 230 kDa to permit possible identification of some fractions that are insufficiently studied. The gel capillary electrophoresis method on the chip was performed using an Agilent 2100 bioassay with the 80 and 230-LabChip Protein kits. Interpretation of electrophoresis was performed using the Protein 2100 expert (Agilent) test software as follows: a) Protein 80 (peak size scale): 1.60, 3.5, 6.50, 15.00, 28.00, 46.00, 63.00, 95.00 kDa; b) Protein 230 (peak size scale): 4.50, 7.00, 15.00, 28.00, 46.00, 63.00, 95.00, 150.00, 240.00 kDa. The screening revealed the presence of compounds with a molecular weight greater than 80 kDa, in the case of Vipera aspis and Vipera xantina palestinae. For V. aspis, a 125 kDa molecular weight pro-coagulant protein was identified, known as being involved in the reduction of plasma clotting time without any direct activity in the fibrinogen coagulation process. The samples examined on the Protein 230-LabChip electrophoresis chip can be considered as a novelty with possible uses in medicine, requiring further approaches by advanced proteomics techniques to confirm the intimate structural features and biological properties of snake venoms.

Protein structure of the venom in nine species of snake: from bio-compounds to possible healing agents

Due to its various structures in bio-compounds, snake venom is the indisputable result of evolutionary stages of molecules with an increasingly complex structure, high specificity, and of great importance for medicine because of their potential. The present study proposed an underpinning examination of venom composition from nine species of venomous snakes using a useful and replicable methodology. The objective was the extension of the evaluation of protein fractions in the field up to 230 kDa to permit possible identification of some fractions that are insufficiently studied. The gel capillary electrophoresis method on the chip was performed using an Agilent 2100 bioassay with the 80 and 230-LabChip Protein kits. Interpretation of electrophoresis was performed using the Protein 2100 expert (Agilent) test software as follows: a) Protein 80 (peak size scale): 1.60, 3.5, 6.50, 15.00, 28.00, 46.00, 63.00, 95.00 kDa; b) Protein 230 (peak size scale): 4.50, 7.00, 15.00, 28.00, 46.00, 63.00, 95.00, 150.00, 240.00 kDa. The screening revealed the presence of compounds with a molecular weight greater than 80 kDa, in the case of Vipera aspis and Vipera xantina palestinae. For V. aspis, a 125 kDa molecular weight pro-coagulant protein was identified, known as being involved in the reduction of plasma clotting time without any direct activity in the fibrinogen coagulation process. The samples examined on the Protein 230-LabChip electrophoresis chip can be considered as a novelty with possible uses in medicine, requiring further approaches by advanced proteomics techniques to confirm the intimate structural features and biological properties of snake venoms.

Calibrating the tree of vipers under the fossilized birth-death model.

Scaling evolutionary trees to time is essential for understanding the origins of clades. Recently developed methods allow including the entire fossil record known for the group of interest and eliminated the need for specifying prior distributions for node ages. Here we apply the fossilized birth-death (FBD) approach to reconstruct the diversification timeline of the viperines (subfamily Viperinae). Viperinae are an Old World snake subfamily comprising 102 species from 13 genera. The fossil record of vipers is fairly rich and well assignable to clades due to the unique vertebral and fang morphology. We use an unprecedented sampling of 83 modern species and 13 genetic markers in combination with 197 fossils representing 28 extinct taxa to reconstruct a time-calibrated phylogeny of the Viperinae. Our results suggest a late Eocene-early Oligocene origin with several diversification events following soon after the group's establishment. The age estimates inferred with the FBD model correspond to those from previous studies that were based on node dating but FBD provides notably narrower credible intervals around the node ages. Viperines comprise two African and an Eurasian clade, but the ancestral origin of the subfamily is ambiguous. The most parsimonious scenarios require two transoceanic dispersals over the Tethys Sea during the Oligocene.

Evolutionary history of burrowing asps (Lamprophiidae: Atractaspidinae) with emphasis on fang evolution and prey selection.

Atractaspidines are poorly studied, fossorial snakes that are found throughout Africa and western Asia, including the Middle East. We employed concatenated gene-tree analyses and divergence dating approaches to investigate evolutionary relationships and biogeographic patterns of atractaspidines with a multi-locus data set consisting of three mitochondrial (16S, cyt b, and ND4) and two nuclear genes (c-mos and RAG1). We sampled 91 individuals from both atractaspidine genera (Atractaspis and Homoroselaps). Additionally, we used ancestral-state reconstructions to investigate fang and diet evolution within Atractaspidinae and its sister lineage (Aparallactinae). Our results indicated that current classification of atractaspidines underestimates diversity within the group. Diversification occurred predominantly between the Miocene and Pliocene. Ancestral-state reconstructions suggest that snake dentition in these taxa might be highly plastic within relatively short periods of time to facilitate adaptations to dynamic foraging and life-history strategies.

Retrospective Documentation of a Confirmed White-Lipped Green Pit Viper (Trimeresurus albolabris Gray, 1842) Bite in the South-Central Hills of Nepal.

This case report documents envenomation by an arboreal white-lipped green pit viper (Trimeresurus albolabris), a species found in South and Southeast Asia that causes the majority of venomous snakebites among Southeast Asian pit vipers. Clinical features vary from asymptomatic to serious coagulopathy that may progress into life-threatening or fatal hemorrhage. The proven life-threatening cases described in published literature, however, are sparse. Practically, no specific antivenom targeted to pit viper bites is available in Nepal. We report a case (managed with several non-evidence-based interventions) of noticeable coagulopathic envenomation due to confirmed T albolabris bite in Nepal. This is the first known reported case of such a bite in Nepal. This case highlights the urgent need to improve diagnosis, monitoring, and supportive care for bite victims and to study the effectiveness of Thai pit viper antivenoms for the treatment of T albolabris envenomations.

Thrombotic Microangiopathy, Hemolytic Uremic Syndrome, and Thrombotic Thrombocytopenic Purpura Following Hump-nosed Pit Viper (Genus: Hypnale) Envenoming in Sri Lanka.

Thrombotic microangiopathy (TMA), which includes the spectrum of hemolytic uremic syndrome and thrombotic thrombocytopenic purpura, is an uncommon complication of hump-nosed pit viper envenomation. We describe 4 cases of TMA following hump-nosed pit viper (Hypnale spp) bites in Sri Lanka. The first case is a typical TMA that spontaneously resolved with supportive treatments. The second and third cases are related to hemolytic uremic syndrome complicated with acute kidney injury that required hemodialysis. The fourth case is thrombotic thrombocytopenic purpura associated with acute kidney injury that required hemodialysis and therapeutic plasma exchange. For each patient we describe the circumstances of the bite, clinical features, laboratory findings, and management.

Multilocus phylogeography of the brown-spotted pitviper Protobothrops mucrosquamatus (Reptilia: Serpentes: Viperidae) sheds a new light on the diversification pattern in Asia.

Understanding the influence of geographical events and climate changes on genetic diversity is essential in explaining current patterns of genetic structure and geographic distribution of organisms. We inferred phylogenetic relationships, investigated historical demography, explored the evolutionary history, and clarified intraspecific taxonomy of Protobothrops mucrosquamatus, which is one of the commonest and most wide-ranging Asian pitvipers. A total of 184 samples from 54 localities were sequenced and analyzed for two mitochondrial gene fragments and two nuclear genes. Phylogenetic reconstruction based on mtDNA sequences revealed the existence of a minimum of five geographically structured and well-supported lineages within P. mucrosquamatus. Based on the mtDNA gene tree, and the geographic relationship between populations allied by matrilineal lineages, a complex longitudinal and latitudinal diversification pattern was uncovered in P. mucrosquamatus. The estimated date of the origin of the species (about 5.3 Ma) and divergence of the intraspecific lineages match the rapid uplifting of Qinghai-Xizang Plateau, and is also consistent with those of some other co-distributed Asian pitvipers. Formation of the two island lineages (Taiwan and Hainan) was generally congruent with the first isolation of the islands, but the two lineages showed different relationships with the continental Asian populations in comparison with some other pitvipers. Population historical demographic analyses, based on three methods, showed that all lineages have experienced slight population expansion in and around the Dali Glacial. Tests of intraspecific taxonomy indicated that no cryptic taxon is present within this widely distributed snake.

Vipers of the Middle East: A Rich Source of Bioactive Molecules.

Snake venom serves as a tool of defense against threat and helps in prey digestion. It consists of a mixture of enzymes, such as phospholipase A2, metalloproteases, and l-amino acid oxidase, and toxins, including neurotoxins and cytotoxins. Beside their toxicity, venom components possess many pharmacological effects and have been used to design drugs and as biomarkers of diseases. Viperidae is one family of venomous snakes that is found nearly worldwide. However, three main vipers exist in the Middle Eastern region: Montivipera bornmuelleri, Macrovipera lebetina, and Vipera (Daboia) palaestinae. The venoms of these vipers have been the subject of many studies and are considered as a promising source of bioactive molecules. In this review, we present an overview of these three vipers, with a special focus on their venom composition as well as their biological activities, and we discuss further frameworks for the exploration of each venom.

Variability in venom composition of European viper subspecies limits the cross-effectiveness of antivenoms.

Medically relevant cases of snakebite in Europe are predominately caused by European vipers of the genus Vipera. Systemic envenoming by European vipers can cause severe pathology in humans and different clinical manifestations are associated with different members of this genus. The most representative vipers in Europe are V. aspis and V. berus and neurological symptoms have been reported in humans envenomed by the former but not by the latter species. In this study we determined the toxicological profile of V. aspis and V. berus venoms in vivo in mice and we tested the effectiveness of two antivenoms, commonly used as antidotes, in counteracting the specific activities of the two venoms. We found that V. aspis, but not V. berus, is neurotoxic and that this effect is due to the degeneration of peripheral nerve terminals at the NMJ and is not neutralized by the two tested antisera. Differently, V. berus causes a haemorrhagic effect, which is efficiently contrasted by the same antivenoms. These results indicate that the effectiveness of different antisera is strongly influenced by the variable composition of the venoms and reinforce the arguments supporting the use polyvalent antivenoms.

Comparative Profiling of Three Atheris Snake Venoms: A. squamigera, A. nitschei and A. chlorechis.

A proteomic and transcriptomic comparative analysis of the venoms of three Atheris species (A. squamigera, A. nitschei and A. chlorechis) was carried out by size exclusion liquid chromatography, gel electrophoresis, mass spectrometry, and mRNA sequencing. The improved proteomic profiling utilised in this work was combined with transcript studies, advancing our insights into venom composition, protein distribution and inter-species variation among the three bush vipers. Crude venoms of all three samples contained at least 10-20 protein components, ranging in size from ≤ 3 to > 98 kDa. Both approaches yielded converging overall information, pointing to phospholipases, disintegrins, serine proteases and metalloproteases as the major toxin classes, which are likely to explain the local and systemic symptoms observed in envenomation by Atheris genus. Being considered as the main factors involved in the distinct venom-induced pathologies, these identified snake venom proteins are of particular interest in terms of understanding their physiological and biological function as well as for their contribution in potential medical treatments.

Conservation Below the Species Level: Suitable Evolutionarily Significant Units among Mountain Vipers (the Montivipera raddei complex) in Iran.

Northern and western mountains of Iran are among the most important biodiversity and endemism hot spots for reptiles in the Middle East. Among herpetofauna, the montivipers represent an emblematic and fragmented endemic group for which estimating their level of genetic differentiation and defining conservation priorities is urgently needed. Here, we present the most comprehensive phylogenetic study on the Montivipera raddei species group comprising all 5 known taxa, among which 3 are endemic to Iran. Based on 2 mitochondrial genes, phylogenetic and phylogeographic analyses revealed 3 major lineages each presenting very contrasting distribution areas. The Iranian montivipers are highly structured in clades showing low genetic diversity and corresponding to high altitude summits. Molecular dating revealed the role of Quaternary paleo-climatic oscillations and altitudinal movements of montivipers in shaping genetic diversity and differentiation of these sky-island taxa. In addition, the best scenario of historical biogeography allowed identifying 3 possible refugial areas in Iran most likely arising by vicariance. Based on our mitochondrial results and pending additional data, we recognize 3 candidate species among the M. raddei complex: M. raddei, Montivipera latifii, and Montivipera kuhrangica that are coherent with their geographical distribution. We propose that the most appropriate evolutionary significant units for conservation of the montivipers are represented by 13 units among which 6 are recognized as high priority. Finally, we suggest some recommendations to the IUCN as well as to the Iranian conservation policies with respect to conservation prioritization.

Diversification in vipers: Phylogenetic relationships, time of divergence and shifts in speciation rates.

Snakes of the cosmopolitan family Viperidae comprise around 329 venomous species showing a striking heterogeneity in species richness among lineages. While the subfamily Azemiopinae comprises only two species, 70% of all viper species are arranged in the subfamily Crotalinae or the "pit vipers". The radiation of the pit vipers was marked by the evolution of the heat-sensing pits, which has been suggested to be a key innovation for the successful diversification of the group. Additionally, only crotalines were able to successfully colonize the New World. Here, we present the most complete molecular phylogeny for the family to date that comprises sequences from nuclear and mitochondrial genes representing 79% of all living vipers. We also investigated the time of divergence between lineages, using six fossils to calibrate the tree, and explored the hypothesis that crotalines have undergone an explosive radiation. Our phylogenetic analyses retrieved high support values for the monophyly of the family Viperidae, subfamilies Viperinae and Crotalinae, and 22 out of 27 genera, as well as well-supported intergeneric relationships throughout the family. We were able to recover a strongly supported sister clade to the New World pit vipers that comprises Gloydius, Ovophis, Protobothrops and Trimeresurus gracilis. Our results agree in many aspects with other studies focusing on the phylogenetics of vipers, but we recover new relationships as well. Despite the addition of new sequences we were not able to resolve some of the poor supported relationships previously suggested. Time of divergence estimates suggested that vipers started to radiate around the late Paleocene to middle Eocene with subfamilies most likely dating back to the Eocene. The invasion of the New World might have taken place sometime close to the Oligocene/Miocene boundary. Diversification analyses suggested a shift in speciation rates during the radiation of a sub-clade of pit vipers where speciation rates rapidly increased but slowed down toward the present. Thus, the evolution of the loreal pits alone does not seem to explain their explosive speciation rates. We suggest that climatic and geological changes in Asia and the invasion of the New World may have also contributed to the speciation shift found in vipers.

Respiratory Effects of Sarafotoxins from the Venom of Different Atractaspis Genus Snake Species.

Sarafotoxins (SRTX) are endothelin-like peptides extracted from the venom of snakes belonging to the Atractaspididae family. A recent in vivo study on anesthetized and ventilated animals showed that sarafotoxin-b (SRTX-b), extracted from the venom of Atractaspis engaddensis, decreases cardiac output by inducing left ventricular dysfunction while sarafotoxin-m (SRTX-m), extracted from the venom of Atractaspis microlepidota microlepidota, induces right ventricular dysfunction with increased airway pressure. The aim of the present experimental study was to compare the respiratory effects of SRTX-m and SRTX-b. Male Wistar rats were anesthetized, tracheotomized and mechanically ventilated. They received either a 1 LD50 IV bolus of SRTX-b (n = 5) or 1 LD50 of SRTX-m (n = 5). The low-frequency forced oscillation technique was used to measure respiratory impedance. Airway resistance (Raw), parenchymal damping (G) and elastance (H) were determined from impedance data, before and 5 min after SRTX injection. SRTX-m and SRTX-b injections induced acute hypoxia and metabolic acidosis with an increased anion gap. Both toxins markedly increased Raw, G and H, but with a much greater effect of SRTX-b on H, which may have been due to pulmonary edema in addition to bronchoconstriction. Therefore, despite their structural analogy, these two toxins exert different effects on respiratory function. These results emphasize the role of the C-terminal extension in the in vivo effect of these toxins.

A new large species of Bitis Gray, 1842 (Serpentes: Viperidae) from the Bale Mountains of Ethiopia.

A new species of viperine viperid snake is described, Bitis harenna sp. nov. The new species is a member of the subgenus Macrocerastes based on it having three scales separating the nasal and rostral shields, and on the combination of 'divisions' of dorsal scale rows on the upper flanks and 'fusions' of rows on the lower flanks. Bitis harenna sp. nov. is distinguished from other members of the subgenus by its unique colour pattern, posterior parietal flange on the lateral wall of the braincase, and possibly by differences in scalation and head proportions. Only a single museum specimen is known, a female collected from 'Dodola' in Ethiopia probably in the late 1960s and previously identified as a possibly unusually coloured and patterned B. parviocula. A live, presumably male, specimen very closely resembling the holotype of Bitis harenna sp. nov. was photographed on the Harenna escarpment of the Bale Mountains National Park, Ethiopia in 2013, providing secure occurrence data and evidence that the holotype is not a uniquely aberrant specimen. A revised key to the species of Bitis in Ethiopia is presented. Aspects of body scalation are compared among species of the subgenus Macrocerastes and between species of Macrocerastes and Bitis, and several systematic characters are highlighted and clarified.

Multilocus phylogeny of the Asian Lance-headed pitvipers (Squamata, Viperidae, Protobothrops).

Despite the public health significance and potential applications to medical research, the evolutionary history of the Asian Lance-headed pitvipers (genus Protobothrops) remains inadequately studied. Previous research generally focused on a few selected species with no comprehensive molecular study of Protobothrops. Here, we conduct the first study to infer the phylogenetic relationships of all currently recognized species in this genus based on four mitochondrial DNA fragments and four nuclear genes. Bayesian Inference and Maximum Likelihood analyses show that within Protobothrops there are four strongly supported clades forming distinct subgroups. The first subgroup, which is sister to the other three, consists of three species, P. himalayanus, P. kaulbacki, and P. sieversorum. The second subgroup contains only P. mangshanensis. The final two subgroups, which are sister groups, include the other four and six Protobothrops species. Although our findings provide additional information on the phylogenetic relationships of the genus Protobothrops, the placement of P. dabieshanensis and P. elegans remains problematic. In addition, our molecular results indicate that P. trungkhanhensis, currently considered endemic to Vietnam, should be added to the species known from China. Our ancestral area estimation indicated that Protobothrops likely originated in southwestern China. This study improves our understanding of the evolutionary relationships among species of Asian Lance-headed pitvipers, providing a greater framework for future studies.