Evolutionary relationships among the true vipers (Reptilia: Viperidae) inferred from mitochondrial DNA sequences

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.

Vertebral morphology and intracolumnar variation of the iconic African viperid snake Atheris (Serpentes, Viperidae)

We here provide a detailed description of the vertebral morphology of the African arboreal viperid snakes of the genus Atheris. Vertebrae of three different species of the genus, i.e., Atheris desaixi, Atheris hispida, and Atheris katangensis, were investigated via the aid of μCT (micro-computed tomography) scanning. We describe several vertebrae from different regions of the vertebral column for all three species, starting from the atlas-axis complex to the caudal tip, in order to demonstrate important differences regarding the intracolumnar variation. Comparison of these three species shows an overall similar general morphology of the trunk vertebrae among the Atheris species. We extensively compare Atheris with other known viperids. As the sole arboreal genus of Viperinae the prehensile nature of the tail of Atheris is reflected in its caudal vertebral morphology, which is characterized by a high number of caudal vertebrae but also robust and anteroventrally oriented pleurapophyses as a skeletal adaptation, linked with the myology of the tail, to an arboreal lifestyle. We anticipate that the extensive figuring of these viperid specimens will also aid identifications in paleontology.

Phylogeography and Ecological Niche Modeling of the Alashan Pit Viper (Gloydius cognatus; Reptilia, Viperidae) in Northwest China and Adjacent Areas

The joint impacts of historical geological events and Quaternary climatic oscillations in Northwest China on species evolution have been examined extensively in plant under a phylogeographic perspective. However, animal phylogeographic analyses in this region are still limited. The Alashan pit viper, Gloydius cognatus, occurs primarily in arid Northwest China and adjacent areas. Based on variation at two mtDNA genes (ND4 and Cytb) in 27 individuals representing 24 populations, the spatial genetic structure and demographic history of G. cognatus were examined across its geographic range. Phylogenetic analyses revealed two well-supported allopatric clades (each with two distinct subclades/lineages), distributed across the southern (Qaidam Basin, Lanzhou Basin, and Zoige Basin [S1]; Loess Plateau [S2]) and northern (Ily Basin [N1]; Junggar Basin and Mongolian Plateau [N2]) regions. AMOVA analysis demonstrated that over 76% of the observed genetic variation was related to these lineage splits, indicating substantial genetic differentiation among the four lineages. A strong pattern of isolation-by-distance across the sampling populations suggested that geographic distance principally shaped the genetic structure. The four lineages diverged by 0.9-2.2% for the concatenated data, which were estimated to have coalesced ~1.17 million years ago (Mya), suggesting that the expansions of the Badain Jaran, Tengger, and Mu Us deserts during the Xixiabangma glaciation likely interrupted gene flow and triggered the observed divergence in the southern and northern regions. Subsequently, the early Pleistocene integration of the Yellow River and associated deserts expansion promoted the differentiation of S1 and S2 lineages (~0.9 Mya). Both mitochondrial evidence and ecological niche modeling (ENM) reject the signature of demographic and range contractions during the LGM for G. cognatus. In addition, ENM predicts that the suitable habitat of G. cognatus will contract in the future. As such, the conservation and management of ESUs should be a priority. Our findings provide the first insights on the lineage diversification and population dynamics of the Alashan pit viper in relation to geological history and Pleistocene climatic oscillations in arid Northwest China.

Venomics of the milos viper (Macrovipera schweizeri) unveils patterns of venom composition and exochemistry across blunt-nosed viper venoms

Introduction: Snakebite is a neglected tropical disease and a globally important driver of death and morbidity. Vipers of the genus Macrovipera (Viperidae: Viperinae) are among the snakes of higher medical importance in the Old World. Despite the medical relevance of Macrovipera venoms, the knowledge regarding them is heterogeneously distributed with virtually all works conducted so far focusing on subspecies of Macrovipera lebetinus, while other species within the genus are largely overlooked. Here we present the first proteomic evaluation of the venom from the Greek endemic Milos viper (Macrovipera schweizeri). In line with clinical symptoms typically elicited by Macrovipera envenomations, Milos viper venom primarily comprises coagulotoxic and cytotoxic protein families, such as metalloproteinases (svMP) and serine proteases (svSP). Methods: We conducted comparative bioactivity assays on venoms from M. schweizeri and the M. lebetinus subspecies M. lebetinus cernovi, M. lebetinus obtusa, and M. lebetinus turanica, and showed that they all exhibit similarities in levels of cytotoxicity proteolytic activity, and inhibition of prokaryotic growth. Lastly, we compared Macrovipera venom profiles by 1D-SDS-PAGE and RP-HPLC, as well as our proteomic data with previously published Macrovipera venom proteomes. Results and discussion: The analyzes performed to reveal that a general venom profile seems to be conserved across blunt-nosed vipers, and that, M. schweizeri envenomations, similarly to those caused by other blunt-nosed vipers, are able to cause significant tissue damage. The present work represents an important starting point for the development of comparative studies across the full taxonomic range of the genus Macrovipera and can potentially help optimize the treatment of envenomations caused by M. schweizeri.

Total-evidence phylogeny and evolutionary morphology of New World pitvipers (Serpentes: Viperidae: Crotalinae)

Crotalines (pitvipers) in the Americas are distributed from southern Canada to southern Argentina, and are represented by 13 genera and 163 species that constitute a monophyletic group. Their phylogenetic relationships have been assessed mostly based on DNA sequences, while morphological data have scarcely been used for phylogenetic inquiry. We present a total-evidence phylogeny of New World pitvipers, the most taxon/character comprehensive phylogeny to date. Our analysis includes all genera, morphological data from external morphology, cranial osteology and hemipenial morphology, and DNA sequences from mitochondrial and nuclear genes. We performed analyses with parsimony as an optimality criterion, using different schemes for character weighting. We evaluated the contribution of the different sources of characters to the phylogeny through analyses of reduced datasets and calculation of weighted homoplasy and retention indexes. We performed a morphological character analysis to identify synapomorphies for the main clades. In terms of biogeography, our results support a single colonization event of the Americas by pitvipers, and a cladogenetic event into a Neotropical clade and a North American/Neotropical clade. The results also shed light on the previously unstable position of some taxa, although they could not sufficiently resolve the position of Bothrops lojanus, which may lead to the paraphyly of either Bothrops or Bothrocophias. The morphological character analyses demonstrated that an important phylogenetic signal is contained in characters related to head scalation, the jaws and the dorsum of the skull, and allowed us to detect morphological convergences in external morphology associated with arboreality.

Anti-Metalloproteases: Production and Characterization of Polyclonal IgG Anti-F2 Fraction Antibodies Purified from the Venom of the Snake Bitis arietans

Bitis arietans is a medically important snake found in Sub-Saharan Africa. The envenomation is characterized by local and systemic effects, and the lack of antivenoms aggravates the treatment. This study aimed to identify venom toxins and develop antitoxins. The F2 fraction obtained from Bitis arietans venom (BaV) demonstrated the presence of several proteins in its composition, including metalloproteases. Titration assays carried out together with the immunization of mice demonstrated the development of anti-F2 fraction antibodies by the animals. The determination of the affinity of antibodies against different Bitis venoms was evaluated, revealing that only BaV had peptides recognized by anti-F2 fraction antibodies. In vivo analyses demonstrated the hemorrhagic capacity of the venom and the effectiveness of the antibodies in inhibiting up to 80% of the hemorrhage and 0% of the lethality caused by BaV. Together, the data indicate: (1) the prevalence of proteins that influence hemostasis and envenomation; (2) the effectiveness of antibodies in inhibiting specific activities of BaV; and (3) isolation and characterization of toxins can become crucial steps in the development of new alternative treatments. Thus, the results obtained help in understanding the envenoming mechanism and may be useful for the study of new complementary therapies.

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

We conducted a comprehensive analysis of the phylogenetic, phylogeographic, and demographic relationships of Caspian cobra (Naja oxiana; Eichwald, 1831) populations based on a concatenated dataset of two mtDNA genes (cyt b and ND4) across the species' range in Iran, Afghanistan, and Turkmenistan, along with other members of Asian cobras (i.e., subgenus Naja Laurenti, 1768). Our results robustly supported that the Asiatic Naja are monophyletic, as previously suggested by other studies. Furthermore, N. kaouthia and N. sagittifera were recovered as sister taxa to each other, and in turn sister clades to N. oxiana. Our results also highlighted the existence of a single major evolutionary lineage for populations of N. oxiana in the Trans-Caspian region, suggesting a rapid expansion of this cobra from eastern to western Asia, coupled with a rapid range expansion from east of Iran toward the northeast. However, across the Iranian range of N. oxiana, subdivision of populations was not supported, and thus, a single evolutionary significant unit is proposed for inclusion in future conservation plans in this region.

Macroevolutionary variation and environmental correlates of scalation traits in Eurasian vipers (Serpentes: Viperinae)

Understanding how phenotypic variation across species is shaped by the combination of shared evolutionary history and environmental factors is key to elucidating the processes that underlie biodiversity. In reptiles, morphological traits have traditionally been used to delimit species and make systematic inferences. Recent studies highlight the possibility that phenotypic variation, particularly in scalation traits, might instead be driven by environmental factors and therefore not reflect the phylogenetic relationships among species. In this study, we combined morphological and ecological data in a macroevolutionary framework, in order to describe the morphological variation across species of Eurasian vipers (Serpentes: Viperinae), investigate the phylogenetic structure of scalation traits and test the contribution of environmental factors in shaping morphological patterns. We found considerable variation in all examined traits, which, in most cases, agreed with the phylogenetic relationships among species, reinforcing their usefulness for taxonomic inferences. Interestingly, however, the number of ventral scales exhibited lower phylogenetic signal and a tight association with environmental factors of geographical ranges, suggesting potential adaptive or developmental sources of variation in the trait. This is the first comparative study of macroevolutionary variation in scalation traits in Eurasian vipers, validating the use of most of them for systematic inferences, but also indicating possible environmental factors that might shape phenotypic variation across species.

Dagestan blunt-nosed viper, Macrovipera lebetina obtusa (Dwigubsky, 1832), venom. Venomics, antivenomics, and neutralization assays of the lethal and toxic venom activities by anti-Macrovipera lebetina turanica and anti-Vipera berus berus antivenoms

We have applied a combination of venomics, in vivo neutralization assays, and in vitro third-generation antivenomics analysis to assess the preclinical efficacy of the monospecific anti-Macrovipera lebetina turanica (anti-Mlt) antivenom manufactured by Uzbiopharm® (Uzbekistan) and the monospecific anti-Vipera berus berus antivenom from Microgen® (Russia) against the venom of Dagestan blunt-nosed viper, Macrovipera lebetina obtusa (Mlo). Despite their low content of homologous (anti-Mlt, 5–10%) or para-specific (anti-Vbb, 4–9%) F(ab')₂ antibody fragments against M. l. obtusa venom toxins, both antivenoms efficiently recognized most components of the complex venom proteome's arsenal, which is made up of toxins derived from 11 different gene families and neutralized, albeit at different doses, key toxic effects of M. l. obtusa venom, i.e., in vivo lethal and hemorrhagic effects in a murine model, and in vitro phospholipase A₂, proteolytic and coagulant activities. The calculated lethality neutralization potencies for Uzbiopharm® anti-Mlt and anti-Vbb Microgen® antivenoms were 1.46 and 1.77 mg/mL, indicating that 1 mL of Uzbiopharm® and Microgen® antivenoms may protect mice from 41 to 50 LD₅₀s of Mlo venom, respectively. The remarkable degree of conservation of immunogenic determinants between species of the clades of European and Oriental viper, which evolved geographically segregated since the early Miocene, suggests an eventual window of opportunity for the treatment of envenomings by Eurasian snakes. Clearly, the rational use of heterologous antivenoms requires establishing their para-specificity landscapes. This paper illustrates the analytical power of combining in vitro and in vivo preclinical quantitative assays toward this goal.

•

Efficacy against M. l. obtusa venom by two antivenoms was investigated.

•

Quantification of lethality neutralizing antibodies was assessed.

•

Anti-Vipera berus antivenom showed paraspecificity against M. l. obtusa venom.

•

This study provides hints as how to improve the potency of the antivenoms sampled.

Genetic evidence for occurrence of Macrovipera razii (Squamata, Viperidae) in the central Zagros region, Iran

This study presents the first molecular evidence of Macrovipera razii from central Zagros, more than 300 km north-west of its prior records in southern Iran. Molecular analyses based on mitochondrial cytochrome b sequences identified the individuals from central Zagros as a lineage of M. razii. Specimens from the new localities are separated by a genetic distance of 1.46% from the known populations of M. razii. The results extend the known distribution range of M. razii as an endemic species of Iran.

Venomous Landmines: Clinical Implications of Extreme Coagulotoxic Diversification and Differential Neutralization by Antivenom of Venoms within the Viperid Snake Genus Bitis

The genus Bitis comprises 17 snake species that inhabit Africa and the Arabian Peninsula. They are responsible for a significant proportion of snakebites in the region. The venoms of the two independent lineages of giant Bitis (B. arietans and again in the common ancestor of the clade consisting of B. gabonica, B. nasicornis, B. parviocula and B. rhinoceros) induce an array of debilitating effects including anticoagulation, hemorrhagic shock and cytotoxicity, whilst the dwarf species B. atropos is known to have strong neurotoxic effects. However, the venom effects of the other species within the genus have not been explored in detail. A series of coagulation assays were implemented to assess the coagulotoxic venom effects of fourteen species within the genus. This study identified procoagulant venom as the ancestral condition, retained only by the basal dwarf species B. worthingtoni, suggesting anticoagulant venom is a derived trait within the Bitis genus and has been secondarily amplified on at least four occasions. A wide range of anticoagulant mechanisms were identified, such as coagulant and destructive activities upon fibrinogen in both giant and dwarf Bitis and the action of inhibiting the prothrombinase complex, which is present in a clade of dwarf Bitis. Antivenom studies revealed that while the procoagulant effects of B. worthingtoni were poorly neutralized, and thus a cause for concern, the differential mechanisms of anticoagulation in other species were all well neutralized. Thus, this study concludes there is a wide range of coagulotoxic mechanisms which have evolved within the Bitis genus and that clinical management strategies are limited for the procoagulant effects of B. worthingtoni, but that anticoagulant effects of other species are readily treated by the South African polyvalent antivenom. These results therefore have direct, real-work implications for the treatment of envenomed patients.

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.

Vipera berus berus Venom from Russia: Venomics, Bioactivities and Preclinical Assessment of Microgen Antivenom

The common European adder, Vipera berus berus, is a medically relevant species, which is widely distributed in Russia and thus, is responsible for most snakebite accidents in Russia. We have investigated the toxic and enzymatic activities and have determined the proteomic composition of its venom. Phospholipases A₂ (PLA₂, 25.3% of the venom proteome), serine proteinases (SVSP, 16.2%), metalloproteinases (SVMP, 17.2%), vasoactive peptides (bradykinin-potentiating peptides (BPPs), 9.5% and C-type natriuretic peptides (C-NAP, 7.8%), cysteine-rich secretory protein (CRISP, 8%) and L-amino acid oxidase (LAO, 7.3%) represent the major toxin classes found in V. b. berus (Russia) venom. This study was also designed to assess the in vivo and in vitro preclinical efficacy of the Russian Microgen antivenom in neutralizing the main effects of V. b. berus venom. The results show that this antivenom is capable of neutralizing the lethal, hemorrhagic and PLA₂ activities. Third-generation antivenomics was applied to quantify the toxin-recognition landscape and the maximal binding capacity of the antivenom for each component of the venom. The antivenomics analysis revealed that 6.24% of the anti-V. b. berus F(ab’)₂ molecules fraction are toxin-binding antibodies, 60% of which represent clinically relevant antivenom molecules.

Next-generation sequencing yields the complete mitochondrial genome of the endangered Milos viper Macrovipera schweizeri (Reptilia, Viperidae)

The Milos viper, Macrovipera schweizeri, is an endangered viperid snake found on four Aegean islands (Greece). Its complete mitochondrial genome, the first reported for the genus Macrovipera, was assembled through next-generation sequencing. Its total length is 17,152 bp and includes 22 tRNAs, two ribosomal RNA genes, 13 protein-coding genes and two control regions, showing the typical gene-arrangement for Viperidae. Eight tRNAs and ND3 are encoded on the light strand, while all other genes are encoded on the heavy strand. A mitogenomic phylogeny that included Macrovipera schweizeri and 13 other viperid genera returned an unresolved relationship among the genera Macrovipera, Daboia and Vipera.

Daboia (Vipera) palaestinae envenomation in horses: Clinical and hematological signs, risk factors for mortality and construction of a novel severity scoring system

Daboia palaestinae is the most common venomous snake in Israel and an important cause of envenomations in humans and animals. Although specific antivenom is produced from horses, little documentation exists regarding the characteristics of envenomed horses. This survey was constructed to describe the clinical, hematological and biochemical characteristics of D. palaestinae envenomation in horses, to identify risk factors for fatality and to construct a prognostic snakebite severity scale (SSS) to be used by veterinarians in the field. Data regarding 123 equine snakebite cases were collected over four years. Clinical signs most frequently described were local swelling (100%), tachycardia (67%), increased respiratory effort (62%), bite marks (44%) and reduced appetite (50%). Blood samples were collected at initial presentation from 23 of these horses. Hematological and biochemical abnormalities described in the majority of the sampled cases included leukocytosis (52%), neutrophilia (78%), increased creatine phosphokinase (77%), hypoproteinemia (50%), decreased total carbon dioxide (54%), and hyperphosphatemia (54%). The overall case fatality rate was 20.3% (25/123 horses). Risk factors associated with increased mortality (P < 0.05) were body weight below 200 kg, hind limb envenomation, severity of local swelling, mental status, respiratory effort, bleeding disorders, dehydration, basophilia, monocytosis, bilirubinemia, hypocalcemia and hyperphosphatemia. These risk factors were used to construct an SSS adapted for the assessment of envenomed horses in the field, consisting of seven dichotomous parameters. A cutoff value of 1.5 predicted mortality with 95% sensitivity and 68% specificity. In the future, this tool could assist clinicians with case management and communication with owners.

An integrative systematic revision and biogeography of Rhynchocalamus snakes (Reptilia, Colubridae) with a description of a new species from Israel

BACKGROUND

The colubrid snakes of the genus Rhynchocalamus are seldom studied and knowledge of their ecology and life history is scarce. Three species of Rhynchocalamus are currently recognized, R. satunini (from Turkey eastwards to Iran), R. arabicus (Yemen and Oman), and R. melanocephalus (from the Sinai Peninsula northwards to Turkey). All are slender, secretive, mainly nocturnal and rare fossorial snakes. This comprehensive study is the first to sample all known Rhynchocalamus species in order to review the intra-generic phylogenetic relationships and historical biogeography of the genus.

METHODS

We revised the systematics of Rhynchocalamus using an integrative approach and evaluated its phylogeography. The phylogenetic position within the Colubridae and the phylogenetic relationships within the genus were inferred using 29 individuals belonging to the three known species, with additional sampling of two other closely-related genera, Muhtarophis and Lytorhynchus. We analysed three mitochondrial (12S, 16S, cytb) and one nuclear (c-mos) gene fragments. Phylogenetic trees were reconstructed using maximum likelihood and Bayesian inference methods; the latter method also used to provide the first time-calibrated molecular phylogeny of the genus. We generated a nuclear network and carried out a topology test and species delimitation analysis. Morphological comparisons were used to differentiate among species and to describe a new species from Israel. The studied material was comprised of 108 alcohol-preserved specimens, 15 photographs, and data from the literature for the examination of 17 mensural, 14 meristic, and two categorical characters.

RESULTS

The molecular results support Rhynchocalamus as monophyletic, and as having split from its sister genus Lytorhynchus during the Late Oligocene. The three recognized species of Rhynchocalamus comprise four independently evolving groups. The molecular results reveal that the genus began to diverge during the Middle Miocene. We revealed that the best-studied species, R. melanocephalus, is paraphyletic. A population, formally ascribed to this species, from the Negev Mountain area in southern Israel is phylogenetically closer to R. arabicus from Oman than to the northern populations of the species from Israel, Syria and Turkey. Herein we describe this population as a new species: Rhynchocalamus dayanaesp. nov.

DISCUSSION

We identify four species within Rhynchocalamus: R. satunini, R. arabicus, R. melanocephalus, and R. dayanaesp. nov., the latter, to the best of our knowledge, is endemic to southern Israel. The onset of Rhynchocalamus diversification is very old and estimated to have occurred during the Middle Miocene, possibly originating in the Levant region. Radiation probably resulted from vicariance and dispersal events caused by continuous geological instability, sea-level fluctuations and climatic changes within the Levant region.

Morphological and functional alteration of erythrocyte ghosts and giant unilamellar vesicles caused by Vipera latifi venom

Snake bites are an endemic public health problem in Iran, both in rural and urban area. Viper venom as a hemolytic biochemical "cocktail" of toxins, primarily cause to the systemic alteration of blood cells. In the sixties and seventies, human erythrocytes were extensively studied, but the mechanical and chemical stresses commonly exerted on red blood cells continue to attract interest of scientists for the study of membrane structure and function. Here, we monitor the effect of Vipera latifi venom on human erythrocytes ghost membranes using phase contrast and fluorescent microscopy and changes in ATPase activity under snake venom influence in vitro. The ion pumps [Na⁺,K⁺]-ATPase and (Ca²⁺+Mg²⁺)-ATPase plays a pivotal role in the active transport of certain cations and maintenance of intracellular electrolyte homeostasis. We also describe the interaction of Vipera latifi (VL) venom with giant unilamellar vesicles (GUVs) composed of the native phospholipid mixtures visualized by the membrane fluorescence probe, ANS, used to assess the state of membrane and specifically mark the phospholipid domains.

A Species-Level Phylogeny of Extant Snakes with Description of a New Colubrid Subfamily and Genus

BACKGROUND

With over 3,500 species encompassing a diverse range of morphologies and ecologies, snakes make up 36% of squamate diversity. Despite several attempts at estimating higher-level snake relationships and numerous assessments of generic- or species-level phylogenies, a large-scale species-level phylogeny solely focusing on snakes has not been completed. Here, we provide the largest-yet estimate of the snake tree of life using maximum likelihood on a supermatrix of 1745 taxa (1652 snake species + 7 outgroup taxa) and 9,523 base pairs from 10 loci (5 nuclear, 5 mitochondrial), including previously unsequenced genera (2) and species (61).

RESULTS

Increased taxon sampling resulted in a phylogeny with a new higher-level topology and corroborate many lower-level relationships, strengthened by high nodal support values (> 85%) down to the species level (73.69% of nodes). Although the majority of families and subfamilies were strongly supported as monophyletic with > 88% support values, some families and numerous genera were paraphyletic, primarily due to limited taxon and loci sampling leading to a sparse supermatrix and minimal sequence overlap between some closely-related taxa. With all rogue taxa and incertae sedis species eliminated, higher-level relationships and support values remained relatively unchanged, except in five problematic clades.

CONCLUSION

Our analyses resulted in new topologies at higher- and lower-levels; resolved several previous topological issues; established novel paraphyletic affiliations; designated a new subfamily, Ahaetuliinae, for the genera Ahaetulla, Chrysopelea, Dendrelaphis, and Dryophiops; and appointed Hemerophis (Coluber) zebrinus to a new genus, Mopanveldophis. Although we provide insight into some distinguished problematic nodes, at the deeper phylogenetic scale, resolution of these nodes may require sampling of more slowly-evolving nuclear genes.

Phylogeny and diversification of mountain vipers (Montivipera, Nilson et al., 2001) triggered by multiple Plio-Pleistocene refugia and high-mountain topography in the Near and Middle East

The Near and Middle East is a hotspot of biodiversity, but the region remains underexplored at the level of genetic biodiversity. Here, we present an extensive molecular phylogeny of the viperid snake genus Montivipera, including all known taxa. Based on nuclear and mitochondrial data, we present novel insights into the phylogeny of the genus and review the status of its constituent species. Maximum likelihood methods revealed a montane origin of Montivipera at 12.3Mya. We then analyzed factors of mountain viper diversity. Our data support substantial changes in effective population size through Plio-Pleistocene periods. We conclude that climatic oscillations were drivers of allopatric speciation, and that mountain systems of the Near and Middle East have strongly influenced the evolution and survival of taxa, because climatic and topographical heterogeneities induced by mountains have played a crucial role as filters for dispersal and as multiple refugia. The wide diversity of montane microhabitats enabled mountain vipers to retain their ecological niche during climatic pessima. In consequence the varied geological and topographical conditions between refugia favoured genetic isolation and created patterns of species richness resulting in the formation of neoendemic taxa. Our data support high concordance between geographic distributions of Montivipera haplotypes with putative plant refugia.

African adders: partial characterization of snake venoms from three Bitis species of medical importance and their neutralization by experimental equine antivenoms

BACKGROUND

An alarming number of fatal accidents involving snakes are annually reported in Africa and most of the victims suffer from permanent local tissue damage and chronic disabilities. Envenomation by snakes belonging to the genus Bitis, Viperidae family, are common in Sub-Saharan Africa. The accidents are severe and the victims often have a poor prognosis due to the lack of effective specific therapies. In this study we have biochemically characterized venoms from three different species of Bitis, i.e., Bitis arietans, Bitis gabonica rhinoceros and Bitis nasicornis, involved in the majority of the human accidents in Africa, and analyzed the in vitro neutralizing ability of two experimental antivenoms.

METHODOLOGY/PRINCIPAL FINDINGS

The data indicate that all venoms presented phospholipase, hyaluronidase and fibrinogenolytic activities and cleaved efficiently the FRET substrate Abz-RPPGFSPFRQ-EDDnp and angiotensin I, generating angiotensin 1-7. Gelatinolytic activity was only observed in the venoms of B. arietans and B. nasicornis. The treatment of the venoms with protease inhibitors indicated that Bitis venoms possess metallo and serinoproteases enzymes, which may be involved in the different biological activities here evaluated. Experimental antivenoms produced against B. arietans venom or Bitis g. rhinoceros plus B. nasicornis venoms cross-reacted with the venoms from the three species and blocked, in different degrees, all the enzymatic activities in which they were tested.

CONCLUSION

These results suggest that the venoms of the three Bitis species, involved in accidents with humans in the Sub-Saharan Africa, contain a mixture of various enzymes that may act in the generation and development of some of the clinical manifestations of the envenomations. We also demonstrated that horse antivenoms produced against B. arietans or B. g. rhinoceros plus B. nasicornis venoms can blocked some of the toxic activities of these venoms.

Identification and phylogeny of Arabian snakes: Comparison of venom chromatographic profiles versus 16S rRNA gene sequences

Identification of snake species is important for various reasons including the emergency treatment of snake bite victims. We present a simple method for identification of six snake species using the gel filtration chromatographic profiles of their venoms. The venoms of Echis coloratus, Echis pyramidum, Cerastes gasperettii, Bitis arietans, Naja arabica, and Walterinnesia aegyptia were milked, lyophilized, diluted and centrifuged to separate the mucus from the venom. The clear supernatants were filtered and chromatographed on fast protein liquid chromatography (FPLC). We obtained the 16S rRNA gene sequences of the above species and performed phylogenetic analysis using the neighbor-joining method. The chromatograms of venoms from different snake species showed peculiar patterns based on the number and location of peaks. The dendrograms generated from similarity matrix based on the presence/absence of particular chromatographic peaks clearly differentiated Elapids from Viperids. Molecular cladistics using 16S rRNA gene sequences resulted in jumping clades while separating the members of these two families. These findings suggest that chromatographic profiles of snake venoms may provide a simple and reproducible chemical fingerprinting method for quick identification of snake species. However, the validation of this methodology requires further studies on large number of specimens from within and across species.

Hybridization at an ecotone: ecological and genetic barriers between three Iberian vipers

The formation of stable genetic boundaries between emerging species is often diagnosed by reduced hybrid fitness relative to parental taxa. This reduced fitness can arise from endogenous and/or exogenous barriers to gene flow. Although detecting exogenous barriers in nature is difficult, we can estimate the role of ecological divergence in driving species boundaries by integrating molecular and ecological niche modelling tools. Here, we focus on a three-way secondary contact zone between three viper species (Vipera aspis, V. latastei and V. seoanei) to test for the contribution of ecological divergence to the development of reproductive barriers at several species traits (morphology, nuclear DNA and mitochondrial DNA). Both the nuclear and mitochondrial data show that all taxa are genetically distinct and that the sister species V. aspis and V. latastei hybridize frequently and backcross over several generations. We find that the three taxa have diverged ecologically and meet at a hybrid zone coincident with a steep ecotone between the Atlantic and Mediterranean biogeographical provinces. Integrating landscape and genetic approaches, we show that hybridization is spatially restricted to habitats that are suboptimal for parental taxa. Together, these results suggest that niche separation and adaptation to an ecological gradient confer an important barrier to gene flow among taxa that have not achieved complete reproductive isolation.

A phylogeny and revised classification of Squamata, including 4161 species of lizards and snakes

BACKGROUND

The extant squamates (>9400 known species of lizards and snakes) are one of the most diverse and conspicuous radiations of terrestrial vertebrates, but no studies have attempted to reconstruct a phylogeny for the group with large-scale taxon sampling. Such an estimate is invaluable for comparative evolutionary studies, and to address their classification. Here, we present the first large-scale phylogenetic estimate for Squamata.

RESULTS

The estimated phylogeny contains 4161 species, representing all currently recognized families and subfamilies. The analysis is based on up to 12896 base pairs of sequence data per species (average = 2497 bp) from 12 genes, including seven nuclear loci (BDNF, c-mos, NT3, PDC, R35, RAG-1, and RAG-2), and five mitochondrial genes (12S, 16S, cytochrome b, ND2, and ND4). The tree provides important confirmation for recent estimates of higher-level squamate phylogeny based on molecular data (but with more limited taxon sampling), estimates that are very different from previous morphology-based hypotheses. The tree also includes many relationships that differ from previous molecular estimates and many that differ from traditional taxonomy.

CONCLUSIONS

We present a new large-scale phylogeny of squamate reptiles that should be a valuable resource for future comparative studies. We also present a revised classification of squamates at the family and subfamily level to bring the taxonomy more in line with the new phylogenetic hypothesis. This classification includes new, resurrected, and modified subfamilies within gymnophthalmid and scincid lizards, and boid, colubrid, and lamprophiid snakes.

Snake velvet black: hierarchical micro- and nanostructure enhances dark colouration in Bitis rhinoceros

The West African Gaboon viper (Bitis rhinoceros) is a master of camouflage due to its colouration pattern. Its skin is geometrically patterned and features black spots that purport an exceptional spatial depth due to their velvety surface texture. Our study shades light on micromorphology, optical characteristics and principles behind such a velvet black appearance. We revealed a unique hierarchical pattern of leaf-like microstructures striated with nanoridges on the snake scales that coincides with the distribution of black colouration. Velvet black sites demonstrate four times lower reflectance and higher absorbance than other scales in the UV-near IR spectral range. The combination of surface structures impeding reflectance and absorbing dark pigments, deposited in the skin material, provides reflecting less than 11% of the light reflected by a polytetrafluoroethylene diffuse reflectance standard in any direction. A view-angle independent black structural colour in snakes is reported here for the first time.

The continuing saga of snake venom disintegrins

Disintegrins, a family of polypeptides released in the venoms of viperid snakes (vipers and rattlesnakes) by the proteolytic processing of multidomain metalloproteinases, selectively block the function of β(1) and β(3) integrin receptors. Few of the proteins isolated and characterized from snake venoms have proven to be more structural and functional versatile than the disintegrins. Not surprisingly, 25 years after their discovery, our knowledge on the evolutionary history and the molecular determinants modulating the integrin inhibitory activity of disintegrins still remain fragmentary. This paper highlights some seminal contributions, including personal accounts of pioneer authors, related to basic and applied research on disintegrins. Investigators have evaluated disintegrin applications in therapies for a number of pathologies in which integrin receptors play relevant roles, particularly myocardial infarction and inappropriate tumor angiogenesis. Completing the continuing story of the disintegrin family by applying novel research approaches may hold the key to learn how to use deadly toxins as therapeutic agents.

Recombinant expression of mutants of the Frankenstein disintegrin, RTS-ocellatusin. Evidence for the independent origin of RGD and KTS/RTS disintegrins

The requirements to transform a short disintegrin of the RGD clade into an RTS disintegrin, were investigated through the generation of recombinant mutants of ocellatusin in which the RGD tripeptide was substituted for RTS in different positions along the integrin-specificity loop. Any attempt to create an active integrin α(1)β(1) inhibitory motif within the specificity loop of ocellatusin was unsuccessful. Replacing the whole RGD-loop of ocellatusin by the RTS-loop of jerdostatin was neither sufficient for confering α(1)β(1) binding specificity to this ocellatusin-RTS Frankenstein(2) mutant. Factors other than the integrin-binding loop sequence per se are thus required to transform a disintegrin scaffold from the RGD clade into another scaffold from the RTS/KTS clade. Moreover, our results provide evidences, that the RTS/KTS short disintegrins have potentially been recruited into the venom gland of Eurasian vipers independently from the canonical neofunctionalization pathway of the RGD disintegrins. PCR-amplifications of jerdostatin-like sequences from a number of taxa across reptiles, including snakes (Crotalinae, Viperinae, and Elapidae taxa) and lizards (Lacertidae and Iguanidae) clearly showed that genes coding for RTS/KTS disintegrins existed long before the split of Lacertidae and Iguania, thus predating the recruitment of the SVMP precursors of disintegrins, providing strong support for the view of an independent evolutionary history of the RTS/KTS and the RGD clades of short disintegrins.

Epidemiology of snakebites in Europe: a systematic review of the literature

Snakebites are rare medical emergency cases in Europe but may sometimes be severe and lead to complications. A better knowledge of snakebite epidemiology may help health authorities to better understand therapeutic requirements, especially concerning antivenoms, and thus improve treatment of snakebite. An extensive literature search for studies and articles published between 1970 and 2010 was performed. Both indexed and non-indexed articles were examined, the analysis of which took into account the heterogeneity between the studies and weighted the studies according to size of the study population covered. Most of the articles involved hospitalized patients who represented more than 90% of snakebites. Incidence, mortality and population at risk were estimated after stratification into three regions (northern, central and southern Europe) based both on viper species distribution and climatic characteristics. There was no significant variation in incidence from the north to the south of Europe. In the whole of Europe, including European Russia and Turkey, the annual number of snakebite cases was estimated at 7992 [CI 95% = 6860-9178] bites, out of which approximately 15% were considered severe (grade 3). These bites usually occurred between May and September, with a more dispersed distribution in southern Europe. The average number of deaths per annum was 4 [0.7-7.7]. Children and male victims are more affected, contrary to what one would expect given their respective proportion in the entire population. Both upper and lower limb bites were recorded at an equal frequency while the bites in other parts of the body were very rare. Immunotherapy was prescribed in one out of three snakebites in Europe, with a very high geographical variability, in spite of excellent tolerance, at least considering highly-purified immunoglobulin fragments. Snakebites are uncommon in Europe but can cause life-threatening envenomation. Fragments of highly-purified immunoglobulins are now very well tolerated and dramatically reduce both severity and mortality of snakebites when used in treatment.

Neutralization of Vipera and Macrovipera venoms by two experimental polyvalent antisera: a study of paraspecificity

We conducted an extensive study of neutralization of lethality of 11 species and one subspecies of snakes of the genus Vipera, and of five species of Macrovipera, by two experimental equine antisera. One antiserum was a trivalent preparation raised against the venoms of Vipera aspis aspis, Vipera berus berus and Vipera ammodytes ammodytes; the other was a pentavalent preparation that also included venoms of Vipera (now Montivipera) xanthina and Macrovipera lebetina obtusa. We measured specific neutralization of lethality against all venoms included in the immunization schemes, and paraspecific neutralization against the venoms of Vipera ammodytes montandoni, Vipera (Montivipera) bornmuelleri, Vipera latastei, Vipera (Mo.) latifii, Vipera (Mo.) lotievi, Vipera (Daboia) palaestinae, Vipera (Mo.) raddei and Vipera seoanei, as well as against Macrovipera (D.) deserti, Macrovipera lebetina cernovi, Macrovipera lebetina turanica and Macrovipera schweitzeri. We found an important degree of paraspecific protection within each genera (omitting recent reclassification) that was quite independent of both the lethal potency of the venoms and their geographic origin. This information may be of use to clinicians charged with the treatment of Vipera or Macrovipera envenomations with non-specific antivenoms.