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.     

Population genetic structure of the bank vole Myodes glareolus within its glacial refugium in peninsular Italy

It is now well established that Southern European peninsulas have been major glacial refugia for temperate species during Pleistocene climatic oscillations. However, substantial environmental changes occurred also within these peninsulas throughout the Pleistocene, raising questions about the role and interplay of various microevolutionary processes in shaping patterns of intraspecific diversity within these areas. Here, we investigate the patterns of genetic variation in the bank vole Myodes glareolus within the Italian peninsula. By using a panel of 13 microsatellite loci, we found more intraspecific variation than expected based on previous assessments. Indeed, both Bayesian and ordination‐based clustering analyses of variation recovered five main geographic/genetic clusters along the peninsula, with three clusters geographically restricted to the southern portion of the peninsula. This clustering is supported by previous evidences of some morphological distinctiveness among these populations. This pattern can be explained by a refugia‐within‐refugia scenario, with the occurrence of multiple sub‐refugia for the bank vole within the Italian peninsula, likely promoted by the major palaeo‐environmental changes which affected forested habitats within this area during the Pleistocene. Moreover, our results support a scenario whereby the high levels of intraspecific diversity observed within major Pleistocene refugia are better explained by dynamic microevolutionary processes occurred within these areas, rather than by long‐term demographic stability of refugial population. Finally, the narrow and isolated distribution of some of the identified lineages suggests the need for future assessments of their conservation and taxonomic status.     

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.     

Phylogeography of the widespread African puff adder (Bitis arietans) reveals multiple Pleistocene refugia in southern Africa

Evidence from numerous Pan‐African savannah mammals indicates that open‐habitat refugia existed in Africa during the Pleistocene, isolated by expanding tropical forests during warm and humid interglacial periods. However, comparative data from other taxonomic groups are currently lacking. We present a phylogeographic investigation of the African puff adder (Bitis arietans), a snake that occurs in open‐habitat formations throughout sub‐Saharan Africa. Multiple parapatric mitochondrial clades occur across the current distribution of B. arietans, including a widespread southern African clade that is subdivided into four separate clades. We investigated the historical processes responsible for generating these phylogeographic patterns in southern Africa using species distribution modelling and genetic approaches. Our results show that interior regions of South Africa became largely inhospitable for B. arietans during glacial maxima, whereas coastal and more northerly areas remained habitable. This corresponds well with the locations of refugia inferred from mitochondrial data using a continuous phylogeographic diffusion model. Analysis of data from five anonymous nuclear loci revealed broadly similar patterns to mtDNA. Secondary admixture was detected between previously isolated refugial populations. In some cases, this is limited to individuals occurring near mitochondrial clade contact zones, but in other cases, more extensive admixture is evident. Overall, our study reveals a complex history of refugial isolation and secondary expansion for puff adders and a mosaic of isolated refugia in southern Africa. We also identify key differences between the processes that drove isolation in B. arietans and those hypothesized for sympatric savannah mammals.     

Population genetic structure and diversity of the Apennine endemic stream frog, Rana italica--insights on the Pleistocene evolutionary history of the Italian peninsular biota

For most species in the Western Palaearctic region, southern Mediterranean peninsulas have been identified as major Quaternary refugia and hotspots of intraspecific diversity, and thus, as areas of particular relevance for the conservation of the evolutionary potential. We analysed the patterns of geographical variation among 26 populations of the Italian stream frog, using both nuclear (allozymes) and mitochondrial (partial cytochrome b sequences) markers. Phylogenetic, phylogeographical and population genetic analyses suggested that the species survived the last glacial–interglacial cycles in two distinct refugia, one restricted to the tip of the Calabrian peninsula, at the extreme south of the species' range, the other spanning from central Calabria to central Apennines and showing evidences for further population subdivision therein. Historical demographic tests suggested a significant population expansion from the latter, which most likely began around the last pleniglacial. This expansion would have led to the rapid colonization of the northern Apennines to the north, and to a secondary contact and population admixture with the population from the southern refugium in southern central Calabria. A comparison of the evolutionary history inferred for the Italian stream frog with the data emerging for other codistributed species suggests: (i) the generality of a multiple-refugia scenario for the Italian peninsula, (ii) the possible occurrence of at least one suture zone in southern Italy, and (iii) that for most species, this Pleistocene refugium is not only a hotspot, but also a melting pot of intraspecific genetic diversity. Finally, the conservation implications of these results are also briefly highlighted.     

Pleistocene evolutionary history of the Clouded Apollo (Parnassius mnemosyne): genetic signatures of climate cycles and a 'time-dependent' mitochondrial substitution rate

Genetic data are currently providing a large amount of new information on past distribution of species and are contributing to a new vision of Pleistocene ice ages. Nonetheless, an increasing number of studies on the 'time dependency' of mutation rates suggest that date assessments for evolutionary events of the Pleistocene might be overestimated. We analysed mitochondrial (mt) DNA (COI) sequence variation in 225 Parnassius mnemosyne individuals sampled across central and eastern Europe in order to assess (i) the existence of genetic signatures of Pleistocene climate shifts; and (ii) the timescale of demographic and evolutionary events. Our analyses reveal a phylogeographical pattern markedly influenced by the Pleistocene/Holocene climate shifts. Eastern Alpine and Balkan populations display comparatively high mtDNA diversity, suggesting multiple glacial refugia. On the other hand, three widely distributed and spatially segregated lineages occupy most of northern and eastern Europe, indicating postglacial recolonization from different refugial areas. We show that a conventional 'phylogenetic' substitution rate cannot account for the present distribution of genetic variation in this species, and we combine phylogeographical pattern and palaeoecological information in order to determine a suitable intraspecific rate through a Bayesian coalescent approach. We argue that our calibrated 'time-dependent' rate (0.096 substitutions/ million years), offers the most convincing time frame for the evolutionary events inferred from sequence data. When scaled by the new rate, estimates of divergence between Balkan and Alpine lineages point to c. 19 000 years before present (last glacial maximum), and parameters of demographic expansion for northern lineages are consistent with postglacial warming (5-11 000 years before present).     

Mitochondrial DNA phylogeography of Lissotriton boscai (Caudata, Salamandridae): evidence for old, multiple refugia in an Iberian endemic

In Europe, southern peninsulas served as refugia during cold periods in the Pleistocene, acting both as centres of origin of endemisms and as sources from which formerly glaciated areas were recolonized during interglacial periods. Previous studies have revealed that within the main refugial areas, intraspecific lineages often survived in allopatric refugia. We analysed two mitochondrial markers (nad4, control region, approximately 1.4 kb) in 103 individuals representing the entire distribution of Lissotriton boscai, a newt endemic to the western Iberian Peninsula. We inferred the evolutionary history of the species through phylogenetic, phylogeographic and historical demographic analyses. The results revealed unexpected, deep levels of geographically structured genetic variability. We identified two main evolutionary lineages, each containing three well-supported clades. The first historical split involved populations from central-southwestern coastal Portugal and the ancestor of all the remaining populations around 5.8 million years ago. Both lineages were subsequently fragmented into different population groups between 2.5 and 1.2 million years ago. According to nested clade analysis, at lower hierarchical levels the patterns suggest restricted gene flow with isolation by distance, whereas at higher levels the clades exhibit signatures of contiguous range expansion. Bayesian Skyline Plots show recent bottlenecks, followed by demographic expansions in all lineages. The significant genetic structure found is consistent with long-term survival of populations in allopatric refugia, supporting the 'refugia-within-refugia' scenario for southern European peninsulas. The comparison of our results with other co-distributed species highlights the generality of this hypothesis for the Iberian herpetofauna and suggests that Mediterranean refuges had more relevance for the composition and distribution of present biodiversity patterns than currently acknowledged. We briefly discuss the taxonomic and conservation implications of our results.     

On glacial refugia,genetic diversity,and microevolutionary processes: deep phylogeographical structure in the endemic newt Lissotriton italicus1

We investigated the phylogeographical patterns of Lissotriton italcus, a newt endemic to the Italian peninsula, aiming to determine why hotspots of intraspecific diversity so ‘hot’. We found two main mitochindrial DNA lineages (net sequence divergence of 6.8% at two fragments of total length of 1897 bp): one restricted to part of the Calabrian peninsula (i.e the southernmost portion of the species range) and the other widespread throughout the rest of the species range. Both lineages, which had a parapatric distribution, showed evidence of further subdivisions, with an overall number of eight terminal haplogroups, most of whose times to the most recent common ancestors were estimated at the Late Pleistocene. Analysis of molecular variance suggested that partitioning populations according to the geographical distribution of these haplogroups can explain 97% of the observed genetic variation. These results suggest that L. italicus underwent repeated cycles of allopatric fragmentation throughout the Pleistocene, and that it likely survived the Late Pleistocene paleoenvironmental changes within eight separate refugia. Thus, the current hotspot of intraspecific diversity of L. italicus (within the Calabrian peninsula) has not been moulded by long‐term stability of large populations but rather by multiple events of allopatric fragmentation and divergence. When compared with the patterns recently identified in other species, these results suggest that the occurrence of phases of allopatric divergence (eventually followed by secondary admixture) could be a common, albeit probably underrated feature in the history of formation of hotspots of intraspecific diversity. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 105 , 42–55.     

Out of Anatolia: longitudinal gradients in genetic diversity support an eastern origin for a circum-Mediterranean oak gallwasp Andricus quercustozae

Many studies have addressed the latitudinal gradients in intraspecific genetic diversity of European taxa generated during postglacial range expansion from southern refugia. Although Asia Minor is known to be a centre of diversity for many taxa, relatively few studies have considered its potential role as a Pleistocene refugium or a potential source for more ancient westward range expansion into Europe. Here we address these issues for an oak gallwasp, Andricus quercustozae (Hymenoptera: Cynipidae), whose distribution extends from Morocco along the northern coast of the Mediterranean through Turkey to Iran. We use sequence data for a fragment of the mitochondrial gene cytochrome b and allele frequency data for 12 polymorphic allozyme loci to answer the following questions: (1) which regions represent current centres of genetic diversity for A. quercustozae? Do eastern populations represent one refuge or several discrete glacial refugia? (2) Can we infer the timescale and sequence of the colonization processes linking current centres of diversity? Our results suggest that A. quercustozae was present in five distinct refugia (Iberia, Italy, the Balkans, southwestern Turkey and northeastern Turkey) with recent genetic exchange between Italy and Hungary. Genetic diversity is greatest in the Turkish refugia, suggesting that European populations are either (a) derived from Asia Minor, or (b) subject to more frequent population bottlenecks. Although Iberian populations show the lowest diversity for putatively selectively neutral markers, they have colonized a new oak host and represent a genetically and biologically discrete entity within the species.     

Genetic diversity and phylogeography of the Apennine yellow-bellied toad Bombina pachypus, with implications for conservation

Genetic variation was investigated in 17 populations of the Italian endemic Apennine yellow-bellied toad using both mitochondrial (598 bp of the cytochrome b gene) and nuclear (21 allozyme loci) markers. Populations from central Calabria (southern Italy) showed the highest levels of intrapopulation genetic variation, whereas samples located north of this region were nearly lacking in variation. This appears to be a typical pattern of 'southern richness and northern purity', usually attributed to the prolonged population stability within southern refugia coupled with the loss of variation during postglacial northward expansion. However, the overall pattern of genetic variation observed has a strong geographical component, suggesting two Calabrian plains, Catanzaro and Crati-Sibari, as historical barriers to dispersal separating three population groups. These findings cannot be explained by the prolonged stability of southern populations alone, and suggest that the southern richness has been at least in part shaped by allopatric differentiation within the refugial range, followed by intermixing of previously differentiated lineages. From a conservation standpoint, Calabria is the major genetic diversity reservoir for this species, thus deserving particular conservation efforts. Furthermore, although the low intrapopulation genetic variation outside Calabria appears to be of clear historical origin, evidence of a current reduction of gene flow suggests that human disturbance has also played a part, particularly in the anthropogenic impacted Volturno river drainage basin.     

Phylogeography of Potamon ibericum (Brachyura: Potamidae) identifies Quaternary glacial refugia within the Caucasus biodiversity hot spot

Refugia are critical for the maintenance of biodiversity during the periods of Quaternary climatic oscillations. The long‐term persistence of refugial populations in a large continuous refugium has resulted in a homogenous pattern of genetic structure among populations, while highly structured evolutionary lineages characterize the restriction of refugial populations to smaller subrefugia. These mechanisms have resulted in the identification of hot spots of biodiversity within putative glacial refugia. We studied phylogeography of Potamon ibericum (Brachyura: Potamidae) in the drainages of the western Caucasus biodiversity hot spot (i.e., Colchis and the Caucasus) to infer spatial genetic structure and potential refugia for a freshwater crab in this region. These areas have traditionally considered as a refugium due to the presence of Tertiary relict species. We integrated population genetic data and historical demographic analysis from cytochrome oxidase subunit I sequences and paleoclimatic data from species distribution modeling (SDM). The results revealed the lack of phylogeographic structure and provided evidence for demographic expansion. The SDM presented a rather homogenous and large refugium that extended from northeast Turkey to Colchis during the last glacial period. In contrast to these findings, previous phylogeographic study on P. ibericum of the eastern Caucasus biodiversity hot spot (i.e., Hyrcania) identified multiple independent refugia. By combining these results, we explain the significance of this important western Palearctic hot spot of biological diversity in shaping the geographic distribution of intraspecific genetic diversity in a freshwater taxon.     

Stronger spatial genetic structure in recolonized areas than in refugia in the European beech

Extant rear‐edge populations located in former glacial refugia remain understudied despite their high conservation value. These populations should have experienced strong genetic drift due to their small size and long isolation. Moreover, the prolonged action of isolation by distance in refugial areas should result in stronger regional spatial genetic structure (SGS) than in recolonized areas, but empirical tests of this prediction are scarce. To fill this gap, we first used a set of 16 microsatellite markers to investigate the genetic structure of European beech in France in 65 populations from three refugial areas and one control recolonized (nonrefugial) area. Then, using the same approach, we reanalysed published isozyme data from 375 populations distributed across the entire species range. We found stronger genetic differentiation among populations in refugia than in recolonized areas. However, contrary to expectations, regional SGS was lower within refugia than within recolonized areas. Published studies presenting similar analyses suggest that our results could have generality across different biogeographical settings and types of organisms. Strong and prolonged genetic drift in refugial areas could have erased the signature of range expansions that is still visible in recolonized areas. Our results therefore suggest that Pleistocene population isolation has played a key role in increasing the genetic complexity of extant rear‐edge populations.     

Phylogeography and historical demography of the Italian treefrog, Hyla intermedia, reveals multiple refugia, population expansions and secondary contacts within peninsular Italy

We investigated the geographical patterns of genetic diversity in the Italian treefrog through sequence analysis of a mitochondrial cytochrome b gene fragment. Three main mitochondrial lineages were identified, distributed in northern, central and southern Italy, respectively. Their divergence appears indicative of a split time largely predating Late Pleistocene climatic oscillations, and syntopy between them was only observed in the geographically intermediate populations. The historical demographic reconstructions suggest that in both northern and central Italy, an expansion occurred during the last major glacial phase, when a vast widening of the lowland habitats followed the glaciation-induced fall of the sea level. Instead, in southern Italy an expansion event likely followed the end of the last glaciation, although the inference of expansion appears less reliable for the southern clade than for the others. Within this geographical area, a sharp phylogeographic discontinuity separated peninsular from Sicilian populations, and the overall pattern of diversity suggests that the latter derived from a recent colonization of the island, probably through a Late Pleistocene land bridge. Phylogenetic, phylogeographic and historical demographic analyses thus concur in delineating a scenario of multiple refugia, with four groups of populations which survived the last glacial-interglacial cycles in at least three distinct refugia arranged along peninsular Italy, and have recently come into contact following range expansions. Therefore, these results support the hypothesis that a plethora of microevolutionary processes, rather than the prolonged stability of populations, were mainly responsible for shaping the patterns of diversity within this major biodiversity hotspot.     

Evolutionary melting pots: a biodiversity hotspot shaped by ring diversifications around the Black Sea in the Eastern tree frog (Hyla orientalis)

Hotspots of intraspecific genetic diversity, which are of primary importance for the conservation of species, have been associated with glacial refugia, that is areas where species survived the Quaternary climatic oscillations. However, the proximate mechanisms generating these hotspots remain an open issue. Hotspots may reflect the long‐term persistence of large refugial populations; alternatively, they may result from allopatric differentiation between small and isolated populations, that later admixed. Here, we test these two scenarios in a widely distributed species of tree frog, Hyla orientalis, which inhabits Asia Minor and southeastern Europe. We apply a fine‐scale phylogeographic survey, combining fast‐evolving mitochondrial and nuclear markers, with a dense sampling throughout the range, as well as ecological niche modelling, to understand what shaped the genetic variation of this species. We documented an important diversity centre around the Black Sea, composed of multiple allopatric and/or parapatric diversifications, likely driven by a combination of Pleistocene climatic fluctuations and complex regional topography. Remarkably, this diversification forms a ring around the Black Sea, from the Caucasus through Anatolia and eastern Europe, with terminal forms coming into contact and partially admixing in Crimea. Our results support the view that glacial refugia generate rather than host genetic diversity and can also function as evolutionary melting pots of biodiversity. Moreover, we report a new case of ring diversification, triggered by a large, yet cohesive dispersal barrier, a very rare situation in nature. Finally, we emphasize the Black Sea region as an important centre of intraspecific diversity in the Palearctic with implications for conservation.     

Molecular phylogeography of Dryas integrifolia: glacial refugia and postglacial recolonization

Chloroplast DNA variation in the Arctic plant species Dryas integrifolia (Rosaceae) was analysed in relation to both the present-day geographical distribution of populations and to Pleistocene fossil records of this species. The phylogeographical structure was weak but the analysis of haplotype diversity revealed several groups of haplotypes having present-day geographical ranges that overlap locations postulated from geographical and fossil evidence to have been glacial refugia. Based on this information we infer that two important refugial sources of Arctic recolonization by this species were Beringia and the High Arctic. Two other putative refugia, located southeast of the ice sheet and along coastal regions of the eastern Arctic may have served as sources for recolonization of smaller portions of the Arctic. The genetic substructure in the species is mostly due to variation among populations regardless of the ecogeographical region in which they are found. Spatial autocorrelation at the regional scale was also detected. High levels of diversity both within populations and ecogeographical regions are probably indicative of population establishment from several sources possibly combined with recent gene flow.     

Endemic species may have complex histories: within‐refugium phylogeography of an endangered Iberian vole

Glacial refugia protected and promoted biodiversity during the Pleistocene, not only at a broader scale, but also for many endemics that contracted and expanded their ranges within refugial areas. Understanding the evolutionary history of refugial endemics is especially important in the case of endangered species to recognize the origins of their genetic structure and thus produce better informed conservation practices. The Iberian Peninsula is an important European glacial refugium, rich in endemics of conservation concern, including small mammals, such as the Cabrera vole (Microtus cabrerae). This near‐threatened rodent is characterized by an unusual suite of genetic, life history and ecological traits, being restricted to isolated geographic nuclei in fast‐disappearing Mediterranean subhumid herbaceous habitats. To reconstruct the evolutionary history of the Cabrera vole, we studied sequence variation at mitochondrial, autosomal and sex‐linked loci, using invasive and noninvasive samples. Despite low overall mitochondrial and nuclear nucleotide diversities, we observed two main well‐supported mitochondrial lineages, west and east. Phylogeographic modelling in the context of the Cabrera vole's detailed fossil record supports a demographic scenario of isolation of two populations during the Last Glacial Maximum from a single focus in the southern part of the Iberian Peninsula. In addition, our data suggest subsequent divergence within the east, and secondary contact and introgression of the expanding western population, during the late Holocene. This work emphasizes that refugial endemics may have a phylogeographic history as rich as that of more widespread species, and conservation of such endemics includes the preservation of that genetic legacy.     

Comparative phylogeography reveals a shared impact of pleistocene environmental change in shaping genetic diversity within nine Anopheles mosquito species across the Indo-Burma biodiversity hotspot

South-East Asia is one of the world's richest regions in terms of biodiversity. An understanding of the distribution of diversity and the factors shaping it is lacking, yet essential for identifying conservation priorities for the region's highly threatened biodiversity. Here, we take a large-scale comparative approach, combining data from nine forest-associated Anopheles mosquito species and using statistical phylogeographical methods to disentangle the effects of environmental history, species-specific ecology and random coalescent effects. Spatially explicit modelling of Pleistocene demographic history supports a common influence of environmental events in shaping the genetic diversity of all species examined, despite differences in species' mtDNA gene trees. Populations were periodically restricted to allopatric northeastern and northwestern refugia, most likely due to Pleistocene forest fragmentation. Subsequent southwards post-glacial recolonization is supported by a north-south gradient of decreasing genetic diversity. Repeated allopatric fragmentation and recolonization have led to the formation of deeply divergent geographical lineages within four species and a suture zone where these intraspecific lineages meet along the Thai-Myanmar border. A common environmental influence for this divergence was further indicated by strong support for simultaneous divergence within the same four species, dating to approximately 900 thousand years ago (kya). Differences in the geographical structuring of genetic diversity between species are probably the result of varying species' biology. The findings have important implications for conservation planning; if the refugial regions and suture zone identified here are shared by other forest taxa, the unique and high levels of genetic diversity they house will make these areas conservation priorities.     

Population genomic evidence for multiple Pliocene refugia in a montane‐restricted harvestman (Arachnida,Opiliones, Sclerobunus robustus) from the southwestern United States

The integration of ecological niche modelling into phylogeographic analyses has allowed for the identification and testing of potential refugia under a hypothesis‐based framework, where the expected patterns of higher genetic diversity in refugial populations and evidence of range expansion of nonrefugial populations are corroborated with empirical data. In this study, we focus on a montane‐restricted cryophilic harvestman, Sclerobunus robustus, distributed throughout the heterogeneous Southern Rocky Mountains and Intermontane Plateau of southwestern North America. We identified hypothetical refugia using ecological niche models (ENMs) across three time periods, corroborated these refugia with population genetic methods using double‐digest RAD‐seq data and conducted population‐level phylogenetic and divergence dating analyses. ENMs identify two large temporally persistent regions in the mid‐latitude highlands. Genetic patterns support these two hypothesized refugia with higher genetic diversity within refugial populations and evidence for range expansion in populations found outside hypothesized refugia. Phylogenetic analyses identify five to six genetically divergent, geographically cohesive clades of S. robustus. Divergence dating analyses suggest that these separate refugia date to the Pliocene and that divergence between clades pre‐dates the late Pleistocene glacial cycles, while diversification within clades was likely driven by these cycles. Population genetic analyses reveal effects of both isolation by distance (IBD) and isolation by environment (IBE), with IBD more important in the continuous mountainous portion of the distribution, while IBE was stronger in the populations inhabiting the isolated sky islands of the south. Using model‐based coalescent approaches, we find support for postdivergence migration between clades from separate refugia.     

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.     

Are glacial refugia hotspots of speciation and cytonuclear discordances? Answers from the genomic phylogeography of Spanish common frogs

Subdivided Pleistocene glacial refugia, best known as “refugia within refugia”, provided opportunities for diverging populations to evolve into incipient species and/or to hybridize and merge following range shifts tracking the climatic fluctuations, potentially promoting extensive cytonuclear discordances and “ghost” mtDNA lineages. Here, we tested which of these opposing evolutionary outcomes prevails in northern Iberian areas hosting multiple historical refugia of common frogs (Rana cf. temporaria), based on a genomic phylogeography approach (mtDNA barcoding and RAD‐sequencing). We found evidence for both incipient speciation events and massive cytonuclear discordances. On the one hand, populations from northwestern Spain (Galicia and Asturias, assigned to the regional endemic R. parvipalmata), are deeply‐diverged at mitochondrial and nuclear genomes (~4 My of independent evolution), and barely admix with northeastern populations (assigned to R. temporaria sensu stricto) across a narrow hybrid zone (~25 km) located in the Cantabrian Mountains, suggesting that they represent distinct species. On the other hand, the most divergent mtDNA clade, widespread in Cantabria and the Basque country, shares its nuclear genome with other R. temporaria s. s. lineages. Patterns of population expansions and isolation‐by‐distance among these populations are consistent with past mitochondrial capture and/or drift in generating and maintaining this ghost mitochondrial lineage. This remarkable case study emphasizes the complex evolutionary history that shaped the present genetic diversity of refugial populations, and stresses the need to revisit their phylogeography by genomic approaches, in order to make informed taxonomic inferences.