History or ecology? Substrate type as a major driver of patial genetic structure in Alpine plants

Ice holes microrefugia harbor genetically and functionally distinct populations of Vaccinium vitis-idaea (Ericaceae)

In the mountain terrain, ice holes are little depressions between rock boulders that are characterized by the exit of cold air able to cool down the rock surface even in summer. This cold air creates cold microrefugia in warmer surroundings that preserve plant species probably over thousands of years under extra-zonal climatic conditions. We hypothesized that ice hole populations of the model species Vaccinium vitis-idaea (Ericaceae) show genetic differentiation from nearby zonal subalpine populations, and high functional trait distinctiveness, in agreement with genetic patterns. We genotyped almost 30,000 single nucleotide polymorphisms using restriction site-associated DNA sequencing and measured eight functional traits indicative of individual performance and ecological strategies. Genetic results showed high differentiation among the six populations suggesting isolation. On siliceous bedrock, ice hole individuals exhibited higher levels of admixture than those from subalpine populations which could have experienced more bottlenecks during demographic fluctuations related to glacial cycles. Ice hole and subalpine calcareous populations clearly separated from siliceous populations, indicating a possible effect of bedrock in shaping genetic patterns. Trait analysis reflected the bedrock effect on populations' differentiation. The significant correlation between trait and genetic distances suggests the genetic contribution in shaping intraspecific functional differentiation. In conclusion, extra-zonal populations reveal a prominent genetic and phenotypic differentiation determined by history and ecological contingency. Therefore, microrefugia populations can contribute to the overall variability of the species and lead to intraspecific-driven responses to upcoming environmental changes.

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.

Postglacial range expansion of high-elevation plants is restricted by dispersal ability and habitat specialization

AIM

Species' ecological traits influence their spatial genetic patterns. Bedrock preference strongly shapes the phylogeography of alpine plants, but its interactions with other ecological traits have rarely been disentangled. Here, we explore whether dispersal ability and degree of habitat specialization account for divergent postglacial expansion patterns of high-elevation plants in spite of similar bedrock preference.

LOCATION

The Pyrenees, southwestern Europe.

TAXON

Cirsium glabrum (Asteraceae), Salix pyrenaica (Salicaceae) and Silene borderei (Caryophyllaceae).

METHODS

Phylogenetic, genetic structure and demographic modelling analyses based on restriction-site-associated DNA sequencing (RADseq) data from a range-wide populational sampling were conducted. Occurrence data and environmental variables were used to construct species distribution models, which were projected under current and Last Glacial Maximum conditions, and were combined with RADseq data to reconstruct the postglacial history of the study species. The degree of habitat specialization of each species was estimated based on the plant communities within which they occur, and their climatic niche breadth.

RESULTS

Salix pyrenaica, which occupies a broad range of habitats, shows a high level of range filling, a blurred genetic structure and an admixture cline between the two main genetic groups, congruent with rapid postglacial expansion. The microsite specialists C. glabrum and S. borderei exhibit a strong genetic structure and low levels of range filling, indicative of slow postglacial expansion. The good disperser C. glabrum shows higher levels of admixture between genetic groups and weaker population differentiation than the poor disperser S. borderei.

MAIN CONCLUSIONS

Factors other than bedrock preference have a strong impact on the postglacial range dynamics of high-elevation species. Habitat specialization plays an important role, allowing species occupying a broad range of habitats to more rapidly expand their ranges after environmental change. The effect of dispersal ability is lower than expected for the study species.

Tempo and drivers of plant diversification in the European mountain system

There is still limited consensus on the evolutionary history of species-rich temperate alpine floras due to a lack of comparable and high-quality phylogenetic data covering multiple plant lineages. Here we reconstructed when and how European alpine plant lineages diversified, i.e., the tempo and drivers of speciation events. We performed full-plastome phylogenomics and used multi-clade comparative models applied to six representative angiosperm lineages that have diversified in European mountains (212 sampled species, 251 ingroup species total). Diversification rates remained surprisingly steady for most clades, even during the Pleistocene, with speciation events being mostly driven by geographic divergence and bedrock shifts. Interestingly, we inferred asymmetrical historical migration rates from siliceous to calcareous bedrocks, and from higher to lower elevations, likely due to repeated shrinkage and expansion of high elevation habitats during the Pleistocene. This may have buffered climate-related extinctions, but prevented speciation along elevation gradients as often documented for tropical alpine floras.

Here, the authors use full-plastome phylogenomics and multiclade comparative models to reconstruct the tempo and drivers of six European Alpine angiosperm lineages before and during the Pleistocene. They find that geographic divergence and bedrock shifts drive speciation events, while diversification rates remained steady.

Cytogenetic, Morphometric, and Ecological Characterization of Festuca indigesta Boiss. in the Southeast of Spain

Festuca indigesta subsp. indigesta (Poaceae) is endemic to the southeast of Spain, and until recently, it was considered that its range of distribution was restricted to the siliceous core of Sierra Nevada. However, it has been recently extended in the territory to others calcareous mountains. This study investigates the cytogenetic variability throughout the geographic range of this taxon, the possible edaphic preferences of each cytotype, and the morphological variation of cytotypes. Genome sizes and ploidy levels were estimated using flow cytometry and chromosome count. Soil samples were collected to test the nature of the substrate, i.e., pH, and calcium and magnesium contents. Finally, morphological characters were measured in herbarium specimens. This study provides the first genome size data for the species. Hidden cytogenetic diversity was detected in the taxon, comprising hexaploid (2n = 6x = 42), octoploid (2n = 8x = 56) and dodecaploid (2n = 12x = 84) individuals. No relationship between substrate nature and cytotype was observed. Morphological differences were detected for the size of floral parts and stomata among cytotypes, but these were blurred if the entire morphological variation range was considered. Our results suggest that each mountain range could act as a reservoir of morphologically cryptic genetic diversity regarding this taxon.

Life in the desert: The impact of geographic and environmental gradients on genetic diversity and population structure of Ivesia webberi

For range-restricted species with disjunct populations, it is critical to characterize population genetic structure, gene flow, and factors that influence functional connectivity among populations in order to design effective conservation programs. In this study, we genotyped 314 individuals from 16 extant populations of Ivesia webberi, a United States federally threatened Great Basin Desert using six microsatellite loci. We assessed the effects of Euclidean distance, landscape features, and ecological dissimilarity on the pairwise genetic distance of the sampled populations, while also testing for a potential relationship between I. webberi genetic diversity and diversity in the vegetative communities. The results show low levels of genetic diversity overall (H e = 0.200-0.441; H o = 0.192-0.605) and high genetic differentiation among populations. Genetic diversity was structured along a geographic gradient, congruent with patterns of isolation by distance. Populations near the species' range core have relatively high genetic diversity, supporting in part a central-marginal pattern, while also showing some evidence for a metapopulation dynamic. Peripheral populations have lower genetic diversity, significantly higher genetic distances, and higher relatedness. Genotype cluster admixture results suggest a complex dispersal pattern among populations with dispersal direction and distance varying on the landscape. Pairwise genetic distance strongly correlates with elevation, actual evapotranspiration, and summer seasonal precipitation, indicating a role for isolation by environment, which the observed phenological mismatches among the populations also support. The significant correlation between pairwise genetic distance and floristic dissimilarity in the germinated soil seed bank suggests that annual regeneration in the plant communities contribute to the maintenance of genetic diversity in I. webberi.

History or demography? Determining the drivers of genetic variation in North American plants

Understanding the impact of historical and demographic processes on genetic variation is essential for devising conservation strategies and predicting responses to climate change. Recolonization after Pleistocene glaciations is expected to leave distinct genetic signatures, characterised by lower genetic diversity in previously glaciated regions. Populations’ positions within species ranges also shape genetic variation, following the central‐marginal paradigm dictating that peripheral populations are depauperate, sparse and isolated. However, the general applicability of these patterns and relative importance of historical and demographic factors remains unknown. Here, we analysed the distribution of genetic variation in 91 native species of North American plants by coupling microsatellite data and species distribution modelling. We tested the contributions of historical climatic shifts and the central‐marginal hypothesis on genetic diversity and structure on the whole data set and across subsets based on taxonomic groups and growth forms. Decreased diversity was found with increased distance from potential glacial refugia, coinciding with the expected make‐up of postglacially colonised localities. At the range periphery, lower genetic diversity, higher inbreeding levels and genetic differentiation were reported, following the assumptions of the central‐marginal hypothesis. History and demography were found to have approximately equal importance in shaping genetic variation.

Arabis alpina: A perennial model plant for ecological genomics and life-history evolution

Many model organisms were chosen and achieved prominence because of an advantageous combination of their life‐history characteristics, genetic properties and also practical considerations. Discoveries made in Arabidopsis thaliana, the most renowned noncrop plant model species, have markedly stimulated studies in other species with different biology. Within the family Brassicaceae, the arctic–alpine Arabis alpina has become a model complementary to Arabidopsis thaliana to study the evolution of life‐history traits, such as perenniality, and ecological genomics in harsh environments. In this review, we provide an overview of the properties that facilitated the rapid emergence of A. alpina as a plant model. We summarize the evolutionary history of A. alpina, including genomic aspects, the diversification of its mating system and demographic properties, and we discuss recent progress in the molecular dissection of developmental traits that are related to its perennial life history and environmental adaptation. From this published knowledge, we derive open questions that might inspire future research in A. alpina, other Brassicaceae species or more distantly related plant families.

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.

Discovery of cryptic plant diversity on the rooftops of the Alps

High elevation temperate mountains have long been considered species poor owing to high extinction or low speciation rates during the Pleistocene. We performed a phylogenetic and population genomic investigation of an emblematic high-elevation plant clade (Androsace sect. Aretia, 31 currently recognized species), based on plant surveys conducted during alpinism expeditions. We inferred that this clade originated in the Miocene and continued diversifying through Pleistocene glaciations, and discovered three novel species of Androsace dwelling on different bedrock types on the rooftops of the Alps. This highlights that temperate high mountains have been cradles of plant diversity even during the Pleistocene, with in-situ speciation driven by the combined action of geography and geology. Our findings have an unexpected historical relevance: H.-B. de Saussure likely observed one of these species during his 1788 expedition to the Mont Blanc and we describe it here, over two hundred years after its first sighting.

Range and niche expansion through multiple interspecific hybridization: a genotyping by sequencing analysis of Cherleria (Caryophyllaceae)

BACKGROUND

Cherleria (Caryophyllaceae) is a circumboreal genus that also occurs in the high mountains of the northern hemisphere. In this study, we focus on a clade that diversified in the European High Mountains, which was identified using nuclear ribosomal (nrDNA) sequence data in a previous study. With the nrDNA data, all but one species was monophyletic, with little sequence variation within most species. Here, we use genotyping by sequencing (GBS) data to determine whether the nrDNA data showed the full picture of the evolution in the genomes of these species.

RESULTS

The overall relationships found with the GBS data were congruent with those from the nrDNA study. Most of the species were still monophyletic and many of the same subclades were recovered, including a clade of three narrow endemic species from Greece and a clade of largely calcifuge species. The GBS data provided additional resolution within the two species with the best sampling, C. langii and C. laricifolia, with structure that was congruent with geography. In addition, the GBS data showed significant hybridization between several species, including species whose ranges did not currently overlap.

CONCLUSIONS

The hybridization led us to hypothesize that lineages came in contact on the Balkan Peninsula after they diverged, even when those lineages are no longer present on the Balkan Peninsula. Hybridization may also have helped lineages expand their niches to colonize new substrates and different areas. Not only do genome-wide data provide increased phylogenetic resolution of difficult nodes, they also give evidence for a more complex evolutionary history than what can be depicted by a simple, branching phylogeny.

Slow and steady wins the race: contrasted phylogeographic signatures in two Alpine amphibians

A deeper phylogeographic structure is expected for slow-dispersing habitat specialists compared to widespread adaptable species, especially in topographically complex regions. We tested this classic assumption by comparing the genomic (RAD-sequencing) phylogeographies of two amphibians inhabiting the Swiss Alps: the mobile, cosmopolitan common frog (Rana temporaria) against the stationary, mountain endemic Alpine salamander (Salamandra atra). Our results ran opposite of predictions: the frog displayed significantly higher genetic divergences and lower within-population variation compared to the salamander. This implies a prominent role for their distinctive glacial histories in shaping intraspecific diversity and structure: diversification and recolonization from several circum-Alpine micro-refugia for the frog versus a single refugium for the salamander, potentially combined with better population connectivity and stability. These striking differences emphasize the great variability of phylogeographic responses to the Quaternary glaciations, hence the complexity to predict general patterns of genetic diversity at the regional scale, and the forces that underlie them.

Refugia within refugium of Geranium yesoense (Geraniaceae) in Japan were driven by recolonization into the southern interglacial refugium

Recent studies have found that geographical fragmentation and recurrent colonization result in complex genetic structures in refugial areas. This phenomenon, known as ‘refugia within refugium’, has been identified from many geographical locations. In Japan, the high-elevation mountains of central Honshu provided an interglacial refugium for alpine plants. Here we focused on the Geranium yesoense complex, which exhibits increased morphological variation in the refugial area, to determine whether this variation was shaped by recurrent colonization, range fragmentation or phenotypic changes independent of population history. We analysed single nucleotide polymorphism data and chloroplast genome sequences. Diversification in the G. yesoense species complex occurred in the mid-Pleistocene. The varieties are distinct entities and suggest the presence of a genetic cluster with highly disjunct distributions, occurring both in northern Japan and in southern refugial areas in central Honshu. Demographic analysis suggests that a single ancestral variety (var. nipponicum) evolved in the alpine region of central Honshu, and that subsequent migration from one of the two diverged northern varieties (var. pseudopratense) led to secondary contact with var. nipponicum during the last glacial period. Recolonization into refugial populations in central Honshu and hybridization between diverged populations have resulted in complex genetic structures among refugial populations.

As old as the hills: Pliocene palaeogeographical processes influence patterns of genetic structure in the widespread, common shrub Banksia sessilis

The impact of Quaternary glaciation on the development of phylogeographic structure in plant species is well documented. In unglaciated landscapes, phylogeographic patterns tend to reflect processes relating to persistence and stochasticity, yet other factors, associated with the palaeogeographical history of the landscape, including geomorphological events, can also have a significant influence. The unglaciated landscape of south-western Western Australia is an ideal location to observe these ancient drivers of lineage diversification, with tectonic activity associated with the Darling Fault in the late Pliocene attributed to patterns of deep phylogeographic divergence in a widespread tree from this region. Interestingly, other species within this region have not shown this pattern and this palaeogeographical boundary therefore presents an opportunity to examine age and historical distribution of plant species endemic to this region. In this study, we assess patterns of genetic diversity and structure across 28 populations of the widespread shrub Banksia sessilis using three cpDNA markers and nine nuclear microsatellite markers. Sixteen cpDNA haplotypes were identified, comprising two major chloroplast DNA lineages that are estimated to have diverged in the Pliocene, approximately 3.3 million years ago. This timing coincides with major geomorphological processes in the landscape, including the separation of the Darling Plateau from the adjacent Swan Coastal Plain, as well as eustatic changes on the Swan Coastal Plain that are likely to have resulted in the physical isolation of historical plant lineages. Chloroplast lineages were broadly aligned with populations associated with older lateritic soils of the Darling Plateau and Geraldton sandplains or the younger sandy soils associated with the Swan Coastal Plain and Southern Coastline. This structural pattern of lateritic versus non-lateritic division was not observed in the nuclear microsatellite data that identified three genetic clades that roughly corresponded to populations in the North, South, and Central portions of the distributions.

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.

Comparative studies on population genetic structure of two closely related selfing and outcrossing Zingiber species in Hainan Island

How mating system impacts the genetic diversity of plants has long fascinated and puzzled evolutionary biologists. Numerous studies have shown that self-fertilising plants have less genetic diversity at both the population and species levels than outcrossers. However, the phylogenetic relationships between species and correlated ecological traits have not been accounted for in these previous studies. Here, we conduct a comparative population genetic study of two closely related selfing and outcrossing Zingiber species, with sympatric distribution in Hainan Island, and obtain a result contrary to previous studies. The results indicate that selfing Z. corallinum can maintain high genetic diversity through differentiation intensified by local adaptation in populations across the species’ range. In contrast, outcrossing Z. nudicarpum preserves high genetic diversity through gene exchange by frequent export of pollen within or among populations. Contrary to expectations, the major portion of genetic variation of outcrossing Z. nudicarpum may exist among populations, depending on the dispersal ability of pollen and seed. Our results also reveal that the main factor affecting population structure of selfing Z. corallinum is mountain ranges, followed by a moist climate, while that of outcrossing Z. nudicarpum is likely moisture, but not mountain ranges, due to gene flow via pollen.

Plant evolution in alkaline magnesium-rich soils: A phylogenetic study of the Mediterranean genus Hormathophylla (Cruciferae: Alysseae) based on nuclear and plastid sequences

Habitats with alkaline edaphic substrates are often associated with plant speciation and diversification. The tribe Alysseae, in the family Brassicaceae, epitomizes this evolutionary trend. In this lineage, some genera, like Hormathophylla, can serve as a good case for testing the evolutionary framework. This genus is centered in the western Mediterranean. It grows on different substrates, but mostly on alkaline soils. It has been suggested that diversification in many lineages of the tribe Alysseae and in the genus Hormathophylla is linked to a tolerance for high levels of Mg+2 in xeric environments. In this study, we investigated the controversial phylogenetic placement of Hormathophylla in the tribe, the generic limits and the evolutionary relationships between the species using ribosomal and plastid DNA sequences. We also examined the putative association between the evolution of different ploidy levels, trichome morphology and the type of substrates. Our analyses demonstrated the monophyly of the genus Hormathophylla including all previously described species. Nuclear sequences revealed two lineages that differ in basic chromosome numbers (x = 7 and x = 8 or derived 11, 15) and in their trichome morphology. Contrasting results with plastid genes indicates more complex relationships between these two lineages involving recent hybridization processes. We also found an association between chloroplast haplotypes and substrate, especially in populations growing on dolomites. Finally, our dated phylogeny demonstrates that the origin of the genus took place in the mid-Miocene, during the establishment of temporal land bridges between the Tethys and Paratethys seas, with a later diversification during the upper Pliocene.

Glacial refugia and speciation in a group of wind-pollinated and -dispersed, endemic Alpine species of Helictotrichon (Poaceae)

In the Alps phylogeographic studies indicate for small insect-pollinated herbs that climatic fluctuations caused significant population migrations and fragmentations into glacial refugia at the periphery of the Alps. Here we investigate whether this holds also for wind-pollinated and -dispersed species. We therefore analysed the phylogeographic pattern (nuclear and chloroplast dataset) of a clade of the four species of the Helictotrichon parlatorei species group (Poaceae) endemic to the Alps. In contrast to earlier findings for small insect-pollinated herbs no clear barriers to gene flow could be detected in this species group. Instead a few haplotypes are widespread across the entire Alpine region. While the complete absence of a phylogeographic structure in the plastid dataset hints towards very efficient long distance seed dispersal, the moderate phylogeographic structure in the nuclear dataset indicates at least some spatial restriction to pollen dispersal. Rare haplotypes cluster solely in the Western and Southern central Alps and thereby suggest this to be the area of origin for the H. parlatorei species group from where expansion occurred following the presence of calcareous bedrock into the Eastern Alps. We thus conclude that the inclusion of taxa with complementary life-history traits is vital in understanding the glacial history of the Alpine flora.

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.

Subsampling reveals that unbalanced sampling affects STRUCTURE results in a multi-species dataset

Studying the genetic population structure of species can reveal important insights into several key evolutionary, historical, demographic, and anthropogenic processes. One of the most important statistical tools for inferring genetic clusters is the program STRUCTURE. Recently, several papers have pointed out that STRUCTURE may show a bias when the sampling design is unbalanced, resulting in spurious joining of underrepresented populations and spurious separation of overrepresented populations. Suggestions to overcome this bias include subsampling and changing the ancestry model, but the performance of these two methods has not yet been tested on actual data. Here, I use a data set of 12 high-alpine plant species to test whether unbalanced sampling affects the STRUCTURE inference of population differentiation between the European Alps and the Carpathians. For four of the 12 species, subsampling of the Alpine populations-to match the sample size between the Alps and the Carpathians-resulted in a drastically different clustering than the full data set. On the other hand, STRUCTURE results with the alternative ancestry model were indistinguishable from the results with the default model. Based on these results, the subsampling strategy seems a more viable approach to overcome the bias than the alternative ancestry model. However, subsampling is only possible when there is an a priori expectation of what constitute the main clusters. Though these results do not mean that the use of STRUCTURE should be discarded, it does indicate that users of the software should be cautious about the interpretation of the results when sampling is unbalanced.

Environmental effects on fine-scale spatial genetic structure in four Alpine keystone forest tree species

Genetic responses to environmental changes take place at different spatial scales. While the effect of environment on the distribution of species' genetic diversity at large geographical scales has been the focus of several recent studies, its potential effects on genetic structure at local scales are understudied. Environmental effects on fine-scale spatial genetic structure (FSGS) were investigated in four Alpine conifer species (five to eight populations per species) from the eastern Italian Alps. Significant FSGS was found for 11 of 25 populations. Interestingly, we found no significant differences in FSGS across species but great variation among populations within species, highlighting the importance of local environmental factors. Interannual variability in spring temperature had a small but significant effect on FSGS of Larix decidua, probably related to species-specific life history traits. For Abies alba, Picea abies and Pinus cembra, linear models identified spring precipitation as a potentially relevant climate factor associated with differences in FSGS across populations; however, models had low explanatory power and were strongly influenced by a P. cembra outlier population from a very dry site. Overall, the direction of the identified effects is according to expectations, with drier and more variable environments increasing FSGS. Underlying mechanisms may include climate-related changes in the variance of reproductive success and/or environmental selection of specific families. This study provides new insights on potential changes in local genetic structure of four Alpine conifers in the face of environmental changes, suggesting that new climates, through altering FSGS, may also have relevant impacts on plant microevolution.

Environmental and geographic variables are effective surrogates for genetic variation in conservation planning

Protected areas buffer species from anthropogenic threats and provide places for the processes that generate and maintain biodiversity to continue. However, genetic variation, the raw material for evolution, is difficult to capture in conservation planning, not least because genetic data require considerable resources to obtain and analyze. Here we show that freely available environmental and geographic distance variables can be highly effective surrogates in conservation planning for representing adaptive and neutral intraspecific genetic variation. We obtained occurrence and genetic data from the IntraBioDiv project for 27 plant species collected over the European Alps using a gridded sampling scheme. For each species, we identified loci that were potentially under selection using outlier loci methods, and mapped their main gradients of adaptive and neutral genetic variation across the grid cells. We then used the cells as planning units to prioritize protected area acquisitions. First, we verified that the spatial patterns of environmental and geographic variation were correlated, respectively, with adaptive and neutral genetic variation. Second, we showed that these surrogates can predict the proportion of genetic variation secured in randomly generated solutions. Finally, we discovered that solutions based only on surrogate information secured substantial amounts of adaptive and neutral genetic variation. Our work paves the way for widespread integration of surrogates for genetic variation into conservation planning.

Simultaneous speciation in the European high mountain flowering plant genus Facchinia (Minuartia s.l., Caryophyllaceae) revealed by genotyping-by-sequencing

Understanding the relative importance of different mechanisms of speciation in a given lineage requires fully resolved interspecific relationships. Using Facchinia, a genus of seven species centred in the European Alps, we explore whether the polytomy found by Sanger sequencing analyses of standard nuclear (ITS) and plastid markers (trnQ-rps16) is a hard or soft polytomy by substantially increasing the amount of DNA sequence data, generated by genotyping-by-sequencing. In comparison to 142 phylogenetically informative sites in the Sanger sequences the GBS sequences yielded 3363 phylogenetically informative sites after exclusion of apparently oversaturated SNPs. Maximum parsimony, maximum likelihood, NeighborNet, SVDquartets and Astral-II analyses all resulted in phylogenetic trees (and networks) in which interspecific relationships were largely unresolved. After excluding incomplete lineage sorting, hybridisation and oversaturation of characters as possible causes for lack of phylogenetic resolution, we conclude that the polytomy obtained most likely represents a hard polytomy. We hypothesize that diversification of Facchinia is best interpreted as the result of multiple simultaneous vicariance in response to climatic changes during the Early Quaternary.

The DNA of coral reef biodiversity: predicting and protecting genetic diversity of reef assemblages

Conservation of ecological communities requires deepening our understanding of genetic diversity patterns and drivers at community-wide scales. Here, we use seascape genetic analysis of a diversity metric, allelic richness (AR), for 47 reef species sampled across 13 Hawaiian Islands to empirically demonstrate that large reefs high in coral cover harbour the greatest genetic diversity on average. We found that a species's life history (e.g. depth range and herbivory) mediates response of genetic diversity to seascape drivers in logical ways. Furthermore, a metric of combined multi-species AR showed strong coupling to species richness and habitat area, quality and stability that few species showed individually. We hypothesize that macro-ecological forces and species interactions, by mediating species turnover and occupancy (and thus a site's mean effective population size), influence the aggregate genetic diversity of a site, potentially allowing it to behave as an apparent emergent trait that is shaped by the dominant seascape drivers. The results highlight inherent feedbacks between ecology and genetics, raise concern that genetic resilience of entire reef communities is compromised by factors that reduce coral cover or available habitat, including thermal stress, and provide a foundation for new strategies for monitoring and preserving biodiversity of entire reef ecosystems.

Niche shifts and range expansions along cordilleras drove diversification in a high-elevation endemic plant genus in the tropical Andes

The tropical Andes represent one of the world's biodiversity hot spots, but the evolutionary drivers generating their striking species diversity still remain poorly understood. In the treeless high-elevation Andean environments, Pleistocene glacial oscillations and niche differentiation are frequently hypothesized diversification mechanisms; however, sufficiently densely sampled population genetic data supporting this are still lacking. Here, we reconstruct the evolutionary history of Loricaria (Asteraceae), a plant genus endemic to the Andean treeless alpine zone, based on comprehensive population-level sampling of 289 individuals from 67 populations across the entire distribution ranges of its northern Andean species. Partly incongruent AFLP and plastid DNA markers reveal that the distinct genetic structure was shaped by a complex interplay of biogeography (spread along and across the cordilleras), history (Pleistocene glacial oscillations) and local ecological conditions. While plastid variation documents an early split or colonization of the northern Andes by at least two lineages, one of which further diversified, a major split in the AFLP data correlate with altitudinal ecological differentiation. This suggests that niche shifts may be important drivers of Andean diversification not only in forest-alpine transitions, but also within the treeless alpine zone itself. The patterns of genetic differentiation at the intraspecific level reject the hypothesized separation in spatially isolated cordilleras and instead suggest extensive gene flow among populations from distinct mountain chains. Our study highlights that leveraging highly variable markers against extensive population-level sampling is a promising approach to address mechanisms of rapid species diversifications.

Toward a paradigm shift in comparative phylogeography driven by trait-based hypotheses

For three decades, comparative phylogeography has conceptually and methodologically relied on the concordance criterion for providing insights into the historical/biogeographic processes driving population genetic structure and divergence. Here we discuss how this emphasis, and the corresponding lack of methods for extracting information about biotic/intrinsic contributions to patterns of genetic variation, may bias our general understanding of the factors driving genetic structure. Specifically, this emphasis has promoted a tendency to attribute discordant phylogeographic patterns to the idiosyncracies of history, as well as an adherence to generic null expectations of concordance with reduced predictive power. We advocate that it is time for a paradigm shift in comparative phylogeography, especially given the limited utility of the concordance criterion as genomic data provide ever-increasing levels of resolution. Instead of adhering to the concordance-discordance dichotomy, comparative phylogeography needs to emphasize the contribution of taxon-specific traits that will determine whether concordance is a meaningful criterion for evaluating hypotheses or may predict discordant phylogeographic structure. Through reference to some case studies we illustrate how refined hypotheses based on taxon-specific traits can provide improved predictive frameworks to forecast species responses to climatic change or biogeographic barriers while gaining unique insights about the taxa themselves and their interactions with their environment. We outline a potential avenue toward a synthetic comparative phylogeographic paradigm that includes addressing some important conceptual and methodological challenges related to study design and application of model-based approaches for evaluating support of trait-based hypotheses under the proposed paradigm.

Genomic and phenotypic differentiation of Arabidopsis thaliana along altitudinal gradients in the North Italian Alps

Altitudinal gradients in mountain regions are short-range clines of different environmental parameters such as temperature or radiation. We investigated genomic and phenotypic signatures of adaptation to such gradients in five Arabidopsis thaliana populations from the North Italian Alps that originated from 580 to 2350 m altitude by resequencing pools of 19-29 individuals from each population. The sample includes two pairs of low- and high-altitude populations from two different valleys. High-altitude populations showed a lower nucleotide diversity and negative Tajima's D values and were more closely related to each other than to low-altitude populations from the same valley. Despite their close geographic proximity, demographic analysis revealed that low- and high-altitude populations split between 260 000 and 15 000 years before present. Single nucleotide polymorphisms whose allele frequencies were highly differentiated between low- and high-altitude populations identified genomic regions of up to 50 kb length where patterns of genetic diversity are consistent with signatures of local selective sweeps. These regions harbour multiple genes involved in stress response. Variation among populations in two putative adaptive phenotypic traits, frost tolerance and response to light/UV stress was not correlated with altitude. Taken together, the spatial distribution of genetic diversity reflects a potentially adaptive differentiation between low- and high-altitude populations, whereas the phenotypic differentiation in the two traits investigated does not. It may resemble an interaction between adaptation to the local microhabitat and demographic history influenced by historical glaciation cycles, recent seed dispersal and genetic drift in local populations.

Phylogeography of Arabidopsis halleri (Brassicaceae) in mountain regions of Central Europe inferred from cpDNA variation and ecological niche modelling

The present study aimed to investigate phylogeographical patterns present within A. halleri in Central Europe. 1,281 accessions sampled from 52 populations within the investigated area were used in the study of genetic variation based on chloroplast DNA. Over 500 high-quality species occurrence records were used in ecological niche modelling experiments. We evidenced the presence of a clear phylogeographic structure within A. halleri in Central Europe. Our results showed that two genetically different groups of populations are present in western and eastern part of the Carpathians. The hypothesis of the existence of a glacial refugium in the Western Carpathians adn the Bohemian Forest cannot be rejected from our data. It seems, however, that the evidence collected during the present study is not conclusive. The area of Sudetes was colonised after LGM probably by migrants from the Bohemian Forest.

Evolution of the Selfing Syndrome in Arabis alpina (Brassicaceae)

Introduction

The transition from cross-fertilisation (outcrossing) to self-fertilisation (selfing) frequently coincides with changes towards a floral morphology that optimises self-pollination, the selfing syndrome. Population genetic studies have reported the existence of both outcrossing and selfing populations in Arabis alpina (Brassicaceae), which is an emerging model species for studying the molecular basis of perenniality and local adaptation. It is unknown whether its selfing populations have evolved a selfing syndrome.

Methods

Using macro-photography, microscopy and automated cell counting, we compared floral syndromes (size, herkogamy, pollen and ovule numbers) between three outcrossing populations from the Apuan Alps and three selfing populations from the Western and Central Alps (Maritime Alps and Dolomites). In addition, we genotyped the plants for 12 microsatellite loci to confirm previous measures of diversity and inbreeding coefficients based on allozymes, and performed Bayesian clustering.

Results and Discussion

Plants from the three selfing populations had markedly smaller flowers, less herkogamy and lower pollen production than plants from the three outcrossing populations, whereas pistil length and ovule number have remained constant. Compared to allozymes, microsatellite variation was higher, but revealed similar patterns of low diversity and high Fis in selfing populations. Bayesian clustering revealed two clusters. The first cluster contained the three outcrossing populations from the Apuan Alps, the second contained the three selfing populations from the Maritime Alps and Dolomites.

Conclusion

We conclude that in comparison to three outcrossing populations, three populations with high selfing rates are characterised by a flower morphology that is closer to the selfing syndrome. The presence of outcrossing and selfing floral syndromes within a single species will facilitate unravelling the genetic basis of the selfing syndrome, and addressing which selective forces drive its evolution.

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.

Validation of outlier loci through replication in independent data sets: a test on Arabis alpina

Outlier detection and environmental association analysis are common methods to search for loci or genomic regions exhibiting signals of adaptation to environmental factors. However, a validation of outlier loci and corresponding allele distribution models through functional molecular biology or transplant/common garden experiments is rarely carried out. Here, we employ another method for validation, namely testing outlier loci in specifically designed, independent data sets. Previously, an outlier locus associated with three different habitat types had been detected in Arabis alpina. For the independent validation data set, we sampled 30 populations occurring in these three habitat types across five biogeographic regions of the Swiss Alps. The allele distribution model found in the original study could not be validated in the independent test data set: The outlier locus was no longer indicative of habitat-mediated selection. We propose several potential causes of this failure of validation, of which unaccounted genetic structure and technical issues in the original data set used to detect the outlier locus were most probable. Thus, our study shows that validating outlier loci and allele distribution models in independent data sets is a helpful tool in ecological genomics which, in the case of positive validation, adds confidence to outlier loci and their association with environmental factors or, in the case of failure of validation, helps to explain inconsistencies.

Genetic Drift Linked to Heterogeneous Landscape and Ecological Specialization Drives Diversification in the Alpine Endemic Columbine Aquilegia thalictrifolia

The European Alpine system is an extensive mountain range, whose heterogeneous landscape together with Quaternary climatic oscillations significantly affected organismal diversity and distribution in Europe. The model genus Aquilegia represents a textbook example of a rapid and recent radiation through the Northern hemisphere, with the majority of the European taxa occuring in the Alpine system. However, the processes governing genetic differentiation of the genus in this complex geographic area are still widely unexplored. In this work, we used 9 microsatellite loci to study the genetic structure and diversity of 11 populations of Aquilegia thalictrifolia Schott & Kotschy, an alpine taxon characterized by a marked ecological specificity. We found that, despite the endemic and fragmented distribution, A. thalictrifolia has overall high levels of heterozygosity, which is consistent to the substantial inbreeding depression that characterizes the genus. Strong spatial genetic structuring of populations suggests a historical prevalence of genetic drift over gene flow, with natural barriers and ecological niche hindering migration. An analytical comparison of fixation and population differentiation indexes allowed us to infer hypotheses of the postglacial history and more recent demographic events that have influenced the genetics of the species. Overall, our results indicate allopatry as a major force of differentiation in the European scenario, likely to underlie the development of taxonomic boundaries in a broader geographic context. This adds to previous notions on the primary evolutionary forces shaping the Aquilegia radiation in Europe.