A genome-wide genetic association study reveals SNPs significantly associated with environmental variables and specific leaf area in European beech.

European beech is negatively affected by climate change and a further growth decline is predicted for large parts of its distribution range. Despite the importance of this species, little is known about its genetic adaptation and especially the genetic basis of its physiological traits. Here, we used genotyping by sequencing to identify SNPs in 43 German European beech populations growing under different environmental conditions. In total, 28 of these populations were located along a precipitation and temperature gradient in northern Germany, and single tree-based hydraulic and morphological traits were available. We obtained a set of 13,493 high-quality SNPs that were used for environmental and SNP-trait association analysis. In total, 22 SNPs were identified that were significantly associated with environmental variables or specific leaf area (SLA). Several SNPs were located in genes related to stress response. The majority of the significant SNPs were located in non-coding (intergenic and intronic) regions. These may be in linkage disequilibrium with the causative coding or regulatory regions. Our study gives insights into the genetic basis of abiotic adaptation in European beech, and provides genetic resources that can be used in future studies on this species. Besides clear patterns of local adaptation to environmental conditions of the investigated populations, the analyzed morphological and hydraulic traits explained most of the explainable genetic variation. Thus, they could successfully be altered in tree breeding programs, which may help to increase the adaptation of European beech to changing environmental conditions in the future.

Leaf gene expression trajectories during the growing season are consistent between sites and years in American beech.

Transcriptomics provides a versatile tool for ecological monitoring. Here, through genome-guided profiling of transcripts mapping to 33 042 gene models, expression differences can be discerned among multi-year and seasonal leaf samples collected from American beech trees at two latitudinally separated sites. Despite a bottleneck due to post-Columbian deforestation, the single nucleotide polymorphism-based population genetic background analysis has yielded sufficient variation to account for differences between populations and among individuals. Our expression analyses during spring-summer and summer-autumn transitions for two consecutive years involved 4197 differentially expressed protein coding genes. Using Populus orthologues we reconstructed a protein-protein interactome representing leaf physiological states of trees during the seasonal transitions. Gene set enrichment analysis revealed gene ontology terms that highlight molecular functions and biological processes possibly influenced by abiotic forcings such as recovery from drought and response to excess precipitation. Further, based on 324 co-regulated transcripts, we focused on a subset of GO terms that could be putatively attributed to late spring phenological shifts. Our conservative results indicate that extended transcriptome-based monitoring of forests can capture diverse ranges of responses including air quality, chronic disease, as well as herbivore outbreaks that require activation and/or downregulation of genes collectively tuning reaction norms maintaining the survival of long living trees such as the American beech.

Development of polymorphic microsatellite markers for Fagus pashanica (Fagaceae) using next-generation sequencing.

Fagus pashanica is an endangered and endemic tree species in China. To understand its genetic diversity and structure for effective conservation, we used next-generation sequencing data to develop a set of microsatellite markers. Twenty-three of the 68 designed loci were successfully amplified. Fifteen polymorphic loci with clear peaks were selected for further analyses in three F. pashanica populations sampled from Nanjiang, Wangcang and Pingwu counties in Sichuan Province, China. The number of alleles per locus ranged from two to 11. The levels of observed and expected heterozygosity ranged from 0.033-0.852 and 0.033-0.787, respectively. All 23 loci were also successfully amplified in F. longipetiolata and F. lucida, and 19 were successfully amplified in F. engleriana. These microsatellite markers will be useful for population genetic studies of F. pashanica and other Fagus species.

First genome-wide data from Italian European beech (Fagus sylvatica L.): Strong and ancient differentiation between Alps and Apennines.

The European beech (Fagus sylvatica L.) is one of the most widespread forest trees in Europe whose distribution and intraspecific diversity has been largely shaped by repeated glacial cycles. Previous studies, mainly based on palaeobotanical evidence and a limited set of chloroplast and nuclear genetic markers, highlighted a complex phylogeographic scenario, with southern and western Europe characterized by a rather heterogeneous genetic structure, as a result of recolonization from different glacial refugia. Despite its ecological and economic importance, the genome of this broad-leaved tree has only recently been assembled, and its intra-species genomic diversity is still largely unexplored. Here, we performed whole-genome resequencing of nine Italian beech individuals sampled from two stands located in the Alpine and Apennine mountain ranges. We investigated patterns of genetic diversity at chloroplast, mitochondrial and nuclear genomes and we used chloroplast genomes to reconstruct a temporally-resolved phylogeny. Results allowed us to test European beech differentiation on a whole-genome level and to accurately date their divergence time. Our results showed comparable, relatively high levels of genomic diversity in the two populations and highlighted a clear differentiation at chloroplast, mitochondrial and nuclear genomes. The molecular clock analysis indicated an ancient split between the Alpine and Apennine populations, occurred between the Günz and the Riss glaciations (approximately 660 kyrs ago), suggesting a long history of separation for the two gene pools. This information has important conservation implications in the context of adaptation to ongoing climate changes.

Development of Primers Specific for Detection of Litylenchus crenatae, the Causal Agent of Beech Leaf Disease, in Plant Tissue.

Beech leaf disease (BLD), an emerging threat to American beech (Fagus grandifolia) in the northern United States and Canada, was recently confirmed to be caused by the nematode Litylenchus crenatae subsp. mccannii (hereafter L. crenatae). Consequently, there is a need for a rapid, sensitive, and accurate method for detecting L. crenatae for both diagnostic as well as control purposes. This research developed a new set of DNA primers that specifically amplify L. crenatae and allow for accurate detection of the nematode in plant tissue. These primers have also been used in quantitative PCR (qPCR) to determine relative differences in gene copy number between samples. This primer set provides an improved, effective tool for monitoring and detecting L. crenatae in temperate tree leaf tissue which is necessary to understand the spread of this emerging forest pest and to develop management strategies.

Both abundant and rare fungi colonizing Fagus sylvatica ectomycorrhizal root-tips shape associated bacterial communities.

Ectomycorrhizal fungi live in close association with their host plants and form complex interactions with bacterial/archaeal communities in soil. We investigated whether abundant or rare ectomycorrhizal fungi on root-tips of young beech trees (Fagus sylvatica) shape bacterial/archaeal communities. We sequenced 16S rRNA genes and fungal internal transcribed spacer regions of individual root-tips and used ecological networks to detect the tendency of certain assemblies of fungal and bacterial/archaeal taxa to inhabit the same root-tip (i.e. modularity). Individual ectomycorrhizal root-tips hosted distinct fungal communities associated with unique bacterial/archaeal communities. The structure of the fungal-bacterial/archaeal association was determined by both, dominant and rare fungi. Integrating our data in a conceptual framework suggests that the effect of rare fungi on the bacterial/archaeal communities of ectomycorrhizal root-tips contributes to assemblages of bacteria/archaea on root-tips. This highlights the potential impact of complex fine-scale interactions between root-tip associated fungi and other soil microorganisms for the ectomycorrhizal symbiosis.

5S-IGS rDNA in wind-pollinated trees (Fagus L.) encapsulates 55 million years of reticulate evolution and hybrid origins of modern species.

Standard models of plant speciation assume strictly dichotomous genealogies in which a species, the ancestor, is replaced by two offspring species. The reality in wind-pollinated trees with long evolutionary histories is more complex: species evolve from other species through isolation when genetic drift exceeds gene flow; lineage mixing can give rise to new species (hybrid taxa such as nothospecies and allopolyploids). The multi-copy, potentially multi-locus 5S rDNA is one of few gene regions conserving signal from dichotomous and reticulate evolutionary processes down to the level of intra-genomic recombination. Therefore, it can provide unique insights into the dynamic speciation processes of lineages that diversified tens of millions of years ago. Here, we provide the first high-throughput sequencing (HTS) of the 5S intergenic spacers (5S-IGS) for a lineage of wind-pollinated subtropical to temperate trees, the Fagus crenata - F. sylvatica s.l. lineage, and its distant relative F. japonica. The observed 4963 unique 5S-IGS variants reflect a complex history of hybrid origins, lineage sorting, mixing via secondary gene flow, and intra-genomic competition between two or more paralogous-homoeologous 5S rDNA lineages. We show that modern species are genetic mosaics and represent a striking case of ongoing reticulate evolution during the past 55 million years.

Complete Chloroplast Genomes of Fagus sylvatica L. Reveal Sequence Conservation in the Inverted Repeat and the Presence of Allelic Variation in NUPTs.

Growing amounts of genomic data and more efficient assembly tools advance organelle genomics at an unprecedented scale. Genomic resources are increasingly used for phylogenetic analyses of many plant species, but are less frequently used to investigate within-species variability and phylogeography. In this study, we investigated genetic diversity of Fagus sylvatica, an important broadleaved tree species of European forests, based on complete chloroplast genomes of 18 individuals sampled widely across the species distribution. Our results confirm the hypothesis of a low cpDNA diversity in European beech. The chloroplast genome size was remarkably stable (158,428 ± 37 bp). The polymorphic markers, 12 microsatellites (SSR), four SNPs and one indel, were found only in the single copy regions, while inverted repeat regions were monomorphic both in terms of length and sequence, suggesting highly efficient suppression of mutation. The within-individual analysis of polymorphisms showed >9k of markers which were proportionally present in gene and non-gene areas. However, an investigation of the frequency of alternate alleles revealed that the source of this diversity originated likely from nuclear-encoded plastome remnants (NUPTs). Phylogeographic and Mantel correlation analysis based on the complete chloroplast genomes exhibited clustering of individuals according to geographic distance in the first distance class, suggesting that the novel markers and in particular the cpSSRs could provide a more detailed picture of beech population structure in Central Europe.

Genetic signatures of divergent selection in European beech (Fagus sylvatica L.) are associated with the variation in temperature and precipitation across its distribution range.

High genetic variation and extensive gene flow may help forest trees with adapting to ongoing climate change, yet the genetic bases underlying their adaptive potential remain largely unknown. We investigated range-wide patterns of potentially adaptive genetic variation in 64 populations of European beech (Fagus sylvatica L.) using 270 SNPs from 139 candidate genes involved either in phenology or in stress responses. We inferred neutral genetic structure and processes (drift and gene flow) and performed differentiation outlier analyses and gene-environment association (GEA) analyses to detect signatures of divergent selection. Beech range-wide genetic structure was consistent with the species' previously identified postglacial expansion scenario and recolonization routes. Populations showed high diversity and low differentiation along the major expansion routes. A total of 52 loci were found to be putatively under selection and 15 of them turned up in multiple GEA analyses. Temperature and precipitation related variables were equally represented in significant genotype-climate associations. Signatures of divergent selection were detected in the same proportion for stress response and phenology-related genes. The range-wide adaptive genetic structure of beech appears highly integrated, suggesting a balanced contribution of phenology and stress-related genes to local adaptation, and of temperature and precipitation regimes to genetic clines. Our results imply a best-case scenario for the maintenance of high genetic diversity during range shifts in beech (and putatively other forest trees) with a combination of gene flow maintaining within-population neutral diversity and selection maintaining between-population adaptive differentiation.

Genomic signatures of natural selection at phenology-related genes in a widely distributed tree species Fagus sylvatica L.

BACKGROUND: Diversity among phenology-related genes is predicted to be a contributing factor in local adaptations seen in widely distributed plant species that grow in climatically variable geographic areas, such as forest trees. European beech (Fagus sylvatica L.) is widespread, and is one of the most important broadleaved tree species in Europe; however, its potential for adaptation to climate change is a matter of uncertainty, and little is known about the molecular basis of climate change-relevant traits like bud burst. RESULTS: We explored single nucleotide polymorphisms (SNP) at candidate genes related to bud burst in beech individuals sampled across 47 populations from Europe. SNP diversity was monitored for 380 candidate genes using a sequence capture approach, providing 2909 unlinked SNP loci. We used two complementary analytical methods to find loci significantly associated with geographic variables, climatic variables (expressed as principal components), or phenotypic variables (spring and autumn phenology, height, survival). Redundancy analysis (RDA) was used to detect candidate markers across two spatial scales (entire study area and within subregions). We revealed 201 candidate SNPs at the broadest scale, 53.2% of which were associated with phenotypic variables. Additive polygenic scores, which provide a measure of the cumulative signal across significant candidate SNPs, were correlated with a climate variable (first principal component, PC1) related to temperature and precipitation availability, and spring phenology. However, different genotype-environment associations were identified within Southeastern Europe as compared to the entire geographic range of European beech. CONCLUSIONS: Environmental conditions play important roles as drivers of genetic diversity of phenology-related genes that could influence local adaptation in European beech. Selection in beech favors genotypes with earlier bud burst under warmer and wetter habitats within its range; however, selection pressures may differ across spatial scales.

Genomic basis for drought resistance in European beech forests threatened by climate change.

In the course of global climate change, Central Europe is experiencing more frequent and prolonged periods of drought. The drought years 2018 and 2019 affected European beeches (Fagus sylvatica L.) differently: even in the same stand, drought-damaged trees neighboured healthy trees, suggesting that the genotype rather than the environment was responsible for this conspicuous pattern. We used this natural experiment to study the genomic basis of drought resistance with Pool-GWAS. Contrasting the extreme phenotypes identified 106 significantly associated single-nucleotide polymorphisms (SNPs) throughout the genome. Most annotated genes with associated SNPs (>70%) were previously implicated in the drought reaction of plants. Non-synonymous substitutions led either to a functional amino acid exchange or premature termination. An SNP assay with 70 loci allowed predicting drought phenotype in 98.6% of a validation sample of 92 trees. Drought resistance in European beech is a moderately polygenic trait that should respond well to natural selection, selective management, and breeding.

Climate change is having a serious impact on many ecosystems. In the summer of 2018 and 2019, around two thirds of European beech trees were damaged or killed by extreme drought. It is critical to keep these beech woods healthy, as they are central to the survival of over 6,000 other species of animals and plants. The level of damage caused by the drought varied between forests. However, not all the trees in each forest responded in the same way, with severely damaged trees often sitting next to fully healthy ones. This suggests that the genetic make-up of each tree determines how well it can adapt to drought rather than its local environment. To investigate this further, Pfenninger et al. studied the genome of over 400 European beech trees from the Hesse region in Germany. The samples came from pairs of neighbouring trees that had responded differently to the droughts. The analysis found more than 80 parts of the genome that differed between healthy and damaged trees. Pfenninger et al. then used this information to create a genetic test which can quickly and inexpensively predict how well an individual beech tree might survive in a drought. Applying this test to another 92 trees revealed that it can reliably detect which ones were healthy and which ones were damaged. Beech forests are typically managed by private owners, agencies or breeders that could use this genetic test to select and reproduce trees that are better adapted to drought. The goal now is to develop the test so that it can be used more widely to manage European beech trees and potentially other species.

Integrated biological and chemical characterisation of a pair of leonardesque canal lock gates.

The Museo Nazionale della Scienza e della Tecnologia "Leonardo da Vinci" in Milan is exposing two pairs of canal lock gates, used to control the water flow in Milan canal system, whose design appears in the Leonardo's Codex Atlanticus. The wood present in the gates has been deeply characterised by mean of a multidisciplinary investigation involving i) DNA barcoding of wood fragments; ii) microbial community characterisation, and iii) chemical analyses. DNA barcoding revealed that two fragments of the gates belonged to wood species widely used in the middle age: Fagus sylvatica and Picea abies. The chemical characterisations were based on the use of ionic liquid as dissolving medium in order to analyse the entire cell wall material by means of Gel Permeation Chromatography (GPC) and 2D-NMR-HSQC techniques. This multidisciplinary analytical approach was able to highlight the complex nature of the degradation occurred during the gate operation (XVI-XVIII centuries): an intricate interplay between microbial populations (i.e. Shewanella), inorganic factors (i.e. iron from nails), physical factors and the lignocellulosic material.

Comparative genome analyses suggest a hemibiotrophic lifestyle and virulence differences for the beech bark disease fungal pathogens Neonectria faginata and Neonectria coccinea.

Neonectria faginata and Neonectria coccinea are the causal agents of the insect-fungus disease complex known as beech bark disease (BBD), known to cause mortality in beech forest stands in North America and Europe. These fungal species have been the focus of extensive ecological and disease management studies, yet less progress has been made toward generating genomic resources for both micro- and macro-evolutionary studies. Here, we report a 42.1 and 42.7 mb highly contiguous genome assemblies of N. faginata and N. coccinea, respectively, obtained using Illumina technology. These species share similar gene number counts (12,941 and 12,991) and percentages of predicted genes with assigned functional categories (64 and 65%). Approximately 32% of the predicted proteomes of both species are homologous to proteins involved in pathogenicity, yet N. coccinea shows a higher number of predicted mitogen-activated protein kinase genes, virulence determinants possibly contributing to differences in disease severity between N. faginata and N. coccinea. A wide range of genes encoding for carbohydrate-active enzymes capable of degradation of complex plant polysaccharides and a small number of predicted secretory effector proteins, secondary metabolite biosynthesis clusters and cytochrome oxidase P450 genes were also found. This arsenal of enzymes and effectors correlates with, and reflects, the hemibiotrophic lifestyle of these two fungal pathogens. Phylogenomic analysis and timetree estimations indicated that the N. faginata and N. coccinea species divergence may have occurred at ∼4.1 million years ago. Differences were also observed in the annotated mitochondrial genomes as they were found to be 81.7 kb (N. faginata) and 43.2 kb (N. coccinea) in size. The mitochondrial DNA expansion observed in N. faginata is attributed to the invasion of introns into diverse intra- and intergenic locations. These first draft genomes of N. faginata and N. coccinea serve as valuable tools to increase our understanding of basic genetics, evolutionary mechanisms and molecular physiology of these two nectriaceous plant pathogenic species.

Double-digest RAD-sequencing: do pre- and post-sequencing protocol parameters impact biological results?

Next-generation sequencing technologies have opened a new era of research in population genetics. Following these new sequencing opportunities, the use of restriction enzyme-based genotyping techniques, such as restriction site-associated DNA sequencing (RAD-seq) or double-digest RAD-sequencing (ddRAD-seq), has dramatically increased in the last decade. From DNA sampling to SNP calling, the laboratory and bioinformatic parameters of enzyme-based techniques have been investigated in the literature. However, the impact of those parameters on downstream analyses and biological results remains less documented. In this study, we investigated the effects of sevral pre- and post-sequencing settings on ddRAD-seq results for two biological systems: a complex of butterfly species (Coenonympha sp.) and several populations of common beech (Fagus sylvatica). Our results suggest that pre-sequencing parameters (i.e., DNA quantity, number of PCR cycles during library preparation) have a significant impact on the number of recovered reads and SNPs, on the number of unique alleles and on individual heterozygosity. In the same way, we found that post-sequencing settings (i.e., clustering and minimum coverage thresholds) influenced loci reconstruction (e.g., number of loci, mean coverage) and SNP calling (e.g., number of SNPs; heterozygosity) but had only a marginal impact on downstream analyses (e.g., measure of genetic differentiation, estimation of individual admixture, and demographic inferences). In addition, replication analyses confirmed the reproducibility of the ddRAD-seq procedure. Overall, this study assesses the degree of sensitivity of ddRAD-seq data to pre- and post-sequencing protocols, and illustrates its robustness when studying population genetics.

Weak founder effects but significant spatial genetic imprint of recent contraction and expansion of European beech populations.

Understanding the ecological and evolutionary processes occurring during species range shifts is important in the current context of global change. Here, we investigate the interplay between recent expansion, gene flow and genetic drift, and their consequences for genetic diversity and structure at landscape and local scales in European beech (Fagus sylvatica L.) On Mont Ventoux, South-Eastern France, we located beech forest refugia at the time of the most recent population minimum, ~150 years ago, and sampled 71 populations (2042 trees) in both refugia and expanding populations over an area of 15,000 ha. We inferred patterns of gene flow and genetic structure using 12 microsatellite markers. We identified six plots as originating from planting, rather than natural establishment, mostly from local genetic material. Comparing genetic diversity and structure in refugia versus recent populations did not support the existence of founder effects: heterozygosity (He = 0.667) and allelic richness (Ar = 4.298) were similar, and FST was low (0.031 overall). Still, significant spatial evidence of colonization was detected, with He increasing along the expansion front, while genetic differentiation from the entire pool (ßWT) decreased. Isolation by distance was found in refugia but not in recently expanding populations. Our study indicates that beech capacities for colonization and gene flow were sufficient to preserve genetic diversity despite recent forest contraction and expansion. Because beech has long distance pollen and seed dispersal, these results illustrate a 'best case scenario' for the maintenance of high genetic diversity and adaptive potential under climate-change-related range change.

Plastome Characterization and Phylogenomics of East Asian Beeches with a Special Emphasis on Fagus multinervis on Ulleung Island, Korea.

Beech trees of the genus Fagus (Fagaceae) are monoecious and distributed in the Northern Hemisphere. They represent an important component of mixed broad-leaved evergreen-deciduous forests and are an economically important source of timber. Despite their ecological and economical importance, however, little is known regarding the overall plastome evolution among Fagus species in East Asia. In particular, the taxonomic position and status of F. multinervis, a beech species endemic to Ulleung Island of Korea, remains unclear even today. Therefore, in this study, we characterized four newly completed plastomes of East Asian Fagus species (one accession each of F. crenata and F. multinervis and two accessions of F. japonica). Moreover, we performed phylogenomic analyses comparing these four plastomes with F. sylvatica (European beech) plastome. The four plastomes were highly conserved, and their size ranged from 158,163 to 158,348 base pair (bp). The overall GC content was 37.1%, and the sequence similarity ranged from 99.8% to 99.99%. Codon usage patterns were similar among species, and 7 of 77 common protein-coding genes were under positive selection. Furthermore, we identified five highly variable hotspot regions of the Fagus plastomes (ccsA/ndhD, ndhD/psaC, ndhF/rpl32, trnS-GCU/trnG-UCC, and ycf1). Phylogenetic analysis revealed the monophyly of Fagus as well as early divergence of the subgenus Fagus and monophyletic Engleriana. Finally, phylogenetic results supported the taxonomic distinction of F. multinervis from its close relatives F. engleriana and F. japonica. However, the sister species and geographic origin of F. multinervis on Ulleung Island could not be determined.

Plant functional traits are correlated with species persistence in the herb layer of old-growth beech forests.

This paper explores which traits are correlated with fine-scale (0.25 m2) species persistence patterns in the herb layer of old-growth forests. Four old-growth beech forests representing different climatic contexts (presence or absence of summer drought period) were selected along a north-south gradient in Italy. Eight surveys were conducted in each of the sites during the period spanning 1999-2011. We found that fine-scale species persistence was correlated with different sets of plant functional traits, depending on local ecological context. Seed mass was found to be as important for the fine-scale species persistence in the northern sites, while clonal and bud-bank traits were markedly correlated with the southern sites characterised by summer drought. Leaf traits appeared to correlate with species persistence in the drier and wetter sites. However, we found that different attributes, i.e. helomorphic vs scleromorphic leaves, were correlated to species persistence in the northernmost and southernmost sites, respectively. These differences appear to be dependent on local trait adaptation rather than plant phylogenetic history. Our findings suggest that the persistent species in the old-growth forests might adopt an acquisitive resource-use strategy (i.e. helomorphic leaves with high SLA) with higher seed mass in sites without summer drought, while under water-stressed conditions persistent species have a conservative resource-use strategy (i.e. scleromorphic leaves with low SLA) with an increased importance of clonal and resprouting ability.

Assessments of fine-scale spatial patterns of SNPs in an old-growth beech forest.

The spatial patterns of non-neutral genetic variations at fine spatial scales and their possible associations with microenvironments have not been well-documented for tree populations. Based on 25-32 SNP markers, we examine whether non-neutral SNPs and their associations with microenvironments can be detected in FcMYB1603, a gene homologous to that encoding a protein induced by drought stress in Arabidopsis thaliana for the 166 adult trees in a 1-ha plot in a mature population of Fagus crenata. In the 83 individuals of a younger cohort of below canopy trees, the nonsynonymous SNP at locus FcMYB1603_684 exhibited a spatial signature representing a departure from the expected spatial patterns of neutral genetic variation. Evaluations of non-neutrality for this locus were robust against the potential risks of false positives due to the low number of SNP loci, a low criterion set for minor allele frequency, and any edge effect on the trees' spatial structure. An older cohort exhibited no signal of the existence of non-neutral genetic variation, suggesting that temporal fluctuation in the microenvironmental conditions on the forest floor may have exposed different cohorts to different magnitudes of selection pressure. Although genotypes of the locus showed a spatial association with a microenvironmental variable potentially related to soil moisture, the present study was subject to a limitation due to the generally low polymorphism of nonsynonymous loci within the single plot, which suggests that it will be important to replicate the study design in order to carry out research on fine-scale non-neutral genetic variations.

Climate-associated genetic variation in Fagus sylvatica and potential responses to climate change in the French Alps.

Local adaptation patterns have been found in many plants and animals, highlighting the genetic heterogeneity of species along their range of distribution. In the next decades, global warming is predicted to induce a change in the selective pressures that drive this adaptive variation, forcing a reshuffling of the underlying adaptive allele distributions. For species with low dispersion capacity and long generation time such as trees, the rapidity of the change could impede the migration of beneficial alleles and lower their capacity to track the changing environment. Identifying the main selective pressures driving the adaptive genetic variation is thus necessary when investigating species capacity to respond to global warming. In this study, we investigate the adaptive landscape of Fagus sylvatica along a gradient of populations in the French Alps. Using a double-digest restriction-site-associated DNA (ddRAD) sequencing approach, we identified 7,000 SNPs from 570 individuals across 36 different sites. A redundancy analysis (RDA)-derived method allowed us to identify several SNPs that were strongly associated with climatic gradients; moreover, we defined the primary selective gradients along the natural populations of F. sylvatica in the Alps. Strong effects of elevation and humidity, which contrast north-western and south-eastern site, were found and were believed to be important drivers of genetic adaptation. Finally, simulations of future genetic landscapes that used these findings allowed identifying populations at risk for F. sylvatica in the Alps, which could be helpful for future management plans.

Climate warming disrupts mast seeding and its fitness benefits in European beech.

Many plants benefit from synchronous year-to-year variation in seed production, called masting. Masting benefits plants because it increases the efficiency of pollination and satiates predators, which reduces seed loss. Here, using a 39-year-long dataset, we show that climate warming over recent decades has increased seed production of European beech but decreased the year-to-year variability of seed production and the reproductive synchrony among individuals. Consequently, the benefit that the plants gained from masting has declined. While climate warming was associated with increased reproductive effort, we demonstrate that less effective pollination and greater losses of seeds to predators offset any benefits to the plants. This shows that an apparently simple benefit of climate warming unravels because of complex ecological interactions. Our results indicate that in masting systems, the main beneficiaries of climate-driven increases in seed production are seed predators, not plants.