A new scenario for the quaternary history of European beech populations: palaeobotanical evidence and genetic consequences

Here, palaeobotanical and genetic data for common beech (Fagus sylvatica) in Europe are used to evaluate the genetic consequences of long-term survival in refuge areas and postglacial spread. Four large datasets are presented, including over 400 fossil-pollen sites, 80 plant-macrofossil sites, and 450 and 600 modern beech populations for chloroplast and nuclear markers, respectively. The largely complementary palaeobotanical and genetic data indicate that: (i) beech survived the last glacial period in multiple refuge areas; (ii) the central European refugia were separated from the Mediterranean refugia; (iii) the Mediterranean refuges did not contribute to the colonization of central and northern Europe; (iv) some populations expanded considerably during the postglacial period, while others experienced only a limited expansion; (v) the mountain chains were not geographical barriers for beech but rather facilitated its diffusion; and (vi) the modern genetic diversity was shaped over multiple glacial-interglacial cycles. This scenario differs from many recent treatments of tree phylogeography in Europe that largely focus on the last ice age and the postglacial period to interpret genetic structure and argue that the southern peninsulas (Iberian, Italian and Balkan) were the main source areas for trees in central and northern Europe.     

High chloroplast haplotype diversity in Greek populations of beech (Fagus sylvatica L.)

The distribution of chloroplast DNA (cpDNA) variations in Greek beech ( Fagus sylvatica L.) populations was studied using chloroplast microsatellite markers. Thirteen haplotypes were identified from 40 populations by combining three different primers. Most of the cpDNA variation was distributed among populations, but a considerable variation was also observed within populations. The total diversity was very high for all regions. The N_st/G_st comparison was significant, indicating phylogenetic subdivision, but no strong spatial structure was detected, suggesting complex post-glacial migration patterns. Possible scenarios explaining this diversity pattern include the existence of several separated refugia in the region, the recolonisation of mountains by different beech lineages and the formation of an introgression zone between two different beech subspecies in the eastern part of the country.     

Post-glacial migration of silver fir (Abies alba Mill.) in the south-western Alps

Aim  Previous studies have failed to reconstruct the regional post-glacial migration pattern of Abies alba in southern France. Based on the first exhaustive compilation of palaeoecological data in this region, we present the state-of-the-art and attempt to synthesize the available information concerning glacial refugia and post-glacial migration, and analyse the information with regard to climate and orography.
Location  South-western Alps and adjacent areas, southern France.
Methods  The work compiles the available palaeoecological data in the south-western Alps (52 sites, 290 radiocarbon dates). The post-glacial migration pattern of Abies alba is reconstructed based on 22 selected palynological analyses (11 well-dated reference sites and 11 supplementary ones).
Results  The geographical patterns of approaching area limit, immigration and expansion are reconstructed at the scale of the southern French Alps.
Main conclusions  Despite previous assertions, the evidence of refugia in southern France is non-existent. The late-glacial records of fir pollen, previously interpreted in French Mediterranean regions and on adjacent foot-hills as possibly reflecting regional refugia, most probably correspond to reworking phenomena or long-distance pollen transport. Fir migration, originating in the Apennine refugia and through the south-western extremity of the Alps, was extremely rapid in the southern French Alps, only spanning a few centuries between 10,100 and 9800 cal. yr bp . The subsequent spread of fir populations was controlled by local parameters, such as the aridity of the inner valleys, which resulted in a delayed expansion in comparison to other regions. Abies almost disappeared from the south-western Alps during the Roman era, around 2000 cal. yr bp .     

Complex fine‐scale phylogeographical patterns in a putative refugial region for Fagus sylvatica (Fagaceae)

Broad‐scale plastid (chloroplast) DNA studies of beech (Fagus sylvatica) populations suggest the existence of glacial refugia and introgression zones in south‐eastern Europe. We choose a possible refugium of beech in northern Greece, Mt. Paggeo, which hosts a private plastid haplotype for beech, to conduct a fine‐scale genetic study. We attempt to confirm or reject the hypothesis of the existence of a small‐scale refugium and to gain an understanding of the ecological and topographical factors affecting the spatial distribution of plastid haplotypes in the area. Our results reveal a high haplotype diversity on Mt. Paggeo, but the overall distribution of haplotypes shows no significant correlation with the ecological characteristics of the beech forests. However, the private haplotype is found at high frequencies in beech forests located in or near ravines, having a high spatial overlap with a relict vegetation type occurring in ecological conditions found mainly in ravines. This result emphasizes the importance of topography in the existence of glacial refugia in the wider area. Furthermore, haplotypes originating from two more widespread beech lineages in Greece are found on Mt. Paggeo, indicating a possible mixing of populations originating from a local refugium with populations from remote refugia that possibly migrated into the area after the last glaciation. © 2014 The Linnean Society of London, Botanical Journal of the Linnean Society, 2014, 174 , 516–528.     

Phylogeny of European Pyrochroa (Coleoptera,Pyrochroidae) reveals cryptic taxa and different glacial histories

Only three saproxylic species of Pyrochroinae (Coleoptera: Pyrochroidae) are distributed in Europe, two of which belonging to Pyrochroa: P. coccinea and P. serraticornis. However, P. serraticornis is polytypic, for the presence of the endemic subspecies P. s. kiesenwetteri in southern Italy. Using both molecular and morphological data, we explored the phylogeny of the European Pyrochroa species. A multilocus (COI, CAD, 28S) phylogenetic analysis helped highlight different evolutionary histories for the two examined species. First, P. coccinea, distributed throughout Europe, showed a high differentiation among Italian and European populations. Furthermore, three different taxonomic entities were identified within P. serraticornis, among which the cryptic species Pyrochroa bifoveata sp. n. from central Europe is described and illustrated. A comprehensive identification key to the European Pyrochroinae is also provided. Our results also suggested an historical survival of P. coccinea and P. s. kiesenwetteri in glacial refugia in Italy, and a subsequent post-glacial spread of the former species throughout the Peninsula. In contrast, the current distribution of P. s. serraticornis likely originated from a post-glacial colonization of western European relict populations, while the survival of P. bifoveata plausibly occurred in more eastern glacial refugia (e.g. Carpathian or Balkan regions). Similarly, the European populations of P. coccinea could have originated from relict populations in glacial refugia out from the Italian Peninsula. More comprehensive data on the taxonomy, ecology and biogeography of Pyrochroa are needed to learn more about these species and to help preserve the European saproxylic fauna.     

Inferring long-distance dispersal and topographic barriers during post-glacial colonization from the genetic structure of red maple (Acer rubrum L.) in New England

Aim  This study aims to assess the role of long-distance seed dispersal and topographic barriers in the post-glacial colonization of red maple ( Acer rubrum L.) using chloroplast DNA (cpDNA) variation, and to understand whether this explains the relatively higher northern diversity found in eastern North American tree species compared with that in Europe.
Location  North-eastern United States.
Methods  The distribution of intraspecific cpDNA variation in temperate tree populations has been used to identify aspects of post-glacial population spread, including topographic barriers to population expansion and spread by long-distance seed dispersal. We sequenced c.  370 cpDNA base pairs from 221 individuals in 100 populations throughout the north-eastern United States, and analysed spatial patterns of diversity and differentiation.
Results  Red maple has high genetic diversity near its northern range limit, but this diversity is not partitioned by topographic barriers, suggesting that the northern Appalachian Mountains were not a barrier to the colonization of red maple. We also found no evidence of the patchy genetic structure that has been associated with spread by rare long-distance seed dispersal in previous studies.
Main conclusions  Constraints on post-glacial colonization in eastern North America seem to have been less stringent than those in northern Europe, where bottlenecks arising from long-distance colonization and topographic barriers appear to have strongly reduced genetic diversity. In eastern North America, high northern genetic diversity may have been maintained by a combination of frequent long-distance dispersal, minor topographic obstacles and diffuse northern refugia near the ice sheet.     

Multiple differentiation centres of a non-Mediterranean butterfly species in south-eastern Europe

Aim  The analysis of the phylogeographical structures of many European species reveals the importance of Mediterranean glacial refugia for many thermophilic species, but also underlines the relevance of extra-Mediterranean glacial differentiation centres for a number of temperate species. In this context, phylogeographical analyses of species from south-eastern Europe are highly important for a comprehensive understanding of Europe as a whole.
Location  Romania and Bulgaria.
Methods  We analysed 19 allozyme loci for 615 individuals of the temperate butterfly species Erebia medusa from 28 populations.
Results  These populations had an intermediate genetic diversity, but the Bulgarian populations were significantly more diverse than the ones north of the Danube in Romania. The differentiation among populations was strong, and 52.1% of the genetic variance among populations was distributed between these two countries. The genetic differentiation was considerably stronger in Romania than in Bulgaria, but several sublineages were distinguished within each of these countries.
Main conclusions  The observed genetic structure is so strong that it is most probably the result of glacial differentiation processes in south-eastern Europe and not a post-glacial structure. The strong differentiation into the two groups north and south of the Danube suggests a separating effect by this river valley. The strong differentiation accompanied with genetic impoverishment in Romania suggests the existence of several differentiation centres: at least two small ones on the southern slopes of the southern Carpathians and one in the eastern Carpathian Basin. The considerably weaker differentiation among the Bulgarian samples and their significantly higher genetic diversity imply that gene flow occurred among different regions of Bulgaria during the last ice age.     

Glacial refugia: sanctuaries for allelic richness, but not for gene diversity

Glacial refugia are generally expected to harbor higher levels of genetic diversity than are areas that have been colonized after the retreat of the glaciers because colonization often involves only a few individuals. A new paper by Comps et al. challenges this expectation by demonstrating a more complex situation in the European beech Fagus sylvatica, for which some measures of genetic diversity are higher in newly colonized areas than in refugia. The key to understanding this counter-intuitive result rests both in the estimators used to measure genetic diversity and in the processes affecting these estimators during postglacial recolonization.     

Range-wide phylogeography of the European temperate-montane herbaceous plant Meum athamanticum Jacq.: evidence for periglacial persistence

Aim  The aim of this study is to analyse the genetic population structure of Meum athamanticum Jacq. in order to explore the alternative hypotheses (1) that the central and northern highland populations are the result of post-glacial recolonization from southern refugia, and the disjunct distribution of M. athamanticum can be explained by modern ecological conditions, or (2) that extant populations north of the Alps and Pyrenees persisted in situ during glacial periods.
Location  Europe.
Methods  Variation of amplified fragment length polymorphisms (AFLPs) was analysed for 23 populations from the entire range of the species. We used band-based approaches and methods based on allele frequencies to measure genetic diversity and to identify population structure.
Results  Our analyses reveal a north–south differentiation within M. athamanticum . High levels of genetic diversity, as well as private fragments, are found in populations both north and south of the Alps. Differentiation among populations is lower in the northern than in the southern population group, and significant isolation-by-distance is found only in the latter group.
Main conclusions  Our results indicate that M. athamanticum survived the last ice age in multiple refugia throughout its contemporary range and did not expand into areas north of the Alps from southern refugia. We found evidence that regional-scale migration in northern, formerly periglacial, parts of the species range has resulted in the intermingling of populations. In contrast, southern populations are characterized by long-term isolation. The south-west Alps represent an area where immigration and mixing of populations from northern and southern refugia appears to have taken place.     

Some considerations about the postglacial history and recent distribution ofFagus sylvatica in the NW Iberian Peninsula

This paper presents a re-assessment of the evidence of beech (Fagus sylvatica) presence on the NW Iberian Peninsula in the different stages of the Late Quaternary. For this purpose, numerous pollen diagrams were re-examined and the present-day distribution of the species was documented. Beech pollen has been recorded in different periods, before the Würm Pleni-Glacial, during the Late Glacial and the Holocene, leading to the hypothesis that there existed local refugia forFagus sylvatica in the NW Iberian Peninsula, comparable to those in more eastern areas of southern Europe, the existence of which has been recognized.     

The Holocene history of low altitude Mediterranean Fagus sylvatica forests in southern France

Questions

As the dominant tree in many European forests, Fagus sylvatica functions as an ecosystem engineer, yet its istory remains little understood. Here we ask: (a) are there indications for its presence in southeast France during the last Glacial period; (b) what was the timing of the expansion and decline of F. sylvatica dominated forests; (c) which factors influenced their dynamics and in particular to what extent did past precipitation changes impact upon them; and (d) at which altitudes did these beech forests occur within the region?

Location

Languedoc, the French Mediterranean area.

Method

This article presents a well dated and high‐resolution pollen sequence covering the last 7,800 years from the Palavas Lagoon in the Languedoc together with a review of Fagus charcoal occurrences in the Languedoc and the lower Rhône Valley, and a review of pollen data from a compilation of 69 sites in southeast France.

Results

The Palavas pollen sequence provides a regional summary of F. sylvatica abundance changes near the Mediterranean coast. Around 6,000 years cal BP , an abrupt transition from small beech populations to well‐developed forests is recorded. The maximum development of beech forests occurred between 4,000 and 3,000 years cal BP , while F. sylvatica started to regress after 3,000 years cal BP .

Conclusion

Scattered F. sylvatica populations probably survived throughout southern France during the last Glacial period. F. sylvatica started to spread around 8,000 years cal BP while beech forests never expanded before 6,000 years cal BP . The complex patterns of F. sylvatica expansion in southern France after 6,000 years cal BP suggests that a combination of global (climate change) and local (human impact) factors were responsible for this major change. Recurrent abrupt climate changes, the aridity trend and human deforestation caused beech forests to decline after 3,000 years cal BP .

     

Phylogeographic pattern of the European forest grass species Hordelymus europaeus: cpDNA evidence

The phylogeographical structure of the European forest grass Hordelymus europaeus (Poaceae) was studied by sequencing three non-coding regions of chloroplast DNA. Forty-three populations within the entire natural distribution area were analysed. The greatest haplotype variation and divergence were revealed on Balkan and Apennine Peninsula, suggesting main glacial refugia in these regions. Among southern refugia, probably only the Balkans could have remarkably contributed to postglacial re-colonisation of the northern parts of Europe. A distinctly different haplotype group found on the Crimean Peninsula and in central Italy may represent either relicts of a previously more widespread ancestor or result of long-distance dispersal. The phylogeographical pattern found in H. europaeus is to certain extent similar with that found in Fagus sylvatica. This might imply a partly common postglacial colonisation history of these ecologically narrowly tied species. Nevertheless, unlike in the case of F. sylvatica, we did not found convincing evidence for the existence of Central European glacial refugia for H. europaeus.     

The spread of<Emphasis Type="Italic"> Abies</Emphasis> throughout Europe since the last glacial period: combined macrofossil and pollen data

In this paper we reconstruct the location of the last glacial refugia and postglacial spread of Abies throughout Europe based on combined pollen and macrofossil data. More than 208 pollen sequences available in the European Pollen Database (EPD) and 38 macrofossil sites are used to produce distribution maps encompassing a time span between 38000 and 5500 years B.P. The investigation excludes more recent periods, because these could be strongly influenced by human impact. The pollen data presented here confirm long-lasting refugial areas such as southern Italy and Greece already described in previous studies. The combined pollen/macrofossil dataset identifies the Pyrenees as a further important refugium. In addition the pollen data indicate potential refugia in south-east France and north-west Italy. Possible migration tracks of Abies are discussed by comparing the palaeobotanical evidence with isozyme studies on gene markers of recent fir populations.Electronic Supplementary Material Supplementary material is available in the online version of this article at http://dx.doi.org/10.1007/s00334-004-0049-4.     

Imprints of glacial refugia in the modern genetic diversity of Pinus sylvestris

Aim To understand the impact of glacial refugia and migration pathways on the modern genetic diversity of Pinus sylvestris. Location The study was carried out throughout Europe. Methods An extended set of data of pollen and macrofossil remains was used to locate the glacial refugia and reconstruct the migrating routes of P. sylvestris throughout Europe. A vegetation model was used to simulate the extent of the potential refugia during the last glacial period. At the same time a genetic survey was carried out on this species. Results The simulated distribution of P. sylvestris during the last glacial period is coherent with the observed fossil data, which showed a patchy distribution of the refugia between c. 40° N and 50° N. Several migrational fronts were detected within the Iberian and the Italian peninsulas, and outside the Hungarian plain and around the Alps. The modern mitochondrial DNA depicted three different haplotypes for P. sylvestris. Two distinct haplotypes were restricted to northern Spain and Italy, and the third haplotype dominated most of the present‐day remaining distribution range of P. sylvestris in Europe. Main conclusions During the last glacial period P. sylvestris was constrained under severe climatic conditions to survive in scattered and restricted refugial areas. Combining palaeoenvironmental data, vegetation modelling and the genetic data, we have shown that the long‐term isolation in the glacial refugia and the migrational process during the Holocene have played a major role in shaping the modern genetic diversity of P. sylvestris in Europe.     

Multiple Glacial Refugia of the Low-Dispersal Ground Beetle Carabus irregularis: Molecular Data Support Predictions of Species Distribution Models

Classical glacial refugia such as the southern European peninsulas were important for species survival during glacial periods and acted as sources of post-glacial colonisation processes. Only recently, some studies have provided evidence for glacial refugia north of the southern European peninsulas. In the present study, we combined species distribution models (SDMs) with phylogeographic analyses (using mitochondrial DNA = mtDNA) to investigate if the cold-adapted, stenotopic and flightless ground beetle species, Carabus irregularis, survived the Last Glacial Maximum (LGM) in classical and/or other refugia. SDMs (for both a western European and for a Carpathian subgroup) were calculated with MAXENT on the basis of 645 species records to predict current and past distribution patterns. Two mtDNA loci (CO1 and ND5, concatenated sequence length: 1785 bp) were analyzed from 91 C. irregularis specimens to reconstruct the phylogeography of Central and eastern European populations and to estimate divergence times of the given lineages. Strong intra-specific genetic differentiation (inter-clade Φ_ST values ranged from 0.92 to 0.99) implied long-term isolation of major clades and subsclades. The high divergence between the nominate subspecies and the Carpathian subspecies C. i. montandoni points to two independent species rather than subspecies (K-2P distance 0.042 ± 0.004; supposed divergence of the maternal lineages dated back 1.6 to 2.5 million years BP) differing not only morphologically but also genetically and ecologically from each other. The SDMs also inferred classical as well as other refugia for C. irregularis, especially north of the Alps, in southeastern Europe and in the Carpathians. The coincidences between the results of both methods confirm the assumption of multiple glacial refugia for the studied species and the usefulness of combining methodological approaches for the understanding of the history of low-dispersal insect species.     

Beyond the Mediterranean peninsulas: evidence of central European glacial refugia for a temperate forest mammal species, the bank vole (Clethrionomys glareolus)

This study details the phylogeographic pattern of the bank vole, Clethrionomys glareolus, a European rodent species strongly associated with forest habitat. We used sequences of 1011 base pairs of the mitochondrial DNA cytochrome b gene from 207 bank voles collected in 62 localities spread throughout its distribution area. Our results reveal the presence of three Mediterranean (Spanish, Italian and Balkan) and three continental (western, eastern and 'Ural') phylogroups. The endemic Mediterranean phylogroups did not contribute to the post-glacial recolonization of much of the Palaearctic range of species. Instead, the major part of this region was apparently recolonized by bank voles that survived in glacial refugia in central Europe. Moreover, our phylogeographic analyses also reveal differentiated populations of bank voles in the Ural mountains and elsewhere, which carry the mitochondrial DNA of another related vole species, the ruddy vole (Clethrionomys rutilus). In conclusion, this study demonstrates a complex phylogeographic history for a forest species in Europe which is sufficiently adaptable that, facing climate change, survives in relict southern and northern habitats. The high level of genetic diversity characterizing vole populations from parts of central Europe also highlights the importance of such regions as a source of intraspecific genetic biodiversity.     

Postglacial colonization of Europe by the barbastelle bat: agreement between molecular data and past predictive modelling

The barbastelle (Barbastella barbastellus) is a rare forest bat with a wide distribution in Europe. Here, we combine results from the analysis of two mtDNA fragments with species distribution modelling to determine glacial refugia and postglacial colonization routes. We also investigated whether niche conservatism occurs in this species. Glacial refugia were identified in the three southern European peninsulas: Iberia, Italy and the Balkans. These latter two refugia played a major role in the postglacial colonization process, with their populations expanding to England and central Europe, respectively. Palaeo‐distribution models predicted that suitable climatic conditions existed in the inferred refugia during the last glacial maximum (LGM). Nevertheless, the overlap between the current and the LGM distributions was almost inexistent in Italy and in the Balkans, meaning that B. barbastellus populations were forced to shift range between glacial and interglacial periods, a process that probably caused some local extinctions. In contrast, Iberian populations showed a ‘refugia within refugium’ pattern, with two unconnected areas containing stable populations (populations that subsisted during both glacial and interglacial phases). Moreover, the match between LGM models and the refugial areas determined by molecular analysis supported the hypothesis of niche conservatism in B. barbastellus. We argue that geographic patterns of genetic structuring, altogether with the modelling results, indicate the existence of four management units for conservation: Morocco, Iberia, Italy and UK, and Balkans and central Europe. In addition, all countries sampled possessed unique gene pools, thus stressing the need for the conservation of local populations.     

Long-range transport of beech (<Emphasis Type="Italic">Fagus sylvatica</Emphasis> L.) pollen to Catalonia (north-eastern Spain)

Local and long-range transport of beech (Fagus sylvatica) pollen was analysed by using 23-year data (1983-2007) at six stations in Catalonia, Spain, and numerical simulations. Back trajectories and synoptic meteorology indicated a consistent north European provenance during beech pollen peak days. Specifically, the area from northern Italy to central Germany was the most probable source, as indicated by a source-receptor model based on back trajectories. For the event with the highest pollen levels (17 May 2004), back trajectories indicated a source in the Vosges (NE France) and the Schwarzwald (SW Germany) regions. By applying a mesoscale model (MM5) to this event, pollen transport could be further refined, allowing its entrance to Catalonia through the lower easternmost pass of the Pyrenees (the Alberes pass, 500 m a.s.l.) to be described. Hourly counts of Fagus pollen allowed the timing of pollen arrival during this episode to be matched with the model results regarding the above-mentioned passage. This study may help to interpret some results of modern beech genetic diversity and contribute to the understanding of paleopalynological records by taking long-range transport into consideration.     

A complex scenario of glacial survival in Mediterranean and continental refugia of a temperate continental vole species (Microtus arvalis) in Europe

The role of glacial refugia in shaping contemporary species distribution is a long-standing question in phylogeography and evolutionary ecology. Recent studies are questioning previous paradigms on glacial refugia and postglacial recolonization pathways in Europe, and more flexible phylogeographic scenarios have been proposed. We used the widespread common vole Microtus arvalis as a model to investigate the origin, locations of glacial refugia, and dispersal pathways, in the group of “Continental” species in Europe. We used a Bayesian spatiotemporal diffusion analysis (relaxed random walk model) of cytochrome b sequences across the species range, including newly collected individuals from 10 Iberian localities and published sequences from 68 localities across 22 European countries. Our data suggest that the species originated in Central Europe, and we revealed the location of multiple refugia (in both southern peninsulas and continental regions) for this continental model species. Our results confirm the monophyly of Iberian voles and the pre-LGM divergence between Iberian and European voles. We found evidence of restricted postglacial dispersal from refugia in Mediterranean peninsulas. We inferred a complex evolutionary and demographic history of M. arvalis in Europe over the last 50,000 years that does not adequately fit previous glacial refugial scenarios. The phylogeography of M. arvalis provides a paradigm of ice-age survival of a temperate continental species in western and eastern Mediterranean peninsulas (sources of endemism) and multiple continental regions (sources of postglacial spread). Our findings also provide support for a major role of large European river systems in shaping geographic boundaries of M. arvalis in Europe.     

Multiplexed microsatellite markers for genetic studies of beech

European beech (Fagus sylvatica L.) is one of the economically most important broadleaved tree species in Europe and has become a model for studying climate change effects on forests. Multiplex PCR of microsatellites is a fast and cost-effective technique allowing high-throughput genotyping. Here we present the procedure used to develop two multiplex kits (8-plexes) for European beech. We paid particular attention to quality control throughout all steps of the multiplex kits development (null allele detection, error rate measurements, linkage disequilibrium). Preliminary assays suggest that the 16 amplified loci are largely devoid of null alleles and allow rapid and cost-effective genotyping of beech with low error rates. The two kits, which differ in their levels of polymorphism, most likely due to marker origin, were also informative in seven other beech species tested.