Less Pollen-Mediated Gene Flow for More Signatures of Glacial Lineages: Congruent Evidence from Balsam Fir cpDNA and mtDNA for Multiple Refugia in Eastern and Central North America

The phylogeographic structure and postglacial history of balsam fir (Abies balsamea), a transcontinental North American boreal conifer, was inferred using mitochondrial DNA (mtDNA) and chloroplast DNA (cpDNA) markers. Genetic structure among 107 populations (mtDNA data) and 75 populations (cpDNA data) was analyzed using Bayesian and genetic distance approaches. Population differentiation was high for mtDNA (dispersed by seeds only), but also for cpDNA (dispersed by seeds and pollen), indicating that pollen gene flow is more restricted in balsam fir than in other boreal conifers. Low cpDNA gene flow in balsam fir may relate to low pollen production due to the inherent biology of the species and populations being decimated by recurrent spruce budworm epidemics, and/or to low dispersal of pollen grains due to their peculiar structural properties. Accordingly, a phylogeographic structure was detected using both mtDNA and cpDNA markers and population structure analyses supported the existence of at least five genetically distinct glacial lineages in central and eastern North America. Four of these would originate from glacial refugia located south of the Laurentide ice sheet, while the last one would have persisted in the northern Labrador region. As expected due to reduced pollen-mediated gene flow, congruence between the geographic distribution of mtDNA and cpDNA lineages was higher than in other North American conifers. However, concordance was not complete, reflecting that restricted but nonetheless detectable cpDNA gene flow among glacial lineages occurred during the Holocene. As a result, new cpDNA and mtDNA genome combinations indicative of cytoplasmic genome capture were observed.     

Phylogeography of Douglas‐fir based on mitochondrial and chloroplast DNA sequences: testing hypotheses from the fossil record

The integration of fossil and molecular data can provide a synthetic understanding of the ecological and evolutionary history of an organism. We analysed range‐wide maternally inherited mitochondrial DNA and paternally inherited chloroplast DNA sequence data with coalescent simulations and traditional population genetic methods to test hypotheses of population divergence generated from the fossil record of Douglas‐fir (Pseudotsuga menziesii), an ecologically and economically important western North American conifer. Specifically, we tested (i) the hypothesis that the Pliocene orogeny of the Cascades and Sierra Nevada caused the divergence of coastal and Rocky Mountain Douglas‐fir varieties; and (ii) the hypothesis that multiple glacial refugia existed on the coast and in the Rocky Mountains. We found that Douglas‐fir varieties diverged about 2.11 Ma (4.37 Ma–755 ka), which could be consistent with a Pliocene divergence. Rocky Mountain Douglas‐fir probably resided in three or more glacial refugia. More variable molecular markers would be required to detect the two coastal refugia suggested in the fossil record. Comparison of mitochondrial DNA and chloroplast DNA variation revealed that gene flow via pollen linked populations isolated from seed exchange. Postglacial colonization of Canada from coastal and Rocky Mountain refugia near the ice margin at the Last Glacial Maximum produced a wide hybrid zone among varieties that formed almost exclusively by pollen exchange and chloroplast DNA introgression, not seed exchange. Postglacial migration rates were 50–165 m/year, insufficient to track projected 21st century warming in some regions. Although fossil and genetic data largely agree, each provides unique insights.     

Multiple glacial refugia for cool‐temperate deciduous trees in northern East Asia: the Mongolian oak as a case study

In East Asia, temperate forests are predicted to have retracted southward to c. 30° N during the last glacial maximum (LGM) based on fossil pollen data, whereas phylogeographic studies have often suggested glacial in situ survival of cool‐temperate deciduous trees in their modern northern ranges. Here we report a study of the genetic diversity and structure of 29 natural Mongolian oak (Quercus mongolica) populations using 19 nuclear simple sequence repeat (nSSR) loci and four chloroplast DNA fragments. Bayesian clustering analysis with nSSRs revealed five groups, which were inferred by approximate Bayesian computation (ABC) to have diverged in multiple refugia through multiple glacial–interglacial cycles. Analysis of chloroplast DNA variation revealed four lineages that were largely but incompletely geographically disjunct. Ecological niche modelling (ENMs) indicated a southward range shift of the oak's distribution at the LGM, although high suitability scores were also evident in the Changbai Mts. (Northeast China), the Korean Peninsula, areas surrounding the Bohai Sea, and along the coast of the Russian Far East. In addition, endemic chloroplast DNA haplotypes and nuclear lineages occurred in high‐latitude northern areas where the ENM predicted no suitable habitat. The combined evidence from nuclear and chloroplast DNA, and the results of the ENM clearly demonstrate that multiple northern refugia, including cryptic ones, were maintained across the current distributional range of the Mongolian oak during the LGM or earlier glacial periods. Though spatially limited, postglacial expansions from these refugia have led to a pattern of decreased genetic diversity with increasing latitude.     

When new human‐modified habitats favour the expansion of an amphibian pioneer species: Evolutionary history of the natterjack toad (Bufo calamita) in a coal basin

Human activities affect microevolutionary dynamics by inducing environmental changes. In particular, land cover conversion and loss of native habitats decrease genetic diversity and jeopardize the adaptive ability of populations. Nonetheless, new anthropogenic habitats can also promote the successful establishment of emblematic pioneer species. We investigated this issue by examining the population genetic features and evolutionary history of the natterjack toad (Bufo [Epidalea] calamita) in northern France, where populations can be found in native coastal habitats and coalfield habitats shaped by European industrial history, along with an additional set of European populations located outside this focal area. We predicted contrasting patterns of genetic structure, with newly settled coalfield populations departing from migration–drift equilibrium. As expected, coalfield populations showed a mosaic of genetically divergent populations with short‐range patterns of gene flow, and native coastal populations indicated an equilibrium state with an isolation‐by‐distance pattern suggestive of postglacial range expansion. However, coalfield populations exhibited (i) high levels of genetic diversity, (ii) no evidence of local inbreeding or reduced effective population size and (iii) multiple maternal mitochondrial lineages, a genetic footprint depicting independent colonization events. Furthermore, approximate Bayesian computations suggested several evolutionary trajectories from ancient isolation in glacial refugia during the Pleistocene, with biogeographical signatures of recent expansion probably confounded by human‐mediated mixing of different lineages. From an evolutionary and conservation perspective, this study highlights the ecological value of industrial areas, provided that ongoing regional gene flow is ensured within the existing lineage boundaries.     

Phylogeography of spruce beetles (Dendroctonus rufipennis Kirby) (Curculionidae: Scolytinae) in North America

Tree-feeding insects that are widespread in north temperate regions are excellent models for studying how past glaciations have impacted differentiation and speciation. We used mitochondrial DNA (mtDNA) sequences and allele frequencies at nine microsatellite loci to examine genetic population structure across the current range of the spruce beetle (Dendroctonus rufipennis), an economically important insect in North America. Two major haplotype groups occur across northern North America, from Newfoundland to Alaska, on white spruce (Picea glauca), and a third distinctive haplotype group occurs throughout the Rocky Mountains on Engelmann spruce (Picea engelmannii). The two mtDNA lineages found in northern populations are 3-4% divergent from each other and from the lineages found in the Rocky Mountains. Analyses of microsatellite data also suggest the existence of major population groupings associated with different geographical regions. In the Pacific Northwest, concordant contact zones for genetically distinct populations of spruce beetles and their principal hosts appear to reflect recent secondary contact. Although we could detect no evidence of historical mtDNA gene flow between allopatric population groups, patterns of variation in the Pacific Northwest suggest recent hybridization and introgression. Together with the pollen record for spruce, they also suggest that beetles have spread from at least three glacial refugia. A minimum estimate of divergence time between the Rocky Mountain and northern populations was 1.7 Myr (million years), presumably reflecting the combined effects of isolation during multiple glacial cycles.     

Mitochondrial phylogeography of a Beringian relict: the endemic freshwater genus of blackfish Dallia (Esociformes)

Mitochondrial genetic variability among populations of the blackfish genus Dallia (Esociformes) across Beringia was examined. Levels of divergence and patterns of geographic distribution of mitochondrial DNA lineages were characterized using phylogenetic inference, median‐joining haplotype networks, Bayesian skyline plots, mismatch analysis and spatial analysis of molecular variance (SAMOVA) to infer genealogical relationships and to assess patterns of phylogeography among extant mitochondrial lineages in populations of species of Dallia. The observed variation includes extensive standing mitochondrial genetic diversity and patterns of distinct spatial segregation corresponding to historical and contemporary barriers with minimal or no mixing of mitochondrial haplotypes between geographic areas. Mitochondrial diversity is highest in the common delta formed by the Yukon and Kuskokwim Rivers where they meet the Bering Sea. Other regions sampled in this study host comparatively low levels of mitochondrial diversity. The observed levels of mitochondrial diversity and the spatial distribution of that diversity are consistent with persistence of mitochondrial lineages in multiple refugia through the last glacial maximum.     

Identification of glacial refugia in south-eastern North America by phylogeographical analyses of a forest understorey plant, Trillium cuneatum

Aim We examine several hypotheses emerging from biogeographical and fossil records regarding glacial refugia of a southern thermophilic plant species. Specifically, we investigated the glacial history and post‐glacial colonization of a forest understorey species, Trillium cuneatum. We focused on the following questions: (1) Did T. cuneatum survive the Last Glacial Maximum (LGM) in multiple refugia, and (if so) where were they located, and is the modern genetic structure congruent with the fossil record‐based reconstruction of refugia for mesic deciduous forests? (2) What are the post‐glacial colonization patterns in the present geographical range? Location South‐eastern North America. Methods We sampled 45 populations of T. cuneatum throughout its current range. We conducted phylogeographical analyses based on maternally inherited chloroplast DNA (cpDNA haplotypes) and used TCS software to reconstruct intraspecific phylogeny. Results We detected six cpDNA haplotypes, geographically highly structured into non‐overlapping areas. With one exception, none of the populations had mixed haplotype composition. TCS analysis resulted in two intraspecific cpDNA lineages, with one clade subdivided further by shallower diversification. Main conclusions Our investigation revealed that T. cuneatum survived the LGM in multiple refugia, belonging to two (western, eastern) genealogical lineages geographically structured across south‐eastern North America. The western clade is confined to the south‐western corner of T. cuneatum’s modern range along the Lower Mississippi Valley, where fossil records document a major refugium of mesic deciduous forest. For the eastern clade, modern patterns of cpDNA haplotype distribution suggest cryptic vicariance, in the form of forest contractions and subsequent expansions associated with Pleistocene glacial cycles, rather than simple southern survival and subsequent northward colonization. The north–south partitioning of cpDNA haplotypes was unexpected, suggesting that populations of this rather southern thermophilic species may have survived in more northern locations than initially expected based on LGM climate reconstruction, and that the Appalachian Mountains functioned as a barrier to the dispersal of propagules originating in more southern refugia. Furthermore, our results reveal south‐west to north‐east directionality in historical migration through the Valley and Ridge region of north‐west Georgia.     

Genetic consequences of glacial survival and postglacial colonization in Norway spruce: combined analysis of mitochondrial DNA and fossil pollen

Norway spruce (Picea abies [L.] Karst.) is a broadly distributed European conifer tree whose history has been intensively studied by means of fossil records to infer the location of full‐glacial refugia and the main routes of postglacial colonization. Here we use recently compiled fossil pollen data as a template to examine how past demographic events have influenced the species’ modern genetic diversity. Variation was assessed in the mitochondrial nad1 gene containing two minisatellite regions. Among the 369 populations (4876 trees) assayed, 28 mitochondrial variants were identified. The patterns of population subdivision superimposed on interpolated fossil pollen distributions indicate that survival in separate refugia and postglacial colonization has led to significant structuring of genetic variation in the southern range of the species. The populations in the northern range, on the other hand, showed a shallow genetic structure consistent with the fossil pollen data, suggesting that the vast northern range was colonized from a single refugium. Although the genetic diversity decreased away from the putative refugia, there were large differences between different colonization routes. In the Alps, the diversity decreased over short distances, probably as a result of population bottlenecks caused by the presence of competing tree species. In northern Europe, the diversity was maintained across large areas, corroborating fossil pollen data in suggesting that colonization took place at high population densities. The genetic diversity increased north of the Carpathians, probably as a result of admixture of expanding populations from two separate refugia.     

Widespread secondary contact and new glacial refugia in the halophilic rotifer Brachionus plicatilis in the Iberian Peninsula

Small aquatic organisms harbour deep phylogeographic patterns and highly structured populations even at local scales. These patterns indicate restricted gene flow, despite these organisms' high dispersal abilities, and have been explained by a combination of (1) strong founder effects due to rapidly growing populations and very large population sizes, and (2) the development of diapausing egg banks and local adaptation, resulting in low effective gene flow, what is known as the Monopolization hypothesis. In this study, we build up on our understanding of the mitochondrial phylogeography of the halophilic rotifer Brachionus plicatilis in the Iberian Peninsula by both increasing the number of sampled ponds in areas where secondary contact is likely and doubling sample sizes. We analyzed partial mitochondrial sequences of 252 individuals. We found two deep mitochondrial DNA lineages differing in both their genetic diversity and the complexity of their phylogeographic structure. Our analyses suggest that several events of secondary contact between clades occurred after their expansion from glacial refugia. We found a pattern of isolation-by-distance, which we interpret as being the result of historical colonization events. We propose the existence of at least one glacial refugium in the SE of the Iberian Peninsula. Our findings challenge predictions of the Monopolization hypothesis, since coexistence (i.e., secondary contact) of divergent lineages in some ponds in the Iberian Peninsula is common. Our results indicate that phylogeographic structures in small organisms can be very complex and that gene flow between diverse lineages after population establishment can indeed occur.     

A new cytochrome b phylogroup of the common vole (Microtus arvalis) endemic to the Balkans and its implications for the evolutionary history of the species

The phylogeographic architecture of the common vole, Microtus arvalis, has been well‐studied using mitochondrial DNA and used to test hypotheses relating to glacial refugia. The distribution of the five described cytochrome b (cyt b) lineages in Europe west of Russia has been interpreted as a consequence of postglacial expansion from both southern and central European refugia. A recently proposed competing model suggests that the ‘cradle’ of the M. arvalis lineages is in western central Europe from where they dispersed in different directions after the Last Glacial Maximum. In the present study, we report a new cyt b lineage of the common vole from the Balkans that is not closely related to any other lineage and whose presence might help resolve these issues of glacial refugia. The Balkan phylogroup occurs along the southern distributional border of M. arvalis in central and eastern Bosnia and Herzegovina, Montenegro, and eastern Serbia. Further north and west in Slovenia, Bosnia and Herzegovina, and Serbia, common voles belong to the previously‐described Eastern lineage, whereas both lineages are sympatric in one site in Bosnia and Herzegovina. The Balkan phylogroup most reasonably occupied a glacial refugium already known for various Balkan endemic species, in contrast to the recently proposed model. South‐east Europe is an absolutely crucial area for understanding the postglacial colonization history of small mammals in Europe and the present study adds to the very few previous detailed phylogeographic studies of this region. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 788–796.     

Natural hybridization in eastern-Mediterranean firs: The case of Abies borisii-regis

The genus Abies is represented in southern Balkans by A. alba, A. cephalonica, and A. borisii-regis. To infer the status of southern-Balkans firs, as well as the extent and patterns of introgression within this taxonomical complex, we analyzed genetic variation patterns of 29 indigenous fir populations in Bulgaria, Macedonia, Greece, and Calabria using a combination of maternally and biparentally inherited markers. Three mitochondrial lineages were observed, one comprising Calabrian populations and two distributed in the Balkans, coinciding with A. alba and A. cephalonica. The boundary between lineages is sharp; only two populations containing a mixture of haplotypes were found. Bayesian analysis of population structure based on seven nuclear microsatellite (nSSR) loci revealed the existence of two clusters whose proportions exhibited a latitudinal cline with a width of 2.3° ( ≈ 255 km). Populations in the center of the latitudinal cline exhibit the most symmetrical, the flattest, and the broadest distributions of cluster proportions within individual tree genomes. A neighbor-net network reflects the cline resulting from the Bayesian analysis. The observed variation patterns are not consistent with the hypothesis of A. borisii-regis as a monophyletic taxon or a stabilized hybridogenous species resulted from ancient hybridization; the taxon rather is a product of recent introgression.     

Late Pleistocene divergence between eastern and western populations of wood ducks (Aix sponsa) inferred by the 'isolation with migration' coalescent method

During the Late Pleistocene, glaciers sundered many species into multiple glacial refugia where populations diverged in allopatry. Although deeply divergent mitochondrial DNA (mtDNA) lineages often reflect the number of refugia occupied, it is unlikely that populations that split during the recent Wisconsin glaciations will have reached reciprocal monophyly. We examined mtDNA control region sequences from eastern and western populations of wood ducks (Aix sponsa) to determine whether their current, disjunct distribution is consistent with the occupancy of two glacial refugia. We used the 'isolation with migration' coalescent method (im) to simultaneously estimate effective population sizes, maternal gene flow, and time since divergence. We found 24 unique haplotypes, none of which were shared between the eastern and western populations, but we did not find diagnostic monophyletic lineages suggestive of long-term isolation in multiple glacial refugia. However, a high Phi ST (0.31) indicates that eastern and western populations are well differentiated in mtDNA, and results from im suggest that these populations have been diverging, without extensive gene flow, for 10,000 to 124,000 years. Results from im further suggest that these populations most likely split about 34,000 years ago, and this time of divergence is consistent with the occupancy of multiple glacial refugia during the Late Wisconsin glaciation. Eastern wood ducks are characterized by high genetic diversity, a large effective population size, and a recent population expansion, while western wood ducks have much less genetic diversity, a smaller population size, and have not undergone a recent population expansion.     

Colonization routes of Pinus sylvestris inferred from distribution of mitochondrial DNA variation

Understanding the present-day distribution of molecular variation requires knowledge about the history of the species. Past colonization routes and locations of refugia of Scots pine (Pinus sylvestris) were inferred from variation in mitochondrial DNA in material collected from 37 populations located in countries within, and immediately adjacent to the continent of Europe. Two mitochondrial regions, nad1 intron (exon B/C) and nad7 intron 1, were included in the study. Differentiation in maternally inherited mitochondria was high (G _ST′ = 0.824). Two new haplotypes were found at the nad7 intron 1. The occurrence of a 5-bp indel variant was restricted to the Turkish Kalabak population and a 32 bp only found in Central, Eastern, and Northern Europe. The complete absence of the 32-bp indel from the Mediterranean peninsulas supports the view that coniferous forests existed outside these areas during the last glacial maximum, and these populations contributed to the subsequent colonization of the northern parts of Europe. P. sylvestris shares features of its glacial and postglacial history with two other northern, cold-tolerant tree species, Picea abies and Betula sp. These three species differ from many other European trees for which pollen core and molecular evidence indicate colonization from southern refugia after the last glacial period.     

Landscape genetics reveal broad and fine‐scale population structure due to landscape features and climate history in the northern leopard frog (Rana pipiens) in North Dakota

Prehistoric climate and landscape features play large roles structuring wildlife populations. The amphibians of the northern Great Plains of North America present an opportunity to investigate how these factors affect colonization, migration, and current population genetic structure. This study used 11 microsatellite loci to genotype 1,230 northern leopard frogs (Rana pipiens) from 41 wetlands (30 samples/wetland) across North Dakota. Genetic structure of the sampled frogs was evaluated using Bayesian and multivariate clustering methods. All analyses produced concordant results, identifying a major east–west split between two R. pipiens population clusters separated by the Missouri River. Substructuring within the two major identified population clusters was also found. Spatial principal component analysis (sPCA) and variance partitioning analysis identified distance, river basins, and the Missouri River as the most important landscape factors differentiating R. pipiens populations across the state. Bayesian reconstruction of coalescence times suggested the major east–west split occurred ~13–18 kya during a period of glacial retreat in the northern Great Plains and substructuring largely occurred ~5–11 kya during a period of extreme drought cycles. A range‐wide species distribution model (SDM) for R. pipiens was developed and applied to prehistoric climate conditions during the Last Glacial Maximum (21 kya) and the mid‐Holocene (6 kya) from the CCSM4 climate model to identify potential refugia. The SDM indicated potential refugia existed in South Dakota or further south in Nebraska. The ancestral populations of R. pipiens in North Dakota may have inhabited these refugia, but more sampling outside the state is needed to reconstruct the route of colonization. Using microsatellite genotype data, this study determined that colonization from glacial refugia, drought dynamics in the northern Great Plains, and major rivers acting as barriers to gene flow were the defining forces shaping the regional population structure of R. pipiens in North Dakota.     

Rethinking long‐term vegetation dynamics: multiple glacial refugia and local expansion of a species complex

Evidence is accumulating that some arcto‐boreal plant taxa persisted through the last glacial maximum (LGM) in Alaska and adjacent Canada. However, the spatial patterns of glacial persistence and associated postglacial colonization remain largely unknown. In this study, we investigated the LGM refugia of an alder (Alnus) species complex (n = 3 taxa) and assess the spatiotemporal dynamics of Alnus in this vast region. Specifically, we conducted high‐throughput DNA sequencing (ddRADseq) on Alnus foliar samples collected from a dense population network to investigate patterns of genetic structure and infer the presence of glacial lineages. Species distribution modeling (SDM) was used to investigate the probability and possible locations of glacial persistence. These analyses were integrated and then compared with fossil pollen data to identify the locations of refugial populations and spatial patterns of postglacial colonization. Our genetic analyses revealed two glacial lineages with separate geographic origins for each Alnus taxon, suggesting that the genus persisted in multiple LGM refugia. Non‐overlapping hindcast distributions based on SDMs further support the presence of multiple, spatially distinct refugia. These ddRADseq and SDM results, in conjunction with reassessment of fossil pollen records, suggest that Alnus expanded from several population nuclei that existed during the LGM and coalesced during the Holocene to form its present range. These results challenge the unidirectional model for postglacial vegetation expansion, implying that climate buffering associated with landscape heterogeneity and adaptation to millennial‐scale environmental variability played important roles in driving late‐Quaternary population dynamics.     

Genetic diversity of North American populations of Cristatella mucedo, inferred from microsatellite and mitochondrial DNA

Research over the past 20 years has shown, with the help of molecular markers, that the population genetics and distribution patterns of freshwater invertebrates in North America are often more complex than was previously believed. Here we extend this research to an, as yet, unstudied but widespread and common group, the freshwater bryozoans. Colonies of the bryozoan Cristatella mucedo were collected from a number of lakes across central North America, and were characterized genetically by analysis of microsatellite loci and mitochondrial DNA (mtDNA) cytochrome b sequences. The microsatellites illustrate a pattern of generally diverse and highly differentiated populations that contain little evidence of recent gene flow. The mtDNA sequences yielded highly variable levels of divergence, ranging from 0.0 to 8.8% within populations, and 0.0 to 9.8% among populations. The multiple divergent mtDNA lineages within populations provide evidence for repeated colonization events. The lack of clustering of haplotypes by site suggests that there has been widespread dispersal of multiple genetic lineages since the last ice age. While some of the haplotype lineages may have evolved in disjunct glacial refugia, the maximum levels of divergence predate the time since the last glacial-interglacial cycles. It is likely that multiple factors including vicariance events, patterns of dispersal, localized extinction, and an unusual life history, explain the unique phylogeographic patterns evident today in populations of C. mucedo.     

Nuclear and mitochondrial phylogeography of the European fire-bellied toads Bombina bombina and Bombina variegata supports their independent histories

Exact location and number of glacial refugia still remain unclear for many European cold‐blooded terrestrial vertebrates. We performed a fine‐scaled multilocus phylogeographic analysis of two Bombina species combining mitochondrial variation of 950 toads from 385 sites and nuclear genes (Rag‐1, Ncx‐1) from a subset of samples to reconstruct their colonization and contemporary variation patterns. We identified the lowlands northwest of the Black Sea and the Carpathians to be important refugial areas for B. bombina and B. variegata, respectively. This result emphasizes the importance of Central European refugia for ectothermic terrestrial species, far north of the Mediterranean areas regarded as exclusive glacial refugia for the animals. Additional refugia for B. variegata have been located in the southern Apennines and Balkans. In contrast, no evidence for the importance of other east European plains as refugial regions has been found. The distribution of mtDNA and Ncx‐1 variation suggests the presence of local refugia near the Black Sea for B. bombina; however, coalescent simulations did not allow to distinguish whether one or two refugia were present in the region. Strong genetic drift apparently accompanied postglacial expansions reducing diversity in the colonization areas. Extended sampling, coupled with the multilocus isolation with migration analysis, revealed a limited and geographically restricted gene flow from the Balkan to Carpathian populations of B. variegata. However, despite proximity of inferred B. bombina and B. variegata refugia, gene exchange between them was not detected.     

Phylogeography and Pleistocene refugia of the adder (Vipera berus) as inferred from mitochondrial DNA sequence data

In order to contribute to the debate about southern glacial refugia used by temperate species and more northern refugia used by boreal or cold-temperate species, we examined the phylogeography of a widespread snake species (Vipera berus) inhabiting Europe up to the Arctic Circle. The analysis of the mitochondrial DNA (mtDNA) sequence variation in 1043 bp of the cytochrome b gene and in 918 bp of the noncoding control region was performed with phylogenetic approaches. Our results suggest that both the duplicated control region and cytochrome b evolve at a similar rate in this species. Phylogenetic analysis showed that V. berus is divided into three major mitochondrial lineages, probably resulting from an Italian, a Balkan and a Northern (from France to Russia) refugial area in Eastern Europe, near the Carpathian Mountains. In addition, the Northern clade presents an important substructure, suggesting two sequential colonization events in Europe. First, the continent was colonized from the three main refugial areas mentioned above during the Lower-Mid Pleistocene. Second, recolonization of most of Europe most likely originated from several refugia located outside of the Mediterranean peninsulas (Carpathian region, east of the Carpathians, France and possibly Hungary) during the Mid-Late Pleistocene, while populations within the Italian and Balkan Peninsulas fluctuated only slightly in distribution range, with larger lowland populations during glacial times and with refugial mountain populations during interglacials, as in the present time. The phylogeographical structure revealed in our study suggests complex recolonization dynamics of the European continent by V. berus, characterized by latitudinal as well as altitudinal range shifts, driven by both climatic changes and competition with related species.     

Colonization of Northern Europe by Zygaena filipendulae (Lepidoptera)

Northern and mountainous ice sheets have expanded and contracted many times due to ice ages. Consequently, temperate species have been confined to refugia during the glacial periods wherefrom they have recolonized warming northern habitats between ice ages. In this study, we compare the gene CYP405A2 between different populations of the common burnet moth Zygaena filipendulae from across the Western Palearctic region to illuminate the colonization history of this species. These data show two major clusters of Z. filipendulae populations possibly reflecting two different refugial populations during the last ice age. The two types of Z. filipendulae only co‐occur in Denmark, Sweden, and Scotland indicating that Northern Europe comprise the hybridization zone where individuals from two different refugia met after the last ice age. Bayesian phylogeographic and ecological clustering analyses show that one cluster probably derives from an Alpe Maritime refugium in Southern France with ancestral expansive tendencies to the British Isles in the west, touching Northern Europe up to Denmark and Sweden, and extending throughout Central Europe into the Balkans, the Peleponnes, and South East Europe. The second cluster encompasses East Anatolia as the source area, from where multiple independent dispersal events to Armenia, to the Alborz mountains in north‐western Iran, and to the Zagros mountains in western Iran are suggested. Consequently, the classical theory of refugia for European temperate species in the Iberian, Italian, and Balkan peninsulas does not fit with the data from Z. filipendulae populations, which instead support more Northerly, mountainous refugia.     

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

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