Molecular systematics and evolution of Regina and the thamnophiine snakes

Molecular phylogeny reveals cryptic diversity in Sibynophis from China (Serpentes: Sibynophiidae)

The elucidation of species diversity and distribution is critical within the fields of evolution, genetics, and conservation. The genus Sibynophis contains rare snakes that have historically received little attention. In this study, we conducted comprehensive sampling and use both mitochondrial and nuclear genetic markers to explore Sibynophis species diversity within China. Our findings revealed that S. c. miyiensis should be considered synonymous with S. c. grahami, and S. c. grahami should be gave a specific rank as S. grahami. In addition, we discovered S. triangularis was new to China and Myanmar. Based on the specimens and molecular phylogeny results, we redefined the species distribution boundaries of each Chinese species.

The Phylogenomic and Biogeographic History of the Gartersnakes, Watersnakes, and Allies (Natricidae: Thamnophiini)

North American Thamnophiini (gartersnakes, watersnakes, brownsnakes, and swampsnakes) are an ecologically and phenotypically diverse temperate clade of snakes representing 61 species across 10 genera. In this study, we estimate phylogenetic trees using ∼3,700 ultraconserved elements (UCEs) for 76 specimens representing 75% of all Thamnophiini species. We infer phylogenies using multispecies coalescent methods and time calibrate them using the fossil record. We also conducted ancestral area estimation to identify how major biogeographic boundaries in North America affect broadscale diversification in the group. While most nodes exhibited strong statistical support, analysis of concordant data across gene trees reveals substantial heterogeneity. Ancestral area estimation demonstrated that the genus Thamnophis was the only taxon in this subfamily to cross the Western Continental Divide, even as other taxa dispersed southward toward the tropics. Additionally, levels of gene tree discordance are overall higher in transition zones between bioregions, including the Rocky Mountains. Therefore, the Western Continental Divide may be a significant transition zone structuring the diversification of Thamnophiini during the Neogene and Pleistocene. Here we show that despite high levels of discordance across gene trees, we were able to infer a highly resolved and well-supported phylogeny for Thamnophiini, which allows us to understand broadscale patterns of diversity and biogeography.

Maximizing Molecular Data From Low-Quality Fluid-Preserved Specimens in Natural History Collections

Over the past decade, museum genomics studies have focused on obtaining DNA of sufficient quality and quantity for sequencing from fluid-preserved natural history specimens, primarily to be used in systematic studies. While these studies have opened windows to evolutionary and biodiversity knowledge of many species worldwide, published works often focus on the success of these DNA sequencing efforts, which is undoubtedly less common than obtaining minimal or sometimes no DNA or unusable sequence data from specimens in natural history collections. Here, we attempt to obtain and sequence DNA extracts from 115 fresh and 41 degraded samples of homalopsid snakes, as well as from two degraded samples of a poorly known snake, Hydrablabes periops. Hydrablabes has been suggested to belong to at least two different families (Natricidae and Homalopsidae) and with no fresh tissues known to be available, intractable museum specimens currently provide the only opportunity to determine this snake’s taxonomic affinity. Although our aim was to generate a target-capture dataset for these samples, to be included in a broader phylogenetic study, results were less than ideal due to large amounts of missing data, especially using the same downstream methods as with standard, high-quality samples. However, rather than discount results entirely, we used mapping methods with references and pseudoreferences, along with phylogenetic analyses, to maximize any usable molecular data from our sequencing efforts, identify the taxonomic affinity of H. periops, and compare sequencing success between fresh and degraded tissue samples. This resulted in largely complete mitochondrial genomes for five specimens and hundreds to thousands of nuclear loci (ultra-conserved loci, anchored-hybrid enrichment loci, and a variety of loci frequently used in squamate phylogenetic studies) from fluid-preserved snakes, including a specimen of H. periops from the Field Museum of Natural History collection. We combined our H. periops data with previously published genomic and Sanger-sequenced datasets to confirm the familial designation of this taxon, reject previous taxonomic hypotheses, and make biogeographic inferences for Hydrablabes. A second H. periops specimen, despite being seemingly similar for initial raw sequencing results and after being put through the same protocols, resulted in little usable molecular data. We discuss the successes and failures of using different pipelines and methods to maximize the products from these data and provide expectations for others who are looking to use DNA sequencing efforts on specimens that likely have degraded DNA.

Life Science Identifier (Hydrablabes periops)

urn:lsid:zoobank.org:pub:F2AA44 E2-D2EF-4747-972A-652C34C2C09D.

Phylogenomic analyses resolve relationships among garter snakes (Thamnophis: Natricinae: Colubridae) and elucidate biogeographic history and morphological evolution

Garter snakes (Thamnophis) are a successful group of natricines endemic to North America. They have become important natural models for ecological and evolutionary research, yet prior efforts to resolve phylogenetic relationships have resulted in conflicting topologies and weak support for certain relationships. Here, we use genomic data generated with a reduced representation double-digest RADseq approach to reassess evolutionary relationships across Thamnophis. We then use the resulting phylogeny to better understand how biogeography and feeding ecology have influenced lineage diversification and morphological evolution. We recovered highly congruent and strongly supported topologies from maximum likelihood and Bayesian analyses, but some discordance with a multispecies coalescent approach. All phylogenomic estimates split Thamnophis into two clades largely defined by northern and southern North American species. Divergence time estimates and biogeographic analyses indicate a mid-Miocene origin of Thamnophis in Mexico. In addition, historic vicariant events thought to explain biogeographic patterns in other lineages (e.g., Isthmus of Tehuantepec, Rocky Mountain Range, and Trans-Mexican Volcanic Belt) appear to have influenced patterns of diversification in Thamnophis as well. Analyses of morphological traits associated with feeding ecology showed moderate to strong phylogenetic signal. Nevertheless, phylogenetic ANOVA suggested significant differences in certain cranial morphologies between aquatic specialists and garter snakes that are terrestrial-aquatic generalists, independent of evolutionary history. Our new estimate of Thamnophis phylogeny yields an improved understanding of the biogeographic history and morphological evolution of garter snakes, and provides a robust framework for future research on these snakes.

Multilocus phylogeny, natural history traits and classification of natricine snakes (Serpentes: Natricinae)

Natricine snakes are geographically widespread, species rich (with ~250 extant species) and both morphologically and ecologically diverse. We present a multilocus DNA sequence phylogeny for 249 natricine specimens representing 189 named species, including 69 specimens and 21 species not previously sampled. Our inferred Bayesian and maximum likelihood trees form the basis for evaluations of genus-level classification, historical biogeography, lineage diversification, and dietary, habit and reproductive-mode diversity and evolution, although several, mostly deeper, relationships remain poorly resolved. The optimal trees support natricine origins in Asia, with dispersals to Australo-Melanesia, sub-Saharan Africa (including Seychelles Archipelago, excluding Aldabra), Europe and North Africa and into North and Central America. Viviparity appears to have evolved independently three times in Natricinae but was not significantly associated with an aquatic habit. We found limited associations between habit and diet categories. We propose generic reallocations for four natricine species and highlight other points of uncertainty in natricine classification.

Cryptic diversity and phylogeography of the Rhabdophis nuchalis group (Squamata: Colubridae)

Previous studies, have found that the rapid uplift of the Tibetan plateau accelerated the diversification of species. However, there are few relevant biogeographic data for the Colubridae in this region. We conducted a comprehensive study of the Rhabdophis nuchalis Group, which presently contains four nominal species, R. nuchalis, R. pentasupralabialis, R. leonardi, and R. chiwen. Building upon previous studies with specimens we have recently examined, greater interspecific and intraspecific diversity has been revealed. Here we address three questions: (1) Do the intraspecific differences represent only geographic variation within lineages, or are there cryptic species? (2) What are the interspecific relationships among members of the R. nuchalis Group? (3) What has been the biogeographic history of this species group? To resolve these questions we used four mitochondrial gene sequences and one nuclear sequence to investigate the molecular phylogenetic and geographic relationships among populations. Our molecular analysis reveals cryptic species diversity within the R. nuchalis Group, and seven clades were identified in the analysis. Ancestral area estimation suggests that the R. nuchalis Group originated in the Hengduan Mountains approximately 6.24 Mya and expanded its range northward to the Qinling-Daba Mountains. The Sichuan Basin appears to have been a barrier to migration. Species divergence seems to have been related to the rapid uplift of the Qinghai-Tibet Plateau.

Phylogenetics of mud snakes (Squamata: Serpentes: Homalopsidae): A paradox of both undescribed diversity and taxonomic inflation

Mud snakes (Serpentes: Homalopsidae) are a family of 55 described, mainly aquatic, species primarily distributed throughout mainland Southeast Asia and the Indo-Australian Archipelago. Although they have been the focus of prior research, the basic relationships amongst genera and species remain poorly known. We used a combined mitochondrial and nuclear gene dataset to infer their phylogenetic relationships, using the highest levels of taxon and geographic sampling for any homalopsid phylogeny to date (62% generic and 62% species coverage; 140 individuals). Our results recover two reciprocally monophyletic groups: the fangless Brachyorrhos and its sister clade comprised of all rear-fanged homalopsids. Most genera and interspecific relationships were monophyletic and strongly supported, but intergeneric relationships and intraspecific population structure lack support. We find evidence of both undescribed diversity as well as cases of taxonomic inflation within several species. Tree-based species delimitation approaches (mPTP) support potential new candidate species as distinct from their conspecifics and also suggest that many named taxa may not be distinct species. Divergence date estimation and lineage-through-time analyses indicate lower levels of speciation in the Eocene, with a subsequent burst in diversification in the Miocene. Homalopsids may have diversified most rapidly during the Pliocene and Pleistocene, possibly in relation to tectonic shifts and sea-level fluctuations that took place in Sundaland and the Sahul Shelf. Our analyses provide new insights on homalopsid taxonomy, a baseline phylogeny for the family, and further biogeographic implications demonstrating how dynamic tectonics and Quaternary sea level changes may have shaped a widespread, diverse family of snakes.

Eye size in North American watersnakes (genus Nerodia) correlates with variation in feeding ecology

Visual acuity and sensitivity positively correlate to eye size in vertebrates, and eye size relates to the ecology of colubrid snakes. We investigated whether eye morphology of North American colubrids of the genus Nerodia correlates with ecology as well. Although all members of the genus utilize aquatic habits, they differ widely in the proportion of anurans they eat. We specifically tested whether eye size and placement is associated with the proportion of frogs in the diet to determine whether these two aspects of eye morphology relate to feeding ecology. Using phylogenetic comparative methods, we found a significantly positive association between eye size and the proportion of anurans eaten by Nerodia species. Although the evidence is equivocal, the anterior placement of relatively small eyes in one species may also enhance anurophagy. Although eye size may improve a snake’s ability to feed on frogs, eye size must compete with other selective forces on head shape in trade-offs that may also influence eye size.

Out of the Hengduan Mountains: Molecular phylogeny and historical biogeography of the Asian water snake genus Trimerodytes (Squamata: Colubridae)

The Asian water snake genus Trimerodytes (formerly Sinonatrix) is endemic to East and Southeast Asia. Although several species have been included in various phylogenetic studies previously, the evolution and relationships among members of this genus as a whole remain unexplored. In this study, we report the sequencing two protein-coding mitochondrial gene fragments (MTCYB and ND2) and three nuclear genes (c-mos, NT3, and Rag1), reconstruct interspecific phylogeny, and explore biogeography for the genus Trimerodytes. Both Bayesian inference and maximum likelihood analyses consistently recover the monophyly of Trimerodytes with strong support, with T. yapingi the sister-group to the remaining species. The divergence date and ancestral area estimation suggest that Trimerodytes likely originated in Hengduan Mountains (eastern Tibetan Plateau) in western China at 23.93 Ma (95% HPD: 17.09-31.30), and intraspecific divergence began at about 4.23 Ma (95% HPD: 2.74-6.10). Analyses support the validity of T. yunnanensis.

Congenital feeding response to a novel prey in a Mexican gartersnake

In this study, we explored chemosensory, ingestive and prey-catching responses of neonate Mexican Black-bellied Gartersnakes (Thamnophis melanogaster) to crayfish (Cambarellus montezumae). By comparing snakes from a recently discovered crayfish-eating population and a typical non-crayfish-eating population, we asked which behavioral components change as a species enlarges its feeding niche. In the crayfish-eating population chemosensory responsiveness to crayfish was not enhanced but its heritability was higher. Neonates of both populations showed similar preference for freshly-molted versus unmolted crayfish, and whereas the tendency to ingest both crayfish stages remained stable between ages 15 and 90 days in the non-crayfish-eating population, in the crayfish-eating population it actually decreased. Techniques to catch and manipulate molted crayfish were similar in the two populations. We discuss the possibility that there is no increase in the behavioral response to eat crayfish by the neonates of the crayfish-eating populations, possibly due to the absence of ecological and spatial isolation between the two T. melanogaster populations. The crayfish ingestion in some population of T. melanogaster can be explained by environmental differences between populations, or by recent origin of crayfish ingestion in T. melanogaster.

The relationship of head morphology and diet among three sympatric watersnake species

Investigating dietary resource utilization and head morphology of gape-limited predators can provide a basis for understanding interspecific competition and species coexistence. For sympatric species, convergence of head morphology can indicate similar prey or foraging strategies while divergence can suggest competition, resource partitioning or expansion into new habitats. Sexual dimorphism can further complicate the head morphology-diet relationship by potentially reducing intraspecific competition. To understand the relationship between head morphology and diet, we studied three sympatric species, the plain-bellied (Nerodia erythrogaster), diamondback (N. rhombifer), and northern (N. sipedon) watersnakes in western Kentucky. All three species inhabit similar wetlands and feed on amphibians and fishes. The anurophagous plain-bellied watersnakes had longer, narrower heads that likely allow snakes to capture anuran metamorphs and froglets while facilitating movement through dense vegetation. The piscivorous diamondback watersnake had a wider head, which would enhance contact with prey in this open-mouth forager, and smaller interocular distance, resulting in dorsally placed eyes that facilitate seeing fish prey from below. The mostly piscivorous northern watersnake had an intermediate, generalized head shape, which may reflect the typically wide diet range of this species. Head shape did not differ between snake sexes but diamondback and northern females had larger heads than males and fed more frequently on fishes. The relationship between head morphology and diet highlights potential interspecific and intersexual competition and aids in the understanding of species coexistence patterns.

Evolution of nuchal glands, unusual defensive organs of Asian natricine snakes (Serpentes: Colubridae), inferred from a molecular phylogeny

A large body of evidence indicates that evolutionary innovations of novel organs have facilitated the subsequent diversification of species. Investigation of the evolutionary history of such organs should provide important clues for understanding the basis for species diversification. An Asian natricine snake, Rhabdophis tigrinus, possesses a series of unusual organs, called nuchal glands, which contain cardiotonic steroid toxins known as bufadienolides. Rhabdophis tigrinus sequesters bufadienolides from its toad prey and stores them in the nuchal glands as a defensive mechanism. Among more than 3,500 species of snakes, only 17 Asian natricine species are known to possess nuchal glands or their homologues. These 17 species belong to three nominal genera, Balanophis, Macropisthodon, and Rhabdophis. In Macropisthodon and Rhabdophis, however, species without nuchal glands also exist. To infer the evolutionary history of the nuchal glands, we investigated the molecular phylogenetic relationships among Asian natricine species with and without nuchal glands, based on variations in partial sequences of Mt-CYB, Cmos, and RAG1 (total 2,767 bp). Results show that all species with nuchal glands belong to a single clade (NGC). Therefore, we infer that the common ancestor of this clade possessed nuchal glands with no independent origins of the glands within the members. Our results also imply that some species have secondarily lost the glands. Given the estimated divergence time of related species, the ancestor of the nuchal gland clade emerged 19.18 mya. Our study shows that nuchal glands are fruitful subjects for exploring the evolution of novel organs. In addition, our analysis indicates that reevaluation of the taxonomic status of the genera Balanophis and Macropisthodon is required. We propose to assign all species belonging to the NGC to the genus Rhabdophis, pending further study.

Morphological convergence in a Mexican garter snake associated with the ingestion of a novel prey

Morphological convergence is expected when organisms which differ in phenotype experience similar functional demands, which lead to similar associations between resource utilization and performance. To consume prey with hard exoskeletons, snakes require either specialized head morphology, or to deal with them when they are vulnerable, for example, during molting. Such attributes may in turn reduce the efficiency with which they prey on soft‐bodied, slippery animals such as fish. Snakes which consume a range of prey may present intermediate morphology, such as that of Thamnophiine (Natricinae), which may be classified morphometrically across the soft–hard prey dietary boundary. In this study, we compared the dentition and head structure of populations of Thamnophis melanogaster that have entered the arthropod–crustacean (crayfish)‐eating niche and those that have not, and tested for convergence between the former and two distantly related crayfish specialists of the genus Regina (R. septemvittata and R. grahamii). As a control, we included the congener T. eques. Multivariate analysis of jaw length, head length, head width, and number of maxillary teeth yielded three significant canonical variables that together explained 98.8% of the variance in the size‐corrected morphological data. The first canonical variable significantly discriminated between the three species. The results show that head dimensions and number of teeth of the two Regina species are more similar to those of crayfish‐eating T. melanogaster than to non‐crayfish‐eating snakes or of T. eques. It is unclear how particular head proportions or teeth number facilitates capture of crayfish, but our results and the rarity of soft crayfish ingestion by T. melanogaster may reflect the novelty of this niche expansion, and are consistent with the hypothesis that some populations of T. melanogaster have converged in their head morphology with the two soft crayfish‐eating Regina species, although we cannot rule out the possibility of a morphological pre‐adaptation to ingest crayfish.

A Species-Level Phylogeny of Extant Snakes with Description of a New Colubrid Subfamily and Genus

BACKGROUND

With over 3,500 species encompassing a diverse range of morphologies and ecologies, snakes make up 36% of squamate diversity. Despite several attempts at estimating higher-level snake relationships and numerous assessments of generic- or species-level phylogenies, a large-scale species-level phylogeny solely focusing on snakes has not been completed. Here, we provide the largest-yet estimate of the snake tree of life using maximum likelihood on a supermatrix of 1745 taxa (1652 snake species + 7 outgroup taxa) and 9,523 base pairs from 10 loci (5 nuclear, 5 mitochondrial), including previously unsequenced genera (2) and species (61).

RESULTS

Increased taxon sampling resulted in a phylogeny with a new higher-level topology and corroborate many lower-level relationships, strengthened by high nodal support values (> 85%) down to the species level (73.69% of nodes). Although the majority of families and subfamilies were strongly supported as monophyletic with > 88% support values, some families and numerous genera were paraphyletic, primarily due to limited taxon and loci sampling leading to a sparse supermatrix and minimal sequence overlap between some closely-related taxa. With all rogue taxa and incertae sedis species eliminated, higher-level relationships and support values remained relatively unchanged, except in five problematic clades.

CONCLUSION

Our analyses resulted in new topologies at higher- and lower-levels; resolved several previous topological issues; established novel paraphyletic affiliations; designated a new subfamily, Ahaetuliinae, for the genera Ahaetulla, Chrysopelea, Dendrelaphis, and Dryophiops; and appointed Hemerophis (Coluber) zebrinus to a new genus, Mopanveldophis. Although we provide insight into some distinguished problematic nodes, at the deeper phylogenetic scale, resolution of these nodes may require sampling of more slowly-evolving nuclear genes.

Mitochondrial divergence between slow- and fast-aging garter snakes

Mitochondrial function has long been hypothesized to be intimately involved in aging processes--either directly through declining efficiency of mitochondrial respiration and ATP production with advancing age, or indirectly, e.g., through increased mitochondrial production of damaging free radicals with age. Yet we lack a comprehensive understanding of the evolution of mitochondrial genotypes and phenotypes across diverse animal models, particularly in species that have extremely labile physiology. Here, we measure mitochondrial genome-types and transcription in ecotypes of garter snakes (Thamnophis elegans) that are adapted to disparate habitats and have diverged in aging rates and lifespans despite residing in close proximity. Using two RNA-seq datasets, we (1) reconstruct the garter snake mitochondrial genome sequence and bioinformatically identify regulatory elements, (2) test for divergence of mitochondrial gene expression between the ecotypes and in response to heat stress, and (3) test for sequence divergence in mitochondrial protein-coding regions in these slow-aging (SA) and fast-aging (FA) naturally occurring ecotypes. At the nucleotide sequence level, we confirmed two (duplicated) mitochondrial control regions one of which contains a glucocorticoid response element (GRE). Gene expression of protein-coding genes was higher in FA snakes relative to SA snakes for most genes, but was neither affected by heat stress nor an interaction between heat stress and ecotype. SA and FA ecotypes had unique mitochondrial haplotypes with amino acid substitutions in both CYTB and ND5. The CYTB amino acid change (Isoleucine → Threonine) was highly segregated between ecotypes. This divergence of mitochondrial haplotypes between SA and FA snakes contrasts with nuclear gene-flow estimates, but correlates with previously reported divergence in mitochondrial function (mitochondrial oxygen consumption, ATP production, and reactive oxygen species consequences).

Projecting invasion risk of non-native watersnakes (Nerodia fasciata and Nerodia sipedon) in the western United States

Species distribution models (SDMs) are increasingly used to project the potential distribution of introduced species outside their native range. Such studies rarely explicitly evaluate potential conflicts with native species should the range of introduced species expand. Two snake species native to eastern North America, Nerodia fasciata and Nerodia sipedon, have been introduced to California where they represent a new stressor to declining native amphibians, fish, and reptiles. To project the potential distributions of these non-native watersnakes in western North America, we built ensemble SDMs using MaxEnt, Boosted Regression Trees, and Random Forests and habitat and climatic variables. We then compared the overlap between the projected distribution of invasive watersnakes and the distributions of imperiled native amphibians, fish, and reptiles that can serve as prey or competitors for the invaders, to estimate the risk to native species posed by non-native watersnakes. Large areas of western North America were projected to be climatically suitable for both species of Nerodia according to our ensemble SDMs, including much of central California. The potential distributions of both N. fasciata and N. sipedon overlap extensively with the federally threatened Giant Gartersnake, Thamnophis gigas, which inhabits a similar ecological niche. N. fasciata also poses risk to the federally threatened California Tiger Salamander, Ambystoma californiense, whereas N. sipedon poses risk to some amphibians of conservation concern, including the Foothill Yellow-legged Frog, Rana boylii. We conclude that non-native watersnakes in California can likely inhabit ranges of several native species of conservation concern that are expected to suffer as prey or competing species for these invaders. Action should be taken now to eradicate or control these invasions before detrimental impacts on native species are widespread. Our methods can be applied broadly to quantify the risk posed by incipient invasions to native biodiversity.

A phylogeny and revised classification of Squamata, including 4161 species of lizards and snakes

BACKGROUND

The extant squamates (>9400 known species of lizards and snakes) are one of the most diverse and conspicuous radiations of terrestrial vertebrates, but no studies have attempted to reconstruct a phylogeny for the group with large-scale taxon sampling. Such an estimate is invaluable for comparative evolutionary studies, and to address their classification. Here, we present the first large-scale phylogenetic estimate for Squamata.

RESULTS

The estimated phylogeny contains 4161 species, representing all currently recognized families and subfamilies. The analysis is based on up to 12896 base pairs of sequence data per species (average = 2497 bp) from 12 genes, including seven nuclear loci (BDNF, c-mos, NT3, PDC, R35, RAG-1, and RAG-2), and five mitochondrial genes (12S, 16S, cytochrome b, ND2, and ND4). The tree provides important confirmation for recent estimates of higher-level squamate phylogeny based on molecular data (but with more limited taxon sampling), estimates that are very different from previous morphology-based hypotheses. The tree also includes many relationships that differ from previous molecular estimates and many that differ from traditional taxonomy.

CONCLUSIONS

We present a new large-scale phylogeny of squamate reptiles that should be a valuable resource for future comparative studies. We also present a revised classification of squamates at the family and subfamily level to bring the taxonomy more in line with the new phylogenetic hypothesis. This classification includes new, resurrected, and modified subfamilies within gymnophthalmid and scincid lizards, and boid, colubrid, and lamprophiid snakes.

Gravity and the evolution of cardiopulmonary morphology in snakes

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

Refugial isolation and divergence in the Narrowheaded Gartersnake species complex (Thamnophis rufipunctatus) as revealed by multilocus DNA sequence data

Glacial-interglacial cycles of the Pleistocene are hypothesized as one of the foremost contributors to biological diversification. This is especially true for cold-adapted montane species, where range shifts have had a pronounced effect on population-level divergence. Gartersnakes of the Thamnophis rufipunctatus species complex are restricted to cold headwater streams in the highlands of the Sierra Madre Occidental and southwestern USA. We used coalescent and multilocus phylogenetic approaches to test whether genetic diversification of this montane-restricted species complex is consistent with two prevailing models of range fluctuation for species affected by Pleistocene climate changes. Our concatenated nuDNA and multilocus species analyses recovered evidence for the persistence of multiple lineages that are restricted geographically, despite a mtDNA signature consistent with either more recent connectivity (and introgression) or recent expansion (and incomplete lineage sorting). Divergence times estimated using a relaxed molecular clock and fossil calibrations fall within the Late Pleistocene, and zero gene flow scenarios among current geographically isolated lineages could not be rejected. These results suggest that increased climate shifts in the Late Pleistocene have driven diversification and current range retraction patterns and that the differences between markers reflect the stochasticity of gene lineages (i.e. ancestral polymorphism) rather than gene flow and introgression. These results have important implications for the conservation of T. rufipunctatus (sensu novo), which is restricted to two drainage systems in the southwestern US and has undergone a recent and dramatic decline.

A proposal to sequence the genome of a garter snake (Thamnophis sirtalis)

Here we develop an argument in support of sequencing a garter snake (Thamnophis sirtalis) genome, and outline a plan to accomplish this. This snake is a common, widespread, nonvenomous North American species that has served as a model for diverse studies in evolutionary biology, physiology, genomics, behavior and coevolution. The anole lizard is currently the only genome sequence available for a non-avian reptile. Thus, the garter snake at this time would be the first available snake genome sequence and as such would provide much needed comparative representation of non-avian reptilian genomes, and would also allow critical new insights for vertebrate comparative genomic studies. We outline the major areas of discovery that the availability of the garter snake genome would enable, and describe a plan for whole-genome sequencing.

Phylogeographic analysis and environmental niche modeling of the plain-bellied watersnake (Nerodia erythrogaster) reveals low levels of genetic and ecological differentiation

Species that exhibit geographically defined phenotypic variation traditionally have been divided into subspecies. Subspecies based on phenotypic features may not comprise monophyletic groups due to selection, gene flow, and/or convergent evolution. In many taxonomic groups the number of species once designated as widespread is dwindling rapidly, and many workers reject the concept of subspecies altogether. We tested whether currently recognized subspecies in the plain-bellied watersnake Nerodia erythrogaster are concordant with relationships based on mitochondrial markers, and whether it represents a single widespread species. The range of this taxon spans multiple potential biogeographic barriers (especially the Mississippi and Apalachicola Rivers) that correspond with lineage breaks in many species, including other snakes. We sequenced three mitochondrial genes (NADH-II, Cyt-b, Cox-I) from 156 geo-referenced specimens and developed ecological niche models using Maxent and spatially explicit climate data to examine historical and ecological factors affecting variation in N. erythrogaster across its range. Overall, we found little support for the recognized subspecies as either independent evolutionary lineages or geographically circumscribed units and conclude that although some genetic and niche differentiation has occurred, most populations assigned to N. erythrogaster appear to represent a single, widespread species. However, additional sampling and application of nuclear markers are necessary to clarify the status of the easternmost populations.