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

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.

Phylogeography of the smooth greensnake, Opheodrys vernalis (Squamata: Colubridae): divergent lineages and variable demographics in a widely distributed yet enigmatic species

Phylogeographic studies can uncover robust details about the population structure, demographics, and diversity of species. The smooth greensnake, Opheodrys vernalis, is a small, cryptic snake occupying mesic grassland and sparsely wooded habitats. Although O. vernalis has a wide geographical range, many metapopulations are patchy and some are declining. We used mitochondrial DNA and double digest restriction-site associated DNA sequencing to construct the first phylogeographic assessment of O. vernalis. Genomic analysis of 119 individuals (mitochondrial DNA) and a subset of another 45 smooth greensnakes (nuclear DNA; N = 3031 single nucleotide polymorphisms) strongly supports two longitudinally separated lineages, with admixture in the Great Lakes region. Post-Pleistocene secondary contact best explains admixture from populations advancing northwards. Overall, populations expressed low heterozygosity, variable inbreeding rates, and moderate to high differentiation. Disjunct populations in the Rocky Mountains and central Great Plains regions might be contracting relicts, whereas northerly populations in more continuous mesic habitats (e.g., Prairie Pothole region, southern Canada) had signals of population expansion. Broadly, conservation management efforts should be focused on local populations, because habitat connectivity may facilitate gene flow and genetic diversity.

Genome-wide data reveal extensive gene flow during the diversification of the western rattlesnakes (Viperidae: Crotalinae: Crotalus)

Hybridization and introgression are important, but often overlooked processes when inferring phylogenies. When these processes are not accounted for and a strictly diverging phylogenetic model is applied to groups with a history of hybridization, phylogenetic inference and parameter estimation can be inaccurate. Recent developments in phylogenetic network approaches coupled with the increasing availability of genomic data allow inferences of reticulate evolutionary histories across the tree of life. The western rattlesnake species group (C. viridis species complex, C. mitchellii species complex, C. scutulutas, and C. tigris) is an iconic snake lineage that is widespread across western North America. This group is composed of several species complexes with unclear species limits, likely the result of ongoing gene flow among nascent lineages. Here I generate reduced representation genomic data and test for a history of reticulation within this group. I demonstrate that all species have undergone hybridization with at least one other lineage, suggesting introgression is widespread in this group. Topologies differ between phylogenies estimated under the multispecies coalescent and multispecies network coalescent methods, indicating that gene flow has obscured phylogenetic relationships within this group. These past introgression events are predominantly restricted to species that co-occur geographically. However, within species that have a history of introgression, this signature is detected regardless of specimen sampling across geography. Overall, my results suggest the accumulation of reproductive isolating barriers occurs slowly in rattlesnakes which likely leads to the difficulty in delimiting species, furthermore, the results of this study have implications for trait evolution in this group.

Biological and Biochemical Characterization of Coronado Island Rattlesnake (Crotalus helleri caliginis) Venom and Antivenom Neutralization

The Baja California Peninsula has over 250 islands and islets with many endemic species. Among them, rattlesnakes are the most numerous but also one of the least studied groups. The study of island rattlesnake venom could guide us to a better understanding of evolutionary processes and the description of novel toxins. Crotalus helleri caliginis venom samples were analyzed to determine possible ontogenetic variation with SDS-PAGE in one and two dimensions and with RP-HPLC. Western Blot, ELISA, and amino-terminal sequencing were used to determine the main components of the venom. The biological and biochemical activities demonstrate the similarity of C. helleri caliginis venom to the continental species C. helleri helleri, with both having low proteolytic and phospholipase A₂ (PLA₂) activity but differing due to the absence of neurotoxin (crotoxin-like) in the insular species. The main components of the snake venom were metalloproteases, serine proteases, and crotamine, which was the most abundant toxin group (30–35% of full venom). The crotamine was isolated using size-exclusion chromatography where its functional effects were tested on mouse phrenic nerve–hemidiaphragm preparations in which a significant reduction in muscle twitch contractions were observed. The two Mexican antivenoms could neutralize the lethality of C. helleri caliginis venom but not the crotamine effects.

A symphony of destruction: Dynamic differential fibrinogenolytic toxicity by rattlesnake (Crotalus and Sistrurus) venoms

What factors influence the evolution of a heavily selected functional trait in a diverse clade? This study adopts rattlesnakes as a model group to investigate the evolutionary history of venom coagulotoxicity in the wider context of phylogenetics, natural history, and biology. Venom-induced clotting of human plasma and fibrinogen was determined and mapped onto the rattlesnake phylogenetic tree to reconstruct the evolution of coagulotoxicity across the group. Our results indicate that venom phenotype is often independent of phylogenetic relationships in rattlesnakes, suggesting the importance of diet and/or other environmental variables in driving venom evolution. Moreover, the striking inter- and intraspecific variability in venom activity on human blood highlights the considerable variability faced by physicians treating envenomation. This study is the most comprehensive effort to date to describe and characterize the evolutionary and biological aspects of coagulotoxins in rattlesnake venom. Further research at finer taxonomic levels is recommended to elucidate patterns of variation within species and lineages.

Thermal ecology and baseline energetic requirements of a large-bodied ectotherm suggest resilience to climate change

Most studies on how rising temperatures will impact terrestrial ectotherms have focused on single populations or multiple sympatric species. Addressing the thermal and energetic implications of climatic variation on multiple allopatric populations of a species will help us better understand how a species may be impacted by altered climates.We used eight years of thermal and behavioral data collected from four populations of Pacific rattlesnakes (Crotalus oreganus) living in climatically distinct habitat types (inland and coastal) to determine the field-active and laboratory-preferred body temperatures, thermoregulatory metrics, and maintenance energetic requirements of snakes from each population.Physical models showed that thermal quality was best at coastal sites, but inland snakes thermoregulated more accurately despite being in more thermally constrained environments. Projected increases of 1 and 2°C in ambient temperature result in an increase in overall thermal quality at both coastal and inland sites.Population differences in modeled standard metabolic rate estimates were driven by body size and not field-active body temperature, with inland snakes requiring 1.6× more food annually than coastal snakes.All snakes thermoregulated with high accuracy, suggesting that small increases in ambient temperature are unlikely to impact the maintenance energetic requirements of individual snakes and that some species of large-bodied reptiles may be robust to modest thermal perturbations under conservative climate change predictions. ​.

A Clot Twist: Extreme Variation in Coagulotoxicity Mechanisms in Mexican Neotropical Rattlesnake Venoms

Rattlesnakes are a diverse clade of pit vipers (snake family Viperidae, subfamily Crotalinae) that consists of numerous medically significant species. We used validated in vitro assays measuring venom-induced clotting time and strength of any clots formed in human plasma and fibrinogen to assess the coagulotoxic activity of the four medically relevant Mexican rattlesnake species Crotalus culminatus, C. mictlantecuhtli, C. molossus, and C. tzabcan. We report the first evidence of true procoagulant activity by Neotropical rattlesnake venom in Crotalus culminatus. This species presented a strong ontogenetic coagulotoxicity dichotomy: neonates were strongly procoagulant via Factor X activation, whereas adults were pseudo-procoagulant in that they converted fibrinogen into weak, unstable fibrin clots that rapidly broke down, thereby likely contributing to net anticoagulation through fibrinogen depletion. The other species did not activate clotting factors or display an ontogenetic dichotomy, but depleted fibrinogen levels by cleaving fibrinogen either in a destructive (non-clotting) manner or via a pseudo-procoagulant mechanism. We also assessed the neutralization of these venoms by available antivenom and enzyme-inhibitors to provide knowledge for the design of evidence-based treatment strategies for envenomated patients. One of the most frequently used Mexican antivenoms (Bioclon Antivipmyn®) failed to neutralize the potent procoagulant toxic action of neonate C. culminatus venom, highlighting limitations in snakebite treatment for this species. However, the metalloprotease inhibitor Prinomastat substantially thwarted the procoagulant venom activity, while 2,3-dimercapto-1-propanesulfonic acid (DMPS) was much less effective. These results confirm that venom-induced Factor X activation (a procoagulant action) is driven by metalloproteases, while also suggesting Prinomastat as a more promising potential adjunct treatment than DMPS for this species (with the caveat that in vivo studies are necessary to confirm this potential clinical use). Conversely, the serine protease inhibitor 4-(2-aminoethyl)benzenesulfonyl fluoride hydrochloride (AEBSF) inhibited the direct fibrinogen cleaving actions of C. mictlantecuhtli venom, thereby revealing that the pseudo-procoagulant action is driven by kallikrein-type serine proteases. Thus, this differential ontogenetic variation in coagulotoxicity patterns poses intriguing questions. Our results underscore the need for further research into Mexican rattlesnake venom activity, and also highlights potential limitations of current antivenom treatments.

The scales of coevolution: comparative phylogeography and genetic demography of a locally adapted venomous predator and its prey

Coevolutionary theory predicts that differences in the genetic demography of interacting species can influence patterns of local adaptation by affecting the potential of local populations to respond to selection. We conducted a comparative phylogeographical study of venomous rattlesnakes and their venom-resistant ground squirrel prey across California, and assessed how effective population size (Ne) estimates correspond with a previously documented pattern of rattlesnake local adaptation. Using RAD sequencing markers, we detected lineage relationships among both the rattlesnakes (Crotalus oreganus ssp.) and ground squirrels (Otospermophilus sp.) that are incongruent with previous phylogenetic hypotheses. Both rattlesnakes and squirrels share a deep divergence at the Sacramento–San Joaquin River Delta. At this broad phylogeographical scale, we found that the locally adapted rattlesnakes had higher Ne than squirrels. At the population scale, snakes also had larger Ne accompanied by larger values of several metrics of population genetic diversity. However, the specific magnitude of local adaptation of venom activity to ground squirrel venom resistance was not significantly correlated with local differences in Ne or other diversity statistics between predator and prey populations, suggesting that other factors in the geographic mosaic of coevolution contribute to the specific local-scale outcomes of this interaction. These results suggest an evolutionary mechanism that may explain some (but clearly not all) of rattlesnake local adaptation in this coevolutionary interaction – larger population sizes raise the adaptive potential of rattlesnakes compared to ground squirrels.

Rattlesnake migrations and the implications of thermal landscapes

BACKGROUND

The importance of thermal resources to terrestrial ectotherms has been well documented but less often considered in larger-scale analyses of habitat use and selection, such as those routinely conducted using standard habitat features such as vegetation and physical structure. Selection of habitat based on thermal attributes may be of particular importance for ectothermic species, especially in colder climates. In Canada, Western Rattlesnakes (Crotalus oreganus) reach their northern limits, with limited time to conduct annual migratory movements between hibernacula and summer habitat. We radio-tracked 35 male snakes departing from 10 different hibernacula. We examined coarse-scale differences in migratory movements across the region, and then compared the route of each snake with thermal landscapes and ruggedness GIS maps generated for different periods of the animals' active season.

RESULTS

We observed dichotomous habitat use (grasslands versus upland forests) throughout most of the species' northern range, reflected in different migratory movements of male snakes emanating from different hibernacula. Snakes utilizing higher-elevation forests moved further during the course of their annual migrations, and these snakes were more likely to use warmer areas of the landscape.

CONCLUSION

In addition to thermal benefits, advantages gained from selective migratory patterns may include prey availability and outbreeding. Testing these alternative hypotheses was beyond the scope of this study, and to collect the data to do so will require overcoming certain challenges. Still, insight into migratory differences between rattlesnake populations and the causal mechanism(s) of migrations will improve our ability to assess the implications of landscape change, management, and efficacy of conservation planning. Our findings suggest that such assessments may need to be tailored to individual dens and the migration strategies of their inhabitants. Additionally, local and landscape-scale migration patterns, as detected in this study, will have repercussions for snakes under climate-induced shifts in ecosystem boundaries and thermal regimes.

Allopatric divergence and secondary contact with gene flow: a recurring theme in rattlesnake speciation

The study of recently diverged lineages whose geographical ranges come into contact can provide insight into the early stages of speciation and the potential roles of reproductive isolation in generating and maintaining species. Such insight can also be important for understanding the strategies and challenges for delimiting species within recently diverged species complexes. Here, we use mitochondrial and nuclear genetic data to study population structure, gene flow and demographic history across a geographically widespread rattlesnake clade, the western rattlesnake species complex (Crotalus cerberus, Crotalus viridis, Crotalus oreganus and relatives), which contains multiple lineages with ranges that overlap geographically or contact one another. We find evidence that the evolutionary history of this group does not conform to a bifurcating tree model and that pervasive gene flow has broadly influenced patterns of present-day genetic diversity. Our results suggest that lineage diversity has been shaped largely by drift and divergent selection in isolation, followed by secondary contact, in which reproductive isolating mechanisms appear weak and insufficient to prevent introgression, even between anciently diverged lineages. The complexity of divergence and secondary contact with gene flow among lineages also provides new context for why delimiting species within this complex has been difficult and contentious historically.

DNA barcodes from snake venom: a broadly applicable method for extraction of DNA from snake venoms

DNA barcoding is a simple technique used to develop a large-scale system of classification that is broadly applicable across a wide variety of taxa. DNA-based analysis of snake venoms can provide a system of classification independent of currently accepted taxonomic relationships by generating DNA barcodes specific to each venom sample. DNA purification from dried snake venoms has previously required large amounts of starting material, has resulted in low yields and inconsistent amplification, and was possible with front-fanged snakes only. Here, we present a modified DNA extraction protocol applied to venoms of both front- and rear-fanged snakes that requires significantly less starting material (1 mg) and yields sufficient amounts of DNA for successful PCR amplification of regions commonly used for DNA barcoding. [Formula: see text].

Geographical features are the predominant driver of molecular diversification in widely distributed North American whipsnakes

Allopatric divergence following the formation of geographical features has been implicated as a major driver of evolutionary diversification. Widespread species complexes provide opportunities to examine allopatric divergence across varying degrees of isolation in both time and space. In North America, several geographical features may play such a role in diversification, including the Mississippi River, Pecos River, Rocky Mountains, Cochise Filter Barrier, Gulf of California and Isthmus of Tehuantepec. We used thousands of nuclear single nucleotide polymorphisms (SNPs) and mitochondrial DNA from several species of whipsnakes (genera Masticophis and Coluber) distributed across North and Central America to investigate the role that these geographical features have played on lineage divergence. We hypothesize that these features restrict gene flow and separate whipsnakes into diagnosable genomic clusters. We performed genomic clustering and phylogenetic reconstructions at the species and population levels using Bayesian and likelihood analyses and quantified migration levels across geographical features to assess the degree of genetic isolation due to allopatry. Our analyses suggest that (i) major genetic divisions are often consistent with isolation by geographical features, (ii) migration rates between clusters are asymmetrical across major geographical features, and (iii) areas that receive proportionally more migrants possess higher levels of genetic diversity. Collectively, our findings suggest that multiple features of the North American landscape contributed to allopatric divergence in this widely distributed snake group.

Phylogenetic analysis of the king cobra, Ophiophagus hannah in Thailand based on mitochondrial DNA sequences

Background: Snakes possess adaptive characteristics of morphology that may result in incorrect reconstruction of phylogeny. Molecular approaches have become the major source of new information for advancing our understanding of evolutionary, genetic relationships, and species identification.

Objective: We studied DNA sequences of Ophiophagus hannah in different parts of Thailand and compared them with those of O. hannah from other countries.

Materials and Methods: We studied the mitochondrial NADH dehydrogenase subunit 2 (ND2) and control region (CR) sequences of 12 individuals O. hannah from different localities across Thailand. Phylogenetic analysis of DNA sequences were compared to the published sequences of O. hannah deposited in NCBI GenBank database from other countries.

Results: O. hannah could be categorized into 2 Clades, 5 haplotypes, and 4 localities based on 43 different nucleotide positions from the 736 bp of ND2 and 673 bp of CR. Clade A was haplotype A from southern Thailand. Clade B consisted of haplotypes B, C, D, and E. Haplotype B and C came from northern Thailand. Haplotype D came from western Thailand, while haplotype E was O. hannah from central Thailand. The DNA sequences of Clade B were similar to the sequences of O. hannah in Myanmar and southern China that are already deposited in NCBI GenBank database.

Conclusion: We found a different genotype of O. hannah from southern Thailand and suggest that this may be a new species of O. hannah.

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.

Incipient speciation with biased gene flow between two lineages of the Western Diamondback Rattlesnake (Crotalus atrox)

We used mitochondrial DNA sequence data from 151 individuals to estimate population genetic structure across the range of the Western Diamondback Rattlesnake (Crotalus atrox), a widely distributed North American pit viper. We also tested hypotheses of population structure using double-digest restriction site associated DNA (ddRADseq) data, incorporating thousands of nuclear genome-wide SNPs from 42 individuals. We found strong mitochondrial support for a deep divergence between eastern and western C. atrox populations, and subsequent intermixing of these populations in the Inter-Pecos region of the United States and Mexico. Our nuclear RADseq data also identify these two distinct lineages of C. atrox, and provide evidence for nuclear admixture of eastern and western alleles across a broad geographic region. We identified contrasting patterns of mitochondrial and nuclear genetic variation across this genetic fusion zone that indicate partially restricted patterns of gene flow, which may be due to either pre- or post-zygotic isolating mechanisms. The failure of these two lineages to maintain complete genetic isolation, and evidence for partially-restricted gene flow, imply that these lineages were in the early stages of speciation prior to secondary contact.

A novel phospholipase A2 (D49) from the venom of the Crotalus oreganus abyssus (North American Grand canyon rattlesnake)

Currently, Crotalus viridis was divided into two species: Crotalus viridis and Crotalus oreganus. The current classification divides "the old" Crotalus viridis into two new and independent species: Crotalus viridis (subspecies: viridis and nuntius) and Crotalus oreganus (subspecies: abyssus, lutosus, concolor, oreganus, helleri, cerberus, and caliginis). The analysis of a product from cDNA (E6d), derived from the gland of a specie Crotalus viridis viridis, was found to produce an acid phospholipase A2. In this study we isolated and characterized a PLA2 (D49) from Crotalus oreganus abyssus venom. Our studies show that the PLA2 produced from the cDNA of Crotalus viridis viridis (named E6d) is exactly the same PLA2 primary sequence of amino acids isolated from the venom of Crotalus oreganus abyssus. Thus, the PLA2 from E6d cDNA is actually the same PLA2 presented in the venom of Crotalus oreganus abyssus and does not correspond to the venom from Crotalus viridis viridis. These facts highlight the importance of performing more studies on subspecies of Crotalus oreganus and Crotalus viridis, since the old classification may have led to mixed results or mistaken data.

Population genetic structure and geographic differentiation in the hot spring snake Thermophis baileyi (Serpentes, Colubridae): indications for glacial refuges in southern-central Tibet

The endangered hot spring snake (Thermophis baileyi) is endemic to the Tibetan Plateau and restricted to a few sites at high altitudes (>3500 m a.s.l.). Its strong preference for habitats with hot springs might be an early adaptation to the cooling climate during the plateau uplift. Some of these thermal sites may have been free of ice during the last glacial maximum (LGM) serving as refuges for the snake. To test this and to examine the population genetic structure data of 12 microsatellites and mitochondrial DNA were obtained from 153 individuals from 12 sites across the plateau. Four of the microsatellite loci were newly developed for this study. The results suggest that T. baileyi has at least two genetically diverse clades in Tibet, which developed during the Pleistocene and expanded after the LGM. Accordingly, the existence of separate glacial refuges on the central plateau can be assumed. Analyses of the genetic variation indicated a high level of geographic differentiation and population structure on a regional as well as on a range-wide scale. The study shows that, apart from the phylogeographic signatures, the diversification of current Thermophis populations is caused by (a) a limited dispersal due to mountain ranges, (b) a strong preference for hot springs and (c) the insular distribution of suitable habitats on the plateau. The present data should be taken into account in necessary monitoring and management plans for T. baileyi.

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).

The relationship between plasma steroid hormone concentrations and the reproductive cycle in the Northern Pacific rattlesnake, Crotalus oreganus

We describe the reproductive cycle of Northern Pacific rattlesnakes (Crotalus oreganus) by quantifying steroid hormone concentrations and observing reproductive behaviors in free-ranging individuals. Additionally, we examined reproductive tissues from museum specimens. Plasma steroid hormone concentrations were quantified for both male and female snakes throughout the active season (March-October). We measured testosterone (T), 5alpha-dihydrotestosterone (DHT), and corticosterone (B) concentrations in both sexes and 17beta-estradiol (E2) and progesterone (P) in females only. We observed reproductive behaviors (e.g., consortship, courtship, and copulation) in the field and measured testis and follicle size in male and female snakes from museum collections to relate steroid hormone concentrations to the timing of reproductive events. Our study revealed that C. oreganus in central California exhibits a bimodal pattern of breeding, with most mating behavior occurring in the spring and some incidences of mating behavior observed in late summer/fall. Each breeding period corresponded with elevated androgen (T or DHT) levels in males. Testes were regressed in the spring when the majority of reproductive behavior was observed in this population, and they reached peak volume in August and September during spermatogenesis. Although we did not detect seasonal variation in female hormone concentrations, some females had high E2 in the spring and fall, coincident with mating and with increased follicle size (indicating vitellogenesis) in museum specimens. Females with high E2 concentrations also had high T and DHT concentrations. Corticosterone concentrations in males and females were not related either to time of year or to concentrations of any other hormones quantified. Progesterone concentrations in females also did not vary seasonally, but this likely reflected sampling bias as females tended to be underground, and thus unobtainable, in summer months when P would be expected to be elevated during gestation. In females, P was positively correlated with T and DHT, and E2 was positively correlated with T.

Evolutionary trends in venom composition in the western rattlesnakes (Crotalus viridis sensu lato): toxicity vs. tenderizers

The Western Rattlesnake (Crotalus viridis sensu lato, now including Crotalus oreganus) is broadly distributed across the western half of the United States, northwestern Mexico and southwestern Canada, and eight subspecies are currently recognized. Although some venom characteristics have been noted for most subspecies, a systematic study of venoms from all subspecies has not been reported. Venom was extracted from snakes collected from approximate geographic range centers for all subspecies and analyzed using SDS-PAGE, MALDI-TOF mass spectrometry, enzyme and toxicity assays. Electrophoretic and mass spectrometric analyses demonstrated that small myotoxins, disintegrins and PLA(2) were abundant in most venoms. PIII and PI metalloproteinases ( approximately 54 kDa and 23 kDa, respectively) were common to all venoms except C. o. concolor, C. o. caliginis and C.o. helleri. Metalloproteinase activity was highest in C. o. cerberus and lowest in C. o. concolor venoms ( approximately 100-fold difference). Conversely, C. o. concolor venom was the most toxic and C. o. cerberus venom was least toxic (15-fold difference). In general, venoms with high metalloproteinase activity were less toxic (type I venoms), while venoms which were highly toxic showed low protease activity (type II venoms). Within the C. viridis/oreganus complex, these two extremes of venom compositional phenotypes are observed, and it appears that high metalloproteinase activity and high toxicity are incompatible qualities of these venoms. The functional significance of these biochemical characteristics likely relates to characteristics of prey consumed, and venoms with low metalloproteinase activity may constrain snake prey selection or foraging activity patterns.

Conservation genetics of the Chinese cobra (Naja atra) investigated with mitochondrial DNA sequences

We collected Chinese cobras (Naja atra) from one island (Dinghai) and four mainland (Huangshan, Lishui, Quanzhou, and Baise) populations in southeastern China, and used sequence data derived from the ND2 (1032 bp) and cytochrome b (1117 bp) genes and molecular variance estimates to investigate the population genetic structure of the species. Our sequence data show that: (1) the three eastern (Dinghai, Huangshan, and Lishui) populations are genetically segregated from the two southern (Quanzhou and Baise) populations; (2) the Quanzhou and Baise populations consist of two well-defined subclades, suggesting that the two populations have been well differentiated; (3) N. atra from the Huangshan population do not differ from those from the Lishui population, and lineage sorting in the northeastern part of the cobra's distributional range has not yet been completed because of the young age of Zhoushan Islands. The three eastern populations, the Quanzhou population, and the Baise population should be regarded as different management units (MUs). For these MUs, we suggest that in-situ protection measures should be taken because of their genetic uniqueness. Re-introductions or translocations are required to protect or re-establish natural populations of N. atra , but great care should be taken to enhance or retain local genetic variation.

Evolutionary morphology of the rattlesnake style

BACKGROUND

The rattlesnake rattling system is an evolutionary novelty that includes anatomical, behavioral, and physiological modifications of the generalized pitviper tail. One such modification, the formation of a bony clublike style at the terminal region of the caudal vertebrae, has not previously been examined in a phylogenetic context. Here we used skeletal material, cleared and stained preparations, and radiographs of whole preserved specimens to examine interspecific variation in style morphology among 34 rattlesnake species.

RESULTS

Evolutionary Principal Components Analysis revealed an inverse relationship between caudal segmental counts and style size, supporting the hypothesis that bone from caudal vertebral elements was reallocated to style formation during the evolution of this structure. Most of the basal rattlesnake species have small styles consisting of few compacted vertebral elements; however, early in the rattlesnake radiation there appears to have been two independent transitions to relatively large, pronged styles consisting of multiple coalesced vertebrae (once in Sistrurus catenatus, and once in Crotalus following the divergence of the Mexican long-tailed rattlesnakes). In terms of style shape, the two most divergent species, C. catalinensis and C. ericsmithi, provide insight into the possible relationship between style and rattle matrix morphology and lineage-specific evolutionary strategies for retaining rattle segments.

CONCLUSION

The considerable interspecific variation in rattle morphology appears to correspond to variation in the bony style. We hypothesize that style morphology evolves indirectly as an integrated module responding to adaptive evolution on matrix morphology.

Twenty years of phylogeography: the state of the field and the challenges for the Southern Hemisphere

Phylogeography is a young, vigorous and integrative field of study that uses genetic data to understand the history of populations. This field has recently expanded into many areas of biology and also into several historical disciplines of Earth sciences. In this review, I present a numerical synthesis of the phylogeography literature based on an examination of over 3000 articles published during the first 20 years of the field (i.e. from 1987 to 2006). Information from several topics needed to evaluate the progress, tendencies and deficiencies of the field is summarized for 10 major groups of organisms and at a global scale. The topics include the geography of phylogeographic surveys, comparative nature of studies, temporal scales and major environments investigated, and genetic markers used. I also identify disparities in research productivity between the developing and the developed world, and propose ways to reduce some of the challenges faced by phylogeographers from less affluent countries. Phylogeography has experienced explosive growth in recent years fuelled by developments in DNA technology, theory and statistical analysis. I argue that the intellectual maturation of the field will eventually depend not only on these recent developments, but also on syntheses of comparative information across different regions of the globe. For this to become a reality, many empirical phylogeographic surveys in regions of the Southern Hemisphere (and in developing countries of the Northern Hemisphere) are needed. I expect the information and views presented here will assist in promoting international collaborative work in phylogeography and in guiding research efforts at both regional and global levels.

Feeding ecology of the Great Basin Rattlesnake (Crotalus lutosus, Viperidae)

Documenting variation in organismal traits is essential to understanding the ecology of natural populations. We relied on stomach contents of preserved specimens and literature records to assess ontogenetic, intersexual, temporal, and geographic variations in the feeding ecology of the North American Great Basin Rattlesnake ( Crotalus lutosus Klauber, 1930). Snakes preyed mainly on rodents, occasionally on lizards, and less frequently on birds; squamate eggs and frogs were rarely eaten. There was a positive relationship between predator and prey size. The best predictors of this relationship were prey diameter as a function of snake body length and head size, underscoring the importance of prey diameter for gape-limited predators such as snakes. Crotalus lutosus displayed ontogenetic, sexual, and seasonal variations in diet. Smaller rattlesnakes fed predominantly on lizards, whereas larger individuals mostly fed on mammals. Females fed on lizards more often than males. The proportion of mammals in the diet was highest during the summer, a temporal variation that may be related to behavioral shifts in the diel activity and prey selectivity of C. lutosus, and (or) to differential abundance of rodents between seasons. Great Basin Rattlesnakes also displayed geographic variation in feeding habits, with snakes from the Great Basin Desert eating a higher proportion of lizards than serpents from the more northern Columbia Plateau.

Phylogeny, divergence times and species limits of spiny lizards (Sceloporus magister species group) in western North American deserts and Baja California

The broad distribution of the Sceloporus magister species group (squamata: phrynosomatidae) throughout western North America provides an appropriate model for testing biogeographical hypotheses explaining the timing and origins of diversity across mainland deserts and the Baja California Peninsula. We inferred concordant phylogenetic trees describing the higher-level relationships within the magister group using 1.6 kb of mitochondrial DNA (mtDNA) and 1.7 kb of nuclear DNA data. These data provide strong support for the parallel divergence of lineages endemic to the Baja California Peninsula (S. zosteromus and the orcutti complex) in the form of two sequential divergence events at the base of the magister group phylogeny. A relaxed phylogenetic analysis of the mtDNA data using one fossil and one biogeographical constraint provides a chronology of these divergence events and evidence that further diversification within the Baja California clades occurred simultaneously, although patterns of geographical variation and speciation between clades differ. We resolved four major phylogeographical clades within S. magister that (i) provide a novel phylogenetic placement of the Chihuahuan Desert populations sister to the Mojave Desert; (ii) illustrate a mixed history for the Colorado Plateau that includes Mojave and Sonoran Desert components; and (iii) identify an area of overlap between the Mojave and Sonoran Desert clades near Yuma, Arizona. Estimates of bidirectional migration rates among populations of S. magister using four nuclear loci support strong asymmetries in gene flow among the major mtDNA clades. Based on the nonexclusivity of mtDNA haplotypes, nuclear gene flow among populations and wide zones of phenotypic intergradation, S. magister appears to represent a single geographically variable and widespread species.

Phylogeography of the olive sea snake, Aipysurus laevis (Hydrophiinae) indicates Pleistocene range expansion around northern Australia but low contemporary gene flow

Pleistocene sea-level fluctuations profoundly changed landmass configurations around northern Australia. The cyclic emergence of the Torres Strait land bridge and concomitant shifts in the distribution of shallow-water marine habitats repeatedly sundered east and west coast populations. These biogeographical perturbations invoke three possible scenarios regarding the directions of interglacial range expansion: west to east, east to west, or bidirectional. We evaluated these scenarios for the olive sea snake, Aipysurus laevis, by exploring its genetic structure around northern Australia based on 354 individuals from 14 locations in three regions (Western Australia, WA; Gulf of Carpentaria, GoC; Great Barrier Reef, GBR). A 726-bp fragment of the mitochondrial DNA ND4 region revealed 41 variable sites and 38 haplotypes, with no shared haplotypes among the three regions. Population genetic structure was strong overall, phiST=0.78, P<0.001, and coalescent analyses revealed no migration between regions. Genetic diversity was low in the GBR and GoC and the genetic signatures of these regions indicated range or population expansions consistent with their recent marine transgressions around 7000 years ago. By contrast, genetic diversity on most WA reefs was higher and there were no signals of recent expansion events on these reefs. Phylogenetic analyses indicated that GBR and GoC haplotypes were derived from WA haplotypes; however, statistical parsimony suggested that recent range expansion in the GBR-GoC probably occurred from east coast populations, possibly in the Coral Sea. Levels of contemporary female-mediated gene flow varied within regions and reflected potential connectivity among populations afforded by the different regional habitat types.

Comparative mitochondrial genomics of snakes: extraordinary substitution rate dynamics and functionality of the duplicate control region

BACKGROUND

The mitochondrial genomes of snakes are characterized by an overall evolutionary rate that appears to be one of the most accelerated among vertebrates. They also possess other unusual features, including short tRNAs and other genes, and a duplicated control region that has been stably maintained since it originated more than 70 million years ago. Here, we provide a detailed analysis of evolutionary dynamics in snake mitochondrial genomes to better understand the basis of these extreme characteristics, and to explore the relationship between mitochondrial genome molecular evolution, genome architecture, and molecular function. We sequenced complete mitochondrial genomes from Slowinski's corn snake (Pantherophis slowinskii) and two cottonmouths (Agkistrodon piscivorus) to complement previously existing mitochondrial genomes, and to provide an improved comparative view of how genome architecture affects molecular evolution at contrasting levels of divergence.

RESULTS

We present a Bayesian genetic approach that suggests that the duplicated control region can function as an additional origin of heavy strand replication. The two control regions also appear to have different intra-specific versus inter-specific evolutionary dynamics that may be associated with complex modes of concerted evolution. We find that different genomic regions have experienced substantial accelerated evolution along early branches in snakes, with different genes having experienced dramatic accelerations along specific branches. Some of these accelerations appear to coincide with, or subsequent to, the shortening of various mitochondrial genes and the duplication of the control region and flanking tRNAs.

CONCLUSION

Fluctuations in the strength and pattern of selection during snake evolution have had widely varying gene-specific effects on substitution rates, and these rate accelerations may have been functionally related to unusual changes in genomic architecture. The among-lineage and among-gene variation in rate dynamics observed in snakes is the most extreme thus far observed in animal genomes, and provides an important study system for further evaluating the biochemical and physiological basis of evolutionary pressures in vertebrate mitochondria.