Lobo A, Hansen JK, Hansen LN, Kjær ED. Differences among six woody perennials native to Northern Europe in their level of genetic differentiation and adaptive potential at fine local scale. Ecol Evol. 2018;8:2231--2239. <https://doi.org/10.1002/ece3.3824>

1. INTRODUCTION {#ece33824-sec-0001}
===============

Phenology is important for fitness of perennial species. Higher spring temperatures are expected to prolong the growing seasons due to earlier bud burst (Dragoni et al., [2011](#ece33824-bib-0010){ref-type="ref"}; Menzel & Fabian, [1999](#ece33824-bib-0034){ref-type="ref"}; Menzel et al., [2006](#ece33824-bib-0035){ref-type="ref"}; Richardson et al., [2010](#ece33824-bib-0045){ref-type="ref"}, [2013](#ece33824-bib-0046){ref-type="ref"}; Vitasse, Delzon, Dufrêne, et al., [2009](#ece33824-bib-0059){ref-type="ref"}), while lack of synchrony between phenology and occurrence of spring frost events increases risk of damage to early flushing plants (Duputié, Rutschmann, Ronce, & Chuine, [2015](#ece33824-bib-0012){ref-type="ref"}). Yet again, raised temperatures with no frost events and changes in daily temperature fluctuations could also result in later bud burst and influence time between bud development and growth cessation (Chmura et al., [2011](#ece33824-bib-0008){ref-type="ref"}; Rohde, Bastien, & Boerjan, [2011](#ece33824-bib-0048){ref-type="ref"}; Way, [2011](#ece33824-bib-0060){ref-type="ref"}). Several studies have documented substantial variation among populations in their phenology reflecting their geographic origin (Alberto, Derory, et al., [2013](#ece33824-bib-0004){ref-type="ref"}; Chuine & Beaubien, [2001](#ece33824-bib-0009){ref-type="ref"}; Salmela, Cavers, Cottrell, Iason, & Ennos, [2013](#ece33824-bib-0049){ref-type="ref"}). It is therefore fair to assume that genetic variation within species combined with divergent selection has played an important role for the ability of many tree species to thrive in a large distribution range across strong environmental gradients (Brousseau et al., [2016](#ece33824-bib-0006){ref-type="ref"}; Kawecki & Ebert, [2004](#ece33824-bib-0023){ref-type="ref"}; Pluess et al., [2016](#ece33824-bib-0042){ref-type="ref"}). However, species differ in their response and spring and autumn phenology may show different patterns. A study based on 59 tree species from similar climatic clines thus showed a relative clear pattern with respect to bud set in autumn but less clear pattern in spring bud burst (Alberto, Aitken, et al., [2013](#ece33824-bib-0003){ref-type="ref"}). The expectation of local adaptation often leads restoration and conservation programs to focus on gene pools at the local scale rather than the regional scale (Stanturf et al., [2015](#ece33824-bib-0056){ref-type="ref"}), but only few studies have actually assessed the variation at the local scale (within \~100 km).

Fast climate change calls for high phenotypic plasticity (Valladares et al., [2014](#ece33824-bib-0057){ref-type="ref"}) and genetic variation in relevant traits in plant populations (Franks, Weber, & Aitken, [2014](#ece33824-bib-0014){ref-type="ref"}) in order to ensure continued adaptation. Many tree species have revealed phenotypic plasticity and demonstrated their ability to thrive and reproduce when transferred to new areas (Kjær, Lobo, & Myking, [2014](#ece33824-bib-0025){ref-type="ref"}; Lascoux, Glémin, & Savolainen, [2016](#ece33824-bib-0033){ref-type="ref"}; Myking, Rusanen, Steffenrem, Kjær, & Jansson, [2016](#ece33824-bib-0036){ref-type="ref"}; Rehfeldt, Ying, Spittlehouse, & Hamilton, [1999](#ece33824-bib-0043){ref-type="ref"}). However, there are indications of declining plasticity in phenology in spring with continued global warming (Fu et al., [2015](#ece33824-bib-0015){ref-type="ref"}) and increased frequency of warm periods in the early spring may lead to frost backlashes (Jönsson & Bärring, [2011](#ece33824-bib-0022){ref-type="ref"}). The temperatures are likely to change in the future, but the seasonal variations in day length remain the same. Species will therefore respond differently as regards phenology to temperature changes depending on the species‐specific importance of day length versus heat sum in controlling flushing time (Alberto, Aitken, et al., [2013](#ece33824-bib-0003){ref-type="ref"}; Caffarra & Donnelly, [2011](#ece33824-bib-0007){ref-type="ref"}; Fu et al., [2015](#ece33824-bib-0015){ref-type="ref"}; Vitasse, Delzon, Bresson, Michalet, & Kremer, [2009](#ece33824-bib-0058){ref-type="ref"}; Vitasse, Delzon, Dufrêne, et al., [2009](#ece33824-bib-0059){ref-type="ref"}).

Variation among individual genotypes in their phenology triggers selection at the population level in favor of the new conditions, if the variation is in fitness traits and expressed with moderate or high heritability (Alberto, Aitken, et al., [2013](#ece33824-bib-0003){ref-type="ref"}; Falconer & Mackay, [1996](#ece33824-bib-0013){ref-type="ref"}). Genetic differentiation within and among population is therefore important for the species ability to adapt to climate change (Kubisch, Holt, Poethke, & Fronhofer, [2014](#ece33824-bib-0029){ref-type="ref"}). However, natural selection is not the only evolutionary force, because the combined actions of natural selection, gene flow and genetic drift drive the level of genetic variation within and among populations in fitness traits (Nadeau, Meirmans, Aitken, Ritland, & Isabel, [2016](#ece33824-bib-0037){ref-type="ref"}). The realized patterns therefore depend on effective population sizes, force of natural selection and limitations in pollen and seed dispersal (Aitken & Whitlock, [2013](#ece33824-bib-0002){ref-type="ref"}; Savolainen, Pyhäjärvi, & Knürr, [2007](#ece33824-bib-0052){ref-type="ref"}). Natural selection tends to increase the level of genetic variation among populations, but reduce the variation within populations leading to local adaptation. Fragmentation into small populations can create drift that also increase variation among population variation, but without supporting local adaptation. Small populations can rather reduce fitness due to inbreeding, at least in outcrossing species. Pollination and seed dispersal across the landscape among populations will counteract the fragmentation, but seed and pollen dispersal across environmental gradients will also counteract natural selection and thereby reduce the ability of populations to adapt to specific local conditions (Kubisch et al., [2014](#ece33824-bib-0029){ref-type="ref"}). The level of gene flow is expected to differ among woody species depending on occurrence, landscape features, mating system, pollen, and seed dispersal vectors. We therefore expect that the degree of local differentiation---and the degree to which this reflect local adaptation---will differ among woody plant species depending on their major life‐history traits.

1.1. Objectives of this study {#ece33824-sec-0002}
-----------------------------

On the above background, the objectives of this study were to use data on spring phenology to Assess the degree to which woody perennials reveal fine‐scale genetic differences among geographically and climatic close (differences in mean annual temperature around 1°C) populations, and test whether differences (if any) influence the fitness revealing signs of local adaptation.Estimate levels of genetic variation within and among populations in spring phenology and compare the levels of population differentiation (*Q* ~ST~) and additive genetic variation (*V* ~A~) among species.

Based on the results, we discuss the implications for dynamic conservation of genetic resources and the ability of woody plant species to adapt to future climate conditions.

2. MATERIALS AND METHODS {#ece33824-sec-0003}
========================

Six broadleaved species in Denmark were studied, of which three are insect pollinated---Common dogwood (*Cornus sanguinea* L.), European Crab Apple (*Malus sylvestris* (L.) Mill.) and Glaucous Dog Rose (*Rosa dumalis* (Bechst) Boulay), while three are wind pollinated---Downy birch (*Betula pubescens* Ehrh.), Common Hazel (*Corylus avellana* L.) and Sessile oak (*Quercus petraea* (Matt.) Liebl.). We used data from progeny trials with offspring from open pollinated trees based on seed collected from putative natural populations across Denmark (Figure [1](#ece33824-fig-0001){ref-type="fig"}). Description of the field trials, location of populations, and climate information at the population site and information on families within each population are provided in Appendix [S1](#ece33824-sup-0001){ref-type="supplementary-material"}.

![Locations of original populations for the different species in the study are spread across different climate zones (here represented by temperature minimum in May \[*T* ~min~\])](ECE3-8-2231-g001){#ece33824-fig-0001}

2.1. Assessment of phenotypic data and estimation of variance components {#ece33824-sec-0004}
------------------------------------------------------------------------

We used spring phenology assessed as bud burst scores using a scale from 0 to 8, where class 0 was closed winter buds and class 8 was fully opened leaves. The time of assessment was adapted for each species so that the trees in the trials were approximately in the middle of their flushing time (average values of bud burst score between 3 and 5). This was to ensure maximum variation among individuals in their bud burst ranking. Bud burst for *B. pubescens, C. sanguinea, M. sylvestris,* and *R. dumalis* were assessed in spring 2010, while assessments of *C. avellana* and *Q. petraea* were from spring 2013. Height was measured in the year of phenology assessment except for *R. dumalis*.

All analyzes were made per species, because they were tested in different common garden trials. The variation among trees in their phenological stage was separated into variance and covariance components reflecting genetic differences among and within populations as implemented in the software ASReml (Gilmour, Gogel, Cullis, & Thompson, [2009](#ece33824-bib-0017){ref-type="ref"}) from the model described below.

Population values were predicted by application of mixed model (1), and genetic variance component corresponding to each of the random effects in the model was estimated from the analysis; $$Y_{\mathit{ijkl}} = \mu + B_{i} + P_{j} + \lambda_{\mathit{ij}} + f_{k(j)} + \tau_{\mathit{ikj}} + \varepsilon_{\mathit{ijkl}}$$

where *Y* ~*ijkl*~ is the trait measured for tree *l*, μ is the overall mean of the trait, *B* ~*i*~ is the fixed block effect, *P* ~*j*~ is the fixed population effect, λ~*ij*~ is the fixed population by block interaction, *f* ~*k*(*j*)~ is the random effect of family within population, τ~*ijk*~ is the random effect of plot *k*, and ε~*ijkl*~ is the residual. *Q* ~ST~ values were estimated from model (1) as well having population effects as random.

The significance of the genetic variances of traits within sites and genotype by environment interaction across sites was tested using the log likelihood ratio (Gilmour et al., [2009](#ece33824-bib-0017){ref-type="ref"}). The significance of the populations was tested using the Satterthwaite approximation (Satterthwaite, [1946](#ece33824-bib-0051){ref-type="ref"}) in the procedure GLM in the statistical software program SAS (SAS Institute Inc. [2011](#ece33824-bib-0050){ref-type="ref"}).

2.2. Quantitative genetic analysis {#ece33824-sec-0005}
----------------------------------

We estimated the additive genetic variance (*V* ~A~) as 4σ~f~ ^2^ and narrow sense heritability within sites as; $\binom{\frown}{h}^{2} = 4\sigma_{f}^{2}/{(\sigma_{f}^{2} + \sigma_{\tau}^{2} + \sigma_{\varepsilon}^{2})}$, where σ~f~ ^2^ is the estimated family variance, σ~τ~ ^2^ is the estimated plot variance, and σ~ε~ ^2^ is the estimated within plot variance. The phenotypic variance was estimated as $\binom{\frown}{V}_{P} = \sigma_{f}^{2} + \sigma_{\tau}^{2} + \sigma_{\varepsilon}^{2}$. Families were considered to represent groups of half‐sibs. This assumption will overestimate the additive genetic variance, heritability, and expected response to selection, if the offspring within progeny groups on average are more related than half‐sibs (Falconer & Mackay, [1996](#ece33824-bib-0013){ref-type="ref"}). While the assumption of half‐sibs may provide a good fit for most trees species (Kjær, McKinney, Nielsen, Hansen, & Hansen, [2012](#ece33824-bib-0026){ref-type="ref"}; Larsen & Kjær, [2009](#ece33824-bib-0032){ref-type="ref"}), the situation may be more complicated for *R. dumalis,* where polyploidy may be present as the variance between families is biased due to a fraction of the dominance genetic variance, even if the families consist of pure half‐sibs (Roberts, Gladis, & Brumme, [2009](#ece33824-bib-0047){ref-type="ref"}).

The degree of differentiation among populations was estimated as the *Q* ~ST~ values (Spitze, [1993](#ece33824-bib-0055){ref-type="ref"}) using the formula;$$\binom{\frown}{Q}_{ST} = \binom{\frown}{V}_{POP}/{(\binom{\frown}{V}_{POP} + 2\binom{\frown}{V}_{A})}$$

where $\binom{\frown}{V}_{POP}$ is the estimated variance between populations, and $\binom{\frown}{V}_{A}$ is the estimated additive genetic variance.

Variance components, narrow sense heritability estimates, *Q* ~ST~ estimates, and their standard errors were estimated from model (1) using the software ASReml (Gilmour et al., [2009](#ece33824-bib-0017){ref-type="ref"}).

2.3. Support to the hypothesis of location adaptation {#ece33824-sec-0006}
-----------------------------------------------------

Climate data for the locations of origins were estimated with the ClimateEU v4.63 software package following the methodology described by (Hamann, Wang, Spittlehouse, & Murdock, [2013](#ece33824-bib-0019){ref-type="ref"}). Weighted regression analysis to test relationship between phenology of populations and climate at the population sites was carried out using procedure REG in SAS (SAS Institute Inc. [2011](#ece33824-bib-0050){ref-type="ref"}) having the predicted values of populations as dependent variable and climate variables as explanatory variables. One divided by the error variance for the predicted population values were used as weights. The climate variables tested for bud burst were monthly minimum temperatures and differences between monthly maximum and minimum temperatures in March, April, May and June, because these data are expected to provide good proxies for the risk of frost occurring after flushing. A backward selection approach (5% level) was used for the selection of climate variables that could best explain the variation in flushing time. The relationship between geographic distances and difference in average budburst was tested by a Mantel test as implemented in R version 3.2.2 (Dray & Dufour, [2007](#ece33824-bib-0011){ref-type="ref"}).

2.4. Fitness effects of spring phenology {#ece33824-sec-0007}
----------------------------------------

Genetic correlations between bud burst and height growth were estimated based on individual tree data to assess whether the height as a proxy for realized fitness in the present climate varied with the phenology. Genetic correlations could be estimated between survival and phenology based on plot means, because the plants were grown in family plots in all trials (see Appendix [S1](#ece33824-sup-0001){ref-type="supplementary-material"}). Genetic correlations were estimated according to Falconer and Mackay ([1996](#ece33824-bib-0013){ref-type="ref"}) as $\binom{\frown}{r}_{g} = \frac{\sigma_{\mathit{fxy}}}{\sqrt{\sigma_{\mathit{fx}}^{2}\sigma_{\mathit{fy}}^{2}}}$ , where σ~*fxy*~ is the family (within population) covariance between trait *x* and *y*, σ~*fx*~ ^2^ is the family (within population) variance for trait *x*, and σ~*fy*~ ^2^ is the family (within population) variance for trait *y*.

3. RESULTS {#ece33824-sec-0008}
==========

3.1. Presence of genetic variation within and among populations in bud burst {#ece33824-sec-0009}
----------------------------------------------------------------------------

Genetic variation both within and among populations was significant in all six species, but with large differences among the species. The additive genetic variance for bud burst (*V* ~A~) thus ranged from 0.07 in *C. sanguinea* to 0.34 in *Q. petraea* (Table [1](#ece33824-tbl-0001){ref-type="table-wrap"}), while the level of differentiation among populations (*Q* ~ST~) ranged from 0.04 in *M. sylvestris* to 0.27 in *C. sanguinea* (Table [2](#ece33824-tbl-0002){ref-type="table-wrap"}). The three wind‐pollinated species *B. pubescens, C. avellana* and *Q. petraea* showed the highest level of additive variance, while the insect‐pollinated species *C. sanguinea* and *R. dumalis* showed the highest level of population differentiation as measured by the *Q* ~ST~ (Figure [2](#ece33824-fig-0002){ref-type="fig"}).

###### 

Family variation and genetic parameters for bud burst in different species

  Species              Trait                     Mean   *p*‐Value family   *V* ~A~   *V* ~P~   *h* ^2^   *SE*
  -------------------- ------------------------- ------ ------------------ --------- --------- --------- ------
  *Betula pubescens*   Bud burst 24 April 2010   3.27   .003               0.15      0.30      0.50      0.12
  *Cornus sanguinea*   Bud burst 6 May 2010      3.81   \<.0001            0.07      0.27      0.24      0.11
  *Corylus avellana*   Budburst 24 April 2013    3.09   \<.0001            0.13      0.27      0.47      0.09
  *Malus sylvestris*   Bud burst 28 April 2010   4.05   \<.0001            0.03      0.10      0.31      0.08
  *Quercus petraea*    Bud burst 24 May 2013     4.64   \<.0001            0.34      0.58      0.58      0.26
  *Rosa dumalis*       Bud burst 26 April 2010   3.58   \<.0001            0.09      0.21      0.41      0.14

*V* ~A~ = additive genetic variance, *V* ~P~ = *p* variance, *h* ^2^ = narrow sense heritability, *SE* = standard error for *h* ^2^.

John Wiley & Sons, Ltd

###### 

Level of population differentiation for bud burst in different species

  Species              Trait                     *Q* ~st~   *SE*     *p*‐Value population   Max. pop. diff.[a](#ece33824-note-0002){ref-type="fn"}
  -------------------- ------------------------- ---------- -------- ---------------------- --------------------------------------------------------
  *Betula pubescens*   Bud burst 24 April 2010   0.13       0.08     \<.001                 0.71
  *Cornus sanguinea*   Bud burst 6 May 2010      0.27       *0.14*   \<.001                 0.88
  *Corylus avellana*   Budburst 24 April 2013    *0.12*     *0.05*   \<.0001                0.86
  *Malus sylvestris*   Bud burst 28 April 2010   0.04       *0.03*   .02                    0.28
  *Quercus petraea*    Bud burst 24 May 2013     0.13       0.09     .0003                  2.45
  *Rosa dumalis*       Bud burst 26 April 2010   0.18       0.09     \<.001                 1.52

Maximum difference between populations in bud burst score.

QST = quantitative genetic differentiation among populations in bud burst, SE = standard error for QST

John Wiley & Sons, Ltd

![Difference among species in genetic variation within populations (*V*~A~ = additive genetic variance) and population differentiation (*Q*~ST~ value) for bud ![](ECE3-8-2231-g005.jpg "image") burst insect‐pollinated species, ![](ECE3-8-2231-g006.jpg "image") wind‐pollinated species](ECE3-8-2231-g002){#ece33824-fig-0002}

3.2. Relationship of phenology with growth, fitness, distance between populations and climate at original population site {#ece33824-sec-0010}
-------------------------------------------------------------------------------------------------------------------------

The regression between bud burst of populations and population site minimum temperatures in March were significant (*p \< .05*), or close to significant (*p = .05*) for *B. pubescens* and *C. avellana* (Figure [3](#ece33824-fig-0003){ref-type="fig"}), and the regression between bud burst of populations and population site minimum temperatures in May was close to significant (*p = .05*) for *Q. petraea* (Figure [3](#ece33824-fig-0003){ref-type="fig"}). This corresponds to the general time for bud burst, where *B. pubescens* and *C. avellana* flush before *Q. petraea* in Denmark and indicate that spring frost event has played an important role in developing the pattern of variation among population in these species. Similar significant correlations were not observed in any of the three insect‐pollinated species. Mantel tests for correlation between geographic distance and difference in bud burst at population were nonsignificant in all the species (Figure [4](#ece33824-fig-0004){ref-type="fig"}). Survival was not significantly correlated to bud burst data in any of the species evaluated, and a significant additive genetic correlation between height and budburst was only observed in *B. pubescens (rg = *0.35; *SE 0.23*).

![Weighted regression between population means of bud burst score and *T* ~min~ (temperature minimum) of different months in spring at original population site in different species](ECE3-8-2231-g003){#ece33824-fig-0003}

![Population genetic differentiation in bud burst show no relationship with geographic distance (in 100 km) between them in different species](ECE3-8-2231-g004){#ece33824-fig-0004}

4. DISCUSSION {#ece33824-sec-0011}
=============

Our results show that fine‐scale local genetic differentiation in fitness traits such as phenology indeed can be present and to a larger extent than previously anticipated in the Danish gene conservation strategy (Graudal, Kjær, & Canger, [1995](#ece33824-bib-0018){ref-type="ref"}). However, the analysis across the six different species did not provide unique patterns of local adaptation, as the absence of significant regression between bud burst time and spring temperature in the three insect‐pollinated species suggests that the variation among population is not always simply reflecting local adaptation. It is rather more likely a result of natural selection and neutral processes simultaneously acting as predicted by theory (Nadeau et al., [2016](#ece33824-bib-0037){ref-type="ref"}).The fact that populations in all studied six species were significantly differentiated in their spring phenology suggests that local populations of woody perennials more often than not are genetically differentiated, and this can be the case even if they only are separated by 10 to 35 km (Figure [4](#ece33824-fig-0004){ref-type="fig"}) and growing in areas with low variation in altitude and where the spring temperature varies only between 1°C and 2.5°C (Figure [3](#ece33824-fig-0003){ref-type="fig"}). This finding should be compared to the prediction of Northern Europe becoming 2--4°C warmer by the year 2100 (Nikulin, Kjellström, Hansson, Strandberg, & Ullerstig, [2011](#ece33824-bib-0039){ref-type="ref"}). However, the observed high level of additive genetic variation in the studied species indicates that the native populations of woody species do possess substantial evolutionary potential. Divergent natural selection is possible within few generations, if the timing of phenology creates large differences in fitness. This is supported by the fact that the observed patterns in our study for *B. pubescens*,*C. avellana,* and *Q. petraea* (Figure [3](#ece33824-fig-0003){ref-type="fig"}) are at least partly due to past natural selection in favor of local adaption.

The presence of local adaptation and evolutionary potential in woody species are often attributed to their life‐history traits including pollination mechanisms and associated genetic processes, and it is therefore interesting that the patterns were different for the wind‐ and insect‐pollinated species in our study. The result that the three wind‐pollinated, continuous distributed species (*B. pubescens, C. avellana, Q. petraea)* were characterized by having high levels of genetic variation within populations (right circle of Figure [2](#ece33824-fig-0002){ref-type="fig"}) supports the general assumption that wind‐pollinated, widely distributed species have allocated genetic variation within populations to a larger extent than insect‐pollinated species with small and scattered populations (Petit & Hampe, [2006](#ece33824-bib-0041){ref-type="ref"}; Smith & Donoghue, [2008](#ece33824-bib-0053){ref-type="ref"}). The results correspond to the expectation that long‐distance gene flow through wind pollination will maintain connectivity among populations (lowering *Q* ~ST~) and counteract loss of genetic variation (maintaining high *h* ^2^) within populations (Sork, [2016](#ece33824-bib-0054){ref-type="ref"}). High levels of gene flow among populations across landscapes have been reported in the wind‐pollinated *Quercus* (Gerber et al., [2014](#ece33824-bib-0016){ref-type="ref"}). The insect‐pollinated *C. sanguinea* was in the opposite corner of *Q. petraea* in Figure [2](#ece33824-fig-0002){ref-type="fig"}, reflecting a much higher proportion of genetic variation located among populations. *C. sanguinea* in general occurs in scattered and small populations in Denmark (Ødum, [1968](#ece33824-bib-0040){ref-type="ref"}). Seed dispersal by birds (Krüsi & Debussche, [1988](#ece33824-bib-0028){ref-type="ref"}) may generate long‐distance gene flow across the landscape, but the observed level and distribution of genetic diversity in the present study point toward genetic drift rather than natural selection as important driver behind the observed population differentiation. The true difference between *Q. petraea* and *C. sanguinea* in level of genetic variation may be even higher, because the assumption of half‐sib relationship within seed from single tree collection is likely to create a bias toward overestimation of *V* ~A~ in the rose as discussed in M&M above. *M. sylvestris* is again a different case in Figure [2](#ece33824-fig-0002){ref-type="fig"}, with low level of genetic differentiation but also relatively low additive genetic variance. Pollinating insects have been reported to visit primarily trees in immediate vicinity, but long‐distance pollination events are also likely to occur (Larsen & Kjær, [2009](#ece33824-bib-0032){ref-type="ref"}; Reim et al., [2015](#ece33824-bib-0044){ref-type="ref"}) and seeds moved by birds and deer that feed on the wild apples also create long‐distance gene flow (Larsen & Kjær, [2009](#ece33824-bib-0032){ref-type="ref"}). Studies based on neutral SSR markers revealed low level of population differences within Denmark (*F* ~ST~ = 0.03) (Larsen, Asmussen, Coart, Olrik, & Kjær, [2006](#ece33824-bib-0031){ref-type="ref"}), which is very close to the estimated *Q* ~ST~ = 0.04 for flushing time in the present. Therefore, the observed variation among populations in *M. sylvestris* may have been generated by random processes without significant effects of divergent selection for flushing time.

The studied populations have in general maintained a high level of genetic variation and thereby possess a high ability to adapt to new conditions if exposed to natural selection in favor of---or against---early flushing. The question is how fast species can adapt under the anticipated rapid increase in temperature over the next 100 years, and how important such changes are for the fitness of the species. In *B. pubescens*, the maximum difference in minimum temperature in March of the populations tested in this study was 2.4°C (Figure [3](#ece33824-fig-0003){ref-type="fig"}). Even with a conservative estimate of a generation time of 50 years (Hynynen et al., [2010](#ece33824-bib-0021){ref-type="ref"}), the results suggest that *B. pubescens* possesses the resilience at the population level to change phenology corresponding to the predicted increase in temperature. A similar conclusion can probably be drawn for *Q. petraea* also (Figure [3](#ece33824-fig-0003){ref-type="fig"}). For many of the other studied species, the generation time is possibly less than 20 years (e.g., *R. dumalis*,*M. sylvestris*,*C. avellana,* and *C. sanguinea*), and for these species, it therefore also seems reasonable to assume that populations can evolve in spring phenology with a speed that can match the predicted increase in temperature. The response will only be realized if early bud burst has substantial influence on fitness (survival and reproduction) or because of directional selection implemented in domestication programmes. The estimates in the present study refer to the latter situation, because the heritability was estimated in managed progeny test and the realized heritability in natural population may be substantially lower due to higher environmental and developmental heterogeneity. In Denmark, all major woody plant species are included in domestication programs based on breeding seed orchards in order to support effective and rapid selection for continued fitness (Kjær et al., [2009](#ece33824-bib-0024){ref-type="ref"}).

Our finding on genetic differences among geographically close populations (cf. Table [2](#ece33824-tbl-0002){ref-type="table-wrap"}) supports that selectively divergent phenotypes can be maintained even if there are options for gene flow among these populations as suggested by Aitken and Bemmels ([2016](#ece33824-bib-0001){ref-type="ref"}). Kremer et al. ([2012](#ece33824-bib-0027){ref-type="ref"}) argue that the positive effect of gene flow on adaptation of trees to climate change is larger than the negative effect due to the introduction of genetic variation in adaptive traits. Species are different with respect to their genetic differentiation within and among their populations as revealed in our study. Assisted migration for genetic enrichment can serve as a supplement to genetic management based on native populations; especially, the species associated with small population sizes and limited gene flow. Populations at the low latitudinal limit of the species range have in general maintained high biological diversity over time as shown for, for example, *Abies alba* (Bergmann, Gregorius, & Larsen, [1990](#ece33824-bib-0005){ref-type="ref"}; Larsen, [1986](#ece33824-bib-0030){ref-type="ref"}) and southern populations should therefore be considered as a source of gene pool for assisted migration of genotypes under the future climate predictions (Hampe & Petit, [2005](#ece33824-bib-0020){ref-type="ref"}). The populations in the present study are close to the northern distribution range of the species (except for *B. pubescens*) and genetic diversity may therefore be lower compared to populations closer to the refugial areas. But we still observed substantial level of genetic variation in the fitness‐related trait. For species occurring in extremely small fragmented populations assisted migration at a more localized landscape level may still be desirable to counteract inbreeding within these populations due to random drift. Adaptive potential of a species is not only determined by the presence of genetic variation within and among populations as studied here, but also by the presence of phenotypic plasticity in adaptive traits (Nicotra et al., [2010](#ece33824-bib-0038){ref-type="ref"}). Final conclusions on the adaptive potential of species to climate change will therefore also depend on the role of plasticity in local adaptation in the studied species.

DATA ACCESSIBILITY {#ece33824-sec-0013}
==================

Data for this study is available at University of Copenhagen---Electronic Research Data Archive (UCPH‐ERDA).

CONFLICT OF INTEREST {#ece33824-sec-0014}
====================

None declared.

AUTHOR CONTRIBUTION {#ece33824-sec-0015}
===================

Albin Lobo was responsible for the design of study, analysis and interpretation of data, and writing of the manuscript. Jon Kehlet Hansen was responsible for supervision of the work, design of study, and writing and approving of manuscript. Lars Nørgaard Hansen was responsible for field data collection, and revising and approving the final version of the manuscript before submission. Erik Dahl Kjær was responsible for overall supervision of the study, conception and design of work, data interpretation, and writing and final approval of the manuscript.

Supporting information
======================

###### 

 

###### 

Click here for additional data file.

We would like to thank the Danish Nature Agency for providing access to all trials and supporting the huge effort of data collection. The authors would also like to thank the Villum Foundation for financial support for the data collection and the analysis as part of the Trees for future forests project (VKR‐023063). We also thank Dr. Lene Rostgaard Nielsen for the critical reading of this manuscript.
