Resistance or emigration: response of the high-alpine plant Eritrichium nanum (L.) Gaudin to the ice age within the Central Alps

Addressing alpine plant phylogeography using integrative distributional, demographic and coalescent modeling

Phylogeographic studies of alpine plants have evolved considerably in the last two decades from ad hoc interpretations of genetic data to statistical model-based approaches. In this review we outline the developments in alpine plant phylogeography focusing on the recent approach of integrative distributional, demographic and coalescent (iDDC) modeling. By integrating distributional data with spatially explicit demographic modeling and subsequent coalescent simulations, the history of alpine species can be inferred and long-standing hypotheses, such as species-specific responses to climate change or survival on nunataks during the last glacial maximum, can be efficiently tested as exemplified by available case studies. We also discuss future prospects and improvements of iDDC.

An explicit test of Pleistocene survival in peripheral versus nunatak refugia in two high mountain plant species

Pleistocene climate fluctuations had profound influence on the biogeographical history of many biota. As large areas in high mountain ranges were covered by glaciers, biota were forced either to peripheral refugia (and possibly beyond to lowland refugia) or to interior refugia (nunataks). However, nunatak survival remains controversial as it relies solely on correlative genetic evidence. Here, we test hypotheses of glacial survival using two high alpine plant species (the insect‐pollinated Pedicularis asplenifolia and wind‐pollinated Carex fuliginosa) in the European Alps. Employing the iDDC (integrative Distributional, Demographic and Coalescent) approach, which couples species distribution modelling, spatial and temporal demographic simulation and Approximate Bayesian Computation, we explicitly test three hypotheses of glacial survival: (a) peripheral survival only, (b) nunatak survival only and (c) peripheral plus nunatak survival. In P. asplenifolia the peripheral plus nunatak survival hypothesis was supported by Bayes factors (BF> 100), whereas in C. fuliginosa the peripheral survival only hypothesis, although best supported, could not be unambiguously distinguished from the peripheral plus nunatak survival hypothesis (BF = 5.58). These results are consistent with current habitat preferences (P. asplenifolia extends to higher elevations) and the potential for genetic swamping (i.e., replacement of local genotypes via hybridization with immigrating genotypes [expected to be higher in the wind‐pollinated C. fuliginosa]). Although the persistence of plants on nunataks during glacial periods has been debated and studied over decades, this is one of the first studies to explicitly test the hypothesis instead of solely using correlative evidence.

No evidence of contemporary interploidy gene flow between the closely related European woodland violets Viola reichenbachiana and V. riviniana (sect. Viola, Violaceae)

Viola reichenbachiana (2n = 4x = 20) and V. riviniana (2n = 8x = 40) are closely related species widely distributed in Europe, often sharing the same habitat throughout their overlapping ranges. It has been suggested in numerous studies that their high intraspecific morphological variability and plasticity might have been further increased by interspecific hybridisation in contact zones, given the sympatry of the species and the incomplete sterility of their hybrid. The aims of this study were to: (i) confirm that V. reichenbachiana and V. riviniana have one 4x genome in common, and (ii) determine the impact of hybridisation and introgression on genetic variation of these two species in selected European populations. For our study, we used 31 Viola populations from four European countries, which were analysed using AFLP and sequencing of a variable plastid intergenic spacer, trnH-psbA. Our analysis revealed that V. reichenbachiana exhibited larger haplotype diversity, having three species-specific haplotypes versus one in V. riviniana. The relationships among haplotypes suggest transfer of common haplotypes into V. riviniana from both V. reichenbachiana and hypothetically the other, now extinct, parental species. AFLP analysis showed low overall genetic diversity of both species, with V. riviniana showing higher among-population diversity. None of the morphologically designated hybrid populations had additive AFLP polymorphisms that would have indicated recent hybridisation. Also, kinship coefficients between both species did not indicate gene flow. V. riviniana showed significant population subdivision and significant isolation by distance, in contrast to V. reichenbachiana. The results indicate lack of gene flow between species, high influence of selfing on genetic variability, as well as probably only localised introgression toward V. riviniana.

Long-term demographic decline and late glacial divergence in a Californian paleoendemic: Sequoiadendron giganteum (giant sequoia)

Mediterranean ecosystems comprise a high proportion of endemic taxa whose response to climate change will depend on their evolutionary origins. In the California flora, relatively little attention has been given to the evolutionary history of paleoendemics from a molecular perspective, yet they number among some of the world's most iconic plant species. Here, we address questions of demographic change in Sequoiadendron giganteum (giant sequoia) that is restricted to a narrow belt of groves in the Sierra Nevada Mountains. We ask whether the current distribution is a result of northward colonization since the last glacial maximum (LGM), restriction of a broader range in the recent past (LGM) or independent colonizations in the deeper past. Genetic diversity at eleven microsatellite loci decreased with increasing latitude, but partial regressions suggested this was a function of smaller population sizes in the north. Disjunct populations north of the Kings River were divergent from those south of the Kings River that formed a single cluster in Bayesian assignment tests. Demographic inferences supported a demographic contraction just prior to the LGM as the most likely scenario for the current disjunct range of the species. This contraction appeared to be superimposed upon a long‐term decline in giant sequoia over the last 2 million years, associated with increasing aridity due to the Mediterranean climate. Overall, low genetic diversity, together with competition in an environment to which giant sequoia is likely already poorly adapted, will pose major constraints on its success in the face of increasing aridity.

The genetic diversity and spatial genetic structure of the Corso-Sardinian endemic Ferula arrigonii Bocchieri (Apiaceae)

Corsica and Sardinia represent major hotspots of plant diversity in the Mediterranean area and are priority regions for conservation due to their high number of endemic plant species. However, information supporting human decision-making on the conservation of these species is still scarce, especially at the genetic level. In this work, the first assessment is reported of the species-wide spatial genetic structure and diversity of Ferula arrigonii Bocchieri, a Corso-Sardinian endemic located in a few coastal sites and on small islands. Nine populations covering the entire natural range of the species were investigated by means of AFLP (amplified fragment length polymorphism) markers. Results indicate that this species is characterised by high levels of genetic polymorphism (92% polymorphic fragments) and of genetic diversity (H(w) = 0.317) and by relatively low differentiation among populations (F(st) = 0.057). PCoA, Bayesian analysis and neighbour-joining clustering were also employed to investigate the genetic structure of this species. Three genetically distinct groups were detected, although with considerable overlap between populations.

Climate rather than geography separates two European honeybee subspecies

Both climatic and geographical factors play an important role for the biogeographical distribution of species. The Carpathian mountain ridge has been suggested as a natural geographical divide between the two honeybee subspecies Apis mellifera carnica and A. m. macedonica. We sampled one worker from one colony each at 138 traditional apiaries located across the Carpathians spanning from the Hungarian plains to the Danube delta. All samples were sequenced at the mitochondrial tRNA(Leu)-cox2 intergenic region and genotyped at twelve microsatellite loci. The Carpathians had only limited impact on the biogeography because both subspecies were abundant on either side of the mountain ridge. In contrast, subspecies differentiation strongly correlated with the various temperature zones in Romania. A. m. carnica is more abundant in regions with the mean average temperature below 9 °C, whereas A. m. macedonica honeybees are more frequent in regions with mean temperatures above 9 °C. This range selection may have impact on the future biogeography in the light of anticipated global climatic changes.

Population genetics of the westernmost distribution of the glaciations-surviving black truffle Tuber melanosporum

The black truffle (Tuber melanosporum Vittad.) is an important natural resource due to its relevance as a delicacy in gastronomy. Different aspects of this hypogeous fungus species have been studied, including population genetics of French and Italian distribution ranges. Although those studies include some Spanish populations, this is the first time that the genetic diversity and genetic structure of the wide geographical range of the natural Spanish populations have been analysed. To achieve this goal, 23 natural populations were sampled across the Spanish geographical distribution. ISSR technique demonstrated its reliability and capability to detect high levels of polymorphism in the species. Studied populations showed high levels of genetic diversity (h N  = 0.393, h S  = 0.678, Hs = 0.418), indicating a non threatened genetic conservation status. These high levels may be a consequence of the wide distribution range of the species, of its spore dispersion by animals, and by its evolutionary history. AMOVA analysis showed a high degree of genetic structure among populations (47.89%) and other partitions as geographical ranges. Bayesian genetic structure analyses differentiated two main Spanish groups separated by the Iberian Mountain System, and showed the genetic uniqueness of some populations. Our results suggest the survival of some of these populations during the last glaciation, the Spanish southern distribution range perhaps surviving as had occurred in France and Italy, but it is also likely that specific northern areas may have acted as a refugia for the later dispersion to other calcareous areas in the Iberian Peninsula and probably France.

Genetic roadmap of the Arctic: plant dispersal highways, traffic barriers and capitals of diversity

We provide the first comparative multispecies analysis of spatial genetic structure and diversity in the circumpolar Arctic using a common strategy for sampling and genetic analyses. We aimed to identify and explain potential general patterns of genetic discontinuity/connectivity and diversity, and to compare our findings with previously published hypotheses. We collected and analyzed 7707 samples of 17 widespread arctic-alpine plant species for amplified fragment length polymorphisms (AFLPs). Genetic structure, diversity and distinctiveness were analyzed for each species, and extrapolated to cover the geographic range of each species. The resulting maps were overlaid to produce metamaps. The Arctic and Atlantic Oceans, the Greenlandic ice cap, the Urals, and lowland areas between southern mountain ranges and the Arctic were the strongest barriers against gene flow. Diversity was highest in Beringia and gradually decreased into formerly glaciated areas. The highest degrees of distinctiveness were observed in Siberia. We conclude that large-scale general patterns exist in the Arctic, shaped by the Pleistocene glaciations combined with long-standing physical barriers against gene flow. Beringia served as both refugium and source for interglacial (re)colonization, whereas areas further west in Siberia served as refugia, but less as sources for (re)colonization.

Evolutionary history of almond tree domestication in the Mediterranean basin

Genetic diversity of contemporary domesticated species is shaped by both natural and human-driven processes. However, until now, little is known about how domestication has imprinted the variation of fruit tree species. In this study, we reconstruct the recent evolutionary history of the domesticated almond tree, Prunus dulcis, around the Mediterranean basin, using a combination of nuclear and chloroplast microsatellites [i.e. simple sequence repeat (SSRs)] to investigate patterns of genetic diversity. Whereas conservative chloroplast SSRs show a widespread haplotype and rare locally distributed variants, nuclear SSRs show a pattern of isolation by distance with clines of diversity from the East to the West of the Mediterranean basin, while Bayesian genetic clustering reveals a substantial longitudinal genetic structure. Both kinds of markers thus support a single domestication event, in the eastern side of the Mediterranean basin. In addition, model-based estimation of the timing of genetic divergence among those clusters is estimated sometime during the Holocene, a result that is compatible with human-mediated dispersal of almond tree out of its centre of origin. Still, the detection of region-specific alleles suggests that gene flow from relictual wild preglacial populations (in North Africa) or from wild counterparts (in the Near East) could account for a fraction of the diversity observed.

Pleistocene survival on central Alpine nunataks: genetic evidence from the jumping bristletail Machilis pallida

Mechanisms of survival during the Pleistocene glaciation periods have been studied for more than a century. Until now, molecular studies that confirmed animal survival on Alpine nunataks, that is, ice-free summits surrounded by glaciers, were restricted to peripheral areas. Here, we search for molecular signatures of inner-Alpine survival of the narrow endemic and putatively parthenogenetic Alpine jumping bristletail Machilis pallida combining mitochondrial and AFLP data from its three known populations. The mitochondrial data indicate survival on both peripheral and central nunataks, the latter suggesting that refugia in the centre of the Alpine main ridge were more widespread than previously recognized. Incongruences between mitochondrial and AFLP patterns suggest a complex evolutionary history of the species and may be explained via parallel fixation of parthenogenesis of different origins during the last glacial maximum. We suggest that the inferred parthenogenesis may have been essential for central nunatak survival, but may pose a serious threat for M. pallida in consideration of the present climatic changes.

Integrating species distribution models (SDMs) and phylogeography for two species of Alpine Primula

The major intention of the present study was to investigate whether an approach combining the use of niche-based palaeodistribution modeling and phylo-geography would support or modify hypotheses about the Quaternary distributional history derived from phylogeographic methods alone. Our study system comprised two closely related species of Alpine Primula. We used species distribution models based on the extant distribution of the species and last glacial maximum (LGM) climate models to predict the distribution of the two species during the LGM. Phylogeographic data were generated using amplified fragment length polymorphisms (AFLPs). In Primula hirsuta, models of past distribution and phylogeographic data are partly congruent and support the hypothesis of widespread nunatak survival in the Central Alps. Species distribution models (SDMs) allowed us to differentiate between alpine regions that harbor potential nunatak areas and regions that have been colonized from other areas. SDMs revealed that diversity is a good indicator for nunataks, while rarity is a good indicator for peripheral relict populations that were not source for the recolonization of the inner Alps. In P. daonensis, palaeo-distribution models and phylogeographic data are incongruent. Besides the uncertainty inherent to this type of modeling approach (e.g., relatively coarse 1-km grain size), disagreement of models and data may partly be caused by shifts of ecological niche in both species. Nevertheless, we demonstrate that the combination of palaeo-distribution modeling with phylogeographical approaches provides a more differentiated picture of the distributional history of species and partly supports (P. hirsuta) and partly modifies (P. daonensis and P. hirsuta) hypotheses of Quaternary distributional history. Some of the refugial area indicated by palaeodistribution models could not have been identified with phylogeographic data.

Extensive range persistence in peripheral and interior refugia characterizes Pleistocene range dynamics in a widespread Alpine plant species (Senecio carniolicus, Asteraceae)

Recent evidence suggests that survival of arctic-alpine organisms in peripheral or interior glacial refugia are not mutually exclusive and may both be involved in shaping an organism’s Pleistocene history, yet potentially at different time levels. Here, we test this hypothesis in a high-mountain plant (diploid lineage of Senecio carniolicus, Asteraceae) from the Eastern European Alps, in which patterns of morphological variation and current habitat requirements suggest survival in both types of refugia. To this end, we used AFLPs, nuclear and plastid DNA sequences and analysed them, among others, within a graph theoretic framework and using novel Bayesian methods of phylogeographic inference. On the basis of patterns of genetic diversity, occurrence of rare markers, distribution of distinct genetic lineages and patterns of range connectivity both interior refugia in the formerly strongly glaciated central Alps and peripheral refugia along the southern margin of the Alps were identified. The presence of refugia congruently inferred by markers resolving at different time levels suggests that these refugia acted as such throughout several glacial cycles. The high degree of range persistence together with gradual range expansion, which contrasts with the extent of range shifts implied for other Alpine species, is likely responsible for incipient lineage differentiation evident from the genetic data. Replacing a simplistic peripheral vs. interior refugia dualism by more complex models involving both types of refugia and considering different time levels will help identifying common phylogeographic patterns with respect to, for instance, location of refugia and colonization routes and elucidating their underlying genetic and/or ecological causes.

The phylogeographical history of the Iberian steppe plant Ferula loscosii (Apiaceae): a test of the abundant-centre hypothesis

The geology and climate of the western Mediterranean area were strongly modified during the Late Tertiary and the Quaternary. These geological and climatic events are thought to have induced changes in the population histories of plants in the Iberian Peninsula. However, fine-scale genetic spatial architecture across western Mediterranean steppe plant refugia has rarely been investigated. A population genetic analysis of amplified fragment length polymorphism variation was conducted on present-day, relict populations of Ferula loscosii (Apiaceae). This species exhibits high individual/population numbers in the middle Ebro river valley and, according to the hypothesis of an abundant-centre distribution, these northern populations might represent a long-standing/ancestral distribution centre. However, our results suggest that the decimated southern and central Iberian populations are more variable and structured than the northeastern ones, representing the likely vestiges of an ancestral distribution centre of the species. Phylogeographical analysis suggests that F. loscosii likely originated in southern Spain and then migrated towards the central and northeastern ranges, further supporting a Late Miocene southern-bound Mediterranean migratory way for its oriental steppe ancestors. In addition, different glacial-induced conditions affected the southern and northern steppe Iberian refugia during the Quaternary. The contrasting genetic homogeneity of the Ebro valley range populations compared to the southern Iberian ones possibly reflects more severe bottlenecks and subsequent genetic drift experienced by populations of the northern Iberia refugium during the Pleistocene, followed by successful postglacial expansion from only a few founder plants.

Pleistocene distribution range shifts were accompanied by breeding system divergence within Hornungia alpina (Brassicaceae) in the Alps

Impact of glacial history on the phylogeography of silicate-dwelling plants of the European Alps has been particularly well studied, whereas virtually no data are available for species growing on different bedrock types, as for Hornungia alpina. Bayesian clustering of AFLP data only partly support the distinction of three subspecies as morphologically defined. Whereas the phylogeographical N-group corresponds to subsp. alpina, the congruence of the SW-group and SE-group with subsp. brevicaulis, and subsp. austroalpina, respectively, is limited. High levels of rarity and genetic diversity in the N-group suggest Pleistocene survival along the outer margin of the Alpine arc. For subsp. brevicaulis we suggest a single origin from a refugium in the Southwestern Alps, whereas subsp. austroalpina might have originated twice in the Southern and Southeastern Alps. Different levels of genetic diversity and partitioning of genetic variation indicate a divergence in breeding system, which is corroborated by pollinator exclusion experiments revealing self-incompatibility in the N-group and autonomous selfing in the SE-group.

Genetic structure and evolution of Alpine polyploid complexes: Ranunculus kuepferi (Ranunculaceae) as a case study

The alpine white-flowered buttercup, Ranunculus kuepferi Greuter & Burdet, is a polyploid complex with diploids endemic to the southwestern Alps and polyploids - which have been previously described as apomictic - widespread throughout European mountains. Due to the polymorphic status of both its ploidy level and its reproductive mode, R. kuepferi represents a key species for understanding the evolution of polyploid lineages in alpine habitats. To disentangle the phylogeography of this polyploid taxon, we used cpDNA sequences and AFLP (amplified fragment length polymorphism) markers in 33 populations of R. kuepferi representative of its ploidy level and distribution area. Polyploid individuals were shown to be the result of at least two polyploidization events that may have taken place in the southwestern Alps. From this region, one single main migration of tetraploids colonized the entire Alpine range, the Apennines and Corsica. Genetic recombination among tetraploids was also observed, revealing the facultative nature of the apomictic reproductive mode in R. kuepferi polyploids. Our study shows the contrasting role played by diploid lineages mostly restricted to persistent refugia and by tetraploids, whose dispersal abilities have permitted their range extension all over the previously glaciated Alpine area and throughout neighbouring mountain massifs.

Multiple glacial refugia and complex postglacial range shifts of the obligatory woodland plant Polygonatum verticillatum (Convallariaceae)

The phylogeography of typical alpine plant species is well understood in Europe. However, the genetic patterns of boreo-montane species are mostly unstudied. Therefore, we analysed the AFLPs of 198 individuals of Polygonatum verticillatum over a major part of its European distribution. We obtained a total of 402 reproducible fragments, of which 96.8% were polymorphic. The average Phi(ST) over all samples was high (73.0%). The highest number of private fragments was observed in the Cantabrian Mountains; the highest genetic diversities of the populations were detected in populations from the Alps. BAPS, Principal Coordinates and Cluster analyses revealed a deep split between the Cantabrian population and all other samples. The latter further distinguished two major groups in western and eastern Europe. These results suggest a complex biogeographical history of P. verticillatum. The Cantabrian population was most probably isolated for the longest time. Furthermore, putative glacial survival centres might have existed in the western group around the glaciated Alps and in the eastern group in the foothills of the Carpathian and Balkan mountain systems. The origin of the Scandinavian populations is still unresolved, but an origin from the southeastern Alps or the western Balkans appears the most likely scenario.

The homogenous genetic structure and inferred unique history of range shifts during the Pleistocene climatic oscillations of Arcterica nana (Maxim.) Makino (Ericaceae)

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

Intraspecific variation in Viola suavis in Europe: parallel evolution of white-flowered morphotypes

BACKGROUND AND AIMS

Viola species are commonly grown for their ornamental flowers, but their evolutionary history and taxonomy are often complicated and have been poorly explored so far. This is a study of the polymorphic, typically blue-flowered species Viola suavis, concentrating on the white-flowered populations of uncertain taxonomic assignment that occur in Spain and central and south-eastern Europe. The aim was to resolve their origin and taxonomic status and to study the intraspecific structure and (post)glacial history of this species.

METHODS

Viola suavis and five close relatives were sampled from multiple locations and subjected to molecular (AFLP, sequencing of nrDNA ITS) and morphometric analyses. Data on ploidy level and pollen fertility were also obtained, to address an assumed hybrid origin of the white-flowered populations.

KEY RESULTS

In V. suavis a strong intraspecific genetic split into two groups was observed, indicating that there has been a long-term isolation and survival in distinct glacial refugia. The white-flowered populations could be placed within the variation range of this species, and it is clear that they evolved independently in two distant areas. Their parallel evolution is supported by both morphological and genetic differentiation. The strongly reduced genetic variation and absence of unique AFLP fragments suggest their derived status and origin from the typical, blue-flowered populations.

CONCLUSIONS

These results suggest that intraspecific variation in V. suavis has been largely shaped by population isolations during the last glaciation and subsequent recolonizations, although cultivation and vegetative spread by humans have affected the present picture as well.

Conservation genetics of the endangered Iberian steppe plant Ferula loscosii (Apiaceae)

Ferula loscosii (Lange) Willk (Apiaceae) is a threatened endemic species native to the Iberian Peninsula. The plant has a narrow and disjunct distribution in three regions, NE, C and SE Spain. Genetic variability within and among 11 populations from its natural distribution was assessed using allozymes. Intermediate levels of genetic diversity were detected in F. loscosii (P(99%) = 36.83; H(E) = 0.125; H(T) = 0.152). However, the highest genetic diversity (58%) corresponded to the threatened populations from SE and C Spain (H(T) = 0.169) rather than the more abundant and larger populations from NE Spain (Ebro valley) (H(T) = 0.122). Low to moderate levels of genetic structure were found among regional ranges (G(ST) = 0.134), and several statistical spatial correlation analyses corroborated substantial genetic differentiation among the three main regional ranges. However, no significant genetic differentiation was found among the NE Spain populations, except for a northernmost population that is geographically isolated. Outcrossing mating and other biological traits of the species could account for the maintenance of the present values of genetic diversity within populations. The existence of an ancestral late Tertiary wider distribution of the species in SE and C Spain, followed by the maintenance of different Quaternary refugia in these warmer areas, together with a more recent and rapid post-glacial expansion towards NE Spain, are arguments that could explain the low genetic variability and structure found in the Ebro valley and the higher levels of diversity in the southern Iberian populations.

The genetic structure of the mountain forest butterfly Erebia euryale unravels the late Pleistocene and postglacial history of the mountain coniferous forest biome in Europe

The distribution of the mountain coniferous forest biome in Europe throughout time is not sufficiently understood. One character species of this habitat type is the large ringlet, Erebia euryale well reflecting the extension of this biome today, and the genetic differentiation of this species among and within mountain systems may unravel the late Pleistocene history of this habitat type. We therefore analysed the allozyme pattern of 381 E. euryale individuals from 11 populations in four different European mountain systems (Pyrenees, Alps, Carpathians, Rila). All loci analysed were polymorphic. The mean F(ST) over all samples was high (20%). Furthermore, the mean genetic distance among samples was quite high (0.049). We found four different groups well supported by cluster analyses, bootstraps and hierarchical variance analyses: Pyrenees, western Alps, eastern Alps and southeastern Europe (Carpathians and Rila). The genetic diversity of the populations was highest in the southeastern European group and stepwise decreased westwards. Interestingly, the populations from Bulgaria and Romania were almost identical; therefore, we assume that they were not separated by the Danube Valley, at least during the last ice age. On the contrary, the differentiation among the three western Alps populations was considerable. For all these reasons, we assume that (i) the most important refugial area for the coniferous mountain forest biome in Europe has been located in southeastern Europe including at least parts of the Carpathians and the Bulgarian mountains; (ii) important refugial areas for this biome existed at the southeastern edge of the Alps; (iii) fragments of this habitat types survived along the southwestern Alps, but in a more scattered distribution; and (iv) relatively small relicts have persisted somewhere at the foothills of the Pyrenees.

Cryptic differentiation in alpine-endemic, high-altitude butterflies reveals down-slope glacial refugia

The influence of cyclic climate fluctuations and their impact on high-altitude species is still insufficiently understood. We therefore analysed in this study the genetic structure of cold-adapted animals and their coherence with geographical distributions throughout the Late Quaternary. We analysed 588 individuals from 23 populations of the alpine-endemic lesser mountain ringlet, Erebia melampus, by allozyme electrophoresis to detect its intraspecific differentiation. As an outgroup, we added one population of Erebia sudetica inalpina from Grindelwald (Swiss Alps). Seventeen of 18 loci were polymorphic. The mean F(ST) over all samples was 37%. We detected strong differentiation into three lineages with the genetic distances between the two E. melampus groups being larger than between each of the two E. melampus groups and E. sudetica. The mean genetic distance among these three groups was 0.17. These results give evidence for the existence of a species complex within the E. melampus/sudetica group and indicate a discontinuous distribution within this group during at least the last ice age. One of them, E. sudetica inalpina, is found in the northern Alps and most probably had its Würm glacial refugium north of the glaciated Alps. The western E. melampus group might have had a refugium at the southwestern Alps margin, the eastern group in the lower altitudes of the southeastern and/or eastern Alps. In the latter, a further subdivision within this relict area is possible.

Glacial in situ survival in the Western Alps and polytopic autopolyploidy in Biscutella laevigata L. (Brassicaceae)

Past climatic changes and especially the ice ages have had a great impact on both the distribution and the genetic composition of plant populations, but whether they promoted speciation is still controversial. The autopolyploid complex Biscutella laevigata is a classical example of polyploidy linked to glaciations and is an interesting model to explore migration and speciation driven by climate changes in a complex alpine landscape. Diploid taxa survived the last glacial maximum in several never-glaciated areas and autotetraploids are clearly dominant in the central parts of the Alps; however, previous range-wide studies failed to identify their diploid ancestor(s). This study highlights the phylogeographical relationships of maternal lineages in the Western Alps and investigates the polyploidy process using plastid DNA sequences (trnS-trnG and trnK-intron) combined with plastid DNA length polymorphism markers, which were transferable among Brassicaceae species. Twenty-one distinct plastid DNA haplotypes were distinguished in 67 populations densely sampled in the Western Alps and main lineages were identified by a median-joining network. The external Alps harboured high levels of genetic diversity, while the Central Alps contained only a subset of haplotypes due to postglacial recolonization. Several haplotypes were restricted to local peripheral refugia and evidence of in situ survival in central nunataks was detected by the presence of highly differentiated haplotypes swamped by frequent ones. As hierarchical genetic structure pointed to an independent evolution of the species in different biogeographical districts, and since tetraploids displayed haplotypes belonging to different lineages restricted to either the northern or the southern parts of the Alpine chain, polytopic autopolyploidy was also apparent in the Western Alps.

Phylogeography of the high alpine plant Senecio halleri (Asteraceae) in the European Alps: in situ glacial survival with postglacial stepwise dispersal into peripheral areas

Whether alpine plant species survived Pleistocene glaciations in situ on high alpine nunatak mountains is still under debate. To test this hypothesis, Senecio halleri, a high alpine and endemic species with a narrow distribution range in the European Alps, was chosen as a model organism. Polymerase chain reaction-restriction fragment length polymorphisms of chloroplast DNA (cpDNA PCR-RFLPs) were used in a phylogeographic analysis of 14 populations of S. halleri, covering its total distribution area. The results of haplotype diversity and distribution gave evidence of in situ glacial survival on siliceous central-alpine nunatak mountains in two areas, southwest and northeast of the Aosta valley. According to the absence of genetic differentiation between these two nunatak areas (based on amova), nested clade analysis implied a history of preglacial gene flow, in situ survival and extinction of intermediate populations during glaciation and postglacial stepwise recolonization of peripheral and intermediate areas.

Genetic diversity and spatial correlation patterns unravel the biogeographical history of the European sweet vernal grasses (Anthoxanthum L., Poaceae)

Different processes have contributed to shaping the present distribution of the European biotas. Up to three different tertiary- to quaternary-time-scale evolutionary scenarios have been proposed to interpret the divergence and genetic structuring of plant species in Europe. In the present study, the Amplified Fragment Length Polymorphisms technique has been used to unravel the species and regional phylogeography of the European sweet vernal grasses (Anthoxanthum L. Poaceae). Forty-six populations belonging to all seven European species of Anthoxanthum and covering a broad geographical and ecological range were selected. Different phylogeography and population genetics diversity and structure estimates indicated a clear divergence of old Messinian Mediterranean lineages, followed by a pre-Pliocene split between Mediterranean annuals and Eurosiberian perennials and a more recent Pleistocene differentiation of Arctic-Alpine, Atlantic and Submediterranean diploid to polyploid landraces. Regional and population correlation tests between geographical and genetic distances allowed to postulate distinct pre- and post-glacial colonization pathways across Europe for the taxa of this widespread genus.

Molecular biogeography of Europe: Pleistocene cycles and postglacial trends

The climatic cycles with subsequent glacial and intergalcial periods have had a great impact on the distribution and evolution of species. Using genetic analytical tools considerably increased our understanding of these processes. In this review I therefore give an overview of the molecular biogeography of Europe. For means of simplification, I distinguish between three major biogeographical entities: (i) "Mediterranean" with Mediterranean differentiation and dispersal centres, (ii) "Continental" with extra-Mediterranean centres and (iii) "Alpine" and/or "Arctic" with recent alpine and/or arctic distribution patterns. These different molecular biogeographical patterns are presented using actual examples. Many "Mediterranean" species are differentiated into three major European genetic lineages, which are due to glacial isolation in the three major Mediterranean peninsulas. Postglacial expansion in this group of species is mostly influenced by the barriers of the Pyrenees and the Alps with four resulting main patterns of postglacial range expansions. However, some cases are known with less than one genetic lineage per Mediterranean peninsula on the one hand, and others with a considerable genetic substructure within each of the Mediterranean peninsulas, Asia Minor and the Maghreb. These structures within the Mediterranean sub-centres are often rather strong and in several cases even predate the Pleistocene. For the "Continental" species, it could be shown that the formerly supposed postglacial spread from eastern Palearctic expansion centres is mostly not applicable. Quite the contrary, most of these species apparently had extra-Mediterranean centres of survival in Europe with special importance of the perialpine regions, the Carpathian Basin and parts of the Balkan Peninsula. In the group of "Alpine" and/or "Arctic" species, several molecular biogeographical patterns have been found, which support and improve the postulates based on distribution patterns and pollen records. Thus, genetic studies support the strong linkage between southwestern Alps and Pyrenees, northeastern Alps and Carpathians as well as southeastern Alps and the Dinaric mountain systems, hereby allowing conclusions on the glacial distribution patterns of these species. Furthermore, genetic analyses of arctic-alpine disjunct species support their broad distribution in the periglacial areas at least during the last glacial period. The detailed understanding of the different phylogeographical structures is essential for the management of the different evolutionary significant units of species and the conservation of their entire genetic diversity. Furthermore, the distribution of genetic diversity due to biogeographical reasons helps understanding the differing regional vulnerabilities of extant populations.

Molecular biogeography of Europe: Pleistocene cycles and postglacial trends

The climatic cycles with subsequent glacial and intergalcial periods have had a great impact on the distribution and evolution of species. Using genetic analytical tools considerably increased our understanding of these processes. In this review I therefore give an overview of the molecular biogeography of Europe. For means of simplification, I distinguish between three major biogeographical entities: (i) "Mediterranean" with Mediterranean differentiation and dispersal centres, (ii) "Continental" with extra-Mediterranean centres and (iii) "Alpine" and/or "Arctic" with recent alpine and/or arctic distribution patterns. These different molecular biogeographical patterns are presented using actual examples.

Many "Mediterranean" species are differentiated into three major European genetic lineages, which are due to glacial isolation in the three major Mediterranean peninsulas. Postglacial expansion in this group of species is mostly influenced by the barriers of the Pyrenees and the Alps with four resulting main patterns of postglacial range expansions. However, some cases are known with less than one genetic lineage per Mediterranean peninsula on the one hand, and others with a considerable genetic substructure within each of the Mediterranean peninsulas, Asia Minor and the Maghreb. These structures within the Mediterranean sub-centres are often rather strong and in several cases even predate the Pleistocene.

For the "Continental" species, it could be shown that the formerly supposed postglacial spread from eastern Palearctic expansion centres is mostly not applicable. Quite the contrary, most of these species apparently had extra-Mediterranean centres of survival in Europe with special importance of the perialpine regions, the Carpathian Basin and parts of the Balkan Peninsula.

In the group of "Alpine" and/or "Arctic" species, several molecular biogeographical patterns have been found, which support and improve the postulates based on distribution patterns and pollen records. Thus, genetic studies support the strong linkage between southwestern Alps and Pyrenees, northeastern Alps and Carpathians as well as southeastern Alps and the Dinaric mountain systems, hereby allowing conclusions on the glacial distribution patterns of these species. Furthermore, genetic analyses of arctic-alpine disjunct species support their broad distribution in the periglacial areas at least during the last glacial period.

The detailed understanding of the different phylogeographical structures is essential for the management of the different evolutionary significant units of species and the conservation of their entire genetic diversity. Furthermore, the distribution of genetic diversity due to biogeographical reasons helps understanding the differing regional vulnerabilities of extant populations.

Reticulate phylogenetics and phytogeographical structure of Heliosperma (Sileneae, Caryophyllaceae) inferred from chloroplast and nuclear DNA sequences

The Balkan Peninsula is known to be one of the most diverse and species-rich parts of Europe, but its biota has gained much less attention in phylogenetic and evolutionary studies compared to other southern European mountain systems. We used nuclear ribosomal internal transcribed spacer (ITS) sequences and intron sequences of the chloroplast gene rps16 to examine phylogenetic and biogeographical patterns within the genus Heliosperma (Sileneae, Caryophyllaceae). The ITS and rps16 intron sequences both support monophyly of Heliosperma, but the data are not conclusive with regard to its exact origin. Three strongly supported clades are found in both data sets, corresponding to Heliosperma alpestre, Heliosperma macranthum and the Heliosperma pusillum clade, including all other taxa. The interrelationships among these three differ between the nuclear and the plastid data sets. Hierarchical relationships within the H. pusillum clade are poorly resolved by the ITS data, but the rps16 intron sequences form two well-supported clades which are geographically, rather than taxonomically, correlated. A similar geographical structure is found in the ITS data, when analyzed with the NeighbourNet method. The apparent rate of change within Heliosperma is slightly higher for rps16 as compared to ITS. In contrast, in the Sileneae outgroup, ITS substitution rates are more than twice as high as those for rps16, a situation more in agreement with what has been found in other rate comparisons of noncoding cpDNA and ITS. Unlike most other Sileneae ITS sequences, the H. pusillum group sequences display extensive polymorphism. A possible explanation to these patterns is extensive hybridization and gene flow within Heliosperma, which together with concerted evolution may have eradicated the ancient divergence suggested by the rps16 data. The morphological differentiation into high elevation, mainly widely distributed taxa, and low elevation narrow endemics is not correlated with the molecular data, and is possibly a result of ecological differentiation.

Population genetic consequences of geographic disjunction: a prairie plant isolated on Great Lakes alvars

Species may often exhibit geographic variation in population genetic structure due to contemporary and historical variation in population size and gene flow. Here, we test the predictions that populations on the margins of a species' distribution contain less genetic variation and are more differentiated than populations towards the core of the range by comparing patterns of genetic variation at five microsatellite loci between disjunct and core populations of the perennial, allohexaploid herb Geum triflorum. We sampled nine populations isolated on alvar habitat within the eastern Great Lakes region in North America, habitats that include disjunct populations of several plant species, and compared these to 16 populations sampled from prairie habitat throughout the core of the species' distribution in midwestern Canada and the USA. Alvar populations exhibited much lower within-population diversity and contained only a subset of alleles found in prairie populations. We detected isolation by distance across the species' range and within alvar and prairie regions separately. As predicted, genetic differentiation was higher among alvar populations than among prairie populations, even after controlling for the geographic distance between sampled populations. Low diversity and high differentiation can be accounted for by the greater contemporary spatial isolation of alvar populations. However, the genetic structure of alvar populations may also have been influenced by postglacial range expansion and contraction. Our results are consistent with alvar populations being founded during an expansion of prairie habitat during the warmer, hypsithermal period approximately 5000 bp and subsequently becoming stranded on isolated alvar habitat as the climate grew cooler and wetter.

Refugia of Potentilla matsumurae (Rosaceae) located at high mountains in the Japanese archipelago

Molecular phylogeographic studies have revealed the genetic patterns and glacial-interglacial history of many plant and animal species. To infer the Quaternary history of alpine plants in the Japanese archipelago, which is poorly known, we investigated 203 individuals of Potentilla matsumurae and its varieties from 22 populations. We found 11 haplotypes based on approximately 1400 bp of two intergenetic spacers in chloroplast DNA (trnT-L and rpl20-rps20). The distribution of these haplotypes was geographically structured, which was supported by haplotype composition, principal component analysis, and unweighted pair group method with arithmetic mean (UPGMA), and N(ST) (0.71) was significantly greater than G(ST) (0.68). In addition to the positive correlation between genetic and geographic distance (Mantel test, r = 0.497, P < 0.001), an abrupt genetic change was detected between mountains in central Honshu and the Tohoku region. This genetic boundary was further supported by analysis of molecular variance (AMOVA), and high variation (54.0%) was explained by differences on either side of this boundary. Moreover, haplotypes in central Honshu were thought to have diverged, based on an outgroup comparison. These results suggest that mountains in central Honshu served as refugia during the Quaternary climatic oscillation, although the results could not reveal the history of most mountains in the Tohoku region and Hokkaido. Nevertheless, following floristic studies, our results indicate that alpine plants in Japan experienced a history different from that in Europe; i.e. they retreated into refugia during warm periods to avoid forest development, rather than glaciers.

Palaeopolyploidy, spatial structure and conservation genetics of the narrow steppe plant Vella pseudocytisus subsp. paui (Vellinae, Cruciferae)

BACKGROUND AND AIMS

Vella pseudocytisus subsp. paui (Cruciferae) is a narrow endemic plant to the Teruel province (eastern Spain), which is listed in the National Catalogue of Endangered Species. Two distinct ploidy levels (diploid, 2n = 34, and tetraploid, 2n = 68) have been reported for this taxon that belongs to the core subtribe Vellinae, a western Mediterranean group of shrubby taxa with a chromosome base number of x = 17. Allozyme and AFLP analyses were conducted (a) to test for the ploidy and putative palaeo-allopolyploid origin of this taxon, (b) to explore levels of genetic diversity and spatial structure of its populations, and (c) to address in-situ and ex-situ strategies for its conservation.

METHODS

Six populations that covered the entire geographical range of this taxon were sampled and examined for 19 allozyme loci and three AFLP primer pair combinations. In addition, the gametic progenies of five individuals were analysed for two allozyme loci that showed fixed heterozygosity.

KEY RESULTS

Multiple banded allozyme profiles for most of the surveyed loci indicated the polyploidy of this taxon. Co-inherited fixed heterozygous patterns were exhibited by the gametophytic tissues of the mother plants. Both allozyme and AFLP markers detected high levels of genetic diversity, and a strong micro-spatial genetic structure was recovered from AFLP phenetic analyses and Mantel correlograms.

CONCLUSIONS

Allozyme data support the hypothesis of an allotetraploid origin of Vella pseudocytisus subsp. paui that could be representative of other taxa of the core Vellinae group. AFLP data distinguished three geographically distinct groups with no genetic interaction among them. Allotetraploidy and outcrossing reproduction have probably contributed to maintenance of high levels of genetic variability of the populations, whereas habitat fragmentation may have enhanced the high genetic isolation observed among groups. In-situ microgenetic reserves and a selective sampling of germplasm stocks for ex-situ conservation of this taxon are proposed.

Molecular evidence for glacial refugia of mountain plants in the European Alps

Many mountain ranges have been strongly glaciated during the Quaternary ice ages, and the locations of glacial refugia of mountain plants have been debated for a long time. A series of detailed molecular studies, investigating intraspecific genetic variation of mountain plants in the European Alps, now allows for a first synopsis. A comparison of the phylogeographic patterns with geological and palaeoenvironmental data demonstrates that glacial refugia were located along the southwestern, southern, eastern and northern border of the Alps. Additional glacial refugia were present in central Alpine areas, where high-elevation plants survived the last glaciation on ice-free mountain tops. The observed intraspecific phylogeographies suggest general patterns of glacial survival, which conform to well-known centres of Alpine species diversity and endemism. This implies that evolutionary or biogeographic processes induced by climatic fluctuations act on gene and species diversity in a similar way.

RAPD Profiling in Detecting Genetic Variation in Endemic Coelonema (Brassicaceae) of Qinghai-Tibet Plateau of China

Random amplified polymorphic DNA (RAPD) markers were used to measure genetic diversity of Coelonema draboides (Brassicaceae), a genus endemic to the Qilian Mountains of the Qinghai-Tibet Plateau. We sampled 90 individuals in 30 populations of Coelonema draboides from Datong and Huzhu counties of Qinghai Province in P.R. China. A total of 186 amplified bands were scored from the 14 RAPD primers, with a mean of 13.3 amplified bands per primer, and 87% (161 bands) polymorphic bands (PPB) was found. Analysis of molecular variance (AMOVA) shows that a large proportion of genetic variation (84.2%) resides among individuals within populations, while only 15.8% resides among populations. The species shows higher genetic diversity between individuals than other endemic and endangered plants. The RAPDs provide a useful tool for assessing genetic diversity of rare, endemic species and for resolving relationships among populations. The results show that the genetic diversity of this species is high, possibly allowing it to adapt more easily to environmental variations. The main factor responsible for the high level of differentiation within populations and the low level of diversity among populations is probably the outcrossing and long-lived nature of this species. Some long-distance dispersal, even among far separated populations, is also a crucial determinant for the pattern of genetic variation in the species. This distributive pattern of genetic variation of C. draboides populations provides important baseline data for conservation and collection strategies for the species. It is suggested that only populations in different habitats should be studied and protected, not all populations, so as to retain as much genetic diversity as possible.

Sex-specific SCAR markers in the dioecious plant Rumex nivalis (Polygonaceae) and implications for the evolution of sex chromosomes

We developed SCAR primers based on isolated and sequenced male-specific fragments as identified in an AFLP analysis of the dioecious plant Rumex nivalis. PCR amplification using these primers on females and males resulted in fragments exclusively present in males. Co-amplification of the nuclear rDNA internal transcribed spacer 2 together with the male-specific fragment was applied as an internal control for successful PCR reactions to avoid false-negative sex scoring. With a length of about 164 bp, the AFLP fragment was of a similar size as the tandemly arranged, repetitive sequences of 180 bp located on the Y chromosomes of Rumex acetosa. The genetic distances between the Y-chromosomal sequences of R. nivalis and R. acetosa, both members of the section Acetosa, were substantial. We found intra-individual divergence among cloned sequences of the male-specific fragment in R. nivalis. The patterns of interspecific and intra-individual sequence variation found are in accordance with proposed modes of the evolution of sex chromosomes. Y chromosomes possibly arose only once in the genus Rumex and consist mainly of heterochromatic DNA. Due to the almost complete absence of selection on them, Y chromosomes are likely to accumulate large numbers of mutations.

Comparative biogeography of the cytotypes of annual Microthlaspi perfoliatum (Brassicaceae) in Europe using isozymes and cpDNA data: refugia, diversity centers, and postglacial colonization

In the last few years, the biogeography of many European plant species has been analyzed using molecular markers, and some consistent patterns of Pleistocenic differentiation and range fluctuations have been established. These studies mostly focused on perennial herbs or woody species, rarely considering annual taxa. This study focused on the annual Microthlaspi perfoliatum, which is distributed all over Europe and comprises three cytotypes. Morphologically, these cytotypes are hard to distinguish, although, based on molecular markers, they should be treated as two different species. Diploid and polyploid cytotypes had a different biogeographical history, with distinct glacial refugia. For the polyploids, a well-known distribution pattern of relict areas was confirmed, with Iberia, Italy, and the Balkans serving as primary Pleistocene refugia. Secondary refuge areas have been detected in southeastern France and neighboring Switzerland, with closer affinity to the Iberian refugium than to any other region based on allozyme and cpDNA haplotype data. For the diploids, two refugia have been characterized, one of which is congruent to the secondary refugium of the polyploids in France and Switzerland. The second refuge of diploid populations is located in unglaciated lowland areas of East Austria and Croatia. Isozyme and cpDNA haplotype data favor a postglacial colonization of diploid populations into Germany from Austrian lowland areas along the Danube River as well as from Switzerland. This scenario is also true for polyploids in Germany, Belgium, and Sweden.

The Pleistocene glaciations and the evolutionary history of the polytypic snail species Arianta arbustorum (Gastropoda, Pulmonata, Helicidae)

The evolutionary history of the snail Arianta arbustorum is controversial. This diverse, polytypic species has two distinct forms: one, with a globular shell and closed umbilicus, is found from lowland to high altitudes; the other, with a depressed shell and open umbilicus, is found at a few scattered, high altitude localities. What is the origin of these two forms? Some believe that the depressed shell is a recent, local, ecotypic adaptation to alpine environments. Others believe that this form is a relic of an ancestral condition that may have survived the Pleistocene glaciations on nunatak-like montane refugia, while the globular shell is a derived condition and its presence at high altitudes follows post-Pleistocene recolonisation. We analysed a portion of the mitochondrial gene cytochrome oxidase I for 100 snails of the species A. arbustorum, three additional Arianta species, and nine outgroup taxa from five genera, in order to understand the phylogeographic history of the species. Despite some confounding artefacts that are likely due to introgression among the morphological forms, the resulting phylogeny shows that the depressed shell is plesiomorphic, while the globular shell is derived. Moreover, their disparate histories suggest that the depressed shell variety survived the glaciations in pockets of alpine refugia, while the globular shell variety recolonised the alpine environment post-glacially.

Mitochondrial differentiation in a polymorphic land snail: evidence for Pleistocene survival within the boundaries of permafrost

The genetic differentiation of populations having colonized formerly unsuitable habitats after the Pleistocene glaciations depends to a great extent on the speed of expansion. Slow dispersers maintain their refugial diversity whereas fast dispersal leads to a reduction of diversity in the newly colonized areas. During the Pleistocene, almost the entire current range of the land snail Arianta arbustorum has repeatedly been covered with ice or been subjected to permafrost. Owing to the low potential for dispersal of land snails, slow (re)colonization of the wide range from southern refugia can be excluded. Alternatively, fast, passive dispersal from southern refugia or survival in and expansion from multiple refugia within the area subjected to permafrost may account for the current distribution. To distinguish between these scenarios we reconstructed a phylogeography based on the sequences of a fragment of the cytochrome oxidase I from 133 individuals collected at 45 localities and analysed the molecular variance. Seventy-five haplotypes were found that diverged on average at 7.52% of positions. This high degree of diversity suggests that A. arbustorum is an old species in which the population structure, isolation and the hermaphroditic nature have reduced the probability of lineage extinction. The genetic structure was highly significant with the highest variance partition found among regions. Geographic distance and mitochondrial differentiation were not congruent. Lineages had overlapping ranges. The clear genetic differentiation and the patchy pattern of haplotype distribution suggest that colonization of formerly unsuitable habitats was mainly achieved from multiple populations from within the permafrost area.