Parallel plumage colour evolution and introgressive hybridization in wheatears

Genetic and phenotypic mosaics, in which various phenotypes and different genomic regions show discordant patterns of species or population divergence, offer unique opportunities to study the role of ancestral and introgressed genetic variation in phenotypic evolution. Here, we investigated the evolution of discordant phenotypic and genetic divergence in a monophyletic clade of four songbird taxa—pied wheatear (O. pleschanka), Cyprus wheatear (Oenanthe cypriaca), and western and eastern subspecies of black‐eared wheatear (O. h. hispanica and O. h. melanoleuca). Phenotypically, black back and neck sides distinguish pied and Cyprus wheatears from the white‐backed/necked black‐eared wheatears. Meanwhile, mitochondrial variation only distinguishes western black‐eared wheatear. In the absence of nuclear genetic data, and given frequent hybridization among eastern black‐eared and pied wheatear, it remains unclear whether introgression is responsible for discordance between mitochondrial divergence patterns and phenotypic similarities, or whether plumage coloration evolved in parallel. Multispecies coalescent analyses of about 20,000 SNPs obtained from RAD data mapped to a draft genome assembly resolve the species tree, provide evidence for the parallel evolution of colour phenotypes and establish western and eastern black‐eared wheatears as independent taxa that should be recognized as full species. The presence of the entire admixture spectrum in the Iranian hybrid zone and the detection of footprints of introgression from pied into eastern black‐eared wheatear beyond the hybrid zone despite strong geographic structure of ancestry proportions furthermore suggest a potential role for introgression in parallel plumage colour evolution. Our results support the importance of standing heterospecific and/or ancestral variation in phenotypic evolution.     

Comparative phylogeography and demographic history of two marsupials of the Atlantic Forest in eastern Brazil

The Atlantic Forest biodiversity hotspot in eastern South America has been the focus of several phylogeographic studies concerning relationships between populations and areas and how taxa respond to environmental changes. We infer and compare the demographic and biogeographic histories of two didelphid marsupial species, Gracilinanus microtarsus and Marmosops incanus, from the Atlantic Forest of eastern Brazil to determine how these species responded to environmental changes over time, using mitochondrial and nuclear DNA sequences. We found great intraspecific genetic divergence in both species and a strong geographic structure related to similar and spatially cohesive groups within each species. These groups are consistent with the same topographical barriers, such as mountains and river valleys. Intraspecific clades are very old, dating back to a period of tectonic activities in the Neogene (5.39–8.57 Mya). Changes in the environment over the last 7 million years lead to fairly concordant demographic changes in both marsupial species, including population expansion during the last glacial maximum (ca. 21,000 years ago) or last interglacial (ca. 120,000 years ago) or both. These results do not fit the Pleistocene refuge hypothesis as an explanation of the historical biogeography and diversification of both species in the Atlantic Forest, but are compatible with the Atlantis Forest hypothesis.     

Crater lake cichlids individually specialize along the benthic–limnetic axis

A common pattern of adaptive diversification in freshwater fishes is the repeated evolution of elongated open water (limnetic) species and high‐bodied shore (benthic) species from generalist ancestors. Studies on phenotype‐diet correlations have suggested that population‐wide individual specialization occurs at an early evolutionary and ecological stage of divergence and niche partitioning. This variable restricted niche use across individuals can provide the raw material for earliest stages of sympatric divergence. We investigated variation in morphology and diet as well as their correlations along the benthic‐limnetic axis in an extremely young Midas cichlid species, Amphilophus tolteca, endemic to the Nicaraguan crater lake Asososca Managua. We found that A. tolteca varied continuously in ecologically relevant traits such as body shape and lower pharyngeal jaw morphology. The correlation of these phenotypes with niche suggested that individuals are specialized along the benthic‐limnetic axis. No genetic differentiation within the crater lake was detected based on genotypes from 13 microsatellite loci. Overall, we found that individual specialization in this young crater lake species encompasses the limnetic‐ as well as the benthic macro‐habitat. Yet there is no evidence for any diversification within the species, making this a candidate system for studying what might be the early stages preceding sympatric divergence.     

Long‐term survival and diversification of an endemic Melitaea species in mountains of Iran and adjacent areas

Disjunct distribution patterns regularly resulted in the separation of different genetic lineages in glacial refugia. Recent patterns of survival and expansion have been often revealed by climatic niche modelling. We used the combination of genetic markers, geometric morphometry and climatic niche modelling to clear up the taxonomy and reconstruct the potential range of an endemic Iranian, taxonomically disputed Melitaea species in climatically different epochs. Our results show that this species (Melitaea abbas Gross and Ebert 1975, comb. n. = M. zagrosi Tóth and Varga syn. nova) is clearly separated from all taxa of the Melitaea phoebe species group and only occurs in Iran and Azerbaijan but were also predicted for some adjacent regions. Molecular markers and distribution modelling show consistently that this species should have had a long‐term survival in this area, and its range could have been slightly larger during the last glacial maximum than currently. Based on the studied molecular markers, three main groups in M. abbas can be recognized: those of steppic area of Azerbaijan, western Iran and northeastern Iran. Each group is characterized by own mitochondrial haplotypes, but also a high level of genetic diversity appears in the central part of the distribution area (Zagros Mts.).     

Refugia,colonization and diversification of an arid‐adapted bird: coincident patterns between genetic data and ecological niche modelling

Phylogeographical studies are common in boreal and temperate species from the Palaearctic, but scarce in arid‐adapted species. We used nuclear and mitochondrial markers to investigate phylogeography and to estimate chronology of colonization events of the trumpeter finch Bucanetes githagineus, an arid‐adapted bird. We used 271 samples from 16 populations, most of which were fresh samples but including some museum specimens. Microsatellite data showed no clear grouping according to the sampling locations. Microsatellite and mitochondrial data showed the clearest differentiation between Maghreb and Canary Islands and between Maghreb and Western Sahara. Mitochondrial data suggest differentiation between different Maghreb populations and among Maghreb and Near East populations, between Iberian Peninsula and Canary Islands, as well as between Western Sahara and Maghreb. Our coalescence analyses indicate that the trumpeter finch colonized North Africa during the humid Marine Isotope Stage 5 (MIS5) period of the Sahara region 125 000 years ago. We constructed an ecological niche model (ENM) to estimate the geographical distribution of climatically suitable habitats for the trumpeter finch. We tested whether changes in the species range in relation to glacial–interglacial cycles could be responsible for observed patterns of genetic diversity and structure. Modelling results matched with those from genetic data as the species' potential range increases in interglacial scenarios (in the present climatic scenario and during MIS5) and decreases in glacial climates (during the last glacial maximum, LGM, 21 000 years ago). Our results suggest that the trumpeter finch responded to Pleistocene climatic changes by expanding and contracting its range.     

The integration of multiple independent data reveals an unusual response to Pleistocene climatic changes in the hard tick Ixodes ricinus

In the last few years, improved analytical tools and the integration of genetic data with multiple sources of information have shown that temperate species exhibited more complex responses to ice ages than previously thought. In this study, we investigated how Pleistocene climatic changes affected the current distribution and genetic diversity of European populations of the tick Ixodes ricinus, an ectoparasite with high ecological plasticity. We first used mitochondrial and nuclear genetic markers to investigate the phylogeographic structure of the species and its Pleistocene history using coalescent‐based methods; then we used species distribution modelling to infer the climatic niche of the species at last glacial maximum; finally, we reviewed the literature on the I. ricinus hosts to identify the locations of their glacial refugia. Our results support the scenario that during the last glacial phase, I. ricinus never experienced a prolonged allopatric divergence in separate glacial refugia, but persisted with interconnected populations across Southern and Central Europe. The generalist behaviour in host choice of I. ricinus would have played a major role in maintaining connections between its populations. Although most of the hosts persisted in separate refugia, from the point of view of I. ricinus, they represented a continuity of ‘bridges’ among populations. Our study highlights the importance of species‐specific ecology in affecting responses to Pleistocene glacial–interglacial cycles. Together with other cases in Europe and elsewhere, it contributes to setting new hypotheses on how species with wide ecological plasticity coped with Pleistocene climatic changes.     

The Quaternary evolutionary history,potential distribution dynamics,and conservation implications for a Qinghai–Tibet Plateau endemic herbaceous perennial,Anisodus tanguticus (Solanaceae)

Various hypotheses have been proposed about the Quaternary evolutionary history of plant species on the Qinghai–Tibet Plateau (QTP), yet only a handful of studies have considered both population genetics and ecological niche context. In this study, we proposed and compared climate refugia hypotheses based on the phylogeographic pattern of Anisodus tanguticus (three plastid DNA fragments and nuclear internal transcribed spacer regions from 32 populations) and present and past species distribution models (SDMs). We detected six plastid haplotypes in two well‐differentiated lineages. Although all haplotypes could be found in its western (sampling) area, only haplotypes from one lineage occurred in its eastern area. Meanwhile, most genetic variations existed between populations (F_ST = 0.822). The SDMs during the last glacial maximum and last interglacial periods showed range fragmentation in the western area and significant range contraction in the eastern area, respectively, in comparison with current potential distribution. This species may have undergone intraspecific divergence during the early Quaternary, which may have been caused by survival in different refugia during the earliest known glacial in the QTP, rather than geological isolation due to orogenesis events. Subsequently, climate oscillations during the Quaternary resulted in a dynamic distribution range for this species as well as the distribution pattern of its plastid haplotypes and nuclear genotypes. The interglacial periods may have had a greater effect on A. tanguticus than the glacial periods. Most importantly, neither genetic data nor SDM alone can fully reveal the climate refugia history of this species. We also discuss the conservation implications for this important Tibetan folk medicine plant in light of these findings and SDMs under future climate models. Together, our results underline the necessity to combine phylogeographic and SDM approaches in future investigations of the Quaternary evolutionary history of species in topographically complex areas, such as the QTP.     

Late‐glacial recolonization and phylogeography of European red deer (Cervus elaphus L.)

The Pleistocene was an epoch of extreme climatic and environmental changes. How individual species responded to the repeated cycles of warm and cold stages is a major topic of debate. For the European fauna and flora, an expansion–contraction model has been suggested, whereby temperate species were restricted to southern refugia during glacial times and expanded northwards during interglacials, including the present interglacial (Holocene). Here, we test this model on the red deer (Cervus elaphus) a large and highly mobile herbivore, using both modern and ancient mitochondrial DNA from the entire European range of the species over the last c. 40 000 years. Our results indicate that this species was sensitive to the effects of climate change. Prior to the Last Glacial Maximum (LGM) haplogroups restricted today to South‐East Europe and Western Asia reached as far west as the UK. During the LGM, red deer was mainly restricted to southern refugia, in Iberia, the Balkans and possibly in Italy and South‐Western Asia. At the end of the LGM, red deer expanded from the Iberian refugium, to Central and Northern Europe, including the UK, Belgium, Scandinavia, Germany, Poland and Belarus. Ancient DNA data cannot rule out refugial survival of red deer in North‐West Europe through the LGM. Had such deer survived, though, they were replaced by deer migrating from Iberia at the end of the glacial. The Balkans served as a separate LGM refugium and were probably connected to Western Asia with genetic exchange between the two areas.     

Niche divergence promotes rapid diversification of East African sky island white‐eyes (Aves: Zosteropidae)

The Eastern Afromontane biodiversity hotspot composed of highly fragmented forested highlands (sky islands) harbours exceptional diversity and endemicity, particularly within birds. To explain their elevated diversity within this region, models founded on niche conservatism have been offered, although detailed phylogeographic studies are limited to a few avian lineages. Here, we focus on the recent songbird genus Zosterops, represented by montane and lowland members, to test the roles of niche conservatism versus niche divergence in the diversification and colonization of East Africa's sky islands. The species‐rich white‐eyes are a typically homogeneous family with an exceptional colonizing ability, but in contrast to their diversity on oceanic islands, continental diversity is considered depauperate and has been largely neglected. Molecular phylogenetic analysis of ~140 taxa reveals extensive polyphyly among different montane populations of Z. poliogastrus. These larger endemic birds are shown to be more closely related to taxa with divergent habitat types, altitudinal distributions and dispersal abilities than they are to populations of restricted endemics that occur in neighbouring montane forest fragments. This repeated transition between lowland and highland habitats over time demonstrate that diversification of the focal group is explained by niche divergence. Our results also highlight an underestimation of diversity compared to morphological studies that has implications for their taxonomy and conservation. Molecular dating suggests that the spatially extensive African radiation arose exceptionally rapidly (1–2.5 Ma) during the fluctuating Plio‐Pleistocene climate, which may have provided the primary driver for lineage diversification.     

Niche conservatism and phylogenetic clustering in a tribe of arid‐adapted marsupial mice,the Sminthopsini

The progressive expansion of the Australian arid zone during the last 20 Ma appears to have spurred the diversification of several families of plants, vertebrates and invertebrates, yet such taxonomic groups appear to show limited niche radiation. Here, we test whether speciation is associated with niche conservatism (constraints on ecological divergence) or niche divergence in a tribe of marsupial mice (Sminthopsini; 23 taxa) that includes the most speciose genus of living dasyurids, the sminthopsins. To that end, we integrated phylogenetic data with ecological niche modelling, to enable us to reconstruct the evolution of climatic suitability within Sminthopsini. Niche overlap among species was low‐moderate (but generally higher than expected given environmental background similarity), and the degree of phylogenetic clustering increased with aridity. Climatic niche reconstruction illustrates that there has been little apparent evolution of climatic tolerance within clades. Accordingly, climatic disparity tends to be accumulated among clades, suggesting considerable niche conservatism. Our results also indicate that evolution of climatic tolerances has been heterogeneous across different dimensions of climate (temperature vs. precipitation) and across phylogenetic clusters (Sminthopsis murina group vs. other groups). Although some results point to the existence of shifts in climatic niches during the speciation of sminthopsins, our study provides evidence for substantial phylogenetic niche conservatism in the group. We conclude that niche diversification had a low impact on the speciation of this tribe of small, but highly mobile marsupials.     

Molecular data and ecological niche modelling reveal the phylogeographic pattern of Cotinus coggygria (Anacardiaceae) in China's warm‐temperate zone

The phylogeography of common and widespread species helps to elucidate the history of local flora and vegetation. In this study, we selected Cotinus coggygria, a species widely distributed in China's warm‐temperate zone. One chloroplast DNA (cpDNA) region and ecological niche modelling were used to examine the phylogeographic pattern of C. coggygria. The cpDNA data revealed two phylogeographic groups (Southern and Northern) corresponding to the geographic regions. Divergence time analyses revealed that divergence of the two groups occurred at approximately 147,000 years before the present (BP), which coincided with the formation of the downstream area of the Yellow River, indicating that the Yellow River was a weak phylogeographic divide for C. coggygria. The molecular data and ecological niche modelling also indicated that C. coggyria did not experience population expansion after glaciations. This study thus supports the fact that Pleistocene glacial cycles only slightly affected C. coggygria, which survived in situ and occupied multiple localised glacial refugia during glaciations. This finding is contrary to the hypothesis of large‐scale range habitat contraction and retreat into a few main refugia.     

Diversification across the Palaearctic desert belt throughout the Pleistocene: phylogeographic history of the Houbara–Macqueen's bustard complex (Otididae: Chlamydotis) as revealed by mitochondrial DNA

Studies on the influence of Pleistocene climatic fluctuations and associated habitat changes on arid‐adapted bird species living in the Holarctic region are comparatively rare. In contrast to temperate species, the populations of arid‐adapted avian species might be characterized by low genetic differentiation because periods of population isolation were associated with the short interglacial periods, while population expansion events might have occurred during the longer glacial periods when steppe‐like vegetation might have been prevalent. In this study, we tested this hypothesis in a widespread arid‐adapted taxon of the Palaearctic desert belt, the Houbara–Macqueen's bustard complex. The later includes the Houbara bustard Chlamydotis undulata, comprising the North African subspecies Chlamydotis u. undulata and Chlamydotis u. fuertaventurae from the Canary Islands, and the Asian Macqueen's bustard Chlamydotis macqueenii. A long fragment (1042 bp) of the Cyt‐b gene was investigated in 39 representatives of the two species to assess phylogenetic and phylogeographic patterns, and demographic history and to compute divergence time estimates using a Bayesian relaxed molecular clock approach based on different coalescent priors. While the two species are genetically distinct, we found little intraspecific genetic differentiation. The divergence time of the two species falls within a period of extreme aridity at around 0.9 million years ago, which most likely resulted in an east–west vicariance along the Arabo‐Saharan deserts. Differentiation within Houbara and Macqueen's bustard occurred later during the Middle to Upper Pleistocene, and as we have predicted, periods of range expansion were associated to the last glacial period at least in the Macqueen's bustard.     

Phylogeographic structure in long‐tailed voles (Rodentia: Arvicolinae) belies the complex Pleistocene history of isolation,divergence, and recolonization of Northwest North America's fauna

Quaternary climate fluctuations restructured biodiversity across North American high latitudes through repeated episodes of range contraction, population isolation and divergence, and subsequent expansion. Identifying how species responded to changing environmental conditions not only allows us to explore the mode and tempo of evolution in northern taxa, but also provides a basis for forecasting future biotic response across the highly variable topography of western North America. Using a multilocus approach under a Bayesian coalescent framework, we investigated the phylogeography of a wide‐ranging mammal, the long‐tailed vole, Microtus longicaudus. We focused on populations along the North Pacific Coast to refine our understanding of diversification by exploring the potentially compounding roles of multiple glacial refugia and more recent fragmentation of an extensive coastal archipelago. Through a combination of genetic data and species distribution models (SDMs), we found that historical climate variability influenced contemporary genetic structure, with multiple isolated locations of persistence (refugia) producing multiple divergent lineages (Beringian or northern, southeast Alaska or coastal, and southern or continental) during glacial advances. These vole lineages all occur along the North Pacific Coast where the confluence of numerous independent lineages in other species has produced overlapping zones of secondary contact, collectively a suture zone. Finally, we detected high levels of neoendemism due to complex island geography that developed in the last 10,000 years with the rising sea levels of the Holocene.     

Postglacial climate changes and rise of three ecotypes of harbour porpoises,Phocoena phocoena,in western Palearctic waters

Despite no obvious barriers to gene flow in the marine realm, environmental variation and ecological specializations can lead to genetic differentiation in highly mobile predators. Here, we investigated the genetic structure of the harbour porpoise over the entire species distribution range in western Palearctic waters. Combined analyses of 10 microsatellite loci and a 5085 base‐pair portion of the mitochondrial genome revealed the existence of three ecotypes, equally divergent at the mitochondrial genome, distributed in the Black Sea (BS), the European continental shelf waters, and a previously overlooked ecotype in the upwelling zones of Iberia and Mauritania. Historical demographic inferences using approximate Bayesian computation (ABC) suggest that these ecotypes diverged during the last glacial maximum (c. 23–19 kilo‐years ago, kyrbp ). ABC supports the hypothesis that the BS and upwelling ecotypes share a more recent common ancestor (c. 14 kyrbp ) than either does with the European continental shelf ecotype (c. 28 kyrbp ), suggesting they probably descended from the extinct populations that once inhabited the Mediterranean during the glacial and post‐glacial period. We showed that the two Atlantic ecotypes established a narrow admixture zone in the Bay of Biscay during the last millennium, with highly asymmetric gene flow. This study highlights the impacts that climate change may have on the distribution and speciation process in pelagic predators and shows that allopatric divergence can occur in these highly mobile species and be a source of genetic diversity.     

Evolutionary history of the Persian Jird,Meriones persicus,based on genetics,species distribution modelling and morphometric data

The Persian Jird, Meriones persicus, is distributed from Eastern Anatolia to Afghanistan and western Pakistan. Six subspecies were described based on skull features and coat colours, but the validity of these subspecies is uncertain, and no molecular work has ever been conducted on this species. Iran appears to be a key geographical region in which to revise the systematic and evolutionary history of this species, because five of the six subspecies are present in this country. To evaluate the phylogeographical history and taxonomy of this species in Iran, we used a combination of genetic (cytochrome b gene sequences of 70 specimens) and geometric morphometric (2D landmarks on the ventral side of skull of 258 specimens) analyses. We also used ecological niche modelling to make inferences about the evolutionary history of these lineages. Our molecular data highlight the existence of four genetic lineages, but they only partly correspond to the previously described subspecies. Our molecular and morphometric data confirm the validity of M. p. rossicus and show that it has a wider geographical range than previously thought. M. p. gurganensis and M. p. baptistae are genetically very close. The skull of M. p. gurganensis is morphologically distinguishable from other subspecies. The subspecies M. p. persicus and M. p. baptistae are genetically distinct, but morphologically close. Meriones p. ambrosius is genetically close to M. p. persicus, and additional analyses with more specimens are needed to validate its subspecific status. The genetic structure observed in Iran seems to fit the topography and biogeography of the country and emphasize the role of the Abarkooh, Central and Lut deserts as barriers to gene flow. All intraspecific divergent events within the Persian Jird occurred during the last 1.4 My, suggesting that climatic changes probably trigger diversification within this species. Our genetic and species niche modelling results suggest that potential refugial areas persisted during glacial periods for this species in north‐western Zagros Mountains, north‐eastern Alborz Mountains and Kohrud Mountains.     

Divergent parasite infections in sympatric cichlid species in Lake Victoria

Parasitism has been proposed as a factor in host speciation, as an agent affecting coexistence of host species in species‐rich communities and as a driver of post‐speciation diversification. Young adaptive radiations of closely related host species of varying ecological and genomic differentiation provide interesting opportunities to explore interactions between patterns of parasitism, divergence and coexistence of sympatric host species. Here, we explored patterns in ectoparasitism in a community of 16 fully sympatric cichlid species at Makobe Island in Lake Victoria, a model system of vertebrate adaptive radiation. We asked whether host niche, host abundance or host genetic differentiation explains variation in infection patterns. We found significant differences in infections, the magnitude of which was weakly correlated with the extent of genomic divergence between the host species, but more strongly with the main ecological gradient, water depth. These effects were most evident with infections of Cichlidogyrus monogeneans, whereas the only host species with a strictly crevice‐dwelling niche, Pundamilia pundamilia, deviated from the general negative relationship between depth and parasitism. In accordance with the Janzen–Connell hypothesis, we also found that host abundance tended to be positively associated with infections in some parasite taxa. Data on the Pundamilia sister species pairs from three other islands with variable degrees of habitat (crevice) specialization suggested that the lower parasite abundance of P. pundamilia at Makobe could result from both habitat specialization and the evolution of specific resistance. Our results support influences of host genetic differentiation and host ecology in determining infections in this diverse community of sympatric cichlid species.     

Reevaluation of a classic phylogeographic barrier: new techniques reveal the influence of microgeographic climate variation on population divergence

We evaluated the mtDNA divergence and relationships within Geomys pinetis to assess the status of formerly recognized Geomys taxa. Additionally, we integrated new hypothesis‐based tests in ecological niche models (ENM) to provide greater insight into causes for divergence and potential barriers to gene flow in Southeastern United States (Alabama, Florida, and Georgia). Our DNA sequence dataset confirmed and strongly supported two distinct lineages within G. pinetis occurring east and west of the ARD. Divergence date estimates showed that eastern and western lineages diverged about 1.37 Ma (1.9 Ma–830 ka). Predicted distributions from ENMs were consistent with molecular data and defined each population east and west of the ARD with little overlap. Niche identity and background similarity tests were statistically significant suggesting that ENMs from eastern and western lineages are not identical or more similar than expected based on random localities drawn from the environmental background. ENMs also support the hypothesis that the ARD represents a ribbon of unsuitable climate between more suitable areas where these populations are distributed. The estimated age of divergence between eastern and western lineages of G. pinetis suggests that the divergence was driven by climatic conditions during Pleistocene glacial–interglacial cycles. The ARD at the contact zone of eastern and western lineages of G. pinetis forms a significant barrier promoting microgeographic isolation that helps maintain ecological and genetic divergence.     

Lineage-specific late pleistocene expansion of an endemic subtropical gossamer-wing damselfly, Euphaea formosa, in Taiwan

Background  

Pleistocene glacial oscillations have significantly affected the historical population dynamics of temperate taxa. However, the general effects of recent climatic changes on the evolutionary history and genetic structure of extant subtropical species remain poorly understood. In the present study, phylogeographic and historical demographic analyses based on mitochondrial and nuclear DNA sequences were used. The aim was to investigate whether Pleistocene climatic cycles, paleo-drainages or mountain vicariance of Taiwan shaped the evolutionary diversification of a subtropical gossamer-wing damselfly, Euphaea formosa.     

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.     

Glacial refugia,recolonization patterns and diversification forces in Alpine‐endemic Megabunus harvestmen

The Pleistocene climatic fluctuations had a huge impact on all life forms, and various hypotheses regarding the survival of organisms during glacial periods have been postulated. In the European Alps, evidence has been found in support of refugia outside the ice shield (massifs de refuge) acting as sources for postglacial recolonization of inner‐Alpine areas. In contrast, evidence for survival on nunataks, ice‐free areas above the glacier, remains scarce. Here, we combine multivariate genetic analyses with ecological niche models (ENMs) through multiple timescales to elucidate the history of Alpine Megabunus harvestmen throughout the ice ages, a genus that comprises eight high‐altitude endemics. ENMs suggest two types of refugia throughout the last glacial maximum, inner‐Alpine survival on nunataks for four species and peripheral refugia for further four species. In some geographic regions, the patterns of genetic variation are consistent with long‐distance dispersal out of massifs de refuge, repeatedly coupled with geographic parthenogenesis. In other regions, long‐term persistence in nunataks may dominate the patterns of genetic divergence. Overall, our results suggest that glacial cycles contributed to allopatric diversification in Alpine Megabunus, both within and at the margins of the ice shield. These findings exemplify the power of ENM projections coupled with genetic analyses to identify hypotheses about the position and the number of glacial refugia and thus to evaluate the role of Pleistocene glaciations in driving species‐specific responses of recolonization or persistence that may have contributed to observed patterns of biodiversity.