Phylogeographic diversification of antelope squirrels (Ammospermophilus) across North American deserts

We investigated the biogeographic history of antelope squirrels, genus Ammospermophilus, which are widely distributed across the deserts and other arid lands of western North America. We combined range‐wide sampling of all currently recognized species of Ammospermophilus with a multilocus data set to infer phylogenetic relationships. We then estimated divergence times within identified clades of Ammospermophilus using fossil‐calibrated and rate‐calibrated molecular clocks. Lastly, we explored generalized distributional changes of Ammospermophilus since the last glacial maximum using species distribution models, and assessed responses to Quaternary climate change by generating demographic parameter estimates for the three wide‐ranging clades of A. leucurus. From our phylogenetic estimates we inferred strong phylogeographic structure within Ammospermophilus and the presence of three well‐supported major clades. Initial patterns of historical divergence were coincident with dynamic alterations in the landscape of western North America, and the formation of regional deserts during the Late Miocene and Pliocene. Species distribution models and demographic parameter estimates revealed patterns of recent population expansion in response to glacial retreat. When combined with evidence from co‐distributed taxa, the historical biogeography of Ammospermophilus provides additional insight into the mechanisms that impacted diversification of arid‐adapted taxa across the arid lands of western North America. We propose species recognition of populations of the southern Baja California peninsula to best represent our current understanding of evolutionary relationships among genetic units of Ammospermophilus. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 109 , 949–967.     

Genetic differentiation and demographic history in Drosophila pachea from the Sonoran Desert

Genetic variation at six microsatellite DNA loci and a segment of the mitochondrial cytochrome oxidase subunit I (COI) locus was used to estimate gene flow, population structure, and demographic history in the cactophilic Drosophila pachea from the Sonoran Desert of North America, a species that shows a strict association with its senita host cactus (genus Lophocereus). For microsatellite analyses, thirteen populations of D. pachea were sampled, five in mainland Mexico and the southwestern USA, and eight on the Baja California (Baja) peninsula, covering essentially the entire range of the species. Analysis of molecular variance (AMOVA) of microsatellite data revealed that populations from both the mainland and the Baja peninsula generally showed little structure, although there were a few exceptions, suggesting some local differentiation and restriction of gene flow within both regions. Pairwise comparisons of F(ST) among each of the mainland and Baja populations showed evidence of both panmixia and population subdivision. AMOVA performed on grouped populations from both the mainland and Baja, however, revealed significant partitioning of genetic variation among the two regions, but no partitioning among localities within each region. Bayesian skyline analyses of the COI data set, consisting of four mainland and seven peninsular populations, revealed population expansions dating to the Pleistocene or late Pliocene in D. pachea from both regions, although regional differences were seen in the estimated timing of the expansions and in changes in effective population size over time.     

Geographic genetic architecture of pocket gopher (Thomomys bottae) populations in Baja California, Mexico

Phylogenetic analyses of complete mitochondrial cytochrome b sequences support the monophyly of pocket gopher (Thomomys bottae) populations from the 1000 km length of the Baja California peninsula of Mexico, relative to other geographical segments of the species range in western North America. The Baja California peninsula is an area that encompasses considerable ecomorphological and infraspecific diversity within this pocket gopher species. However, detailed population analyses encompassing 35 localities distributed over the southern half of the peninsula reveal only trivial phylogeographical structure. Rather, most of the 72 unique 500-base pair haplotypes examined from 142 individuals is restricted to single populations, although a few haplotypes are shared broadly across geography. Individual populations are typically comprised of haplotype sets from different branches in a network of relationships. Analysis of molecular variance (amova) indicates that approximately half of the total pool of variation is contained among individuals within local populations, and that only about 25% can be explained by the regional subdivisions of current subspecies distributions or physiographic realms. A hypothesized historical vicariant event that has been causally linked to the phylogeographical structure of other, codistributed species has had little influence on these pocket gopher populations, explaining only 13% of the total variation. The temporal depth, estimated by coalescence parameters, of the haplotype lineage in Baja California is relatively recent, approximately 300,000 generations; both the mismatch distribution of pairwise comparisons and a significantly positive exponential growth estimate support a recent history of expanding populations; but current, or recent past, migration estimates have remained small, are largely unidirectional from north to south, and weak isolation by distance is present. All data suggest that pocket gophers have relatively recently invaded the southern half of peninsular Baja California, with the genetic signature of expansion still evident but with sufficient time having lapsed to result in a weak isolation by distance pattern. The geographical assemblage of sampled populations thus appears as a meta-population, with limited gene flow contrasting with random haplotype loss due to drift in small, localized populations.     

Phylogeography of the mountain chickadee (Poecile gambeli): diversification, introgression, and expansion in response to Quaternary climate change

Since the late 1990s, molecular techniques have fuelled debate about the role of Pleistocene glacial cycles in structuring contemporary avian diversity in North America. The debate is still heated; however, there is widespread agreement that the Pleistocene glacial cycles forced the repeated contraction, fragmentation, and expansion of the North American biota. These demographic processes should leave genetic 'footprints' in modern descendants, suggesting that detailed population genetic studies of contemporary species provide the key to elucidating the impact of the late Quaternary (late Pleistocene-Holocene). We present an analysis of mitochondrial DNA (mtDNA) variation in the mountain chickadee (Poecile gambeli) in an attempt to examine the genetic evidence of the impact of the late Quaternary glacial cycles. Phylogenetic analyses reveal two strongly supported clades of P. gambeli: an Eastern Clade (Rocky Mountains and Great Basin) and a Western Clade (Sierra Nevada and Cascades). Post-glacial introgression is apparent between these two clades in the Mono Lake region of Central California. Within the Eastern Clade there is evidence of isolation-by-distance in the Rocky Mountain populations, and of limited gene flow into and around the Great Basin. Coalescent analysis of genetic variation in the Western Clade indicates that northern (Sierra Nevada/Cascades) and southern (Transverse/Peninsular Ranges) populations have been isolated and evolving independently for nearly 60,000 years.     

Beringian origins and cryptic speciation events in the fly agaric (Amanita muscaria)

Amanita muscaria sensu lato has a wide geographic distribution, occurring in Europe, Asia, Africa, Australia, New Zealand, and North, Central and South America. Previous phylogenetic work by others indicates three geographic clades (i.e. 'Eurasian', 'Eurasian-alpine' and 'North American' groups) within A. muscaria. However, the historical dispersal patterns of A. muscaria remained unclear. In our project, we collected specimens from arctic, boreal and humid temperate regions in Alaska, and generated DNA sequence data from the protein-coding beta-tubulin gene and the internal transcribed spacer (ITS) and large subunit (LSU) regions of the ribosomal DNA repeat. Homologous sequences from additional A. muscaria isolates were downloaded from GenBank. We conducted phylogenetic and nested clade analyses (NCA) to reveal the phylogeographic history of the species complex. Although phylogenetic analyses confirmed the existence of the three above-mentioned clades, representatives of all three groups were found to occur sympatrically in Alaska, suggesting that they represent cryptic phylogenetic species with partially overlapping geographic distributions rather than being allopatric populations. All phylogenetic species share at least two morphological varieties with other species, suggesting ancestral polymorphism in pileus and wart colour pre-dating their speciations. The ancestral population of A. muscaria likely evolved in the Siberian-Beringian region and underwent fragmentation as inferred from NCA and the coalescent analyses. The data suggest that these populations later evolved into species, expanded their range in North America and Eurasia. In addition to range expansions, populations of all three species remained in Beringia and adapted to the cooling climate.     

Phylogeography of Ramalina menziesii,a widely distributed lichen‐forming fungus in western North America

The complex topography and climate history of western North America offer a setting where lineage formation, accumulation and migration have led to elevated inter‐ and intraspecific biodiversity in many taxa. Here, we study Ramalina menziesii, an epiphytic lichenized fungus with a range encompassing major ecosystems from Baja California to Alaska to explore the predictions of two hypotheses: (i) that the widespread distribution of R. menziesii is due to a single migration episode from a single lineage and (ii) that the widespread distribution is due to the formation and persistence of multiple lineages structured throughout the species' range. To obtain evidence for these predictions, we first construct a phylogenetic tree and identify multiple lineages structured throughout the species' range – some ancient ones that are localized and other more recent lineages that are widely distributed. Second, we use an isolation with migration model to show that sets of ecoregion populations diverged from each other at different times, demonstrating the importance of historical and current barriers to gene flow. Third, we estimated migration rates among ecoregions and find that Baja California populations are relatively isolated, that inland California ecoregion populations do not send out emigrants and that migration out of California coastal and Pacific Northwest populations into inland California ecoregions is high. Such intraspecific geographical patterns of population persistence and dispersal both contribute to the wide range of this genetically diverse lichen fungus and provide insight into the evolutionary processes that enhance species diversity of the California Floristic Province.     

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

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

Geological barriers and restricted gene flow in the holarctic skipper Hesperia comma (Hesperiidae)

Patterns of genetic variation within a species may be a consequence of historical factors, such as past fragmentation, as well as current barriers to gene flow. Using sequence data from the mitochondrial cytochrome oxidase subunit II region (COII) and the nuclear gene wingless, we conducted a phylogeographical study of the holarctic skipper Hesperia comma to elucidate patterns of genetic diversity and to infer historical and contemporary processes maintaining genetic variation. One hundred and fifty-one individuals were sampled from throughout North America and Eurasia, focusing on California and adjacent regions in the western United States where morphological diversity is highest compared to the rest of the range. Analyses of sequence data obtained from both genes revealed a well-supported division between the Old and New World. Within western North America, wingless shows little geographical structure, while a hierarchical analysis of genetic diversity of COII sequences indicates three major clades: a western clade in Oregon and Northern California, an eastern clade including the Great Basin, Rocky Mountains and British Columbia, and a third clade in southern California. The Sierra Nevada and the Transverse Ranges appear to be the major barriers to gene flow for H. comma in the western United States. Relatively reduced haplotype diversity in Eurasia compared to North America suggests that populations on the two continents have been affected by different historical processes.     

Phylogeography and the origins of range disjunctions in a north temperate fish,the pygmy whitefish (Prosopium coulterii), inferred from mitochondrial and nuclear DNA sequence analysis

Aim To investigate the degree of phylogeographical divergence within pygmy whitefish (Prosopium coulterii) and to test hypotheses concerning the origin of disjunct populations within North America. Location North America from western Alaska to Lake Superior. Methods Mitochondrial (ATPase subunit VI) and nuclear (ITS‐1, ITS‐2) DNA sequence variation was assessed across the species’ North American range to test for the existence of distinct phylogeographical groupings of pygmy whitefish associated with known glacial refugia. Coalescent simulations of the mitochondrial DNA (mtDNA) data were used to test hypotheses of population structure. Results This species is composed of two monophyletic mitochondrial clades across its North American range. The two mtDNA clades differed by an average 3.3% nucleotide sequence divergence. These clades were also distinguished by ITS‐2, but the relationships among lineages were not resolved by the ITS‐1 analysis. Coalescent analyses rejected the null hypothesis of the current disjunct distributions being a result of fragmentation of a single widespread ancestral lineage across a variety of effective population sizes and divergence times. Main conclusions The current range disjunctions of pygmy whitefish in North America probably resulted from isolation, genetic divergence, and selective dispersal from at least two major Pleistocene glacial refugia: Beringia and Cascadia. More recent isolation and dispersal from an upper Mississippi refugium is suggested by relationships within one of the clades and by distributional evidence from co‐distributed species. The Beringian and Cascadian refugia have played major roles in the zoogeography of Nearctic temperate aquatics, but the roles of smaller refugia appear more variable among other species.     

Genetic structure of desert ground squirrels over a 20-degree-latitude transect from Oregon through the Baja California peninsula

The genetic structure of populations over a wide geographical area should reflect the demographic and evolutionary processes that have shaped a species across its range. We examined the population genetic structure of antelope ground squirrels (Ammospermophilus leucurus) across the complex of North American deserts from the Great Basin of Oregon to the cape region of the Baja California peninsula. We sampled 73 individuals from 13 major localities over this 2500-km transect, from 43 to 22 degrees north. Our molecular phylogeographical analysis of 555 bp of the mitochondrial cytochrome b gene and 510 bp of the control region revealed great genetic uniformity in a single clade that extends from Oregon to central Baja California. A second distinct clade occupies the southern half of the peninsula. The minimal geographical structure of the northern clade, its low haplotype diversity and the distribution of pairwise differences between haplotypes suggest a rapid northward expansion of the population that must have followed a northward desert habitat shift associated with the most recent Quaternary climate warming and glacial retreat. The higher haplotype diversity within the southern clade and distribution of pairwise differences between haplotypes suggest that the southern clade has a longer, more stable history associated with a southern peninsular refugium. This system, as observed, reflects both historical and contemporary ecological and evolutionary responses to physical environmental gradients within genetically homogeneous populations.     

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.     

Conservation Genetics of Kangaroo Mice, Genus Microdipodops

The two currently recognized species of kangaroo mice, Microdipodops megacephalus and M. pallidus, inhabit sandy soils of the Great Basin Desert in western North America. Given their habitat specificity and the fluctuating climate throughout the Pleistocene, kangaroo mice likely endured a turbulent biogeographic history that resulted in disjunct distributions and isolation of genetic lineages. Recent phylogenetic investigations using mitochondrial data have revealed several mitochondrial clades within this genus that may represent cryptic species. These mitochondrial clades are genetically unique, occupy relatively small distributions, and, as such, may be at an increased risk of extinction due to climate change and extensive recent habitat alteration. Herein, we apply haplotype network, population genetic, and historical demographic analyses to mitochondrial data of each Micropdipodops species and mitochondrial clade to assess conservation genetics within kangaroo mice. Results indicate that each mitochondrial clade is a distinct lineage with little to no gene flow occurring among clades. Additionally, historical demographic analyses support past population expansions and identify locations of past refugium for each distinct lineage. Although mitochondrial data indicate that the clades appear to be in approximate genetic equilibrium and have not suffered any extreme bottlenecks over time, there is still concern for the survival of smaller and more vulnerable Microdipodops subpopulations due to impending habitat threats in the Great Basin Desert.     

Dispersal Pathways and Genetic Differentiation among Worldwide Populations of the Invasive Weed Centaurea solstitialis L. (Asteraceae)

The natural history of introduced species is often unclear due to a lack of historical records. Even when historical information is readily available, important factors of the invasions such as genetic bottlenecks, hybridization, historical relationships among populations and adaptive changes are left unknown. In this study, we developed a set of nuclear, simple sequence repeat markers and used these to characterize the genetic diversity and population structure among native (Eurasian) and non-native (North and South American) populations of Centaurea solstitialis L., (yellow starthistle). We used these data to test hypotheses about the invasion pathways of the species that were based on historical and geographical records, and we make inferences about historical relationships among populations and demographic processes following invasion. We confirm that the center of diversity and the native range of the species is likely the eastern Mediterranean region in the vicinity of Turkey. From this region, the species likely proceeded to colonize other parts of Europe and Asia via a slow, stepwise range expansion. Spanish populations were the primary source of seed to invade South America via human-mediated events, as was evident from historical records, but populations from the eastern Mediterranean region were also important. North American populations were largely derived from South America, but had secondary contributors. We suggest that the introduction history of non-native populations from disparate parts of the native range have allowed not just one, but multiple opportunities first in South America then again in North America for the creation of novel genotypes via intraspecific hybridization. We propose that multiple intraspecific hybridization events may have created especially potent conditions for the selection of a noxious invader, and may explain differences in genetic patterns among North and South America populations, inferred differences in demographic processes, as well as morphological differences previously reported from common garden experiments.     

Model-based comparisons of phylogeographic scenarios resolve the intraspecific divergence of cactophilic Drosophila mojavensis

The cactophilic fly Drosophila mojavensis exhibits considerable intraspecific genetic structure across allopatric geographic regions and shows associations with different host cactus species across its range. The divergence between these populations has been studied for more than 60years, yet their exact historical relationships have not been resolved. We analysed sequence data from 15 intronic X-linked loci across populations from Baja California, mainland Sonora-Arizona and Mojave Desert regions under an isolation-with-migration model to assess multiple scenarios of divergence. We also compared the results with a pre-existing sequence data set of eight autosomal loci. We derived a population tree with Baja California placed at its base and link their isolation to Pleistocene climatic oscillations. Our estimates suggest the Baja California population diverged from an ancestral Mojave Desert/mainland Sonora-Arizona group around 230,000-270,000years ago, while the split between the Mojave Desert and mainland Sonora-Arizona populations occurred one glacial cycle later, 117,000-135,000years ago. Although we found these three populations to be effectively allopatric, model ranking could not rule out the possibility of a low level of gene flow between two of them. Finally, the Mojave Desert population showed a small effective population size, consistent with a historical population bottleneck. We show that model-based inference from multiple loci can provide accurate information on the historical relationships of closely related groups allowing us to set into historical context a classic system of incipient ecological speciation.     

Chloroplast DNA phylogeography of a distylous shrub (Palicourea padifolia, Rubiaceae) reveals past fragmentation and demographic expansion in Mexican cloud forests

Several phylogeographic studies in northern Mesoamerica have examined the influence of Pleistocene glaciations on the genetic structure of temperate tree species with their southern limit by the contact zone between species otherwise characteristic of North or South America, but few have featured plant species that presumably colonized northern Mesoamerica from South America. A phylogeographical study of Palicourea padifolia, a fleshy-fruited, bird dispersed distylous shrub, was conducted to investigate genetic variation at two chloroplast regions (trnS-trnG and rpl32-trnL) across cloud forest areas to determine if such patterns are consistent with the presence of Pleistocene refugia and/or with the historical fragmentation of the Mexican cloud forests. We conducted population and spatial genetic analyses as well as phylogenetic and isolation with migration analyses on 122 individuals from 22 populations comprising the distribution of P. padifolia in Mexico to gain insight of the evolutionary history of these populations. Twenty-six haplotypes were identified after sequencing 1389 bp of chloroplast DNA. These haplotypes showed phylogeographic structure (N(ST) = 0.508, G(ST) = 0.337, N(ST) > G(ST), P < 0.05), including a phylogeographic break at the Isthmus of Tehuantepec, with private haplotypes at either side of the isthmus, and a divergence time of the split in the absence of gene flow dating back c. 309,000-103,000 years ago. The patterns of geographic structure found in this study are consistent with past fragmentation and demographic range expansion, supporting the role of the Isthmus of Tehuantepec as a biogeographical barrier in the dispersal of P. padifolia. Our data suggest that P. padifolia populations were isolated throughout glacial cycles by the Isthmus of Tehuantepec, accumulating genetic differences due to the lack of migration across the isthmus in either direction, but the results of our study are not consistent with the existence of the previously proposed Pleistocene refugia for rain forest plant species in the region.     

Montane refugia predict population genetic structure in the Large‐blotched Ensatina salamander

Understanding the biotic consequences of Pleistocene range shifts and fragmentation remains a fundamental goal in historical biogeography and evolutionary biology. Here, we combine species distribution models (SDM) from the present and two late Quaternary time periods with multilocus genetic data (mitochondrial DNA and microsatellites) to evaluate the effect of climate‐induced habitat shifts on population genetic structure in the Large‐blotched Ensatina (Ensatina eschscholtzii klauberi), a plethodontid salamander endemic to middle and high‐elevation conifer forest in the Transverse and Peninsular Ranges of southern California and northern Baja California. A composite SDM representing the range through time predicts two disjunct refugia, one in southern California encompassing the core of the species range and the other in the Sierra San Pedro Mártir of northern Baja California at the southern limit of the species range. Based on our spatial model, we would expect a pattern of high connectivity among populations within the northern refugium and, conversely, a pattern of isolation due to long‐term persistence of the Sierra San Pedro Mártir population. Our genetic results are consistent with these predictions based on the hypothetical refugia in that (i) historical measures of population connectivity among stable areas are correlated with gene flow estimates; and (ii) there is strong geographical structure between separate refugia. These results provide evidence for the role of recent climatic change in shaping patterns of population persistence and connectivity within the Transverse and Peninsular Ranges, an evolutionary hotspot.     

Range-wide phylogeography of a temperate lizard, the five-lined skink (Eumeces fasciatus)

We used mitochondrial DNA and microsatellite loci to examine the phylogeographic patterns of the most broadly distributed lizard in eastern North America, the five-lined skink (Eumeces fasciatus). We infer that longitudinal phylogeographic patterns in E. fasciatus are consistent with fragmentation due to refugial and post-glacial dynamics, but that deep divergences within the species imply historical fragmentation that predates the Pleistocene. The effect of multiple refugia is implied from our nested clade analyses, including a northern refugium in Wisconsin. Analysis of population structure using nuclear microsatellite data within the species suggests the importance of glacial dynamics in shaping more recent genetic structuring within one widely distributed lineage that ranges from the Mississippi River to the Atlantic Ocean in longitude and from southern Ontario to the Gulf of Mexico in latitude. Results shed light on the historical processes that have influenced current population structure of a temperate lizard, support the striking similarity of longitudinal phylogeographic structure across many herpetofaunal species in eastern North America, and illustrate the utility of employing multiple markers in phylogeographic studies.     

Historical demography and genetic structure of sister species: deermice (Peromyscus) in the North American temperate rain forest

To investigate the evolutionary and biogeographical history of Peromyscus keeni and P. maniculatus within the coastal forest ecosystem of the Pacific Northwest of North America, we sampled 128 individuals from 43 localities from southeastern Alaska through Oregon. We analysed mitochondrial DNA variation using DNA sequence data from the mitochondrial cytochrome-b (cyt-b) gene and control region, and we found two distinct clades consistent with the morphological designation of the two species. The sequence divergence between the two clades was 0.0484 substitutions per site for cyt-b and 0.0396 for the control region, suggesting that divergence of the two clades occurred during the middle to late Pleistocene. We also examined the historical demography of the two clades using stepwise and exponential expansion models, both of which indicated recent rapid population growth. Furthermore, using the program migrate we found evidence of migration from populations north of the Fraser River (British Columbia) to the south in both clades. This study demonstrates the utility of these model-based demographic methods in illuminating the evolutionary and biogegographic history of natural systems.     

Phylogeography of Diadophis punctatus: extensive lineage diversity and repeated patterns of historical demography in a trans-continental snake

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.     

Phylogeography of a widespread small carnivore,the western spotted skunk (Spilogale gracilis) reveals temporally variable signatures of isolation across western North America

We analyzed phylogeographic patterns in the western spotted skunk, Spilogale gracilis Merriam, 1890 (Carnivora: Mephitidae) in relation to historical events associated with Pre‐Pleistocene Divergence (PPD) and Quaternary climate change (QCC) using mitochondrial DNA from 97 individuals distributed across Western North America. Divergence times were generated using BEAST to estimate when isolation in putative refugia occurred. Patterns and timing of demographic expansion was performed using Bayesian skyline plot. Putative climatic refugia resulting from Quaternary climate change were identified using paleoecological niche modeling and divergence dates compared to major vicariant events associated with Pre‐Pleistocene conditions. We recovered three major mitochondrial clades corresponding to western North America (California, Baja, and across the Great Basin), east‐central North America (Texas, central Mexico, New Mexico), and southwestern Arizona/northwestern Mexico. Time to most recent common ancestor for S. gracilis occurred ~1.36 Ma. Divergence times for each major clade occurred between 0.25 and 0.12 Ma, with signature of population expansion occurring 0.15 and 0.10 Ma. Ecological niche models identified three potential climatic refugia during the Last Interglacial, (1) west coast of California and Oregon, (2) northwestern Mexico, and (3) southern Texas/northeastern Mexico as well as two refugia during the Last Glacial Maximum, (1) western USA and (2) southern Texas/northeastern Mexico. This study supports PPD in shaping species‐level diversity compared to QCC‐driven changes at the intraspecific level for Spilogale, similar to the patterns reported for other small mammals (e.g., rodents and bats). Phylogeographic patterns also appear to have been shaped by both habitat and river vicariance, especially across the desert southwest. Further, continuing climate change during the Holocene coupled with anthropogenic modifications during the Anthropocene appears to be removing both of these barriers to current dispersal of western spotted skunks.