Alpine biogeography of Parnassian butterflies during Quaternary climate cycles in North America

Growth of alpine glaciers during the Pleistocene had profound effects on montane landscapes in North America and the organisms now inhabiting alpine ecosystems. Biogeography of this region has often been viewed as a system of sky islands despite the fact that species richness patterns deviate from a strict island biogeographic model. One explanation is that alpine species are not in equilibrium because of late Quaternary geographic range shifts. Genetic data can provide evidence of nonequilibrium dynamics and the distributional shifts that occur during glaciation events in alpine landscapes. Using mitochondrial and nuclear sequence data, we examine the evolutionary history of butterflies in the Parnassius phoebus complex. We test explicit, alternative models of the biogeographic history of Parnassius smintheus and Parnassius behrii , including an equilibrium island model, ancestral radiation and fragmentation, an expanding alpine archipelago and an alpine archipelago refuge model. Our results support the alpine archipelago refuge model, in which alpine butterflies undergo population contraction during glacial climates followed by population expansion during interglacial phases. While butterflies can disperse between distant mountain ranges during glacial periods, gene flow is rare. We find evidence of recent connectivity between California and Colorado, population expansion events following deglaciation ∼20 000 years B.P., and small population sizes during the last glacial period. An analysis of lineage splitting suggests that morphological differences in P. smintheus and P. behrii are the result of late Pleistocene divergence (∼48 000 years B.P.) with limited gene flow. Our results demonstrate that spatially complex and nonequilibrium population dynamics influence alpine diversity patterns.     

The homogenous genetic structure and inferred unique history of range shifts during the Pleistocene climatic oscillations of <Emphasis Type="Italic">Arcterica nana</Emphasis> (Maxim.) Makino (Ericaceae)

Previous phylogeographic studies of alpine plants in Japan have inferred that populations in central Honshu persisted during the Pleistocene climatic oscillations and suggested interglacial survival in high mountains. However, Arcterica nana (Maxim.) Makino (Ericaceae) exhibits a homogenous genetic structure throughout Japan and may therefore have a unique phylogeographic history. This inconsistency could have resulted from insufficient resolution of previously analyzed chloroplast DNA sequences. Therefore, we conducted a phylogeographic investigation based on amplified fragment length polymorphisms. Using 176 individuals from 21 populations, the relationships among individuals and populations were determined by principal coordinate analysis and a neighbor-joining tree, respectively. In addition, genetic differentiation was estimated using analysis of molecular variance and spatial autocorrelation analysis. These analyses demonstrate a homogenous structure throughout the entire Japanese range, supporting the previous cpDNA phylogeography. Although this genetic structure is inconsistent with those of other alpine plants, it is difficult to postulate that pre-existing genetic differentiation was swamped exclusively within A. nana. Therefore, this homogenous genetic structure may have been caused by the distinct history of populations of A. nana. Specifically, the southern-ward migration and the subsequent continuous populations enabled gene flow throughout the Japanese archipelago during the last glacial period. Thus, our data suggest that alpine plants in the Japanese archipelago did not always experience a shared distribution change following climatic oscillations.     

A range-wide genetic bottleneck overwhelms contemporary landscape factors and local abundance in shaping genetic patterns of an alpine butterfly (Lepidoptera: Pieridae: Colias behrii)

Spatial and environmental heterogeneity are major factors in structuring species distributions in alpine landscapes. These landscapes have also been affected by glacial advances and retreats, causing alpine taxa to undergo range shifts and demographic changes. These nonequilibrium population dynamics have the potential to obscure the effects of environmental factors on the distribution of genetic variation. Here, we investigate how demographic change and environmental factors influence genetic variation in the alpine butterfly Colias behrii. Data from 14 microsatellite loci provide evidence of bottlenecks in all population samples. We test several alternative models of demography using approximate Bayesian computation (ABC), with the results favouring a model in which a recent bottleneck precedes rapid population growth. Applying independent calibrations to microsatellite loci and a nuclear gene, we estimate that this bottleneck affected both northern and southern populations 531–281 years ago, coinciding with a period of global cooling. Using regression approaches, we attempt to separate the effects of population structure, geographical distance and landscape on patterns of population genetic differentiation. Only 40% of the variation in F_ST is explained by these models, with geographical distance and least‐cost distance among meadow patches selected as the best predictors. Various measures of genetic diversity within populations are also decoupled from estimates of local abundance and habitat patch characteristics. Our results demonstrate that demographic change can have a disproportionate influence on genetic diversity in alpine species, contrasting with other studies that suggest landscape features control contemporary demographic processes in high‐elevation environments.     

A comparative phylogeographic study reveals discordant evolutionary histories of alpine ground beetles (Coleoptera,Carabidae)

Taiwan, an island with three major mountain ranges, provides an ideal topography to study mountain–island effect on organisms that would be diversified in the isolation areas. Glaciations, however, might drive these organisms to lower elevations, causing gene flow among previously isolated populations. Two hypotheses have been proposed to depict the possible refugia for alpine organisms during glaciations. Nunatak hypothesis suggests that alpine species might have stayed in situ in high mountain areas during glaciations. Massif de refuge, on the other hand, proposes that alpine species might have migrated to lower ice‐free areas. By sampling five sympatric carabid species of Nebria and Leistus, and using two mitochondrial genes and two nuclear genes, we evaluated the mountain–island effect on alpine carabids and tested the two proposed hypotheses with comparative phylogeographic method. Results from the phylogenetic relationships, network analysis, lineage calibration, and genetic structure indicate that the deep divergence among populations in all L. smetanai, N. formosana, and N. niitakana was subjected to long‐term isolation, a phenomenon in agreement with the nunatak hypothesis. However, genetic admixture among populations of N. uenoiana and some populations of L. nokoensis complex suggests that gene flow occurred during glaciations, as a massif de refuge depicts. The speciation event in N. niitakana is estimated to have occurred before 1.89 million years ago (Mya), while differentiation among isolated populations in N. niitakana, N. formosana, L. smetanai, and L. nokoensis complex might have taken place during 0.65–1.65 Mya. While each of the alpine carabids arriving in Taiwan during different glaciation events acquired its evolutionary history, all of them had confronted the existing mountain ranges.     

Integrating environmental,molecular, and morphological data to unravel an ice‐age radiation of arctic‐alpine Campanula in western North America

Many arctic‐alpine plant genera have undergone speciation during the Quaternary. The bases for these radiations have been ascribed to geographic isolation, abiotic and biotic differences between populations, and/or hybridization and polyploidization. The Cordilleran Campanula L. (Campanulaceae Juss.), a monophyletic clade of mostly endemic arctic‐alpine taxa from western North America, experienced a recent and rapid radiation. We set out to unravel the factors that likely influenced speciation in this group. To do so, we integrated environmental, genetic, and morphological datasets, tested biogeographic hypotheses, and analyzed the potential consequences of the various factors on the evolutionary history of the clade. We created paleodistribution models to identify potential Pleistocene refugia for the clade and estimated niche space for individual taxa using geographic and climatic data. Using 11 nuclear loci, we reconstructed a species tree and tested biogeographic hypotheses derived from the paleodistribution models. Finally, we tested 28 morphological characters, including floral, vegetative, and seed characteristics, for their capacity to differentiate taxa. Our results show that the combined effect of Quaternary climatic variation, isolation among differing environments in the mountains in western North America, and biotic factors influencing floral morphology contributed to speciation in this group during the mid‐Pleistocene. Furthermore, our biogeographic analyses uncovered asynchronous consequences of interglacial and glacial periods for the timing of refugial isolation within the southern and northwestern mountains, respectively. These findings have broad implications for understanding the processes promoting speciation in arctic‐alpine plants and the rise of numerous endemic taxa across the region.     

Historical biogeography of two alpine butterflies in the Rocky Mountains: broad-scale concordance and local-scale discordance

Aim We inferred the phylogeography of the alpine butterfly Colias meadii Edwards (Pieridae) and compared its genetic structure with that of another high elevation, co‐distributed butterfly, Parnassius smintheus Doubleday (Papilionidae), to test if the two Rocky Mountain butterflies responded similarly to the palaeoclimatic cycles of the Quaternary. Location Specimens were collected from 18 alpine sites in the Rocky Mountains of North America, from southern Colorado to northern Montana. Methods We sequenced 867 and 789 nucleotides of cytochrome oxidase I from an average of 19 and 20 individuals for C. meadii and P. smintheus, respectively, from each of the same 18 localities. From the sequence data, we calculated measures of genetic diversity within each population (H, θ), genetic divergence among populations (F_ST), and tested for geographic structure through an analysis of molecular variance (amova ). Population estimates were compared against latitude and between species using a variety of statistical tests. Furthermore, nested clade analysis was implemented to infer historic events underlying the geographic distribution of genetic variation in each species. Then, we compared the number of inferred population events between species using a nonparametric Spearman's rank correlation test. Finally, we ran coalescent simulations on each species’ genealogy to test whether the two species of Lepidoptera fit the same model of population divergence. Results Our analyses revealed that: (1) measures of within‐population diversity were not correlated with latitude for either species, (2) within‐site diversity was not correlated between species, (3) within a species, nearly all populations were genetically isolated, (4) both species exhibited significant and nearly identical partitioning of genetic variation at all hierarchical levels of the amova , including a strong break between populations across the Wyoming Basin, (5) both species experienced similar cycles of expansion and contraction, although fewer were inferred for C. meadii, and (6) data from both species fit a model of three refugia diverging during the Pleistocene. Main conclusions While our findings supported a shared response of the two butterfly species to historic climate change across coarse spatial scales, a common pattern was not evident at finer spatial and temporal scales. The shared demographic history of the two species is consistent with an expanding–contracting archipelago model, suggesting that populations persisted across the geographic range throughout the climate cycles, experiencing isolation on ‘sky islands’ during interglacial periods and becoming connected as they migrated down‐slope during cool, wet climates.     

Population fragmentation causes randomly fixed genotypes in populations of <Emphasis Type="Italic">Arabidopsis kamchatica</Emphasis> in the Japanese Archipelago

Populations of arctic alpine plants likely disappeared and re-colonised several times at the southern edge of their distributions during glacial and interglacial cycles throughout the Quaternary. Range shift and population fragmentation after a glacial period would affect the genetic structure of such plants in southernmost populations. We aimed to elucidate how climatic oscillations influenced the population subsistence of alpine plants in the Japanese Archipelago as one of the southernmost populations, by inferring the genetic structure of Arabidopsis kamchatica subsp. kamchatica and the intraspecific littoral taxon, subsp. kawasakiana. We identified genotypes based on the haplotypes of five nuclear genes and two chloroplast DNA spacers for 164 individuals from 24 populations. Most populations harboured only one private genotype, whereas few polymorphisms were found in each population. Two genetic genealogies were found, suggesting that northern Japanese populations of alpine subsp. kamchatica, subsp. kawasakiana and the northerly subsp. kamchatica in eastern Russia and Alaska clustered and differentiated from populations in central Honshu, western Japan and Taiwan. During climatic oscillations, the genetic structure of extant southernmost populations would have been shaped by strong genetic drift under population fragmentation and randomly fixed to a single genotype among their ancestral polymorphisms.     

Long‐term in situ persistence of biodiversity in tropical sky islands revealed by landscape genomics

Tropical mountains are areas of high species richness and endemism. Two historical phenomena may have contributed to this: (i) fragmentation and isolation of habitats may have promoted the genetic differentiation of populations and increased the possibility of allopatric divergence and speciation and (ii) the mountain areas may have allowed long‐term population persistence during global climate fluctuations. These two phenomena have been studied using either species occurrence data or estimating species divergence times. However, only few studies have used intraspecific genetic data to analyse the mechanisms by which endemism may emerge at the microevolutionary scale. Here, we use landscape analysis of genomic SNP data sampled from two high‐elevation plant species from an archipelago of tropical sky islands (the Trans‐Mexican Volcanic Belt) to test for population genetic differentiation, synchronous demographic changes and habitat persistence. We show that genetic differentiation can be explained by the degree of glacial habitat connectivity among mountains and that mountains have facilitated the persistence of populations throughout glacial/interglacial cycles. Our results support the ongoing role of tropical mountains as cradles for biodiversity by uncovering cryptic differentiation and limits to gene flow.     

Speciation with gene flow: Evidence from a complex of alpine butterflies (Coenonympha,Satyridae)

Until complete reproductive isolation is achieved, the extent of differentiation between two diverging lineages is the result of a dynamic equilibrium between genetic isolation and mixing. This is especially true for hybrid taxa, for which the degree of isolation in regard to their parental species is decisive in their capacity to rise as a new and stable entity. In this work, we explored the past and current patterns of hybridization and divergence within a complex of closely related butterflies in the genus Coenonympha in which two alpine species, C. darwiniana and C. macromma, have been shown to result from hybridization between the also alpine C. gardetta and the lowland C. arcania. By testing alternative scenarios of divergence among species, we show that gene flow has been uninterrupted throughout the speciation process, although leading to different degrees of current genetic isolation between species in contact zones depending on the pair considered. Nonetheless, at broader geographic scale, analyses reveal a clear genetic differentiation between hybrid lineages and their parental species, pointing out to an advanced stage of the hybrid speciation process. Finally, the positive correlation observed between ecological divergence and genetic isolation among these butterflies suggests a potential role for ecological drivers during their speciation processes.     

Post-glacial range fragmentation is responsible for the current distribution of Potentilla matsumurae Th. Wolf (Rosaceae) in the Japanese archipelago

Aim We aimed to elucidate how the current geographic distribution of alpine plants in the Japanese archipelago was shaped during Quaternary climatic oscillations, using Potentilla matsumurae as a case study. According to previous phylogeographic studies, post‐glacial range fragmentation (vicariance scenario) and stepwise migration (dispersal scenario) are both possible. We thus aimed to assess which scenario is more probable for the distribution changes of alpine plants in the Japanese archipelago. Location The alpine zone in the Japanese archipelago. Methods Using amplified fragment length polymorphism we determined the genotype of 161 individuals of P. matsumurae from 22 populations. Relationships among individuals and populations were examined using principal coordinates analysis and a neighbour‐joining (NJ) tree, respectively. To examine the genetic population structure, we performed analysis of molecular variance (amova ) and structure analysis. Results Differentiation between central Honshu and northern Japan was not very strong based on the principal coordinates analysis among individuals, the NJ tree of populations (59% bootstrap support), or amova (12% of genetic variation). Moreover, structure analysis did not detect clear geographic differentiation across populations. Although the populations in central Honshu were structured geographically (Mantel test: r = 0.45, P < 0.005; NJ tree), those in northern Japan did not exhibit geographic structure regardless of geographic distance (Mantel test: r = 0.26, P = 0.03; NJ tree). Population relationships in the NJ tree did not always reflect the geographic location. Main conclusions The current geographic structure of P. matsumurae could not be explained by stepwise migration. This suggests that a single continuous distribution during the last glacial period was later fragmented, perhaps by recovering forest, during the post‐glacial period, resulting in the current distribution and phylogeographic structure of P. matsumurae. Our data support the vicariance scenario.     

Importance of demographic history for phylogeographic inference on the arctic–alpine plant Phyllodoce caerulea in East Asia

Arctic–alpine plants have enormous ranges in the Northern Hemisphere. Phylogeographic studies have provided insights into their glacial survival as well as their postglacial colonization history. However, our understanding of the population dynamics of disjunct alpine populations in temperate regions remains limited. During Pleistocene cold periods, alpine populations of arctic–alpine species in East Asia were either connected to an ice-free Beringia refugium or they persisted with prolonged isolation after their establishment. To estimate which of these scenarios is more likely, we elucidated the genetic structure of Phyllodoce caerulea (Ericaceae) in Beringia and northern Japan, East Asia. Sequence variation in multiple nuclear loci revealed that P. caerulea can be distinguished into northern and southern groups. A demographic analysis demonstrated that the north–south divergence did not predate the last glacial period and detected introgression from Phyllodoce aleutica, relative widely distributed in East Asia, exclusively into the southern group. Therefore, although there has been genetic divergence between northern Japan and Beringia in P. caerulea, the divergence is unlikely to have resulted from their prolonged geographic separation throughout several cycles of glacial and interglacial periods. Instead, our study suggests that the introgression contributed to the genetic divergence of P. caerulea and that the range of P. caerulea was plausibly connected between northern Japan and Beringia during the last glacial period. Overall, our study not only provides a biogeographic insight into alpine populations of arctic–alpine plants in East Asia but also emphasizes the importance of careful interpretation of genetic structure for inferring phylogeographic history.     

Gene flow from multiple sources maintains high genetic diversity and stable population history of Common Moorhen Gallinula chloropus in China

In this study, we used mitochondrial control sequences and microsatellite data from 231 Common Moorhen Gallinula chloropus individuals sampled from 19 sites in China to analyse their genetic structure and evolutionary history. High genetic diversity was found for all populations, although microsatellite analysis showed that the genetic diversity in non‐migratory populations was significantly higher than in migratory populations. High gene flow occurred between neighbouring populations, although long‐distance gene flow also occurred. The Huazhong population was the single greatest genetic source for other populations. High gene flow probably led to the shallow genetic structure that we observed. Demographic expansion was found in migratory populations, non‐migratory populations and with all individuals combined. The expansion time for all populations combined was estimated to be 221 000 years ago. The Common Moorhen population grew rapidly during the interglacial before the last glacial maximum (LGM), then remained generally stable from the LGM to the present.     

Repeated landmass reformation limits diversification in the widespread littoral zone mosquito Anopheles sundaicus sensu lato in the Indo‐Oriental Region

Southeast Asia harbours abundant biodiversity, hypothesized to have been generated by Pliocene and Pleistocene climatic and environmental change. Vicariance between the island of Borneo, the remaining Indonesian archipelago and mainland Southeast Asia caused by elevated sea levels during interglacial periods has been proposed to lead to diversification in the littoral zone mosquito Anopheles (Cellia) sundaicus (Rodenwaldt) sensu lato. To test this biogeographical hypothesis, we inferred the population history and assessed gene flow of A. sundaicus s.l. sampled from 18 populations across its pan‐Asian species range, using sequences from mitochondrial cytochrome c oxidase subunit 1 (CO1), the internal transcribed spacer 2 (ITS2) and the mannose phosphate isomerase (Mpi) gene. A hypothesis of ecological speciation for A. sundaicus involving divergent adaptation to brackish and freshwater larval habitats was also previously proposed, based on a deficiency of heterozygotes for Mpi allozyme alleles in sympatry. This hypothesis was not supported by Mpi sequence data, which exhibited no fixed differences between brackish and freshwater larval habitats. Mpi and CO1 supported the presence of up to eight genetically distinct population groupings. Counter to the hypothesis of three allopatric species, divergence was often no greater between Borneo, Sumatra/Java and the Southeast Asian mainland than it was between genetic groupings within these landmasses. An isolation‐with‐migration (IM) model indicates recurrent gene flow between the current major landmasses. Such gene flow would have been possible during glacial periods when the current landmasses merged, presenting opportunities for dispersal along expanding and contracting coastlines. Consequently, Pleistocene climatic variation has proved a homogenizing, rather than diversifying, force for A. sundaicus diversity.     

Genetic consequences of interglacial isolation in a steppe bird

In response to climate changes that have occurred during Pleistocene glacial cycles, taxa associated to steppe vegetation might have followed a pattern of historical evolution in which isolation and fragmentation of populations occurred during the short interglacials and expansion events occurred during the long glacial periods, in contrast to the pattern described for temperate species. Here, we use molecular genetic data to evaluate this idea in a steppe bird with Palaearctic distribution, the little bustard (Tetrax tetrax). Overall, extremely low genetic diversity and differentiation was observed among eight little bustard populations distributed in Spain and France. Mismatch distribution analyses showed that most little bustard populations expanded during cooling periods previous to, and just after, the last interglacial period (127,000-111,000 years before present), when steppe habitats were widespread across Europe. Coalescent-based methods suggested that glacial expansions have resulted in substantial admixture in Western Europe due to the existence of different interglacial refugia. Our results are consistent with a model of evolution and genetic consequences of Pleistocene cycles with low between-population genetic differentiation as a result of short-term isolation periods during interglacials and long-term exchange during glacial periods.     

Historic cycles of fragmentation and expansion in Parnassius smintheus (papilionidae) inferred using mitochondrial DNA

Climate oscillations of the Quaternary drove the repeated expansion and contraction of ecosystems. Alpine organisms were probably isolated in sky island refugia during warm interglacials, such as now, and expanded their range by migrating down-slope during glacial periods. We used population genetic and phylogenetic approaches to infer how paleoclimatic events influenced the distribution of genetic variation in the predominantly alpine butterfly Parnassius smintheus. We sequenced a 789 bp region of cytochrome oxidase I for 385 individuals from 20 locations throughout the Rocky Mountains, ranging from southern Colorado to northern Montana. Analyses revealed at lease two centers of diversity in the northern and southern Rocky Mountains and strong population structure. Nested clade analysis suggested that the species experienced repeated cycles of population expansion and fragmentation. The estimated ages of these events, assuming a molecular clock, corresponded with paleoclimatic data on habitat expansion and contraction over the past 400,000 years. We propose that alpine butterflies persisted in an archipelago of isolated sky islands during interglacials and that populations expanded and became more connected during cold glacial periods. An archipelago model implies that the effects of genetic drift and selection varied among populations, depending on their latitude, area, and local environment. Alpine organisms are sensitive indicators of climate change and their history can be used to predict how high-elevation ecosystems might respond to further climate warming.     

ANALYSIS OF MECHANISMS OF MICROEVOLUTIONARY CHANGE IN CEPPHUS GUILLEMOTS USING PATTERNS OF CONTROL REGION VARIATION

We surveyed population-level sequence variation in part of the mitochondrial control region for three species including eight subspecies of Cepphus guillemots (Charadriiformes: Alcidae) to test specific predictions about mechanisms of population differentiation. We found that sequences of spectacled guillemots (C. carbo) were more closely related to those of pigeon guillemots (C. columba; both found in the Pacific Ocean) than to those of black guillemots (C. grylle; Arctic and Atlantic Oceans), despite dissimilarities in plumage between spectacled guillemots and the other species. Distributions of species and timing of divergence events suggest that speciation involved allopatric and microallopatric populations isolated by Pleistocene glaciers. Control region sequences were significantly differentiated among populations within species and suggest that gene flow is low; however, populations are probably not in genetic equilibrium, so these results probably reflect historical isolation of colonies. In contrast, phylogenetic relationships among sequences within species were poorly resolved, probably because of a combination of incomplete lineage sorting and contemporary gene flow. Indices of genetic diversity provided no suggestion of recent bottlenecks in most populations, although two populations apparently underwent recent severe bottlenecks. Genetic divergence among populations was not correlated with geographic distance, which argues against isolation by distance. Results of these analyses, combined with breeding distributions and timing of divergence events, suggest that populations diverged during isolation in glacial refugia. Our results are consistent with earlier hypotheses posed by Storer and Udvardy.     

Evolution of a hybrid zone of two willow species (Salix L.) in the European Alps analyzed by RAD-seq and morphometrics

Natural hybridization of plants can result in many outcomes with several evolutionary consequences, such as hybrid speciation and introgression. Natural hybrid zones can arise in mountain systems as a result of fluctuating climate during the exchange of glacial and interglacial periods, where species retract and expand their territories, resulting in secondary contacts. Willows are a large genus of woody plants with an immense capability of interspecific crossing. In this study, the sympatric area of two diploid sister species, S. foetida and S. waldsteiniana in the eastern European Alps, was investigated to study the genomic structure of populations within and outside their contact zone and to analyze congruence of morphological phenotypes with genetic data. Eleven populations of the two species were sampled across the Alps and examined using phylogenetic network and population genetic structure analyses of RAD Seq data and morphometric analyses of leaves. The results showed that a homoploid hybrid zone between the two species was established within their sympatric area. Patterns of genetic admixture in homoploid hybrids indicated introgression with asymmetric backcrossing to not only one of the parental species but also one hybrid population forming a separate lineage. The lack of F1 hybrids indicated a long-term persistence of the hybrid populations. Insignificant isolation by distance suggests that gene flow can act over large geographical scales. Morphometric characteristics of hybrids supported the molecular data and clearly separated populations of the parental species, but showed intermediacy in the hybrid zone populations with a bias toward S. waldsteiniana. The homoploid hybrid zone might have been established via secondary contact hybridization, and its establishment was fostered by the low genetic divergence of parental species and a lack of strong intrinsic crossing barriers. Incomplete ecological separation and the ability of long-distance dispersal of willows could have contributed to the spatial expansion of the hybrid zone.     

Interglacial genetic diversification of Moussonia deppeana (Gesneriaceae), a hummingbird‐pollinated,cloud forest shrub in northern Mesoamerica

Recent empirical work on cloud forest‐adapted species supports the role of both old divergences across major geographical areas and more recent divergences attributed to Pleistocene climate changes. The shrub Moussonia deppeana is distributed in northern Mesoamerica, with geographically disjunct populations. Based on sampling throughout the species range and employing plastid and nuclear markers, we (i) test whether the fragmented distribution is correlated with main evolutionary lineages, (ii) reconstruct its phylogeographical history to infer the history of cloud forest in northern Mesoamerica and (iii) evaluate a set of refugia/vicariance scenarios for the region and demographic patterns of the populations whose ranges expanded and tracked cloud forest conditions during the Last Glacial Maximum. We found a deep evolutionary split in M. deppeana about 6–3 Ma, which could be consistent with a Pliocene divergence. Comparison of variation in plastid and nuclear markers revealed several lineages mostly congruent with their isolated geographical distribution and restricted gene flow among groups. Results of species distribution modelling and coalescent simulations fit a model of multiple refugia diverging during interglacial cycles. The demographic history of M. deppeana is not consistent with an expanding–contracting cloud forest archipelago model during the Last Glacial Maximum. Instead, our data suggest that populations persisted across the geographical range throughout the glacial cycles, and experienced isolation and divergence during interglacial periods.     

Potential warm-stage microrefugia for alpine plants: Feedback between geomorphological and biological processes

During interglacial stages, microrefugia are sites that support locally favorable climates within larger areas with unfavorable warmer climates. Despite recent theoretical representations of microrefugia, an appropriate ecological characterization is still lacking, mostly for warm periods. Across mountain/alpine areas, cold-adapted plant species could adopt different strategies to manage the effects of climate warming: (A) migration toward higher elevations and summits; (B) in situ resilience of communities and species populations within microrefugia; and C) adaptation and evolution by genetic differentiation. This review aims to distinguish and characterize from an ecological perspective glacial, nival, periglacial and composite landforms and deposits that may function as potential microrefugia during interglacial warm periods.We conducted a literature screening related to the geomorphological processes and landforms associated with vegetation and plant communities in alpine/mountain environments of Europe. They include glacial deposits rock glaciers, debris-covered glaciers, composite cones and channels. In Alpine regions, geomorphologic niches that constantly maintain cold-air pooling and temperature inversions are the main candidates for microrefugia. Within such microrefugia, microhabitat diversity modulates the responses of plants to disturbances caused by geomorphologic processes and supports their aptitude for surviving under extreme conditions on unstable surfaces in isolated patches. Currently, European marginal mountain chains may be considered as examples of macrorefugia where relict boreo-alpine species persist within peculiar geomorphological niches that act as microrefugia.This review contributes to identifying potential warm-stage microrefugia areas across alpine and mountain regions and determining certain landforms that play or may play such role under global-change scenarios. The occurrence of warm-stage microrefugia within these locations may be of great importance for the modeling of future distributions of species and assessing the risk of extinction for alpine species. Microrefugia may have important implications in micro-evolutionary processes that occur across alternating climatic phases.     

Population history and genetic structure of a circumpolar species: the arctic fox

The circumpolar arctic fox Alopex lagopus thrives in cold climates and has a high migration rate involving long-distance movements. Thus, it differs from many temperate taxa that were subjected to cyclical restriction in glacial refugia during the Ice Ages. We investigated population history and genetic structure through mitochondrial control region variation in 191 arctic foxes from throughout the arctic. Several haplotypes had a Holarctic distribution and no phylogeographical structure was found. Furthermore, there was no difference in haplotype diversity between populations inhabiting previously glaciated and unglaciated regions. This suggests current gene flow among the studied populations, with the exception of those in Iceland, which is surrounded by year-round open water. Arctic foxes have often been separated into two ecotypes: 'lemming' and 'coastal'. An analysis of molecular variance suggested particularly high gene flow among populations of the 'lemming' ecotype. This could be explained by their higher migration rate and reduced fitness in migrants between ecotypes. A mismatch analysis indicated a sudden expansion in population size around 118 000 BP, which coincides with the last interglacial. We propose that glacial cycles affected the arctic fox in a way opposite to their effect on temperate species, with interglacials leading to short-term isolation in northern refugia.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 79–89.