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

Ontogenetic Variation in Biological Activities of Venoms from Hybrids between Bothrops erythromelas and Bothrops neuwiedi Snakes

Lance-headed snakes are found in Central and South America, and they account for most snakebites in Brazil. The phylogeny of South American pitvipers has been reviewed, and the presence of natural and non-natural hybrids between different species of Bothrops snakes demonstrates that reproductive isolation of several species is still incomplete. The present study aimed to analyze the biological features, particularly the thrombin-like activity, of venoms from hybrids born in captivity, from the mating of a female Bothrops erythromelas and a male Bothrops neuwiedi, two species whose venoms are known to display ontogenetic variation. Proteolytic activity on azocoll and amidolytic activity on N-benzoyl-DL-arginine-p-nitroanilide hydrochloride (BAPNA) were lowest when hybrids were 3 months old, and increased over body growth, reaching values similar to those of the father when hybrids were 12 months old. The clotting activity on plasma diminished as hybrids grew; venoms from 3- and 6-months old hybrids showed low clotting activity on fibrinogen (i.e., thrombin-like activity), like the mother venom, and such activity was detected only when hybrids were older than 1 year of age. Altogether, these results point out that venom features in hybrid snakes are genetically controlled during the ontogenetic development. Despite the presence of the thrombin-like enzyme gene(s) in hybrid snakes, they are silenced during the first six months of life.

Mitochondrial genome of Ogmocotyle sikae and implications for phylogenetic studies of the Notocotylidae trematodes

Ogmocotyle spp. (Trematoda: Digenea: Notocotylidae) are neglected but important trematodes that can infect numerous mammal species, causing significant economic losses to livestock industries. However, there have been few studies on the molecular ecology of these trematodes. We amplified and sequenced the complete mitochondrial (mt) genome of Ogmocotyle sikae (14,307 bp). The gene content and gene arrangement of O. sikae mt genome was similar to that of Dicrocoelium chinensis, except that trnE and trnG were reversed. Phylogenetic analysis of O. sikae and selected parasites using Bayesian inference was performed based on concatenated amino acid sequence datasets conceptually translated from the 12 protein-coding genes. The results indicated that the family Notocotylidae is related to the family Paramphistomatidae. Our description of O. sikae mt genome provides a significant resource of molecular markers for future comparative studies of the Notocotylidae and other trematodes.

Eye and pit size are inversely correlated in crotalinae: Implications for selection pressure relaxation

Mate, prey, and predator recognition often depend on the integration of information from multiple sensory modalities including visual, auditory, and/or olfactory inputs. In Crotalinae, the eyes sense visible light while the pit organs detect infrared (IR) radiation. Previous studies indicate that there is significant overlap between the eye and pit sensory fields and that both senses are involved in recognition processes. This study investigated the relationships between eye and pit sizes in this taxonomic group as a function of phylogeny and habitat. In view of the fact that pit orientation depends largely on snout shape, pit vipers were grouped as follows: 1) arboreal, 2) terrestrial with rounded snout, and 3) terrestrial with pointed snout. The pit orientations and habitant patterns were fully independent of the Crotalinae phylogenetic tree. The phylogenetic generalized least squares model showed that both eye and pit areas were not of significantly phylogenetic relatedness, implying alternatively a strong effect of adaptation on eye and pit sizes. Negative correlations between relative eye and pit areas in terrestrial (both pointed and rounded snouts) and arboreal species were statistically significant. Our results suggest that the eyes and pits function in a complementary fashion such that selection for IR-perception relaxes selection pressures on the visual system and selection for visual discrimination relaxes selection pressures acting on the IR-system.

Combined venomics, venom gland transcriptomics, bioactivities, and antivenomics of two Bothrops jararaca populations from geographic isolated regions within the Brazilian Atlantic rainforest

Bothrops jararaca is a slender and semi-arboreal medically relevant pit viper species endemic to tropical and subtropical forests in southern Brazil, Paraguay, and northern Argentina (Misiones). Within its geographic range, it is often abundant and is an important cause of snakebite. Although no subspecies are currently recognized, geographic analyses have revealed the existence of two well-supported B. jararaca clades that diverged during the Pliocene ~3.8Mya and currently display a southeastern (SE) and a southern (S) Atlantic rainforest (Mata Atlântica) distribution. The spectrum, geographic variability, and ontogenetic changes of the venom proteomes of snakes from these two B. jararaca phylogroups were investigated applying a combined venom gland transcriptomic and venomic analysis. Comparisons of the venom proteomes and transcriptomes of B. jararaca from the SE and S geographic regions revealed notable interpopulational variability that may be due to the different levels of population-specific transcriptional regulation, including, in the case of the southern population, a marked ontogenetic venom compositional change involving the upregulation of the myotoxic PLA2 homolog, bothropstoxin-I. This population-specific marker can be used to estimate the proportion of venom from the southern population present in the B. jararaca venom pool used for the Brazilian soro antibotrópico (SAB) antivenom production. On the other hand, the southeastern population-specific D49-PLA2 molecules, BinTX-I and BinTX-II, lend support to the notion that the mainland ancestor of Bothrops insularis was originated within the same population that gave rise to the current SE B. jararaca phylogroup, and that this insular species endemic to Queimada Grande Island (Brazil) expresses a pedomorphic venom phenotype. Mirroring their compositional divergence, the two geographic B. jararaca venom pools showed distinct bioactivity profiles. However, the SAB antivenom manufactured in Vital Brazil Institute neutralized the lethal effect of both venoms to a similar extent. In addition, immobilized SAB antivenom immunocaptured most of the venom components of the venoms of both B. jararaca populations, but did not show immunoreactivity against vasoactive peptides. The Costa Rican bothropic-crotalic-lachesic (BCL) antivenom showed the same lack of reactivity against vasoactive peptides but, in addition, was less efficient immunocapturing PI- and PIII-SVMPs from the SE venom, and bothropstoxin-I, a CRISP molecule, and a D49-PLA2 from the venom of the southern B. jararaca phylogroup. The remarkable paraspecificity exhibited by the Brazilian and the Costa Rican antivenoms indicates large immunoreactive epitope conservation across the natural history of Bothrops, a genus that has its roots in the middle Miocene. This article is part of a Special Issue entitled: Omics Evolutionary Ecolog.

Venomics and antivenomics of Bothrops erythromelas from five geographic populations within the Caatinga ecoregion of northeastern Brazil

The Caatinga lancehead, Bothrops erythromelas, is a medically relevant species, responsible for most of the snakebite accidents in most parts of its distribution range in northeastern Brazil. The spectrum and geographic variability of its venom toxins were investigated applying a venomics approach to venom pools from five geographic areas within the Caatinga ecoregion. Despite its wide habitat, populations of B. erythromelas from Ceará, Pernambuco, Juazeiro, Paraiba, and Ilha de Itaparica exhibit highly conserved venom proteomes. Mirroring their compositional conservation, the five geographic venom pools also showed qualitatively and quantitatively overlapping antivenomic profiles against antivenoms generated in Vital Brazil (BR) and Clodomiro Picado (CR) Institutes, using different venoms in the immunization mixtures. The paraspecificity exhibited by the Brazilian SAB and the Costa Rican BCL antivenoms against venom toxins from B. erythromelas indicates large immunoreactive epitope conservation across genus Bothrops during the last ~14 million years, thus offering promise for the possibility of generating a broad-spectrum bothropic antivenom. Biological Significance Accidental snakebite envenomings represent an important public health hazard in Brazil. Ninety per cent of the yearly estimated 20-30,000 snakebite accidents are caused by species of the Bothrops genus. Bothrops erythromelas, a small, moderately stocky terrestrial venomous snake, is responsible for most of the snakebite accidents in its broad distribution range in the Caatinga, a large ecoregion in northeastern Brazil. To gain a deeper insight into the spectrum of medically important toxins present in the venom of the Caatinga lancehead, we applied a venomics approach to define the proteome and geographic variability of adult B. erythromelas venoms from five geographic regions. Although intraspecific compositional variation between venoms among specimens from different geographic regions has long been appreciated by herpetologists and toxinologists as a general feature of highly adaptable and widely distributed snake species, the five B. erythromelas populations investigated exhibit highly conserved venom proteomes. The overall toxin profile of the Caatinga lancehead's venom explains the local and systemic effects observed in envenomations by B. erythromelas. The five geographic venom pools sampled also showed qualitatively and quantitatively overlapping antivenomic profiles against antivenoms generated using different bothropic venoms in the immunization mixtures. The large immunoreactive epitope conservation across genus Bothrops offers promise for the generation of a broad-spectrum bothropic antivenom.

Phylogeny of courtship and male-male combat behavior in snakes

Background

Behaviors involved in courtship and male-male combat have been recorded in a taxonomically broad sample (76 species in five families) of snakes in the clade Boidae + Colubroidea, but before now no one has attempted to find phylogenetic patterns in such behaviors. Here, we present a study of phylogenetic patterns in such behaviors in snakes.

Methodology/Principal Findings

From the literature on courtship and male-male combat in snakes we chose 33 behaviors to analyze. We plotted the 33 behaviors onto a phylogenetic tree to determine whether phylogenetic patterns were discernible. We found that phylogenetic patterns are discernible for some behaviors but not for others. For behaviors with discernible phylogenetic patterns, we used the fossil record to determine minimum ages for the addition of each behavior to the courtship and combat behavioral repertoire of each snake clade.

Conclusions/Significance

The phylogenetic patterns of behavior reveal that male-male combat in the Late Cretaceous common ancestors of Boidae and Colubridae involved combatants raising the head and neck and attempting to topple each other. Poking with spurs was added in Boidae. In Lampropeltini the toppling behavior was replaced by coiling without neck-raising, and body-bridging was added. Phylogenetic patterns reveal that courtship ancestrally involved rubbing with spurs in Boidae. In Colubroidea, courtship ancestrally involved chin-rubbing and head- or body-jerking. Various colubroid clades subsequently added other behaviors, e.g. moving undulations in Natricinae and Lampropeltini, coital neck biting in the Eurasian ratsnake clade, and tail quivering in Pantherophis. The appearance of each group in the fossil record provides a minimum age of the addition of each behavior to combat and courtship repertoires. Although many gaps in the story of the evolution of courtship and combat in snakes remain, this study is an important first step in the reconstruction of the evolution of these behaviors in snakes.

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.

Phylogenetic relationships of yellowjackets inferred from nine loci (Hymenoptera: Vespidae, Vespinae, Vespula and Dolichovespula)

Eusociality has arisen repeatedly and independently in the history of insects, often leading to evolutionary success and ecological dominance. Eusocial wasps of the genera Vespula and Dolichovespula, or yellowjackets, have developed advanced social traits in a relatively small number of species. The origin of traits such as effective paternity and colony size has been interpreted with reference to an established phylogenetic hypothesis that is based on phenotypic data, while the application of molecular evidence to phylogenetic analysis within yellowjackets has been limited. Here, we investigate the evolutionary history of yellowjackets on the basis of mitochondrial and nuclear markers (nuclear: 28S, EF1α, Pol II, and wg; mitochondrial: 12S, 16S, COI, COII, and Cytb). We use these data to test the monophyly of yellowjackets and species groups, and resolve species-level relationships within each genus using parsimony and Bayesian inference. Our results indicate that a yellowjacket clade is either weakly supported (parsimony) or rejected (Bayesian inference). However, the monophyly of each yellowjacket genus as well as species groups are strongly supported and concordant between methods. Our results agree with previous studies regarding the monophyly of the Vespula vulgaris group and the sister relationship between the V. rufa and V. squamosa groups. This suggests convergence of large colony size and high effective paternity in the vulgaris group and V. squamosa, or a single origin of both traits in the most recent common ancestor of all Vespula species and their evolutionary reversal in the rufa group.

Venomics of New World pit vipers: genus-wide comparisons of venom proteomes across Agkistrodon

We report a genus-wide comparison of venom proteome variation across New World pit vipers in the genus Agkistrodon. Despite the wide variety of habitats occupied by this genus and that all its taxa feed on diverse species of vertebrates and invertebrate prey, the venom proteomes of copperheads, cottonmouths, and cantils are remarkably similar, both in the type and relative abundance of their different toxin families. The venoms from all the eleven species and subspecies sampled showed relatively similar proteolytic and PLA2 activities. In contrast, quantitative differences were observed in hemorrhagic and myotoxic activities in mice. The highest myotoxic activity was observed with the venoms of A. b. bilineatus, followed by A. p. piscivorus, whereas the venoms of A. c. contortrix and A. p. leucostoma induced the lowest myotoxic activity. The venoms of Agkistrodon bilineatus subspecies showed the highest hemorrhagic activity and A. c. contortrix the lowest. Compositional and toxicological analyses agree with clinical observations of envenomations by Agkistrodon in the USA and Central America. A comparative analysis of Agkistrodon shows that venom divergence tracks phylogeny of this genus to a greater extent than in Sistrurus rattlesnakes, suggesting that the distinct natural histories of Agkistrodon and Sistrurus clades may have played a key role in molding the patterns of evolution of their venom protein genes.

BIOLOGICAL SIGNIFICANCE

A deep understanding of the structural and functional profiles of venoms and of the principles governing the evolution of venomous systems is a goal of venomics. Isolated proteomics analyses have been conducted on venoms from many species of vipers and pit vipers. However, making sense of these large inventories of data requires the integration of this information across multiple species to identify evolutionary and ecological trends. Our genus-wide venomics study provides a comprehensive overview of the toxic arsenal across Agkistrodon and a ground for understanding the natural histories of, and clinical observations of envenomations by, species of this genus.

Comparison of phylogeny, venom composition and neutralization by antivenom in diverse species of bothrops complex

In Latin America, Bothrops snakes account for most snake bites in humans, and the recommended treatment is administration of multispecific Bothrops antivenom (SAB – soro antibotrópico). However, Bothrops snakes are very diverse with regard to their venom composition, which raises the issue of which venoms should be used as immunizing antigens for the production of pan-specific Bothrops antivenoms. In this study, we simultaneously compared the composition and reactivity with SAB of venoms collected from six species of snakes, distributed in pairs from three distinct phylogenetic clades: Bothrops, Bothropoides and Rhinocerophis. We also evaluated the neutralization of Bothrops atrox venom, which is the species responsible for most snake bites in the Amazon region, but not included in the immunization antigen mixture used to produce SAB. Using mass spectrometric and chromatographic approaches, we observed a lack of similarity in protein composition between the venoms from closely related snakes and a high similarity between the venoms of phylogenetically more distant snakes, suggesting little connection between taxonomic position and venom composition. P-III snake venom metalloproteinases (SVMPs) are the most antigenic toxins in the venoms of snakes from the Bothrops complex, whereas class P-I SVMPs, snake venom serine proteinases and phospholipases A₂ reacted with antibodies in lower levels. Low molecular size toxins, such as disintegrins and bradykinin-potentiating peptides, were poorly antigenic. Toxins from the same protein family showed antigenic cross-reactivity among venoms from different species; SAB was efficient in neutralizing the B. atrox venom major toxins. Thus, we suggest that it is possible to obtain pan-specific effective antivenoms for Bothrops envenomations through immunization with venoms from only a few species of snakes, if these venoms contain protein classes that are representative of all species to which the antivenom is targeted.

Snakebite envenomation is a serious health issue in Latin America, particularly in the Amazon, where antivenom administration may be delayed due to logistic constraints. Bothrops snakes are involved in most of the snakebite-related accidents in Brazil. This work reports a comparative study of the toxin composition and antigenicity of the Bothrops venoms used to prepare the commercial antivenom and its effectiveness against the venom from Bothrops atrox, a prevalent Amazon species that is not included in the pool. Our data show a lack of connection between Bothrops taxonomic identity and venom composition. We also show that different toxins display distinct reactivity with the tested antivenom. However, the antivenom reacted similarly with each class of toxin present in the venoms of the different snakes studied. Important evidence was the neutralization of the major toxic effects of B. atrox venom, not included in the mixture of antigens used to produce the antivenom. Based on the observed antigenicity of the distinct protein classes of toxins, we suggest that it is possible to obtain pan-specific and efficient Bothrops antivenoms via immunization with venoms from a few species of snakes that are representative of the protein composition of a large number of targeted species.

Phylogenetic relationships of the enigmatic longtailed rattlesnakes (Crotalus ericsmithi, C. lannomi, and C. stejnegeri)

The longtailed rattlesnakes of western Mexico represent an enigmatic group of poorly known venomous snake species: Crotalus ericsmithi, C. lannomi, and C. stejnegeri. In the 120 years since their discovery, fewer than twenty individuals have been deposited in natural history collections worldwide. These three species share similar morphological traits, including a particularly long tail that has been interpreted as either an ancestral condition among rattlesnakes or as derived within the longtailed group. An understanding of the phylogenetic distinctiveness and relationships among the longtailed rattlesnakes, and their relationships to other rattlesnake groups, has previously been hampered by a dearth of comparative material and tissues for collection of DNA sequence data. Facilitated by the recent availability of tissue samples from multiple individuals of each species, we estimate the phylogenetic relationships among the longtailed rattlesnakes and their placement among other rattlesnake groups, using DNA sequence data from three mitochondrial and three nuclear gene fragments. We explore phylogenetic signal in our data using Bayesian and maximum likelihood methods, species tree analyses and hypothesis testing. Our results strongly support the monophyly of longtailed rattlesnakes and suggest the three species diverged from each other during the mid to late Pliocene or early Pleistocene (~1.5-5.6 mya). Contrary to prevailing hypotheses, we find no evidence for an early or basal divergence of the longtailed clade within the rattlesnake tree, and instead estimate that it diverged relatively recently (~6.8 mya) from its sister lineage, composed of the diamondback rattlesnakes (C. atrox group) and the prairie rattlesnakes (C. viridis group). With our added sampling of lineages and identification of previously used problematic sequences, we provide a revised hypothesis for relationships among Crotalus species, yet underscore the need for future studies and new data to deliver a well-supported robust estimate of rattlesnake relationships.

A relict lineage and new species of green palm-pitviper (Squamata, Viperidae, Bothriechis) from the Chortís Highlands of Mesoamerica

A new species of palm-pitviper of the genus Bothriechis is described from Refugio de Vida Silvestre Texíguat in northern Honduras. The new species differs from congeners by having 19 dorsal scale rows at midbody, a bright green dorsal coloration in adults, the prelacunal scale fused to the second supralabial, and in representing a northern lineage that is sister to Bothriechis lateralis, which is distributed in Costa Rica and western Panama and is isolated from the new taxon by the Nicaraguan Depression. This represents the 15th endemic species occurring in Refugio de Vida Silvestre Texíguat, one of the richest herpetofaunal sites in Honduras, itself being the country with the highest degree of herpetofaunal endemism in Central America. We name this new species in honor of a Honduran conservationist slain in fighting against illegal logging, highlighting the sacrifices of rural activists in battling these issues and the critical importance of conservation in these areas.

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.

Using rattlesnake microsatellites to determine paternity in captive bushmasters (Lachesis muta)

We used 10 microsatellite DNA markers originally described for two Crotalus and one Sistrurus species to infer paternity in a captive-hatched clutch of Lachesis muta. Although the dam was known, records listed two potential sires, which prevented the inclusion of those offspring in a captive breeding program. Samples were collected from both possible sires, the presumed dam and two offspring. Five of the ten markers were paternity informative and either unambiguously identified one of the two males as the sire of the offspring or excluded the other male. A sixth marker identified the sire in one of the two offspring. These results reveal rattlesnake DNA markers as a useful tool in paternity testing and captive breeding management of Lachesis muta. It also indicates that, within the American crotalid radiation, a selection of microsatellite DNA markers cross-amplify and provide useful genotypic information across species and genera.

Molecular systematics and historical biogeography of tree boas (Corallus spp.)

Inferring the evolutionary and biogeographic history of taxa occurring in a particular region is one way to determine the processes by which the biodiversity of that region originated. Tree boas of the genus Corallus are an ancient clade and occur throughout Central and South America and the Lesser Antilles, making it an excellent group for investigating Neotropical biogeography. Using sequenced portions of two mitochondrial and three nuclear loci for individuals of all recognized species of Corallus, we infer phylogenetic relationships, present the first molecular analysis of the phylogenetic placement of the enigmatic C. cropanii, develop a time-calibrated phylogeny, and explore the biogeographic history of the genus. We found that Corallus diversified within mainland South America, via over-water dispersals to the Lesser Antilles and Central America, and via the traditionally recognized Panamanian land bridge. Divergence time estimates reject the South American Caribbean-Track as a general biogeographic model for Corallus and implicate a role for events during the Oligocene and Miocene in diversification such as marine incursions and the uplift of the Andes. Our findings also suggest that recognition of the island endemic species, C. grenadensis and C. cookii, is questionable as they are nested within the widely distributed species, C. hortulanus. Our results highlight the importance of using widespread taxa when forming and testing biogeographic hypotheses in complex regions and further illustrate the difficulty of forming broadly applicable hypotheses regarding patterns of diversification in the Neotropical region.

Phylogenetic relationships among Staphylococcus species and refinement of cluster groups based on multilocus data

Background

Estimates of relationships among Staphylococcus species have been hampered by poor and inconsistent resolution of phylogenies based largely on single gene analyses incorporating only a limited taxon sample. As such, the evolutionary relationships and hierarchical classification schemes among species have not been confidently established. Here, we address these points through analyses of DNA sequence data from multiple loci (16S rRNA gene, dnaJ, rpoB, and tuf gene fragments) using multiple Bayesian and maximum likelihood phylogenetic approaches that incorporate nearly all recognized Staphylococcus taxa.

Results

We estimated the phylogeny of fifty-seven Staphylococcus taxa using partitioned-model Bayesian and maximum likelihood analysis, as well as Bayesian gene-tree species-tree methods. Regardless of methodology, we found broad agreement among methods that the current cluster groups require revision, although there was some disagreement among methods in resolution of higher order relationships. Based on our phylogenetic estimates, we propose a refined classification for Staphylococcus with species being classified into 15 cluster groups (based on molecular data) that adhere to six species groups (based on phenotypic properties).

Conclusions

Our findings are in general agreement with gene tree-based reports of the staphylococcal phylogeny, although we identify multiple previously unreported relationships among species. Our results support the general importance of such multilocus assessments as a standard in microbial studies to more robustly infer relationships among recognized and newly discovered lineages.

Gravity and the evolution of cardiopulmonary morphology in snakes

Physiological investigations of snakes have established the importance of heart position and pulmonary structure in contexts of gravity effects on blood circulation. Here we investigate morphological correlates of cardiopulmonary physiology in contexts related to ecology, behavior and evolution. We analyze data for heart position and length of vascular lung in 154 species of snakes that exhibit a broad range of characteristic behaviors and habitat associations. We construct a composite phylogeny for these species, and we codify gravitational stress according to species habitat and behavior. We use conventional regression and phylogenetically independent contrasts to evaluate whether trait diversity is correlated with gravitational habitat related to evolutionary transitions within the composite tree topology. We demonstrate that snake species living in arboreal habitats, or which express strongly climbing behaviors, possess relatively short blood columns between the heart and the head, as well as relatively short vascular lungs, compared to terrestrial species. Aquatic species, which experience little or no gravity stress in water, show the reverse - significantly longer heart-head distance and longer vascular lungs. These phylogenetic differences complement the results of physiological studies and are reflected in multiple habitat transitions during the evolutionary histories of these snake lineages, providing strong evidence that heart-to-head distance and length of vascular lung are co-adaptive cardiopulmonary features of snakes.

Snake venomics of Crotalus tigris: the minimalist toxin arsenal of the deadliest Nearctic rattlesnake venom. Evolutionary Clues for generating a pan-specific antivenom against crotalid type II venoms [corrected]

We report the proteomic and antivenomic characterization of Crotalus tigris venom. This venom exhibits the highest lethality for mice among rattlesnakes and the simplest toxin proteome reported to date. The venom proteome of C. tigris comprises 7-8 gene products from 6 toxin families; the presynaptic β-neurotoxic heterodimeric PLA(2), Mojave toxin, and two serine proteinases comprise, respectively, 66 and 27% of the C. tigris toxin arsenal, whereas a VEGF-like protein, a CRISP molecule, a medium-sized disintegrin, and 1-2 PIII-SVMPs each represent 0.1-5% of the total venom proteome. This toxin profile really explains the systemic neuro- and myotoxic effects observed in envenomated animals. In addition, we found that venom lethality of C. tigris and other North American rattlesnake type II venoms correlates with the concentration of Mojave toxin A-subunit, supporting the view that the neurotoxic venom phenotype of crotalid type II venoms may be described as a single-allele adaptation. Our data suggest that the evolutionary trend toward neurotoxicity, which has been also reported for the South American rattlesnakes, may have resulted by pedomorphism. The ability of an experimental antivenom to effectively immunodeplete proteins from the type II venoms of C. tigris, Crotalus horridus , Crotalus oreganus helleri, Crotalus scutulatus scutulatus, and Sistrurus catenatus catenatus indicated the feasibility of generating a pan-American anti-Crotalus type II antivenom, suggested by the identification of shared evolutionary trends among South and North American Crotalus species.

The Great American Biotic Interchange in frogs: multiple and early colonization of Central America by the South American genus Pristimantis (Anura: Craugastoridae)

The completion of the land bridge between North and South America approximately 3.5-3.1 million years ago (Ma) initiated a tremendous biogeographic event called the Great American Biotic Interchange (GABI), described principally from the mammalian fossil record. The history of biotic interchange between continents for taxonomic groups with poor fossil records, however, is not well understood. Molecular and fossil data suggest that a number of plant and animal lineages crossed the Isthmus of Panama well before 3.5 Ma, leading biologists to speculate about trans-oceanic dispersal mechanisms. Here we present a molecular phylogenetic analysis of the frog genus Pristimantis based on 189 individuals of 137 species, including 71 individuals of 31 species from Panama and Colombia. DNA sequence data were obtained from three mitochondrial (COI, 12S, 16S) and two nuclear (RAG-1 and Tyr) genes, for a total of 4074 base pairs. The resulting phylogenetic hypothesis showed statistically significant conflict with most recognized taxonomic groups within Pristimantis, supporting only the rubicundus Species Series, and the Pristimantis myersi and Pristimantis pardalis Species Groups as monophyletic. Inference of ancestral areas based on a likelihood model of geographic range evolution via dispersal, local extinction, and cladogenesis (DEC) suggested that the colonization of Central America by South American Pristimantis involved at least 11 independent events. Relaxed-clock analyses of divergence times suggested that at least eight of these invasions into Central America took place prior to 4 Ma, mainly in the Miocene. These findings contribute to a growing list of molecular-based biogeographic studies presenting apparent temporal conflicts with the traditional GABI model.

Pleistocene climate, phylogeny, and climate envelope models: an integrative approach to better understand species' response to climate change

Mean annual temperature reported by the Intergovernmental Panel on Climate Change increases at least 1.1°C to 6.4°C over the next 90 years. In context, a change in climate of 6°C is approximately the difference between the mean annual temperature of the Last Glacial Maximum (LGM) and our current warm interglacial. Species have been responding to changing climate throughout Earth's history and their previous biological responses can inform our expectations for future climate change. Here we synthesize geological evidence in the form of stable oxygen isotopes, general circulation paleoclimate models, species' evolutionary relatedness, and species' geographic distributions. We use the stable oxygen isotope record to develop a series of temporally high-resolution paleoclimate reconstructions spanning the Middle Pleistocene to Recent, which we use to map ancestral climatic envelope reconstructions for North American rattlesnakes. A simple linear interpolation between current climate and a general circulation paleoclimate model of the LGM using stable oxygen isotope ratios provides good estimates of paleoclimate at other time periods. We use geologically informed rates of change derived from these reconstructions to predict magnitudes and rates of change in species' suitable habitat over the next century. Our approach to modeling the past suitable habitat of species is general and can be adopted by others. We use multiple lines of evidence of past climate (isotopes and climate models), phylogenetic topology (to correct the models for long-term changes in the suitable habitat of a species), and the fossil record, however sparse, to cross check the models. Our models indicate the annual rate of displacement in a clade of rattlesnakes over the next century will be 2 to 3 orders of magnitude greater (430-2,420 m/yr) than it has been on average for the past 320 ky (2.3 m/yr).

Profiling the venom gland transcriptomes of Costa Rican snakes by 454 pyrosequencing

BACKGROUND

A long term research goal of venomics, of applied importance for improving current antivenom therapy, but also for drug discovery, is to understand the pharmacological potential of venoms. Individually or combined, proteomic and transcriptomic studies have demonstrated their feasibility to explore in depth the molecular diversity of venoms. In the absence of genome sequence, transcriptomes represent also valuable searchable databases for proteomic projects.

RESULTS

The venom gland transcriptomes of 8 Costa Rican taxa from 5 genera (Crotalus, Bothrops, Atropoides, Cerrophidion, and Bothriechis) of pitvipers were investigated using high-throughput 454 pyrosequencing. 100,394 out of 330,010 masked reads produced significant hits in the available databases. 5.165,220 nucleotides (8.27%) were masked by RepeatMasker, the vast majority of which corresponding to class I (retroelements) and class II (DNA transposons) mobile elements. BLAST hits included 79,991 matches to entries of the taxonomic suborder Serpentes, of which 62,433 displayed similarity to documented venom proteins. Strong discrepancies between the transcriptome-computed and the proteome-gathered toxin compositions were obvious at first sight. Although the reasons underlaying this discrepancy are elusive, since no clear trend within or between species is apparent, the data indicate that individual mRNA species may be translationally controlled in a species-dependent manner. The minimum number of genes from each toxin family transcribed into the venom gland transcriptome of each species was calculated from multiple alignments of reads matched to a full-length reference sequence of each toxin family. Reads encoding ORF regions of Kazal-type inhibitor-like proteins were uniquely found in Bothriechis schlegelii and B. lateralis transcriptomes, suggesting a genus-specific recruitment event during the early-Middle Miocene. A transcriptome-based cladogram supports the large divergence between A. mexicanus and A. picadoi, and a closer kinship between A. mexicanus and C. godmani.

CONCLUSIONS

Our comparative next-generation sequencing (NGS) analysis reveals taxon-specific trends governing the formulation of the venom arsenal. Knowledge of the venom proteome provides hints on the translation efficiency of toxin-coding transcripts, contributing thereby to a more accurate interpretation of the transcriptome. The application of NGS to the analysis of snake venom transcriptomes, may represent the tool for opening the door to systems venomics.

Ending the drought: new strategies for improving the flow of affordable, effective antivenoms in Asia and Africa

The development of snake antivenoms more than a century ago should have heralded effective treatment of the scourge of snakebite envenoming in impoverished, mostly rural populations around the world. That snakebite still exists today, as a widely untreated illness that maims, kills and terrifies men, women and children in vulnerable communities, is a cruel anachronism. Antivenom can be an effective, safe and affordable treatment for snakebites, but apathy, inaction and the politicisation of public health have marginalised both the problem (making snakebite arguably the most neglected of all neglected tropical diseases) and its solution. For lack of any coordinated approach, provision of antivenoms has been pushed off the public health agenda, leading to an incongruous decline in demand for these crucial antidotes, excused and fed by new priorities, an absence of epidemiological data, and a poor regulatory framework. These factors facilitated the infiltration of poor quality products that degrade user confidence and undermine legitimate producers. The result is that tens of thousands are denied an essential life-saving medicine, allowing a toll of human suffering that is a summation of many individual catastrophes. No strategy has been developed to address this problem and to overcome the intransigence and inaction responsible for the global tragedy of snakebite. Attempts to engage with the broader public health community through the World Health Organisation (WHO), GAVI, and other agencies have failed. Consequently, the toxinology community has taken on a leadership role in a new approach, the Global Snakebite Initiative, which seeks to mobilise the resources, skills and experience of scientists and clinicians for whom venoms, toxins, antivenoms, snakes and snakebites are already fields of interest. Proteomics is one such discipline, which has embraced the potential of using venoms in bio-discovery and systems biology. The fields of venomics and antivenomics have recently evolved from this discipline, offering fresh hope for the victims of snakebites by providing an exciting insight into the complexities, nature, fundamental properties and significance of venom constituents. Such a rational approach brings with it the potential to design new immunising mixtures from which to raise potent antivenoms with wider therapeutic ranges. This addresses a major practical limitation in antivenom use recognised since the beginning of the 20th century: the restriction of therapeutic effectiveness to the specific venom immunogen used in production. Antivenomic techniques enable the interactions between venoms and antivenoms to be examined in detail, and if combined with functional assays of specific activity and followed up by clinical trials of effectiveness and safety, can be powerful tools with which to evaluate the suitability of current and new antivenoms for meeting urgent regional needs. We propose two mechanisms through which the Global Snakebite Initiative might seek to end the antivenom drought in Africa and Asia: first by establishing a multidisciplinary, multicentre, international collaboration to evaluate currently available antivenoms against the venoms of medically important snakes from specific nations in Africa and Asia using a combination of proteomic, antivenomic and WHO-endorsed preclinical assessment protocols, to provide a validated evidence base for either recommending or rejecting individual products; and secondly by bringing the power of proteomics to bear on the design of new immunising mixtures to raise Pan-African and Pan-Asian polyvalent antivenoms of improved potency and quality. These products will be subject to rigorous clinical assessment. We propose radically to change the basis upon which antivenoms are produced and supplied for the developing world. Donor funding and strategic public health alliances will be sought to make it possible not only to sustain the financial viability of antivenom production partnerships, but also to ensure that patients are relieved of the costs of antivenom so that poverty is no longer a barrier to the treatment of this important, but grossly neglected public health emergency.

Enzymatic toxins from snake venom: structural characterization and mechanism of catalysis

Snake venoms are cocktails of enzymes and non-enzymatic proteins used for both the immobilization and digestion of prey. The most common snake venom enzymes include acetylcholinesterases, l-amino acid oxidases, serine proteinases, metalloproteinases and phospholipases A(2) . Higher catalytic efficiency, thermal stability and resistance to proteolysis make these enzymes attractive models for biochemists, enzymologists and structural biologists. Here, we review the structures of these enzymes and describe their structure-based mechanisms of catalysis and inhibition. Some of the enzymes exist as protein complexes in the venom. Thus we also discuss the functional role of non-enzymatic subunits and the pharmacological effects of such protein complexes. The structures of inhibitor-enzyme complexes provide ideal platforms for the design of potent inhibitors which are useful in the development of prototypes and lead compounds with potential therapeutic applications.