Comparative venomics of the Prairie Rattlesnake (Crotalus viridis viridis) from Colorado: Identification of a novel pattern of ontogenetic changes in venom composition and assessment of the immunoreactivity of the commercial antivenom CroFab®.

Here we describe and compare the venomic and antivenomic characteristics of both neonate and adult Prairie Rattlesnake (Crotalus viridis viridis) venoms. Although both neonate and adult venoms contain unique components, similarities among protein family content were seen. Both neonate and adult venoms consisted of myotoxin, bradykinin-potentiating peptide (BPP), phospholipase A2 (PLA2), Zn(2+)-dependent metalloproteinase (SVMP), serine proteinase, L-amino acid oxidase (LAAO), cysteine-rich secretory protein (CRISP) and disintegrin families. Quantitative differences, however, were observed, with venoms of adults containing significantly higher concentrations of the non-enzymatic toxic compounds and venoms of neonates containing higher concentrations of pre-digestive enzymatic proteins such as SVMPs. To assess the relevance of this venom variation in the context of snakebite and snakebite treatment, we tested the efficacy of the common antivenom CroFab® for recognition of both adult and neonate venoms in vitro. This comparison revealed that many of the major protein families (SVMPs, CRISP, PLA2, serine proteases, and LAAO) in both neonate and adult venoms were immunodepleted by the antivenom, whereas myotoxins, one of the major toxic components of C. v. viridis venom, in addition to many of the small peptides, were not efficiently depleted by CroFab®. These results therefore provide a comprehensive catalog of the venom compounds present in C. v. viridis venom and new molecular insight into the potential efficacy of CroFab® against human envenomations by one of the most widely distributed rattlesnake species in North America. BIOLOGICAL SIGNIFICANCE: Comparative proteomic analysis of venoms of neonate and adult Prairie Rattlesnake (Crotalus viridis viridis) from a discrete population in Colorado revealed a novel pattern of ontogenetic shifts in toxin composition for viperid snakes. The observed stage-dependent decrease of the relative content of disintegrins, catalytically active D49-PLA2s, L-amino acid oxidase, and SVMPs, and the concomitant increase of the relative abundance of paralytic small basic myotoxins and ohanin-like toxin, and hemostasis-disrupting serine proteinases, may represent an age-dependent strategy for securing prey and avoiding injury as the snake switches from small ectothermic prey and newborn rodents to larger endothermic prey. Such age-dependent shifts in venom composition may be relevant for antivenom efficacy and treatment of snakebite. However, applying a second-generation antivenomics approach, we show that CroFab®, developed against venom of three Crotalus and one Agkistrodon species, efficiently immunodepleted many, but not all, of the major compounds present in neonate and adult C. v. viridis venoms.

Phylogenetics, likelihood, evolution and complexity.

SUMMARY: Phylogenetics, likelihood, evolution and complexity (PLEX) is a flexible and fast Bayesian Markov chain Monte Carlo software program for large-scale analysis of nucleotide and amino acid data using complex evolutionary models in a phylogenetic framework. The program gains large speed improvements over standard approaches by implementing 'partial sampling of substitution histories', a data augmentation approach that can reduce data analysis times from months to minutes on large comparative datasets. A variety of nucleotide and amino acid substitution models are currently implemented, including non-reversible and site-heterogeneous mixture models. Due to efficient algorithms that scale well with data size and model complexity, PLEX can be used to make inferences from hundreds to thousands of taxa in only minutes on a desktop computer. It also performs probabilistic ancestral sequence reconstruction. Future versions will support detection of co-evolutionary interactions between sites, probabilistic tests of convergent evolution and rigorous testing of evolutionary hypotheses in a Bayesian framework. AVAILABILITY AND IMPLEMENTATION: PLEX v1.0 is licensed under GPL. Source code and documentation will be available for download at www.evolutionarygenomics.com/ProgramsData/PLEX. PLEX is implemented in C++ and supported on Linux, Mac OS X and other platforms supporting standard C++ compilers. Example data, control files, documentation and accessory Perl scripts are available from the website. CONTACT: David.Pollock@UCDenver.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Rapid microsatellite identification from Illumina paired-end genomic sequencing in two birds and a snake.

Identification of microsatellites, or simple sequence repeats (SSRs), can be a time-consuming and costly investment requiring enrichment, cloning, and sequencing of candidate loci. Recently, however, high throughput sequencing (with or without prior enrichment for specific SSR loci) has been utilized to identify SSR loci. The direct "Seq-to-SSR" approach has an advantage over enrichment-based strategies in that it does not require a priori selection of particular motifs, or prior knowledge of genomic SSR content. It has been more expensive per SSR locus recovered, however, particularly for genomes with few SSR loci, such as bird genomes. The longer but relatively more expensive 454 reads have been preferred over less expensive Illumina reads. Here, we use Illumina paired-end sequence data to identify potentially amplifiable SSR loci (PALs) from a snake (the Burmese python, Python molurus bivittatus), and directly compare these results to those from 454 data. We also compare the python results to results from Illumina sequencing of two bird genomes (Gunnison Sage-grouse, Centrocercus minimus, and Clark's Nutcracker, Nucifraga columbiana), which have considerably fewer SSRs than the python. We show that direct Illumina Seq-to-SSR can identify and characterize thousands of potentially amplifiable SSR loci for as little as $10 per sample--a fraction of the cost of 454 sequencing. Given that Illumina Seq-to-SSR is effective, inexpensive, and reliable even for species such as birds that have few SSR loci, it seems that there are now few situations for which prior hybridization is justifiable.

Bayesian analysis of high-throughput quantitative measurement of protein-DNA interactions.

Transcriptional regulation depends upon the binding of transcription factor (TF) proteins to DNA in a sequence-dependent manner. Although many experimental methods address the interaction between DNA and proteins, they generally do not comprehensively and accurately assess the full binding repertoire (the complete set of sequences that might be bound with at least moderate strength). Here, we develop and evaluate through simulation an experimental approach that allows simultaneous high-throughput quantitative analysis of TF binding affinity to thousands of potential DNA ligands. Tens of thousands of putative binding targets can be mixed with a TF, and both the pre-bound and bound target pools sequenced. A hierarchical Bayesian Markov chain Monte Carlo approach determines posterior estimates for the dissociation constants, sequence-specific binding energies, and free TF concentrations. A unique feature of our approach is that dissociation constants are jointly estimated from their inferred degree of binding and from a model of binding energetics, depending on how many sequence reads are available and the explanatory power of the energy model. Careful experimental design is necessary to obtain accurate results over a wide range of dissociation constants. This approach, which we call Simultaneous Ultra high-throughput Ligand Dissociation EXperiment (SULDEX), is theoretically capable of rapid and accurate elucidation of an entire TF-binding repertoire.