Historical divergence vs. contemporary gene flow: evolutionary history of the calcicoleRanunculus alpestrisgroup (Ranunculaceae) in the European Alps and the Carpathians

Perhaps there were northern refugia in LGM? The phylogeographic structure of the thermophilic tree Carpinus betulus in the Carpathian region

Carpinus betulus L., the hornbeam, is a component of lowland and highland forests in Europe. By examining the postglacial migratory history of thermophilic tree species, the study aimed to unravel their putative glacial microrefugia in the Carpathian region. The present study points to the two distinct genetic AFLP groups of C. betulus in the Carpathian region that represent different genetic lineages based on Bayesian analysis. They differed in Nei's gene diversity index h, and the analysis of molecular variance AMOVA showed a percentage variation of the populations between the groups of 13.74 %. Principal coordinate analysis (PCoA) of 368 AFLP tree samples confirmed the presence of two genetic groups. Ninety-five populations underwent principal component analysis (PCA) to show the main correlations between genetic diversity indices and bioclimatic/climate variables (WorldClim and Carpatclim). The generalized logistic model (GLM) showed the significance of Nei's genetic index h in delimiting genetic groups. The results of population-genetic and multivariate analyses determined that the two genetic groups nowadays are spatially diffused and do not show a clear geographic pattern, pointing to a genetic melting pot. We found ecological links between genetic diversity and bioclimatic characteristics, especially the precipitation in the coldest quarter - Bio19. The refugial Maxent model indicates a significant contribution of the Bio7 variable (both linked with a continental type of climate) to the occurrence of the species during the LGM in Europe. We suggest the relict character of hornbeam populations in a specific climatic-terrain niche in the northern part of the Carpathian Basin.

Evolution of the Sabulina verna group (Caryophyllaceae) in Europe: A deep split, followed by secondary contacts, multiple allopolyploidization and colonization of challenging substrates

One of the major goals of contemporary evolutionary biology is to elucidate the relative roles of allopatric and ecological differentiation and polyploidy in speciation. In this study, we address the taxonomically intricate Sabulina verna group, which has a disjunct Arctic-alpine postglacial range in Europe and occupies a broad range of ecological niches, including substrates toxic to plants. Using genome-wide ddRAD sequencing combined with morphometric analyses based on extensive sampling of 111 natural populations, we aimed to disentangle internal evolutionary relationships and examine their correspondence with the pronounced edaphic and ploidy diversity within the group. We identified two spatially distinct groups of diploids: a widespread Arctic-alpine group and a spatially restricted yet diverse Balkan group. Most tetraploids exhibited a considerably admixed ancestry derived from both these groups, suggesting their allopolyploid origin. Four genetic clusters in congruence with geography and mostly supported by morphological traits were recognized in the diploid Arctic-alpine group. Tetraploids are split into two distinct and geographically vicariant groups, indicating their repeated polytopic origin. Furthermore, our results also revealed at least five-fold parallel colonization of toxic substrates (serpentine and metalliferous), altogether demonstrating a complex interaction between geography, challenging substrates and polyploidy in the evolution of the group. Finally, we propose a new taxonomic treatment of this complex.

Parallel adaptation to lower altitudes is associated with enhanced plasticity in Heliosperma pusillum (Caryophyllaceae)

High levels of phenotypic plasticity are thought to be inherently costly in stable or extreme environments, but enhanced plasticity may evolve as a response to new environments and foster novel phenotypes. Heliosperma pusillum forms glabrous alpine and pubescent montane ecotypes that diverged recurrently and polytopically (parallel evolution) and can serve as evolutionary replicates. The specific alpine and montane localities are characterized by distinct temperature conditions, available moisture, and light. Noteworthy, the ecotypes show a home-site fitness advantage in reciprocal transplantations. To disentangle the relative contribution of constitutive versus plastic gene expression to altitudinal divergence, we analyze the transcriptomic profiles of two parallely evolved ecotype pairs, grown in reciprocal transplantations at native altitudinal sites. In this incipient stage of divergence, only a minor proportion of genes appear constitutively differentially expressed between the ecotypes in both pairs, regardless of the growing environment. Both derived, montane populations bear comparatively higher plasticity of gene expression than the alpine populations. Genes that change expression plastically or constitutively underlie similar ecologically relevant pathways, related to response to drought and trichome formation. Other relevant processes, such as photosynthesis, rely mainly on plastic changes. The enhanced plasticity consistently observed in the montane ecotype likely evolved as a response to the newly colonized, drier, and warmer niche. We report a striking parallelism of directional changes in gene expression plasticity. Thus, plasticity appears to be a key mechanism shaping the initial stages of phenotypic evolution, likely fostering adaptation to novel environments.

Pervasive Introgression During Rapid Diversification of the European Mountain Genus Soldanella (L.) (Primulaceae)

Hybridization is a key mechanism involved in lineage diversification and speciation, especially in ecosystems that experienced repeated environmental oscillations. Recently radiated plant groups, which have evolved in mountain ecosystems impacted by historical climate change provide an excellent model system for studying the impact of gene flow on speciation. We combined organellar (whole-plastome) and nuclear genomic data (RAD-seq) with a cytogenetic approach (rDNA FISH) to investigate the effects of hybridization and introgression on evolution and speciation in the genus Soldanella (snowbells, Primulaceae). Pervasive introgression has already occurred among ancestral lineages of snowbells and has persisted throughout the entire evolutionary history of the genus, regardless of the ecology, cytotype, or distribution range size of the affected species. The highest extent of introgression has been detected in the Carpathian species, which is also reflected in their extensive karyotype variation. Introgression occurred even between species with dysploid and euploid cytotypes, which were considered to be reproductively isolated. The magnitude of introgression detected in snowbells is unprecedented in other mountain genera of the European Alpine System investigated hitherto. Our study stresses the prominent evolutionary role of hybridization in facilitating speciation and diversification on the one hand, but also enriching previously isolated genetic pools. [chloroplast capture; diversification; dysploidy; European Alpine system; introgression; nuclear-cytoplasmic discordance; ribosomal DNA.].

Genome-environment associations along elevation gradients in two snowbed species of the North-Eastern Calcareous Alps

BACKGROUND

Anthropogenic climate change leads to increasing temperatures and altered precipitation and snowmelt patterns, especially in alpine ecosystems. To understand species' responses to climate change, assessment of genetic structure and diversity is crucial as the basis for the evaluation of migration patterns, genetic adaptation potential as well as the identification of adaptive alleles.

RESULTS

We studied genetic structure, diversity and genome-environment associations of two snowbed species endemic to the Eastern Alps with a large elevational range, Achillea clusiana Tausch and Campanula pulla L. Genotyping-by-sequencing was employed to assemble loci de novo, call variants and perform population genetic analyses. Populations of either species were distinguishable by mountain, and to some extent by elevation. We found evidence for gene flow between elevations. Results of genome-environment associations suggested similar selective pressures acting on both species, emanating mainly from precipitation and exposition rather than temperature.

CONCLUSIONS

Given their genetic structure and amount of gene flow among populations the two study species are suitable to serve as a model for genetic monitoring of climate change adaptation along an elevation gradient. Consequences of climate change will predominantly manifest via changes in precipitation and, thus, duration of snow cover in the snowbeds and indirectly via shrub encroachment accompanied by increasing shading of snowbeds at lower range margins. Assembling genomes of the study species and studying larger sample sizes and time series will be necessary to functionally characterize and validate the herein identified genomic loci putatively involved in adaptive processes.

Phenotypic, chemical component and molecular assessment of genetic diversity and population structure of Morinda officinalis germplasm

Background

Morinda officinalis How (MO) is a perennial herb distributed in tropical and subtropical regions, which known as one of the “Four Southern Herbal Medicines”. The extent of genetic variability and the population structure of MO are presently little understood. Here, nine morphological traits, six chemical components and Single nucleotide polymorphism (SNP) markers were used in integrative research of MO germplasm variation among 88 individuals collected from ten populations across four geographical provinces of China.

Results

Both phenotype and chemical composition have significant genetic variation, and there is a certain correlation between them such as root diameter and the nystose content, as well as geographical distribution. The principal component analysis (PCA) showed the leaf length, leaf width, nystose, 1^F-furanosaccharide nystose, and the section color were the major contributors to diversity. The cluster analysis based on phenotypic and oligosaccharide data distinguished three significant groups, which was consistent with the result of a corresponding analysis with 228,615 SNP markers, and importantly, they all showed a significant correlation with geographical origin. However, there was little similarity between two cluster results. The Shannon’s information index (I) varied from 0.17 to 0.53 with a mean of 0.37, suggesting a high level of genetic diversity in MO populations, which mainly existed among individuals within populations, accounting for 99.66% of the total according to the analysis of molecular variance (AMOVA) results. Each population also maintains the connection because of certain gene communication, so that the genetic differentiation between populations was not very significant. The STRUCTURE software was used to analyse the population structure and the result showed that 88 accessions were clustered into three groups, and 67% of them were pure type, which was also confirmed through PCA.

Conclusions

The comprehensive study of phenotypic, chemical and molecular markers will provide valuable information for future breeding plans and understanding the phylogenetic relationship of MO population.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08817-w.

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.

Genetic diversity of Salixlapponum populations in Central Europe

Salixlapponum is a cold-tolerant relict species in Europe that occurs in several sites, probably reflecting previous migration routes of S.lapponum during the Pleistocene. However, only a few data are available on the genetic structures of populations of S.lapponum. In this study, we use PCR-ISSR markers to investigate genetic variation in 19 European populations of S.lapponum L. AMOVA analysis shows that most of the variation (55.8%) occurs within populations; variability among groups accounts for 19.7%. An AMOVA analysis based on four groups determined by STRUCTURE analysis shows similar results: variability of 54.1% within the population and variability of 18.9% between the four population groups, based on geographic regions. Within individual geographic groups, which are characterised by the studied populations, the lowest variability (as well as the highest homogeneity) was found in populations located in Belarus. The obtained results are consistent with our expectations that the European Lowland could be a significant geographic barrier for gene flow over large geographic distances for S.lapponum. Both the Scandinavian and Belarusian populations, as well as those coming from NE Poland, are characterised by significant genetic distinctiveness. However, some populations from NE Poland and the Sudetes show similarities with populations from other geographic regions, indicating existing genetic relationships between them. Moreover, the results suggest a fairly clear division of the population into 4 emerging geographic regions, although separated by a geographical barrier: the Polish lowland, which forms part of the larger geographic unit known as the European Lowland.

Genetic melting pot and importance of long-distance dispersal indicated in the Gladiolus imbricatus L. populations in the Polish Carpathians

The genetic diversity in 11 populations of Gladiolus imbricatus in five mountain ranges, including the Tatra, Pieniny, Gorce, Beskid Niski (Western Carpathians) and Bieszczady Mts (Eastern Carpathians), was studied with inter-simple sequence repeat (ISSR) markers. The species is a perennial plant occurring in open and semi-open sites of anthropogenic origin (meadows and forest margins). We checked a hypothesis on the microrefugial character of the plant populations in the Pieniny Mts, a small calcareous Carpathian range of complicated relief that has never been glaciated. Plant populations in the Tatra and Pieniny Mts had the highest genetic diversity indices, pointing to their long-term persistence. The refugial vs. the non-refugial mountain ranges accounted for a relatively high value of total genetic variation [analysis of molecular variance (AMOVA), 14.12%, p = 0.003]. One of the Pieniny populations was of hybridogenous origin and shared genetic stock with the Tatra population, indicating there is a local genetic melting pot. A weak genetic structuring of populations among particular regions was found (AMOVA, 4.5%, p > 0.05). This could be an effect of the frequent short-distance and sporadic long-distance gene flow. The dispersal of diaspores between the remote populations in the Western Carpathians and Eastern Carpathians could be affected by the historical transportation of flocks of sheep from the Tatra to Bieszczady Mts.

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.

Assisted species migration and hybridization to conserve cold-adapted plants under climate change

Temperature rise due to climate change is putting many arctic and alpine plants at risk of extinction because their ability to react is outpaced by the speed of climate change. We considered assisted species migration (ASM) and hybridization as methods to conserve cold-adapted species (or the genes thereof) and to minimize the potential perturbation of ecosystems due to climate change. Assisted species migration is the deliberate movement of individuals from their current location to where the species' ecological requirements will be matched under climate projections. Hybridization refers to crossbreeding of closely related species, where for arctic and alpine plants, 1 parent is the threatened cold-adapted and the other its reproductively compatible, warm-adapted sibling. Traditionally, hybridization is viewed as negative and leading to a loss of biodiversity, even though hybridization has increased biodiversity over geological times. Furthermore, the incorporation of warm-adapted genes into a hybrid may be the only means for the persistence of increasingly more maladapted, cold-adapted species. If approached with thorough consideration of fitness-related parameters of the source population and acknowledgement of the important role hybridization has played in shaping current biodiversity, ASM and hybridization could help save partial or whole genomes of key cold-adapted species at risk due to climate change with minimal negative effects on ecosystem functioning.

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.

Lack of gene flow between Phytophthora infestans populations of two neighboring countries with the largest potato production

Gene flow is an important evolutionary force that enables adaptive responses of plant pathogens in response to changes in the environment and plant disease management strategies. In this study, we made a direct inference concerning gene flow in the Irish famine pathogen Phytophthora infestans between two of its hosts (potato and tomato) as well as between China and India. This was done by comparing sequence characteristics of the eukaryotic translation elongation factor 1 alpha (eEF-1α) gene, generated from 245 P. infestans isolates sampled from two countries and hosts. Consistent with previous results, we found that eEF-1α gene was highly conserved and point mutation was the only mechanism generating any sequence variation. Higher genetic variation was found in the eEF-1α sequences in the P. infestans populations sampled from tomato compared to those sampled from potato. We also found the P. infestans population from India displayed a higher genetic variation in the eEF-1α sequences compared to China. No gene flow was detected between the pathogen populations from the two countries, which is possibly attributed to the geographic barrier caused by Himalaya Plateau and the minimum cross-border trade of potato and tomato products. The implications of these results for a sustainable management of late blight diseases are discussed.

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.

Diversification and independent polyploid origins in the disjunct species Alyssum repens from the Southeastern Alps and the Carpathians

PREMISE

Disjunct distributions have been commonly observed in mountain plant species and have stimulated phylogeographic and phylogenetic research. Here we studied Alyssum repens, a member of the polyploid species complex A. montanum-A. repens, which exhibits SE Alpine-Carpathian disjunctions with a large elevational span and consists of diploid and tetraploid populations. We aimed to investigate the species' genetic and cytotype structure in the context of its distribution patterns, to elucidate the polyploid origins and to propose an appropriate taxonomic treatment.

METHODS

We combined AFLP fingerprinting markers, sequence variation of the highly repetitive ITS region of rDNA and the low-copy DET1 nuclear gene, genome size, and morphometric data.

RESULTS

We identified four geographically structured genetic lineages. One consisted of diploid populations from the foothills of the Southeastern Alps and neighboring regions, and the three others were allopatric montane to alpine groups comprising diploids and tetraploids growing in the Southeastern Carpathians and the Apuseni Mts. in Romania.

CONCLUSIONS

We inferred a vicariance scenario associated with Quaternary climatic oscillations, accompanied by one auto- and two allopolyploidization events most likely involving a northern Balkan relative. Whereas genetic differentiation and allopatric distribution would favor the taxonomic splitting of this species, the genetic lineages largely lack morphological distinguishability, and their ecological, cytotype and genome size divergence is only partial. Even though we probably face here a case of incipient speciation, we propose to maintain the current taxonomic treatment of Alyssum repens as a single, albeit variable, species.

Divergent selection and genetic structure of Sideritis scardica populations from southern Balkan Peninsula as revealed by AFLP fingerprinting

Sideritis scardica Giseb. is a subalpine/alpine plant species endemic to the central part of the Balkan Peninsula. In this study, we combined Amplified Fragment Length Polymorphism (AFLP) and environmental data to examine the adaptive genetic variations in S. scardica natural populations sampled in contrasting environments. A total of 226 AFLP loci were genotyped in 166 individuals from nine populations. The results demonstrated low gene diversity, ranging from 0.095 to 0.133 and significant genetic differentiation ranging from 0.115 to 0.408. Seven genetic clusters were revealed by Bayesian clustering methods as well as by Discriminant Analysis of Principal Components and each population formed its respective cluster. The exception were populations P02 Mt. Shara and P07 Mt. Vermio, that were admixed between two clusters. Both landscape genetic methods Mcheza and BayeScan identified a total of seven (3.10%) markers exhibiting higher levels of genetic differentiation among populations. The spatial analysis method Samβada detected 50 individual markers (22.12%) associated with bioclimatic variables, among them seven were identified by both Mcheza and BayeScan as being under directional selection. Four bioclimatic variables associated with five out of seven outliers were related to precipitation, suggesting that this variable is the key factor affecting the adaptive variation of S. scardica.

Contrasting evolutionary origins of two mountain endemics: Saxifraga wahlenbergii (Western Carpathians) and S. styriaca (Eastern Alps)

BACKGROUND

The Carpathians and the Alps are the largest mountain ranges of the European Alpine System and important centres of endemism. Among the distinctive endemic species of this area is Saxifraga wahlenbergii, a Western Carpathians member of the speciose genus Saxifraga. It was frequently considered a taxonomically isolated Tertiary palaeopolyploid and palaeoendemic, for which the closest relatives could not yet be traced. A recently described narrow endemic of the Eastern Alps, S. styriaca, was hypothesized to be closely related to S. wahlenbergii based on shared presence of peculiar glandular hairs. To elucidate the origin and phylogenetic relationships of both species we studied nuclear and plastid DNA markers based on multiple accessions and analysed the data in a wide taxonomic context. We applied Sanger sequencing, followed by targeted next-generation sequencing (NGS) for a refined analysis of nrITS variants to detect signatures of ancient hybridization. The ITS data were used to estimate divergence times of different lineages using a relaxed molecular clock.

RESULTS

We demonstrate divergent evolutionary histories for the two mountain endemics. For S. wahlenbergii we revealed a complicated hybrid origin. Its maternal parent belongs to a Western Eurasian lineage of high mountain taxa grouped in subsect. Androsaceae and is most likely the widespread S. androsacea. The putative second parent was most likely S. adscendens, which belongs to the distantly related subsect. Tridactylites. While Sanger sequencing of nrITS only showed S. adscendens-related variants in S. wahlenbergii, our NGS screening revealed presence of sequences from both lineages with clear predominance of the paternal over the maternal lineage.

CONCLUSIONS

Saxifraga styriaca was unambiguously assigned to subsect. Androsaceae and is not the sister taxon of S. wahlenbergii. Accordingly, the similarity of the glandular hairs observed in both taxa rests on parallelism and both species do not constitute an example of a close evolutionary link between the floras of the Western Carpathians and Eastern Alps. With the origin of its paternal, S. adscendens-like ITS DNA estimated to ca. 4.7 Ma, S. wahlenbergii is not a relict of the mid-Tertiary climate optimum. Its hybrid origin is much younger and most likely took place in the Pleistocene.

Phylogeography of Swertia perennis in Europe based on cpDNA markers

Background

Swertia perennis(Gentianaceae) is a perennial diploid and clonal plant species that is discontinuously distributed in peat bogs in the mountains of Europe, Asia and North America as well as in the lowlands of Europe. The current geographical dispersion of S. perennis is probably the result of quaternary climatic changes that have played an important role in determining the distribution of Swertia and other plant and animal species.

Methods

In this study we used molecular techniques and combined data from chloroplast DNA markers (trnLF region and trnH-psbA spacer) to elucidate the phylogeography of S. perennis in Europe. Plants were collected from 28 populations in different locations in the lowlands and mountainous areas of Europe (e.g., the Carpathians, Sudetes, Bohemian Forest and Alps). cDNA was analysed to detect the genetic relationship between specimens from different locations.

Results

A total of 20 haplotypes were identified across the dataset. They were characterised by a high level of genetic variability but showed a lack of phylogeographical structure. This pattern may be the result of repeated recolonization and expansion from several areas. Such genetic differentiation may also be attributed to the relatively long-term isolation of S. perennis in Pleistocene refugia in Europe, which resulted in independent separation of different cpDNA phylogenetic lineages and variation in the nucleotide composition of cpDNA.

Discussion

The lack of strong phylogeographical structure makes it impossible to indicate the centre of haplotype diversity; however, refugia located in the Carpathians, Sudetes or Alps are the most probable sites where S. perennis existed in Europe. This lack of structure may also indicate a high level of gene flow in times when the landscape and fen systems were not fragmented in numerous geographically-isolated populations. This makes it difficult to speculate about the relationships between Asiatic and European plant populations and the origin and distribution of this species in Europe. Today, it seems to be restricted due to the occurrence of plants which clearly reflects the genetic variability from the ancient period.

Molecular markers provide evidence for a broad-fronted recolonisation of the widespread calcareous grassland species Sanguisorba minor from southern and cryptic northern refugia

Calcareous grasslands belong to the most species-rich and endangered habitats in Europe. However, little is known about the origin of the species typically occurring in these grasslands. In this study we analysed the glacial and post-glacial history of Sanguisorba minor, a typical plant species frequently occurring in calcareous grasslands. The study comprised 38 populations throughout the whole distribution range of the species across Europe. We used molecular markers (AFLP) and applied Bayesian cluster analysis as well as spatial principal components analysis (sPCA) to identify glacial refugia and post-glacial migration routes to Central Europe. Our study revealed significant differences in the level of genetic variation and the occurrence of rare fragments within populations of S. minor and a distinct separation of eastern and western lineages. The analyses uncovered traditional southern but also cryptic northern refugia and point towards a broad fronted post-glacial recolonisation. Based on these results we postulate that incomplete lineage sorting may have contributed to the detected pattern of genetic variation and that S. minor recolonised Central Europe post-glacially from Iberia and northern glacial refugia in France, Belgium or Germany. Our results highlight the importance of refugial areas for the conservation of intraspecific variation in calcareous grassland species.

Phylogeography of a tough rock survivor in European dry grasslands

Phylogeographic analyses of plants in Europe have revealed common glacial refugia and migration routes for several trees and herbs with arctic-alpine distributions. The postglacial histories of dry grassland species in central Europe have rarely been analyzed, even though the extremely species-rich habitat is threatened. Sedum album (Crassulaceae) is a common inhabitant of rocky sites in central European dry grasslands. We inferred the phylogeographic history of S. album over its distribution range in Europe. Genetic diversity within and differentiation between 34 S. album populations was examined using AFLP markers. Population isolation was indicated based on the rarity of the fragments and by isolation-by-distance effects. We sequenced the trnL-trnF region in 32 populations and used chloroplast microsatellites to analyze chloroplast haplotype distributions. Two distinct S. album lineages were detected. One lineage was comprised of populations from eastern and central parts of central Europe, and the Apennine Peninsula. A second lineage was comprised of populations from the Iberian Peninsula and western and northern parts of central Europe. Glacial refugia were identified based on the accumulation of ancient chloroplast haplotypes, high diversity of AFLP fragments within populations, and high levels of rare fragments in Liguria, Serbia, the Apennine and Iberian peninsulas. Cryptic refugia were detected in the Czech Republic and Slovakia. Isolation by distance was present all over the distribution range, and it was separately detected in southwestern and central Europe. In western Europe, where a contact zone between the two lineages can be expected, no isolation by distance was detected. Our results suggest migration routes of S. album northeastward from glacial refugia in southern Iberia, northward from the Apennine Peninsula, and northward and westward from the southeastern parts of central Europe. Therefore, central European grasslands were recently colonized by northern cryptic populations and source populations originating in the east and the Apennine Peninsula.

Species distribution modeling and molecular markers suggest longitudinal range shifts and cryptic northern refugia of the typical calcareous grassland species Hippocrepis comosa (horseshoe vetch)

Calcareous grasslands belong to the most diverse, endangered habitats in Europe, but there is still insufficient information about the origin of the plant species related to these grasslands. In order to illuminate this question, we chose for our study the representative grassland species Hippocrepis comosa (Horseshoe vetch). Based on species distribution modeling and molecular markers, we identified the glacial refugia and the postglacial migration routes of the species to Central Europe. We clearly demonstrate that H. comosa followed a latitudinal and due to its oceanity also a longitudinal gradient during the last glacial maximum (LGM), restricting the species to southern refugia situated on the Peninsulas of Iberia, the Balkans, and Italy during the last glaciation. However, we also found evidence for cryptic northern refugia in the UK, the Alps, and Central Germany. Both species distribution modeling and molecular markers underline that refugia of temperate, oceanic species such as H. comosa must not be exclusively located in southern but also in western of parts of Europe. The analysis showed a distinct separation of the southern refugia into a western cluster embracing Iberia and an eastern group including the Balkans and Italy, which determined the postglacial recolonization of Central Europe. At the end of the LGM, H. comosa seems to have expanded from the Iberian refugium, to Central and Northern Europe, including the UK, Belgium, and Germany.

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.

Echoes of a distant time: effects of historical processes on contemporary genetic patterns in Galaxias platei in Patagonia

Interpreting the genetic structure of a metapopulation as the outcome of gene flow over a variety of timescales is essential for the proper understanding of how changes in landscape affect biological connectivity. Here we contrast historical and contemporary connectivity in two metapopulations of the freshwater fish Galaxias platei in northern and southernmost Patagonia where paleolakes existed during the Holocene and Pleistocene, respectively. Contemporary gene flow was mostly high and asymmetrical in the northern system while extremely reduced in the southernmost system. Historical migration patterns were high and symmetric in the northern system and high and largely asymmetric in the southern system. Both systems showed a moderate structure with a clear pattern of isolation by distance (IBD). Effective population sizes were smaller in populations with low contemporary gene flow. An approximate Bayesian computation (ABC) approach suggests a late Holocene colonization of the lakes in the northern system and recent divergence of the populations from refugial populations from east and west of the Andes. For the southern system, the ABC approach reveals that some of the extant G. platei populations most likely derive from an ancestral population inhabiting a large Pleistocene paleolake while the rest derive from a higher-altitude lake. Our results suggest that neither historical nor contemporary processes individually fully explain the observed structure and geneflow patterns and both are necessary for a proper understanding of the factors that affect diversity and its distribution. Our study highlights the importance of a temporal perspective on connectivity to analyse the diversity of spatially complex metapopulations.

Polyploidisation and geographic differentiation drive diversification in a European High Mountain Plant Group (Doronicum clusii Aggregate, Asteraceae)

Range shifts (especially during the Pleistocene), polyploidisation and hybridization are major factors affecting high-mountain biodiversity. A good system to study their role in the European high mountains is the Doronicum clusii aggregate (Asteraceae), whose four taxa (D. clusii s.s., D. stiriacum, D. glaciale subsp. glaciale and D. glaciale subsp. calcareum) are differentiated geographically, ecologically (basiphilous versus silicicolous) and/or via their ploidy levels (diploid versus tetraploid). Here, we use DNA sequences (three plastid and one nuclear spacer) and AFLP fingerprinting data generated for 58 populations to infer phylogenetic relationships, origin of polyploids-whose ploidy level was confirmed by chromosomally calibrated DNA ploidy level estimates-and phylogeographic history. Taxonomic conclusions were informed, among others, by a Gaussian clustering method for species delimitation using dominant multilocus data. Based on molecular data we identified three lineages: (i) silicicolous diploid D. clusii s.s. in the Alps, (ii) silicicolous tetraploid D. stiriacum in the eastern Alps (outside the range of D. clusii s.s.) and the Carpathians and (iii) the basiphilous diploids D. glaciale subsp. glaciale (eastern Alps) and D. glaciale subsp. calcareum (northeastern Alps); each taxon was identified as distinct by the Gaussian clustering, but the separation of D. glaciale subsp. calcareum and D. glaciale subsp. glaciale was not stable, supporting their taxonomic treatment as subspecies. Carpathian and Alpine populations of D. stiriacum were genetically differentiated suggesting phases of vicariance, probably during the Pleistocene. The origin (autopolyploid versus allopolyploid) of D. stiriacum remained unclear. Doronicum glaciale subsp. calcareum was genetically and morphologically weakly separated from D. glaciale subsp. glaciale but exhibited significantly higher genetic diversity and rarity. This suggests that the more widespread D. glaciale subsp. glaciale originated from D. glaciale subsp. calcareum, which is restricted to a prominent Pleistocene refugium previously identified in other alpine plant species.

Key processes for Cheirolophus (Asteraceae) diversification on oceanic islands inferred from AFLP data

The radiation of the genus Cheirolophus (Asteraceae) in Macaronesia constitutes a spectacular case of rapid diversification on oceanic islands. Twenty species - nine of them included in the IUCN Red List of Threatened Species - have been described to date inhabiting the Madeiran and Canarian archipelagos. A previous phylogenetic study revealed that the diversification of Cheirolophus in Macaronesia started less than 2 Ma. As a result of such an explosive speciation process, limited phylogenetic resolution was reported, mainly due to the low variability of the employed molecular markers. In the present study, we used highly polymorphic AFLP markers to i) evaluate species' boundaries, ii) infer their evolutionary relationships and iii) investigate the patterns of genetic diversity in relation to the potential processes likely involved in the radiation of Cheirolophus. One hundred and seventy-two individuals representing all Macaronesian Cheirolophus species were analysed using 249 AFLP loci. Our results suggest that geographic isolation played an important role in this radiation process. This was likely driven by the combination of poor gene flow capacity and a good ability for sporadic long-distance colonisations. In addition, we also found some traces of introgression and incipient ecological adaptation, which could have further enhanced the extraordinary diversification of Cheirolophus in Macaronesia. Last, we hypothesize that current threat categories assigned to Macaronesian Cheirolophus species do not reflect their respective evolutionary relevance, so future evaluations of their conservation status should take into account the results presented here.

Genetic diversity and structure of Dalmatian pyrethrum (Tanacetum cinerariifolium Trevir./Sch./Bip., Asteraceae) within the Balkan refugium

Dalmatian pyrethrum (Tanacetum cinerariifolium Trevir. /Sch./ Bip.) is an outcrossing, perennial insecticidal plant, restricted to the eastern Adriatic coast (Mediterranean). Amplified fragment-length polymorphisms (AFLP) were used to investigate the genetic diversity and structure within and among 20 natural plant populations. The highest level of gene diversity, the number of private alleles and the frequency down-weighted marker values (DW) were found in northern Adriatic populations and gradually decreased towards the southern boundary of the species range. Genetic impoverishment of these southern populations is most likely the result of human-related activities. An analysis of molecular variance (AMOVA) indicated that most of the genetic diversity was attributed to differences among individuals within populations (85.78%), which are expected due to the outcrossing nature of the species. A Bayesian analysis of the population structure identified two dominant genetic clusters. A spatial analysis of the genetic diversity indicated that 5.6% of the genetic differentiation resulted from isolation by distance (IBD), while 12.3% of the genetic differentiation among populations followed the pattern of isolation by environmental distance (IBED). Knowledge of the genetic diversity patterns of the natural populations and the mechanism behind these patterns is required for the exploitation and possible conservation management of this endemic and economically important species.

Analyses of amplified fragment length polymorphisms (AFLP) indicate rapid radiation of Diospyros species (Ebenaceae) endemic to New Caledonia

Background

Radiation in some plant groups has occurred on islands and due to the characteristic rapid pace of phenotypic evolution, standard molecular markers often provide insufficient variation for phylogenetic reconstruction. To resolve relationships within a clade of 21 closely related New Caledonian Diospyros species and evaluate species boundaries we analysed genome-wide DNA variation via amplified fragment length polymorphisms (AFLP).

Results

A neighbour-joining (NJ) dendrogram based on Dice distances shows all species except D. minimifolia, D. parviflora and D. vieillardii to form unique clusters of genetically similar accessions. However, there was little variation between these species clusters, resulting in unresolved species relationships and a star-like general NJ topology. Correspondingly, analyses of molecular variance showed more variation within species than between them. A Bayesian analysis with BEAST produced a similar result. Another Bayesian method, this time a clustering method, Structure, demonstrated the presence of two groups, highly congruent with those observed in a principal coordinate analysis (PCO). Molecular divergence between the two groups is low and does not correspond to any hypothesised taxonomic, ecological or geographical patterns.

Conclusions

We hypothesise that such a pattern could have been produced by rapid and complex evolution involving a widespread progenitor for which an initial split into two groups was followed by subsequent fragmentation into many diverging populations, which was followed by range expansion of then divergent entities. Overall, this process resulted in an opportunistic pattern of phenotypic diversification. The time since divergence was probably insufficient for some species to become genetically well-differentiated, resulting in progenitor/derivative relationships being exhibited in a few cases. In other cases, our analyses may have revealed evidence for the existence of cryptic species, for which more study of morphology and ecology are now required.

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.

Degradation of sexual reproduction in Veronica filiformis after introduction to Europe

BACKGROUND

Baker's law predicts that self-incompatible plant species are generally poor colonizers because their mating system requires a high diversity of genetically differentiated individuals and thus self-compatibility should develop after long-distance dispersal. However, cases like the introduction of the self-incompatible Veronica filiformis (Plantaginaceae) to Europe constitute an often overlooked alternative to this rule. This species was introduced from subalpine areas of the Pontic-Caucasian Mountains and colonized many parts of Central and Western Europe in the last century, apparently without producing seeds. To investigate the consequences of the absence of sexual reproduction in this obligate outcrosser since its introduction, AFLP fingerprints, flower morphology, pollen and ovule production and seed vitality were studied in introduced and native populations.

RESULTS

Interpopulation crossings of 19 introduced German populations performed in the greenhouse demonstrated that introduced populations are often unable to reproduce sexually. These results were similar to intrapopulation crossings, but this depended on the populations used for crossings. Results from AFLP fingerprinting confirmed a lack of genetic diversity in the area of introduction, which is best explained by the dispersal of clones. Flower morphology revealed the frequent presence of mutations affecting the androecium of the flower and decreasing pollen production in introduced populations. The seeds produced in our experiments were smaller, had a lower germination rate and had lower viability than seeds from the native area.

CONCLUSIONS

Taken together, our results demonstrate that V. filiformis was able to spread by vegetative means in the absence of sexual reproduction. This came at the cost of an accumulation of phenotypically observable mutations in reproductive characters, i.e. Muller's ratchet.

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.

A European phylogeography of Rhinanthus minor compared to Rhinanthus angustifolius: unexpected splits and signs of hybridization

Rhinanthus minor and Rhinanthus angustifolius (Orobanchaceae) are annual hemiparasites, which occur sympatrically in Europe and are known to hybridize. We studied chloroplast and nuclear (amplified fragment length polymorphism [AFLP]) diversity in R. minor and compared genetic structuring in this species with R. angustifolius by analyzing the AFLP data for both species simultaneously. The AFLP data revealed that populations in Italy, Greece, and southeast Russia initially identified as R. minor were so distant from the other R. minor populations that they probably belong to another, yet unidentified taxon, and we refer to them as Rhinanthus sp. R. minor s.s. showed a clear geographic genetic structure in both the chloroplast DNA (cpDNA) and nuclear genome. The simultaneous analysis of both species shed new light on the previously published findings for R. angustifolius, because some populations now turned out to belong to R. minor. The admixture analysis revealed very few individuals of mixed R. minor–R.angustifolius ancestry in the natural populations in the west of Europe, while admixture levels were higher in the east. The combined haplotype network showed that haplotype H1 was shared among all species and is likely to be ancestral. H2 was more abundant in R. angustifolius and H3 in R. minor, and the latter probably arose from H1 in this species in the east of Europe. The occurrence of H3 in R. angustifolius may be explained by introgression from R. minor, but without interspecific admixture, these are likely to have been old hybridization events. Our study underlines the importance of including related species in phylogeographic studies.

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.

Bringing together evolution on serpentine and polyploidy: spatiotemporal history of the diploid-tetraploid complex of Knautia arvensis (Dipsacaceae)

Polyploidization is one of the leading forces in the evolution of land plants, providing opportunities for instant speciation and rapid gain of evolutionary novelties. Highly selective conditions of serpentine environments act as an important evolutionary trigger that can be involved in various speciation processes. Whereas the significance of both edaphic speciation on serpentine and polyploidy is widely acknowledged in plant evolution, the links between polyploid evolution and serpentine differentiation have not yet been examined. To fill this gap, we investigated the evolutionary history of the perennial herb Knautia arvensis (Dipsacaceae), a diploid-tetraploid complex that exhibits an intriguing pattern of eco-geographic differentiation. Using plastid DNA sequencing and AFLP genotyping of 336 previously cytotyped individuals from 40 populations from central Europe, we unravelled the patterns of genetic variation among the cytotypes and the edaphic types. Diploids showed the highest levels of genetic differentiation, likely as a result of long term persistence of several lineages in ecologically distinct refugia and/or independent immigration. Recurrent polyploidization, recorded in one serpentine island, seems to have opened new possibilities for the local serpentine genotype. Unlike diploids, the serpentine tetraploids were able to escape from the serpentine refugium and spread further; this was also attributable to hybridization with the neighbouring non-serpentine tetraploid lineages. The spatiotemporal history of K. arvensis allows tracing the interplay of polyploid evolution and ecological divergence on serpentine, resulting in a complex evolutionary pattern. Isolated serpentine outcrops can act as evolutionary capacitors, preserving distinct karyological and genetic diversity. The serpentine lineages, however, may not represent evolutionary 'dead-ends' but rather dynamic systems with a potential to further influence the surrounding populations, e.g., via independent polyplodization and hybridization. The complex eco-geographical pattern together with the incidence of both primary and secondary diploid-tetraploid contact zones makes K. arvensis a unique system for addressing general questions of polyploid research.

Contemporary gene flow and mating system of Arabis alpina in a Central European alpine landscape

BACKGROUND AND AIMS

Gene flow is important in counteracting the divergence of populations but also in spreading genes among populations. However, contemporary gene flow is not well understood across alpine landscapes. The aim of this study was to estimate contemporary gene flow through pollen and to examine the realized mating system in the alpine perennial plant, Arabis alpina (Brassicaceae).

METHODS

An entire sub-alpine to alpine landscape of 2 km(2) was exhaustively sampled in the Swiss Alps. Eighteen nuclear microsatellite loci were used to genotype 595 individuals and 499 offspring from 49 maternal plants. Contemporary gene flow by pollen was estimated from paternity analysis, matching the genotypes of maternal plants and offspring to the pool of likely father plants. Realized mating patterns and genetic structure were also estimated.

KEY RESULTS

Paternity analysis revealed several long-distance gene flow events (≤1 km). However, most outcrossing pollen was dispersed close to the mother plants, and 84 % of all offspring were selfed. Individuals that were spatially close were more related than by chance and were also more likely to be connected by pollen dispersal.

CONCLUSIONS

In the alpine landscape studied, genetic structure occurred on small spatial scales as expected for alpine plants. However, gene flow also covered large distances. This makes it plausible for alpine plants to spread beneficial alleles at least via pollen across landscapes at a short time scale. Thus, gene flow potentially facilitates rapid adaptation in A. alpina likely to be required under ongoing climate change.

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.

Temperature-related natural selection in a wind-pollinated tree across regional and continental scales

Adaptive genetic variation is a key factor in evolutionary biology, but the detection of signatures of natural selection remains challenging in nonmodel organisms. We used a genome scan approach to detect signals of natural selection in the Black alder (Alnus glutinosa), a widespread wind-pollinated tree. Gene flow through pollen dispersal is believed to be high in this species, and we therefore expected to find a clear response to natural selection. In combination with two different landscape genetic approaches, we determined which environmental variables were most associated with the inferred selection. This analysis was performed on a regional scale (northern Belgium) and on a continental scale (Europe). Because climate-related differences are much more pronounced at the continental scale, we expected to find more selection-sensitive genetic markers across Europe than across northern Belgium. At both spatial scales, a substantial number of genetic loci were considered outliers, with respect to neutral expectations, and were therefore identified as selective. Based on results from our combined approach, four putative selective loci (or 2.5%) were recovered with high statistical support. Although these loci seemed to be associated with different environmental variables, they were mainly temperature-related. Our study demonstrates that the use of complementary methods in landscape genetics allows the discovery of selective loci which otherwise might stay hidden. In combination with a genome scan, the selective loci can be verified and the nature of the selection pressure can be identified.