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.

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.

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.

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.

Discovery of a novel accessory structure of the pitviper infrared receptor organ (serpentes: viperidae)

The facial pits of rattlesnakes, copperheads, lanceheads, bushmasters and other American and Asian pitvipers (Crotalinae) are highly innervated and densely vascularized infrared (IR) receptor organs. For over a century, studies have focused on a small sample of model species from North America and Asia. Based on an expanded survey of Central and South American crotalines, we report a conspicuous accessory structure composed of well-defined papillae that project from the anterior orbital adnexa. The papillae are continuous with the inner chamber of the IR receptor organ and our histological and ultrastructural data suggest that they possess a well-developed nervous network and extensive vascularization; however, they lack the characteristic IR-sensitive terminal nerve masses found in the IR-receptive pit membrane. The function of the IR receptor organ papillae is unknown.

Clinico-epidemiological features of viper bite envenomation: a study from Manipal, South India

INTRODUCTION

Snakebite is an important and preventable health hazard. Viper bites are more common than other poisonous snakebites in human beings. The present study aimed to collate the victim profile of viper bite cases in the region and to determine the pattern, manifestations, complications and the associated risk factors of these bites.

METHODS

This was a prospective study of viper bite cases admitted to Kasturba Hospital, Manipal, India between August 2003 and November 2005. The demographic and clinical details of each case were obtained from the patients, their relatives and the patients' hospital records, and analysed.

RESULTS

A total of 31 viper bite cases were reported during the study period. The victims were predominantly male and aged 25-55 years. The highest number of cases occurred during daytime in the months of September and October, which coincided with the harvesting season, and involved the lower limbs. 94% of the snakebite victims were farmers, suggesting that this was an occupational hazard. Envenomation was observed in patients with scratch marks, suggesting the importance of keeping the victim under observation in all alleged snakebite cases, even in the absence of clear fang marks. The mortality rate in our study was 6.5%.

CONCLUSION

Immobilising and transporting snakebite victims to the hospital and prompt administration of anti-snake venom remain the best way to reduce morbidity and mortality. It is also important to practise correct first aid measures, as otherwise they may cause more harm than good.

Phospholipase A2 diversity and polymorphism in European viper venoms: Paradoxical molecular evolution in Viperinae

We report the diversity and polymorphism of phospholipase A(2) (PLA(2)) transcripts from snakes belonging to nine European viper subspecies. This diversity results in the expression of a combination of six PLA(2) species--ammodytin I1, ammodytin I2, ammodytin L, ammodytoxin, vaspin A and vaspin B--with 19 known isoforms of the first five of these species. Most of the European viper venoms studied contained either a myotoxin or a neurotoxin, and all contained ammodytin I1 and ammodytin I2. There is no evidence that a given pattern of PLA(2) species constitutes a taxonomic criterion, and isoform analysis would be required for such discrimination. Analysis of the phylogenetic relationships between PLA(2) species from European vipers and those of other members of the Viperinae revealed a strong correlation between the geographical source of the viper and the clustering seen for the different isoforms, for each PLA(2) species. The K(a)/K(s) values calculated for the mature protein-coding region of paralogous genes showed that ratios for pairs including vaspin B or one ammodytoxin isoform were greater than 1.09, whereas those for most of the remaining pairs were less than 1. Different patterns of mutation were observed in comparisons of the different PLA(2) isoforms. The mechanisms directing a mutation toward a precise exon remain unresolved.

Snake Venomics ofBitisSpecies Reveals Large Intragenus Venom Toxin Composition Variation: Application to Taxonomy of Congeneric Taxa

The protein composition of the venoms of the West African Gaboon viper (Bitis gabonica rhinoceros), the rhinoceros viper (Bitis nasicornis), and the horned puff adder (Bitis caudalis) were analyzed by RP-HPLC, N-terminal sequencing, SDS-PAGE, MALDI-TOF peptide mass fingerprinting, and CID-MS/MS. In line with previous proteomic and transcriptomic analyses showing that snake venom proteins belong to only a few major protein families, the venom proteomes of Bitis gabonica rhinoceros, Bitis nasicornis, and Bitis caudalis comprise, respectively, toxins from 11, 9, and 8 toxin families. Dimeric disintegrins, PLA2 molecules, serine proteinases, a CRISP, C-type lectin-like proteins, L-amino acid oxidases, and snake venom metalloproteases are present in the three Bitis snake venoms, though they depart from each other in the composition and the relative abundance of their toxins. The venom composition appears to keep information on the evolutionary history of congeneric taxa. Protein similarity coefficients used to estimate the similarity of venom proteins of the Bitis taxa sampled here and in previous studies (eg. Bitis arietans and Bitis gabonica gabonica) support the monophyly of the three West African taxa (B.g. gabonica, B.g. rhinoceros, and B. nasicornis) based on genetic distance reconstructions, the lack of alliances between B. arietans and any other Bitis species, and are consistent with the taxonomic association of Bitis caudalis within the differentiated group of small Bitis species. The low level of venom toxin composition similarity between the two conventionally recognized subspecies of Bitis gabonica, B. g. gabonica and B. g. rhinoceros, supports the consideration by some authors of B. g. rhinoceros as a separate species, Bitis rhinoceros. Moreover, our proteomic data fit better to a weighted phylogram based on overall genetic distances than to an unweighted maximum-parsimony tree.

Molecular phylogeography of endangered sharp-snouted pitviper (Deinagkistrodon acutus; Reptilia, Viperidae) in Mainland China

Using phylogenetic and population genetic approaches, the present study reports the phylogeographic structure of the sharp-snouted pitviper (Deinagkistrodon acutus), a threatened snake species with commercial and medicinal importance in China. The entire mitochondrial ND2 gene (NADH dehydrogenase subunit 2) sequences of 86 individuals of D. acutus from 14 localities across its range in China were determined. Based on the results of phylogenetic analyses, distribution of diagnostic sites, haplotype network, and AMOVA hierarchical analysis, an east-west division of the whole D. acutus population could be observed. Geographically, a line formed by a lake, river, and mountain chain (the Poyang Lake, Gan River to the southern end of the Wuyi Mountains), results in vicariance and approximately vertically splits the range into two and the whole population into two main lineages (western and eastern). The bifurcating tree suggested generally west to east dispersal trend. The data fit the isolation by distance (IBD) model well. Star-like clusters in haplotype network, significantly negative values of Fs statistics, and unimodal mismatch distributions all suggest recent demographic expansions in four areas. The results show that isolation, dispersal, bottleneck, and expansion jointly constitute the history of D. acutus. In a haplotype network, the excessive predominance of central haplotypes, few medium-frequency haplotypes, predominance (73.1%) of the singletons among the derived haplotypes, most of which are connected to the central haplotype by only one mutational step, unsymmetrical campanulate unimodal curve of mismatch distributions and leftwards shift of the peaks, all suggest that the whole D. acutus population is a young population with low genetic diversity. Based on the data, the first priority for conservation action should be given to the Huangshan unit.

Loss of introns along the evolutionary diversification pathway of snake venom disintegrins evidenced by sequence analysis of genomic DNA from Macrovipera lebetina transmediterranea and Echis ocellatus

Analysis of cDNAs from Macrovipera lebetina transmediterranea (Mlt) and Echis ocellatus (Eo) venom gland libraries encoding disintegrins argued strongly for a common ancestry of the messengers of short disintegrins and those for precursors of dimeric disintegrin chains. We now report the sequence analysis of disintegrin-coding genes from these two vipers. Genomic DNAs for dimeric disintegrin subunits Ml_G1 and Ml_G2 (Mlt) and Eo_D3 (Eo) contain single 1-kb introns exhibiting the 5'-GTAAG (donor)/3'-AG (acceptor) consensus intron splicing signature. On the other hand, the short RTS-disintegrins Ml_G3 (Mlt) and Eo_RTS (Eo) and the short RGD-disintegrin ocellatusin (Eo) are transcribed from intronless genomic DNA sequences, indicating that the evolutionary pathway leading to the emergence of short disintegrins involved the removal of all intronic sequences. The insertion position of the intron within Ml_G1, Ml_G2, and Eo_D3 is conserved in the genes for vertebrate ADAM (A disintegrin and metalloproteinase) protein disintegrin-like domains and within the gene for the medium-size snake disintegrins halystatins 2 and 3. However, a comparative analysis of currently available disintegrin(-like) genes outlines the view that a minimization of both the gene organization and the protein structure underlies the evolution of the snake venom disintegrin family.

Venomous Adversaries: A Reference to Snake Identification, Field Safety, and Bite-Victim First Aid for Disaster-Response Personnel Deploying Into the Hurricane-Prone Regions of North America

Each hurricane season, emergency-preparedness deployment teams including but not limited to the Office of Force Readiness and Deployment of the US Public Health Service, Federal Emergency Management Agency, Deployment Medical Assistance Teams, Veterinary Medical Assistance Teams, and the US Army and Air Force National Guard are at risk for deploying into hurricane-stricken areas that harbor indigenous hazards, including those posed by venomous snakes. North America is home to 2 distinct families of venomous snakes: 1) Viperidae, which includes the rattlesnakes, copperheads, and cottonmouths; and 2) Elapidae, in which the only native species are the coral snakes. Although some of these snakes are easily identified, some are not, and many rank among the most feared and misunderstood animals. This article specifically addresses all the native species of venomous snakes that inhabit the hurricane-prone regions of North America and is intended to serve as a reference to snake identification, basic field safety procedures, and the currently recommended first-aid measures for snakebite casualties.

Identifying the biting species in snakebite by clinical features: an epidemiological tool for community surveys

The outcome of snakebite is related to the biting species but it is often difficult to identify the biting snake, particularly in community settings. We have developed a clinical scoring system suitable for use in epidemiological surveys, with the main aim of identifying the presumed biting species in those with systemic envenoming who require treatment. The score took into account ten features relating to bites of the five medically important snakes in Sri Lanka, and an algorithm was developed applying different weightings for each feature for different species. A systematically developed artificial data set was used to fine tune the score and to develop criteria for definitive identification. The score was prospectively validated using 134 species-confirmed snakebites. It correctly differentiated the bites caused by the three snakes that commonly cause major clinical problems (Russell's viper (RV), kraits and cobra) from other snakes (hump-nosed viper (HNV) and saw-scaled viper (SSV)) with 80% sensitivity and 100% specificity. For individual species, sensitivity and specificity were, respectively: cobra 76%, 99%; kraits 85%, 99%; and RV 70%, 99%. As anticipated, the score was insensitive in the identification of bites due to HNV and SSV.

Evidence for a Müllerian mimetic radiation in Asian pitvipers

Müllerian mimicry, in which toxic species gain mutual protection from shared warning signals, is poorly understood in vertebrates, reflecting a paucity of examples. Indirect evidence for mimicry is found if monophyletic species or clades show parallel geographic variation in warning patterns. Here, we evaluate a hypothesis of Müllerian mimicry for the pitvipers in Southeast Asia using a phylogeny derived from DNA sequences from four combined mitochondrial regions. Mantel matrix correlation tests show that conspicuous red colour pattern elements are significantly associated with sympatric and parapatric populations in four genera. To our knowledge, this represents the first evidence of a Müllerian mimetic radiation in vipers. The putative mimetic patterns are rarely found in females. This appears paradoxical in light of the Müllerian prediction of monomorphism, but may be explained by divergent selection pressures on the sexes, which have different behaviours. We suggest that biased predation on active males causes selection for protective warning coloration, whereas crypsis is favoured in relatively sedentary females.

Some aspects of the venom proteome of the Colubridae snakePhilodryas olfersiirevealed from a Duvernoy's (venom) gland transcriptome

We investigated the putative toxins of Philodryas olfersii (Colubridae), a representative of a family of snakes neglected in venom studies despite their growing medical importance. Transcriptomic data of the venom gland complemented by proteomic analysis of the gland secretion revealed the presence of major toxin classes from the Viperidae family, including serine proteases, metalloproteases, C-type lectins, Crisps, and a C-type natriuretic peptide (CNP). Interestingly, the phylogenetic analysis of the CNP precursor showed it as a linker between two related precursors found in Viperidae and Elapidae snakes. We suggest that these precursors constitute a monophyletic group derived from the vertebrate CNPs.

Lachesis muta (Viperidae) cDNAs reveal diverging pit viper molecules and scaffolds typical of cobra (Elapidae) venoms: implications for snake toxin repertoire evolution

Efforts to describe toxins from the two major families of venomous snakes (Viperidae and Elapidae) usually reveal proteins belonging to few structural types, particular of each family. Here we carried on an effort to determine uncommon cDNAs that represent possible new toxins from Lachesis muta (Viperidae). In addition to nine classes of typical toxins, atypical molecules never observed in the hundreds of Viperidae snakes studied so far are highly expressed: a diverging C-type lectin that is related to Viperidae toxins but appears to be independently originated; an ohanin-like toxin, which would be the third member of the most recently described class of Elapidae toxins, related to human butyrophilin and B30.2 proteins; and a 3FTx-like toxin, a new member of the widely studied three-finger family of proteins, which includes major Elapidae neurotoxins and CD59 antigen. The presence of these common and uncommon molecules suggests that the repertoire of toxins could be more conserved between families than has been considered, and their features indicate a dynamic process of venom evolution through molecular mechanisms, such as multiple recruitments of important scaffolds and domain exchange between paralogs, always keeping a minimalist nature in most toxin structures in opposition to their nontoxin counterparts.

Bayesian mixed models and the phylogeny of pitvipers (Viperidae: Serpentes)

The subfamily Crotalinae (pitvipers) contains over 190 species of venomous snakes distributed in both the Old and New World. We incorporated an extensive sampling of taxa (including 28 of 29 genera), and sequences of four mitochondrial gene fragments (2.3kb) per individual, to estimate the phylogeny of pitvipers based on maximum parsimony and Bayesian phylogenetic methods. Our Bayesian analyses incorporated complex mixed models of nucleotide evolution that allocated independent models to various partitions of the dataset within combined analyses. We compared results of unpartitioned versus partitioned Bayesian analyses to investigate how much unpartitioned (versus partitioned) models were forced to compromise estimates of model parameters, and whether complex models substantially alter phylogenetic conclusions to the extent that they appear to extract more phylogenetic signal than simple models. Our results indicate that complex models do extract more phylogenetic signal from the data. We also address how differences in phylogenetic results (e.g., bipartition posterior probabilities) obtained from simple versus complex models may be interpreted in terms of relative credibility. Our estimates of pitviper phylogeny suggest that nearly all recently proposed generic reallocations appear valid, although certain Old and New World genera (Ovophis, Trimeresurus, and Bothrops) remain poly- or paraphyletic and require further taxonomic revision. While a majority of nodes were resolved, we could not confidently estimate the basal relationships among New World genera and which lineage of Old World species is most closely related to this New World group.

Optimal Intron Analyses in the Trimeresurus Radiation of Asian Pitvipers

Nuclear introns are commonly used as phylogenetic markers, but a number of issues related to alignment strategies, indel treatments, and the incorporation of length-variant heterozygotes (LVHs) are not routinely addressed when generating phylogenetic hypotheses. Topological congruence in relation to an extensive mitochondrial DNA multigene phylogeny (derived from 2,423 bp of 12S, 16S, ND4, and CYTB genes) of the Asian pitviper Trimeresurus radiation was used to compare combinations of "by eye" and edited and unedited ClustalX 1.8 alignments of two nuclear introns. Indels were treated as missing data, fifth character states, and assigned simple and multistate codes. Upon recovery of the optimal alignment and indel treatment strategy, a total evidence approach was used to investigate the phylogenetic utility of the indels and test new generic arrangements within Trimeresurus. Approximately one third of the intron data partitions exhibited LVHs, suggesting that they are common in introns. Furthermore, a simple concatenation approach can facilitate the incorporation of LVHs into phylogenetic analyses to make use of all available data and investigate mechanisms of molecular evolution. Analyses of ClustalX 1.8-assisted alignments were generally more congruent than the "by eye" alignment and the analysis of a simple coded, edited ClustalX 1.8 (gap opening cost 5, gap extension cost 1) alignment revealed the most congruent tree. The total evidence approach supported the new arrangements within Trimeresurus, suggesting that the phylogeny should be considered as a working benchmark in Asian pitviper systematics. Finally, a critical appraisal of the diverse array of indels (56 to 57 per intron, ranging from 1 to 151 bp in length) suggested that they are a combination of Hennigian and homoplasious events unrelated to indel size or location within the intron. [Alignment; indels; intron analysis; length-variant heterozygotes; Trimeresurus.].

Modeling nucleotide evolution at the mesoscale: the phylogeny of the neotropical pitvipers of the Porthidium group (viperidae: crotalinae)

We analyzed the phylogeny of the Neotropical pitvipers within the Porthidium group (including intra-specific through inter-generic relationships) using 1.4 kb of DNA sequences from two mitochondrial protein-coding genes (ND4 and cyt-b). We investigated how Bayesian Markov chain Monte-Carlo (MCMC) phylogenetic hypotheses based on this 'mesoscale' dataset were affected by analysis under various complex models of nucleotide evolution that partition models across the dataset. We develop an approach, employing three statistics (Akaike weights, Bayes factors, and relative Bayes factors), for examining the performance of complex models in order to identify the best-fit model for data analysis. Our results suggest that: (1) model choice may have important practical effects on phylogenetic conclusions even for mesoscale datasets, (2) the use of a complex partitioned model did not produce widespread increases or decreases in nodal posterior probability support, and (3) most differences in resolution resulting from model choice were concentrated at deeper nodes. Our phylogenetic estimates of relationships among members of the Porthidium group (genera: Atropoides, Cerrophidion, and Porthidium) resolve the monophyly of the three genera. Bayesian MCMC results suggest that Cerrophidion and Porthidium form a clade that is the sister taxon to Atropoides. In addition to resolving the intra-specific relationships among a majority of Porthidium group taxa, our results highlight phylogeographic patterns across Middle and South America and suggest that each of the three genera may harbor undescribed species diversity.

Adenosine in the venoms from viperinae snakes of the genus Bitis: Identification and quantitation using LC/MS and CE/MS

Snake venoms are rich sources of toxic proteins and small molecules. This study was directed at molecules of molecular mass below 1 kDa. Thirty different venoms, of either neurotoxic or haemorrhagic type, were fractionated using size-exclusion chromatography. Only venoms of the Puff adder (Bitis arietans), Gaboon viper (Bitis gabonica), and Rhinoceros viper (Bitis nasicornis) exhibited large absorbance peaks at lambda(280 nm) in the total volume range of the chromatographic column indicating the presence of abundant low molecular mass material. Analysis of fractions containing this material using both HPLC and capillary electrophoresis interfaced with electrospray ion-trap mass spectrometry unequivocally established that the bioactive nucleoside, adenosine, was the major component. The concentrations of adenosine found (Puff adder--97.7 x 10(-6) mol L(-1); Gaboon viper--28.0 x 10(-6) mol L(-1); and Rhinoceros viper-56.8 x 10(-6) mol L(-1)) were above those required to activate all known sub-types of adenosine receptors. Adenosine may thus act at the site of envenomation causing local vasodilatation and may play a role in the subsequent systemic hypotension observed.

Estimates of disease burden due to land-snake bite in Sri Lankan hospitals

Snake bite is a common cause of hospital admission in Sri Lanka. Despite this, there have been no countrywide studies or national estimates of disease burden due to snake bites in Sri Lankan hospitals. We assessed the disease burden due to snake bite in our hospitals and estimated the frequency of admissions due to bites by different snake species. Sri Lanka was divided into four zones based on climate and topography. Hospital morbidity and mortality data, which are available on an administrative district basis, were collated for the four zones. A survey of opinion among specialist physicians (the Delphi technique) was used to estimate the proportion of bites by different species, and requirements for anti-venom (AV) and intensive care facilities for management of snake bites in hospitals in each of the four zones. A study of hospital admissions due to snake bites in seven selected hospitals was also performed to validate the opinion survey. There was a clear difference in the incidence of hospital admissions due to snake bites in the different zones. Estimates of hospital admissions due to bites by different species also varied considerably between zones. These trends corresponded to estimates of requirements of AV and other supportive health care. Health care planning using data based on environmental information, rather than merely on political boundaries, could lead to targeted distribution of AV and intensive care requirements to manage snake bites.

Molecular phylogeny of Vipera Laurenti, 1768 and the related genera Macrovipera (Reuss, 1927) and Daboia (Gray, 1842), with comments about neurotoxic Vipera aspis aspis populations

We used mtDNA sequences (cytochrome b and NADH dehydrogenase subunit 2) to reconstruct molecular phylogenies of Vipera sensu lato, Vipera sensu stricto, and Vipera aspis. Three major clades were identified within the Vipera s.l. group: (1) the European vipers, (2) the oriental vipers, consisting of Montivipera (Vipera 2) plus Macrovipera lebetina, and (3) a group of Asian and North African vipers consisting of Daboia russelii, V. palaestinae, and Macrovipera mauritanica. We also distinguished three clades within the monophyletic European Vipera group: V. ammodytes, V. aspis, and V. latastei, and Pelias with monophyly of Vipera 1 uncertain. Within V. aspis, the specimens collected in France formed the sister group of an Italian clade. The "neurotoxic" French population of V. aspis, which has a specific venom profile, separated from other French V. aspis early in the history of this group.

Infrared imaging in vipers: differential responses of crotaline and viperine snakes to paired thermal targets

Pit vipers use infrared-sensitive pit organs to accurately target homeothermic prey even in the absence of visual cues. It has been suggested that other vipers, including large ambush predators of the genus Bitis, also may use radiant infrared information for predatory targeting. We compared behavioral responses of pit vipers and snakes of the viperine genus Bitis to paired targets of different temperatures (i.e. cool and warm balloons), some scented with rodent odors. The rates of tongue flicking, head turning and approaches by pit vipers (Agkistrodon contortrix and Crotalus atrox) were significantly higher toward warm targets than toward cool ones. Moreover, they all were significantly higher in pit vipers than in the vipers Bitis arietans, Bitis gabonica, and Bitis nasicornis. Bitis sp. exhibited no significant differences in their behaviors toward warm versus cool targets. Pit vipers often struck at targets (always the warm target even when paired with a prey-scented cool target), but Bitis sp. never struck at either warm or cool targets. These results show that the behavioral correlates of infrared-based predatory and/or defensive targeting differ significantly between these two groups of viperid snakes, and suggest that the neural substrates of infrared imaging in crotaline snakes are either absent or functionally distinct in viperine snakes of the genus Bitis.

Novel sequences encoding venom C-type lectins are conserved in phylogenetically and geographically distinct Echis and Bitis viper species

Envenoming by Echis saw scaled vipers and Bitis arietans puff adders is the leading cause of death and morbidity in Africa due to snake bite. Despite their medical importance, the composition and constituent functionality of venoms from these vipers remains poorly understood. Here, we report the cloning of cDNA sequences encoding seven clusters or isoforms of the haemostasis-disruptive C-type lectin (CTL) proteins from the venom glands of Echis ocellatus, E. pyramidum leakeyi, E. carinatus sochureki and B. arietans. All these CTL sequences encoded the cysteine scaffold that defines the carbohydrate-recognition domain of mammalian CTLs. All but one of the Echis and Bitis CTL sequences showed greater sequence similarity to the beta than alpha CTL subunits in venoms of related Asian and American vipers. Four of the new CTL clusters showed marked inter-cluster sequence conservation across all four viper species which were significantly different from that of previously published viper CTLs. The other three Echis and Bitis CTL clusters showed varying degrees of sequence similarity to published viper venom CTLs. Because viper venom CTLs exhibit a high degree of sequence similarity and yet exert profoundly different effects on the mammalian haemostatic system, no attempt was made to assign functionality to the new Echis and Bitis CTLs on the basis of sequence alone. The extraordinary level of inter-specific and inter-generic sequence conservation exhibited by the Echis and Bitis CTLs leads us to speculate that antibodies to representative molecules should neutralise the biological function of this important group of venom toxins in vipers that are distributed throughout Africa, the Middle East and the Indian subcontinent.

A new species of hognose pitviper, genus Porthidium, from the southwestern Pacific of Costa Rica (Serpentes: Viperidae)

A new species of terrestrial pitviper, Porthidium porrasi, is described from mesophytic forests of the Península de Osa and surrounding area of the Pacific versant of southwestern Costa Rica. It is most similar to P. nasutum and is characterized by a pattern of bands, persistence of the juvenile tail color in adults, and a high number of dorsal scales. Analysis of mtDNA sequences confirms its distinction from P. nasutum. The existence of this species reinforces the notion of elevated herpetofaunal endemism in southwestern Costa Rica.

Inhibition of hemorrhagic and edematogenic activities of snake venoms by a broad-spectrum protease inhibitor, murinoglobulin; the effect on venoms from five different genera in Viperidae family

In order to obtain basic data on the effect of broad-spectrum protease inhibitor against local symptoms of Viperidae snake envenomation, inhibitory capacity of rat murinoglobulin on local hemorrhagic and edematogenic activities of venoms from Crotalus atrox, Bothrops jararaca, Lachesis muta muta, Trimeresurus flavoviridis and Echis carinatus sochureki were examined. Murinoglobulin, pre-incubated with the crude venoms at 37 degrees C for 15 min, inhibited hemorrhagic activity of all five venoms to various extents. The activity of C. atrox was almost completely inhibited at the murinoglobulin/venom ratio (w/w) of 20. The activity of B. jararaca, Lachesis muta muta and T. flavoviridis venoms was considerably inhibited at the ratio of 20 (77.2, 80.0 and 86.2% inhibition, respectively), however some of the activity still remained even at the ratio of 40 (84.2, 79.8 and 86.2% inhibition, respectively). Among the five venoms, E. c. sochureki venom is quite resistant to murinoglobulin treatment and statistically significant inhibition was only found at the ratio of 40 (64.1% inhibition). Fibrinolytic and gelatinase activities were more susceptible to murinoglobulin inhibition. The treatment at the ratios of 10 and 20 almost completely inhibited respectively the fibrinolytic and the gelatinase activities of all the venoms. Murinoglobulin treatment also significantly inhibited the edematogenic activity of L. muta muta, T. flavoviridis and Echis carinatus sochureki. The treatment of murinoglobulin at the ratio of 40 considerably suppressed the swelling up to 60 min after subcutaneous injection of L. muta muta and E. c. sochureki venoms, and up to 30 min after T. flavoviridis venom injection. Murinoglobulin is a potent inhibitor against local effects of multiple snake venoms in Viperidae family.

Assessing the phylogenetic utility of four mitochondrial genes and a nuclear intron in the asian pit viper genus, Trimeresurus: separate, simultaneous, and conditional data combination analyses

A number of methods have been proposed for addressing how to optimize the analysis of multiple data sets from diverse mitochondrial and nuclear gene partitions in the pursuit of robust organismal phylogenies. The present study used separate, simultaneous, and conditional data combination methods to analyze 3,135 bp of data from four mitochondrial partitions and the seventh intron of the beta fibrinogen gene in the Asian pit viper genus, Trimeresurus sensu stricto. The phylogenetic utility and homogeneity of all partitions were estimated via a combination of homogeneity partition tests, homoplasy indices, and partitioned Bremer support. Despite the detection of significant heterogeneity of phylogenetic signal between the mitochondrial and nuclear partitions, the simultaneous analysis represented the best-supported topology of all the data. The relatively slow rate (approximately one quarter of the rate of mtDNA) and functionally unconstrained molecular evolution of the intron resulted in much lower levels of homoplasy compared with the mitochondrial partitions. This was further shown via partitioned Bremer support, which, when considered throughout hierarchical clade levels, highlighted the phylogenetic strength and limitations of the intron at deeper and shallower phylogenetic levels, respectively. The simultaneous analysis helped to resolve the phylogenetic relationships of taxa that were unresolved throughout all individual gene trees and tentatively supports the existence of morphologically and genetically distinct clades within the genus. Topological appraisals of the mitochondrial gene partitions suggest that the cytochrome b and the NADH subunit 4 gene partitions are better estimators of phylogenetic relationships than are the 12S and 16S ribosomal RNA partitions at the taxonomic levels under consideration.

Adaptive evolution in the snake venom Kunitz/BPTI protein family

Snake venoms are rich sources of serine proteinase inhibitors that are members of the Kunitz/BPTI (bovine pancreatic trypsin inhibitor) family. However, only a few of their gene sequences have been determined from snakes. We therefore cloned the cDNAs for the trypsin and chymotrypsin inhibitors from a Vipera ammodytes venom gland cDNA library. Phylogenetic analysis of these and other snake Kunitz/BPTI homologs shows the presence of three clusters, where sequences cluster by functional role. Analysis of the nucleotide sequences from the snake Kunitz/BPTI family shows that positive Darwinian selection was operating on the highly conserved BPTI fold, indicating that this family evolved by gene duplication and rapid diversification.

Sequences and structural organization of phospholipase A2 genes from Vipera aspis aspis, V. aspis zinnikeri and Vipera berus berus venom. Identification of the origin of a new viper population based on ammodytin I1 heterogeneity

We used a PCR-based method to determine the genomic DNA sequences encoding phospholipases A2 (PLA2s) from the venoms of Vipera aspis aspis (V. a. aspis), Vipera aspis zinnikeri (V. a. zinnikeri), Vipera berus berus (V. b. berus) and a neurotoxic V. a. aspis snake (neurotoxic V. a. aspis) from a population responsible for unusual neurotoxic envenomations in south-east France. We sequenced five groups of genes, each corresponding to a different PLA2. The genes encoding the A and B chains of vaspin from the neurotoxic V. a. aspis, PLA2-I from V. a. zinnikeri, and the anticoagulant PLA2 from V. b. berus are described here. Single nucleotide differences leading to amino-acid substitutions were observed both between genes encoding the same PLA2 and between genes encoding different PLA2s. These differences were clustered in exons 3 and 5, potentially altering the biological activities of PLA2. The distribution and characteristics of the PLA2 genes differed according to the species or subspecies. We characterized for the first time genes encoding neurotoxins from the V. a. aspis and V. b. berus snakes of central France. Genes encoding ammodytins I1 and I2, described previously in Vipera ammodytes ammodytes (V. am. ammodytes), were also present in V. a. aspis and V. b. berus. Three different ammodytin I1 gene sequences were characterized: one from V. b. berus, the second from V. a. aspis, V. a. zinnikeri and the neurotoxic V. a. aspis, and the third from the neurotoxic V. a. aspis. This third sequence was identical with the reported sequence of the V. am. ammodytes ammodytin I1 gene. Genes encoding monomeric neurotoxins of V. am. ammodytes venom, ammodytoxins A, B and C, and the Bov-B LINE retroposon, a phylogenetic marker found in V. am. ammodytes genome, were identified in the genome of the neurotoxic V. a. aspis. These results suggest that the population of neurotoxic V. a. aspis snakes from south-east France may have resulted from interbreeding between V. a. aspis and V. am. ammodytes.

[The snakes of Senegal: an annotated species list]

Between 1990 and 2001, the laboratory of Paludologie of IRD at Dakar collected over 5,500 snakes from all over Senegal. By studying this collection, an entirely new species was discovered for science and eight new species for Senegal. The presence of many specimens of rare species and the great number of localities from which snakes were collected also allowed us to solve several delicate taxonomic problems and better to specify the biogeography of Senegalese snakes. Currently, the presence in Senegal of at least 55 different species has been established with certainty. The data collected simultaneously on the burden of snakebites for public health shows considerable differences according to geographic area, with populations from south-eastern Senegal being more exposed to the risk of death by snakebite.

Venomous snake husbandry in Thailand

A captive breeding program for venomous Thai snakes was established at the Queen Saovabha Memorial Institute at Bangkok, Thailand. This was necessary to secure a stable, healthy, and species-confirmed source of snake venom for antivenom production. In 1994, wild-caught specimens were collected, sexed, quarantined, and housed appropriately. All data in this report, with the exclusion of Table 6, were collected from 1994 to 1997. Two species were bred successfully in captivity to date during this study period. Although captive breeding has not yet been achieved with all species and subspecies, our early success was encouraging.

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

Nucleotide sequences of mitochondrial cytochrome b and 16S rRNA genes, totaling 946 bp, were used to reconstruct a molecular phylogeny of 42 species of the subfamily Viperinae representing 12 of the 13 recognized genera. Maximum-parsimony and maximum-likelihood were used as methods for phylogeny reconstruction with and without a posteriori weighting. When representatives of the Causinae were taken as outgroup, five major monophyletic groups were consistently identified: Bitis, Cerastes, Echis, the Atherini (Atheris s.l.), and the Eurasian viperines. Proatheris was affiliated with Atheris, and Adenorhinos clustered within Atheris. The African Bitis consisted of at least three monophyletic groups: (i) the B. gabonica group, (ii) the B. caudalis group, and (iii) the B. cornuta group. B. worthingtoni and B. arietans are not included in any of these lineages. Eurasian viperines could be unambiguously devided into four monophyletic groups: (i) Pseudocerastes and Eristicophis, (ii) European vipers (Vipera s.str.), (iii) Middle East Macrovipera plus Montivipera (Vipera xanthina group), and (iv) North African Macrovipera plus Vipera palaestinae and Daboia russelii. These evolutionary lineages are consistent with historical biogeographical patterns. According to our analyses, the viperines originated in the Oligocene in Africa and successively underwent a first radiation leading to the five basal groups. The radiation might have been driven by the possession of an effective venom apparatus and a foraging startegy (sit-wait-strike) superior in most African biomes and might have been adaptive. The next diversifications led to the Proatheris-Atheris furcation, the basal Bitis splitting, and the emergence of the basal lineages within the Eurasian stock. Thereafter, lineages within Echis, Atheris, and Cerastes evolved. The emergence of three groups within Vipera s.l. might have been forced by the existence of three land masses during the early Miocene in the area of the Paratethys and the Mediterranean Seas. Taxonomic consequences of these findings are discussed.

Comparative biochemistry of disintegrins isolated from snake venom: consideration of the taxonomy and geographical distribution of snakes in the genus Echis

Species in the genus Echis have been classified mainly based on their morphological appearance and the analytical patterns of their serum. However, re-classification of the genus Echis has recently been suggested by taxonomists, toxicologists, and clinicians, since there have been problems with the current classification, such as the efficacy of antivenoms used for treating bites and the broad geographical distribution of Echis snakes. In this study, we purified five novel disintegrins, the platelet aggregation inhibitors pyramidin A and B from the venom of Echis pyramidum, ocellatin from the venom of Echis ocellatus, and leucogastin A and B from the venom of Echis leucogaster, to compare their sequences and allow us to re-evaluate the classification of various species in the genus Echis. Comparison of the amino acid sequences of five new and four known isolated disintegrins from snake venoms of six Echis species and their distribution strongly support the recent re-classification of the genus Echis.

Generalized frequency coding: a method of preparing polymorphic multistate characters for phylogenetic analysis

A new method of coding polymorphic multiistate characters for phylogenetic analysis is presented. By dividing such characters into subcharacters, their frequency distributions can be represented with discrete states. Differential weighting is used to counter the effect of representing one character with multiple characters. The new method, generalized frequency coding (GFC), is potentially superior to previously used methods in that it incorporates more information and is applicable to both qualitative and quantitative characters. When applied to a previously published data set that includes both types of polymorphic multistate characters, the method performed well, as assessed with g1 and nonparametric bootstrap statistics and giving results congruent with those of other studies. The data set was also used to compare GFC with both gap-weighting and Manhattan distance step matrix coding. On these grounds and for philosophical reasons, we consider GFC to be a better estimator of phylogeny.

A new species of arboreal pitviper from the Atlantic versant of northern Central America

A new species of green, prehensile-tailed pitviper of the genus Bothriechis is described from the Atlantic slopes of eastern Guatemala and western Honduras. This species appears to be most closely related to B. bicolor of the Pacific versant of Chiapas (Mexico) and Guatemala. Several other species of Bothriechis occur on the Atlantic versant of northern Central America, including two montane species, B. aurifer and B. marchi but, with one possible exception, these are not known to be sympatric with the new species and occur in different mountain ranges. The widespread B. schlegelii occurs up to at least 900 m on the Sierra de Caral, where the lowest elevation recorded for the new species is 885 m.

Weighting and congruence: a case study based on three mitochondrial genes in pitvipers

Weighting strategies in a total evidence approach are often conducted in order to remove homoplasy, with the implicit hope to increase congruence between data partitions. Incongruence was measured using the test of Farris et al. (J. S. Farris, M. Källersjö, A. G. Kluge, and C. Bult, 1995, Cladistics 10, 315-319) within and between three mitochondrial genes from pitvipers (Crotalinae) in partitioning each codon position for the coding genes. Incongruence between codon positions of a gene can be worse than incongruence between genes at analogous positions. Third positions of codons behave very differently in terms of incongruence from one gene to another while showing similar patterns in saturation tests. Instead of removing characters in order to discard homoplasy, which is hopeless and does not increase general congruence, we advocate for the removal of those substitutions that are incongruent with the rest. The genus Calloselasma and its sister group the genus Hypnale are the most basal Crotalinae. Asiatic pitvipers are paraphyletic, while American pitvipers are monophyletic.

Venomous snakebite in Thailand. I: Medically important snakes

Thailand has an abundance of venomous snakes. Among the neurotoxic family Elapidae, there are three species of the genus Naja (cobras), three of the genus Bungarus (kraits), and the king cobra of the genus Ophiophagus. Other Elapidae snakes in Thailand include sea snakes and Asian coral snakes of the genus Calliophis. They have potent venoms but rarely bite humans. Tissue and hemotoxic snakes are represented by family Viperidae, subfamilies Viperinae and Crotalinae. They remain an occupational hazard for farmers and rubber tappers, causing serious morbidity but only rare deaths, since competent treatment is now widely available throughout Thailand. Purified equine antivenin is manufactured locally for the monocled and Siamese spitting cobras (Naja kaouthia and N. siamensis), king cobra (Ophiophagus hannah), banded krait (Bungarus fasciatus), most green pit vipers (Trimeresurus sp.), Malayan pit viper (Calloselasma rhodostoma), and the Siamese Russell's viper (Daboia russelli siamensis).

[Poisonous snakes of Europe]

Characteristics of poisonous snakes and their toxines are described. The appearance and biology of all European poisonous snakes, eight vipers (family Viperidae) and one opisthoglyph colubride snake (family Colubridae) are given in detail.

[Thrombin-like serine proteases in Cerasted venoms (Cerasted cerastes and Cerastes vipera)]

Cerastes cerastes and Cerastes vipera snake venoms are rich of thrombin-like, serine-protease polypeptides. Many proteins have been isolated, purified and characterized from these vipers. These proteins act in various way on blood coagulation pathway and platelet function.

Venomous bites and stings in the tropical world

Snakes of the families Viperidae and Elapidae are responsible for the high incidence of morbidity and mortality after snake bites in countries of West Africa, the Indian subcontinent, South-East Asia, New Guinea and Latin America. Envenoming can cause local effects, notably tissue necrosis; and systemic effects, including paralysis, haemostatic disturbances, shock, increased capillary permeability, myocardial damage, rhabdomyolysis and acute renal failure. Specific hyperimmune serum (antivenom) is the mainstay of medical treatment for severe envenoming. Ancillary treatments such as assisted ventilation, repletion of circulating volume, renal dialysis and surgical debridement of necrotic tissues are needed in some cases. Scorpion stings are a common medical problem in middle and southern America, North Africa and the Middle East. Vasodilator drugs are important to counter the effects of massive catecholamine release. Bites by spiders and stings by hymenoptera and marine animals are responsible for deaths and morbidity in some tropical countries.