The genetic regulatory architecture and epigenomic basis for age-related changes in rattlesnake venom

Developmental phenotypic changes can evolve under selection imposed by age- and size-related ecological differences. Many of these changes occur through programmed alterations to gene expression patterns, but the molecular mechanisms and gene-regulatory networks underlying these adaptive changes remain poorly understood. Many venomous snakes, including the eastern diamondback rattlesnake (Crotalus adamanteus), undergo correlated changes in diet and venom expression as snakes grow larger with age, providing models for identifying mechanisms of timed expression changes that underlie adaptive life history traits. By combining a highly contiguous, chromosome-level genome assembly with measures of expression, chromatin accessibility, and histone modifications, we identified cis-regulatory elements and trans-regulatory factors controlling venom ontogeny in the venom glands of C. adamanteus. Ontogenetic expression changes were significantly correlated with epigenomic changes within genes, immediately adjacent to genes (e.g., promoters), and more distant from genes (e.g., enhancers). We identified 37 candidate transcription factors (TFs), with the vast majority being up-regulated in adults. The ontogenetic change is largely driven by an increase in the expression of TFs associated with growth signaling, transcriptional activation, and circadian rhythm/biological timing systems in adults with corresponding epigenomic changes near the differentially expressed venom genes. However, both expression activation and repression contributed to the composition of both adult and juvenile venoms, demonstrating the complexity and potential evolvability of gene regulation for this trait. Overall, given that age-based trait variation is common across the tree of life, we provide a framework for understanding gene-regulatory-network-driven life-history evolution more broadly.

A genome assembly for the southern Pacific rattlesnake, Crotalus oreganus helleri, in the western rattlesnake species complex

Rattlesnakes play important roles in their ecosystems by regulating prey populations, are involved in complex coevolutionary dynamics with their prey, and exhibit a variety of unusual adaptations, including maternal care, heat-sensing pit organs, hinged fangs, and medically-significant venoms. The western rattlesnake (Crotalus oreganus) is one of the widest ranging rattlesnake species, with a distribution from British Columbia, where it is listed as threatened, to Baja California and east across the Great Basin to western Wyoming, Colorado and New Mexico. Here, we report a new reference genome assembly for one of six currently recognized subspecies, C. oreganus helleri, as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genomic sequencing strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises a total of 698 scaffolds spanning 1,564,812,557 base pairs, has a contig N50 of 64.7 Mb, a scaffold N50 of 110.8 Mb, and BUSCO complete score of 90.5%. This reference genome will be valuable for studies on the genomic basis of venom evolution and variation within Crotalus, in resolving the taxonomy of C. oreganus and its relatives, and for the conservation and management of rattlesnakes in general.

Functional genomic diversity is correlated with neutral genomic diversity in populations of an endangered rattlesnake

Theory predicts that genetic erosion in small, isolated populations of endangered species can be assessed using estimates of neutral genetic variation, yet this widely used approach has recently been questioned in the genomics era. Here, we leverage a chromosome-level genome assembly of an endangered rattlesnake (Sistrurus catenatus) combined with whole genome resequencing data (N = 110 individuals) to evaluate the relationship between levels of genome-wide neutral and functional diversity over historical and future timescales. As predicted, we found positive correlations between genome-wide estimates of neutral genetic diversity (π) and inferred levels of adaptive variation and an estimate of inbreeding mutation load, and a negative relationship between neutral diversity and an estimate of drift mutation load. However, these correlations were half as strong for projected future levels of neutral diversity based on contemporary effective population sizes. Broadly, our results confirm that estimates of neutral genetic diversity provide an accurate measure of genetic erosion in populations of a threatened vertebrate. They also provide nuance to the neutral-functional diversity controversy by suggesting that while these correlations exist, anthropogenetic impacts may have weakened these associations in the recent past and into the future.

The best of both worlds? Rattlesnake hybrid zones generate complex combinations of divergent venom phenotypes that retain high toxicity

Studying the consequences of hybridization between closely related species with divergent traits can reveal patterns of evolution that shape and maintain extreme trophic adaptations. Snake venoms are an excellent model system for examining the evolutionary and ecological patterns that underlie highly selected polymorphic traits. Here we investigate hybrid venom phenotypes that result from natural introgression between two rattlesnake species that express highly divergent venom phenotypes: Crotalus o. concolor and C. v. viridis. Though not yet documented, interbreeding between these species may lead to novel venom phenotypes with unique activities that break the typical trends of venom composition in rattlesnakes. The characteristics of these unusual phenotypes could unveil the roles of introgression in maintaining patterns of venom composition and variation, including the near ubiquitous dichotomy between neurotoxic or degradative venoms observed across rattlesnakes. We use RADseq data to infer patterns of gene flow and hybrid ancestry between these diverged lineages and link these genetic data with analyses of venom composition, biological activity, and whole animal model toxicity tests to understand the impacts of introgression on venom composition. We find that introgressed populations express admixed venom phenotypes that do not sacrifice biological activity (lethal toxicity) or overall abundance of dominant toxins compared to parental venoms. These hybridized venoms therefore do not represent a trade-off in functionality between the typical phenotypic extremes but instead represent a unique combination of characters whose expression appears limited to the hybrid zone.

Inferring population connectivity in eastern massasauga rattlesnakes (Sistrurus catenatus) using landscape genetics

Assessing the environmental factors that influence the ability of a threatened species to move through a landscape can be used to identify conservation actions that connect isolated populations. However, direct observations of species' movement are often limited, making the development of alternate approaches necessary. Here we use landscape genetic analyses to assess the impact of landscape features on the movement of individuals between local populations of a threatened snake, the eastern massasauga rattlesnake (Sistrurus catenatus). We linked connectivity data with habitat information from two landscapes of similar size: a large region of unfragmented habitat and a previously studied fragmented landscape consisting of isolated patches of habitat. We used this analysis to identify features of the landscape where modification or acquisition would enhance population connectivity in the fragmented region. We found evidence that current connectivity was impacted by both contemporary land-cover features, especially roads, and inherent landscape features such as elevation. Next, we derived estimates of expected movement ability using a recently developed pedigree-based approach and least-cost paths through the unfragmented landscape. We then used our pedigree and resistance map to estimate resistance polygons of the potential extent for S. catenatus movement in the fragmented landscape. These polygons identify possible sites for future corridors connecting currently isolated populations in this landscape by linking the impact of future habitat modification or land acquisition to dispersal ability in this species. Overall, our study shows how modeling landscape resistance across differently fragmented landscapes can identify habitat features that affect contemporary movement in threatened species in fragmented landscapes and how this information can be used to guide mitigation actions whose goal is to connect isolated populations.

Montane Rattlesnakes in México: Venoms of Crotalus tancitarensis and Related Species within the Crotalus intermedius Group

The Crotalus intermedius group is a clade of rattlesnakes consisting of several species adapted to a high elevation habitat, primarily in México. Crotalus tancitarensis was previously classified as C. intermedius, until individuals occurring on Cerro Tancítaro in Michoacán, México, were reevaluated and classified as a new species (C. tancitarensis) based on scale pattern and geographic location. This study aimed to characterize the venom of C. tancitarensis and compare the venom profile to those of other species within the Crotalus intermedius group using gel electrophoresis, biochemical assays, reverse-phase high performance liquid chromatography, mass spectrometry, and lethal toxicity (LD50) assays. Results show that the venom profiles of species within the Crotalus intermedius group are similar, but with distinct differences in phospholipase A2 (PLA2), metalloproteinase PI (SVMP PI), and kallikrein-like serine proteinase (SVSP) activity and relative abundance. Proteomic analysis indicated that the highland forms produce venoms with 50-60 protein isoforms and a composition typical of type I rattlesnake venoms (abundant SVMPs, lack of presynaptic PLA2-based neurotoxins), as well as a diversity of typical Crotalus venom components such as serine proteinases, PLA2s, C-type lectins, and less abundant toxins (LAAOs, CRiSPs, etc.). The overall venom profile of C. tancitarensis appears most similar to C. transversus, which is consistent with a previous mitochondrial DNA analysis of the Crotalus intermedius group. These rattlesnakes of the Mexican highlands represent a radiation of high elevation specialists, and in spite of divergence of species in these Sky Island habitats, venom composition of species analyzed here has remained relatively conserved. The majority of protein family isoforms are conserved in all members of the clade, and as seen in other more broadly distributed rattlesnake species, differences in their venoms are largely due to relative concentrations of specific components.

A Promising Biomolecule Able to Degrade Neutrophil Extracellular Traps: CdcPDE, a Rattlesnake Phosphodiesterase

Neutrophil extracellular traps (NETs) are an important mechanism for defense against pathogens. Their overproduction can be harmful since excessive NET formation promotes inflammation and tissue damage in several diseases. Nucleases are capable to degrade NET on basis of their DNA hydrolysis activity, including the CdcPDE, a nuclease isolated from Crotalus durissus collilineatus snake venom. Here, we report a new finding about CdcPDE activity, demonstrating its efficiency in degrading cell-free DNA from NETs, being a potential candidate to assist in therapies targeting inflammatory diseases.

The roles of balancing selection and recombination in the evolution of rattlesnake venom

The origin of snake venom involved duplication and recruitment of non-venom genes into venom systems. Several studies have predicted that directional positive selection has governed this process. Venom composition varies substantially across snake species and venom phenotypes are locally adapted to prey, leading to coevolutionary interactions between predator and prey. Venom origins and contemporary snake venom evolution may therefore be driven by fundamentally different selection regimes, yet investigations of population-level patterns of selection have been limited. Here, we use whole-genome data from 68 rattlesnakes to test hypotheses about the factors that drive genomic diversity and differentiation in major venom gene regions. We show that selection has resulted in long-term maintenance of genetic diversity within and between species in multiple venom gene families. Our findings are inconsistent with a dominant role of directional positive selection and instead support a role of long-term balancing selection in shaping venom evolution. We also detect rapid decay of linkage disequilibrium due to high recombination rates in venom regions, suggesting that venom genes have reduced selective interference with nearby loci, including other venom paralogues. Our results provide an example of long-term balancing selection that drives trans-species polymorphism and help to explain how snake venom keeps pace with prey resistance.

A Meta-Analysis of the Protein Components in Rattlesnake Venom

The specificity and potency of venom components give them a unique advantage in developing various pharmaceutical drugs. Though venom is a cocktail of proteins, rarely are the synergy and association between various venom components studied. Understanding the relationship between various components of venom is critical in medical research. Using meta-analysis, we observed underlying patterns and associations in the appearance of the toxin families. For Crotalus, Dis has the most associations with the following toxins: PDE; BPP; CRL; CRiSP; LAAO; SVMP P-I and LAAO; SVMP P-III and LAAO. In Sistrurus venom, CTL and NGF have the most associations. These associations can predict the presence of proteins in novel venom and understand synergies between venom components for enhanced bioactivity. Using this approach, the need to revisit the classification of proteins as major components or minor components is highlighted. The revised classification of venom components is based on ubiquity, bioactivity, the number of associations, and synergies. The revised classification can be expected to trigger increased research on venom components, such as NGF, which have high biomedical significance. Using hierarchical clustering, we observed that the genera's venom compositions were similar, based on functional characteristics rather than phylogenetic relationships.

The Tiger Rattlesnake genome reveals a complex genotype underlying a simple venom phenotype

Variation in gene regulation is ubiquitous, yet identifying the mechanisms producing such variation, especially for complex traits, is challenging. Snake venoms provide a model system for studying the phenotypic impacts of regulatory variation in complex traits because of their genetic tractability. Here, we sequence the genome of the Tiger Rattlesnake, which possesses the simplest and most toxic venom of any rattlesnake species, to determine whether the simple venom phenotype is the result of a simple genotype through gene loss or a complex genotype mediated through regulatory mechanisms. We generate the most contiguous snake-genome assembly to date and use this genome to show that gene loss, chromatin accessibility, and methylation levels all contribute to the production of the simplest, most toxic rattlesnake venom. We provide the most complete characterization of the venom gene-regulatory network to date and identify key mechanisms mediating phenotypic variation across a polygenic regulatory network.

Gradual and Discrete Ontogenetic Shifts in Rattlesnake Venom Composition and Assessment of Hormonal and Ecological Correlates

Ontogenetic shifts in venom occur in many snakes but establishing their nature as gradual or discrete processes required additional study. We profiled shifts in venom expression from the neonate to adult sizes of two rattlesnake species, the eastern diamondback and the timber rattlesnake. We used serial sampling and venom chromatographic profiling to test if ontogenetic change occurs gradually or discretely. We found evidence for gradual shifts in overall venom composition in six of eight snakes, which sometimes spanned more than two years. Most chromatographic peaks shift gradually, but one quarter shift in a discrete fashion. Analysis of published diet data showed gradual shifts in overall diet composition across the range of body sizes attained by our eight study animals, while the shifts in abundance of different prey classes varied in form from gradual to discrete. Testosterone concentrations were correlated with the change in venom protein composition, but the relationship is not strong enough to suggest causation. Venom research employing simple juvenile versus adult size thresholds may be failing to account for continuous variation in venom composition lifespan. Our results imply that venom shifts represent adaptive matches to dietary shifts and highlight venom for studies of alternative gene regulatory mechanisms.

The Sequence and a Three-Dimensional Structural Analysis Reveal Substrate Specificity Among Snake Venom Phosphodiesterases

(1) Background. Snake venom phosphodiesterases (SVPDEs) are among the least studied venom enzymes. In envenomation, they display various pathological effects, including induction of hypotension, inhibition of platelet aggregation, edema, and paralysis. Until now, there have been no 3D structural studies of these enzymes, thereby preventing structure-function analysis. To enable such investigations, the present work describes the model-based structural and functional characterization of a phosphodiesterase from Crotalusadamanteus venom, named PDE_Ca. (2) Methods. The PDE_Ca structure model was produced and validated using various software (model building: I-TESSER, MODELLER 9v19, Swiss-Model, and validation tools: PROCHECK, ERRAT, Molecular Dynamic Simulation, and Verif3D). (3) Results. The proposed model of the enzyme indicates that the 3D structure of PDE_Ca comprises four domains, a somatomedin B domain, a somatomedin B-like domain, an ectonucleotide pyrophosphatase domain, and a DNA/RNA non-specific domain. Sequence and structural analyses suggest that differences in length and composition among homologous snake venom sequences may account for their differences in substrate specificity. Other properties that may influence substrate specificity are the average volume and depth of the active site cavity. (4) Conclusion. Sequence comparisons indicate that SVPDEs exhibit high sequence identity but comparatively low identity with mammalian and bacterial PDEs.

Phylogenetic origins for severe acetaminophen toxicity in snake species compared to other vertebrate taxa

While it has been known for a while that some snake species are extremely sensitive to acetaminophen, the underlying mechanism for this toxicity has not been reported. To investigate if essential detoxification enzymes are missing in snake species that are responsible for biotransformation of acetaminophen in other vertebrate species, livers were collected from a variety of snake species, together with samples from alligator, snapping turtle, cat, rat, and cattle. Subcellular fractions were analyzed for enzymatic activities of phenol-type sulfotransferase and UDP‑glucuronosyltransferase, total glutathione S‑transferase, and N‑acetyltransferase. The results showed that none of the snake species, together with the cat samples, had any phenol-type glucuronidation activity, and that this activity was much lower in alligator and turtle samples than in the mammalian species. Combined with the lack of N‑acetyltransferase activity in snakes and cats, this would explain the accumulation of the aminophenol metabolite, which induces methemoglobinemia and subsequent suffocation of snakes and cats after acetaminophen exposure. While previous investigations have concluded that in cats the gene for the phenol-type glucuronosyltransferase isoform has turned into a pseudogene because of several point mutations, evaluation of genomic information for snake species revealed that they have only 2 genes that may code for glucuronosyltransferase isoforms. Similarity of these genes with mammalian genes is <50%, and suggests that the expressed enzymes may act on other types of substrates than aromatic amines. This indicates that the extreme sensitivity for acetaminophen in snakes is based on a different phylogenetic origin than the sensitivity observed in cats.

Genetic diversity and structure of Crotalus triseriatus, a rattlesnake of central Mexico

The isolated and fragmented populations are highly susceptible to stochastic events, increasing the extinction risk because of the decline in putative adaptive potential and individual fitness. The population has high heterozygosity values and a moderate allelic diversity, the heterozygosity values are higher than in most other Crotalus species and snake studies. Possibly these high levels of genetic diversity can be related to a large founder size, high effective population size, multiple paternity and overlapping generations. We did not find the genetic structuring but the effective number of alleles (N_e) was 138.1. We found evidence of bottlenecks and the majority of rattlesnakeswere unrelated, despite the small sample size, endemic status, the isolated and fragmented habitat. The genetic information provided in this study can be useful as a first approach to try to make informed conservation efforts for this species and also, important to preserve the habitat of this species; the endangered Abies-Pinus forest of the Nevado the Toluca Volcano.

Is Hybridization a Source of Adaptive Venom Variation in Rattlesnakes? A Test, Using a Crotalus scutulatus × viridis Hybrid Zone in Southwestern New Mexico

Venomous snakes often display extensive variation in venom composition both between and within species. However, the mechanisms underlying the distribution of different toxins and venom types among populations and taxa remain insufficiently known. Rattlesnakes (Crotalus, Sistrurus) display extreme inter- and intraspecific variation in venom composition, centered particularly on the presence or absence of presynaptically neurotoxic phospholipases A₂ such as Mojave toxin (MTX). Interspecific hybridization has been invoked as a mechanism to explain the distribution of these toxins across rattlesnakes, with the implicit assumption that they are adaptively advantageous. Here, we test the potential of adaptive hybridization as a mechanism for venom evolution by assessing the distribution of genes encoding the acidic and basic subunits of Mojave toxin across a hybrid zone between MTX-positive Crotalus scutulatus and MTX-negative C. viridis in southwestern New Mexico, USA. Analyses of morphology, mitochondrial and single copy-nuclear genes document extensive admixture within a narrow hybrid zone. The genes encoding the two MTX subunits are strictly linked, and found in most hybrids and backcrossed individuals, but not in C. viridis away from the hybrid zone. Presence of the genes is invariably associated with presence of the corresponding toxin in the venom. We conclude that introgression of highly lethal neurotoxins through hybridization is not necessarily favored by natural selection in rattlesnakes, and that even extensive hybridization may not lead to introgression of these genes into another species.

Deconstructing a Species-Complex: Geometric Morphometric and Molecular Analyses Define Species in the Western Rattlesnake (Crotalus viridis)

Morphological data are a conduit for the recognition and description of species, and their acquisition has recently been broadened by geometric morphometric (GM) approaches that co-join the collection of digital data with exploratory 'big data' analytics. We employed this approach to dissect the Western Rattlesnake (Crotalus viridis) species-complex in North America, currently partitioned by mitochondrial (mt)DNA analyses into eastern and western lineages (two and seven subspecies, respectively). The GM data (i.e., 33 dorsal and 50 lateral head landmarks) were gleaned from 2,824 individuals located in 10 museum collections. We also downloaded and concatenated sequences for six mtDNA genes from the NCBI GenBank database. GM analyses revealed significant head shape differences attributable to size and subspecies-designation (but not their interactions). Pairwise shape distances among subspecies were significantly greater than those derived from ancestral character states via squared-change parsimony, with the greatest differences separating those most closely related. This, in turn, suggests the potential for historic character displacement as a diversifying force in the complex. All subspecies, save one, were significantly differentiated in a Bayesian discriminant function analysis (DFA), regardless of whether our priors were uniform or informative (i.e., mtDNA data). Finally, shape differences among sister-clades were significantly greater than expected by chance alone under a Brownian model of evolution, promoting the hypothesis that selection rather than drift was the driving force in the evolution of the complex. Lastly, we combine head shape and mtDNA data so as to derived an integrative taxonomy that produced robust boundaries for six OTUs (operational taxonomic units) of the C. viridis complex. We suggest these boundaries are concomitant with species-status and subsequently provide a relevant nomenclature for its recognition and representation.

Limitations of climatic data for inferring species boundaries: insights from speckled rattlesnakes

Phenotypes, DNA, and measures of ecological differences are widely used in species delimitation. Although rarely defined in such studies, ecological divergence is almost always approximated using multivariate climatic data associated with sets of specimens (i.e., the "climatic niche"); the justification for this approach is that species-specific climatic envelopes act as surrogates for physiological tolerances. Using identical statistical procedures, we evaluated the usefulness and validity of the climate-as-proxy assumption by comparing performance of genetic (nDNA SNPs and mitochondrial DNA), phenotypic, and climatic data for objective species delimitation in the speckled rattlesnake (Crotalus mitchellii) complex. Ordination and clustering patterns were largely congruent among intrinsic (heritable) traits (nDNA, mtDNA, phenotype), and discordance is explained by biological processes (e.g., ontogeny, hybridization). In contrast, climatic data did not produce biologically meaningful clusters that were congruent with any intrinsic dataset, but rather corresponded to regional differences in atmospheric circulation and climate, indicating an absence of inherent taxonomic signal in these data. Surrogating climate for physiological tolerances adds artificial weight to evidence of species boundaries, as these data are irrelevant for that purpose. Based on the evidence from congruent clustering of intrinsic datasets, we recommend that three subspecies of C. mitchellii be recognized as species: C. angelensis, C. mitchellii, and C. Pyrrhus.

Quartet inference from SNP data under the coalescent model

MOTIVATION

Increasing attention has been devoted to estimation of species-level phylogenetic relationships under the coalescent model. However, existing methods either use summary statistics (gene trees) to carry out estimation, ignoring an important source of variability in the estimates, or involve computationally intensive Bayesian Markov chain Monte Carlo algorithms that do not scale well to whole-genome datasets.

RESULTS

We develop a method to infer relationships among quartets of taxa under the coalescent model using techniques from algebraic statistics. Uncertainty in the estimated relationships is quantified using the nonparametric bootstrap. The performance of our method is assessed with simulated data. We then describe how our method could be used for species tree inference in larger taxon samples, and demonstrate its utility using datasets for Sistrurus rattlesnakes and for soybeans.

AVAILABILITY AND IMPLEMENTATION

The method to infer the phylogenetic relationship among quartets is implemented in the software SVDquartets, available at www.stat.osu.edu/∼lkubatko/software/SVDquartets.

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.

Immunohistochemical and functional characterization of nitric oxide signaling pathway in isolated aorta from Crotalus durissus terrificus

We characterized the nitric oxide (NO)-cyclic GMP-phosphodiesterase-5 (PDE5) pathway in Crotalus durissus terrificus aorta. Concentration responses curves to acetylcholine (ACh), sodium nitroprusside (SNP), BAY41-2272 (soluble guanylyl cyclase [sGC] stimulator), BAY60-2770 (sGC activator) and tadalafil (PDE5 inhibitor) were constructed in phenylephrine (10 μM)-precontracted tissues with intact (E(+)) or denuded (E(-)) endothelium. ACh (0.0001-10 μM) and SNP (0.0001-10 μM) relaxed aorta, which were reduced by the NO synthase (L-NAME,100 μM) or the sGC inhibitors (ODQ, 10 μM). Tadalafil (0.0001-10 μM) relaxed E(+) rings with potency (pEC(50)) and maximal response (E(max)) values of 7.34±0.02 and 105±8%, respectively. E(-) or ODQ treatment significantly (P<0.05) reduced tadalafil relaxations (66±18% and 71±7%, respectively). BAY41-2272 (0.0001-300 nM) produced concentration-dependent relaxations in E(+) rings, which were reduced by addition of either ODQ or L-NAME (16.0- and 5.2-fold rightward shifts, respectively). The relaxation of BAY60-2770 was markedly potentiated by ODQ and L-NAME (41.0- and 9.7-fold leftward shifts, respectively), whereas in E(-) the pEC(50) values were shifted by 7-fold to the right. Immunohistochemistry, followed validation by transcriptomic analysis, revealed the presence of eNOS in endothelium, whereas nNOS was observed only in perivascular nerves. sGC and PDE5 were expressed in smooth muscle. Thus, NO-sGC-PDE5 pathway is evolutionarily present in Crotalus sp. vessels, and has a remarkable degree of functional similarity to mammalian vessels.

Predation and the maintenance of color polymorphism in a habitat specialist squamate

Multiple studies have addressed the mechanisms maintaining polymorphism within a population. However, several examples exist where species inhabiting diverse habitats exhibit local population-specific polymorphism. Numerous explanations have been proposed for the maintenance of geographic variation in color patterns. For example, spatial variation in patterns of selection or limited gene flow can cause entire populations to become fixed for a single morph, resulting in separate populations of the same species exhibiting separate and distinct color morphs. The mottled rock rattlesnake (Crotalus lepidus lepidus) is a montane species that exhibits among-population color polymorphism that correlates with substrate color. Habitat substrate in the eastern part of its range is composed primarily of light colored limestone and snakes have light dorsal coloration, whereas in the western region the substrate is primarily dark and snakes exhibit dark dorsal coloration. We hypothesized that predation on high contrast color and blotched patterns maintain these distinct color morphs. To test this we performed a predation experiment in the wild by deploying model snakes at 12 sites evenly distributed within each of the two regions where the different morphs are found. We employed a 2×2 factorial design that included two color and two blotched treatments. Our results showed that models contrasting with substrate coloration suffered significantly more avian attacks relative to models mimicking substrates. Predation attempts on blotched models were similar in each substrate type. These results support the hypothesis that color pattern is maintained by selective predation.

Roads, interrupted dispersal, and genetic diversity in timber rattlesnakes

Anthropogenic habitat modification often creates barriers to animal movement, transforming formerly contiguous habitat into a patchwork of habitat islands with low connectivity. Roadways are a feature of most landscapes that can act as barriers or filters to migration among local populations. Even small and recently constructed roads can have a significant impact on population genetic structure of some species, but not others. We developed a research approach that combines fine-scale molecular genetics with behavioral and ecological data to understand the impacts of roads on population structure and connectivity. We used microsatellite markers to characterize genetic variation within and among populations of timber rattlesnakes (Crotalus horridus) occupying communal hibernacula (dens) in regions bisected by roadways. We examined the impact of roads on seasonal migration, genetic diversity, and gene flow among populations. Snakes in hibernacula isolated by roads had significantly lower genetic diversity and higher genetic differentiation than snakes in hibernacula in contiguous habitat. Genetic-assignment analyses revealed that interruption to seasonal migration was the mechanism underlying these patterns. Our results underscore the sizeable impact of roads on this species, despite their relatively recent construction at our study sites (7 to 10 generations of rattlesnakes), the utility of population genetics for studies of road ecology, and the need for mitigating effects of roads.

Crotalus durissus collilineatus venom gland transcriptome: analysis of gene expression profile

Crotalus durissus rattlesnakes are responsible for the most lethal cases of snakebites in Brazil. Crotalus durissus collilineatus subspecies is related to a great number of accidents in Southeast and Central West regions, but few studies on its venom composition have been carried out to date. In an attempt to describe the transcriptional profile of the C. durissus collilineatus venom gland, we generated a cDNA library and the sequences obtained could be identified by similarity searches on existing databases. Out of 673 expressed sequence tags (ESTs) 489 produced readable sequences comprising 201 singletons and 47 clusters of two or more ESTs. One hundred and fifty reads (60.5%) produced significant hits to known sequences. The results showed a predominance of toxin-coding ESTs instead of transcripts coding for proteins involved in all cellular functions. The most frequent toxin was crotoxin, comprising 88% of toxin-coding sequences. Crotoxin B, a basic phospholipase A(2) (PLA(2)) subunit of crotoxin, was represented in more variable forms comparing to the non-enzymatic subunit (crotoxin A), and most sequences coding this molecule were identified as CB1 isoform from Crotalus durissus terrificus venom. Four percent of toxin-related sequences in this study were identified as growth factors, comprising five sequences for vascular endothelial growth factor (VEGF) and one for nerve growth factor (NGF) that showed 100% of identity with C. durissus terrificus NGF. We also identified two clusters for metalloprotease from PII class comprising 3% of the toxins, and two for serine proteases, including gyroxin (2.5%). The remaining 2.5% of toxin-coding ESTs represent singletons identified as homologue sequences to cardiotoxin, convulxin, angiotensin-converting enzyme inhibitor and C-type natriuretic peptide, Ohanin, crotamin and PLA(2) inhibitor. These results allowed the identification of the most common classes of toxins in C. durissus collilineatus snake venom, also showing some unknown classes for this subspecies and even for C. durissus species, such as cardiotoxins and VEGF.

Soft epidermis of a scaleless snake lacks beta-keratin

Beta-keratins are responsible for the mechanical resistance of scales in reptiles. In a scaleless crotalus snake (Crotalus atrox), large areas of the skin are completely devoid of scales, and the skin appears delicate and wrinkled. The epidermis of this snake has been assessed for the presence of beta-keratin by immunocytochemistry and immunoblotting using an antibody against chicken scale beta-keratin. This antibody recognizes beta-keratins in normal snake scales with molecular weights of 15-18 kDa and isoelectric points at 6.8, 7.5, 8.3 and 9.4. This indicates that beta-keratins of the stratum corneum are mainly basic proteins, so may interact with cytokeratins of the epidermis, most of which appear acidic (isoelectric points 4.5-5.5). A beta-layer and beta-keratin immunoreactivity are completely absent in moults of the scaleless mutant, and the corneous layer comprises a multi-layered alpha-layer covered by a flat oberhautchen. In conclusion, the present study shows that a lack of beta-keratins is correlated with the loss of scales and mechanical protection in the skin of this mutant snake.

Molecular evolution of PIII-SVMP and RGD disintegrin genes from the genus Crotalus

Several types of disintegrins have been isolated from Crotalus spp rattlesnakes, including RGD disintegrins, and PIII-SVMPs. We isolated six cDNAs from snake venom glands using RT-PCR. Three RGD disintegrins (atroxatin, mojastin, and viridistatin) and three PIII-SVMPs (catroriarin, scutiarin, and viristiarin) cDNAs were isolated from the rattlesnakes Crotalus atrox, Crotalus scutulatus scutulatus, and Crotalus viridis viridis, respectively. Atroxatin and Viridistatin shared 90% amino acid identity to each other, and 87% identity to Mojastin. Scutiarin and Viristiarin were identical. All PIII-SVMPs isolated in this study shared the highest amino acid identity with Catrocollastatin. cDNA and protein sequences for RGD disintegrins, one MVD disintegrin, and PIII-SVMPs of the genus Crotalus (present in the NCBI database), were used in phylogenetic analysis. Neighbor-joining analysis of PIII-SVMP and RGD/MVD disintegrin-coding DNA sequences showed that these groups of genes separate into separate clades. A Phi(ST) pairwise comparison and Analysis of Molecular Variance (AMOVA) between PIII-SVMPs and RGD/MVD disintegrins showed significant genetic differences. Mutations observed in ten of the cDNAs analyzed did not affect Cys-coding sequences. Our K(A)/K(S) data suggest that rapid evolution occurred between the genes coding for PIII-SVMPs resulting, in the production of RGD disintegrin-coding genes. However, once these genes diverged, mutations in the PIII-SVMP-coding genes were accumulated less frequently.

A novel peptide from the ACEI/BPP-CNP precursor in the venom of Crotalus durissus collilineatus

In crotaline venoms, angiotensin-converting enzyme inhibitors [ACEIs, also known as bradykinin potentiating peptides (BPPs)], are products of a gene coding for an ACEI/BPP-C-type natriuretic peptide (CNP) precursor. In the genes from Bothrops jararaca and Gloydius blomhoffii, ACEI/BPP sequences are repeated. Sequencing of a cDNA clone from venom glands of Crotalus durissus collilineatus showed that two ACEIs/BPPs are located together at the N-terminus, but without repeats. An additional sequence for CNP was unexpectedly found at the C-terminus. Homologous genes for the ACEI/BPP-CNP precursor suggest that most crotaline venoms contain both ACEIs/BPPs and CNP. The sequence of ACEIs/BPPs is separated from the CNP sequence by a long spacer sequence. Previously, there was no evidence that this spacer actually coded any expressed peptides. Aird and Kaiser (1986, unpublished) previously isolated and sequenced a peptide of 11 residues (TPPAGPDVGPR) from Crotalus viridis viridis venom. In the present study, analysis of the cDNA clone from C. d. collilineatus revealed a nearly identical sequence in the ACEI/BPP-CNP spacer. Fractionation of the crude venom by reverse phase HPLC (C(18)), and analysis of the fractions by mass spectrometry (MS) indicated a component of 1020.5 Da. Amino acid sequencing by MS/MS confirmed that C. d. collilineatus venom contains the peptide TPPAGPDGGPR. Its high proline content and paired proline residues are typical of venom hypotensive peptides, although it lacks the usual N-terminal pyroglutamate. It has no demonstrable hypotensive activity when injected intravenously in rats; however, its occurrence in the venoms of dissimilar species suggests that its presence is not accidental. Evidence suggests that these novel toxins probably activate anaphylatoxin C3a receptors.

Venom Proteomes of Closely RelatedSistrurusRattlesnakes with Divergent Diets

The protein composition of the venoms of the three subspecies of Sistrurus catenatus (S. c. catenatus, tergeminus, and edwardsii) and a basal species, Sistrurus miliarius barbouri, were analyzed by RP-HPLC, N-terminal sequencing, MALDI-TOF peptide mass fingerprinting, and CID-MS/MS. The venoms of the four Sistrurus taxa contain proteins from 11 families. The protein family profile and the relative abundance of each protein group in the different venoms are not conserved. Myotoxins and 2-chain PLA2s were detected only in S.c. catenatus and S.c. tergeminus, whereas C-type BPP and Kunitz-type inhibitors were exclusively found in S.c. edwardsii and Sistrurus miliarius barbouri. Among major protein families, taxa were most similar in their metalloproteases (protein similarity coefficient value: 34%) and most divergent in PLA2s (12%), with values for disintegrins and serine proteases lying between these extremes (25 and 20%, respectively). The patterns of venom diversity points to either a gain in complexity in S. catenatus taxa or a loss of venom diversity occurring early on in the evolution of the group involving the lineage connecting S. milarius to the other taxa. The high degree of differentiation in the venom proteome among recently evolved congeneric taxa emphasizes the uniqueness of the venom composition of even closely related species that have different diets. Comparative proteomic analysis of Sistrurus venoms provides a comprehensive catalog of secreted proteins, which may contribute to a deeper understanding of the biology and ecology of these North American snakes and may also serve as a starting point for studying structure-function correlations of individual toxins.

Evolution of rattlesnakes (Viperidae; Crotalus) in the warm deserts of western North America shaped by Neogene vicariance and Quaternary climate change

During Pleistocene, the Laurentide ice sheet rearranged and diversified biotic distributions in eastern North America, yet had minimal physical impact in western North America where lineage diversification is instead hypothesized to result from climatic changes. If Pleistocene climatic fluctuations impacted desert species, the latter would reflect patterns of restricted gene flow concomitant with indications of demographic bottlenecks. Accordingly, molecular evidence for refugia should be present within these distributions and for subsequent range expansions as conditions improved. We sought answers to these questions by evaluating mitochondrial DNA (mtDNA) sequences from four species of rattlesnakes [Crotalus mitchellii (speckled rattlesnake), Crotalus cerastes (sidewinder), Crotalus tigris (tiger rattlesnake), Crotalus ruber (red diamond rattlesnake)] with distributions restricted to desert regions of southwestern North America. We inferred relationships using parsimony and maximum likelihood, tested intraspecific clades for population expansions, applied an isolation-with-migration model to determine bi-directional migration rates (m) among regions, and inferred divergence times for species and clades by applying a semiparametric penalized likelihood approach to our molecular data. Evidence for significant range expansion was present in two of eight regions in two species (Crotalus mitchellii pyrrhus, C. tigris region north). Two species (C. cerastes, C. mitchellii) showed a distribution concomitant with northward displacement of Baja California from mainland México, followed by vicariant separation into subclades. Effects of Pleistocene climate fluctuations were found in the distributions of all four species. Three regional diversification patterns were identified: (i) shallow genetic diversity that resulted from Pleistocene climatic events (C. tigris, C. ruber); (ii) deep Pleistocene divisions indicating allopatric segregation of subclades within refugia (C. mitchellii, C. cerastes); and (iii) lineage diversifications that extended to Pliocene or Late Miocene (C. mitchellii, C. cerastes). Clade-diversifying and clade-constraining effects impacted the four species of rattlesnakes unequally. We found relatively high levels of molecular diversification in the two most broadly distributed species (C. mitchellii, C. cerastes), and lower levels of genetic diversification in the two species (C. tigris, C. ruber) whose ranges are relatively more restricted. Furthermore, in several cases, the distributions of subspecies were not congruent with our molecular information. We suggest regional conservation perspectives for southwestern deserts cannot rely upon subspecies as biodiversity surrogates, but must instead employ a molecular and deep historical perspective as a primary mechanism to frame biodiversity reserves within this region.

Phylogeographic structure and historical demography of the western diamondback rattlesnake (Crotalus atrox): A perspective on North American desert biogeography

The western diamondback rattlesnake (Crotalus atrox) is a prominent member of North American desert and semi-arid ecosystems, and its importance extends from its impact on the region's ecology and imagery, to its medical relevance as a large deadly venomous snake. We used mtDNA sequences to identify population genetic structure and historical demographic patterns across the range of this species, and relate these to broader patterns of historical biogeography of desert and semi-arid regions of the southwestern USA and adjacent Mexico. We inferred a Late Pliocene divergence between peninsular and continental lineages of Crotalus, followed by an Early Mid Pleistocene divergence across the continental divide within C. atrox. Within desert regions (Sonoran and Chihuahuan Deserts, Southern Plains, and Tamaulipan Plain) we observed population structure indicating isolation of populations in multiple Pleistocene refugia on either side of the continental divide, which we attempt to identify. Evidence of post-glacial population growth and range expansion was inferred, particularly in populations east of the continental divide. We observed clear evidence of (probably recent) gene flow across the continental divide and secondary contact of haplotype lineages. This recent gene flow appears to be particularly strong in the West-to-East direction. Our results also suggest that Crotalus tortugensis (Tortuga Island rattlesnake) and a population of 'C. atrox' inhabiting Santa Cruz Island (in the Gulf of California) previously suggested to be an unnamed species, are in fact deeply phylogenetically nested within continental lineages of C. atrox. Accordingly, we suggest C. tortugensis and 'C. atrox' from Santa Cruz Island be placed in the synonymy of C. atrox.

Mitochondrial DNA sequences from dried snake venom: a DNA barcoding approach to the identification of venom samples

Outdated nomenclature and incorrect taxonomic characterisation of snake venoms in the current toxinological literature have serious implications for the replicability of results from snake venom toxin research. The situation has not improved, despite attempts to supply toxinologists with regular updates on snake systematics. Here, we demonstrate the successful extraction of DNA, and subsequent sequencing of the mitochondrial 12S gene, from dried snake venoms. This approach offers a new and potentially straightforward method for accurate species identification. Mitochondrial DNA (mtDNA) sequences isolated from snake venom can be used to clarify or validate snake species identification through comparison against existing sequences in the GenBank database, and through phylogenetic analyses with other sequences. Pooled venoms can also be screened a priori for the presence of multiple species, and the species names on the labels of commercial venoms verified. Moreover, if the species from which the venom sample has been taken is known, and the specimen is available as a voucher, the mtDNA sequence of the haplotype isolated from that species venom sample could serve as a sequence standard (or 'DNA barcode') for that species. Our new method of DNA barcoding venoms ensures the identification of venoms even after future taxonomic changes.

Tracing an invasion: landbridges, refugia, and the phylogeography of the Neotropical rattlesnake (Serpentes: Viperidae: Crotalus durissus)

Abstract Pleistocene fragmentation of the Amazonian rainforest has been hypothesized to be a major cause of Neotropical speciation and diversity. However, the role and even the reality of Pleistocene forest refugia have attracted much scepticism. In Amazonia, previous phylogeographical studies have focused mostly on organisms found in the forests themselves, and generally found speciation events to have predated the Pleistocene. However, molecular studies of open-formation taxa found both north and south of the Amazonian forests, probably because of vicariance resulting from expansion of the rainforests, may provide novel insights into the age of continuous forest cover across the Amazon basin. Here, we analyse three mitochondrial genes to infer the phylogeography of one such trans-Amazonian vicariant, the Neotropical rattlesnake (Crotalus durissus), which occupies primarily seasonal formations from Mexico to Argentina, but avoids the rainforests of Central and tropical South America. The phylogeographical pattern is consistent with gradual dispersal along the Central American Isthmus, followed by more rapid dispersal into and across South America after the uplift of the Isthmus of Panama. Low sequence divergence between populations from north and south of the Amazon rainforest is consistent with mid-Pleistocene divergence, approximately 1.1 million years ago (Ma). This suggests that the Amazonian rainforests must have become fragmented or at least shrunk considerably during that period, lending support to the Pleistocene refugia theory as an important cause of distribution patterns, if not necessarily speciation, in Amazonian forest organisms. These results highlight the potential of nonforest species to contribute to an understanding of the history of the Amazonian rainforests themselves.

Identification of crotasin, a crotamine-related gene of Crotalus durissus terrificus

Crotamine is a cationic peptide (4.9 kDa, pI 9.5) of South American rattlesnake, Crotalus durissus terrificus' venom. Its presence varies according to the subspecies or the geographical locality of a given species. At the genomic level, we observed the presence of 1.8 kb gene, Crt-p1, in crotamine-positive specimens and its absence in crotamine-negative ones. In this work, we described a crotamine-related 2.5 kb gene, crotasin (Cts-p2), isolated from crotamine-negative specimens. Reverse transcription coupled to polymerase chain reaction indicates that Cts-p2 is abundantly expressed in several snake tissues, but scarcely expressed in the venom gland. The genome of crotamine-positive specimen contains both Crt-p1 and Cts-p2 genes. The present data suggest that both crotamine and crotasin have evolved by duplication of a common ancestor gene, and the conservation of their three disulfide bonds indicates that they might adopt the same fold as beta-defensin. The physiological function of the crotasin is not yet known.

Direct organogenesis of Mandevilla illustris (Vell) Woodson and effects of its aqueous extract on the enzymatic and toxic activities of Crotalus durissus terrificus snake venom

In order to produce explants of Mandevilla illustris (Vell) Woodson for the "Cerrado in vitro", the Germplasm Bank of UNAERP, we carried out a micropropagation protocol using MS or MS/3 medium supplemented with different concentrations of 6-benzyladeninepurine (BA), Zeatin or 2-isopentenyladenine for nodal segment growth, and alpha-naphthaleneacetic acid, indole-3-butyric acid (IBA) or 1,4 dithiothreitol for rooting. For nodal segments, all the cytokinins tested yielded similar results. However, 2.22 micro M BA is more economical to use. MS/3 medium supplemented with 0.49 micro M IBA was the most appropriate medium for rooting, resulting in 29% rooted explants. The crude aqueous extract from the subterranean system (SS) of M. illustris was assayed for its inhibitory action on the enzymatic activity of Crotalus durissus terrificus snake venom, isolated basic phospholipase A2 (CB) and crotoxin. It totally inhibited the phospholipase activity of crude Cdt venom and CB toxin and inhibited the phospholipase activity of crotoxin by 49%. The toxic action of both the crude venom and crotoxin was partially inhibited-there was a prolonged survival time and a 40.0% decrease in lethality.

Geographic variations, cloning, and functional analyses of the venom acidic phospholipases A2 of Crotalus viridis viridis

Geographic venom samples of Crotalus viridis viridis were obtained from South Dakota, Wyoming, Colorado, Oklahoma, Texas, New Mexico, and Arizona. From these samples, the phospholipases A(2) (PLA(2)s) were purified and their N-terminal sequences, precise masses, and in vitro enzymatic activities were determined. We purified two to four distinct acidic PLA(2)s from each sample; some of them displayed different inhibition specificities toward mammalian platelets. One of the acidic PLA(2)s induced edema, but had no anti-platelet activity. There was also a common basic PLA(2) myotoxin in all the samples. We have cloned five acidic PLA(2)s and several hybrid-like nonexpressing PLA(2)s. Molecular masses and N-terminal sequences of the purified PLA(2)s were matched with those deduced from the cDNA sequences, and the complete amino acid sequences of five novel acidic PLA(2)s were thus solved. They share 78% or greater sequence identity, and a cladogram based on the sequences of many venom acidic PLA(2)s of New World pit vipers revealed at least two subtypes. The results contribute to a better understanding of the ecogenetic adaptation of rattlesnakes and the structure-activity relationships and evolution of the acidic PLA(2)s in pit viper venom.

Proteome analysis of the thermoreceptive pit membrane of the western diamondback rattlesnake Crotalus atrox

Rattlesnakes detect their prey's temperature by means of a cavern-like structure, the pit organ. The sensory component of this organ lies within a thin membrane called the pit membrane. Proteome analysis conducted on this neurosensory tissue revealed only a relatively small number of proteins, thereby depicting its high degree of specialization. In addition to containing blood serum and structural proteins, the proteome of this membrane appears to be strikingly similar to that of isolated rattlesnake brain mitochondria. Indeed, our results show that over 80% of the detected tissue proteins are of mitochondrial origin. Fluorescence microscopy studies of these organelles indicate their dense arrangement and accumulation in structures which have been previously reported to be the terminal ends of free nerve fibers of the innervating trigeminal branches. Thus, original ultrastructural observations are paralleled by our findings at the molecular level.

Molecular systematics of the western rattlesnake, Crotalus viridis (Viperidae), with comments on the utility of the D-loop in phylogenetic studies of snakes

Crotalus viridis, the western rattlesnake, ranges throughout western North America and has been divided into at least eight subspecies. However, the validity of and relationships among these subspecies and the monophyly of C. viridis as a whole are questionable. We used mitochondrial DNA sequence data from the D-loop region and ND2 gene to examine the relationships among 26 populations of C. viridis and to test the monophyly of this species. These data were analyzed separately and combined using maximum-likelihood and maximum-parsimony. The C. viridis group was monophyletic in all combined analyses, consisting of two strongly divergent clades. We recommend that these clades be recognized as two distinct evolutionary species: C. viridis and C. oreganus. Crotalus viridis should be restricted to the subspecies viridis and nuntius and the remaining subspecies be assigned to the species C. oreganus. Our data do not allow strong evaluation of the validity of the subspecies. We found that the ND2 gene had greater sequence divergences among closely related individuals than the D-loop region, but this relationship reversed at higher levels of divergence. This pattern is apparently due to: (1) ND2 third positions evolving faster than the D-loop but becoming saturated at higher levels of divergence, and (2) the D-loop evolving faster than ND2 second (and possibly first) positions. Our results suggest that the ND2 gene is preferable for examining intraspecific relationships and the D-loop may better resolve relationships between species of snakes. The latter result is contrary to the common perception of the phylogenetic utility of the D-loop. Another unusual result is that the 145 bp spacer region, adjacent to the 5' end of the light strand of the D-loop, provides greater phylogenetic resolution than the 1030 bp D-loop.

Insights into population ecology and sexual selection in snakes through the application of DNA-based genetic markers

Hypervariable genetic markers have revolutionized studies of kinship, behavioral ecology, and population biology in vertebrate groups such as birds, but their use in snakes remains limited. To illustrate the value of such markers in snakes, we review studies that have used microsatellite DNA loci to analyze local population differentiation and parentage in snakes. Four ecologically distinct species of snakes all show evidence for differentiation at small spatial scales (2-15 km), but with substantial differences among species. This result highlights how genetic analysis can reveal hidden aspects of the natural history of difficult-to-observe taxa, and it raises important questions about the ecological factors that may contribute to restricted gene flow. A 3-year study of genetic parentage in marked populations of the northern water snake showed that (1) participation in mating aggregations was a poor predictor of genetic-based measures of reproductive success; (2) multiple paternity was high, yet there was no detectable fitness advantage to multiple mating by females; and (3) the opportunity for selection was far higher in males than in females due to a larger variance in male reproductive success, and yet this resulted in no detectable selection on morphological variation in males. Thus genetic markers have provided accurate measures of individual reproductive success in this species, an important step toward resolving the adaptive significance of key features including multiple paternity and reversed sexual size dimorphism. Overall these studies illustrate how genetic analyses of snakes provide previously unobtainable information of long-standing interest to behavioral ecologists.

[Phylogenetic relationships among Viperidae, Crotalinae based on mitochondrial 12S rRNA sequence variations]

This paper analyzed the phylogenetic relationships and classification of pit vipers (Viperidae, Crotalinae) based on mitochondrial 12S rRNA gene sequence variations. We have sequenced mtDNA 12S rRNA gene about 370 bp fragment from Gloydius saxatilis Emelianov, Gloydius shedaoensis Zhao, Gloydius ussurriensis (Emelianov), Trimeresurus stejnegeri Schmidt and Deinagkistrodon acutus, Gloydius brevicaudus (Stejneger) from two different localities, respectively. Combined with the sequence of Dinodon semicarinatus from GenBank, we have constructed two molecular phylogenetic trees using both Maximum-parsimony analysis and Neighbor-joining. Our results support the following conclusions: (1) Deinagkistrodon is a valid genus; (2) we also discussed the origin of Gloydius shedaoensis and showed it is a subspecies of Gloydius saxatilis; (3) Gloydius brevicaudsu and Gloydius ussurriensis are classified two species.

A comparison of RAPD versus microsatellite DNA markers in population studies of the massasauga rattlesnake

We compared genetic differentiation among populations of the threatened massasauga rattlesnake (Sistrurus c. catenatus) using two types of nuclear molecular markers: randomly amplified polymorphic DNA (RAPD) markers and microsatellites. Analyses of molecular variance (AMOVA) and G(ST) and F(ST) analyses indicated that levels of among-population differentiation between regional populations (>100 km) were comparable for both markers. However, microsatellites were superior in population assignment tests and at discerning fine-scale genetic differentiation between subpopulations separated by tens of kilometers. These results argue that both types of markers are suitable for defining broad-scale genetic structures in snake populations and can provide important inputs into conservation initiatives of focal taxa. However, our analyses suggest that microsatellites 3re better for detecting structure at limited spatial scales.

Historical biogeography of the Western Rattlesnake (Serpentes: viperidae: Crotalus viridis), inferred from mitochondrial DNA sequence information

We infer the phylogeography of the Western Rattlesnake (Crotalus viridis) using phylogenetic analysis of mitochondrial DNA sequences from 1345 bp of the genes for cytochrome b and NADH dehydrogenase subunit 4. Two main clades are revealed: one includes populations from east and south of the Rocky Mountains (conventionally referred to as Crotalus viridis viridis and C. v. nuntius), and the other consists of populations west of the Rocky Mountains. Within the western clade, a population from southern Arizona (C. v. cerberus) represents the sister taxon to the remaining western populations. The conventional subspecies recognized in this species do not fully correspond to the phylogenetic pattern, and a review of the systematic status of several populations is needed. Our data allow the inferences that small body size evolved twice and that the ability of one population (C. v. concolor) to secrete highly lethal toxins related to Mojave toxin arose within the complex. Our phylogeny should represent the basis for further studies on the causes of geographical variation in this complex.

Molecular cloning and sequence analysis of aggretin, a collagen-like platelet aggregation inducer

A cDNA library derived from the Malayan-pit-viper (Calloselasma rhodostoma) venom gland was constructed in the phagemid vector. Using the information of the N-terminal amino acid sequences of two subunits of aggretin, synthetic mixed-base oligonucleotides were employed as a screening probe for colony hybridization. Separate cDNA clones encoding for the alpha and beta chains of aggretin were isolated and sequenced. The results revealed that mature alpha and beta chains contain 136 and 123 amino acid residues, respectively. Aggretin subunits show high degrees of identity with respective subunits (50-60% for alpha, 49-58% for beta) of C-type lectin-like snake venoms. The identity to rattlesnake lectin is relatively lower (i.e., 39 and 30%). All cysteine residues in each chain of aggretin are well conserved and located at the positions corresponding to those of C-type lectins. Thus, three intracatenary disulfide bridges and an interchain disulfide bond between Cys83(alpha) and Cys75(beta) may be allocated. This is the first report regarding the entire sequence of venom GPIa/IIa agonist. According to the alignment of amino acid sequences, hypervariable regions among these C-type lectin-like proteins were revealed. These hypervariable regions are proposed to be the counterparts directly interacting with different receptors or different domains of a receptor on the surface of platelet.

Positive Darwinian selection in Vipera palaestinae phospholipase A2 genes is unexpectedly limited to the third exon

The venom of Vipera palaestinae contains a two-component toxin, consisting of an acidic phospholipase A2 (PLA2) and a basic protein. Here we report the cloning and sequence analysis of the complete V. palaestinae PLA2 genes. Since in all Viperidae PLA2 multigene families the 5' and 3' flanking regions are highly conserved, we designed oligonucleotide primers that allow amplification of the whole PLA2 multigene family in a single step. The structural organization of both genes is the same as in the Vipera ammodytes PLA2 multigene family, there being five exons separated by four introns. Comparison of V. palaestinae PLA2 genes with other Viperidae PLA2 genes has shown that the structural organization of the genes and the nucleotide sequence of all introns and flanking regions are highly conserved, whereas the third exon clearly shows a higher number of amino acid replacements, an indication of positive Darwinian selection. The positive Darwinian selection is surprisingly limited to the third exon, in contrast to other Viperidae PLA2 genes, where it is present in all mature protein coding exons.

Genetic structure of populations of the threatened eastern massasauga rattlesnake, Sistrurus c. catenatus: evidence from microsatellite DNA markers

Throughout its distribution in North America, the threatened eastern massasauga rattlesnake (Sistrurus c. catenatus) persists in a series of habitat-isolated disjunct populations of varying size. Here, we use six microsatellite DNA loci to generate information on the degree of genetic differentiation between, and the levels of inbreeding within populations to understand how evolutionary processes operate in these populations and aid the development of conservation plans for this species. Samples were collected from 199 individuals from five populations in Ontario, New York and Ohio. Our results show that all sampled populations: (i) differ significantly in allele frequencies even though some populations are < 50 km apart, and may contain genetically distinct subpopulations < 2 km apart; (ii) have an average of 23% of alleles that are population specific; and (iii) have significant FIS values (mean overall FIS = 0.194) probably due to a combination of Wahlund effects resulting from fine-scale genetic differentiation within populations and the presence of null alleles. Our results imply that massasauga populations may be genetically structured on an extremely fine scale even within continuous populations, possibly due to limited dispersal. Additional information is needed to determine if dispersal and mating behaviour within populations can account for this structure and whether the observed differentiation is due to random processes such as drift or to local adaptation. From a conservation perspective, our results imply that these massasauga populations should be managed as demographically independent units and that each has high conservation value in terms of containing unique genetic variation.

Cloning and characterization of a repetitive 1.9 Kb HindIII DNA fragment from Crotalus durissus terrificus genome

Eukaryote genomes are endowed with varying quantities of repeated DNA families. These families show different patterns of conservation among species, copy numbers, chromosomal distribution, and transcription. Characterization of repeated DNA sequences could help to understand the genome anatomy and organization or be used in molecular systematics and molecular evolution studies. We describe here a repetitive DNA sequence of the HindIII family present in the genome of the rattlesnake Crotalus durissus terrificus. In Brazil, the family Crotalus is comprised only by one species durissus, which include several subspecies. The number and distribution of these subspecies are controversial. In the present study, the genomic DNA of a female rattlesnake was digested with HindIII resulting in a strong 1.9 Kb band. A partial genomic library was constructed from the 1.9 Kb DNAs rescued from the agarose gel after HindIII digestion and ligated to the vector pGEM3Zf(+) (Promega). Analysis of 69 clones, 44 hybridized with the 1.9 Kb probe isolated from one of the clones-clone 76, indicating that the DNA isolated from this clone should represent the 1.9 Kb HindIII fragment. This 1.9 Kb HindIII DNA was completely characterized by sequencing.

A phospholipase A2-like pseudogene retaining the highly conserved introns of Mojave toxin and other snake venom group II PLA2s, but having different exons

Mojave toxin is a neurotoxic, heterodimeric phospholipase A2 (PLA2) from the venom of the Mojave rattlesnake (Crotalus scutulatus scutulatus) and is characteristic of all rattlesnake presynaptic neurotoxins. Here, we describe a phospholipase A2 pseudogene (psi-Mtx) located 2,000 nucleotides upstream, and on the opposite DNA strand, from a gene for Mojave toxin acidic subunit (Mtx-a). The pseudogene lacks the first exon and a few segments of noncoding DNA found in functional snake venom PLA2 genes, but does have the coding information for a complete PLA2 protein. psi-Mtx retains the unusual gene sequence similarity pattern found in functional viperid PLA2 genes. When compared to genes from C. s. scutulatus and the Hahn snake (Trimeresurus flavoviridus), psi-Mtx shows strong conservation of nocoding regions and variable protein-coding regions. Although the nocoding regions of psi-Mtx are conserved with respect to other viperid PLA2 genes, the three exons code for a unique PLA2-like protein similar in sequence to ammodytoxin b found in the venom of the western sand viper (Vipera ammodytes ammodytes). The structure of these genes suggests a common ancestor for all viperid PLA2 genes. Phylogenetic analysis of psi-Mtx, Mtx-a, Mtx-b, pgPLA 1a, and pgPLA 1b suggest that psi-Mtx diverged from an ancestral sequence before the presumed gene duplication event leading to Mtx-a and Mtx-b. However, analysis of the basis of coding regions alone gives a conflicting result.

Isolation and characterization of microsatellite loci for use in population genetic analysis in the timber rattlesnake, Crotalus horridus

A Crotalus horridus genomic library was screened for clones containing microsatellite loci by hybridization with oligonucleotides consisting of a (dC x dA)n dinucleotide repeat. Primers designed to amplify six of the microsatellite loci were used to screen 32 unrelated individuals representing populations in eastern Pennsylvania, southern New Jersey, North Carolina, South Carolina, and Alabama. The six microsatellite loci were all polymorphic, with two to nine alleles, and heterozygote frequencies at each locus from 0.1 to 0.69. Allelic frequencies varied among geographically separated populations. Screening of two families produced no evidence of multiple paternity. These microsatellite markers should be useful for the assessment of kinship relationships and genetic diversity within and between populations of C. horridus. The application of this technology could provide a valuable tool for the development of effective conservation and management programs for threatened and endangered populations.

Trimeresurus flavoviridis venom gland phospholipase A2 isozymes genes have evolved via accelerated substitutions

As a step towards understanding the structure and function of phospholipases A2 (PLA2s), five cDNAs encoding Trimeresurus flavoviridis venom gland PLA2 isozymes have been sequenced. They revealed that the 5' and 3' untranslated regions are much more homologous than the protein-coding regions and that base substitutions have occurred at similar rates for the three positions of codons in the protein-coding regions. Such novel findings are of great interest from the viewpoint of molecular evolution. To gain a further insight into this evolutional phenomenon, we have isolated and sequenced six T. flavoviridis PLA2 isozyme genes. Each gene consisted of four exons and three introns and encoded protein of 138 amino-acid residues, including the signal sequence of 16 amino-acid residues. The introns were much more homologous than the protein-coding regions of exons except for the signal peptide-coding region of the first exon. The absence of apparent functional role in the introns suggested that the protein-coding regions, except for the signal peptide-coding domains, have evolved at greater substitution rates than introns. The fact that the numbers of nucleotide substitutions per non-synonymous site are close to or larger than the numbers of nucleotide substitutions per synonymous site for relevant pairs of genes revealed that Darwinian-type accelerated substitutions have occurred in the protein-coding regions of exons. This is compatible with the presence of PLA2 species with diverse physiological activities in the venom.