Search

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., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Inc. All rights reserved.)

The highly diverse snake superfamily Elapoidea is considered to be a classic example of ancient, rapid radiation. Such radiations are challenging to fully resolve phylogenetically, with the highly diverse Elapoidea a case in point. Previous attempts at inferring a phylogeny of elapoids produced highly incongruent estimates of their evolutionary relationships, often with very low statistical support. We sought to resolve this situation by sequencing over 4,500 ultraconserved element loci from multiple representatives of every elapoid family/subfamily level taxon and inferring their phylogenetic relationships with multiple methods. Concatenation and multispecies coalescent based species trees yielded largely congruent and well-supported topologies. Hypotheses of a hard polytomy were not retained for any deep branches. Our phylogenies recovered Cyclocoridae and Elapidae as diverging early within Elapoidea. The Afro-Malagasy radiation of elapoid snakes, classified as multiple subfamilies of an inclusive Lamprophiidae by some earlier authors, was found to be monophyletic in all analyses. The genus Micrelaps was consistently recovered as sister to Lamprophiidae. We establish a new family, Micrelapidae fam. nov., for Micrelaps and assign Brachyophis to this family based on cranial osteological synapomorphy. We estimate that Elapoidea originated in the early Eocene and rapidly diversified into all the major lineages during this epoch. Ecological opportunities presented by the post-Cretaceous-Paleogene mass extinction event may have promoted the explosive radiation of elapoid snakes., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 The Author(s). Published by Elsevier Inc. All rights reserved.)

We evaluated the role of Quaternary climatic fluctuations on the demographic history and population structure of amphibian species endemic to the 'campo rupestre' in the Neotropics, evaluating their distributional shifts, demographic changes, and lineage formation from the end of Pleistocene to present. We chose two anurans endemic to the high-elevation 'campo rupestre' in the Espinhaço Range (ER) in northeastern and southeastern Brazil (Bokermannohyla alvarengai and Bokermannohyla oxente), as models to test the role of Quaternary climatic fluctuations over their distribution range in this region. We collected tissue samples throughout their distribution range and used statistical phylogeography to examine processes of divergence and population demography. We generated spatial-temporal reconstructions using Bayesian inference in a coalescent framework in combination with hind-cast projections of species distribution models (SDMs). We also used the results and literature information to test alternative diversification scenarios via approximate Bayesian computation (ABC). Our results show that Quaternary climatic fluctuations influenced the geographic ranges of both species showing population expansion during the last glacial maximum (LGM) and range contraction during interglacial periods, as inferred from selected ABC models and from past projections of SDMs. We recovered Pleistocene diversification for both species occuring in distinctly unique periods for each taxon. An older and range-restricted lineage was recovered in a geographically isolated geological massif, deserving conservation and further taxonomic study. The diversification and distribution of these amphibian species endemic to the Neotropical 'campo rupestre' were influenced by Quaternary climatic fluctuations. The expansion of cold adapted species restricted to higher elevations during glacial periods and their concomitant retraction during interglacial periods may have been crucial for producing patterns of species richness and endemism along elevation gradients in tropical and subtropical domains. Such processes may influence the evolution of the biota distributed in heterogeneous landscapes with varied topography., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Within many biomes, the cause of phylogeographic structure remains unknown even across regions throughout North America, including within the biodiverse Chihuahuan Desert. For example, little is known about population structure or the timing of diversification of Chihuahuan endemics. This is due largely to the lack of population genomic studies within this region. We generated ultra-conserved element data for the gray-banded kingsnake (Lampropeltis alterna) to investigate lineage divergence and historical demography across the Chihuahuan Desert. We found three unique lineages corresponding to the Trans-Pecos and Mapimian biogeographic regions of the Chihuahuan Desert, and a distinct population in the Sierra Madre Occidental. Using several mutation rates to calibrate the timing of divergence among these lineages, we show that lineage divergence likely occurred during the Pleistocene, which indicates that careful consideration needs to be used when applying mutation rates to ultra-conserved elements. We suggest that biogeographic provinces within the Chihuahuan Desert may have served as allopatric refugia during climatic fluctuations of the Quaternary. This work serves as an important template for further testing biogeographic hypotheses within the region., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Globally distributed groups may show regionally distinct rates of diversification, where speciation is elevated given timing and sources of ecological opportunity. However, for most organisms, nearly complete sampling at genomic-data scales to reduce topological error in all regions is unattainable, thus hampering conclusions related to biogeographic origins and rates of diversification. We explore processes leading to the diversity of global ratsnakes and test several important hypotheses related to areas of origin and enhanced diversification upon colonizing new continents. We estimate species trees inferred from phylogenomic scale data (304 loci) while exploring several strategies that consider topological error from each individual gene tree. With a dated species tree, we examine taxonomy and test previous hypotheses that suggest the ratsnakes originated in the Old World (OW) and dispersed to New World (NW). Furthermore, we determine if dispersal to the NW represented a source of ecological opportunity, which should show elevated rates of species diversification. We show that ratsnakes originated in the OW during the mid-Oligocene and subsequently dispersed to the NW by the mid-Miocene; diversification was also elevated in a subclade of NW taxa. Finally, the optimal biogeographic region-dependent speciation model shows that the uptick in ratsnake diversification was associated with colonization of the NW. We consider several alternative explanations that account for regionally distinct diversification rates., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Several biogeographic barriers in the Eastern Nearctic appear to reduce gene flow among populations of many species in predictable ways, however these patterns used to infer process of divergence may be deceiving if alternative modes of diversification are not considered. By using a multilocus statistical phylogeographic approach to examine diversity within a North American snake, Lampropeltis calligaster, we find that mode and timing of speciation near the Mississippi River embayment and peninsular Florida, two main biodiversity hotspots in eastern North America, challenge previously held notions of strict vicariance as the causal factor behind patterns of divergence seen among taxa at these locations. We found three species inhabiting distinct ecological niches with divergences dating to the mid- and early-Pleistocene with subsequently stable or increasing effective population sizes, further supporting the idea that the Pleistocene was an important driver of diversification in North America. Our results lead to a revised hypothesis that ecological divergence has occurred in this group across environments associated with the Mississippi River and at the Florida peninsula. Importantly, in their western distributions, we show that species divergence is associated with the ecological transition from distinct forested habitats to grasslands, rather than the nearby Mississippi River, a barrier often implicated for many other organisms. Additionally, we stress the importance of examining each delimited lineage with respect to conservation, since ecological niche models suggest that by the end of the century changes in climate may negatively alter habitat suitability and, barring adaptation, substantially reduce the suitable range of two of the three species we identified., (Copyright © 2016 Elsevier Inc. All rights reserved.)

Next-generation genomic sequencing promises to quickly and cheaply resolve remaining contentious nodes in the Tree of Life, and facilitates species-tree estimation while taking into account stochastic genealogical discordance among loci. Recent methods for estimating species trees bypass full likelihood-based estimates of the multi-species coalescent, and approximate the true species-tree using simpler summary metrics. These methods converge on the true species-tree with sufficient genomic sampling, even in the anomaly zone. However, no studies have yet evaluated their efficacy on a large-scale phylogenomic dataset, and compared them to previous concatenation strategies. Here, we generate such a dataset for Caenophidian snakes, a group with >2500 species that contains several rapid radiations that were poorly resolved with fewer loci. We generate sequence data for 333 single-copy nuclear loci with ∼100% coverage (∼0% missing data) for 31 major lineages. We estimate phylogenies using neighbor joining, maximum parsimony, maximum likelihood, and three summary species-tree approaches (NJst, STAR, and MP-EST). All methods yield similar resolution and support for most nodes. However, not all methods support monophyly of Caenophidia, with Acrochordidae placed as the sister taxon to Pythonidae in some analyses. Thus, phylogenomic species-tree estimation may occasionally disagree with well-supported relationships from concatenated analyses of small numbers of nuclear or mitochondrial genes, a consideration for future studies. In contrast for at least two diverse, rapid radiations (Lamprophiidae and Colubridae), phylogenomic data and species-tree inference do little to improve resolution and support. Thus, certain nodes may lack strong signal, and larger datasets and more sophisticated analyses may still fail to resolve them., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Africa is renowned for its biodiversity and endemicity, yet little is known about the factors shaping them across the continent. African Agama lizards (45 species) have a pan-continental distribution, making them an ideal model for investigating biogeography. Many species have evolved conspicuous sexually dimorphic traits, including extravagant breeding coloration in adult males, large adult male body sizes, and variability in social systems among colorful versus drab species. We present a comprehensive time-calibrated species tree for Agama, and their close relatives, using a hybrid phylogenetic-phylogenomic approach that combines traditional Sanger sequence data from five loci for 57 species (146 samples) with anchored phylogenomic data from 215 nuclear genes for 23 species. The Sanger data are analyzed using coalescent-based species tree inference using (*)BEAST, and the resulting posterior distribution of species trees is attenuated using the phylogenomic tree as a backbone constraint. The result is a time-calibrated species tree for Agama that includes 95% of all species, multiple samples for most species, strong support for the major clades, and strong support for most of the initial divergence events. Diversification within Agama began approximately 23 million years ago (Ma), and separate radiations in Southern, East, West, and Northern Africa have been diversifying for >10Myr. A suite of traits (morphological, coloration, and sociality) are tightly correlated and show a strong signal of high morphological disparity within clades, whereby the subsequent evolution of convergent phenotypes has accompanied diversification into new biogeographic areas., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Delimiting species and clarifying phylogenetic relationships are the main goals of systematics. For species with questionable taxonomic status, species delimitation approaches using multi-species coalescent models with multiple loci are recommended if morphological data are unavailable or unhelpful. Moreover, these methods will also reduce subjectivity based on genetic distance or requirement of monophyletic genetic lineages. We determine the validity and phylogenetic position of a rare and long controversial species of Chinese reptile, the Szechwan ratsnake (Euprepiophis perlaceus), using multi-locus data from multiple individuals and coalescent-based approaches. Species were first delimited using Bayesian Phylogenetics & Phylogeography (BP&P), Brownie and Bayes Factor model comparison approaches, while relationships among species were estimated using species tree inference in (*)BEAST. Results indicate that Euprepiophis perlaceus is a distinct species sister to Euprepiophis mandarinus. Despite gene tree discrepancy, the coalescent model-based approaches used here demonstrate the taxonomic validity and the phylogenetic position of Euprepiophis perlaceus. These approaches objectively test the validity of questionable species diagnoses based on morphological characters and determine their phylogenetic position., (Copyright © 2013 Elsevier Inc. All rights reserved.)

Based on a molecular phylogeny and a comparison of maxillary dentition and morphology, the relationship between the genera Lycodon and Dinodon was investigated. Bayesian Inference and Maximum Likelihood analysis of two mitochondrial genes (cyt b and ND4) and two nuclear genes (c-mos and Rag1) suggested that the two genera shared a most recent common ancestor. However, Dinodon was paraphyletic and Lycodon was polyphyletic, each with respect to the other. The results from counts of maxillary teeth indicated that the diagnostic characters used by previous authors to separate Dinodon and Lycodon were not reliable. Taking the molecular and morphological evidence together, we synonymized Dinodon with Lycodon. In addition, the validity of the species L. futsingensis was confirmed to be distinctly different from the other species of Dinodon and Lycodon., (Copyright © 2013 Elsevier Inc. All rights reserved.)

Inferring the evolutionary and biogeographic history of taxa occurring in a particular region is one way to determine the processes by which the biodiversity of that region originated. Tree boas of the genus Corallus are an ancient clade and occur throughout Central and South America and the Lesser Antilles, making it an excellent group for investigating Neotropical biogeography. Using sequenced portions of two mitochondrial and three nuclear loci for individuals of all recognized species of Corallus, we infer phylogenetic relationships, present the first molecular analysis of the phylogenetic placement of the enigmatic C. cropanii, develop a time-calibrated phylogeny, and explore the biogeographic history of the genus. We found that Corallus diversified within mainland South America, via over-water dispersals to the Lesser Antilles and Central America, and via the traditionally recognized Panamanian land bridge. Divergence time estimates reject the South American Caribbean-Track as a general biogeographic model for Corallus and implicate a role for events during the Oligocene and Miocene in diversification such as marine incursions and the uplift of the Andes. Our findings also suggest that recognition of the island endemic species, C. grenadensis and C. cookii, is questionable as they are nested within the widely distributed species, C. hortulanus. Our results highlight the importance of using widespread taxa when forming and testing biogeographic hypotheses in complex regions and further illustrate the difficulty of forming broadly applicable hypotheses regarding patterns of diversification in the Neotropical region., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Snake diversity in the island of Sri Lanka is extremely high, hosting at least 89 inland (i.e., non-marine) snake species, of which at least 49 are endemic. This includes the endemic genera Aspidura, Balanophis, Cercaspis, Haplocercus, and Pseudotyphlops, which are of uncertain phylogenetic affinity. We present phylogenetic evidence from nuclear and mitochondrial loci showing the relationships of 40 snake species from Sri Lanka (22 endemics) to the remaining global snake fauna. To determine the phylogenetic placement of these species, we create a molecular dataset containing 10 genes for all global snake genera, while also sampling all available species for genera with endemic species occurring in Sri Lanka. Our sampling comprises five mitochondrial genes (12S, 16S, cyt-b, ND2, and ND4) and five nuclear genes (BDNF, c-mos, NT3 RAG-1, and RAG-2), for a total of up to 9582bp per taxon. We find that the five endemic genera represent portions of four independent colonizations of Sri Lanka, with Cercaspis nested within Colubrinae, Balanophis in Natricinae, Pseudotyphlops in Uropeltidae, and that Aspidura+Haplocercus represents a distinct, ancient lineage within Natricinae. We synonymize two endemic genera that render other genera paraphyletic (Haplocercus with Aspidura, and Cercaspis with Lycodon), and discover that further endemic radiations may be present on the island, including a new taxon from the blindsnake family Typhlopidae, suggesting a large endemic radiation. Despite its small size relative to other islands such as New Guinea, Borneo, and Madagascar, Sri Lanka has one of the most phylogenetically diverse island snake faunas in the world, and more research is needed to characterize the island's biodiversity, with numerous undescribed species in multiple lineages., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Numerous taxa show ancient intertropical disjunct distributions. Many can be explained by well-known processes of historical vicariance, such as the breakup of Gondwanaland. Others, such as Asian-Neotropical divergences are not as well understood. To clarify the phylogenetic position and understand biogeographic and temporal origins of the geographically disjunct and morphologically unique genera of hinged-teeth snakes, Scaphiodontophis (n=1) and Sibynophis (n=9; Colubridae), we inferred a time-calibrated phylogeny with additional 107 taxa representing the superfamily Colubroidea using four genes (c-mos, cyt-b, ND2, RAG-1; 3085 bp). We used this tree to estimate ancestral areas for the group. The results show that Scaphiodontophis is sister to Sibynophis, both originated in the late Eocene/Oligocene in Asia and likely dispersed through Beringia to the New World, but unlike other snake groups left no extant species in temperate North America. Current recognition of Scaphiodontophiinae renders Colubrinae paraphyletic, and we resurrect the previously named subfamily Sibynophiinae to encompass both genera and use the tribes Sibynophiini (Sibynophis) and Scaphiodontophiini (Scaphiodontophis) to highlight the geographically distinct areas occupied by these taxa. These results suggest that intercontinental dispersal with extinction in intermediate areas can explain puzzling patterns of ancient intertropical disjunct distributions., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Based on a combination of six mitochondrial gene fragments (12S RNA, cyt b, ND1, ND2, ND4 and CO1) and one nuclear gene (c-mos) from 22 genera we infer phylogenetic relationships among natricine snakes and examine the date and area of origin of these snakes. Our phylogenetic results indicate: (1) the subfamily Natricinae is strongly supported as monophyletic including a majority of extant genera, and a poorly known and previously unassigned species Trachischium monticola, (2) two main clades are inferred within Natricinae, one containing solely taxa from the Old World (OW) and the other comprising taxa from a monophyletic New World (NW) group with a small number of OW relatives. Within the first clade, the genera Xenochrophis and Amphiesma are apparently not monophyletic. Divergence dating and ancestral area estimation indicate that the natricines originated in tropical Asia during the later Eocene or the Oligocene. We recover two major dispersals events out of Asia, the first to Africa in the Oligocene (28 Ma) and the second to the Western Palearctic and the New World at 27 Ma. This date is consistent with the dispersal of numerous other OW groups into the NW., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Studies of speciation and taxon delimitation are usually decoupled. Combining these methods provides a stronger theoretical ground for recognizing new taxa and understanding processes of speciation. Using coalescent methods, we examine speciation, post-speciation population demographics, and taxon delimitation in the Arizona Mountain Kingsnake (Lampropeltis pyromelana), a species restricted to high elevations in southwestern United States and northern Mexico (SW). These methods provide a solid foundation for understanding how biogeographic barriers operate at the regional scale in the SW. Bayesian species delimitation methods, using three loci from samples of L. pyromelana taken throughout their range, show strong support for the existence of two species that are separated by low elevation habitats found between the Colorado Plateau/ Mogollon Rim and the Sierra Madre Occidental. Our results suggest an allopatric mode of speciation given the near absence of gene flow over time, which resulted in two lineages of unequal population sizes. Speciation likely occurred prior to the Pleistocene, during the aridification of the SW and/or the uplift of the Colorado Plateau, and while these species occupy similar high-elevation niches, they are isolated by xeric conditions found in the intervening low deserts. Furthermore, post-speciation demographics suggest that populations of both lineages were not negatively impacted by climate change throughout the Pleistocene. Finally, our results suggest that at least for this group, where divergence is old and gene flow is low, Bayesian species delimitation performs well., (Copyright © 2011 Elsevier Inc. All rights reserved.)

Spatial and temporal heterogeneity in environmental factors can have profound effects on diversification in species that are tightly linked to their environments. The Caddo Mountain Salamander (Plethodon caddoensis) inhabits a unique physiographic section of the Ouachita Mountains in central North America, a region in which Pleistocene climatic fluctuations have been implicated in driving lineage diversification in two other closely related salamanders. We examined P. caddoensis to determine whether it was similarly impacted by historic climatic changes and test whether physiographic features unique to the area also contributed to its diversification. We found that P. caddoensis is composed of four highly divergent, geographically distinct lineages that abut one another along an east-west axis. Phylogeographic structure was significantly related to both geographic distance and stream drainages, indicating that connectivity of streams and stream-associated habitats (e.g., talus) influence patterns of interpopulation gene flow. Lineages originated during the Middle Miocene and population size decreased in all lineages during the Pleistocene. Surface Geology and precipitation were the most important variables predicting the species distribution. Our results show that the unique physiographic features of the area coupled with species response to climatic factors have driven lineage diversification and phylogeographic structure in P. caddoensis. Variation in responses to historic climatic fluctuations among salamander species in this region underscore the importance of integrating species ecology with other factors such as geology and hydrology in order to better understand the effects of climate change on species with close associations to their environments., (Copyright © 2011 Elsevier Inc. All rights reserved.)

The superfamily Colubroidea (> 2500 species) includes the majority of snake species and is one of the most conspicuous and well-known radiations of terrestrial vertebrates. However, many aspects of the phylogeny of the group remain contentious, and dozens of genera have yet to be included in molecular phylogenetic analyses. We present a new, large-scale, likelihood-based phylogeny for the colubroids, including 761 species sampled for up to five genes: cytochrome b (93% of 761 species sampled), ND4 (69%), ND2 (28%), c-mos (54%), and RAG-1 (13%), totaling up to 5814bp per species. We also compare likelihood bootstrapping and a recently proposed ultra-fast measure of branch support (Shimodaira-Hasegawa-like [SHL] approximate likelihood ratio), and find that the SHL test shows strong support for several clades that were weakly-supported by bootstrapping in this or previous analyses (e.g., Dipsadinae, Lamprophiidae). We find that SHL values are positively related to branch lengths, but show stronger support for shorter branches than bootstrapping. Despite extensive missing data for many taxa (mean=67% per species), neither bootstrap nor SHL support values for terminal species are related to their incompleteness, and that most highly incomplete taxa are placed in the expected families from previous taxonomy, typically with very strong support. The phylogeny indicates that the Neotropical colubrine genus Scaphiodontophis represents an unexpectedly ancient lineage within Colubridae. We present a revised higher-level classification of Colubroidea, which includes a new subfamily for Scaphiodontophis (Scaphiodontophiinae). Our study provides the most comprehensive phylogeny of Colubroidea to date, and suggests that SHL values may provide a useful complement to bootstrapping for estimating support on likelihood-based trees., (Copyright © 2010 Elsevier Inc. All rights reserved.)

The impact of barriers to dispersal and gene flow is often inferred to be the primary cause of lineage divergence and phylogeographic structure in terrestrial organisms. In particular, the Mississippi River has been implicated as a barrier to gene flow in many species, including aquatic taxa. However, if barriers are permeable to organisms, then phylogeographic structure may be difficult to detect due to gene flow between lineages. Using time-calibrated Bayesian phylogenetic analyses of mtDNA, and phylogeographic coalescent simulations, we determine if the Mississippi River operates as a barrier to gene flow in the aquatic diamond-backed watersnake (Nerodia rhombifer). The phylogenetic analyses support a basal division within N. rhombifer mtDNA lineages that coincides with populations generally east and west of the Mississippi River. These results, and that of the divergence dating analyses, therefore suggest that the river was a significant barrier to gene flow in the Pleistocene ∼ 1.4 million years ago, presumably during an interglacial period when the river was much wider. However, we also detect western haplotypes in the eastern clade, and vice versa, thereby indicating that this barrier has not been complete. Nonetheless, the coalescent simulations that account for limited migration suggest that the Mississippi River was an important feature that shaped the phylogeographic history of this aquatic snake in the USA despite limited gene flow., (Copyright © 2010 Elsevier Inc. All rights reserved.)

Many studies have revealed that lineages currently inhabiting formerly glaciated areas were pushed into southern glacial refugia and have expanded into their modern range since the last glacial maximum. There have been few studies that compare the effects of glacial cycles on lineage diversification and historical demography in closely related species with overlapping ranges. In this study we compare phylogeographic structure, historical demography, and approximate lineage age in two closely related and broadly co-occurring venomous snakes in eastern North America, the cottonmouth (Agkistrodon piscivorus) and copperhead (A. contortrix) using sequences from the mtDNA gene cytochrome b. We inferred three geographic lineages of A. contortrix and two of A. piscivorus with no common geographic or temporal pattern of lineage diversification identified for these species. Lineage diversification occurred in the Late Pliocene for A. piscivorus (approximately 2.5mya) and in the Early Pleistocene for A. contortrix ( approximately 1.5mya). Demographic estimates revealed population expansion following the last glacial maximum (approximately 20,000 years ago) in two lineages of A. contortrix (the Central clade and Eastern clade) and one lineage of A. piscivorus (the Continental clade). The Florida clade of A. piscivorus is the only lineage for which constant population size through time was inferred, possibly due to stable populations persisting in areas unaffected by glacial advances. Our data suggest that unique habitat preferences may have shaped both the phylogeographic and demographic histories of each species.

Most phylogeographic studies examine organisms that do not have transcontinental distributions and therefore the genetic and temporal effects of barriers across an entire continent cannot be assessed with respect to a single species. We examined the phylogeographic structure, lineage age, and historical demography using sequences from the mtDNA cytochrome b gene of the widespread North American racer (Coluber constrictor), one of the few abundant transcontinental snakes that occurs throughout many diverse biomes. Our results indicate that this complex is comprised of six lineages differing greatly in geographic extent, with the largest (a central US clade) being approximately 26 times greater than the smallest (a lineage restricted to the Florida Panhandle and nearby portions of adjacent States). Most of the six lineages appear to be separated at previously identified genetic barriers for several vertebrates with similar ranges. Lineage diversification in this species began in the late Miocene, separating populations in the Florida Peninsula from the remainder of the US. Diversification of lineages continued throughout the Pliocene and early Pleistocene. Four of the six lineages occur east of the Mississippi River, with only two distinctly young ( approximately 1.5 mya) lineages found west of the Mississippi River (one occurs west of Continental Divide). All methods of demographic inference, including the mismatch distribution, Fu and Li's D and Tajima's D , and Bayesian skyline plots revealed population expansion occurring in the mid-to-late Pleistocene for every lineage, regardless of size or proximity to formerly glaciated areas. Population expansion for lineages found east of the Mississippi River occurred earlier and was much greater than those found west of the River.

Dynamic climatic oscillations during the Pleistocene had profound effects on the distributions of species across North America. Although the role of historical climate change on speciation remains controversial, the impact on genetic variation within species has been well documented. We examined mtDNA sequences from the cytochrome b gene (1117 bp) and a portion of the NADH-4 gene (659 bp) for 286 individuals of Diadophis punctatus to infer phylogeographic patterns and population structure and to examine historical demographic patterns in both glaciated and unglaciated regions of North America. We inferred 14 lineages that replace each other geographically across the United States. Several of these lineages appear to be confined to specific habitats (floodplains, grasslands, montane environments) and traverse previously identified genetic barriers for terrestrial vertebrates including the Mississippi and Apalachicola Rivers, the Appalachian Mountains, and the western continental divide. We also observed overlapping ranges between some haplotype groups and several instances of secondary contact associated with ecological transition zones in eastern South Carolina, southern Oklahoma and central California. Within the US, diversification began during the late Miocene and continued into the mid-Pleistocene, suggesting these lineages pre-dated the last glacial maximum. Coalescent and non-coalescent demographic analyses indicate that independent lineages currently occupying previously glaciated or unsuitable areas in eastern, central and western US underwent post-glacial population expansion likely from southern refugia during the late Pleistocene/early Holocene. Conversely, southern lineages display patterns consistent with long-term population stability. Such long-term persistence of genetic structure may be due to the competitive effects between lineages or ecosystem stability in more southern latitudes.

The systematic relationships of snakes in the Lampropeltis mexicana complex (L. mexicana, L. alterna, and L. ruthveni) are poorly known despite several taxonomic studies over the last 80 years. Mitochondrial DNA sequences were used to infer the phylogeny of the L. mexicana complex. At least one representative sample from the nine currently recognized species of Lampropeltis was sequenced. Our results suggest that a deep basal split resulted in the divergence of two groups of Lampropeltis, with one group occupying the upland areas of western United States and most of western and central Mexico, and the other northeastern Mexico and the lowland areas of the southern United States. Results also revealed that the L. mexicana complex and Lampropeltis triangulum are polyphyletic, with taxa from both groups nested together in deeply divergent northern and southern clades. These results are incongruent with previous hypotheses of phylogenetic relationships based on morphology, and suggest that morphological characters shared among the various tri-colored Lampropeltis (e.g., hemipenal structure and tri-colored pattern) may be difficult to interpret phylogenetically.

To examine Holarctic snake dispersal, we inferred a phylogenetic tree from four mtDNA genes and one scnDNA gene for most species of the Old World (OW) and New World (NW) colubrid group known as ratsnakes. Ancestral area distributions are estimated for various clades using divergence-vicariance analysis and maximum likelihood on trees produced using Bayesian inference. Dates of divergence for the same clades are estimated using penalized likelihood with statistically crosschecked calibration references obtained from the Miocene fossil record. With ancestral areas and associated dates estimated, various hypotheses concerning the age and environment associated with the origin of ratsnakes and the dispersal of NW taxa from OW ancestors were tested. Results suggest that the ratsnakes originated in tropical Asia in the late Eocene and subsequently dispersed to the Western and Eastern Palearctic by the early Oligocene. These analyses also suggest that the monophyletic NW ratsnakes (the Lampropeltini) diverged from OW ratsnakes and dispersed through Beringia in the late Oligocene/early Miocene when this land bridge was mostly composed of deciduous and coniferous forests.

The Colubroidea contains over 85% of all the extant species of snakes and is recognized as monophyletic based on morphological and molecular data. Using DNA sequences (cyt b, c-mos) from 100 species we inferred the phylogeny of colubroids with special reference to the largest family, the Colubridae. Tree inference was obtained using Bayesian, likelihood, and parsimony methods. All analyses produced five major groups, the Pareatidae, Viperidae, Homalopsidae, the Elapidae, and the Colubridae. The specific content of the latter two groups has been altered to accommodate evolutionary history and to yield a more stable taxonomy. We propose an updated classification based on the reallocation of species as indicated by our inferred phylogeny.

Snakes of the subfamily Boinae are found in Madagascar, the Papuan-Pacific Islands, and the Neotropics. It has been suggested that genera within each of these particular areas do not form monophyletic groups. Further, it was proposed that the New World Boa constrictor is more closely related to boine genera in Madagascar than to boines in the Neotropics. Along with inferring the relationship of all boine genera using data from the cytochrome b gene and morphology, the placement of Boa was also examined. Phylogenetic inferences using maximum likelihood and Bayesian (BI) methods for combined data analyses and separate analyses of DNA sequence and morphological data were conducted. Priors, parametric bootstraps, and the Shimodaira-Hasegawa test were used to examine the previously proposed placement of Boa with Madagascan taxa using these DNA data. DNA data and combined data analyses strongly reject the hypothesis that Boa is more closely related to Old World genera than to other New World genera. Additionally, strong tree support suggests that all species within Madagascar, the Papuan-Pacific Islands, and the Neotropics each form a monophyletic group with respect to their geographic region.

Most phylogeographic studies have used maximum likelihood or maximum parsimony to infer phylogeny and bootstrap analysis to evaluate support for trees. Recently, Bayesian methods using Marlov chain Monte Carlo to search tree space and simultaneously estimate tree support have become popular due to its fast search speed and ability to create a posterior distribution of parameters of interest. Here, I present a study that utilizes Bayesian methods to infer phylogenetic relationships of the cornsnake (Elaphe guttata) complex using cytochrome b sequences. Examination of the posterior probability distributions confirms the existence of three geographic lineages. Additionally, there is no support for the monophyly of the subspecies of E. guttata. Results suggest the three geographic lineages partially conform to the ranges of previously defined subspecies, although Shimodaira-Hasegawa tests suggest that subspecies-constrained trees produce significantly poorer likelihood estimates than the most likely trees reflecting the evolution of three geographic assemblages. Based on molecular support, these three geographic assemblages are recognized as species using evolutionary species criteria: E. guttata, Elaphe slowinskii, and Elaphe emoryi [phylogeographic, maximum likelihood, maximum parsimony, bootstrap, Bayesian, Markov chain Monte Carlo, cornsnake, Cytochrome b, geographic lineages, E. guttta, E. slowinskii, and E. emoryi].