Nunatak survival of the high Alpine plant Eritrichium nanum (L.) Gaudin in the central Alps during the ice ages

Evolutionary history of an Alpine Archaeognath (Machilis pallida): Insights from different variant

Reconstruction of species histories is a central aspect of evolutionary biology. Patterns of genetic variation within and among populations can be leveraged to elucidate evolutionary processes and demographic histories. However, interpreting genetic signatures and unraveling the contributing processes can be challenging, in particular for non-model organisms with complex reproductive modes and genome organization. One way forward is the combined consideration of patterns revealed by different molecular markers (nuclear vs. mitochondrial) and types of variants (common vs. rare) that differ in their age, mode, and rate of evolution. Here, we applied this approach to RNAseq data generated for Machilis pallida (Archaeognatha), an Alpine jumping bristletail considered parthenogenetic and triploid. We generated de novo transcriptome and mitochondrial assemblies to obtain high-density data to investigate patterns of mitochondrial and common and rare nuclear variation in 17 M. pallida individuals sampled from all known populations. We find that the different variant types capture distinct aspects of the evolutionary history and discuss the observed patterns in the context of parthenogenesis, polyploidy, and survival during glaciation. This study highlights the potential of different variant types to gain insights into evolutionary scenarios even from challenging but often available data and the suitability of M. pallida and the genus Machilis as a study system for the evolution of sexual strategies and polyploidization during environmental change. We also emphasize the need for further research which will be stimulated and facilitated by these newly generated resources and insights.

Fine-scale genome-wide signature of Pleistocene glaciation in Thitarodes moths (Lepidoptera: Hepialidae), host of Ophiocordyceps fungus in the Hengduan Mountains

The Hengduan Mountains region is a biodiversity hotspot known for its topologically complex, deep valleys and high mountains. While landscape and glacial refugia have been evoked to explain patterns of inter-species divergence, the accumulation of intra-species (i.e. population level) genetic divergence across the mountain-valley landscape in this region has received less attention. We used genome-wide restriction site-associated DNA sequencing (RADseq) to reveal signatures of Pleistocene glaciation in populations of Thitarodes shambalaensis (Lepidoptera: Hepialidae), the host moth of parasitic Ophiocordyceps sinensis (Hypocreales: Ophiocordycipitaceae) or "caterpillar fungus" endemic to the glacier of eastern Mt. Gongga. We used moraine history along the glacier valleys to model the distribution and environmental barriers to gene flow across populations of T. shambalaensis. We found that moth populations separated by less than 10 km exhibited valley-based population genetic clustering and isolation-by-distance (IBD), while gene flow among populations was best explained by models using information about their distributions at the local last glacial maximum (LGM_L , 58 kya), not their contemporary distribution. Maximum likelihood lineage history among populations, and among subpopulations as little as 500 meters apart, recapitulated glaciation history across the landscape. We also found signals of isolated population expansion following the retreat of LGM_L glaciers. These results reveal the fine-scale, long-term historical influence of landscape and glaciation on the genetic structuring of populations of an endangered and economically important insect species. Similar mechanisms, given enough time and continued isolation, could explain the contribution of glacier refugia to the generation of species diversity among the Hengduan Mountains.

Extensive introgression at late stages of species formation: Insights from grasshopper hybrid zones

The process of species formation is characterised by the accumulation of multiple reproductive barriers. The evolution of hybrid male sterility, or Haldane's rule, typically characterises later stages of species formation, when reproductive isolation is strongest. Yet, understanding how quickly reproductive barriers evolve and their consequences for maintaining genetic boundaries between emerging species remains a challenging task because it requires studying taxa that hybridise in nature. Here, we address these questions using the meadow grasshopper Pseudochorthippus parallelus, where populations that show multiple reproductive barriers, including hybrid male sterility, hybridise in two natural hybrid zones. Using mitochondrial data, we infer that such populations have diverged some 100,000 years ago, at the beginning of the last glacial cycle in Europe. Nuclear data shows that contractions at multiple glacial refugia, and post-glacial expansions have facilitated genetic differentiation between lineages that today interact in hybrid zones. We find extensive introgression throughout the sampled species range, irrespective of current strength of reproductive isolation. Populations exhibiting hybrid male sterility in two hybrid zones show repeatable patterns of genomic differentiation, consistent with shared genomic constraints affecting ancestral divergence or with the role of those regions in reproductive isolation. Together, our results suggest that reproductive barriers that characterise late stages of species formation can evolve relatively quickly, particularly when associated with strong demographic changes. Moreover, we show that such barriers persist in the face of extensive gene flow, allowing future studies to identify associated genomic regions.

Climate-Change Impacts on the Southernmost Mediterranean Arctic-Alpine Plant Populations

Human-induced climate- and land-use change have been affecting biogeographical and biodiversity patterns for the past two centuries all over the globe, resulting in increased extinction and biotic homogenization rates. High mountain ecosystems are more sensitive to these changes, which have led to physiological and phenological shifts, as well as to ecosystem processes’ deformation. Glacial relicts, such as arctic-alpine taxa, are sensitive indicators of the effects of global warming and their rear-edge populations could include warm-adapted genotypes that might prove—conservation-wise—useful in an era of unprecedented climate regimes. Despite the ongoing thermophilization in European and Mediterranean summits, it still remains unknown how past and future climate-change might affect the distributional patterns of the glacial relict, arctic-alpine taxa occurring in Greece, their European southernmost distributional limit. Using species distribution models, we investigated the impacts of past and future climate changes on the arctic-alpine taxa occurring in Greece and identified the areas comprising arctic-alpine biodiversity hotspots in Greece. Most of these species will be faced with severe range reductions in the near future, despite their innate resilience to a multitude of threats, while the species richness hotspots will experience both altitudinal and latitudinal shifts. Being long-lived perennials means that there might be an extinction-debt present in these taxa, and a prolonged stability phase could be masking the deleterious effects of climate change on them. Several ex situ conservation measures (e.g., seed collection, population augmentation) should be taken to preserve the southernmost populations of these rare arctic-alpine taxa and a better understanding of their population genetics is urgently needed.

Patterns of protist distribution and diversification in alpine lakes across Europe

Biogeography in Europe is known to be crucially influenced by the large mountain ranges serving as biogeographical islands for cold-adapted taxa and geographical barriers for warm-adapted taxa. While biogeographical patterns are well-known for plants and animals in Europe, we here investigated diversity and distribution patterns of protist freshwater communities on a European scale (256 lakes) in the light of the well-studied post-glacial distribution patterns of macroorganisms. Thus, our study compared 43 alpine protist communities of lakes located in the Alps, Carpathians, Pyrenees, and the Sierra Nevada with that of surrounding lowland lakes. We verified altitudinal diversity gradients of freshwater protists with decreasing richness and diversity across altitudes similar to those observed for plants and animals. Alpine specialists and generalists could be identified differing significantly in richness and diversity, but hardly in occurrence and proportions of major taxonomic groups. High proportions of region-specific alpine specialists indicate an increased occurrence of distinct lineages within each mountain range and thus, suggested either separated glacial refugia or post-glacial diversification within mountain ranges. However, a few alpine specialists were shared between mountain ranges suggesting a post-glacial recolonization from a common lowland pool. Our results identified generalists with wide distribution ranges and putatively wide tolerance ranges toward environmental conditions as main drivers of protist diversification (specification) in alpine lakes, while there was hardly any diversification in alpine specialists.

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.

Evidence for Glacial Refugia of the Forest Understorey Species Helleborus niger (Ranunculaceae) in the Southern as Well as in the Northern Limestone Alps

Glacial refugia of alpine and subnival biota have been intensively studied in the European Alps but the fate of forests and their understory species in that area remains largely unclear. In order to fill this gap, we aimed at disentangling the spatiotemporal diversification of disjunctly distributed black hellebore Helleborus niger (Ranunculaceae). We applied a set of phylogeographic analyses based on restriction-site associated DNA sequencing (RADseq) data and plastid DNA sequences to a range-wide sampling of populations. These analyses were supplemented with species distribution models generated for the present and the Last Glacial Maximum (LGM). We used exploratory analyses to delimit genomically coherent groups and then employed demographic modeling to reconstruct the history of these groups. We uncovered a deep split between two major genetic groups with western and eastern distribution within the Southern Limestone Alps, likely reflecting divergent evolution since the mid-Pleistocene in two glacial refugia situated along the unglaciated southern margin of the Alps. Long-term presence in the Southern Limestone Alps is also supported by high numbers of private alleles, elevated levels of nucleotide diversity and the species' modeled distribution at the LGM. The deep genetic divergence, however, is not reflected in leaf shape variation, suggesting that the morphological discrimination of genetically divergent entities within H. niger is questionable. At a shallower level, populations from the Northern Limestone Alps are differentiated from those in the Southern Limestone Alps in both RADseq and plastid DNA data sets, reflecting the North-South disjunction within the Eastern Alps. The underlying split was dated to ca. 0.1 mya, which is well before the LGM. In the same line, explicit tests of demographic models consistently rejected the hypothesis that the partial distribution area in the Northern Limestone Alps is the result of postglacial colonization. Taken together, our results strongly support that forest understory species such as H. niger have survived the LGM in refugia situated along the southern, but also along the northern or northeastern periphery of the Alps. Being a slow migrator, the species has likely survived repeated glacial-interglacial circles in distributional stasis while the composition of the tree canopy changed in the meanwhile.

Drainage basins serve as multiple glacial refugia for alpine habitats in the Sierra Nevada Mountains, California

The evolutionary histories of alpine species are often directly associated with responses to glaciation. Deep divergence among populations and complex patterns of genetic variation have been inferred as consequences of persistence within glacier boundaries (i.e., on nunataks), while shallow divergence and limited genetic variation are assumed to result from expansion from large refugia at the edge of ice shields (i.e., massifs de refuge). However, for some species, dependence on specific microhabitats could profoundly influence their spatial and demographic response to glaciation, and such a simple dichotomy may obscure the localization of actual refugia. In this study, we use the Nebria ingens complex (Coleoptera: Carabidae), a water-affiliated ground beetle lineage, to test how drainage basins are linked to their observed population structure. By analysing mitochondrial COI gene sequences and genome-wide single nucleotide polymorphisms, we find that the major drainage systems of the Sierra Nevada Mountains in California best explain the population structure of the N. ingens complex. In addition, we find that an intermediate morphotype within the N. ingens complex is the product of historical hybridization of N. riversi and N. ingens in the San Joaquin basin during glaciation. This study highlights the importance of considering ecological preferences in how species respond to climate fluctuations and provides an explanation for discordances that are often observed in comparative phylogeographical studies.

Climate Change and Alpine Screes: No Future for Glacial Relict Papaver occidentale (Papaveraceae) in Western Prealps

Glacial relicts, especially those with very narrow habitat requirements, are particularly affected by global warming. We considered Papaver occidentale, a glacial relict endemic to the Western Prealps, belonging to the alpine poppy complex (P. alpinum aggr.), as a model taxon to study the actual status and potential future distribution of species restricted to particular microrefugia. For this study, all known localities were visited, each population was georeferenced and the number of individuals was estimated. Species Distribution Modelling (SDM) was used to evaluate the present and future potential distribution range and habitat suitability, taking into account the specificity of its habitat (calcareous screes). According to our study, there are globally 19 natural populations of P. occidentale, and a total of about 30,000 individuals. The taxon is a highly specialized alpine plant growing in the majority of natural sites between 1900 and 2100 m a.s.l. on north-facing screes. Predictions for the end of the 21st century indicate that a suitable area will significantly decrease (0–30% remaining). Under the most severe climatic scenarios (RCP 8.5), the species risks complete extinction. The long-term in situ conservation of P. occidentale, and all other taxa of the P. alpinum complex, is unlikely to be achieved without slowing global climate change. More generally, our fine-scale study shows that local environmental buffering of large-scale climate change in high-mountain flora may be very limited in specialised taxa of patchy environments such as screes.

Occurrence of apomictic conspecifics and ecological preferences rather than colonization history govern the geographic distribution of sexual Potentilla puberula

The geographic distribution of sexual-apomictic taxa (i.e., comprising individuals usually reproducing either sexually or asexually via seeds) is traditionally thought to be driven by their ecological preferences and colonization histories. Where sexuals and apomicts get into contact with each other, competitive and reproductive interactions can interfere with these factors, an aspect which hitherto received little attention in biogeographic studies. We disentangled and quantified the relative effects of the three factors on the distribution of tetraploid sexuals in Potentilla puberula in a latitudinal transect through the Eastern European Alps, in which they are codistributed with penta-, hepta-, and octoploid apomictic conspecifics. Effects were explored by means of binomial generalized linear regression models combining a single with a multiple predictor approach. Postglacial colonization history was inferred from population genetic variation (AFLPs and cpDNA) and quantified using a cost distance metric. The study was based on 235 populations, which were purely sexual, purely apomictic, or of mixed reproductive mode. The occurrence of apomicts explained most of the variation in the distribution of sexuals (31%). Specifically, the presence of sexual tetraploids was negatively related to the presence of each of the three apomictic cytotypes. Effects of ecological preferences were substantial too (7% and 12% of the total variation explained by ecological preferences alone, or jointly with apomicts' occurrence, respectively). In contrast, colonization history had negligible effects on the occurrence of sexuals. Taken together, our results highlight the potentially high impact of reproductive interactions on the geographic distribution of sexual and apomictic conspecifics and that resultant mutual exclusion interrelates to ecological differentiation, a situation potentially promoting their local coexistence.

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.

Is the incidence of survival in interior Pleistocene refugia (nunataks) underestimated? Phylogeography of the high mountain plant Androsace alpina (Primulaceae) in the European Alps revisited

Temperate mountain ranges such as the European Alps have been strongly affected by the Pleistocene glaciations. Glacial advances forced biota into refugia, which were situated either at the periphery of mountain ranges or in their interior. Whereas in the Alps peripheral refugia have been repeatedly and congruently identified, support for the latter scenario, termed "nunatak hypothesis," is still limited and no general pattern is recognizable yet. Here, we test the hypothesis of nunatak survival for species growing in the high alpine to subnival zones on siliceous substrate using the cushion plant Androsace alpina (Primulaceae), endemic to the European Alps, as our model species. To this end, we analyzed AFLP and plastid DNA sequence data obtained from a dense and range-wide sampling. Both AFLPs and plastid sequence data identified the southwestern-most population as the most divergent one. AFLP data did not allow for discrimination of interior and peripheral populations, but rather identified two to three longitudinally separated major gene pools. In contrast, in the eastern half of the Alps several plastid haplotypes of regional or local distribution in interior ranges-the Alpine periphery mostly harbored a widespread haplotype-were indicative for the presence of interior refugia. Together with evidence from other Alpine plant species, this study shows that in the eastern Alps silicicolous species of open habitats in the alpine and subnival zone survived, also or exclusively so, in interior refugia. As the corresponding genetic structure may be lost in mostly nuclear-derived, rapidly homogenizing marker systems such as AFLPs or RAD sequencing tags, markers not prone to homogenization, as is the case for plastid sequences (Sanger-sequenced or extracted from an NGS data set) will continue to be important for detecting older, biogeographically relevant patterns.

Glacial refugia and postglacial expansion of the alpine-prealpine plant species Polygala chamaebuxus

The shrubby milkwort (Polygala chamaebuxus L.) is widely distributed in the Alps, but occurs also in the lower mountain ranges of Central Europe such as the Franconian Jura or the Bohemian uplands. Populations in these regions may either originate from glacial survival or from postglacial recolonization. In this study, we analyzed 30 populations of P. chamaebuxus from the whole distribution range using AFLP (Amplified Fragment Length Polymorphism) analysis to identify glacial refugia and to illuminate the origin of P. chamaebuxus in the lower mountain ranges of Central Europe. Genetic variation and the number of rare fragments within populations were highest in populations from the central part of the distribution range, especially in the Southern Alps (from the Tessin Alps and the Prealps of Lugano to the Triglav Massiv) and in the middle part of the northern Alps. These regions may have served, in accordance with previous studies, as long-term refugia for the glacial survival of the species. The geographic pattern of genetic variation, as revealed by analysis of molecular variance, Bayesian cluster analysis and a PopGraph genetic network was, however, only weak. Instead of postglacial recolonization from only few long-term refugia, which would have resulted in deeper genetic splits within the data set, broad waves of postglacial expansion from several short-term isolated populations in the center to the actual periphery of the distribution range seem to be the scenario explaining the observed pattern of genetic variation most likely. The populations from the lower mountain ranges in Central Europe were more closely related to the populations from the southwestern and northern than from the nearby eastern Alps. Although glacial survival in the Bohemian uplands cannot fully be excluded, P. chamaebuxus seems to have immigrated postglacially from the southwestern or central-northern parts of the Alps into these regions during the expansion of the pine forests in the early Holocene.

On the origins of Balkan endemics: the complex evolutionary history of the Cyanus napulifer group (Asteraceae)

BACKGROUND AND AIMS

The Balkan Peninsula is one of the most important centres of plant diversity in Europe. Here we aim to fill the gap in the current knowledge of the evolutionary processes and factors modelling this astonishing biological richness by applying multiple approaches to the Cyanus napulifer group.

METHODS

To reconstruct the mode of diversification within the C. napulifer group and to uncover its relationships with potential relatives with x = 10 from Europe and Northern Africa, we examined variation in genetic markers (amplified fragment length polymorphisms [AFLPs]; 460 individuals), relative DNA content (4',6-diamidino-2-phenylindole [DAPI] flow cytometry, 330 individuals) and morphology (multivariate morphometrics, 40 morphological characters, 710 individuals). To elucidate its evolutionary history, we analysed chloroplast DNA (cpDNA) sequences of the genus Cyanus deposited in the GenBank database.

KEY RESULTS

The AFLPs revealed a suite of closely related entities with variable levels of differentiation. The C. napulifer group formed a genetically well-defined unit. Samples outside the group formed strongly diversified and mostly species-specific genetic lineages with no further geographical patterns, often characterized also by a different DNA content. AFLP analysis of the C. napulifer group revealed extensive radiation and split it into nine allopatric (sub)lineages with varying degrees of congruence among genetic, DNA-content and morphological patterns. Genetic admixture was usually detected in contact zones between genetic lineages. Plastid data indicated extensive maintenance of ancestral variation across Cyanus perennials.

CONCLUSION

The C. napulifer group is an example of a rapidly and recently diversified plant group whose genetic lineages have evolved in spatio-temporal isolation on the topographically complex Balkan Peninsula. Adaptive radiation, accompanied in some cases by long-term isolation and hybridization, has contributed to the formation of this species complex and its mosaic pattern.

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.

A glacial survivor of the alpine Mediterranean region: phylogenetic and phylogeographic insights into Silene ciliata Pourr. (Caryophyllaceae)

Silene ciliata Pourr. (Caryophyllaceae) is a species with a highly disjunct distribution which inhabits the alpine mountains of the Mediterranean Basin. We investigated the phylogeny and phylogeography of the species to (a) clarify the long-suggested division of S. ciliata into two subspecies, (b) evaluate its phylogenetic origin and (c) assess whether the species' diversification patterns were affected by the Mediterranean relief. For this purpose, we collected DNA from 25 populations of the species that inhabit the mountains of Portugal, Spain, France, Italy, former Yugoslav Republic of Macedonia, Bulgaria and Greece and studied the plastid regions rbcL, rps16 and trnL. Major intraspecific variation was supported by all analyses, while the possibility of the existence of more varieties or subspecies was not favoured. Plastid DNA (cpDNA) evidence was in accordance with the division of S. ciliata into the two subspecies, one spreading west (Iberian Peninsula and Central Massif) and the other east of the Alps region (Italian and Balkan Peninsula). This study proposes that the species' geographically disconnected distribution has probably derived from vicariance processes and from the Alps acting as a barrier to the species' dispersal. The monophyletic origin of the species is highly supported. cpDNA patterns were shown independent of the chromosome evolution in the populations and could have resulted from a combination of geographic factors providing links and barriers, climatic adversities and evolutionary processes that took place during Quaternary glaciations.

Postglacial recolonization in a cold climate specialist in western Europe: patterns of genetic diversity in the adder (Vipera berus) support the central-marginal hypothesis

Understanding the impact of postglacial recolonization on genetic diversity is essential in explaining current patterns of genetic variation. The central-marginal hypothesis (CMH) predicts a reduction in genetic diversity from the core of the distribution to peripheral populations, as well as reduced connectivity between peripheral populations. While the CMH has received considerable empirical support, its broad applicability is still debated and alternative hypotheses predict different spatial patterns of genetic diversity. Using microsatellite markers, we analysed the genetic diversity of the adder (Vipera berus) in western Europe to reconstruct postglacial recolonization. Approximate Bayesian Computation (ABC) analyses suggested a postglacial recolonization from two routes: a western route from the Atlantic Coast up to Belgium and a central route from the Massif Central to the Alps. This cold-adapted species likely used two isolated glacial refugia in southern France, in permafrost-free areas during the last glacial maximum. Adder populations further from putative glacial refugia had lower genetic diversity and reduced connectivity; therefore, our results support the predictions of the CMH. Our study also illustrates the utility of highly variable nuclear markers, such as microsatellites, and ABC to test competing recolonization hypotheses.

Out of the alps: the biogeography of a disjunctly distributed mountain butterfly, the almond-eyed ringlet Erebia alberganus (Lepidoptera, Satyrinae)

Many studies on the biogeography of thermophilic and arctic-alpine species were performed during the past. Only little is known about species with intermediate characteristics. We analyzed the molecular biogeography of the butterfly Erebia alberganus (30 populations, representing 1106 individuals), sampled over the Alps, Apennines (Italy), and the Stara Planina (Bulgaria) using allozyme electrophoresis (17 loci). Genetic analyses revealed 3 major splits, with the strongest between the Stara Planina populations and all other populations, and a weaker split between the Alps and the Apennines. Individuals from the Apennines were genetically nested within the Alps group. The Alps cluster was segregated into 3 groups: the Southwestern, Western/Central, and Eastern Alps. The genetic diversities were highest for the Alps populations and significantly lower in the 2 isolates (Apennines, Stara Planina). The remarkable genetic split between Stara Planina and all other populations and the genetic distinctiveness of the former cluster might be interpreted as an ancient colonization event of this Balkan mountain range. The Apennines populations derive from a more recent expansion out of the Southwestern Alps. After surviving the Würm ice age most probably in the central Apennines, accompanied by genetic modification of some of these populations, northward expansion might have started from the western parts of the central Apennines reaching the northern Apennines during the early postglacial. The subtle genetic differentiation found among the Alps populations probably reflects 3 geographically disjunct Würm glacial centers located at the western slopes of the Southwestern Alps, at the southern slopes of the Central Alps, and in the Southeastern Alps.

The intraspecific genetic variability of siliceous and calcareous Gentiana species is shaped by contrasting demographic and re-colonization processes

The Ciminalis section of Gentiana comprises seven species, two of them growing on siliceous substrates (G. alpina and G. acaulis), the other ones being calcareous taxa (G. clusii, G. angustifolia, G. ligustica, G. occidentalis and G. dinarica). A total of 515 individuals from 183 populations over the entire Ciminalis distribution range was analyzed using four chloroplast loci (trnH-psbA, matK, rpoB and rpoC1) and the nuclear ribosomal marker ITS2. The siliceous species display only two chloroplast haplotypes each and are both characterized by patterns of range expansions all over the Alps. Conversely, the calcareous species are on average more diverse (two to 13 haplotypes per species) with strong patterns of local structuring. We suggest that the occurrence of many calcareous refugia at the periphery of the Alps must have led to local adaptation and morphological diversification, and helped preserving intraspecific diversities during the last glaciations for the associated taxa. ITS2 was more efficient in delineating species boundaries than the chloroplast markers for which several haplotypes are shared among species. This might be either due to chloroplast capture among species and/or to recent divergence. Species adapted to the same substrate are generally only distantly related when they co-occur in the same place. For both types of markers, G. clusii is found genetically distant from all other species.

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.

Population genetic analysis of Lobelia rhynchopetalum Hemsl. (Campanulaceae) using DNA sequences from ITS and eight chloroplast DNA regions

DNA sequence data from the internal transcribed spacer of nuclear ribosomal DNA and eight chloroplast DNA regions were used to investigate haplotypic variation and population genetic structure of the Afroalpine giant lobelia, Lobelia rhynchopetalum. The study was based on eight populations sampled from two mountain systems in Ethiopia. A total of 20 variable sites were obtained, which resulted in 13 unique haplotypes and an overall nucleotide diversity (ND) of 0.281 ± 0.15 and gene diversity (GD) of 0.85 ± 0.04. Analysis of molecular variance (AMOVA) revealed a highly significant variation (P < 0.001) among populations (F_ST), and phylogenetic analysis revealed that populations from the two mountain systems formed their own distinct clade with >90% bootstrap support. Each population should be regarded as a significant unit for conservation of this species. The primers designed for this study can be applied to any Lobelia and other closely related species for population genetics and phylogenetic studies.

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.

Isolated history of the coastal plant Lathyrus japonicus (Leguminosae) in Lake Biwa, an ancient freshwater lake

BACKGROUND AND AIMS

Lake Biwa is one of the world's few ancient lakes. Formed ∼4 million years ago, the lake harbours many coastal species that commonly inhabit seashores. The beach pea Lathyrus japonicus is a typical coastal species of this freshwater lake, but its inland populations are faced with the threat of extinction. Here, we investigated the phylogeographical and population structures of both inland and coastal populations of L. japonicus. We also elucidated the historical isolation of the Lake Biwa population.

METHODOLOGY

In total, 520 individuals from 50 L. japonicus populations were sampled throughout the species distribution in Japan. Chloroplast haplotyping using intergenic spacers psbA-trnH and atpI-atpH was performed to investigate the phylogeographical structure as well as the genetic diversity of L. japonicus. Six nuclear microsatellite markers were also used to analyse the population structure.

PRINCIPAL RESULTS

Population structure analyses of chloroplast DNA (cpDNA) and nuclear DNA (nDNA) identified inland and coastal groups. Based on the genetic differentiation, inland populations exhibited a single cpDNA haplotype and significantly lower values of H(S), AR and F(IS) than coastal populations. In addition to the presence of a bottleneck, the lack of gene flow among inland populations was supported by estimates of recent migration rates between subpopulations.

CONCLUSIONS

Our data revealed that inland populations have been isolated in Lake Biwa as 'landlocked' populations since the predecessor lake was isolated from sea. This was also seen in a previous study of Calystegia soldanella. However, the high genetic differentiation, accompanied by a lack of gene flow among the Lake Biwa populations (according to the BAYESASS(+) analysis), contradicts the results with C. soldanella. We conclude that because of the presence of a bottleneck and low genetic diversity of the inland populations, self-sustaining population persistence may be difficult in the future. Conservation strategies must consider the genetic properties of such isolated populations.

Genetic population structure of the alpine species Rhododendron pseudochrysanthum sensu lato (Ericaceae) inferred from chloroplast and nuclear DNA

Background

A complex of incipient species with different degrees of morphological or ecological differentiation provides an ideal model for studying species divergence. We examined the phylogeography and the evolutionary history of the Rhododendron pseudochrysanthum s. l.

Results

Systematic inconsistency was detected between gene genealogies of the cpDNA and nrDNA. Rooted at R. hyperythrum and R. formosana, both trees lacked reciprocal monophyly for all members of the complex. For R. pseudochrysanthum s.l., the spatial distribution of the cpDNA had a noteworthy pattern showing high genetic differentiation (F_ST = 0.56-0.72) between populations in the Yushan Mountain Range and populations of the other mountain ranges.

Conclusion

Both incomplete lineage sorting and interspecific hybridization/introgression may have contributed to the lack of monophyly among R. hyperythrum, R. formosana and R. pseudochrysanthum s.l. Independent colonizations, plus low capabilities of seed dispersal in current environments, may have resulted in the genetic differentiation between populations of different mountain ranges. At the population level, the populations of Central, and Sheishan Mountains may have undergone postglacial demographic expansion, while populations of the Yushan Mountain Range are likely to have remained stable ever since the colonization. In contrast, the single population of the Alishan Mountain Range with a fixed cpDNA haplotype may have experienced bottleneck/founder's events.

Phylogeography of Japanese horse chestnut (Aesculus turbinata) in the Japanese Archipelago based on chloroplast DNA haplotypes

Japanese horse chestnut (Aesculus turbinata: Hippocastanaceae) is one of the typical woody plants that grow in temperate riparian forests in the Japanese Archipelago. To analyze the phylogeography of this plant in the Japanese Archipelago, we determined cpDNA haplotypes for 337 samples from 55 populations covering the entire distribution range. Based on 1,313 bp of two spacers, we determined ten haplotypes that are distinguished from adjacent haplotypes by one or two steps. Most of the populations had a single haplotype, suggesting low diversity. Spatial analysis of molecular variance suggested three obvious phylogeographic structures in western Japan, where Japanese horse chestnut is scattered and isolated in mountainous areas. Conversely, no clear phylogeographic structure was observed from the northern to the southern limit of this species, including eastern Japan, where this plant is more common. Rare and private haplotypes were also found in southwestern Japan, where Japanese horse chestnuts are distributed sparsely. These findings imply that western Japan might have maintained a relatively large habitat for A. turbinata during the Quaternary climatic oscillations, while northerly regions could not.

Cross-scale analysis of the region effect on vascular plant species diversity in southern and northern European mountain ranges

BACKGROUND

The divergent glacial histories of southern and northern Europe affect present-day species diversity at coarse-grained scales in these two regions, but do these effects also penetrate to the more fine-grained scales of local communities?

METHODOLOGY/PRINCIPAL FINDINGS

We carried out a cross-scale analysis to address this question for vascular plants in two mountain regions, the Alps in southern Europe and the Scandes in northern Europe, using environmentally paired vegetation plots in the two regions (n = 403 in each region) to quantify four diversity components: (i) total number of species occurring in a region (total γ-diversity), (ii) number of species that could occur in a target plot after environmental filtering (habitat-specific γ-diversity), (iii) pair-wise species compositional turnover between plots (plot-to-plot β-diversity) and (iv) number of species present per plot (plot α-diversity). We found strong region effects on total γ-diversity, habitat-specific γ-diversity and plot-to-plot β-diversity, with a greater diversity in the Alps even towards distances smaller than 50 m between plots. In contrast, there was a slightly greater plot α-diversity in the Scandes, but with a tendency towards contrasting region effects on high and low soil-acidity plots.

CONCLUSIONS/SIGNIFICANCE

We conclude that there are strong regional differences between coarse-grained (landscape- to regional-scale) diversity components of the flora in the Alps and the Scandes mountain ranges, but that these differences do not necessarily penetrate to the finest-grained (plot-scale) diversity component, at least not on acidic soils. Our findings are consistent with the contrasting regional Quaternary histories, but we also consider alternative explanatory models. Notably, ecological sorting and habitat connectivity may play a role in the unexpected limited or reversed region effect on plot α-diversity, and may also affect the larger-scale diversity components. For instance, plot connectivity and/or selection for high dispersal ability may increase plot α-diversity and compensate for low total γ-diversity.

Disjunct distribution of chloroplast DNA haplotypes in the understory perennial Veratrum album ssp. oxysepalum (Melanthiaceae) in Japan as a result of ancient introgression

• The Quaternary climatic changes resulted in range shifts of species, providing chances for hybridization. However, the genetic signatures of such ancient introgression have rarely been reported. To investigate such signatures, we performed a phylogeographical study on the perennial plant Veratrum album ssp. oxysepalum, which may have hybridized long ago with another congeneric species, V. stamineum. • Sequence variations in chloroplast DNA (cpDNA) were examined in 43 populations in Japan and adjacent areas. Phylogenetic analyses of different cpDNA haplotypes were conducted on the basis of cpDNA and nuclear ribosomal internal transcribed spacer (nrITS) variations. • In the Japanese archipelago, two major groups of haplotypes were detected, one of which was distributed in a disjunct pattern. The major haplotype, occupying the central part of the species' distribution, formed a monophyletic group with V. stamineum in phylogenetic trees on the basis of cpDNA variation, although the two species did not form a monophyletic group in phylogenetic trees on the basis of nrITS variation. • Historical hybridization between V. album ssp. oxysepalum and V. stamineum in refugia during the Quaternary climatic oscillations, and the resulting chloroplast capture of V. stamineum by V. album ssp. oxysepalum, are most probably responsible for the disjunct distribution of cpDNA in V. album ssp. oxysepalum.

Whole genome sequencing of enriched chloroplast DNA using the Illumina GAII platform

BACKGROUND

Complete chloroplast genome sequences provide a valuable source of molecular markers for studies in molecular ecology and evolution of plants. To obtain complete genome sequences, recent studies have made use of the polymerase chain reaction to amplify overlapping fragments from conserved gene loci. However, this approach is time consuming and can be more difficult to implement where gene organisation differs among plants. An alternative approach is to first isolate chloroplasts and then use the capacity of high-throughput sequencing to obtain complete genome sequences. We report our findings from studies of the latter approach, which used a simple chloroplast isolation procedure, multiply-primed rolling circle amplification of chloroplast DNA, Illumina Genome Analyzer II sequencing, and de novo assembly of paired-end sequence reads.

RESULTS

A modified rapid chloroplast isolation protocol was used to obtain plant DNA that was enriched for chloroplast DNA, but nevertheless contained nuclear and mitochondrial DNA. Multiply-primed rolling circle amplification of this mixed template produced sufficient quantities of chloroplast DNA, even when the amount of starting material was small, and improved the template quality for Illumina Genome Analyzer II (hereafter Illumina GAII) sequencing. We demonstrate, using independent samples of karaka (Corynocarpus laevigatus), that there is high fidelity in the sequence obtained from this template. Although less than 20% of our sequenced reads could be mapped to chloroplast genome, it was relatively easy to assemble complete chloroplast genome sequences from the mixture of nuclear, mitochondrial and chloroplast reads.

CONCLUSIONS

We report successful whole genome sequencing of chloroplast DNA from karaka, obtained efficiently and with high fidelity.

Genetic diversity, structure, and demographic change in tanoak, Lithocarpus densiflorus (Fagaceae), the most susceptible species to sudden oak death in California

Knowledge of population genetic structure of tanoak (Lithocarpus densiflorus) is of interest to pathologists seeking natural variation in resistance to sudden oak death disease, to resource managers who need indications of conservation priorities in this species now threatened by the introduced pathogen (Phytophthora ramorum), and to biologists with interests in demographic processes that have shaped plant populations. We investigated population genetic structure using nuclear and chloroplast DNA (cpDNA) and inferred the effects of past population demographic processes and contemporary gene flow. Our cpDNA results revealed a strong pattern of differentiation of four regional groups (coastal California, southern Oregon, Klamath mountains, and Sierra Nevada). The chloroplast haplotype phylogeny suggests relatively deep divergence of Sierra Nevada and Klamath populations from those of coastal California and southern Oregon. A widespread coastal California haplotype may have resulted from multiple refugial sites during the Last Glacial Maximum or from rapid recolonization from few refugia. Analysis of nuclear microsatellites suggests two major groups: (1) central coastal California and (2) Sierra Nevada/Klamath/southern Oregon and an area of admixture in north coastal California. The low level of nuclear differentiation is likely to be due to pollen gene flow among populations during postglacial range expansion.

Bayesian hypothesis testing supports long-distance Pleistocene migrations in a European high mountain plant (Androsace vitaliana, Primulaceae)

Colonization of the south-western European mountain ranges is suggested to have predominantly progressed from the Iberian Peninsula eastwards, but this hypothesis has never been tested in a statistical framework. Here, we test this hypothesis using Androsace vitaliana, a high elevation species with eight mostly allopatric subspecies, which is widely but disjunctly distributed across all major south-western European mountain ranges. To this end, we use plastid and nuclear sequence data as well as fingerprint (amplified fragment length polymorphisms) data and employ Bayesian methods, which allow co-estimation of genealogy and divergence times using explicit demographic models, as well as hypothesis testing via Bayes factors. Irrespective of the ambiguity concerning where A. vitaliana started to diversify -- both the Alps and the mountain ranges of the Iberian Peninsula outside the Pyrenees were possible -- colonization routes were not simply unidirectional, but involved Pleistocene connections between the Alps and mountain ranges of the Iberian Peninsula bypassing the interjacent Pyrenees via long-distance dispersal. In contrast, the species' post-glacial history is shaped by regional gene pool homogenization resulting in the genetic pattern showing good congruence with geographical proximity in agreement with a vicariance model, but only partly supporting current taxonomy.

Haplotype richness in refugial areas: phylogeographical structure of Saxifraga callosa

This paper illustrates the phylogeographical structure of Saxifraga callosa in order to describe its genetic richness in refugial areas and to reconstruct its glacial history. S. callosa is a species spread throughout south-east France and Italy with a high distribution in the Maritime Alps. Four chloroplast microsatellite and AFLP markers were analyzed in populations of S. callosa. The size variants of all tested loci amount to 11 different haplotypes. Intrapopulational haplotype variation was found in two of the populations analyzed: on the Mt. Toraggio in the Maritime Alps, and in the Apuan Alps. On the other hand, no intrapopulational variation was found in 25 populations, most of which were sampled from isolated areas. Analysis of the haplotype distribution showed that population subdivision across all populations was high (G (ST) = 0.899). Moreover, its genetic structure was studied using AMOVA and STRUCTURE analysis. The study legitimated inferred conclusions about the phylogeographical structure of the species and identified centers of diversity. Considerations concerning genetic structure and divergence among three major clades (Maritime Alps, Apuan Alps and Apennines), the patchy distribution of haplotypes, and the high number of private haplotypes support the proposal that S. callosa survived in some refugia within the Italian Peninsula refugium, and that mainly northern populations of refugia were involved in postglacial recolonization.

Repeated range expansion and glacial endurance of Potentilla glabra (Rosaceae) in the Qinghai-Tibetan plateau

To date, little is still known about how alpine species occurring in the Qinghai-Tibetan Plateau (QTP) responded to past climatic oscillations. Here, by using variations of the chloroplast trnT-L, we examined the genetic distribution pattern of 101 individuals of Potentilla glabra, comprising both the interior QTP and the plateau edge. Phylogenetic and network analyses of 31 recovered haplotypes identified three tentative clades (A, B and C). Analysis of molecular variance (amova) revealed that most of the genetic variability was found within populations (0.693), while differentiations between populations were obviously distinct (F(st)= 0.307). Two independent range expansions within clades A and B occurring at approximately 316 and 201 thousand years ago (kya) were recovered from the hierarchical mismatch analysis, and these two expansions were also confirmed by Fu's F(S) values and 'g' tests. However, distant distributions of clade C and private haplotypes from clades A and B suggest that they had survived the Last Glacial Maximum (LGM) and previous glaciers in situ since their origins. Our findings based on available limited samples support that multiple refugia of a few cold-enduring species had been maintained in the QTP platform during LGM and/or previous glacial stages.

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.