Determination of biogeographical range: an application of molecular phylogeography to the European pool frog Rana lessonae

Understanding how species are constrained within their biogeographical ranges is a central problem in evolutionary ecology. Essential prerequisites for addressing this question include accurate determinations of range borders and of the genetic structures of component populations. Human translocation of organisms to sites outside their natural range is one factor that increasingly complicates this issue. In areas not far beyond presumed natural range margins it can be particularly difficult to determine whether a species is native or has been introduced. The pool frog (Rana lessonae) in Britain is a specific example of this dilemma . We used variation at six polymorphic microsatellite loci for investigating the phylogeography of R. lessonae and establishing the affinities of specimens from British populations. The existence and distribution of a distinct northern clade of this species in Norway, Sweden and England infer that it is probably a long-standing native of Britain, which should therefore be included within its natural range. This conclusion was further supported by posterior probability estimates using Bayesian clustering. The phylogeographical analysis revealed unexpected patterns of genetic differentiation across the range of R. lessonae that highlighted the importance of historical colonization events in range structuring.     

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.     

Microsatellite analysis of the phylogeography, Pleistocene history and secondary contact hypotheses for the killifish, Fundulus heteroclitus

The mummichog, Fundulus heteroclitus, exhibits extensive latitudinal clinal variation in a number of physiological and biochemical traits, coupled with phylogeographical patterns at mitochondrial and nuclear DNA loci that suggest a complicated history of spatially variable selection and secondary intergradation. This species continues to serve as a model for understanding local and regional adaptation to variable environments. Resolving the influences of historical processes on the distribution of genetic variation within and among extant populations of F. heteroclitus is crucial to a better understanding of how populations evolve in the context of contemporary environments. In this study, we analysed geographical patterns of genetic variation at eight microsatellite loci among 15 populations of F. heteroclitus distributed throughout the North American range of the species from Nova Scotia to Georgia. Genetic variation in Northern populations was lower than in Southern populations and was strongly correlated with latitude throughout the species range. The most common Northern alleles at all eight loci exhibited concordant latitudinal clinal patterns, and the existence of an abrupt transition zone in allele frequencies between Northern and Southern populations was similar to that observed for mitochondrial DNA and allozyme loci. A significant pattern of isolation by distance was observed both within and between northern and southern regions. This pattern was unexpected, particularly for northern populations, given the recent colonization history of post-Pleistocene habitats, and was inconsistent with either a recent northward population expansion or a geographically restricted northern Pleistocene refugium. The data provided no evidence for recent population bottlenecks, and estimates of historical effective population sizes suggest that post-Pleistocene populations have been large throughout the species distribution. These results suggest that F. heteroclitus was broadly distributed throughout most of its current range during the last glacial event and that the abrupt transition in allele frequencies that separate Northern and Southern populations may reflect regional disequilibrium conditions associated with the post-Pleistocene colonization history of habitats in that region.     

Eco‐evolution on the edge during climate change

We urgently need to predict species responses to climate change to minimize future biodiversity loss and ensure we do not waste limited resources on ineffective conservation strategies. Currently, most predictions of species responses to climate change ignore the potential for evolution. However, evolution can alter species ecological responses, and different aspects of evolution and ecology can interact to produce complex eco‐evolutionary dynamics under climate change. Here we review how evolution could alter ecological responses to climate change on species warm and cool range margins, where evolution could be especially important. We discuss different aspects of evolution in isolation, and then synthesize results to consider how multiple evolutionary processes might interact and affect conservation strategies. On species cool range margins, the evolution of dispersal could increase range expansion rates and allow species to adapt to novel conditions in their new range. However, low genetic variation and genetic drift in small range‐front populations could also slow or halt range expansions. Together, these eco‐evolutionary effects could cause a three‐step, stop‐and‐go expansion pattern for many species. On warm range margins, isolation among populations could maintain high genetic variation that facilitates evolution to novel climates and allows species to persist longer than expected without evolution. This ‘evolutionary extinction debt’ could then prevent other species from shifting their ranges. However, as climate change increases isolation among populations, increasing dispersal mortality could select for decreased dispersal and cause rapid range contractions. Some of these eco‐evolutionary dynamics could explain why many species are not responding to climate change as predicted. We conclude by suggesting that resurveying historical studies that measured trait frequencies, the strength of selection, or heritabilities could be an efficient way to increase our eco‐evolutionary knowledge in climate change biology.     

Exploring demographic, physical, and historical explanations for the genetic structure of two lineages of Greater Antillean bats

Observed patterns of genetic structure result from the interactions of demographic, physical, and historical influences on gene flow. The particular strength of various factors in governing gene flow, however, may differ between species in biologically relevant ways. We investigated the role of demographic factors (population size and sex-biased dispersal) and physical features (geographic distance, island size and climatological winds) on patterns of genetic structure and gene flow for two lineages of Greater Antillean bats. We used microsatellite genetic data to estimate demographic characteristics, infer population genetic structure, and estimate gene flow among island populations of Erophylla sezekorni/E. bombifrons and Macrotus waterhousii (Chiroptera: Phyllostomidae). Using a landscape genetics approach, we asked if geographic distance, island size, or climatological winds mediate historical gene flow in this system. Samples from 13 islands spanning Erophylla's range clustered into five genetically distinct populations. Samples of M. waterhousii from eight islands represented eight genetically distinct populations. While we found evidence that a majority of historical gene flow between genetic populations was asymmetric for both lineages, we were not able to entirely rule out incomplete lineage sorting in generating this pattern. We found no evidence of contemporary gene flow except between two genetic populations of Erophylla. Both lineages exhibited significant isolation by geographic distance. Patterns of genetic structure and gene flow, however, were not explained by differences in relative effective population sizes, island area, sex-biased dispersal (tested only for Erophylla), or surface-level climatological winds. Gene flow among islands appears to be highly restricted, particularly for M. waterhousii, and we suggest that this species deserves increased taxonomic attention and conservation concern.     

Admixture and the organization of genetic diversity in a butterfly species complex revealed through common and rare genetic variants

Detailed information about the geographic distribution of genetic and genomic variation is necessary to better understand the organization and structure of biological diversity. In particular, spatial isolation within species and hybridization between them can blur species boundaries and create evolutionary relationships that are inconsistent with a strictly bifurcating tree model. Here, we analyse genome‐wide DNA sequence and genetic ancestry variation in Lycaeides butterflies to quantify the effects of admixture and spatial isolation on how biological diversity is organized in this group. We document geographically widespread and pervasive historical admixture, with more restricted recent hybridization. This includes evidence supporting previously known and unknown instances of admixture. The genome composition of admixed individuals varies much more among than within populations, and tree‐ and genetic ancestry‐based analyses indicate that multiple distinct admixed lineages or populations exist. We find that most genetic variants in Lycaeides are rare (minor allele frequency <0.5%). Because the spatial and taxonomic distributions of alleles reflect demographic and selective processes since mutation, rare alleles, which are presumably younger than common alleles, were spatially and taxonomically restricted compared with common variants. Thus, we show patterns of genetic variation in this group are multifaceted, and we argue that this complexity challenges simplistic notions concerning the organization of biological diversity into discrete, easily delineated and hierarchically structured entities.     

Northern range expansion of European populations of the wasp spider Argiope bruennichi is associated with global warming–correlated genetic admixture and population‐specific temperature adaptations

Poleward range expansions are observed for an increasing number of species, which may be an effect of global warming during the past decades. However, it is still not clear in how far these expansions reflect simple geographical shifts of species ranges, or whether new genetic adaptations play a role as well. Here, we analyse the expansion of the wasp spider Argiope bruennichi into Northern Europe during the last century. We have used a range‐wide sampling of contemporary populations and historical specimens from museums to trace the phylogeography and genetic changes associated with the range shift. Based on the analysis of mitochondrial, microsatellite and SNP markers, we observe a higher level of genetic diversity in the expanding populations, apparently due to admixture of formerly isolated lineages. Using reciprocal transplant experiments for testing overwintering tolerance, as well as temperature preference and tolerance tests in the laboratory, we find that the invading spiders have possibly shifted their temperature niche. This may be a key adaptation for survival in Northern latitudes. The museum samples allow a reconstruction of the invasion's genetic history. A first, small‐scale range shift started around 1930, in parallel with the onset of global warming. A more massive invasion of Northern Europe associated with genetic admixture and morphological changes occurred in later decades. We suggest that the latter range expansion into far Northern latitudes may be a consequence of the admixture that provided the genetic material for adaptations to new environmental regimes. Hence, global warming could have facilitated the initial admixture of populations and this resulted in genetic lineages with new habitat preferences.     

Population structure and loss of genetic diversity in the endangered white-headed duck, <Emphasis Type="Italic">Oxyura leucocephala</Emphasis>

The white-headed duck is a globally threatened species native to the Palaearctic with a range extending from Spain in the west to the western edge of China in the east. Its populations have become fragmented and undergone major declines in recent decades. To study genetic differences between populations across the range and change in genetic diversity over time, we sequenced a portion of the mitochondrial DNA control region from 67 museum specimens (years 1861–1976) as well as 39 contemporary samples from Spain and seven from Greece (years 1992–2003). In the historical sample, we found a lack of significant genetic structure between populations in different areas. We found evidence that the species experienced a rapid expansion in the past, perhaps from glacial refugia centred around the Mediterranean following the last ice age. In Spain, the population went through a dramatic bottleneck in the 1970s and early 1980s, when only a few dozens individuals remained in the wild. Although population size has since recovered to a few thousand individuals, we found a highly significant loss of mitochondrial haplotype diversity between the historical and contemporary samples. Given ongoing declines in other areas, losses in genetic diversity that may reduce the adaptive potential of white-headed ducks in the future are a continuing concern throughout the geographic range of this species.     

Linking genetic and ecological differentiation in an ungulate with a circumpolar distribution

Genetic differentiation among populations may arise from the disruption of gene flow due to local adaptation to distinct environments and/or neutral accumulation of mutations and genetic drift resulted from geographical isolation. Quantifying the role of these processes in determining the genetic structure of natural populations remains challenging. Here, we analyze the relative contribution of isolation‐by‐resistance (IBR), isolation‐by‐environment (IBE), genetic drift and historical isolation in allopatry during Pleistocene glacial cycles on shaping patterns of genetic differentiation in caribou/reindeer populations Rangifer tarandus across the entire distribution range of the species. Our study integrates analyses at range‐wide and regional scales to partial out the effects of historical and contemporary isolation mechanisms. At the circumpolar scale, our results indicate that genetic differentiation is predominantly explained by IBR and historical isolation. At a regional scale, we found that IBR, IBE and population size significantly explained the spatial distribution of genetic variation among populations belonging to the Euro‐Beringian lineage within North America. In contrast, genetic differentiation among populations within the North American lineage was predominantly explained by IBR and population size, but not IBE. We also found discrepancies between genetic and ecotype designation across the Holarctic species distribution range. Overall, these results indicate that multiple isolating mechanisms have played roles in shaping the spatial distribution of genetic variation across the distribution range of a large mammal with high potential for gene flow. Considering multiple spatial scales and simultaneously testing a comprehensive suite of potential isolating mechanisms, our study contributes to understand the ecological and evolutionary processes underlying organism–landscape interactions.     

Testing the role of interspecific competition in the evolutionary origin of elevational zonation: an example with Buarremon brush-finches (Aves, Emberizidae) in the neotropical mountains

Interspecific competition might drive the evolution of ecological niches and result in pairs of formerly competing species segregating along ecological gradients following a process of character displacement. This mechanism has been proposed to account for replacement of related species along gradients of elevation in many areas of the world, but the fundamental issue of whether competition is responsible for the origin of elevational replacements has not been tested. To test hypotheses about the role of interspecific competition in the origin of complementary elevational ranges, I combined molecular phylogenetics, phylogeography, and population genetic analyses on Buarremon torquatus and B. brunneinucha (Aves, Emberizidae), whose patterns of elevational distribution suggest character displacement or ecological release. The hypothesis that elevational distributions in these species changed in opposite directions as a result of competition is untenable because: (1) a historical expansion of the range of B. brunneinucha into areas occupied by B. torquatus was not accompanied by a shift in the elevational range of the former species; (2) when B. brunneinucha colonized the range of B. torquatus, lineages of the latter distributions had already diverged; and (3) historical trends in effective population size do not suggest populations with elevational ranges abutting those of putative competitors have declined as would be expected if competition caused range contractions. However, owing to uncertainty in coalescent estimates of historical population sizes, the hypothesis that some populations of B. torquatus have declined cannot be confidently rejected, which suggests asymmetric character displacement might have occurred. I suggest that the main role of competition in elevational zonation may be to act as a sorting mechanism that allows the coexistence along mountain slopes only of ecologically similar species that differ in elevational distributions prior to attaining sympatry. The contrasting biogeographic histories of B. brunneinucha and B. torquatus illustrate how present-day ecological interactions can have recent origins, and highlights important challenges for testing the hypothesis of character displacement in the absence of data on population history and robust reconstructions of the evolution of traits and geographic ranges.     

Spatial and temporal patterns of genetic diversity in an endangered Hawaiian honeycreeper,the Hawaii Akepa (<Emphasis Type="Italic">Loxops coccineus coccineus</Emphasis>)

As a result of disease, habitat destruction, and other anthropogenic factors, the Hawaii Akepa (Loxops coccineus coccineus) currently occupies <10% of its original range and exists in five disjunct populations, raising concerns about what effect such reduction and fragmentation has had on the connectivity and diversity of Akepa populations. In this study, we used both historical and contemporary samples to assess genetic diversity and structure in this endangered Hawaiian honeycreeper. We generated sequence data from two mtDNA regions (ND2, control region) and two nuclear introns for contemporary samples representing three of the five current populations. We also generated control region sequence data for museum specimens collected over 100 years ago from throughout the historical range of the bird. Results indicate that despite recent declines and fragmentation, genetic diversity has not been lost. We detected a modest level of genetic differentiation, which followed a combined pattern of isolation-by-barriers and isolation-by-distance, across the historical range of Akepa. The similarly low level of differentiation observed between the contemporary populations indicates that not much divergence, if any, has occurred post-fragmentation. Rather, the present structure seen likely reflects the historical pattern of distribution. Ironically, this declining species exhibits the genetic signal of an expanding population, demonstrating that earlier demographic events are outweighing the effects of recent changes in population size, and genetic estimates of N_e, though crude, suggest Hawaii Akepa were at least an order of magnitude more abundant prior to the decline.     

Genetic variation in the vulnerable and endemic Monkey Puzzle tree, detected using RAPDs

Araucaria araucana (Monkey Puzzle), a southern South American tree species of exceptional cultural and economic importance, is of conservation concern owing to extensive historical clearance and current human pressures. Random amplified polymorphic DNA (RAPD) markers were used to characterise genetic heterogeneity within and among 13 populations of this species from throughout its natural range. Extensive genetic variability was detected and partitioned by analysis of molecular variance, with the majority of variation existing within populations (87.2%), but significant differentiation was recorded among populations (12.8%). Estimates of Shannon's genetic diversity and percent polymorphism were relatively high for all populations and provide no evidence for a major reduction in genetic diversity from historical events, such as glaciation. All pairwise genetic distance values derived from analysis of molecular variance (Phi(ST)) were significant when individual pairs of populations were compared. Although populations are geographically divided into Chilean Coastal, Chilean Andes and Argentinean regions, this grouping explained only 1.77% of the total variation. Within Andean groups there was evidence of a trend of genetic distance with increasing latitude, and clustering of populations across the Andes, suggesting postglacial migration routes from multiple refugia. Implications of these results for the conservation and use of the genetic resource of this species are discussed.     

Genetic subdivision, glacial refugia and postglacial recolonization in the golden-striped salamander, Chioglossa lusitanica (Amphibia: urodela)

The golden-striped salamander (Chioglossa lusitanica) is an ecologically specialized species, endemic to north-western Iberia. Patterns of genetic variation were assessed at seven polymorphic enzyme loci and one mitochondrial DNA (mtDNA) marker (cytochrome b) in 17 populations across its range. Estimates of enzyme genetic diversity revealed a high degree of genetic subdivision (FST = 0.68), mainly attributable to the existence of two groups of populations. The groups were located, respectively, north and south of the Mondego River, indicating that this river coincided with a major historical barrier to gene flow. A significant decrease in genetic variability from the Mondego northwards was associated with the Douro and Minho rivers. mtDNA sequence variation revealed a congruent pattern of two haplotype groups (d = 2.2%), with a geographical distribution resembling that of allozymes. The pattern and depth of genetic variation is consistent with the following hypotheses: (i) subdivision of an ancestral range of the species prior to the middle Pleistocene; (ii) secondary contact between populations representing historical refugia; (iii) relatively recent range expansion giving rise to the northern part of the species range; and (iv) loss of genetic variation through founder effects during range expansion across major rivers.     

Genetic diversity loss in a biodiversity hotspot: ancient DNA quantifies genetic decline and former connectivity in a critically endangered marsupial

The extent of genetic diversity loss and former connectivity between fragmented populations are often unknown factors when studying endangered species. While genetic techniques are commonly applied in extant populations to assess temporal and spatial demographic changes, it is no substitute for directly measuring past diversity using ancient DNA (aDNA). We analysed both mitochondrial DNA (mtDNA) and nuclear microsatellite loci from 64 historical fossil and skin samples of the critically endangered Western Australian woylie (Bettongia penicillata ogilbyi), and compared them with 231 (n = 152 for mtDNA) modern samples. In modern woylie populations 15 mitochondrial control region (CR) haplotypes were identified. Interestingly, mtDNA CR data from only 29 historical samples demonstrated 15 previously unknown haplotypes and detected an extinct divergent clade. Through modelling, we estimated the loss of CR mtDNA diversity to be between 46% and 91% and estimated this to have occurred in the past 2000–4000 years in association with a dramatic population decline. In addition, we obtained near‐complete 11‐loci microsatellite profiles from 21 historical samples. In agreement with the mtDNA data, a number of ‘new’ microsatellite alleles was only detected in the historical populations despite extensive modern sampling, indicating a nuclear genetic diversity loss >20%. Calculations of genetic diversity (heterozygosity and allelic rarefaction) showed that these were significantly higher in the past and that there was a high degree of gene flow across the woylie's historical range. These findings have an immediate impact on how the extant populations are managed and we recommend the implementation of an assisted migration programme to prevent further loss of genetic diversity. Our study demonstrates the value of integrating aDNA data into current‐day conservation strategies.     

Can physiographic regions substitute for genetically-determined conservation units? A case study with the threatened plant, <Emphasis Type="Italic">Silene spaldingii</Emphasis>

Protecting genetic diversity throughout the range of a species is important for conservation, as doing so provides for long-term evolutionary potential and persistence under a changing environment. Conservation of diversity at the intraspecific level requires identification of all genetically distinct population segments within species; i.e., conservation units (CUs). Silene spaldingii occurs in grasslands of the Columbia Plateau region of western North America and is listed as threatened under the Federal Endangered Species Act. The recovery plan identified five physiographic regions across the range of the species to use as surrogates for genetic CUs. We collected leaf samples from an average of 26 plants from each of 19 of the largest populations across all five physiographic regions and used variable microsatellite and chloroplast DNA markers to determine how genetic variation is distributed across the range of the species and how well physiographic regions reflect population structure within this species. Results of several multivariate analyses clustered our samples into four genetic groups which did not correspond well with the physiographic regions. We observed little genetic differentiation among populations in the main range of the species which encompasses nearly all of four contiguous physiographic regions. However, three other distinct genetic groups were identified: two in the disjunct northeast corner and one at the southeast edge of the main range. Modification of the CUs to reflect the genetic groups rather than the physiographic regions would result in CUs which better reflect historical patterns of population structure. Moreover, use of the genetic units to inform translocation and genetic rescue efforts could improve our ability to mimic natural patterns of gene flow. Our results suggest that physiographic regions may not always be an accurate reflection of population structure for threatened or endangered species.     

Isolation by distance in the spore-forming soil bacterium Myxococcus xanthus

Genetic differentiation between spatially separated populations within a species is commonly observed in plants and animals, but its existence in microbes has long been a contentious issue. Traditionally, many microbial ecologists have reasoned that microbes are not limited by dispersal as a result of their immense numbers and microscopic size. In this view, the absence of barriers to gene flow between populations would prevent differentiation of populations by genetic drift and hinder local adaptation. Myxococcus xanthus is a globally distributed, spore-forming bacterium that offers a robust test for genetic differentiation among populations because sporulation is expected to enhance dispersal. Using multi-locus sequence data, we show here that both diversity and the degree of differentiation between populations increase as a function of distance in M. xanthus. Populations are consistently differentiated at scales exceeding 10(2)-10(3) km, and isolation by distance, the divergence of populations by genetic drift due to limited dispersal, is responsible. Our results provide new insights into how genetic diversity within species of free-living microbes is distributed from centimeter to global scales.     

MHC variation reflects the bottleneck histories of New Zealand passerines

Most empirical evidence suggests that balancing selection does not counter the effects of genetic drift in shaping postbottleneck major histocompatibility complex (MHC) genetic diversity when population declines are severe or prolonged. However, few studies have been able to include data from historical specimens, or to compare populations/species with different bottleneck histories. In this study, we examined MHC class II B and microsatellite diversity in four New Zealand passerine (songbird) species that experienced moderate to very severe declines. We compared diversity from historical samples (collected c. 1884–1938) to present‐day populations. Using a Bayesian framework, we found that the change in genetic diversity from historical to contemporary samples was affected by three main factors: (i) whether the data were based on MHC or microsatellite markers, (ii) species (as a surrogate for bottleneck severity) and (iii) whether the comparison between historical and contemporary samples was made using historical samples originating from the mainland, or using historical samples originating from islands. The greatest losses in genetic diversity occurred for the most severely bottlenecked species, particularly between historical mainland and contemporary samples. Additionally, where loss of diversity occurred, the change was greater for MHC genes compared to microsatellite loci.     

Comparative genetic analyses of historical and contemporary collections highlight contrasting demographic histories for the bumble bees Bombus pensylvanicus and B. impatiens in Illinois

Direct comparison of genetic patterns between museum specimens and contemporary collections can be a powerful approach for detecting recent demographic changes. Using microsatellite markers, we examined historical and contemporary genetic variation from an apparently declining bumble bee species, Bombus pensylvanicus , and from a stable species, Bombus impatiens , in central Illinois. For each species, we genotyped specimens from the Illinois Natural History Survey collected from three populations between 1969–1972 and from a resurvey of the same areas conducted in 2008. Population structure in B . pensylvanicus increased markedly over the last four decades (from θ_ST = 0.001 to 0.027) while no structure was detected in B . impatiens for either time period (θ_ST = –0.006 to –0.003). Changes in genetic diversity were not significant for either species, although small reductions were observed for B . pensylvanicus in all three populations. Coalescent simulations incorporating both contemporary and historical samples suggest that this small change is not surprising for recent population declines, as large reductions in genetic diversity were only apparent under the most severe bottleneck scenarios. These results demonstrate how comparisons of genetic patterns between temporal periods and species can help elucidate potential threats to population health and suggest several strategies that might be useful in the conservation of B . pensylvanicus in the Midwestern USA.     

Evolutionary history of the northern leopard frog: reconstruction of phylogeny, phylogeography, and historical changes in population demography from mitochondrial DNA

This study uses a combined methodological approach including phylogenetic, phylogeographic, and demographic analyses to understand the evolutionary history of the northern leopard frog, Rana pipiens. We tested hypotheses concerning how (or if) known geological events and key features of the species biology influenced the contemporary geographic and genetic distribution of R. pipiens. We assayed mitochondrial DNA variation from 389 individuals within 35 populations located throughout the species range. Our a priori expectations for patterns and processes influencing the current genetic structure of R. pipiens were supported by the data. However, our analyses revealed specific aspects of R. pipiens evolutionary history that were unexpected. The phylogenetic analysis indicated that R. pipiens is split into populations containing discrete eastern or western haplotypes, with the Mississippi River and Great Lakes region dividing the geographic ranges. Nested clade analysis indicated that the biological process most often invoked to explain the pattern of haplotype position is restricted gene flow with isolation by distance. Demographic analyses showed evidence of both historical bottlenecks and population expansions. Surprisingly, the genetic evidence indicated that the western haplotypes had significantly reduced levels of genetic diversity relative to the eastern haplotypes and that major range expansions occurred in both regions well before the most recent glacial retreat. This study provides a detailed history of how a widespread terrestrial vertebrate responded to episodic Pleistocene glacial events in North America. Moreover, this study illustrates how complementary methods of data analysis can be used to disentangle recent and ancient effects on the genetic structure of a species.