Phylogeographic and population genetic analyses reveal Pleistocene isolation followed by high gene flow in a wide ranging,but endangered,freshwater mussel

Freshwater organisms of North America have had their contemporary genetic structure shaped by vicariant events, especially Pleistocene glaciations. Life history traits promoting dispersal and gene flow continue to shape population genetic structure. Cumberlandia monodonta, a widespread but imperiled (IUCN listed as endangered) freshwater mussel, was examined to determine genetic diversity and population genetic structure throughout its range. Mitochondrial DNA sequences and microsatellite loci were used to measure genetic diversity and simulate demographic events during the Pleistocene using approximate Bayesian computation (ABC) to test explicit hypotheses explaining the evolutionary history of current populations. A phylogeny and molecular clock suggested past isolation created two mtDNA lineages during the Pleistocene that are now widespread. Two distinct groups were also detected with microsatellites. ABC simulations indicated the presence of two glacial refugia and post-glacial admixture of them followed by simultaneous dispersal throughout the current range of the species. The Ouachita population is distinct from others and has the lowest genetic diversity, indicating that this is a peripheral population of the species. Gene flow within this species has maintained high levels of genetic diversity in most populations; however, all populations have experienced fragmentation. Extirpation from the center of its range likely has isolated remaining populations due to the geographic distances among them.     

Identifying glacial refugia in a geographic parthenogen using palaeoclimate modelling and phylogeography: the New Zealand stick insect Argosarchus horridus (White)

We have used phylogeographic analysis of mitochondrial DNA (COI and COII genes) and ecological niche modelling (ENM) to reconstruct the population history of Argosarchus horridus (White), a widespread species of New Zealand stick insect. These data were used to address outstanding questions on the role of glacial refugia in determining the distribution and genetic structure of New Zealand species. Phylogeographic analysis shows a general pattern of high diversity in upper North Island and reduced diversity in lower North Island and South Island. The ENM indicates that during the last glacial maximum, A. horridus was largely restricted to refugia around coastal areas of North Island. The ENM also suggests refugia on the northeast coast of South Island and southeast coast of North Island and this prediction is verified by phylogeographic analysis, which shows a clade restricted to this region. Argosarchus horridus is also most likely a geographic parthenogen where males are much rarer at higher latitudes. The higher levels of genetic variation in northern, bisexual populations suggest southern and largely unisexual populations originated from southwardly expanding parthenogenetic lineages. Bayesian skyline analysis also provides support for a recent population size increase consistent with a large increase in geographic distribution in the late Pleistocene. These results exemplify the utility of integrating ENM and phylogeographic analysis in testing hypotheses on the origin of geographic parthenogenesis and effects of Pleistocene environmental change on biodiversity.     

Multiple Pleistocene refugia and Holocene range expansion of an abundant southwestern American desert plant species (Melampodium leucanthum,Asteraceae)

Pleistocene climatic fluctuations had major impacts on desert biota in southwestern North America. During cooler and wetter periods, drought‐adapted species were isolated into refugia, in contrast to expansion of their ranges during the massive aridification in the Holocene. Here, we use Melampodium leucanthum (Asteraceae), a species of the North American desert and semi‐desert regions, to investigate the impact of major aridification in southwestern North America on phylogeography and evolution in a widespread and abundant drought‐adapted plant species. The evidence for three separate Pleistocene refugia at different time levels suggests that this species responded to the Quaternary climatic oscillations in a cyclic manner. In the Holocene, once differentiated lineages came into secondary contact and intermixed, but these range expansions did not follow the eastwardly progressing aridification, but instead occurred independently out of separate Pleistocene refugia. As found in other desert biota, the Continental Divide has acted as a major migration barrier for M. leucanthum since the Pleistocene. Despite being geographically restricted to the eastern part of the species’ distribution, autotetraploids in M. leucanthum originated multiple times and do not form a genetically cohesive group.     

Pleistocene glacial refugia across the Appalachian Mountains and coastal plain in the millipede genus Narceus : Evidence from population genetic, phylogeographic, and paleoclimatic data

Background  

Species that are widespread throughout historically glaciated and currently non-glaciated areas provide excellent opportunities to investigate the role of Pleistocene climatic change on the distribution of North American biodiversity. Many studies indicate that northern animal populations exhibit low levels of genetic diversity over geographically widespread areas whereas southern populations exhibit relatively high levels. Recently, paleoclimatic data have been combined with niche-based distribution modeling to locate possible refugia during the Last Glacial Maximum. Using phylogeographic, population, and paleoclimatic data, we show that the distribution and mitochondrial data for the millipede genus Narceus are consistent with classical examples of Pleistocene refugia and subsequent post-glacial population expansion seen in other organismal groups.     

Explicit tests of palaeodrainage connections of southeastern North America and the historical biogeography of Orangethroat Darters (Percidae: Etheostoma: Ceasia)

The alteration in palaeodrainage river connections has shaped patterns of speciation, genetic diversity and the geographical distribution of the species‐rich freshwater fauna of North America. The integration of ancestral range reconstruction methods and divergence time estimates provides an opportunity to infer palaeodrainage connectivity and test alternative palaeodrainage hypotheses. Members of the Orangethroat Darter clade, Ceasia, are endemic to southeastern North America and occur north and south of the Pleistocene glacial front, a distributional pattern that makes this clade of closely related species an ideal system to investigate the number and location of glacial refugia and compare alternative hypotheses regarding the proposed evolution of the Teays‐Mahomet palaeodrainage. This study utilized time‐calibrated mitochondrial and nuclear gene phylogenies and present‐day geographical distributions to investigate hypothesized Teays‐Mahomet River connections through time using a dispersal–extinction–cladogenesis (DEC) framework. Results of DEC ancestral area reconstructions indicate that the Teays‐Mahomet River was a key dispersal route between disjunct highland regions connecting the Mississippi River tributaries to the Old‐Ohio Drainage minimally at two separate occasions during the Pleistocene. There was a dynamic interplay between palaeodrainage connections through time and postglacial range expansion from three glacial refugia that shaped the current genetic structure and geographical distributions of the species that comprise Ceasia.     

High levels of population subdivision in a morphologically conserved Mediterranean toad (Alytes cisternasii) result from recent,multiple refugia: evidence from mtDNA,microsatellites and nuclear genealogies

Pleistocene glaciations often resulted in differentiation of taxa in southern European peninsulas, producing the high levels of endemism characteristic of these regions (e.g. the Iberian Peninsula). Despite their small ranges, endemic species often exhibit high levels of intraspecific differentiation as a result of a complex evolutionary history dominated by successive cycles of fragmentation, expansion and subsequent admixture of populations. Most evidence so far has come from the study of species with an Atlantic distribution in northwestern Iberia, and taxa restricted to Mediterranean‐type habitats remain poorly studied. The Iberian Midwife toad (Alytes cisternasii) is a morphologically conserved species endemic to southwestern and central Iberia and a typical inhabitant of Mediterranean habitats. Applying highly variable genetic markers from both mitochondrial and nuclear genomes to samples collected across the species’ range, we found evidence of high population subdivision within A. cisternasii. Mitochondrial haplotypes and microsatellites show geographically concordant patterns of genetic diversity, suggesting population fragmentation into several refugia during Pleistocene glaciations followed by subsequent events of geographical and demographic expansions with secondary contact. In addition, the absence of variation at the nuclear β‐fibint7 and Ppp3caint4 gene fragments suggests that populations of A. cisternasii have been recurrently affected by episodes of extinction and recolonization, and that documented patterns of population subdivision are the outcome of recent and multiple refugia. We discuss the evolutionary history of the species with particular interest in the increasing relevance of Mediterranean refugia for the survival of genetically differentiated populations during the Pleistocene glaciations as revealed by studies in co‐distributed taxa.     

Multilocus phylogeography and population structure of common eiders breeding in North America and Scandinavia

Aim Glacial refugia during the Pleistocene had major impacts on the levels and spatial apportionment of genetic diversity of species in northern latitude ecosystems. We characterized patterns of population subdivision, and tested hypotheses associated with locations of potential Pleistocene refugia and the relative contribution of these refugia to the post‐glacial colonization of North America and Scandinavia by common eiders (Somateria mollissima). Specifically, we evaluated localities hypothesized as ice‐free areas or glacial refugia for other Arctic vertebrates, including Beringia, the High Arctic Canadian Archipelago, Newfoundland Bank, Spitsbergen Bank and north‐west Norway. Location Alaska, Canada, Norway and Sweden. Methods Molecular data from 12 microsatellite loci, the mitochondrial DNA (mtDNA) control region, and two nuclear introns were collected and analysed for 15 populations of common eiders (n = 716) breeding throughout North America and Scandinavia. Population genetic structure, historical population fluctuations and gene flow were inferred using F‐statistics, analyses of molecular variance, and multilocus coalescent analyses. Results Significant inter‐population variation in allelic and haplotypic frequencies were observed (nuclear DNA F_ST = 0.004–0.290; mtDNA Φ_ST = 0.051–0.927). Whereas spatial differentiation in nuclear genes was concordant with subspecific designations, geographic proximity was more predictive of inter‐population variance in mitochondrial DNA haplotype frequency. Inferences of historical population demography were consistent with restriction of common eiders to four geographic areas during the Last Glacial Maximum: Belcher Islands, Newfoundland Bank, northern Alaska and Svalbard. Three of these areas coincide with previously identified glacial refugia: Newfoundland Bank, Beringia and Spitsbergen Bank. Gene‐flow and clustering analyses indicated that the Beringian refugium contributed little to common eider post‐glacial colonization of North America, whereas Canadian, Scandinavian and southern Alaskan post‐glacial colonization is likely to have occurred in a stepwise fashion from the same glacial refugium. Main conclusions Concordance of proposed glacial refugia used by common eiders and other Arctic species indicates that Arctic and subarctic refugia were important reservoirs of genetic diversity during the Pleistocene. Furthermore, suture zones identified at MacKenzie River, western Alaska/Aleutians and Scandinavia coincide with those identified for other Arctic vertebrates, suggesting that these regions were strong geographic barriers limiting dispersal from Pleistocene refugia.     

Phylogeography of a widespread North American migratory songbird (Setophaga ruticilla)

Genetic analyses for many widespread North American species have revealed significant east-west differentiation, indicating that many survived through the Pleistocene in 2 glacial refugia-1 in the eastern and 1 in the western part of the continent. It remains unclear, however, whether other areas may have served as important glacial refugia. Moreover, many such species exhibit widespread genetic similarity within eastern and western regions because of recent expansion from small refugial populations, making it difficult to evaluate current-day levels of gene flow. In this study, we used mitochondrial DNA (mtDNA) control region sequence and amplified fragment length polymorphism markers to survey genetic variation in a widespread migratory bird, the American redstart (Setophaga ruticilla). mtDNA analyses revealed a pattern that contrasts with that found for most other widespread species studied to date: most redstart populations across North America appear to have spread out from a single glacial refugium, possibly located in the southeastern United States, whereas populations in far-eastern Canada may have survived in a second glacial refugium located on the now-submerged Atlantic coastal shelf off the coast of Newfoundland. A pattern of isolation by distance in mtDNA suggested some constraints on current-day gene flow among extant redstart populations. This study thus reveals a recent evolutionary history for this species that differs from that of most other widespread North American passerines and provides evidence for limited gene flow in a species with potentially large dispersal distances.     

Increased genetic structuring of isolated Salamandra salamandra populations (Caudata: Salamandridae) at the margins of the Carpathian Mountains

The traditional southern Pleistocene refugia hypothesis in Europe has lately been challenged for several animal and plant species. The Carpathian Basin, especially at the marginal regions, is one of the recently recognized biodiversity hotspots in Europe. Marginal populations are prone to have lower genetic diversity and higher genetic differentiation than central populations. Here, we examined one mitochondrial DNA fragment (D‐loop) and nine nuclear (microsatellite) loci to describe the genetic diversity and phylogeographical pattern of fire salamander (Salamandra salamandra) populations in the Carpathian Basin with focusing on the southern margins of the Western Carpathians, where isolated populations of this species are present. Analyses of microsatellites indicated reduced genetic diversity for most of the isolated populations. Based on the mitochondrial DNA, only two haplotypes were found, whereas the analyses with the nuclear markers revealed a more recent genetic split between Western (Alpine) and Eastern (Carpathian) populations, and separated the Apuseni Mountains population (part of the Western Carpathians). Using approximate Bayesian computation analyses, we identified the most probable colonization scenario for the isolated North Hungarian Carpathian Basin populations. The split between isolated salamander populations from the central populations in the Carpathian Mountains dates back to the beginning of the Late Pleistocene, while the split between most of the Hungarian populations can be associated with the Last Glacial Maximum. We found evidence for long‐time isolation between the marginal Carpathian Basin and central populations. Our results also show that S. salamandra survived glacial periods in the temperate forests of north‐east Pannonia (North Hungarian Mountains), confirming that the Carpathian Basin served as important northerly refugia during the Pleistocene climatic oscillations.     

PHYLOGEOGRAPHIC PATTERNS INDICATE TRANSATLANTIC MIGRATION FROM EUROPE TO NORTH AMERICA IN THE RED SEAWEED CHONDRUS CRISPUS (GIGARTINALES,RHODOPHYTA)1

Inferring how the Pleistocene climate oscillations have repopulated the extant population structure of Chondrus crispus Stackh. in the North Atlantic Ocean is important both for our understanding of the glacial episode promoting diversification and for the conservation and development of marine organisms. C. crispus is an ecologically and commercially important red seaweed with broad distributions in the North Atlantic. Here, we employed both partial mtDNA Cox1 and nrDNA internal transcribed spacer region 2 (ITS2) sequences to explore the genetic structure of 17 C. crispus populations from this area. Twenty‐eight and 30 haplotypes were inferred from these two markers, respectively. Analysis of molecular variance (AMOVA) and of the population statistic Φ_ST not only revealed significant genetic structure within C. crispus populations but also detected significant levels of genetic subdivision among and within populations in the North Atlantic. On the basis of high haplotype diversity and the presence of endemic haplotypes, we postulate that C. crispus had survived in Pleistocene glacial refugia in the northeast Atlantic, such as the English Channel and the northwestern Iberian Peninsula. We also hypothesize that C. crispus from the English Channel refugium repopulated most of northeastern Europe and recolonized northeastern North America in the Late Pleistocene. The observed phylogeographic pattern of C. crispus populations is in agreement with a scenario in which severe Quaternary glaciations influenced the genetic structure of North Atlantic marine organisms with contiguous population expansion and locally restricted gene flow coupled with a transatlantic dispersal in the Late Pleistocene.     

Refugia in the Cordilleran ice sheet of western North America: chloroplast DNA diversity in the Arctic–alpine plant Oxyria digyna

Aim  Late Pleistocene ice sheets are thought to have covered most of western Canada, including all of British Columbia (BC). We examine patterns of genetic variation in an Arctic–alpine plant to evaluate the possibility of full glacial refugia within the area covered by the Cordilleran ice sheet (CIS) and to uncover post-glacial migration routes.
Location  Western North America.
Methods  We sampled 1030 individuals of the Arctic–alpine plant Oxyria digyna from 117 populations distributed over much of its range in western and northern North America. DNA haplotypes were identified using restriction site analysis of two chloroplast DNA intergene spacer regions, psb A- trn H and trn T-L. We examined the geographical distribution of haplotype diversity in relation to latitude, and also compared various indices of diversity in putatively glaciated and unglaciated regions. Patterns of migration were inferred using nested clade analysis.
Results  We detected a total of 20 haplotypes. High haplotype diversity was found in Beringia, in unglaciated western USA, and in northern BC at 57–59° N, well within the accepted limits of the CIS. Ancestral haplotypes were also centred in northern BC.
Main conclusions  High genetic diversity of Oxyria digyna is expected in unglaciated regions, but unexpected in northern BC if British Columbia was entirely covered by ice during the late Pleistocene. Our observations suggest the presence of unglaciated areas providing late Pleistocene refugia in northern BC. Such refugia would have important implications for the origins and migrations of many plant and animal species in north-western North America.     

Evolutionary history of the Maltese wall lizard Podarcis filfolensis: insights on the ‘Expansion–Contraction’ model of Pleistocene biogeography

The expansion–contraction (EC) model predicts demographic and range contraction of temperate species during Pleistocene glaciations as a consequence of climate‐related habitat changes, and provides a paradigm for explaining the high intraspecific diversity found in refugia in terms of long‐term demographic stability. However, recent evidence has revealed a weak predictive power of this model for terrestrial species in insular and coastal settings. We investigated the Pleistocene EC dynamics and their evolutionary consequences on temperate species using the Maltese archipelago and its endemic lizard Podarcis filfolensis as a model system. The evolutionary and demographic history of P. filfolensis as inferred from mitochondrial and nuclear sequences data does not conform to the EC model predictions, supporting (i) demographic and spatial stability or expansion, rather than contraction, of the northern and southern lineages during the last glacial period; and (ii) a major role for allopatric differentiation primed by sea‐level dynamics, rather than prolonged demographic stability, in the formation of the observed genetic diversity. When combined with evidence from other Mediterranean refugia, this study shows how the incorporation of Pleistocene sea‐level variations in the EC model accounts for a reverse demographic and range response of insular and coastal temperate biotas relative to continental ones. Furthermore, this cross‐archipelago pattern in which allopatric diversity is formed and shaped by EC cycles resembles that seen between isolated populations within mainland refugia and suggests that the EC model, originally developed to explain population fluctuations into and out‐of refugia, may be appropriate for describing the demographic and evolutionary dynamics driving the high genetic diversity observed in these areas.     

Quaternary range dynamics of ecologically divergent species (Edraianthus serpyllifolius and E. tenuifolius,Campanulaceae) within the Balkan refugium

Aim We investigated Quaternary range dynamics of two closely related but ecologically divergent species (cold‐tolerant Edraianthus serpyllifolius and thermophilic Edraianthus tenuifolius) with overlapping distribution ranges endemic to the western Balkan Peninsula, an important yet understudied Pleistocene refugium. Our aims were: to test predictions of the ‘refugia‐within‐refugia’ model of strong genetic subdivisions due to population isolation in separate refugia; to explore whether two ecologically divergent species reacted differently to Pleistocene climatic fluctuations; and to test predictions of the displacement refugia model of stronger differentiation among populations in the thermophilic E. tenuifolius compared with the cold‐tolerant E. serpyllifolius. Location The western Balkan Peninsula. Methods We gathered amplified fragment‐length polymorphism (AFLP) data and plastid DNA sequences from two to five individuals from 10 populations of E. serpyllifolius and 22 populations of E. tenuifolius, spanning their entire respective distribution areas. AFLP data were analysed using a Bayesian clustering approach and a distance‐based network approach. Plastid sequences were used to depict relationships among haplotypes in a statistical parsimony network, and to obtain age estimates in a Bayesian framework. Results In E. serpyllifolius, both AFLP and plastid sequence data showed clear geographic structure. Western populations showed high AFLP diversity and a high number of rare fragments. In E. tenuifolius, both markers congruently identified a major phylogeographic split along the lower Neretva valley in central Dalmatia. The most distinct and earliest diverging chloroplast DNA (cpDNA) haplotypes were found further south in the south‐easternmost populations. North‐western populations, identified as a separate cluster by Bayesian clustering, were characterized by low genetic diversity and a low number of rare AFLP markers. Main conclusions Clear evidence for multiple Pleistocene refugia is found not only in the high‐elevation E. serpyllifolius, but also in the lowland E. tenuifolius, despite the lack of obvious dispersal barriers, in line with the refugia‐within‐refugia model. Genealogical relationships and genetic diversity patterns support the hypothesis that cold‐adapted E. serpyllifolius responded to climatic oscillations mostly by elevational range shifts, whereas thermophilic E. tenuifolius did so mainly by latitudinal range shifts, with different phases (and probably extents) of range expansion. In contrast to the displacement refugia hypothesis, the two elevationally differentiated species do not differ in their genetic diversity.     

Phylogeography of <Emphasis Type="Italic">Leiopelma hochstetteri</Emphasis> reveals strong genetic structure and suggests new conservation priorities

Leiopelma hochstetteri belongs to a singular anuran lineage endemic to New Zealand that diverged from other frogs about 200 million years ago. The species now is reduced to a series of isolated populations in the northern half of the North Island and Great Barrier Island. We have used mitochondrial and nuclear sequence data to examine the genetic affinities of extant populations of L. hochstetteri. Phylogenetic reconstructions reveal that populations are highly structured. Each of the geographically isolated populations harbours independent mtDNA lineages as well as different degrees of nuDNA differentiation. Moreover, molecular dating reveals that this structure originated in the early Pleistocene. This pattern of genetic structure likely results from unfavourable climatic conditions during the Pleistocene combined with the low dispersal ability of the species. Isolated populations in forested refugia existed even in the southern part of the distribution of the species during glacial cycles. Previously published variation in chromosome numbers and isozyme data are consistent with the new evidence. We identified 13 evolutionary significant units (ESUs) that should serve as the focus for future management and conservation of this species.     

Mitochondrial gene trees support persistence of cold tolerant fairy shrimp throughout the Pleistocene glaciations in both southern and more northerly refugia

Fairy shrimp (Crustacea: Anostraca) are specialist inhabitants of temporary aquatic habitats. In many parts of the world and particularly in Western Europe, however, populations are declining while the development of adequate conservation strategies is impeded by a poor knowledge of the genetic structure and taxonomic status of remaining lineages. We reconstructed a phylogeography of the species Chirocephalus diaphanus Prévost, 1803 using partial sequences of the mitochondrial COI gene and discuss the importance of different Pleistocene refugia to explain current diversity patterns. In addition to 20 C. diaphanus populations, we also included populations of six presumably closely related chirocephalids to evaluate their taxonomic status. Based on molecular data, the Eastern European subspecies C. diaphanus romanicus deserves species status while the species status of two Italian chirocephalids, C. salinus and C. ruffoi is questionable. Results indicate European C. diaphanus lineages diverged well before the last glacial maximum and survived the Pleistocene glaciations in multiple (sub)refugia along the Iberian, Italian, and Balkan peninsula. Northern Europe was subsequently recolonized from Southern France, resulting in high levels of cryptic diversity around glacial refugia but also in more widespread haplotypes in mainland Europe.     

Imprints of glacial history and current environment on correlations between endemic plant and invertebrate species richness

Aim We evaluate how closely diversity patterns of endemic species of vascular plants, beetles, butterflies, molluscs and spiders are correlated with each other, and to what extent similar environmental requirements or survival in common glacial refugia and comparable dispersal limitations account for their existing congruence. Location Austria. Methods We calculated pairwise correlations among species numbers of the five taxonomic groups in 1405 cells of a 3′ × 5′ raster (c. 35 km²) using the raw data as well as the residuals of regression models that accounted for: (1) environmental variables, (2) environmental variables and the occurrence of potential refugia during the Last Glacial Maximum, or (3) environmental variables, refugia and spatial filters. Results Pairwise cross‐taxonomic group Spearman’s rank correlations in the raw data were significantly positive in most cases, but only moderate (0.3 < ρ < 0.5) to weak (ρ < 0.3) throughout. Correlations were closest between plants and beetles, plants and butterflies, and plants and snails, respectively, whereas the distribution of endemic spiders was largely uncorrelated with those of the other groups. Environmental variables explained only a moderate proportion of the variance in endemic richness patterns, and the response of individual groups to environmental gradients was only partly consistent. The inclusion of refugium locations and the spatial filters increased the goodness of model fit for all five taxonomic groups. Moreover, removing the effects of environmental conditions reduced congruence in endemic richness patterns to a lesser extent than did filtering the influence of refugium locations and spatial autocorrelation, except for spiders, which are probably the least dispersal‐limited of the five groups. Main conclusions The moderate to weak congruence of endemic richness patterns clearly limits the usefulness of a surrogacy approach for designating areas for the protection of regional endemics. On the other hand, our results suggest that dispersal limitations still shape the distributions of many endemic plant, snail, beetle and butterfly species, even at the regional scale; that is, survival in shared refugia and subsequent restricted spread retain a detectable signal in existing correlations. Concentrating conservation efforts on well‐known Pleistocene refugia hence appears to be a reasonable first step towards a strategy for protecting regional endemics of at least the less mobile invertebrate groups.     

High latitudes and high genetic diversity: phylogeography of a widespread boreal bird, the gray jay (Perisoreus canadensis)

We describe range-wide phylogeographic variation in gray jays (Perisoreus canadensis), a boreal Nearctic corvid that occurs today primarily in recently glaciated regions. Phylogenetic analysis of mitochondrial DNA (1041 base pairs ND2 gene; N=205, 50 localities) revealed four reciprocally monophyletic groups. One widespread clade occurs across the North American boreal zone, from Newfoundland to Alaska and southwest into Utah. Three other clades occur at lower latitudes in the montane West in Colorado, the northern Rocky Mountains, and the Pacific Northwest respectively. The geographic distribution of clades in gray jays corresponds with a general pattern that is emerging for boreal taxa, having one widespread northern clade and one or more geographically restricted southwestern clades. Population genetic analyses indicate that the larger boreal clade is genetically structured and harbors significantly more genetic diversity than those clades occurring at lower latitudes. Species distribution modeling (SDM) revealed multiple putative Pleistocene refugia including several occurring at higher latitudes. We suggest that multiple post-glacial colonization routes, some of which originate from these northern refugia, are responsible for the relatively high genetic diversity at high latitudes. Conversely, lower latitude clades show little variation, probably as a result of historical restriction to smaller geographical areas with smaller long-term population sizes. This 'upside-down' pattern of genetic diversity contrasts with the conventional view that populations of north-temperate species occupying previously glaciated habitats should possess lower levels of diversity than their southern counterparts.     

Multilocus tests of Pleistocene refugia and ancient divergence in a pair of Atlantic Forest antbirds (Myrmeciza)

The Atlantic Forest (AF) harbours one of the most diverse vertebrate faunas of the world, including 199 endemic species of birds. Understanding the evolutionary processes behind such diversity has become the focus of many recent, primarily single locus, phylogeographic studies. These studies suggest that isolation in forest refugia may have been a major mechanism promoting diversification, although there is also support for a role of riverine and geotectonic barriers, two sets of hypotheses that can best be tested with multilocus data. Here we combined multilocus data (one mtDNA marker and eight anonymous nuclear loci) from two species of parapatric antbirds, Myrmeciza loricata and M. squamosa, and Approximate Bayesian Computation to determine whether isolation in refugia explains current patterns of genetic variation and their status as independent evolutionary units. Patterns of population structure, differences in intraspecific levels of divergence and coalescent estimates of historical demography fit the predictions of a recently proposed model of refuge isolation in which climatic stability in the northern AF sustains higher diversity and demographic stability than in the southern AF. However, a pre‐Pleistocene divergence associated with their abutting range limits in a region of past tectonic activity also suggests a role for rivers or geotectonic barriers. Little or no gene flow between these species suggests the development of reproductive barriers or competitive exclusion. Our results suggests that limited marker sampling in recent AF studies may compromise estimates of divergence times and historical demography, and we discuss the effects of such sampling on this and other studies.     

The impact of Pleistocene climate change on an ancient arctic-alpine plant: multiple lineages of disparate history in Oxyria digyna

The ranges of arctic-alpine species have shifted extensively with Pleistocene climate changes and glaciations. Using sequence data from the trnH-psbA and trnT-trnL chloroplast DNA spacer regions, we investigated the phylogeography of the widespread, ancient (>3 million years) arctic-alpine plant Oxyria digyna (Polygonaceae). We identified 45 haplotypes and six highly divergent major lineages; estimated ages of these lineages (time to most recent common ancestor, T(MRCA)) ranged from ～0.5 to 2.5 million years. One lineage is widespread in the arctic, a second is restricted to the southern Rocky Mountains of the western United States, and a third was found only in the Himalayan and Altai regions of Asia. Three other lineages are widespread in western North America, where they overlap extensively. The high genetic diversity and the presence of divergent major cpDNA lineages within Oxyria digyna reflect its age and suggest that it was widespread during much of its history. The distributions of individual lineages indicate repeated spread of Oxyria digyna through North America over multiple glacial cycles. During the Last Glacial Maximum it persisted in multiple refugia in western North America, including Beringia, south of the continental ice, and within the northern limits of the Cordilleran ice sheet. Our data contribute to a growing body of evidence that arctic-alpine species have migrated from different source regions over multiple glacial cycles and that cryptic refugia contributed to persistence through the Last Glacial Maximum.     

Patterns of freshwater biodiversity in Europe: lessons from the spring snail genus Bythinella

Aim We analyse patterns of biodiversity in the spring snail genus Bythinella, a group of highly isolated and stenotopic freshwater species. We aim to test: (1) whether there are European areas of increased diversity (i.e. ‘hotspots’), (2) whether the potential hotspots inferred show qualitative differences in biodiversity characteristics such as endemicity, distinctiveness of taxa, age of lineages or degree of fragmentation, and (3) whether these hotspots match the Pleistocene refugia of Bythinella spp. Location Europe, Asia Minor. Methods The analyses are based on genetic data from 717 Bythinella specimens sampled at 194 sites. We used haplotypes as operational units in all analyses. To test hypothesis 1, mean pairwise genetic distances between Bythinella populations within each 1° × 1° geographical grid cell sampled in Europe were calculated. Within individual mountain ranges, grid cells with high diversity were grouped with neighbouring ones and hotspots were identified based on pre‐defined criteria. Then, to test hypothesis 2, different biodiversity indices of these regions were calculated and compared. Finally, to test hypothesis 3, the spatial distribution of the identified hotspots was compared with the known Pleistocene refugia of Bythinella spp. Results Five areas showed increased levels of genetic diversity: the Massif Central/Pyrenees, the western and eastern Alps, and the western and eastern Carpathians. These regions showed qualitative differences in biodiversity, with the eastern Carpathians holding the highest number of (endemic) haplotypes, the oldest and most distinct lineages and the highest degree of fragmentation. Only three of the five detected hotspots matched previously identified Pleistocene refugia for Bythinella spp. Main conclusions The genetic diversity of Bythinella spp. is not randomly distributed throughout Europe. Some of the hotspots we identify coincide with those found in other freshwater taxa; others have not previously been reported. Thus, spring organisms may reflect a unique evolutionary history that is distinct from lentic and lotic taxa. Our findings may be useful for conservation purposes even though the species‐level taxonomy of the genus is still under discussion.