Divergent selection in trailing- versus leading-edge populations of Biscutella laevigata

Patterns of Leaf and Fruit Morphological Variation in Marginal Populations of Acer tataricum L. subsp. tataricum

Marginal populations are usually smaller and more isolated and grow in less favourable conditions than those at the distribution centre. The variability of these populations is of high importance, as it can support the adaptations needed for the conditions that they grow in. In this research, the morphological variability of eight Tatar maple (Acer tataricum L. subsp. tataricum) populations was analysed. Tatar maple is an insect-pollinated and wind-dispersed shrub/tree, whose northwestern distribution edge is in southeastern Europe. Morphometric methods were used to analyse the variability of the populations using leaf and fruit morphology. The research revealed significant differences between and within populations. Furthermore, differences in the distribution of the total variability were noted, which suggest that different evolutionarily factors affect different plant traits. Correlation analysis confirmed a weak dependency between the vegetative and generative traits. In addition, no evidence was found for the presence of isolation by environment (IBE). However, the Mantel test for isolation by distance (IBD) was significant for the leaf morphometric traits and non-significant for the fruit morphometric traits. Being the marginal leading-edge populations, they are younger and were less likely to have had time for adaptation to local environments, which would have resulted in the development of IBE. Overall, edge populations of Tatar maple were characterised by great morphological variability, which helps these populations in their response to the intensive selective pressures they face in their environment.

Signature of mid-Pleistocene lineages in the European silver fir (Abies alba Mill.) at its geographic distribution margin

In a conservation and sustainable management perspective, we identify the ecological, climatic, and demographic factors responsible for the genetic diversity patterns of the European silver fir (Abies alba Mill.) at its southwestern range margin (Pyrenees Mountains, France, Europe). We sampled 45 populations throughout the French Pyrenees and eight neighboring reference populations in the Massif Central, Alps, and Corsica. We genotyped 1,620 individuals at three chloroplast and ten nuclear microsatellite loci. We analyzed within- and among-population genetic diversity using phylogeographic reconstructions, tests of isolation-by-distance, Bayesian population structure inference, modeling of demographic scenarios, and regression analyses of genetic variables with current and past environmental variables. Genetic diversity decreased from east to west suggesting isolation-by-distance from the Alps to the Pyrenees and from the Eastern to the Western Pyrenees. We identified two Pyrenean lineages that diverged from a third Alpine-Corsica-Massif Central lineage 0.8 to 1.1 M years ago and subsequently formed a secondary contact zone in the Central Pyrenees. Population sizes underwent contrasted changes, with a contraction in the west and an expansion in the east. Glacial climate affected the genetic composition of the populations, with the western genetic cluster only observed in locations corresponding to the coldest past climate and highest elevations. The eastern cluster was observed over a larger range of temperatures and elevations. All demographic events shaping the current spatial structure of genetic diversity took place during the Mid-Pleistocene Transition, long before the onset of the Holocene. The Western Pyrenees lineage may require additional conservation efforts, whereas the eastern lineage is well protected in in situ gene conservation units. Due to past climate oscillations and the likely emergence of independent refugia, east-west oriented mountain ranges may be important reservoir of genetic diversity in a context of past and ongoing climate change in Europe.

Invoking adaptation to decipher the genetic legacy of past climate change

Persistence of natural populations during periods of climate change is likely to depend on migration (range shifts) or adaptation. These responses were traditionally considered discrete processes and conceptually divided into the realms of ecology and evolution. In a milestone paper, Davis and Shaw (2001) Science 292:673 argued that the interplay of adaptation and migration was central to biotic responses to Quaternary climate, but since then there has been no synthesis of efforts made to set up this research program. Here we review some of the salient findings from molecular genetic studies assessing ecological and evolutionary responses to Quaternary climate change. These studies have revolutionized our understanding of population processes associated with past species migration. However, knowledge remains limited about the role of natural selection for local adaptation of populations to Quaternary environmental fluctuations and associated range shifts, and for the footprints this might have left on extant populations. Next-generation sequencing technologies, high-resolution paleoclimate analyses, and advances in population genetic theory offer an unprecedented opportunity to test hypotheses about adaptation through time. Recent population genomics studies have greatly improved our understanding of the role of contemporary adaptation to local environments in shaping spatial patterns of genetic diversity across modern-day landscapes. Advances in this burgeoning field provide important conceptual and methodological bases to decipher the historical role of natural selection and assess adaptation to past environmental variation. We suggest that a process called "temporal conditional neutrality" has taken place: some alleles favored in glacial environments become selectively neutral in modern-day conditions, whereas some alleles that had been neutral during glacial periods become under selection in modern environments. Building on this view, we present a new integrative framework for addressing the interplay of demographic and adaptive evolutionary responses to Quaternary climate dynamics, the research agenda initially envisioned by Davis and Shaw (2001) Science 292:673.

Multiscale landscape genomic models to detect signatures of selection in the alpine plant Biscutella laevigata

Plant species are known to adapt locally to their environment, particularly in mountainous areas where conditions can vary drastically over short distances. The climate of such landscapes being largely influenced by topography, using fine-scale models to evaluate environmental heterogeneity may help detecting adaptation to micro-habitats. Here, we applied a multiscale landscape genomic approach to detect evidence of local adaptation in the alpine plant Biscutella laevigata. The two gene pools identified, experiencing limited gene flow along a 1-km ridge, were different in regard to several habitat features derived from a very high resolution (VHR) digital elevation model (DEM). A correlative approach detected signatures of selection along environmental gradients such as altitude, wind exposure, and solar radiation, indicating adaptive pressures likely driven by fine-scale topography. Using a large panel of DEM-derived variables as ecologically relevant proxies, our results highlighted the critical role of spatial resolution. These high-resolution multiscale variables indeed indicate that the robustness of associations between genetic loci and environmental features depends on spatial parameters that are poorly documented. We argue that the scale issue is critical in landscape genomics and that multiscale ecological variables are key to improve our understanding of local adaptation in highly heterogeneous landscapes.

Genetic divergence and signatures of natural selection in marginal populations of a keystone, long-lived conifer, Eastern White Pine (Pinus strobus) from Northern Ontario

Marginal populations are expected to provide the frontiers for adaptation, evolution and range shifts of plant species under the anticipated climate change conditions. Marginal populations are predicted to show genetic divergence from central populations due to their isolation, and divergent natural selection and genetic drift operating therein. Marginal populations are also expected to have lower genetic diversity and effective population size (Ne) and higher genetic differentiation than central populations. We tested these hypotheses using eastern white pine (Pinus strobus) as a model for keystone, long-lived widely-distributed plants. All 614 eastern white pine trees, in a complete census of two populations each of marginal old-growth, central old-growth, and central second-growth, were genotyped at 11 microsatellite loci. The central populations had significantly higher allelic and genotypic diversity, latent genetic potential (LGP) and Ne than the marginal populations. However, heterozygosity and fixation index were similar between them. The marginal populations were genetically diverged from the central populations. Model testing suggested predominant north to south gene flow in the study area with curtailed gene flow to northern marginal populations. Signatures of natural selection were detected at three loci in the marginal populations; two showing divergent selection with directional change in allele frequencies, and one balancing selection. Contrary to the general belief, no significant differences were observed in genetic diversity, differentiation, LGP, and Ne between old-growth and second-growth populations. Our study provides information on the dynamics of migration, genetic drift and selection in central versus marginal populations of a keystone long-lived plant species and has broad evolutionary, conservation and adaptation significance.

Adaptive divergence along environmental gradients in a climate-change-sensitive mammal

In the face of predicted climate change, a broader understanding of biotic responses to varying environments has become increasingly important within the context of biodiversity conservation. Local adaptation is one potential option, yet remarkably few studies have harnessed genomic tools to evaluate the efficacy of this response within natural populations. Here, we show evidence of selection driving divergence of a climate-change-sensitive mammal, the American pika (Ochotona princeps), distributed along elevation gradients at its northern range margin in the Coast Mountains of British Columbia (BC), Canada. We employed amplified-fragment-length-polymorphism-based genomic scans to conduct genomewide searches for candidate loci among populations inhabiting varying environments from sea level to 1500 m. Using several independent approaches to outlier locus detection, we identified 68 candidate loci putatively under selection (out of a total 1509 screened), 15 of which displayed significant associations with environmental variables including annual precipitation and maximum summer temperature. These candidate loci may represent important targets for predicting pika responses to climate change and informing novel approaches to wildlife conservation in a changing world.

Historical connectivity, contemporary isolation and local adaptation in a widespread but discontinuously distributed species endemic to Taiwan, Rhododendron oldhamii (Ericaceae)

Elucidation of the evolutionary processes that constrain or facilitate adaptive divergence is a central goal in evolutionary biology, especially in non-model organisms. We tested whether changes in dynamics of gene flow (historical vs contemporary) caused population isolation and examined local adaptation in response to environmental selective forces in fragmented Rhododendron oldhamii populations. Variation in 26 expressed sequence tag-simple sequence repeat loci from 18 populations in Taiwan was investigated by examining patterns of genetic diversity, inbreeding, geographic structure, recent bottlenecks, and historical and contemporary gene flow. Selection associated with environmental variables was also examined. Bayesian clustering analysis revealed four regional population groups of north, central, south and southeast with significant genetic differentiation. Historical bottlenecks beginning 9168-13,092 years ago and ending 1584-3504 years ago were revealed by estimates using approximate Bayesian computation for all four regional samples analyzed. Recent migration within and across geographic regions was limited. However, major dispersal sources were found within geographic regions. Altitudinal clines of allelic frequencies of environmentally associated positively selected outliers were found, indicating adaptive divergence. Our results point to a transition from historical population connectivity toward contemporary population isolation and divergence on a regional scale. Spatial and temporal dispersal differences may have resulted in regional population divergence and local adaptation associated with environmental variables, which may have played roles as selective forces at a regional scale.

Use of qualitative environmental and phenotypic variables in the context of allele distribution models: detecting signatures of selection in the genome of Lake Victoria cichlids

When searching for loci possibly under selection in the genome, an alternative to population genetics theoretical models is to establish allele distribution models (ADM) for each locus to directly correlate allelic frequencies and environmental variables such as precipitation, temperature, or sun radiation. Such an approach implementing multiple logistic regression models in parallel was implemented within a computing program named MATSAM: . Recently, this application was improved in order to support qualitative environmental predictors as well as to permit the identification of associations between genomic variation and individual phenotypes, allowing the detection of loci involved in the genetic architecture of polymorphic characters. Here, we present the corresponding methodological developments and compare the results produced by software implementing population genetics theoretical models (DFDIST: and BAYESCAN: ) and ADM (MATSAM: ) in an empirical context to detect signatures of genomic divergence associated with speciation in Lake Victoria cichlid fishes.

Spatial genetic structure in natural populations of Phragmites australis in a mosaic of saline habitats in the Yellow River Delta, China

Determination of spatial genetic structure (SGS) in natural populations is important for both theoretical aspects of evolutionary genetics and their application in species conservation and ecological restoration. In this study, we examined genetic diversity within and among the natural populations of a cosmopolitan grass Phragmites australis (common reed) in the Yellow River Delta (YRD), China, where a mosaic of habitat patches varying in soil salinity was detected. We demonstrated that, despite their close geographic proximity, the common reed populations in the YRD significantly diverged at six microsatellite loci, exhibiting a strong association of genetic variation with habitat heterogeneity. Genetic distances among populations were best explained as a function of environmental difference, rather than geographical distance. Although the level of genetic divergence among populations was relatively low (F’_ST = 0.073), weak but significant genetic differentiation, as well as the concordance between ecological and genetic landscapes, suggests spatial structuring of genotypes in relation to patchy habitats. These findings not only provided insights into the population dynamics of common reed in changing environments, but also demonstrated the feasibility of using habitat patches in a mosaic landscape as test systems to identify appropriate genetic sources for ecological restoration.

AFLP genome scan to detect genetic structure and candidate loci under selection for local adaptation of the invasive weed Mikania micrantha

Why some species become successful invaders is an important issue in invasive biology. However, limited genomic resources make it very difficult for identifying candidate genes involved in invasiveness. Mikania micrantha H.B.K. (Asteraceae), one of the world's most invasive weeds, has adapted rapidly in response to novel environments since its introduction to southern China. In its genome, we expect to find outlier loci under selection for local adaptation, critical to dissecting the molecular mechanisms of invasiveness. An explorative amplified fragment length polymorphism (AFLP) genome scan was used to detect candidate loci under selection in 28 M. micrantha populations across its entire introduced range in southern China. We also estimated population genetic parameters, bottleneck signatures, and linkage disequilibrium. In binary characters, such as presence or absence of AFLP bands, if all four character combinations are present, it is referred to as a character incompatibility. Since character incompatibility is deemed to be rare in populations with extensive asexual reproduction, a character incompatibility analysis was also performed in order to infer the predominant mating system in the introduced M. micrantha populations. Out of 483 AFLP loci examined using stringent significance criteria, 14 highly credible outlier loci were identified by Dfdist and Bayescan. Moreover, remarkable genetic variation, multiple introductions, substantial bottlenecks and character compatibility were found to occur in M. micrantha. Thus local adaptation at the genome level indeed exists in M. micrantha, and may represent a major evolutionary mechanism of successful invasion. Interactions between genetic diversity, multiple introductions, and reproductive modes contribute to increase the capacity of adaptive evolution.

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.

AFLP genome scan to detect genetic structure and candidate loci under selection for local adaptation of the invasive weed Mikania micrantha

Why some species become successful invaders is an important issue in invasive biology. However, limited genomic resources make it very difficult for identifying candidate genes involved in invasiveness. Mikania micrantha H.B.K. (Asteraceae), one of the world's most invasive weeds, has adapted rapidly in response to novel environments since its introduction to southern China. In its genome, we expect to find outlier loci under selection for local adaptation, critical to dissecting the molecular mechanisms of invasiveness. An explorative amplified fragment length polymorphism (AFLP) genome scan was used to detect candidate loci under selection in 28 M. micrantha populations across its entire introduced range in southern China. We also estimated population genetic parameters, bottleneck signatures, and linkage disequilibrium. In binary characters, such as presence or absence of AFLP bands, if all four character combinations are present, it is referred to as a character incompatibility. Since character incompatibility is deemed to be rare in populations with extensive asexual reproduction, a character incompatibility analysis was also performed in order to infer the predominant mating system in the introduced M. micrantha populations. Out of 483 AFLP loci examined using stringent significance criteria, 14 highly credible outlier loci were identified by Dfdist and Bayescan. Moreover, remarkable genetic variation, multiple introductions, substantial bottlenecks and character compatibility were found to occur in M. micrantha. Thus local adaptation at the genome level indeed exists in M. micrantha, and may represent a major evolutionary mechanism of successful invasion. Interactions between genetic diversity, multiple introductions, and reproductive modes contribute to increase the capacity of adaptive evolution.

Evidence for evolutionary change associated with the recent range expansion of the British butterfly, Aricia agestis, in response to climate change

Poleward range expansions are widespread responses to recent climate change and are crucial for the future persistence of many species. However, evolutionary change in traits such as colonization history and habitat preference may also be necessary to track environmental change across a fragmented landscape. Understanding the likelihood and speed of such adaptive change is important in determining the rate of species extinction with ongoing climate change. We conducted an amplified fragment length polymorphism (AFLP)-based genome scan across the recently expanded UK range of the Brown Argus butterfly, Aricia agestis, and used outlier-based (DFDIST and BayeScan) and association-based (Isolation-By-Adaptation) statistical approaches to identify signatures of evolutionary change associated with range expansion and habitat use. We present evidence for (i) limited effects of range expansion on population genetic structure and (ii) strong signatures of selection at approximately 5% AFLP loci associated with both the poleward range expansion of A. agestis and differences in habitat use across long-established and recently colonized sites. Patterns of allele frequency variation at these candidate loci suggest that adaptation to new habitats at the range margin has involved selection on genetic variation in habitat use found across the long-established part of the range. Our results suggest that evolutionary change is likely to affect species' responses to climate change and that genetic variation in ecological traits across species' distributions should be maximized to facilitate range shifts across a fragmented landscape, particularly in species that show strong associations with particular habitats.

Population genetic diversity and structure of a naturally isolated plant species, Rhodiola dumulosa (Crassulaceae)

Aims

Rhodiola dumulosa (Crassulaceae) is a perennial diploid species found in high-montane areas. It is distributed in fragmented populations across northern, central and northwestern China. In this study, we aimed to (i) measure the genetic diversity of this species and that of its populations; (ii) describe the genetic structure of these populations across the entire distribution range in China; and (iii) evaluate the extent of gene flow among the naturally fragmented populations.

Methods

Samples from 1089 individuals within 35 populations of R. dumulosa were collected, covering as much of the entire distribution range of this species within China as possible. Population genetic diversity and structure were analyzed using AFLP molecular markers. Gene flow among populations was estimated according to the level of population differentiation.

Important Findings

The total genetic diversity of R. dumulosa was high but decreased with increasing altitude. Population-structure analysis indicated that the most closely related populations were geographically restricted and occurred in close proximity to each other. A significant isolation-by-distance pattern, caused by the naturally fragmented population distribution, was observed. At least two distinct gene pools were found in the 35 sampled populations, one composed of populations in northern China and the other composed of populations in central and northwestern China. The calculation of Nei's gene diversity index revealed that the genetic diversity in the northern China pool (0.1972) was lower than that in the central and northwestern China pool (0.2216). The populations were significantly isolated, and gene flow was restricted throughout the entire distribution. However, gene flow among populations on the same mountain appears to be unrestricted, as indicated by the weak genetic isolation among these populations.

Southward Pleistocene migration of Douglas-fir into Mexico: phylogeography, ecological niche modeling, and conservation of 'rear edge' populations

• Poleward Pleistocene plant migration has been an important process structuring modern temperate and boreal plant communities, but the contribution of equatorward migration remains poorly understood. Paleobotanical evidence suggests Miocene or Pleistocene origin for temperate 'sky island' plant taxa in Mexico. These 'rear edge' populations situated in a biodiversity hotspot may be an important reserve of genetic diversity in changing climates. • We used mtDNA sequences, cpDNA sequences and chloroplast microsatellites to test hypotheses of Miocene vs Pleistocene colonization of temperate Douglas-fir in Mexico, explore geographic patterns of molecular variation in relation to Pleistocene climate history using ecological niche models, and assess the taxonomic and conservation implications. • We found strong evidence for Pleistocene divergence of Douglas-fir in Mexico (958 thousand yr before present (ka) with the 90% highest posterior density interval ranging from 1.6 million yr before present (Ma) to 491 ka), consistent with the southward Pleistocene migration hypothesis. Genetic diversity was high and strongly partitioned among populations. Spatial patterns of molecular variation and ecological niche models suggest a complex late Pleistocene history involving periods of isolation and expansion along mountain corridors. • These results highlight the importance of southward Pleistocene migration in establishing modern high-diversity plant communities and provide critical insights into proposals to conserve the unique biodiversity of Mexican Douglas-fir and associated taxa.

Landscape genetics of plants

Landscape genetics is the amalgamation of landscape ecology and population genetics to help with understanding microevolutionary processes such as gene flow and adaptation. In this review, we examine why landscape genetics of plants lags behind that of animals, both in number of studies and consideration of landscape elements. The classical landscape distance/resistance approach to study gene flow is challenging in plants, whereas boundary detection and the assessment of contemporary gene flow are more feasible. By contrast, the new field of landscape genetics of adaptive genetic variation, establishing the relationship between adaptive genomic regions and environmental factors in natural populations, is prominent in plant studies. Landscape genetics is ideally suited to study processes such as migration and adaptation under global change.