Microrefugia and microclimate: Unraveling decoupling potential and resistance to heatwaves

Microrefugia, defined as small areas maintaining populations of species outside their range margins during environmental extremes, are increasingly recognized for their role in conserving species in the face of climate change. Understanding their microclimatic dynamics becomes crucial with global warming leading to severe temperature and precipitation changes. This study investigates the phenomenon of short-term climatic decoupling within microrefugia and its implications for plant persistence in the Mediterranean region of southeastern France. We focus on microrefugia's ability to climatically disconnect from macroclimatic trends, examining temperature and Vapor Pressure Deficit (VPD) dynamics in microrefugia, adjacent control plots, and weather stations. Our study encompasses both "normal" conditions and heatwave episodes to explore the role of microrefugia as thermal and moisture insulators during extreme events. Landscape attributes such as relative elevation, solar radiation, distance to streams, and vegetation height are investigated for their contribution to short-term decoupling. Our results demonstrate that microrefugia exhibit notable decoupling from macroclimatic trends. This effect is maintained during heatwaves, underscoring microrefugia's vital role in responding to climatic extremes. Importantly, microrefugia maintain lower VPD levels than their surroundings outside and during heatwaves, potentially mitigating water stress for plants. This study advances our understanding of microclimate dynamics within microrefugia and underscores their ecological importance for plant persistence in a changing climate. As heatwaves become more frequent and severe, our findings provide insights into the role of microrefugia in buffering but also decoupling against extreme climatic events and, more generally, against climate warming. This knowledge emphasizes the need to detect and protect existing microrefugia, as they can be integrated into conservation strategies and climate change adaptation plans.

Scots pine - panmixia and the elusive signal of genetic adaptation

Scots pine is the foundation species of diverse forested ecosystems across Eurasia and displays remarkable ecological breadth, occurring in environments ranging from temperate rainforests to arid tundra margins. Such expansive distributions can be favored by various demographic and adaptive processes and the interactions between them. To understand the impact of neutral and selective forces on genetic structure in Scots pine, we conducted range-wide population genetic analyses on 2321 trees from 202 populations using genotyping-by-sequencing, reconstructed the recent demography of the species and examined signals of genetic adaptation. We found a high and uniform genetic diversity across the entire range (global FST 0.048), no increased genetic load in expanding populations and minor impact of the last glacial maximum on historical population sizes. Genetic-environmental associations identified only a handful of single-nucleotide polymorphisms significantly linked to environmental gradients. The results suggest that extensive gene flow is predominantly responsible for the observed genetic patterns in Scots pine. The apparent missing signal of genetic adaptation is likely attributed to the intricate genetic architecture controlling adaptation to multi-dimensional environments. The panmixia metapopulation of Scots pine offers a good study system for further exploration into how genetic adaptation and plasticity evolve under gene flow and changing environment.

Cliff-edge forests: Xerothermic hotspots of local biodiversity and models for future climate change

Cliffs are remarkable environments that enable the existence of microclimates. These small, isolated sites, decoupled from the regional macroclimate, play a significant role in maintaining species biodiversity, particularly in topographically homogeneous landscapes. Our study investigated the microclimate of south-exposed forests situated at the edge of sandstone cliffs in the western part of the North Alpine Foreland Basin in Switzerland and its role in local forest community composition. Using direct measurements from data loggers, as well as vegetation analyses, it was possible to quantify the microclimate of the cliff-edge forests and compare it with that of the surrounding forests. Our results highlighted the significant xerothermic and more variable nature of the cliff-edge forest microclimate, with a mean soil temperature up to 3.72°C warmer in the summer, higher annual (+28%) and daily (+250%) amplitudes of soil temperature, which frequently expose vegetation to extreme temperatures, and an 83% higher soil drying rate. These differences have a distinct influence on forest communities: cliff-edge forests are significantly different from surrounding forests. The site particularities of cliff edges support the presence of locally rare species and forest types, particularly of Scots pine. Cliff edges must therefore be considered microrefugia with a high conservation value for both xerothermic species and flora adapted to more continental climates. Moreover, the microclimate of cliff-edge forests could resemble the future climate in many ways. We argue that these small areas, which are already experiencing the future climate, can be seen as natural laboratories to better answer the following question: what will our forests look like in a few decades with accelerated climate change?

Molecular characterisation of Pinus sylvestris (L.) in Ireland at the western limit of the species distribution

BACKGROUND

Scots pine (Pinus sylvestris L.) underwent significant population declines across much of northwest Europe during the mid-to-late Holocene and was thought to have become extirpated in Ireland from about 400 AD. However, most extant populations are plantations reintroduced from Scotland. Others are naturalised therefrom and one in Western Ireland is a putative relict. In this paper, Scots pine in Ireland are genetically described for the first time.

RESULTS

Using two mitochondrial (mtDNA) loci, eight chloroplast (cpSSR) and 18 nuclear (nSSR) loci, the genetic composition and diversity of 19 Irish Scots pine populations is described and compared to other European populations. All trees sampled in Ireland were fixed for mitotype a, which is the most common across northwest Europe. By contrast, cpSSR (HCP = 0.967) and nSSR (He = 0.540) variation was high, and comparable with estimates for other regions across the species range. Differentiation at both sets of loci were similarly low (cpSSR FST = 0.019; nSSR FST = 0.018), but populations from continental Europe were significantly differentiated from all Irish populations based on nSSR variation.

CONCLUSIONS

All Irish Scots pine are likely part of a common Irish-Scottish gene pool which diverged from continental Scots pine following post-glacial recolonisation. A high genetic diversity and an absence of evidence of inbreeding suggests the regional decline of Scots pine did not critically reduce allelic variation. The post-glacial relationship between Irish and Scottish pine is discussed, and a suggestion from recent palaeoecological work that reintroduced Scots pine be managed as a native species is now further supported by genetic data.

Phylogeography of Ramalina farinacea (Lichenized Fungi, Ascomycota) in the Mediterranean Basin, Europe, and Macaronesia

Ramalina farinacea is an epiphytic lichen-forming fungus with a broad geographic distribution, especially in the Northern Hemisphere. In the eighties of the last century, it was hypothesized that R. farinacea had originated in the Macaronesian–Mediterranean region, with the Canary Islands as its probable southernmost limit, and thereafter it would have increased its distribution area. In order to explore the phylogeography of this emblematic lichen, we analyzed 120 thalli of R. farinacea collected in 38 localities distributed in temperate and boreal Europe, the Western Mediterranean Basin, and several Macaronesian archipelagos in the Atlantic Ocean. Data from two nuclear markers (nrITS and uid70) of the mycobiont were obtained to calculate genetic diversity indices to infer the phylogenies and haplotype networks and to investigate population structure. In addition, dating analysis was conducted to provide a valuable hypothesis of the timing of the origin and diversification of R. farinacea and its close allies. Our results highlight that phylogenetic species circumscription in the “Ramalina farinacea group” is complex and suggests that incomplete lineage sorting is at the base of conflicting phylogenetic signals. The existence of a high number of haplotypes restricted to the Macaronesian region, together with the diversification of R. farinacea in the Pleistocene, suggests that this species and its closest relatives originated during relatively recent geological times and then expanded its range to higher latitudes. However, our data cannot rule out whether the species originated from the Macaronesian archipelagos exclusively or also from the Mediterranean Basin. In conclusion, the present work provides a valuable biogeographical hypothesis for disentangling the evolution of this epiphytic lichen in space and time.

Bridging the gap between microclimate and microrefugia: A bottom-up approach reveals strong climatic and biological offsets

In the context of global warming, a clear understanding of microrefugia-microsites enabling the survival of species populations outside their main range limits-is crucial. Several studies have identified forcing factors that are thought to favor the existence of microrefugia. However, there is a lack of evidence to conclude whether, and to what extent, the climate encountered within existing microrefugia differs from the surrounding climate. To investigate this, we adopt a "bottom-up" approach, linking marginal disconnected populations to microclimate. We used the southernmost disconnected and abyssal populations of the circumboreal herbaceous plant Oxalis acetosella in Southern France to study whether populations in sites matching the definition of "microrefugia" occur in particularly favorable climatic conditions compared to neighboring control plots located at distances of between 50 to 100 m. Temperatures were recorded in putative microrefugia and in neighboring plots for approximately 2 years to quantify their thermal offsets. Vascular plant inventories were carried out to test whether plant communities also reflect microclimatic offsets. We found that current microclimatic dynamics are genuinely at stake in microrefugia. Microrefugia climates are systematically colder compared to those found in neighboring control plots. This pattern was more noticeable during the summer months. Abyssal populations showed stronger offsets compared to neighboring plots than the putative microrefugia occurring at higher altitudes. Plant communities demonstrate this strong spatial climatic variability, even at such a microscale approach, as species compositions systematically differed between the two plots, with species more adapted to colder and moister conditions in microrefugia compared to the surrounding area.

Environmental DNA as an emerging tool in botanical research

Over the past quarter century, environmental DNA (eDNA) has been ascendant as a tool to detect, measure, and monitor biodiversity (species and communities), as a means of elucidating biological interaction networks, and as a window into understanding past patterns of biodiversity. However, only recently has the potential of eDNA been realized in the botanical world. Here we synthesize the state of eDNA applications in botanical systems with emphases on aquatic, ancient, contemporary sediment, and airborne systems, and focusing on both single-species approaches and multispecies community metabarcoding. Further, we describe how abiotic and biotic factors, taxonomic resolution, primer choice, spatiotemporal scales, and relative abundance influence the utilization and interpretation of airborne eDNA results. Lastly, we explore several areas and opportunities for further development of eDNA tools for plants, advancing our knowledge and understanding of the efficacy, utility, and cost-effectiveness, and ultimately facilitating increased adoption of eDNA analyses in botanical systems.

MetaDamage tool: Examining post-mortem damage in sedaDNA on a metagenomic scale

The use of metagenomic datasets to support ancient sedimentary DNA (sedaDNA) for paleoecological reconstruction has been demonstrated to be a powerful tool to understand multi-organism responses to climatic shifts and events. Authentication remains integral to the ancient DNA discipline, and this extends to sedaDNA analysis. Furthermore, distinguishing authentic sedaDNA from contamination or modern material also allows for a better understanding of broader questions in sedaDNA research, such as formation processes, source and catchment, and post-depositional processes. Existing tools for the detection of damage signals are designed for single-taxon input, require a priori organism specification, and require a significant number of input sequences to establish a signal. It is therefore often difficult to identify an established cytosine deamination rate consistent with ancient DNA across a sediment sample. In this study, we present MetaDamage, a tool that examines cytosine deamination on a metagenomic (all organisms) scale for multiple previously undetermined taxa and can produce a damage profile based on a few hundred reads. We outline the development and testing of the MetaDamage tool using both authentic sedaDNA sequences and simulated data to demonstrate the resolution in which MetaDamage can identify deamination levels consistent with the presence of ancient DNA. The MetaDamage tool offers a method for the initial assessment of the presence of sedaDNA and a better understanding of key questions of preservation for paleoecological reconstruction.

Possible northern persistence of Siebold's beech, Fagus crenata, at its northernmost distribution limit on an island in Japan Sea: Okushiri Island, Hokkaido

Siebold's beech, Fagus crenata, is widely distributed across the Japanese Archipelago and islands in Japan Sea. Similar to the northern limit of the geographical distribution of F. crenata on the mainland of Hokkaido, the northern limit of the distribution of F. crenata on islands in the Japan Sea is observed on Okushiri Island (ca 42°N). To understand the genetic relationships of F. crenata on Okushiri Island, we examined chloroplast (cp) DNA haplotypes and 11 nuclear microsatellite (SSR) loci among 1,838 individuals from 44 populations from Okushiri Island, mainland Hokkaido, and the northern part of the Tohoku region on Honshu Island. We identified 2 cpDNA haplotypes, which represent not only populations on the Japan Sea coast but also those on the Pacific coast and this suggested the Okushiri Island populations might not be formed by single colonization. Genetic diversity of the Okushiri Island populations of nuclear SSR was not lower than the mainland and the STRUCTURE analysis revealed the Okushiri Island individuals were admixed between Hokkaido and Tohoku clusters. Approximate Bayesian computation inferred that divergence between Tohoku and Hokkaido, and admixture between two populations which generated Okushiri populations occurred before the last glacial maximum (LGM), that is, 7,890 (95% hyper probability density (HPD): 3,420 - 9,910) and 3,870 (95% HPD: 431- 8,540) generations ago, respectively. These inferences were well supported by a geological history which suggested an isolation of Okushiri Island from Hokkaido started prior to the Middle Pleistocene. We discuss the possible persistence of F. crenata during the last glacial maximum on northern islands in the Japan Sea such as Okushiri Island.

Teasing apart the joint effect of demography and natural selection in the birth of a contact zone

Vast population movements induced by recurrent climatic cycles have shaped the genetic structure of plant species. During glacial periods species were confined to low-latitude refugia from which they recolonized higher latitudes as the climate improved. This multipronged recolonization led to many lineages that later met and formed large contact zones. We utilize genomic data from 5000 Picea abies trees to test for the presence of natural selection during recolonization and establishment of a contact zone in Scandinavia. Scandinavian P. abies is today made up of a southern genetic cluster originating from the Baltics, and a northern one originating from Northern Russia. The contact zone delineating them closely matches the limit between two major climatic regions. We show that natural selection contributed to its establishment and maintenance. First, an isolation-with-migration model with genome-wide linked selection fits the data better than a purely neutral one. Second, many loci show signatures of selection or are associated with environmental variables. These loci, regrouped in clusters on chromosomes, are often related to phenology. Altogether, our results illustrate how climatic cycles, recolonization and selection can establish strong local adaptation along contact zones and affect the genetic architecture of adaptive traits.

The untapped potential of macrofossils in ancient plant DNA research

The rapid development of ancient DNA analysis in the last decades has induced a paradigm shift in ecology and evolution. Driven by a combination of breakthroughs in DNA isolation techniques, high-throughput sequencing, and bioinformatics, ancient genome-scale data for a rapidly growing variety of taxa are now available, allowing researchers to directly observe demographic and evolutionary processes over time. However, the vast majority of paleogenomic studies still focus on human or animal remains. In this article, we make the case for a vast untapped resource of ancient plant material that is ideally suited for paleogenomic analyses: plant remains, such as needles, leaves, wood, seeds, or fruits, that are deposited in natural archives, such as lake sediments, permafrost, or even ice caves. Such plant remains are commonly found in large numbers and in stratigraphic sequence through time and have so far been used primarily to reconstruct past local species presences and abundances. However, they are also unique repositories of genetic information with the potential to revolutionize the fields of ecology and evolution by directly studying microevolutionary processes over time. Here, we give an overview of the current state-of-the-art, address important challenges, and highlight new research avenues to inspire future research.

Long-term preservation of biomolecules in lake sediments: potential importance of physical shielding by recalcitrant cell walls

Even though lake sediments are globally important organic carbon (OC) sinks, the controls on long-term OC storage in these sediments are unclear. Using a multiproxy approach, we investigate changes in diatom, green algae, and vascular plant biomolecules in sedimentary records from the past centuries across five temperate lakes with different trophic histories. Despite past increases in the input and burial of OC in sediments of eutrophic lakes, biomolecule quantities in sediments of all lakes are primarily controlled by postburial microbial degradation over the time scales studied. We, moreover, observe major differences in biomolecule degradation patterns across diatoms, green algae, and vascular plants. Degradation rates of labile diatom DNA exceed those of chemically more resistant diatom lipids, suggesting that chemical reactivity mainly controls diatom biomolecule degradation rates in the lakes studied. By contrast, degradation rates of green algal and vascular plant DNA are significantly lower than those of diatom DNA, and in a similar range as corresponding, much less reactive lipid biomarkers and structural macromolecules, including lignin. We propose that physical shielding by degradation-resistant cell wall components, such as algaenan in green algae and lignin in vascular plants, contributes to the long-term preservation of labile biomolecules in both groups and significantly influences the long-term burial of OC in lake sediments.

Norway spruce postglacial recolonization of Fennoscandia

Contrasting theories exist regarding how Norway spruce (Picea abies) recolonized Fennoscandia after the last glaciation and both early Holocene establishments from western microrefugia and late Holocene colonization from the east have been postulated. Here, we show that Norway spruce was present in southern Fennoscandia as early as 14.7 ± 0.1 cal. kyr BP and that the millennia-old clonal spruce trees present today in central Sweden likely arrived with an early Holocene migration from the east. Our findings are based on ancient sedimentary DNA from multiple European sites (N = 15) combined with nuclear and mitochondrial DNA analysis of ancient clonal (N = 135) and contemporary spruce forest trees (N = 129) from central Sweden. Our other findings imply that Norway spruce was present shortly after deglaciation at the margins of the Scandinavian Ice Sheet, and support previously disputed finds of pollen in southern Sweden claiming spruce establishment during the Lateglacial.

Contrasting theories exist regarding how Norway spruce recolonized Fennoscandia after the last glaciation. Here, the authors provide evidences from sedimentary ancient DNA and modern population genomics to support that Norway spruce was present in southern Fennoscandia shortly after deglaciation and the early Holocene migration from the east.

Non-linear genetic diversity and notable population differentiation caused by low gene flow of bermudagrass [Cynodon dactylon (L.) Pers.] along longitude gradients

Background

Environmental variation related to ecological habitat is the main driver of plant adaptive divergence. Longitude plays an important role in the formation of plant population structure, indicating that environmental differentiation can significantly shape population structure.

Methods

Genetic diversity and population genetic structure were estimated using 105 expressed sequence tag-derived simple sequence repeat (EST-SSR) loci. A total of 249 C. dactylon (L.) Pers. (common bermudagrass) individuals were sampled from 13 geographic sites along the longitude (105°57′34″–119°27′06″E).

Results

There was no obvious linear trend of intra-population genetic diversity along longitude and the intra-population genetic diversity was not related to climate in this study. Low gene flow (Nm = 0.7701) meant a rich genetic differentiation among populations of C. dactylon along longitude gradients. Significantly positive Mantel correlation (r = 0.438, P = 0.001) was found between genetic distance and geographical interval while no significant partial Mantel correlation after controlling the effect of mean annual precipitation, which indicated geographic distance correlated with mean annual precipitation affect genetic distance. The genetic diversity of C. dactylon with higher ploidy level was higher than that with lower ploidy level and groups of individuals with higher ploidy level were separated further away by genetic distance from the lower ploidy levels. Understanding the different genetic bases of local adaptation comparatively between latitude and longitude is one of the core findings in the adaptive evolution of plants.

Using target capture to address conservation challenges: Population-level tracking of a globally-traded herbal medicine

The promotion of responsible and sustainable trade in biological resources is widely proposed as one solution to mitigate current high levels of global biodiversity loss. Various molecular identification methods have been proposed as appropriate tools for monitoring global supply chains of commercialized animals and plants. Here, we demonstrate the efficacy of target capture genomic barcoding in identifying and establishing the geographic origin of samples traded as Anacyclus pyrethrum, a medicinal plant assessed as globally vulnerable in the IUCN Red List of Threatened Species. Samples collected from national and international supply chains were identified through target capture sequencing of 443 low-copy nuclear makers and compared to results derived from genome skimming of plastome and DNA barcoding of standard plastid regions and ITS. Both target capture and genome skimming provided approximately 3.4 million reads per sample, but target capture largely outperformed standard plant barcodes and entire plastid genome sequences. We were able to discern the geographical origin of Anacyclus samples collected in Moroccan, Indian and Sri Lankan markets, differentiating between plant materials originally harvested from diverse populations in Algeria and Morocco. Dropping costs of analysing samples enables the potential of target capture to routinely identify commercialized plant species and determine their geographic origin. It promises to play an important role in monitoring and regulation of plant species in trade, supporting biodiversity conservation efforts, and in ensuring that plant products are unadulterated, contributing to consumer protection.

Phylogeographic reconstructions can be biased by ancestral shared alleles: The case of the polymorphic lichen Bryoria fuscescens in Europe and North Africa

Large phylogeographic studies on lichens are scarce, and none involves a single species within which different lineages show fixed alternative dispersal strategies. We investigated Bryoria fuscescens (including B. capillaris) in Europe and western North Africa by phenotypically characterizing 1400 specimens from 64 populations and genotyping them with 14 microsatellites. We studied population structure and genetic diversity at the local and continental scales, discussed the post-glacial phylogeography, and compared dispersal capacities of phenotypes with and without soralia. Our main hypothesis is that the estimated phylogeography, migration routes, and dispersal capacities may be strongly biased by ancestral shared alleles. Scandinavia is genetically the richest area, followed by the Iberian Peninsula, the Carpathians, and the Alps. Three gene pools were detected: two partially linked to phenotypic characteristics, and the third one genetically related to the American sister species B. pseudofuscescens. The comparison of one gene pool producing soredia and one not, suggested both as panmictic, with similar levels of isolation by distance (IBD). The migration routes were estimated to span from north to south, in disagreement with the assessed glacial refugia. The presence of ancestral shared alleles in distant populations can explain the similar IBD levels found in both gene pools while producing a false signal of panmixia, and also biasing the phylogeographic reconstruction. The incomplete lineage sorting recorded for DNA sequence loci also supports this hypothesis. Consequently, the high diversity in Scandinavia may rather come from recent immigration into northern populations than from an in situ diversification. Similar patterns of ancestral shared polymorphism may bias the phylogeographical reconstruction of other lichen species.

The Origin of Isolated Populations of the Mountain Weevil, Liparus glabrirostris-The Flagship Species for Riparian Habitats

Liparus glabrirostris is one of the largest European weevils, and it has been recently proposed as the flagship species for threatened riparian habitats in the mountains. For effective conservation of its populations (and associated habitats), it is crucial to understand its history, contemporary distribution, genetic diversity and predict changes in the range, including its highly isolated populations on the Baltic coast. Here, we examined numerous populations of L. glabrirostris across almost the entire species range using phylogeography and species niche modeling (SNM) approaches. Analyses of mtDNA and nucDNA markers revealed the existence of 2 major mitochondrial lineages generally separated between 1) the Alpine region and 2) the Bohemian Massif, the Carpathians, and the Baltic coast areas. Genetic diversity in nuclear genes was more complicated with no clear division between populations. The origin of Baltic populations was derived from the Carpathians, but there were probably multiple expansion events to the north. SNM suggested the existence of glacial refugia for L. glabrirostris, mainly in the Alps and the Southern Carpathians. Current predictions of species range were found to be generally congruent with zoogeographic data; however, the Baltic coast was not really supported as a suitable area for L. glabrirostris. An important prediction of future distribution (2050-2070 CE) suggests a shrinkage of the L. glabrirostris range and extinction of some of its populations (particularly those isolated on lower altitudes). Based on the aforementioned data, proposals for the protection of this species are proposed, including the designation of several evolutionary units of conservation importance.

Environmental palaeogenomic reconstruction of an Ice Age algal population

Palaeogenomics has greatly increased our knowledge of past evolutionary and ecological change, but has been restricted to the study of species that preserve either as or within fossils. Here we show the potential of shotgun metagenomics to reveal population genomic information for a taxon that does not preserve in the body fossil record, the algae Nannochloropsis. We shotgun sequenced two lake sediment samples dated to the Last Glacial Maximum and reconstructed full chloroplast and mitochondrial genomes to explore within-lake population genomic variation. This revealed two major haplogroups for each organellar genome, which could be assigned to known varieties of N. limnetica, although we show that at least three haplotypes were present using our minimum haplotype diversity estimation method. These approaches demonstrate the utility of lake sedimentary ancient DNA (sedaDNA) for population genomic analysis, thereby opening the door to environmental palaeogenomics, which will unlock the full potential of sedaDNA.

Lammers et al. use sedimentary ancient DNA to reconstruct palaeogenomes of Nannochloropsis. This study demonstrates the value of sedaDNA for palaeogenomic reconstructions and population genomic analysis.

Lake Sedimentary DNA Research on Past Terrestrial and Aquatic Biodiversity: Overview and Recommendations

The use of lake sedimentary DNA to track the long-term changes in both terrestrial and aquatic biota is a rapidly advancing field in paleoecological research. Although largely applied nowadays, knowledge gaps remain in this field and there is therefore still research to be conducted to ensure the reliability of the sedimentary DNA signal. Building on the most recent literature and seven original case studies, we synthesize the state-of-the-art analytical procedures for effective sampling, extraction, amplification, quantification and/or generation of DNA inventories from sedimentary ancient DNA (sedaDNA) via high-throughput sequencing technologies. We provide recommendations based on current knowledge and best practises.

Divergent patterns between phenotypic and genetic variation in Scots pine

In boreal forests, autumn frost tolerance in seedlings is a critical fitness component because it determines survival rates during regeneration. To understand the forces that drive local adaptation in this trait, we conducted freezing tests in a common garden setting for 54 Pinus sylvestris (Scots pine) populations (>5000 seedlings) collected across Scandinavia into western Russia, and genotyped 24 of these populations (>900 seedlings) at >10 000 SNPs. Variation in cold hardiness among populations, as measured by Q_ST , was above 80% and followed a distinct cline along latitude and longitude, demonstrating significant adaptation to climate at origin. In contrast, the genetic differentiation was very weak (mean F_ST 0.37%). Despite even allele frequency distribution in the vast majority of SNPs among all populations, a few rare alleles appeared at very high or at fixation in marginal populations restricted to northwestern Fennoscandia. Genotype-environment associations showed that climate variables explained 2.9% of the genetic differentiation, while genotype-phenotype associations revealed a high marker-estimated heritability of frost hardiness of 0.56, but identified no major loci. Very extensive gene flow, strong local adaptation, and signals of complex demographic history across markers are interesting topics of forthcoming studies on this species to better clarify signatures of selection and demography.

Demography and Natural Selection Have Shaped Genetic Variation in the Widely Distributed Conifer Norway Spruce (Picea abies)

Under the neutral theory, species with larger effective population size are expected to harbor higher genetic diversity. However, across a wide variety of organisms, the range of genetic diversity is orders of magnitude more narrow than the range of effective population size. This observation has become known as Lewontin’s paradox and although aspects of this phenomenon have been extensively studied, the underlying causes for the paradox remain unclear. Norway spruce (Picea abies) is a widely distributed conifer species across the northern hemisphere, and it consequently plays a major role in European forestry. Here, we use whole-genome resequencing data from 35 individuals to perform population genomic analyses in P. abies in an effort to understand what drives genome-wide patterns of variation in this species. Despite having a very wide geographic distribution and an corresponding enormous current population size, our analyses find that genetic diversity of P. abies is low across a number of populations (π = 0.0049 in Central-Europe, π = 0.0063 in Sweden-Norway, π = 0.0063 in Finland). To assess the reasons for the low levels of genetic diversity, we infer the demographic history of the species and find that it is characterized by several reoccurring bottlenecks with concomitant decreases in effective population size can, at least partly, provide an explanation for low polymorphism we observe in P. abies. Further analyses suggest that recurrent natural selection, both purifying and positive selection, can also contribute to the loss of genetic diversity in Norway spruce by reducing genetic diversity at linked sites. Finally, the overall low mutation rates seen in conifers can also help explain the low genetic diversity maintained in Norway spruce.

High genomic diversity in the bank vole at the northern apex of a range expansion: The role of multiple colonizations and end-glacial refugia

The history of repeated northern glacial cycling and southern climatic stability has long dominated explanations for how genetic diversity is distributed within temperate species in Eurasia and North America. However, growing evidence indicates the importance of cryptic refugia for northern colonization dynamics. An important geographic region to assess this is Fennoscandia, where recolonization at the end of the last glaciation was restricted to specific routes and temporal windows. We used genomic data to analyse genetic diversity and colonization history of the bank vole (Myodes glareolus) throughout Europe (>800 samples) with Fennoscandia as the northern apex. We inferred that bank voles colonized Fennoscandia multiple times by two different routes; with three separate colonizations via a southern land-bridge route deriving from a "Carpathian" glacial refugium and one via a north-eastern route from an "Eastern" glacial refugium near the Ural Mountains. Clustering of genome-wide SNPs revealed high diversity in Fennoscandia, with eight genomic clusters: three of Carpathian origin and five Eastern. Time estimates revealed that the first of the Carpathian colonizations occurred before the Younger Dryas (YD), meaning that the first colonists survived the YD in Fennoscandia. Results also indicated that introgression between bank and northern red-backed voles (Myodes rutilus) took place in Fennoscandia just after end-glacial colonization. Therefore, multiple colonizations from the same and different cryptic refugia, temporal and spatial separations and interspecific introgression have shaped bank vole genetic variability in Fennoscandia. Together, these processes drive high genetic diversity at the apex of the northern expansion in this emerging model species.

Ancient DNA analysis of food remains in human dental calculus from the Edo period, Japan

Although there are many methods for reconstructing diets of the past, detailed taxon identification is still challenging, and most plants hardly remain at a site. In this study, we applied DNA metabarcoding to dental calculus of premodern Japan for the taxonomic identification of food items. DNA was extracted from 13 human dental calculi from the Unko-in site (18th–19th century) of the Edo period, Japan. Polymerase chain reaction (PCR) and sequencing were performed using a primer set specific to the genus Oryza because rice (Oryza sativa) was a staple food and this was the only member of this genus present in Japan at that time. DNA metabarcoding targeting plants, animals (meat and fish), and fungi were also carried out to investigate dietary diversity. We detected amplified products of the genus Oryza from more than half of the samples using PCR and Sanger sequencing. DNA metabarcoding enabled us to identify taxa of plants and fungi, although taxa of animals were not detected, except human. Most of the plant taxonomic groups (family/genus level) are present in Japan and include candidate species consumed as food at that time, as confirmed by historical literature. The other groups featured in the lifestyle of Edo people, such as for medicinal purposes and tobacco. The results indicate that plant DNA analysis from calculus provides information about food diversity and lifestyle habits from the past and can complement other analytical methods such as microparticle analysis and stable isotope analysis.

Quaternary DNA: A Multidisciplinary Research Field

The purpose of this Milankovitch review is to explain the significance of Quaternary DNA studies and the importance of the recent methodological advances that have enabled the study of late Quaternary remains in more detail, and the testing of new assumptions in evolutionary biology and phylogeography to reconstruct the past. The topic is wide, and this review is not intended to be an exhaustive account of all the aDNA work performed in the last three decades on late-Quaternary remains. Instead, it is a selection of relevant studies aimed at illustrating how aDNA has been used to reconstruct not only environments of the past, but also the history of many species including our own.

Chloroplast and mitochondrial genetic variation of larches at the Siberian tundra-taiga ecotone revealed by de novo assembly

Larix populations at the tundra-taiga ecotone in northern Siberia are highly under-represented in population genetic studies, possibly due to the remoteness of these regions that can only be accessed at extraordinary expense. The genetic signatures of populations in these boundary regions are therefore largely unknown. We aim to generate organelle reference genomes for the detection of single nucleotide polymorphisms (SNPs) that can be used for paleogenetic studies. We present 19 complete chloroplast genomes and mitochondrial genomic sequences of larches from the southern lowlands of the Taymyr Peninsula (northernmost range of Larix gmelinii (Rupr.) Kuzen.), the lower Omoloy River, and the lower Kolyma River (both in the range of Larix cajanderi Mayr). The genomic data reveal 84 chloroplast SNPs and 213 putatively mitochondrial SNPs. Parsimony-based chloroplast haplotype networks show no spatial structure of individuals from different geographic origins, while the mitochondrial haplotype network shows at least a slight spatial structure with haplotypes from the Omoloy and Kolyma populations being more closely related to each other than to most of the haplotypes from the Taymyr populations. Whole genome alignments with publicly available complete chloroplast genomes of different Larix species show that among official plant barcodes only the rcbL gene contains sufficient polymorphisms, but has to be sequenced completely to distinguish the different provenances. We provide 8 novel mitochondrial SNPs that are putatively diagnostic for the separation of L. gmelinii and L. cajanderi, while 4 chloroplast SNPs have the potential to distinguish the L. gmelinii/L. cajanderi group from other Larix species. Our organelle references can be used for a targeted primer and probe design allowing the generation of short amplicons. This is particularly important with regard to future investigations of, for example, the biogeographic history of Larix by screening ancient sedimentary DNA of Larix.

Genetic resources of relict populations of Pinus sylvestris (L.) in Western Carpathians assessed by chloroplast microsatellites

The relict character of the Scots pine populations from the Tatra Mts. (Western Carpathians) was formed on scarcely accessible, ecologically extreme habitats, in areas which were not under human activity. An exhaustive genetic analysis of the Scots pine populations from the refugial locations in the Tatra Mts. (Poland) had not yet been. In this study, we characterize the genetic variation and differentiation of the relict Scots pine populations from the Tatra Mts., to provide information on their genetic resources and the conservation implications. Eight paternally inherited chloroplast microsatellite loci were used to investigate the genetic structure of 4 relict populations, which were compared with the natural populations from the Baltic region (5 populations from Lithuania, Latvia and Estonia) and to a northern range from Russia (5 populations from the Kola Peninsula and Arkhangelsk Oblast). The four relict pine populations were characterized by a lower level of genetic variation when compared to the Baltic and northern populations, which was expressed in lower values of particular genetic parameters: numbers of alleles and haplotypes, haplotype genetic diversity, haplotype differentiation, and mean genetic distance between individuals in the population. Our results revealed a very high and significant genetic differentiation between all the analyzed populations, as well as between the three analyzed regions (PhiPT = 8%). Furthermore, the relict populations within a close geographic location showed higher differentiation (PhiPT = 5%) compared to the northern Russian and Baltic populations that were separated by tens and hundreds of kilometers. The relict populations of P. sylvestris from the Tatra Mts. deserve conservation efforts due to the fact that they are an important component if the species-poor Tatra forests.

Eyes on the future - evidence for trade-offs between growth, storage and defense in Norway spruce

Carbon (C) allocation plays a central role in tree responses to environmental changes. Yet, fundamental questions remain about how trees allocate C to different sinks, for example, growth vs storage and defense. In order to elucidate allocation priorities, we manipulated the whole-tree C balance by modifying atmospheric CO₂ concentrations [CO₂ ] to create two distinct gradients of declining C availability, and compared how C was allocated among fluxes (respiration and volatile monoterpenes) and biomass C pools (total biomass, nonstructural carbohydrates (NSC) and secondary metabolites (SM)) in well-watered Norway spruce (Picea abies) saplings. Continuous isotope labelling was used to trace the fate of newly-assimilated C. Reducing [CO₂ ] to 120 ppm caused an aboveground C compensation point (i.e. net C balance was zero) and resulted in decreases in growth and respiration. By contrast, soluble sugars and SM remained relatively constant in aboveground young organs and were partially maintained with a constant allocation of newly-assimilated C, even at expense of root death from C exhaustion. We conclude that spruce trees have a conservative allocation strategy under source limitation: growth and respiration can be downregulated to maintain 'operational' concentrations of NSC while investing newly-assimilated C into future survival by producing SM.

Postglacial colonization history reflects in the genetic structure of natural populations of Festuca rubra in Europe

We conducted a large-scale population genetic survey of genetic diversity of the host grass Festuca rubra s.l., which fitness can be highly dependent on its symbiotic fungus Epichloë festucae, to evaluate genetic variation and population structure across the European range. The 27 studied populations have previously been found to differ in frequencies of occurrence of the symbiotic fungus E. festucae and ploidy levels. As predicted, we found decreased genetic diversity in previously glaciated areas in comparison with nonglaciated regions and discovered three major maternal genetic groups: southern, northeastern, and northwestern Europe. Interestingly, host populations from Greenland were genetically similar to those from the Faroe Islands and Iceland, suggesting gene flow also between those areas. The level of variation among populations within regions is evidently highly dependent on the postglacial colonization history, in particular on the number of independent long-distance seed colonization events. Yet, also anthropogenic effects may have affected the population structure in F. rubra. We did not observe higher fungal infection rates in grass populations with lower levels of genetic variability. In fact, the fungal infection rates of E. festucae in relation to genetic variability of the host populations varied widely among geographical areas, which indicate differences in population histories due to colonization events and possible costs of systemic fungi in harsh environmental conditions. We found that the plants of different ploidy levels are genetically closely related within geographic areas indicating independent formation of polyploids in different maternal lineages.

Paleogenomics: reconstruction of plant evolutionary trajectories from modern and ancient DNA

How contemporary plant genomes originated and evolved is a fascinating question. One approach uses reference genomes from extant species to reconstruct the sequence and structure of their common ancestors over deep timescales. A second approach focuses on the direct identification of genomic changes at a shorter timescale by sequencing ancient DNA preserved in subfossil remains. Merged within the nascent field of paleogenomics, these complementary approaches provide insights into the evolutionary forces that shaped the organization and regulation of modern genomes and open novel perspectives in fostering genetic gain in breeding programs and establishing tools to predict future population changes in response to anthropogenic pressure and global warming.

Temporal and spatial patterns of mitochondrial haplotype and species distributions in Siberian larches inferred from ancient environmental DNA and modeling

Changes in species' distributions are classically projected based on their climate envelopes. For Siberian forests, which have a tremendous significance for vegetation-climate feedbacks, this implies future shifts of each of the forest-forming larch (Larix) species to the north-east. However, in addition to abiotic factors, reliable projections must assess the role of historical biogeography and biotic interactions. Here, we use sedimentary ancient DNA and individual-based modelling to investigate the distribution of larch species and mitochondrial haplotypes through space and time across the treeline ecotone on the southern Taymyr peninsula, which at the same time presents a boundary area of two larch species. We find spatial and temporal patterns, which suggest that forest density is the most influential driver determining the precise distribution of species and mitochondrial haplotypes. This suggests a strong influence of competition on the species' range shifts. These findings imply possible climate change outcomes that are directly opposed to projections based purely on climate envelopes. Investigations of such fine-scale processes of biodiversity change through time are possible using paleoenvironmental DNA, which is available much more readily than visible fossils and can provide information at a level of resolution that is not reached in classical palaeoecology.

Genetic analyses of brown hare (Lepus europaeus) support limited migration and translocation of Greek populations

Numerous studies have shown that the phylogeography of many species, including European brown hare, has been affected by the climatic oscillations of the Pleistocene. During this period the Balkans acted as a major refugium offering habitable conditions for many species. However, few studies have focused on the specific role of the Greek peninsula in the phylogeographic history of species in this southernmost margin of Balkans. We, therefore analyzed a 528 bp fragment of the D-loop region of mitochondrial DNA in 154 wild brown hare individuals from unsampled areas from both mainland and island Greece and compared it to 310 available brown hare sequences (including 110 Greek samples). Newly identified haplotypes show characteristic distribution in specific Greek areas reinforcing the theory that Greece can be considered as a subrefuge within Balkans for a number of species, with several “refugia within refugia” spots, holding significant genetic diversity. No haplotypes from wild Greek individuals clustered with the Central and Western Europe group revealing a minimal contribution of this area to the colonization of central Europe. One hundred and ten reared brown hares were also analyzed to elucidate the impact of introductions on local populations. Most of these samples presented close genetic affinity with haplotypes from Central and Western Europe indicating that farms in Greece use breeders from those areas. Therefore, despite human translocation of individuals, the genetic structure of brown hare has mostly been influenced by paleoclimatic conditions and minimally by human actions.

Approximate Bayesian computation analysis of EST-associated microsatellites indicates that the broadleaved evergreen tree Castanopsis sieboldii survived the Last Glacial Maximum in multiple refugia in Japan

Climatic changes have played major roles in plants' evolutionary history. Glacial oscillations have been particularly important, but some of their effects on plants' populations are poorly understood, including the numbers and locations of refugia in Asian warm temperate zones. In the present study, we investigated the demographic history of the broadleaved evergreen tree species Castanopsis sieboldii (Fagaceae) during the last glacial period in Japan. We used approximate Bayesian computation (ABC) for model comparison and parameter estimation for the demographic modeling using 27 EST-associated microsatellites. We also performed the species distribution modeling (SDM). The results strongly support a demographic scenario that the Ryukyu Islands and the western parts in the main islands (Kyushu and western Shikoku) were derived from separate refugia and the eastern parts in the main islands and the Japan Sea groups were diverged from the western parts prior to the coldest stage of the Last Glacial Maximum (LGM). Our data indicate that multiple refugia survived at least one in the Ryukyu Islands, and the other three regions of the western and eastern parts and around the Japan Sea of the main islands of Japan during the LGM. The SDM analysis also suggests the potential habitats under LGM climate conditions were mainly located along the Pacific Ocean side of the coastal region. Our ABC-based study helps efforts resolve the demographic history of a dominant species in warm temperate broadleaved forests during and after the last glacial period, which provides a basic model for future phylogeographical studies using this approach.

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.

Keep your feet warm? A cryptic refugium of trees linked to a geothermal spring in an ocean of glaciers

Up to now, the most widely accepted idea of the periglacial environment is that of treeless ecosystems such as the arctic or the alpine tundra, also called the tabula rasa paradigm. However, several palaeoecological studies have recently challenged this idea, that is, treeless environments in periglacial areas where all organisms would have been exterminated near the glacier formed during the Last Glacial Maximum, notably in the Scandinavian mountains. In the Alps, the issue of glacial refugia of trees remains unanswered. Advances in glacier reconstructions show that ice domes did not cover all upper massifs, but glaciers filled valleys. Here, we used fossils of plant and malacofauna from a travertine formation located in a high mountain region to demonstrate that trees (Pinus, Betula) grew with grasses during the Lateglacial-Holocene transition, while the glacier fronts were 200-300 m lower. The geothermal travertine started to accumulate more than 14,500 years ago, but became progressively more meteogene about 11,500 years ago due to a change in groundwater circulation. With trees, land snails (gastropods) associated to woody or open habitats and aquatic mollusc were also present at the onset of the current interglacial, namely the Holocene. The geothermal spring, due to warm water and soil, probably favoured woody glacial ecosystems. This new finding of early tree growth, combined with other scattered proofs of the tree presence before 11,000 years ago in the western Alps, changes our view of the tree distribution in periglacial environments, supporting the notion of tree refugia on nunataks in an ocean of glaciers. Therefore, the tabula rasa paradigm must be revisited because it has important consequences on the global changes, including postglacial plant migrations and biogeochemical cycles.

Plant DNA metabarcoding of lake sediments: How does it represent the contemporary vegetation

Metabarcoding of lake sediments have been shown to reveal current and past biodiversity, but little is known about the degree to which taxa growing in the vegetation are represented in environmental DNA (eDNA) records. We analysed composition of lake and catchment vegetation and vascular plant eDNA at 11 lakes in northern Norway. Out of 489 records of taxa growing within 2 m from the lake shore, 17–49% (mean 31%) of the identifiable taxa recorded were detected with eDNA. Of the 217 eDNA records of 47 plant taxa in the 11 lakes, 73% and 12% matched taxa recorded in vegetation surveys within 2 m and up to about 50 m away from the lakeshore, respectively, whereas 16% were not recorded in the vegetation surveys of the same lake. The latter include taxa likely overlooked in the vegetation surveys or growing outside the survey area. The percentages detected were 61, 47, 25, and 15 for dominant, common, scattered, and rare taxa, respectively. Similar numbers for aquatic plants were 88, 88, 33 and 62%, respectively. Detection rate and taxonomic resolution varied among plant families and functional groups with good detection of e.g. Ericaceae, Roseaceae, deciduous trees, ferns, club mosses and aquatics. The representation of terrestrial taxa in eDNA depends on both their distance from the sampling site and their abundance and is sufficient for recording vegetation types. For aquatic vegetation, eDNA may be comparable with, or even superior to, in-lake vegetation surveys and may therefore be used as an tool for biomonitoring. For reconstruction of terrestrial vegetation, technical improvements and more intensive sampling is needed to detect a higher proportion of rare taxa although DNA of some taxa may never reach the lake sediments due to taphonomical constrains. Nevertheless, eDNA performs similar to conventional methods of pollen and macrofossil analyses and may therefore be an important tool for reconstruction of past vegetation.

High-Throughput DNA sequencing of ancient wood

Reconstructing the colonization and demographic dynamics that gave rise to extant forests is essential to forecasts of forest responses to environmental changes. Classical approaches to map how population of trees changed through space and time largely rely on pollen distribution patterns, with only a limited number of studies exploiting DNA molecules preserved in wooden tree archaeological and subfossil remains. Here, we advance such analyses by applying high-throughput (HTS) DNA sequencing to wood archaeological and subfossil material for the first time, using a comprehensive sample of 167 European white oak waterlogged remains spanning a large temporal (from 550 to 9,800 years) and geographical range across Europe. The successful characterization of the endogenous DNA and exogenous microbial DNA of 140 (~83%) samples helped the identification of environmental conditions favouring long-term DNA preservation in wood remains, and started to unveil the first trends in the DNA decay process in wood material. Additionally, the maternally inherited chloroplast haplotypes of 21 samples from three periods of forest human-induced use (Neolithic, Bronze Age and Middle Ages) were found to be consistent with those of modern populations growing in the same geographic areas. Our work paves the way for further studies aiming at using ancient DNA preserved in wood to reconstruct the micro-evolutionary response of trees to climate change and human forest management.

Evolutionary Patterns and Processes: Lessons from Ancient DNA

Ever since its emergence in 1984, the field of ancient DNA has struggled to overcome the challenges related to the decay of DNA molecules in the fossil record. With the recent development of high-throughput DNA sequencing technologies and molecular techniques tailored to ultra-damaged templates, it has now come of age, merging together approaches in phylogenomics, population genomics, epigenomics, and metagenomics. Leveraging on complete temporal sample series, ancient DNA provides direct access to the most important dimension in evolution—time, allowing a wealth of fundamental evolutionary processes to be addressed at unprecedented resolution. This review taps into the most recent findings in ancient DNA research to present analyses of ancient genomic and metagenomic data.

DNA metabarcoding of orchid-derived products reveals widespread illegal orchid trade

In eastern Mediterranean countries orchids continue to be collected from the wild for the production of salep, a beverage made of dried orchid tubers. In this study we used nrITS1 and nrITS2 DNA metabarcoding to identify orchid and other plant species present in 55 commercial salep products purchased in Iran, Turkey, Greece and Germany. Thirty samples yielded a total of 161 plant taxa, and 13 products (43%) contained orchid species and these belonged to 10 terrestrial species with tuberous roots. Another 70% contained the substitute ingredient Cyamopsis tetraganoloba (Guar). DNA metabarcoding using the barcoding markers nrITS1 and nrITS2 shows the potential of these markers and approach for identification of species used in salep products. The analysis of interspecific genetic distances between sequences of these markers for the most common salep orchid genera shows that species level identifications can be made with a high level of confidence. Understanding the species diversity and provenance of salep orchid tubers will enable the chain of commercialization of endangered species to be traced back to the harvesters and their natural habitats, and thus allow for targeted efforts to protect or sustainably use wild populations of these orchids.

Population genomics of Mesolithic Scandinavia: Investigating early postglacial migration routes and high-latitude adaptation

Scandinavia was one of the last geographic areas in Europe to become habitable for humans after the Last Glacial Maximum (LGM). However, the routes and genetic composition of these postglacial migrants remain unclear. We sequenced the genomes, up to 57× coverage, of seven hunter-gatherers excavated across Scandinavia and dated from 9,500–6,000 years before present (BP). Surprisingly, among the Scandinavian Mesolithic individuals, the genetic data display an east–west genetic gradient that opposes the pattern seen in other parts of Mesolithic Europe. Our results suggest two different early postglacial migrations into Scandinavia: initially from the south, and later, from the northeast. The latter followed the ice-free Norwegian north Atlantic coast, along which novel and advanced pressure-blade stone-tool techniques may have spread. These two groups met and mixed in Scandinavia, creating a genetically diverse population, which shows patterns of genetic adaptation to high latitude environments. These potential adaptations include high frequencies of low pigmentation variants and a gene region associated with physical performance, which shows strong continuity into modern-day northern Europeans.

The Scandinavian peninsula was the last part of Europe to be colonized after the Last Glacial Maximum. The migration routes, cultural networks, and the genetic makeup of the first Scandinavians remain elusive and several hypotheses exist based on archaeology, climate modeling, and genetics. By analyzing the genomes of early Scandinavian hunter-gatherers, we show that their migrations followed two routes: one from the south and another from the northeast along the ice-free Norwegian Atlantic coast. These groups met and mixed in Scandinavia, creating a population more diverse than contemporaneous central and western European hunter-gatherers. As northern Europe is associated with cold and low light conditions, we investigated genomic patterns of adaptation to these conditions and genes known to be involved in skin pigmentation. We demonstrate that Mesolithic Scandinavians had higher levels of light pigmentation variants compared to the respective source populations of the migrations, suggesting adaptation to low light levels and a surprising signal of genetic continuity in TMEM131, a gene that may be involved in long-term adaptation to the cold.

Looking for relationships between the populations of Dothistroma septosporum in northern Europe and Asia

Dothistroma septosporum, a notorious pine needle pathogen with an unknown historical geographic origin and poorly known distribution pathways, is nowadays found almost in all areas inhabited by pines (Pinus spp.). The main aim of this study was to determine the relationship between North European and East Asian populations. In total, 238 Eurasian D. septosporum isolates from 11 countries, including 211 isolates from northern Europe, 16 isolates from Russian Far East and 11 isolates from Bhutan were analysed using 11 species-specific microsatellite and mating type markers. The most diverse populations were found in northern Europe, including the Baltic countries, Finland and European Russia. Notably, D. septosporum has not caused heavy damage to P. sylvestris in northern Europe, which may suggest a long co-existence of the host and the pathogen. No indication was obtained that the Russian Far East or Bhutan could be the indigenous area of D. septosporum, as the genetic diversity of the fungus there was low and evidence suggests gene flow from northern Europe to Russian Far East. On the western coast of Norway, a unique genetic pattern was observed, which differed from haplotypes dominating other Fennoscandian populations. As an agent of dothistroma needle blight, only D. septosporum was documented in northern Europe and Asia, while D. pini was found in Ukraine and Serbia.

Fire ecology of a tree glacial refugium on a nunatak with a view on Alpine glaciers

In paleoecology, the function of biomass as a fire driver has become a focus of attention in cold ecosystems, and concerns have been raised about climate in this context. Little is known about the fire frequency and fire-plant relationships during glaciation when woodlands were limited and the climate was cold. Fire history and tree biomass were reconstructed from sedimentary charcoal and macroremains, respectively, archived in lake sediments from the western Alps. Two nunataks were investigated, both with lacustrine sediments covering the last 21 000 yr at least. During the Last Glacial Maximum (LGM) and the Lateglacial, fires occurred only on the nunatak sheltering woody plants. Cembra pine (Pinus cembra) and larch (Larix decidua) survived above glaciers during the LGM, thus evidencing a biological refugium and supporting the nunatak theory. We highlighted a long-term relationship between fires and dominant trees over the last 21 000 yr, where fire frequencies track the global climate and the local changes in tree biomass. Glacial climate (dry, cold) does not rule out fires. Fuel load and composition were significant fire drivers, with cembra pine dominating during colder periods with rare fires, and larch during the warmer Holocene with frequent fires. These findings increase knowledge of fire ecology in cold environments, and open perspectives in tree population genetics by considering new areas of tree glacial refugia in Europe.

Genome-Wide Analysis of Colonization History and Concomitant Selection in Arabidopsis lyrata

The high climatic variability in the past hundred thousand years has affected the demographic and adaptive processes in many species, especially in boreal and temperate regions undergoing glacial cycles. This has also influenced the patterns of genome-wide nucleotide variation, but the details of these effects are largely unknown. Here we study the patterns of genome-wide variation to infer colonization history and patterns of selection of the perennial herb species Arabidopsis lyrata, in locally adapted populations from different parts of its distribution range (Germany, UK, Norway, Sweden, and USA) representing different environmental conditions. Using site frequency spectra based demographic modeling, we found strong reduction in the effective population size of the species in general within the past 100,000 years, with more pronounced effects in the colonizing populations. We further found that the northwestern European A. lyrata populations (UK and Scandinavian) are more closely related to each other than with the Central European populations, and coalescent based population split modeling suggests that western European and Scandinavian populations became isolated relatively recently after the glacial retreat. We also highlighted loci showing evidence for local selection associated with the Scandinavian colonization. The results presented here give new insights into postglacial Scandinavian colonization history and its genome-wide effects.

Effects of Mountain Uplift and Climatic Oscillations on Phylogeography and Species Divergence in Four Endangered Notopterygium Herbs

Mountain uplift and climatic fluctuations are important driving forces that have affected the geographic distribution and population dynamics history of organisms. However, it is unclear how geological and climatic events might have affected the phylogeographic history and species divergence in high-alpine herbal plants. In this study, we analyzed the population demographic history and species differentiation of four endangered Notopterygium herbs on the high-altitude Qinghai-Tibetan Plateau (QTP) and adjacent areas. We combined phylogeographic analysis with species distribution modeling to detect the genetic variations in four Notopterygium species (N. incisum, N. franchetii, N. oviforme, and N. forrestii). In total, 559 individuals from 74 populations of the four species were analyzed based on three maternally inherited chloroplast fragments (matK, rbcL, and trnS-trnG) and one nuclear DNA region (internal transcribed spacer, ITS). Fifty-five chloroplast DNA (cpDNA) and 48 ITS haplotypes were identified in the four species. All of the cpDNA and ITS haplotypes were species-specific, except N. franchetii and N. oviforme shared one cpDNA haplotype, H32. Phylogenetic analysis suggested that all four species formed a monophyletic clade with high bootstrap support, where N. franchetii and N. oviforme were sisters. In addition, each Notopterygium species generated an individual clade that corresponded to their respective species in the ITS tree. Population dynamics analyses and species distribution modeling showed that the two widely distributed herbs N. incisum and N. franchetii exhibited obvious demographic expansions during the Pleistocene ice ages. Molecular dating suggested that the divergence of the four Notopterygium species occurred approximately between 3.6 and 1.2 Mya, and it was significantly associated with recent extensive uplifts of the QTP. Our results support the hypothesis that mountain uplift and Quaternary climatic oscillations profoundly shaped the population genetic divergence and demographic dynamics of Notopterygium species. The findings of this and previous studies provide important insights into the effects of QTP uplifts and climatic changes on phylogeography and species differentiation in high altitude mountainous areas. Our results may also facilitate the conservation of endangered herbaceous medicinal plants in the genus Notopterygium.

Demographic inference from whole-genome and RAD sequencing data suggests alternating human impacts on goose populations since the last ice age

We investigated how population changes and fluctuations in the pink-footed goose might have been affected by climatic and anthropogenic factors. First, genomic data confirmed the existence of two separate populations: western (Iceland) and eastern (Svalbard/Denmark). Second, demographic inference suggests that the species survived the last glacial period as a single ancestral population with a low population size (100-1,000 individuals) that split into the current populations at the end of the last glacial maximum with Iceland being the most plausible glacial refuge. While population changes during the last glaciation were clearly environmental, we hypothesize that more recent demographic changes are human-related: (1) the inferred population increase in the Neolithic is due to deforestation to establish new lands for agriculture, increasing available habitat for pink-footed geese, (2) the decline inferred during the Middle Ages is due to human persecution, and (3) improved protection explains the increasing demographic trends during the 20th century. Our results suggest both environmental (during glacial cycles) and anthropogenic effects (more recent) can be a threat to species survival.

The paleosymbiosis hypothesis: host plants can be colonised by root symbionts that have been inactive for centuries to millenia

Paleoecologists have speculated that post-glacial migration of tree species could have been facilitated by mycorrhizal symbionts surviving glaciation as resistant propagules belowground. The general premise of this idea, which we call the 'paleosymbiosis hypothesis', is that host plants can access and be colonised by fungal root symbionts that have been inactive for millennia. Here, we explore the plausibility of this hypothesis by synthesising relevant findings from a diverse literature. For example, the paleoecology literature provided evidence of modern roots penetrating paleosols containing ancient (>6000 years) fungal propagules, though these were of unknown condition. With respect to propagule longevity, the available evidence is of mixed quality, but includes convincing examples consistent with the paleosymbiosis hypothesis (i.e. >1000 years viable propagules). We describe symbiont traits and environmental conditions that should favour the development and preservation of an ancient propagule bank, and discuss the implications for our understanding of soil symbiont diversity and ecosystem functioning. We conclude that the paleosymbiosis hypothesis is plausible in locations where propagule deposition and preservation conditions are favourable (e.g. permafrost regions). We encourage future belowground research to consider and explore the potential temporal origins of root symbioses.

Inferring the demographic history of an oligophagous grasshopper: Effects of climatic niche stability and host-plant distribution

Understanding the consequences of past environmental changes on the abiotic and biotic components of the landscape and deciphering their impacts on the demographic trajectories of species is a major issue in evolutionary biogeography. In this study, we combine nuclear and mitochondrial genetic data to study the phylogeographical structure and lineage-specific demographic histories of the scrub-legume grasshopper (Chorthippus binotatus binotatus), a montane taxon distributed in the Iberian Peninsula and France that exclusively feeds on certain scrub-legume species. Genetic data and paleo-distribution modelling indicate the presence of four main lineages that seem to have diverged in allopatry and long-term persisted in Iberian and French refugia since the Mid Pleistocene. Comparisons of different demographic hypotheses in an Approximate Bayesian Computation (ABC) framework supported a population bottleneck in the northwestern French clade and paleo-distribution modelling indicate that the populations of this lineage have experienced more severe environmental fluctuations during the last 21 000 years than those from the Iberian Peninsula. Accordingly, we found that nuclear genetic diversity of the populations of scrub-legume grasshopper is positively associated with local stability of suitable habitats defined by both Pleistocene climate changes and historical distributional shifts of host-plant species. Overall, our study highlights the importance of integrating the potential effects of abiotic (i.e. climate and geography) and biotic components (i.e. inter-specific interactions) into the study of the evolutionary and demographic history of specialist taxa with narrow ecological requirements.