Persistent genetic signatures of colonization in Brachionus manjavacas rotifers in the Iberian Peninsula

Recent phylogeographical assessments have consistently shown that continental zooplankton display high levels of population subdivision, despite the high dispersal capacity of their diapausing propagules. As such, there is an apparent paradox between observed cosmopolitanism in the zooplankton that is associated with long-distance dispersal, and strong phylogeographical structures at a regional scale. Such population dynamics, far from migration-drift equilibrium, have been shown in the rotifer species complex Brachionus plicatilis, a group of over a dozen species inhabiting salt lakes and coastal lagoons worldwide. Here we present the mitochondrial DNA phylogeography of one of these species, Brachionus manjavacas, in the Iberian Peninsula, where it often co-occurs with the morphologically similar species B. plicatilis sensu stricto. We obtained sequences from 233 individuals from diapausing eggs and clonal cultures from 16 lakes in the Iberian Peninsula, and a Tunisian lake. Two strongly supported deep mitochondrial DNA clades were found (A and B). Phylogenetic and nested clade analysis showed that clade A has a strong phylogeographical structure, with a strong similarity of phylogeographical patterns between B. manjavacas clade A and B. plicatilis s.s. These include (i) signatures of allopatric fragmentation between central and southern populations, and (ii) range expansions in the Iberian Peninsula, both likely to have occurred during the Pleistocene. We find evidence for a glacial refugium in the Guadiana basin. Clades A and B co-occurred in several of these lakes because of range expansion and secondary contact between both clades. The co-occurrence between B. plicatilis s.s. and B. manjavacas is not recent, and both species might have experienced similar environmental challenges during the Pleistocene. The strong correlation of genetic and geographical distance found suggests that historical events can lead to such correlation, mirroring the effects of 'isolation by distance' in equilibrium populations.     

Conservation genetics of the rare Iberian endemic Cheirolophus uliginosus (Asteraceae)

Cheirolophus uliginosus is a rare species, endemic to the south‐western Iberian Peninsula, and listed as a characteristic taxon from the temperate Atlantic wet heaths, a priority habitat for conservation by the European Union. The conservation status of this species in most of its distribution area is poorly known, but, in recent times, some populations have disappeared and there has been a reduction in the number of individuals in others. In this context, we analysed the effects of population size on genetic diversity, revealing that genetic erosion and inbreeding depression could be having a significant impact on smaller populations. Furthermore, we studied the patterns of genetic structure and variability at the species level, finding a strikingly low within‐population diversity and high among‐population genetic differentiation. Finally, the genetic structure analyses suggested a long and complex phylogeographical history of C. uliginosus in the region, in agreement with the climate relict status proposed for this species. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 179 , 157–171.     

Past tree range dynamics in the Iberian Peninsula inferred through phylogeography and palaeodistribution modelling: A review

The forests in the Iberian Peninsula have been strongly influenced by past climatic changes, but reconstructing their historical distributions and dynamics is very difficult due to the complex climatic characteristics and relief of the region. Research disciplines such as phylogeography and species distribution modelling can describe the past range dynamics of individual tree species in relatively great detail and help elucidate how these species have reacted to climatic changes. Here we review phylogeographical and modelling studies from species representative of the major Iberian forest types and attempt to extract general trends from the diversity of individual species histories in the Peninsula. To date most studies focus on geographical range dynamics during the Pleistocene, but an increasing body of evidence shows that some species have also retained the genetic imprints of much more ancient processes. Many widespread Iberian species show a deep and often remarkably clear-cut divide between populations from the Mediterranean and from the Atlantic regions of the Iberian Peninsula, suggesting that both areas have independently sustained viable populations over extended periods. In fact, phylogeographical studies commonly find that species had several glacial refugia across the Iberian Peninsula. On the other hand, distribution models help identifying further suitable areas that could have sustained so far undetected refugia. Such studies are of interest for species conservation, because refugium populations are high-priority targets due to their long-term persistence and unique evolutionary trajectory. Overall, we conclude that palaeoecology, phylogeography and species distribution modelling have a great potential to inform each other because of their complementary perspectives and results. A true integration of these approaches is therefore fundamental for further progress in our understanding of past Iberian environments and the organisms they harboured.     

Refugia within refugia – patterns in endemism and genetic divergence are linked to Late Quaternary climate stability in the Iberian Peninsula

In Europe, southern peninsulas served as major refugia during Pleistocene cold periods. However, growing evidence has revealed complex patterns of glacial survival within these southern regions, with multiple glacial refugia within each larger refugial area. We investigated the extent to which patterns of endemism and phylogeographic are concordant across animal species in the Iberian Peninsula, one of the most important unglaciated areas in Europe during the Pleistocene, can be explained in terms of climatic stability. We found that historical climatic stability (notably climate velocity measures integrating macroclimatic shifts with local spatial topoclimate gradients) was often among the most important predictors of endemic species richness for different taxonomic groups using models that also incorporated measures of modern climate. Furthermore, for some taxonomic groups, climatic stability was also correlated with patterns of spatial concordance in interpopulation genetic divergence across multiple taxa, and private haplotypes were more frequently found in relatively stable areas. Overall, our results suggest that both endemism patterns and cross‐taxa concordant phylogeographic patterns across the Iberian Peninsula to some extent are linked to spatial variation in Late Quaternary climate stability, in agreement with the proposed ‘refugia‐within‐refugia’ scenario. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 13–28.     

Phylogeographic patterns and high levels of chloroplast DNA diversity in four Packera (asteraceae) species in southwestern Alberta

Chloroplast DNA (cpDNA) haplotype variation is compared among alpine and prairie/montane species of Packera from a region in southwestern Alberta that straddles the boundary of Pleistocene glaciation. The phylogeny of the 15 haplotypes identified reveals the presence of two groups: one generally found in coastal and northern species and the other from species in drier habitats. The presence of both groups in all four species and most populations from southwestern Alberta is evidence of past hybridization involving species or lineages that may no longer be present in the region. With the exception of the alpine P. subnuda (phiST = 1.0), interpopulational subdivision of haplotype variation is low (phiST < 0.350), suggesting that interpopulational gene flow is high. However, based on haplotype distribution patterns, we propose that Pleistocene hybridization and incomplete lineage sorting have resulted in reduced subdivision of interpopulational variation so that gene flow may not be as high as indicated. Drift has been more important in the alpine species populations, especially P. subnuda.     

The evolutionary history of the allopolyploid Squalius alburnoides (Cyprinidae) complex in the northern Iberian Peninsula

Understanding the population structure, population dynamics and processes that give rise to polyploidy and helps to maintain it is central to our knowledge of the evolution of asexual vertebrates. Previous studies revealed high genetic diversity and several reproductive pathways in the southern populations of the Squalius alburnoides hybrid complex. In contrast, lower genetic variability and the associated limited chance of introducing new genetic combinations may threaten the survival of the northern Mondego populations. We analysed the genetic diversity and structure of nine populations of S. alburnoides in the Iberian Peninsula using microsatellite loci to provide further insights on the evolutionary history of this complex. Special attention was given to the less-studied northern populations (Mondego and Douro basins). Marked population structure, a high frequency of private alleles and a high diversity of some biotypes in the Douro basin indicate that some northern populations may not be at high risk of extinction, contrary to what was expected. The genetic diversity found in the northern Douro populations contradicts the general trend of remarkable genetic impoverishment northwards that occurs in other species and regions. The results indicate the possible existence of a glacial refugium in the Rabaçal River, corroborating findings in other species of this region. Historical events seem to have affected the geographical patterns of genetic variability found among and within the northern and southern populations of this complex and contributed to different patterns of genome composition. Therefore, historical events might have a major role in the long-term persistence of some polyploid hybrid taxa.     

Phylogeographic structure and past history of the circum-Mediterranean species Tomicus destruens Woll. (Coleoptera: Scolytinae)

Phylogeographic studies are often focused on temperate European species with relict footholds in the Mediterranean region. Past climatic oscillations usually induced range contractions and expansions from refugial areas located in southern Europe, and spatial distribution of genetic diversity show that northward expansions were usually pioneer-like. Actually, few studies have focused on circum-Mediterranean species, which probably were not influenced in the same way by climatic oscillations. We present the phylogeography of the bark beetle Tomicus destruens, which is restricted to the whole Mediterranean basin and the Atlantic coasts of North Africa and Portugal. We systematically sequenced 617 bp of the mitochondrial genes COI and COII for 42 populations (N = 219). Analysis revealed 53 haplotypes geographically structured in two clades, namely eastern and western clades, that diverged during the Pleistocene. A contact zone was identified along the Adriatic coast of Italy. Interestingly, we found contrasting levels of genetic structure within each clade. The eastern group was characterized by a significant phylogeographic pattern and low levels of gene flow, whereas the western group barely showed a spatial structure in haplotype distribution. Moreover, the main pine hosts were different between groups, with the Aleppo-brutia complex in the east and the maritime pine in the west. Potential roles of host species, climatic parameters and geographical barriers are discussed and the phylogeographic patterns are compared to classical models of postglacial recolonization in Europe.     

基于叶绿体微卫星研究鄂报春谱系遗传结构

鄂报春Primula obconica作为一种广泛栽培的园艺植物,其野生居群的遗传多样性及遗传结构的研究还少见报道。本文通过叶绿体微卫星分析了17个鄂报春野生居群(共278个个体),共发现4个多态性位点(16个等位基因),得到14个单倍型。结果表明鄂报春具有很高的总基因多样性(HT=0.971)和极低的居群内基因多样性(HS=0.028);分子方差分析(AMOVA)显示98%的变异存在于居群间。这些结果说明早期的生境片断化及有限的种子传播能力是造成当前遗传结构的重要原因。Nst值显著大于Gst值,表明关系相近的单倍型会出现在相同的地区内,同时最小生成树(MST)的分析结果证实了这样的结论。我们在最小生成树的两个组中推断出一些古老单倍型,并推测在冰期时湖北和我国的西南地区可能是该物种的避难所。     

Steppes, savannahs, forests and phytodiversity reservoirs during the Pleistocene in the Iberian Peninsula

A palaeobotanical analysis of the Pleistocene floras and vegetation in the Iberian Peninsula shows the existence of patched landscapes with Pinus woodlands, deciduous and mixed forests, parklands (savannah-like), shrublands, steppes and grasslands. Extinctions of Arctotertiary woody taxa are recorded during the Early and Middle Pleistocene, but glacial refugia facilitated the survival of a number of temperate, Mediterranean and Ibero-North African woody angiosperms. The responses of Iberian vegetation to climatic changes during the Pleistocene have been spatially and temporarily complex, including rapid changes of vegetation in parallel to orbital and suborbital variability, and situations of multi-centennial resilience or accommodation to climatic changes. Regional characteristics emerged as soon as for the Middle Pleistocene, if not earlier: Ericaceae in the Atlantic coast indicating wetter climate, thermo-mediterranean elements in the south as currently, and broad-leaf trees in the northeastern. Overall, steppe landscapes and open Pinus woodlands prevailed over many continental regions during the cold spells of the Late Pleistocene. The maintenance of a high phytodiversity during the glacials was linked to several refuge zones in the coastal shelves of the Mediterranean and intramountainous valleys. Northern Iberia, especially on coastal areas, was also patched with populations of tree species, and this is not only documented by palaeobotanical data (pollen and charcoal) but also postulated by phylogeographical models.     

Phylogeography of the widespread fern Lemmaphyllum in East Asia: species differentiation and population dynamics in response to change in climate and geography

Past climatic oscillations and complex geodynamic processes had tremendous effects on the current distributions of species in East Asia. Previous studies have revealed that spermatophytes experienced different demographic histories and survived in multiple refuges. However, very few studies involving ferns have been conducted over a large geographical area like East Asia. The monophyletic epiphytic fern genus Lemmaphyllum, which is composed of four species, is widespread in East Asia and offers a good model for exploring how geoclimatic oscillations influence the diversification and demographic history of fern species. We studied the phylogeography of Lemmaphyllum based on 115 populations using plastid sequences and ecological niche modeling. A total of 91 haplotypes were found in Lemmaphyllum. Molecular clock estimation revealed that speciation coincided with the three phases of the Qingzang Movement at beginning of the third uplift of Qinghai-Tibetan Plateau. The “Tanaka-Kaiyong Line” demarcated lineages within L. carnosum. The split of the mainland and island lineages of L. rostratum and L. carnosum var. microphyllum may have resulted from ancestral isolation whereby land-bridges acted as a “barrier” rather than as a “corridor” between mainland and island lineages. Multiple glacial refuges such as Sichuan Basin, Jinggangshan region, YGG region, HDM region, and the islands of the China East Sea during the LGM were revealed. The entities of Lemmaphyllum experienced species-specific demographic histories in response to the Pleistocene climate change. The case study of epiphytic ferns may provide evidences for understanding the migration of evergreen broad-leaf forest under climate oscillation.     

Geographical differentiation of the Euchiloglanis fish complex (Teleostei: Siluriformes) in the Hengduan Mountain Region,China: Phylogeographic evidence of altered drainage patterns

The uplift of the Tibetan Plateau caused significant ecogeographical changes that had a major impact on the exchange and isolation of regional fauna and flora. Furthermore, Pleistocene glacial oscillations were linked to temporal large‐scale landmass and drainage system reconfigurations near the Hengduan Mountain Region and might have facilitated speciation and promoted biodiversity in southwestern China. However, strong biotic evidence supporting this role is lacking. Here, we use the Euchiloglanis fish species complex as a model to demonstrate the compound effects of the Tibetan Plateau uplift and Pleistocene glacial oscillations on species formation in this region. The genetic structure and geographical differentiation of the Euchiloglanis complex in four river systems within the Hengduan Mountain Region were deduced using the cytochrome b (cyt b) gene and 10 microsatellite loci from 360 to 192 individuals, respectively. The results indicated that the populations were divided into four independently evolving lineages, in which the populations from the Qingyi River and Jinsha River formed two sub‐lineages. Phylogenetic relationships were structured by geographical isolation, especially near drainage systems. Divergence time estimation analyses showed that the Euchiloglanis complex diverged from its sister clade Pareuchiloglanis sinensis at around 1.3 Million years ago (Ma). Within the Euchiloglanis complex, the divergence time between the Dadu–Yalong and Jinsha–Qingyi River populations occurred at 1.0 Ma. This divergence time was in concordance with recent geological events, including the Kun‐Huang Movement (1.2–0.6 Ma) and the lag time (<2.0 Ma) of river incision in the Hengduan Mountain Region. Population expansion signals were detected from mismatched distribution analyses, and the expansion times were concurrent with Pleistocene glacier fluctuations. Therefore, current phylogeographic patterns of the Euchiloglanis fish complex in the Hengduan Mountain Region were influenced by the uplift event of the Tibetan Plateau and were subsequently altered by paleo‐river transitions during the late Pleistocene glacial oscillations.     

The phylogeographical history of the Iberian steppe plant Ferula loscosii (Apiaceae): a test of the abundant-centre hypothesis

The geology and climate of the western Mediterranean area were strongly modified during the Late Tertiary and the Quaternary. These geological and climatic events are thought to have induced changes in the population histories of plants in the Iberian Peninsula. However, fine-scale genetic spatial architecture across western Mediterranean steppe plant refugia has rarely been investigated. A population genetic analysis of amplified fragment length polymorphism variation was conducted on present-day, relict populations of Ferula loscosii (Apiaceae). This species exhibits high individual/population numbers in the middle Ebro river valley and, according to the hypothesis of an abundant-centre distribution, these northern populations might represent a long-standing/ancestral distribution centre. However, our results suggest that the decimated southern and central Iberian populations are more variable and structured than the northeastern ones, representing the likely vestiges of an ancestral distribution centre of the species. Phylogeographical analysis suggests that F. loscosii likely originated in southern Spain and then migrated towards the central and northeastern ranges, further supporting a Late Miocene southern-bound Mediterranean migratory way for its oriental steppe ancestors. In addition, different glacial-induced conditions affected the southern and northern steppe Iberian refugia during the Quaternary. The contrasting genetic homogeneity of the Ebro valley range populations compared to the southern Iberian ones possibly reflects more severe bottlenecks and subsequent genetic drift experienced by populations of the northern Iberia refugium during the Pleistocene, followed by successful postglacial expansion from only a few founder plants.     

Ecology and life history affect different aspects of the population structure of 27 high-alpine plants

A plant species' genetic population structure is the result of a complex combination of its life history, ecological preferences, position in the ecosystem and historical factors. As a result, many different statistical methods exist that measure different aspects of species' genetic structure. However, little is known about how these methods are interrelated and how they are related to a species' ecology and life history. In this study, we used the IntraBioDiv amplified fragment length polymorphisms data set from 27 high-alpine species to calculate eight genetic summary statistics that we jointly correlate to a set of six ecological and life-history traits. We found that there is a large amount of redundancy among the calculated summary statistics and that there is a significant association with the matrix of species traits. In a multivariate analysis, two main aspects of population structure were visible among the 27 species. The first aspect is related to the species' dispersal capacities and the second is most likely related to the species' postglacial recolonization of the Alps. Furthermore, we found that some summary statistics, most importantly Mantel's r and Jost's D, show different behaviour than expected based on theory. We therefore advise caution in drawing too strong conclusions from these statistics.     

Genetic assessment of population restorations of the critically endangered Silene hifacensis in the Iberian Peninsula

In order to preserve endangered plant populations and recover their evolutionary potential and ecological behavior, some restoration measures generally involve the reinforcement of the population size in existing natural populations or the reintroduction of new populations. Genetic monitoring of both natural and restored populations can provide an assessment of restoration protocol success in establishing populations that maintain levels of genetic diversity similar to those in natural populations. The highly threatened Spanish species Silene hifacensis (Caryophyllaceae) has only three natural reduced mainland populations in the Iberian Peninsula, following decline and extinction that occurred during the late 20th century. Preterit restoration strategies were essentially based on the implantation of new populations and reinforcement of certain existing populations using transplants mostly cultivated in greenhouses. In the present contribution, levels and patterns of genetic variability within natural and restored populations of Silene hifacensis were assessed using the molecular technique AFLP. Our results pointed out significant genetic diversity differences across the three existing natural populations though their population fragmentation and progressive loss of individuals have not had an impact on the global genetic diversity of this species. For restored populations, their levels of genetic diversity were similar and even higher than in natural populations. As a result, the past restoration protocols were successful in capturing similar and even higher levels of genetic diversity than those observed within natural pools. However, inbreeding processes have been detected for two restored populations. Finally, the main source of plant material for the long-time restored transplants appears to be the natural population of Cova de les Cendres. This study demonstrates, once again, how genetic markers are useful tools to be taken in consideration for endangered plant species conservation plans.     

Species persistence in northerly glacial refugia of Europe: a matter of chance or biogeographical traits?

Aim The southern European peninsulas (Iberian, Italian and Balkan) are considered to have been refugia for many European species of plants and animals during the climatic extremes of the Pleistocene ice ages. A number of recent studies (fossil and genetic), however, have provided evidence for full‐glacial survival of some species beyond these peninsulas. Here we explore the biogeographical traits of these species, and ask whether they possessed certain characteristics that enabled them to persist in more northerly refugia. Location Europe. Methods Fossil and genetic evidence for refugial localities of species that survived in Europe during the last full‐glacial was obtained from the literature (totalling 90 species: 34 woody plants and 56 vertebrates). Forty‐seven of these species (23 woody plants and 24 vertebrates) had fossil evidence, whereas the remaining 43 species (11 woody plants and 32 vertebrates) had only genetic evidence. All species were scored according to their present geographical distribution, habitat preference and life‐history traits. The species were classified on the basis of these traits using hierarchical cluster analysis. Analysis of similarities was used to examine differences in vertebrate and woody plant species groups that survived only in southerly refugia and those that also persisted in more northerly locations. Non‐metric multi‐dimensional scaling was used to examine patterns observed between and within groups. Results Results from our analysis of species with fossil and genetic evidence for survival in refugia reveal that species that survived only in southerly refugia were large‐seeded trees or thermophilous vertebrates. In contrast, species that had a full‐glacial distribution, including more northerly locations, were wind‐dispersed, habitat‐generalist trees with the ability to reproduce vegetatively, and habitat‐generalist mammals with present‐day northerly distributions. Main conclusions Analysis of the geographical distribution, habitat preference and life‐history traits of the species studied suggests that underlying biogeographical traits may have determined their response to Pleistocene glaciation. The traits most commonly found in present populations with a northerly distribution in Europe enabled the same species to exist much farther north than the southern European peninsulas during the full‐glacial. It is possible that many of these species are now in restricted populations, within the ‘warm‐stage’ refugia of the current interglacial. The northerly full‐glacial survival of a number of woody plants and vertebrate species has significant implications for understanding migration rates of these species in response to climate change. It also has important implications for understanding current patterns of genetic diversity of European species. We suggest that both fossil and genetic evidence should be used to identify and prioritize for conservation of refugial localities in southern and northern Europe.     

Comparative study of the phylogenetic structure in six Apodemus species (Mammalia,Rodentia) inferred from ISSR-PCR data

The utility of the Inter Simple Sequence Repeat-Polymerase Chain Reaction (ISSR-PCR) was explored in order to determine genetic variation in six species of the genus Apodemus (A. flavicollis, A. sylvaticus, A. uralensis, A. agrarius, A. mystacinus and A. epimelas) at the individual level, population level, in separate geographic samples and in the species as a whole. Six optimized primers produced highly reproducible and polymorphic DNA markers with 98.3% polymorphic bands on a total sample of 91 individuals from 32 localities in Europe and Asia. Moreover, each primer allowed for an exact diagnosis of each of the six Apodemus species and thus provides a simple and reliable tool for the hitherto problematic discrimination of species from the subgenus Sylvaemus. Genetic distances between species ranged from 0.079, among the closely related A. flavicollis and A. sylvaticus, to 0.203 between A. mystacinus and A. agrarius. A. flavicollis, A. uralensis and A. sylvaticus display a strong population substructure. The range of genetic distances among geographic samples within last two species reaches the values obtained for closely related species. ISSR markers proved to be a simple and reliable tool for species diagnosis, as well as for estimating genetic diversity below the species level and for closely related species, but they showed questionable reliability for larger genetic distances.     

Impact of ice ages on circumpolar molecular diversity: insights from an ecological key species

We address the impact of the ice age cycles on intraspecific cpDNA diversity, for the first time on the full circumboreal-circumarctic scale. The bird-dispersed bog bilberry (or arctic blueberry, Vaccinium uliginosum) is a key component of northern ecosystems and is here used to assess diversity in previously glaciated vs. unglaciated areas and the importance of Beringia as a refugium and source for interglacial expansion. Eighteen chloroplast DNA haplotypes were observed in and among 122 populations, grouping into three main lineages which probably diverged before, and thus were affected more or less independently by, all major glaciations. The boreal 'Amphi-Atlantic lineage' included one haplotype occurring throughout northern Europe and one occurring in eastern North America, suggesting expansion from at least two bottlenecked, glacial refugium populations. The boreal 'Beringian lineage' included seven haplotypes restricted to Beringia and the Pacific coast of USA. The 'Arctic-Alpine lineage' included nine haplotypes, one of them fully circumpolar. This lineage was unexpectedly diverse, also in previously glaciated areas, suggesting that it thrived on the vast tundras during the ice ages and recolonized deglaciated terrain over long distances. Its largest area of persistence during glaciations was probably situated in the north, stretching from Beringia and far into Eurasia, and it probably also survived the last glaciation in southern mountain ranges. Although Beringia apparently was important for the initial divergence and expansion of V. uliginosum as well as for continuous survival of both the Beringian and Arctic-Alpine lineages during all ice ages, this region played a minor role as a source for later interglacial expansions.     

Mitochondrial DNA sequence analysis reveals multiple Pleistocene glacial refugia for the Yellow‐spotted mountain newt,Neurergus derjugini (Caudata: Salamandridae) in the mid‐Zagros range in Iran and Iraq

Phylogeography is often used to investigate the effects of glacial cycles on current genetic structure of various plant and animal species. This approach can also identify the number and location of glacial refugia as well as the recolonization routes from those refugia to the current locations. To identify the location of glacial refugia of the Yellow‐spotted mountain newt, Neurergus derjugini, we employed phylogeography patterns and genetic variability of this species by analyzing partial ND4 sequences (867 bp) of 67 specimens from 15 sampling localities from the whole species range in Iran and Iraq. Phylogenetic trees concordant with haplotype networks showed a clear genetic structure among populations as three groups corresponding to the populations in the north, center, and south. Evolutionary ages of clades north and south ranging from 0.15 to 0.17 Myr, while the oldest clade is the central clade, corresponding to 0.32 Myr. Bayesian skyline plots of population size change through time show a relatively slight increase until about 25 kyr (around the last glacial maximum) and a decline of population size about 2.5 kyr. The presence of geographically structured clades in north, center, and south sections of the species range signifies the disjunct populations that have emerged in three different refugium. This study illustrates the importance of the effect of previous glacial cycles in shaping the genetic structure of mountain species in the Zagros range. These areas are important in terms of long‐term species persistence and therefore valuable areas for conservation of biodiversity.