The dataset presented here was collected by the GenTree project (EU-Horizon 2020), which aims to improve the use of forest genetic resources across Europe by better understanding how trees adapt to their local environment. This dataset of individual tree-core characteristics including ring-width series and whole-core wood density was collected for seven ecologically and economically important European tree species: silver birch (Betula pendula), European beech (Fagus sylvatica), Norway spruce (Picea abies), European black poplar (Populus nigra), maritime pine (Pinus pinaster), Scots pine (Pinus sylvestris), and sessile oak (Quercus petraea). Tree-ring width measurements were obtained from 3600 trees in 142 populations and whole-core wood density was measured for 3098 trees in 125 populations. This dataset covers most of the geographical and climatic range occupied by the selected species. The potential use of it will be highly valuable for assessing ecological and evolutionary responses to environmental conditions as well as for model development and parameterization, to predict adaptability under climate change scenarios.

Genetic association studies in forest trees would greatly benefit from information on the response of trees to environmental stressors over time, which can be provided by dendroecological analysis. Here, we jointly analysed dendroecological and genetic data of surviving silver fir trees to explore the genetic basis of their response to the iconic stress episode of the 1970s and 1980s that led to large-scale forest dieback in Central Europe and has been attributed to air pollution. Specifically, we derived dendrophenotypic measures from 190 trees in the Bavarian Forest that characterize the resistance, resilience and recovery during this growth depression, and in the drought year in 1976. By focusing on relative growth changes of trees and by standardizing the dendrophenotypes within stands, we accounted for variation introduced by micro- and macroscale environmental differences. We associated the dendrophenotypes with single nucleotide polymorphisms (SNPs) in candidate genes using general linear models (GLMs) and the machine learning algorithm random forest with subsequent feature selection. Most trees at our study sites experienced a severe growth decline from 1974 until the mid-1980s with minimum values during the drought year. Fifteen genes were associated with the dendrophenotypes, including genes linked to photosynthesis and drought stress. With our study, we show that dendrophenotypes can be a powerful resource for genetic association studies that permit to account for micro- and macroenvironmental variation when data are derived from natural populations. We call for a wider collaboration of dendroecologists and forest geneticists to integrate individual tree-level dendrophenotypes in genetic association studies.

No abstract

Mediterranean refugial areas are generally underrepresented in large-scale genetic surveys of forest trees. In the case of silver fir (Abies alba Mill.), this has led to divergent hypotheses about the exact location of glacial refugia and the trajectory of recolonization routes. Based on the comprehensive sampling of Apennine populations, we aimed to reconcile discrepancies about the number and location of refugia for silver fir in the Apennines and test alternative demographic scenarios developed from palaeobotanical and genetic data. Location Mediterranean Basin; the Apennines and surrounding areas.Methods 1167 individuals from 16 Apennine populations, extensively covering the species’ distribution along the Italian Peninsula, and eight populations from the Alps and Eastern Europe were genotyped at 16 nuclear and three chloroplast microsatellite markers. The geographical distribution of genetic variation was explored using Bayesian clustering and multivariate methods. Based on the inferred genetic structure, the demographic history of A. alba was assessed by the approximate Bayesian computation (ABC) analysis.Results Two unexpected characteristics of genetic structure emerged: a sharp genetic boundary in the central Apennines and a tight genetic connection between southern Apennine and Eastern European gene pools. Two Apennine areas, corresponding precisely with refugial areas hypothesized in most recent palaeobotanical syntheses, have high genetic diversity on a par with Eastern European populations. ABC analysis showed an ancient separation between Apennine and Eastern European gene pools followed by an admixture event that, mainly through directional gene flow via pollen, might have established the genetic similarity between southern Apennine and Eastern European populations. In addition, there was evidence that the central Apennines acted as a small-scale, isolated refugium during the Last Glacial Maximum. Main conclusionsSilver fir rear edge populations have experienced a complex demographic history across several glacial-interglacial cycles, leading to unexpected genetic structure. Our study provides new insights into forest tree dynamics in the Mediterranean, showing the putative presence of multiple refugia for silver fir in the Apennines and a trans-Adriatic connection between silver fir populations in the southern Italy and the Balkans

BackgroundLocal adaptation is a key driver of phenotypic and genetic divergence at loci responsible for adaptive traits variations in forest tree populations. Its experimental assessment requires rigorous sampling strategies such as those involving population pairs replicated across broad spatial scales.MethodsA hierarchical Bayesian model of selection (HBM) that explicitly considers both the replication of the environmental contrast and the hierarchical genetic structure among replicated study sites is introduced. Its power was assessed through simulations and compared to classical ‘within-site’ approaches (FDIST, BAYESCAN) and a simplified, within-site, version of the model introduced here (SBM).ResultsHBM demonstrates that hierarchical approaches are very powerful to detect replicated patterns of adaptive divergence with low false-discovery (FDR) and false-non-discovery (FNR) rates compared to the analysis of different sites separately through within-site approaches. The hypothesis of local adaptation to altitude was further addressed by analyzing replicated Abies alba population pairs (low and high elevations) across the species’ southern distribution range, where the effects of climatic selection are expected to be the strongest. For comparison, a single population pair from the closely related species A. cephalonica was also analyzed. The hierarchical model did not detect any pattern of adaptive divergence to altitude replicated in the different study sites. Instead, idiosyncratic patterns of local adaptation among sites were detected by within-site approaches.ConclusionHierarchical approaches may miss idiosyncratic patterns of adaptation among sites, and we strongly recommend the use of both hierarchical (multi-site) and classical (within-site) approaches when addressing the question of adaptation across broad spatial scales.

The study of the dispersal capability of a species can provide essential information for the management and conservation of its genetic variability. Comparison of gene flow rates among populations characterized by different management and evolutionary histories allows one to decipher the role of factors such as isolation and tree density on gene movements. We used two paternity analysis approaches and different strategies to handle the possible presence of genotyping errors to obtain robust estimates of pollen flow in four European beech (Fagus sylvatica L.) populations from Austria and France. In each country one of the two plots is located in an unmanaged forest; the other plots are managed with a shelterwood system and inside a colonization area (in Austria and France, respectively). The two paternity analysis approaches provided almost identical estimates of gene flow. In general, we found high pollen immigration (B75% of pollen from outside), with the exception of the plot from a highly isolated forest remnant (B50%). In the two unmanaged plots, the average within-population pollen dispersal distances (from 80 to 184 m) were higher than previously estimated for beech. From the comparison between the Austrian managed and unmanaged plots, that are only 500 m apart, we found no evidence that either gene flow or reproductive success distributions were significantly altered by forest management. The investigated phenotypic traits (crown area, height, diameter and flowering phenology) were not significantly related with male reproductive success. Shelterwood seems to have an effect on the distribution of within-population pollen dispersal distances. In the managed plot, pollen dispersal distances were shorter, possibly because adult tree density is three-fold (163 versus 57 trees per hectare) with respect to the unmanaged one.

At treeline, selection by harsh environmental conditions sets an upward limit to arboreal vegetation. Increasing temperatures and the decline of traditional animal raising have favoured an upward shift of treeline in the last decades. These circumstances create a unique opportunity to study the balance of the main forces (selection and gene flow) that drive tree migration. We conducted a parentage analysis sampling and genotyping with five microsatellite markers in all Norway spruce individuals (342 juveniles and 23 adults) found in a recently colonized treeline area (Paneveggio forest, Eastern Alps, Italy). Our goal was to evaluate local reproductive success versus gene flow from the outside. We were able to identify both parents among local adults for only 11.1% of the juveniles. In the gamete pool we sampled, two-thirds were not produced locally. Effective seed dispersal distance distribution was characterized by a peak far from the seed source (mean 344.66 m ± 191.02 s.d.). Reproductive success was skewed, with six local adults that generated almost two-thirds (62.4%) of juveniles with local parents. Our findings indicate that, although a few local adults seem to play an important role in the colonization process at treeline, large levels of gene flow from outside were maintained, suggesting that the potential advantages of local adults (such as local adaptation, proximity to the colonization area, phenological synchrony) did not prevent a large gamete immigration.

Understanding the genetic structure and diversity of edge populations can shed light on the role of peripheral populations and their relevance for conservation strategies. In this study, three fragmented and isolated Apennine populations of Scots pine (Pinus sylvestris L.) belonging to the rear edge of the species' distribution were analyzed using both nuclear (nSSR) and chloroplast (cpSSR) microsatellites and were compared with an Alpine population belonging to the species' main range. Although small population size and considerable isolation have probably reduced the genetic variability of Apennine populations, these fragmented populations maintain a high level of within-population genetic diversity. A signiWcant amongpopulation diVerentiation was found using both nSSR (F ST = 0.08) and cpSSR markers ( = 0.14). Analysis of molecular variance (AMOVA) on the nSSRs attributed all variabilities to the among Apennine populations component supporting the theoretical predictions regarding fragmentation eVects on genetic structure. On the other hand, AMOVA on the cpSSRs attributed all variances to the between-region component and no diVerentiation was found within region, among the Apennine populations. This result suggests the importance of pollen gene Xow in homogenizing populations on this geographical scale. Our results conWrm the genetic distinctiveness of Apennine populations and their possible derivation from diVerent glacial refugia than those of the Alps. Considering their peculiarity and the high level of intrapopulation genetic diversity that they still retain, fragmented Apennine populations should be considered of high priority for conservation.