Molecular systematics and evolution of Regina and the thamnophiine snakes

The evolutionary origins of beneficial alleles during the repeated adaptation of garter snakes to deadly prey

Where do the genetic variants underlying adaptive change come from? Are currently adaptive alleles recruited by selection from standing genetic variation within populations, moved through introgression from other populations, or do they arise as novel mutations? Here, we examine the molecular basis of repeated adaptation to the toxin of deadly prey in 3 species of garter snakes (Thamnophis) to determine whether adaptation has evolved through novel mutations, sieving of existing variation, or transmission of beneficial alleles across species. Functional amino acid substitutions in the skeletal muscle sodium channel (Na(v)1.4) are largely responsible for the physiological resistance of garter snakes to tetrodotoxin found in their newt (Taricha) prey. Phylogenetic analyses reject the hypotheses that the unique resistance alleles observed in multiple Thamnophis species were present before the split of these lineages, or that alleles were shared among species through occasional hybridization events. Our results demonstrate that adaptive evolution has occurred independently multiple times in garter snakes via the de novo acquisition of beneficial mutations.

Distinctiveness in the face of gene flow: hybridization between specialist and generalist gartersnakes

Patterns of divergence and polymorphism across hybrid zones can provide important clues as to their origin and maintenance. Unimodal hybrid zones or hybrid swarms are composed predominantly of recombinant individuals whose genomes are patchworks of alleles derived from each parental lineage. In contrast, bimodal hybrid zones contain few identifiable hybrids; most individuals fall within distinct genetic clusters. Distinguishing between hybrid swarms and bimodal hybrid zones can be important for taxonomic and conservation decisions regarding the status and value of hybrid populations. In addition, the causes of bimodality are important in understanding the generation and maintenance of biological diversity. For example, are distinct clusters mostly reproductively isolated and co-adapted gene complexes, or can distinctiveness be maintained by a few 'genomic islands' despite rampant gene flow across much of the genome? Here we focus on three patterns of distinctiveness in the face of gene flow between gartersnake taxa in the Great Lakes region of North America. Bimodality, the persistence of distinct clusters of genotypes, requires strong barriers to gene flow and supports recognition of distinct specialist (Thamnophis butleri) and generalist (Thamnophis radix) taxa. Concordance of DNA-based clusters with morphometrics supports the hypothesis that trophic morphology is a key component of divergence. Finally, disparity in the level of differentiation across molecular markers (amplified fragment length polymorphisms) indicates that distinctiveness is maintained by strong selection on a few traits despite high gene flow currently or in the recent past.

What are the closest relatives of the hot-spring snakes (Colubridae, Thermophis), the relict species endemic to the Tibetan Plateau?

We conducted phylogenetic analyses to identify the closest related living relatives of the Xizang and Sichuan hot-spring snakes (T. baileyi and T. zhaoermii) endemic to the Tibetan Plateau, using mitochondrial DNA sequences (cyt b, ND4) from eight specimens, together with sequences from 95 additional caenophidian and five henophidian genera that were downloaded from GenBank. Phylogenetic trees were obtained using Bayesian Inference and Maximum likelihood methods. Results suggest that hot-spring snakes, which are adapted to high and cold environments, were clustered in the monophyletic Xenodontinae. Xenodontinae is one of the largest subfamilies of colubrid snakes, with about 90 genera and more than 500 species known, and are primarily tropical snakes previously thought to be restricted to the New World. Our data failed to provide any evidence that the New World xenodontines diverged from Thermophis and dispersed into the New World, also failed to suggest a colonization of Asia by New World xenodontines by dispersal from the New World. An alternative plausible scenario may be that Thermophis and the New World xenodontines evolved independently in Asia and America, respectively, after the divergence of their common ancestor. The divergence of the two species in Thermophis was caused by the barrier of the Hengduan Mountains, and the speciation had almost occurred when Tibetan Plateau attained present elevation.

Good species despite massive hybridization: genetic research on the contact zone between the watersnakes Nerodia sipedon and N. fasciata in the Carolinas, USA

Genomic markers generated with the amplified fragment length polymorphism method revealed extensive, panmictic-like hybridization along the narrow contact zone between the water snakes Nerodia sipedon and Nerodia fasciata in the Carolinas, USA. However, asymmetric distributions of diagnostic markers between both species and low frequencies of backcrossed hybrids with a high value of interspecific mixture infer selection against certain genotypes. This is consistent with a pronounced genetic and morphological preponderance of N. fasciata characters in the hybrid zone. Despite massive hybridization within the contact zone, the existence of nearly fixed genetic markers and the potential inferiority of certain hybrid genotypes support the species status of the two taxa and corroborate known, but nondiagnostic differences in morphology and ecology. This study stretches the applicability of species concepts to cases, where the genetic compatibility between two closely related species is very high, yet, they still evolve and persist as independent entities.

Morphological integration and adaptation in the snake feeding system: a comparative phylogenetic study

A long-standing hypothesis for the adaptive radiation of macrostomatan snakes is that their enlarged gape--compared to both lizards and basal snakes--enables them to consume "large" prey. At first glance, this hypothesis seems plausible, or even likely, given the wealth of studies showing a tight match between maximum consumed prey mass and head size in snakes. However, this hypothesis has never been tested within a comparative framework. We address this issue here by testing this hypothesis in 12 monophyletic clades of macrostomatan snakes using recently published phylogenies, published maximum consumed prey mass data and morphological measurements taken from a large sample of museum specimens. Our nonphylogenetically corrected analysis shows that head width--independent of body size--is significantly related to mean maximum consumed prey mass among these clades, and this relationship becomes even more significant when phylogeny is taken into account. Therefore, these data do support the hypothesis that head shape is adapted to prey size in snakes. Additionally, we calculated a phylogenetically corrected morphological variance-covariance matrix to examine the role of morphological integration during head shape evolution in snakes. This matrix shows that head width strongly covaries with both jaw length and out-lever length of the lower jaw. As a result, selection on head width will likely be associated with concomitant changes in jaw length and lower jaw out-lever length in snakes.

The energetic consequences of dietary specialization in populations of the garter snake, Thamnophis elegans

We investigated the intraspecific variation in digestive energetics between dietary specialist and generalist populations of the Western Terrestrial garter snake (Thamnophis elegans) in northern California. Coastal populations have a specialized diet of slugs and inland populations have a generalized diet of fish, anurans, mice and leeches. The difference in prey preference between the two populations is congenital, heritable and ontogenetically stable. To examine energetic specializations and trade-offs in these populations, we measured the net assimilation efficiency of each snake population on both slug (Ariolimax columbianus) and fish (Rhinichthys osculus) diets. The net assimilation efficiency was measured during digestion of a meal and continued until metabolic rate re-attained prefeeding levels. Coastal snakes were able to utilize 62% more of the ingested energy towards production from slug diets through both increased assimilation of nutrients and reduced digestive costs. For fish, assimilation and digestive costs were the same in both coastal and inland populations. These results support the hypothesis that snakes with specialized diets can evolve physiological traits to extract nutrients more efficiently. However, there was no apparent trade-off on the more generalized diet that was associated with this specialization.

Phylogenetic hypotheses for the turtle family Geoemydidae

The turtle family Geoemydidae represents the largest, most diverse, and most poorly understood family of turtles. Little is known about this group, including intrafamilial systematics. The only complete phylogenetic hypothesis for this family positions geoemydids as paraphyletic with respect to tortoises, but this arrangement has not been accepted by many workers. We compiled a 79-taxon mitochondrial and nuclear DNA data set to reconstruct phylogenetic relationships for 65 species and subspecies representing all 23 genera of the Geoemydidae. Maximum parsimony (MP) and maximum-likelihood (ML) analyses and Bayesian analysis produced similar, well-resolved trees. Our analyses identified three main clades comprising the tortoises (Testudinidae), the old-world Geoemydidae, and the South American geoemydid genus Rhinoclemmys. Within Geoemydidae, many nodes were strongly supported, particularly based on Bayesian posterior probabilities of the combined three-gene dataset. We found that adding data for a subset of taxa improved resolution of some deeper nodes in the tree. Several strongly supported groupings within the Geoemydidae demonstrate non-monophyly of some genera and possible interspecific hybrids, and we recommend several taxonomic revisions based on available evidence.

Evolutionary transformations of the fetal membranes of viviparous reptiles: a case study of two lineages

The reptilian placenta is a composite structure formed by a functional interaction between extraembryonic membranes and the maternal uterus. Study of placental structure of squamate reptiles over the past century has established that each of the multiple independent origins of placentation, which characterize the reproductive diversity of squamates, has resulted from the evolutionary transformation of these homologous structures. Because each evolutionary transformation is an independent novel relationship between maternal and embryonic tissues, the resulting placentae are not homologous, even though the individual components may be. The evolution of reptilian placentation should reveal much about evolutionary patterns and mechanisms because similar structural-functional systems have been transformed along parallel trajectories on multiple occasions. We compared extraembryonic membrane and placental development and pattern of embryonic nutrition in thamnophiine snakes and Pseudemoia lizards in the context of recent hypotheses of phylogenetic relationships. Two primary types of placentation, chorioallantoic and yolk sac, evolved in each lineage. Smooth, highly vascular regions of chorioallantoic placentation are indistinguishable homoplasies that evolved in parallel, likely to facilitate respiratory exchange. The yolk sac placenta of each lineage is specialized for histotrophic nutrient transfer, yet composition of these structures differs because of variation in the ancestral snakes and lizards. In addition, the omphalopleure that contributes to yolk sac placentation persists to later embryonic stages compared to oviparous outgroups, but the two lineages have evolved different structures that prevent replacement of the omphalopleure by the allantois. Each lineage has also evolved unique structural specializations of the chorioallantoic placenta.

Sweeping and striking: a kinematic study of the trunk during prey capture in three thamnophiine snakes

The trunk plays an obvious and important role in the prey capture behavior of many species of snake, yet trunk function during predatory strikes is poorly understood. Axial kinematics of three thamnophiine snakes (Thamnophis couchii, Thamnophis elegans and Nerodia rhombifer) were studied to quantify differences between sideways-directed and forward-directed attacks and to investigate strike diversity at relatively low phylogenetic levels. Feeding strikes were filmed at 60 Hz, and 13 points along the head and body were digitized. These points were used to calculate body segment displacement, rotation and velocity during predatory strikes. Kinematic analysis revealed significant differences in the foraging modes of these aquatic-feeding species. T. couchii displayed a stereotypical pre-strike posture in which the entire body was arranged in a series of loops directed towards the prey. Forward displacement of body segments sometimes occurred over the entire body in T. couchii but was restricted to the anterior one-third of the trunk in T. elegans and N. rhombifer. T. couchii and N. rhombifer both struck rapidly compared with T. elegans, although N. rhombifer typically had a short strike distance. N. rhombifer struck significantly faster than T. elegans. Aquatic prey capture diversity appears to reflect ecological diversity in thamnophiine snakes.

Molecular phylogeography of common garter snakes (Thamnophis sirtalis) in western North America: implications for regional historical forces

Complete ND2 and partial ND4 and cytochrome b mitochondrial DNA (mtDNA) sequences were analysed to evaluate the phylogeographic patterns of common garter snakes (Thamnophis sirtalis) in western North America. This species is widely distributed throughout North America, and exhibits extensive phenotypic variation in the westernmost part of its range. The overall phylogeographic pattern based on mtDNA sequences is concordant with results from studies of other species in this region, implicating historical vicariant processes during the Pleistocene and indicating bottleneck effects of recent dispersal into postglacial habitat. Indeed, the topology is statistically consistent with the hypothesis of both southern (Great Basin and California) and northern (Haida Gwaii) refugia. Specifically, we identified genetic breaks among three major clades: Northwest Coastal populations, Intermountain populations, and all California populations. The California clade contained the only other well-supported branching patterns detected; relationships among populations within the two northern clades were indistinguishable. These molecular splits contrast sharply with all prior geographical analyses of phenotypic variation in T. sirtalis in this region. Our results suggest that the extensive phenotypic variation in western T. sirtalis has been shaped more by local evolutionary forces than by shared common ancestry. Consequently, we consider all morphologically based subspecies designations of T. sirtalis in this region invalid because they do not reflect reciprocal monophyly of the mtDNA sequences.

Phylogenetic relationships of North American garter snakes (Thamnophis) based on four mitochondrial genes: how much DNA sequence is enough?

The clade of garter snakes (Thamnophis) includes some of the most abundant and well-studied snakes in North America. However, phylogenetic relationships within this group have been little studied. We used DNA sequences of four mitochondrial genes (cytochrome b and NADH dehydrogenase subunits 1, 2, and 4) to estimate relationships among 29 of the 31 recognized species of Thamnophis plus the related species Adelophis foxi. Both maximum parsimony (MP) and maximum-likelihood (ML) analyses of all these genes combined produced well-resolved trees with moderate (70-89%) to strong (90-100%) bootstrap support for most clades. MP and ML trees were very similar, with no strongly supported conflict between the two analyses. These analyses identify a clade of 12 species largely restricted to México (the "Mexican clade"), and a clade containing 15 species that collectively range from Central America to southern Canada (the "widespread clade"). These two groups are identified as sister taxa in both MP and ML analyses. A clade consisting of the ribbon snakes (T. sauritus and T. proximus) and the common garter snake (T. sirtalis) is placed as the sister group to all other Thamnophis (i.e., the Mexican + widespread clades) in our analyses. High bootstrap proportions at several levels in the tree support the inclusion of both Thamnophis validus, which has traditionally been placed in the genus Nerodia, and the poorly known species Adelophis foxi within Thamnophis. We used randomly sampled characters (i.e., standard bootstrapping) and randomly sampled contiguous blocks of characters to examine the effect of number of characters on resolution of and support for relationships within Thamnophis using MP. In general, these analyses indicate that we have reached a point of strongly diminishing returns with respect to the effect of adding mtDNA sequence characters for the current set of taxa; our sample of 3809 mtDNA characters is apparently "enough." The next steps to improve the phylogenetic estimate may be to add nuclear DNA sequences, morphology, or behavior, or to sequence additional mtDNA lineages within species.