Northern glacial refugia and altitudinal niche divergence shape genome‐wide differentiation in the emerging plant model Arabidopsis arenosa

Quaternary climatic oscillations profoundly impacted temperate biodiversity. For many diverse yet undersampled areas, however, the consequences of this impact are still poorly known. In Europe, particular uncertainty surrounds the role of Balkans, a major hotspot of European diversity, in postglacial recolonization of more northerly areas, and the Carpathians, a debatable candidate for a northern ‘cryptic’ glacial refugium. Using genome‐wide SNPs and microsatellites, we examined how the interplay of historical processes and niche shifts structured genetic diversity of diploid Arabidopsis arenosa, a little‐known member of the plant model genus that occupies a wide niche range from sea level to alpine peaks across eastern temperate Europe. While the northern Balkans hosted one isolated endemic lineage, most of the genetic diversity was concentrated further north in the Pannonian Basin and the Carpathians, where it likely survived the last glaciation in northern refugia. Finally, a distinct postglacial environment in northern Europe was colonized by populations of admixed origin from the two Carpathian lineages. Niche differentiation along altitude‐related bioclimatic gradients was the main trend in the phylogeny of A. arenosa. The most prominent niche shifts, however, characterized genetically only slightly divergent populations that expanded into narrowly defined alpine and northern coastal postglacial environments. Our study highlights the role of eastern central European mountains not only as refugia for unique temperate diversity but also sources for postglacial expansion into novel high‐altitude and high‐latitude niches. Knowledge of distinct genetic substructure of diploid A. arenosa also opens new opportunities for follow‐up studies of this emerging model of evolutionary biology.     

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.     

Population genomic evidence for plant glacial survival in Scandinavia

Quaternary glaciations have played a major role in shaping the genetic diversity and distribution of plant species. Strong palaeoecological and genetic evidence supports a postglacial recolonization of most plant species to northern Europe from southern, eastern and even western glacial refugia. Although highly controversial, the existence of small in situ glacial refugia in northern Europe has recently gained molecular support. We used genomic analyses to examine the phylogeography of a species that is critical in this debate. Carex scirpoidea Michx subsp. scirpoidea is a dioecious, amphi‐Atlantic arctic–alpine sedge that is widely distributed in North America, but absent from most of Eurasia, apart from three extremely disjunct populations in Norway, all well within the limits of the Weichselian ice sheet. Range‐wide population sampling and variation at 5,307 single nucleotide polymorphisms show that the three Norwegian populations comprise unique evolutionary lineages divergent from Greenland with high between‐population divergence. The Norwegian populations have low within‐population genetic diversity consistent with having experienced genetic bottlenecks in glacial refugia, and host private alleles that probably accumulated in long‐term isolated populations. Demographic analyses support a single, pre‐Weichselian colonization into Norway from East Greenland, and subsequent divergence of the three populations in separate refugia. Other refugial areas are identified in North‐east Greenland, Minnesota/Michigan, Colorado and Alaska. Admixed populations in British Columbia and West Greenland indicate postglacial contact. Taken together, evidence from this study strongly indicates in situ glacial survival in Scandinavia.     

Genetic homogeneity of the Sedgling Nehalennia speciosa (Odonata: Coenagrionidae) indicates a single Würm glacial refugium and trans‐Palaearctic postglacial expansion

The phylogeographic structures of taiga species often support the hypothesis of East Palaearctic refugia for these taxa, but the phylogeographic structures of northern temperate and southern boreal bog species are still poorly understood. Therefore, we analysed the genetic diversity and differentiation of a stenotopic damselfly, Nehalennia speciosa, across its trans‐Palaearctic range by means of sequencing two mitochondrial gene fragments, 16S rRNA‐ND1 and cytochrome c oxidase II. Only four single nucleotide polymorphisms were detected over the 1130 sequenced nucleotides. This low genetic diversity and differentiation and thus the lack of phylogeographic structure imply postglacial expansion from a single Würm Ice Age refugium, most likely located in the Far East of Asia, i.e. Manchurian refugium. From here, the species could have colonized large parts of the Palaearctics, including Europe, during the postglacial.     

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.     

Evolutionary history of an alpine shrub Hippophae tibetana (Elaeagnaceae): allopatric divergence and regional expansion

Increasing evidence suggests that geological or climatic events in the past promoted allopatric speciation of alpine plants in the Qinghai‐Tibetan Plateau and adjacent region. However, few studies have been undertaken to examine whether such allopatric divergences also occurred within a morphologically uniform species. In the present study, we report the evolutionary history of an alpine shrub species, Hippophae tibetana, based on examining chloroplast DNA (cpDNA) and nuclear ribosomal internal transcribed spacer (ITS) DNA variations. We sequenced two cpDNA fragments (trnL‐F and trnS‐G) and the nuclear ITS region in 183 individuals collected from 21 natural populations. Ten chlorotypes and 17 ITS types were identified. Phylogenetic analyses of both chlorotypes and ITS sequence variations suggested two distinct lineages distributed in the eastern and western region, respectively. On the basis of the fast and low plant substitution rates, these two lineages were estimated to have diverged from each other between 1 and 4 million years ago, during the period of the major glaciations and orogenic processes. In addition, ITS has undergone the accelerated evolution in two populations in the southern Himalaya isolated by the high mountains with a surprising accumulation of the private variations. The east–west split was also supported by an analysis of molecular variance, which partitioned around 91% of the total cpDNA variance between these two groups of populations. A single chlorotype was found for most populations in eastern or western region, suggesting a recent postglacial expansion within each region. Star‐phylogeny and mismatch analyses of all chlorotypes within the eastern group of populations suggested an earlier regional expansion before the Last Glacial Maximum (LGM). The local fixture of the different chlorotypes in multiple populations suggested more than one refugia remained for eastern or western region. Coalescent tests rejected the hypothesis that all current populations originated from a single refugium during the LGM. Instead, they supported hypothesis that two lineages diverged before the late Pleistocene. These findings, when taken together, suggested that this species had experienced long allopatric divergence and recent regional range expansions in response to orogenic processes and the climate changes. The evolutionary history of this shrub species highlights importance of geographical isolations to the intraspecific divergence of alpine plants occurring in the world's ruff. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 102 , 37–50.     

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.     

Distribution and population structure of the anther smut Microbotryum silenes‐acaulis parasitizing an arctic–alpine plant

Cold‐adapted organisms with current arctic–alpine distributions have persisted during the last glaciation in multiple ice‐free refugia, leaving footprints in their population structure that contrast with temperate plants and animals. However, pathogens that live within hosts having arctic–alpine distributions have been little studied. Here, we therefore investigated the geographical range and population structure of a fungus parasitizing an arctic–alpine plant. A total of 1437 herbarium specimens of the plant Silene acaulis were examined, and the anther smut pathogen Microbotryum silenes‐acaulis was present throughout the host's geographical range. There was significantly greater incidence of anther smut disease in more northern latitudes and where the host locations were less dense, indicating a major influence of environmental factors and/or host demographic structure on the pathogen distribution. Genetic analyses with seven microsatellite markers on recent collections of 195 M. silenes‐acaulis individuals revealed three main genetic clusters, in North America, northern Europe and southern Europe, likely corresponding to differentiation in distinct refugia during the last glaciation. The lower genetic diversity in northern Europe indicates postglacial recolonization northwards from southern refugia. This study combining herbarium surveys and population genetics thus uniquely reveals the effects of climate and environmental factors on a plant pathogen species with an arctic–alpine distribution.     

The molecular population structure of the tall forb Cicerbita alpina (Asteraceae) supports the idea of cryptic glacial refugia in central Europe

There is an ongoing debate about the glacial history of non‐arctic species in central and northern Europe. The two main hypotheses are: (1) postglacial colonization from refugia outside this region; (2) glacial survival in microclimatically favourable sites within the periglacial areas. In order to clarify the glacial history of a boreo‐montane tall forb, we analysed AFLPs from populations of Cicerbita alpina through most of its range (Scandinavia, the mountains of central Europe, the Alps, the Pyrenees and the Balkan Peninsula). We found a major differentiation between the Pyrenean population and all others, supported by principal coordinate, neighbour joining and STRUCTURE analyses. Furthermore, three populations from the central and north‐eastern Alps were genetically distinct from the bulk of populations from Scandinavia, central Europe, the Alps and the Balkan Peninsula. Most populations, including those from central and northern Europe, had moderate to high levels of genetic diversity (mean Shannon index H_Sh = 0.292, mean percentage of polymorphic loci P = 54.1%, mean Nei's gene diversity H = 0.195). The results indicate separate glacial refugia in the Pyrenean region and the Italian Alps. Furthermore, they provide evidence of glacial persistence in cryptic refugia north of the Alps, from where Scandinavia and most of the Alps are likely to have been colonized following deglaciation. © 2010 The Linnean Society of London, Botanical Journal of the Linnean Society, 2010, 164 , 142–154.     

Higher genetic diversity in recolonized areas than in refugia of Alnus glutinosa triggered by continent‐wide lineage admixture

Genetic admixture is supposed to be an important trigger of species expansions because it can create the potential for selection of genotypes suitable for new climatic conditions. Up until now, however, no continent‐wide population genetic study has performed a detailed reconstruction of admixture events during natural species expansions. To fill this gap, we analysed the postglacial history of Alnus glutinosa, a keystone species of European swamp habitats, across its entire distribution range using two molecular markers, cpDNA and nuclear microsatellites. CpDNA revealed multiple southern refugia located in the Iberian, Apennine, Balkan and Anatolian Peninsulas, Corsica and North Africa. Analysis of microsatellites variation revealed three main directions of postglacial expansion: (i) from the northern part of the Iberian Peninsula to Western and Central Europe and subsequently to the British Isles, (ii) from the Apennine Peninsula to the Alps and (iii) from the eastern part of the Balkan Peninsula to the Carpathians followed by expansion towards the Northern European plains. This challenges the classical paradigm that most European populations originated from refugial areas in the Carpathians. It has been shown that colonizing lineages have met several times and formed secondary contact zones with unexpectedly high population genetic diversity in Central Europe and Scandinavia. On the contrary, limited genetic admixture in southern refugial areas of A. glutinosa renders rear‐edge populations in the Mediterranean region more vulnerable to extinction due to climate change.     

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 origin of unique diversity in deglaciated areas: traces of Pleistocene processes in north‐European endemics from the Galium pusillum polyploid complex (Rubiaceae)

The role of glacial oscillations in shaping plant diversity has been only rarely addressed in endemics of formerly glaciated areas. The Galium pusillum group represents a rare example of an ecologically diverse and ploidy‐variable species complex that exhibits substantial diversity in deglaciated northern Europe. Using AFLP and plastid and nuclear DNA sequences of 67 populations from northern, central, and western Europe with known ecological preferences, we elucidate the evolutionary history of lineages restricted to deglaciated areas and identify the eco‐geographic partitioning of their genetic variation. We reveal three distinct endemic northern lineages: (i) diploids from southern Sweden + the British Isles, (ii) tetraploids from southern Scandinavia and the British Isles that show signs of ancient hybridization between the first lineage and populations from unglaciated central Europe, and (iii) tetraploids from Iceland + central Norway. Available evidence supports a stepwise differentiation of these three lineages that started at least before the last glacial maximum by processes of genome duplication, interlineage hybridization and/or allopatric evolution in distinct periglacial refugia. We reject the hypothesis of more recent postglacial speciation. Ecological characteristics of the populations under study only partly reflect genetic variation and suggest broad niches of postglacial colonizers. Despite their largely allopatric modern distributions, the north‐European lineages of the G. pusillum group do not show signs of rapid postglacial divergence, in contrast to most other northern endemics. Our study suggests that plants inhabiting deglaciated areas outside the Arctic may exhibit very different evolutionary histories compared with their more thoroughly investigated high‐arctic counterparts.     

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.     

Phylogeographic analysis reveals northerly refugia for the riverine amphibian Triturus dobrogicus (Caudata: Salamandridae)

We investigated the recent evolutionary history of the Danube crested newt, Triturus dobrogicus through reconstructions of: (1) the number and position of refugia at the last glacial maximum, (2) the role of major central European rivers in pattern of post‐glacial dispersal, and (3) the present‐day distribution pattern. We analysed sequences of mitochondrial DNA (ND2, 1065 bp) and six microsatellite loci in 363 T. dobrogicus individuals from 58 populations covering the range of the species. Our analyses suggested that T. dobrogicus survived the last glacial maximum in two separate refugia positioned in northwestern Pannonia and in Southern Pannonia from where its range expanded along the Danube and Tisza Rivers. Our findings also confirmed that rivers played an important role in shaping the evolutionary history of amphibian species in Central Europe. We compared the T. dobrogicus range with another lowland amphibian, the fire‐bellied toad Bombina bombina, using species distribution modelling. In line with these models, the isolated mountains inside Pannonia are occupied not by B. bombina, but by its mountain‐dwelling sister‐species B. variegata. However, in contrast to the model, crested newts (including T. dobrogicus) are absent from these mountains. We attribute this biogeographical discrepancy to the positioning of the species’ refugia at the last glacial maximum.     

Local endemism and within‐island diversification of shrews illustrate the importance of speciation in building Sundaland mammal diversity

Island systems are important models for evolutionary biology because they provide convenient, discrete biogeographic units of study. Continental islands with a history of intermittent dry land connections confound the discrete definitions of islands and have led zoologists to predict (i) little differentiation of terrestrial organisms among continental shelf islands and (ii) extinction, rather than speciation, to be the main cause of differences in community composition among islands. However, few continental island systems have been subjected to well‐sampled phylogeographic studies, leaving these biogeographic assumptions of connectivity largely untested. We analysed nine unlinked loci from shrews of the genus Crocidura from seven mountains and two lowland localities on the Sundaic continental shelf islands of Sumatra and Java. Coalescent species delimitation strongly supported all currently recognized Crocidura species from Sumatra (six species) and Java (five species), as well as one undescribed species endemic to each island. We find that nearly all species of Crocidura in the region are endemic to a single island and several of these have their closest relative(s) on the same island. Intra‐island genetic divergence among allopatric, conspecific populations is often substantial, perhaps indicating species‐level diversity remains underestimated. One recent (Pleistocene) speciation event generated two morphologically distinct, syntopic species on Java, further highlighting the prevalence of within‐island diversification. Our results suggest that both between‐ and within‐island speciation processes generated local endemism in Sundaland, supplementing the traditional view that the region's fauna is relictual and primarily governed by extinction.     

Methods for species delimitation in bumblebees (Hymenoptera,Apidae, Bombus): towards an integrative approach

Delimitation of closely related species is often hindered by the lack of discrete diagnostic morphological characters. This is exemplified in bumblebees (genus Bombus). There have been many attempts to clarify bumblebee taxonomy by using alternative features to discrete morphological characters such as wing shape, DNA, or eco‐chemical traits. Nevertheless each approach has its own limitations. Recent studies have used a multisource approach to gather different lines of speciation evidence in order to draw a strongly supported taxonomic hypothesis in bumblebees. Yet, the resulting taxonomic status is not independent of selected evidence and of consensus methodology (i.e. unanimous procedure, majority, different weighting of evidence). In this article, we compare taxonomic conclusions for a group of taxonomically doubtful species (the Bombus lapidarius‐group) obtained from the four commonly used lines of evidence for species delimitation in bumblebees (geometric morphometric of wing shape, genetic differentiation assessment, sequence‐based species delimitation methods and differentiation of cephalic labial gland secretions). We ultimately aim to assess the usefulness of these lines of evidence as components of an integrative decision framework to delimit bumblebee species. Our results show that analyses based on wing shape do not delineate any obvious cluster. In contrast, nuclear/mitochondrial, sequence‐based species delimitation methods, and analyses based on cephalic labial gland secretions are congruent with each other. This allows setting up an integrative decision framework to establish strongly supported species and subspecies status within bumblebees.     

Genomic data reveal deep genetic structure but no support for current taxonomic designation in a grasshopper species complex

Taxonomy has traditionally relied on morphological and ecological traits to interpret and classify biological diversity. Over the last decade, technological advances and conceptual developments in the field of molecular ecology and systematics have eased the generation of genomic data and changed the paradigm of biodiversity analysis. Here we illustrate how traditional taxonomy has led to species designations that are supported neither by high throughput sequencing data nor by the quantitative integration of genomic information with other sources of evidence. Specifically, we focus on Omocestus antigai and Omocestus navasi, two montane grasshoppers from the Pyrenean region that were originally described based on quantitative phenotypic differences and distinct habitat associations (alpine vs. Mediterranean‐montane habitats). To validate current taxonomic designations, test species boundaries, and understand the factors that have contributed to genetic divergence, we obtained phenotypic (geometric morphometrics) and genome‐wide SNP data (ddRADSeq) from populations covering the entire known distribution of the two taxa. Coalescent‐based phylogenetic reconstructions, integrative Bayesian model‐based species delimitation, and landscape genetic analyses revealed that populations assigned to the two taxa show a spatial distribution of genetic variation that do not match with current taxonomic designations and is incompatible with ecological/environmental speciation. Our results support little phenotypic variation among populations and a marked genetic structure that is mostly explained by geographic distances and limited population connectivity across the abrupt landscapes characterizing the study region. Overall, this study highlights the importance of integrative approaches to identify taxonomic units and elucidate the evolutionary history of species.     

Hitchhiking with forests: population genetics of the epiphytic lichen Lobaria pulmonaria in primeval and managed forests in southeastern Europe

Availability of suitable trees is a primary determinant of range contractions and expansions of epiphytic species. However, switches between carrier tree species may blur co‐phylogeographic patterns. We identified glacial refugia in southeastern Europe for the tree‐colonizing lichen Lobaria pulmonaria, studied the importance of primeval forest reserves for the conservation of genetically diverse populations and analyzed differences in spatial genetic structure between primeval and managed forests with fungus‐specific microsatellite markers. Populations belonged to either of two genepools or were admixed. Gene diversity was higher in primeval than in managed forests. At small distances up to 170 m, genotype diversity was lower in managed compared with primeval forests. We found significant associations between groups of tree species and two L. pulmonaria genepools, which may indicate “hitchhiking” of L. pulmonaria on forest communities during postglacial migration. Genepool B of L. pulmonaria was associated with European Beech (Fagus sylvatica) and we can hypothesize that genepool B survived the last glaciation associated within the refuge of European Beech on the Coastal and Central Dinarides. The allelic richness of genepool A was highest in the Alps, which is the evidence for a northern refuge of L. pulmonaria. Vicariant altitudinal distributions of the two genepools suggest intraspecific ecological differentiation.