Varying Levels of Genetic Control and Phenotypic Plasticity in Timing of Bud Burst, Flower Opening, Leaf Senescence and Leaf Fall in Two Common Gardens of Prunus padus L.

: Several phenological phases mark the seasonal growth pattern in temperate woody perennials. To gain further insight into the way these phases react on an altering growth environment, we tested whether vegetative and reproductive phenophases in a shrub species respond di ﬀ erentially among di ﬀ erent genetic entities and between two di ﬀ erent planting sites. We scored leaf bud burst, ﬂower opening, leaf senescence and leaf fall on 267 ramets of Prunus padus L. belonging to 53 genotypes that were sampled in 9 local populations, and that were planted in 2 common gardens in the northern part of Belgium. The data were processed with cumulative logistic regression. The contribution of genetic and non-genetic components to the total variability varied between the four studied seasonal phenophases. The timing of ﬂower opening displayed the smallest relative amount of intragenotypic variance (between ramets), suggesting a stronger genetic control and a lesser need at the individual plant level for plastic ﬁne tuning to the micro-environment. In addition, whereas leaf bud burst showed the highest relative variance at the interpopulation level among all phenophases, probably at least partly attributable to local adaptation, ﬂower opening displayed the highest intergenotypic variance, which may have been promoted more by assortative mating. Spring phenophases were strongly correlated ( r = 0.89) as well as the autumnal phenophases ( r = 0.72). Flower opening was not correlated with the autumnal phenophases. Timing of leaf bud burst and leaf senescence were negatively correlated, demonstrating that the length of the growing season enlarged or diminished among the studied genotypes. Although the two planting sites were only 24 km apart, all phenophases were advanced at the less exposed site, indicating a phenotypic plastic response. Together, our results suggest that in P. padus , ﬂowering is less sensitive to environmental variation than leaf bud break and may show a lesser impact of a changing environment on this reproductive phenophase.


Introduction
Understanding how plants respond to environmental change will help predict the impact of climate change on plant populations. This is especially true for woody species that are sessile and perennial organisms, have a long life-span, support diverse communities and give structure to a wide range of ecosystems on Earth [1]. An individual plant can maintain its fitness by altering its phenotype in response to a changing growth environment [2]. Trees experience considerable variation in growth conditions throughout their life [3,4]. Therefore, the ability of tree species for phenotypic plastic responses and to adapt locally to altering environments is studied [5]. Plasticity will likely play an autumnal leaf senescence is responding both to photoperiod and to temperature [23]. Still, there is no consensus on the exact factors that control leaf senescence [23]. Leaf fall is the next phenophase following leaf senescence. Similar to the spring phenophases bud burst and flowering, and more so as it concerns the same vegetative organ, the timing of leaf senescence and leaf abscission is most likely strongly auto-correlated. This auto-correlation is reported to be to such an extent that leaf fall can be used as a proxy for leaf senescence, taking into account a certain delay in time [26].
The quantitative genetic knowledge of lesser known or rare forest tree and shrub species remains scarce. We focused on Prunus padus L. (European bird cherry), which has a wide natural distribution area in Europe [27]. The species usually forms small populations in the undergrowth of riparian and moist forests [28]. Shrubby thickets are common in meadows. Furthermore, P. padus creates blooming edgings of forests. The species is insect pollinated and is characterized by endozoic dispersal [28]. Birds are essential disseminating agents. Fruits that fall beneath the shrubs may be dispersed by small mammals. To address the question how the shrub P. padus may respond to a changing environment, we combined two advantages related to this species. One is that it takes less time to achieve a reproductive life stage in comparison to tree species. Secondly, P. padus can easily be propagated vegetatively through cuttings, allowing an additional hierarchical level of relatedness among the individuals planted in a common garden, namely, the clonal level (further called ramets) next to the individual level (genotypes) and the populations. Planting of identical genotypes in clonal repetitions permits a detailed analysis of variance, accounting for both within genotype variation (non-genetic) and variation between genotypes and populations (genetic). Traditionally, half-sib offspring of mother trees, called families, are planted in common gardens of tree species. We adopted a more stringent method, making use of vegetatively propagated plants planted in two common gardens. We firstly hypothesized different levels of genetic control for the four phenophases leaf bud burst, flower opening, leaf senescence and leaf fall when compared with each other. Secondly, we questioned whether different responses could be detected for the four phenophases of clonal replicates at the two plantation sites, which indicates phenotypic plasticity. It should be noted that we focused on the relative differences between the four phenophases without any intention to model the absolute timing of the phenophases. The latter needs a different analysis including chilling and forcing temperatures and will be subject of another study.

Source Material and Planting Sites
Cuttings were taken in the spring of 2006 for the clonal propagation of 53 genotypes sampled in 9 local natural populations in the northern part of Belgium ( Figure 1, Table 1). These populations were classified as autochthonous in an inventory locating sites harbouring autochthonous populations of trees and shrubs in the northern part of Belgium [29]. As P. padus is known to propagate vegetatively by natural root suckers, sampling from the same genotype was avoided by evaluating the distance between two shrubs (at least 5 m without any P. padus shrub in between two sampled shrubs), and by evaluating differences in the habitus of the shrubs. Vigorously growing shoots (the season's new growth) on the mother shrubs in the natural populations were cut and taken to the greenhouse. There, shoots were cut with a grafting knife in pieces of 6 to 7 cm length, still having at least one leaf. The cuttings' rooting edges were wounded by stripping of a small part of the bark. No phytohormones were applied as P. padus can propagate easily by cuttings. Cuttings were placed in trays filled with moist standard seeding and cutting soil, with the rooting edge 2 cm into the soil. Thin transparent plastic covered the cuttings until they rooted. Rooted cuttings were transferred to plant pots and grown to 267 two-year-old potted planting stock in the nursery of the Research Institute for Nature and Forest in Geraardsbergen, Belgium, following standard nursery methods. In Dentergem and in Semmerzake, 147 plants (53 genotypes, on average 2.8 ramets per genotype) and 120 plants (52 genotypes, on average 2.3 ramets per genotype) were planted, respectively. Table 1. Genotypes of Prunus padus L. that were sampled in natural populations in the northern part of Belgium and vegetatively propagated, and the number of ramets that were planted in the two common gardens at Dentergem (n D ) and Semmerzake (n S ).

Brakel BR
The two planting sites were 24 km apart from each other ( Figure 1). For a modest approximation of the nutrient status of the soils at both sites, two random bulk soil samples were taken at each site from the Ah horizon of the mineral soil and dried at 40 • C (24 h). Soil pH (H 2 O) and electrical conductivity (EC) were measured at a 1:5 soil-to-water ratio. The site in Semmerzake (altitude of 7.5 m) was alluvial and appeared nutrient richer (pH values of 7.64 and 7.59; EC values of 127.6 µS/cm and 100.3 µS/cm) compared with the site in Dentergem (altitude of 24 m; pH values of 5.95 and 5.00; EC values of 75.8 and 75.1) which had a texture of loamy sand to sand and an average soil moisture content. The planting in Semmerzake was situated at the southwest of a mature forest. In Dentergem, the planted shrubs were more exposed as the site was located in an agricultural landscape lacking forests, wooded banks, hedges or tree rows that could protect against harsh weathers. As the two planting sites were situated in lowland area very close to each other, no essential mean temperature differences between the sites were expected. We checked this by analysing data from the weather stations nearest to the two common gardens (not on-site) [30]. The average mean yearly temperature between 1980 and 2014 in Semmerzake (weather station at 930 m from plantation site) was only 0.15 • C higher than Dentergem (weather station at 1877 m from plantation site, Figure S1). A paired t-test on the mean monthly temperatures between 1980 and 2014 between the two sites was significant (df = 12,418, p value < 0.001). Although significant, the mean temperature difference was considered negligible in terms of influencing factors on the timings of the studied phenophases. The shrubs were planted at both sites with a spacing of 3 × 3 m and individually mingled (single tree plot design).

Scoring of the Phenophases
Scoring of the phenophases was performed in 2017. The evolution of the four phenophases, namely, leaf bud burst, flower opening, leaf senescence and leaf fall, was recorded during the different visits to the plantations using four easy-to-use protocols that were specifically composed for this species ( Table 2). The recordings were done in both plantations on the same day by the same experienced researcher. For each phenophase, the whole shrub was observed (all buds or leaves were observed), and an average score was given.

Statistical Analysis of Phenological Traits
All statistical analyses were performed in the open source software R 3.6.1 [31]. Models were fit with each phenological trait (T ph ) as a response variable. The phenological scorings were in an ordinal scale and therefore modelled using cumulative logistic regression. The function "clmm" in the package ordinal models the probability (p) to have reached maximally a given score level of the response variable [32]. The score levels of bud burst and flower opening were defined in decreasing order, so that the probability to have reached maximally, e.g., a bud burst score of 3, comprised scores 5, 4 and 3. This allowed it to be understood as having reached at least a score of 3. The score levels for leaf senescence and fall were defined in an increasing order. Mixed models were fit. Location (L), designating the two planting sites Dentergem and Semmerzake, was in the fixed part of each model and accounted for putative deviating responses of the genetic entities at the two different planting sites. A significant location effect pointed therefore to a spatial phenotypic plastic response. Day (D) was also in the fixed part of each model as each phenophase was scored on several days. This variable therefore accounted for the fact that a given plant reached a further development of a given phenophase at a later time of observation. The random part (random intercept) consisted of population (P), genotype (G) and a unique shrub identity (ID). Formula of the applied models: α i is an intercept value between two succeeding score levels. β D and β L are the estimated coefficients for the fixed covariates D and L, respectively. R P , r G and r ID are the random effect coefficients for all the levels of the categorical variables P, G and ID, respectively. Variance components attributable to the variation between the populations, further called "interpopulation" (σ 2 P ), between the genotypes within a population, further called "intergenotype" (σ 2 G ) and to the variation within each genotype, which is the variance between ramets of a genotype, further called "intragenotype" (σ 2 ID ) were extracted from the models.
As several ramets per genotype were planted in each common garden, the variance attributable to this component in the models designated the variance due to the micro-environment or to soma-clonal variation, independent of any genetic contribution, whereas intergenotype and interpopulation variance expressed the genetically controlled variation. To examine the variability in timing of the four phenophases among the sampled populations, the four models were adjusted by moving the population variable (P) from the random part to the fixed part and by adding an interaction term between population (P) and location (L). Formula of the applied models: Based on these models, the DOYs (day of the year) were calculated when the probability for having reached at least a bud burst or flower opening score of 4 and having reached a maximally leaf senescence score of 2 and a leaf fall score of 3 attained 50% (D 50%P ). A D 50%P for a given population therefore indicated the day that half of the plants of this population had reached at least (bud burst and flower opening) or maximally (leaf senescence and fall) this stage of the respective phenophase. This calculation was based on the fact that log (p Tph /(1 − p Tph )) = 0 for p = 50%. For the plantation in Dentergem (Dentergem is the standard level in the model): For the plantation in Semmerzake: The formula for Dentergem was adapted for the calculation of D 50% for each genotype: The genotypic D 50%G values were used to calculate Pearson correlation coefficients among the four phenophases.

Variance Partitioning between Populations, Genotypes and Ramets
General linear mixed models were fit to the observational data to examine the timing of the phenophases leaf bud burst, flower opening, leaf senescence and leaf fall in the shrub P. padus in two common gardens ( Table 3). The estimated variances of the variables present in the random parts of the models allowed discriminating variance between ramets within genotypes, variance between genotypes within populations and variance between populations. As the location variable was accounted for in the fixed part of the models, results of the variance partitioning were independent of the different growth environments at the two planting sites. Relative interpopulation and intergenotypic variances represented the genetic contributions to the timing of the phenophases, whereas the intragenotypic variance denoted the non-genetic contribution. Among the different phenophases, flower opening displayed the relatively smallest non-genetic contribution (mainly environmentally induced, possibly also soma-clonal variability), and, thus, the relatively highest genetic control, as the intragenotypic variance, was smallest, followed by bud burst and subsequently by leaf senescence and finally leaf fall ( Figure 2). The relative interpopulation variance was higher than intergenotypic variance in bud burst, whereas the contributions of these two levels of genetically controlled variances differed relatively less between each other for flower opening, leaf senescence and leaf fall.

Figure 2.
Variance components extracted from the mixed models for the four phenophases observed in the two common gardens. Relative variance is partitioned for populations (between populations), genotypes (between genotypes within populations) and ramets (between ramets within genotypes).

Relative Timing of the Phenophases at the Two Planting Sites
A putative phenotypic plastic response in the phenology of the plants was approached by comparing the responses between the two planting sites. All phenophases showed a significant differentiation in timing between the two studied plantation sites (the variable location being significant in the four models in Table 3: p-values < 0,001 for bud burst, leaf senescence and leaf fall, and p = 0.005 for flower opening), with shrubs at the more exposed location in Dentergem being delayed in comparison with Semmerzake ( Figure 3). As the relative genetic (interpopulation and intergenotypic variance) and the non-genetic (intragenotypic variance) contributions within each location were accounted for in the random parts of the models, a significant location effect represented a different phenotype for the same genotype between the two sites. Variance components extracted from the mixed models for the four phenophases observed in the two common gardens. Relative variance is partitioned for populations (between populations), genotypes (between genotypes within populations) and ramets (between ramets within genotypes).

Relative Timing of the Phenophases at the Two Planting Sites
A putative phenotypic plastic response in the phenology of the plants was approached by comparing the responses between the two planting sites. All phenophases showed a significant differentiation in timing between the two studied plantation sites (the variable location being significant in the four models in Table 3: p-values < 0,001 for bud burst, leaf senescence and leaf fall, and p = 0.005 for flower opening), with shrubs at the more exposed location in Dentergem being delayed in comparison with Semmerzake ( Figure 3). As the relative genetic (interpopulation and intergenotypic variance) and the non-genetic (intragenotypic variance) contributions within each location were accounted for in the random parts of the models, a significant location effect represented a different phenotype for the same genotype between the two sites. Modelled timing of bud burst, flower opening, leaf senescence and leaf fall at the two common gardens Dentergem and Semmerzake. For bud burst and flower opening, the probability is displayed for having reached a score of 4 or higher. For leaf senescence, the probability is displayed for having reached a score of 2 or lower. For leaf fall, the probability is displayed for having reached a score of 3 or lower.

Duration and Order in Timing of the Phenophases at the Population and Genotypic Level
To assess the duration and the order in timing of leaf bud burst, flower opening, leaf senescence and leaf fall at the population and genotypic level, the four phenological models were fit with population as covariate in the fixed part in interaction with the plantation site ( Figure S2, Table S1). For each phenophase in each plantation, DOY values were calculated for which each population or genotype had reached the same phenological stage (D50%P and D50%G, respectively). For each phenophase, the time span between the two populations with the lowest and the highest D50%P were calculated by simple subtractions of the D50%P values. Spring phenophases develop more quickly than autumnal phenophases. Flower opening (9.5 days in Dentergem; 9.1 days in Semmerzake) and leaf fall (17.2 days in Dentergem; 22.6 days in Semmerzake) displayed shorter time spans between the first and the last population in comparison with bud burst (14.2 days in Dentergem; 14.8 days in Semmerzake) and leaf senescence (20.4 days in Dentergem; 24.1 days in Semmerzake), respectively. The D50%P values were standardized for each phenophase separately and are subsequently shown together in Figure 4 for Dentergem ( Figure S3 for Semmerzake). Clearly, the order in timing among the different populations was not maintained across the four phenophases, suggesting deviating influencing factors. Pearson correlation coefficients were calculated between the timing of the different phenophases based on the D50%G values ( Table 4). As expected, the spring phenophases, bud burst and flower opening, are strongly (and significantly) correlated (r = 0.89), as well as the autumnal phenophases, leaf senescence and fall (r = 0.72). Bud burst was negatively correlated (significantly) with leaf senescence (r = −0.38), indicating that early flushing genotypes enter leaf senescence later. Flower opening was not correlated with the autumnal phenophases.

Duration and Order in Timing of the Phenophases at the Population and Genotypic Level
To assess the duration and the order in timing of leaf bud burst, flower opening, leaf senescence and leaf fall at the population and genotypic level, the four phenological models were fit with population as covariate in the fixed part in interaction with the plantation site ( Figure S2, Table S1). For each phenophase in each plantation, DOY values were calculated for which each population or genotype had reached the same phenological stage (D 50%P and D 50%G , respectively). For each phenophase, the time span between the two populations with the lowest and the highest D 50%P were calculated by simple subtractions of the D 50%P values. Spring phenophases develop more quickly than autumnal phenophases. Flower opening (9.5 days in Dentergem; 9.1 days in Semmerzake) and leaf fall (17.2 days in Dentergem; 22.6 days in Semmerzake) displayed shorter time spans between the first and the last population in comparison with bud burst (14.2 days in Dentergem; 14.8 days in Semmerzake) and leaf senescence (20.4 days in Dentergem; 24.1 days in Semmerzake), respectively. The D 50%P values were standardized for each phenophase separately and are subsequently shown together in Figure 4 for Dentergem ( Figure S3 for Semmerzake). Clearly, the order in timing among the different populations was not maintained across the four phenophases, suggesting deviating influencing factors. Pearson correlation coefficients were calculated between the timing of the different phenophases based on the D 50%G values ( Table 4). As expected, the spring phenophases, bud burst and flower opening, are strongly (and significantly) correlated (r = 0.89), as well as the autumnal phenophases, leaf senescence and fall (r = 0.72). Bud burst was negatively correlated (significantly) with leaf senescence (r = −0.38), indicating that early flushing genotypes enter leaf senescence later. Flower opening was not correlated with the autumnal phenophases.  Table 1.

Timing of Leaf Bud Burst and Flower Opening
The modelling analysis of the seasonal phenophases observed on clonal replicates of P. padus in the two common garden plantations allowed partitioning of the variance attributable to genetic (between populations and between genotypes) and non-genetic (between ramets within genotypes) components, while accounting for the variability attributable to the two planting sites (location in fixed part of the models).
Bud burst displayed the highest interpopulation variation among all phenophases, indicating a relatively higher population differentiation. Population differentiation for bud burst in woody perennials is often encountered in common garden experiments and may be shaped, at least partly, by local adaptation. Habitat fragmentation may also have contributed to the observed population differentiation for bud burst in our study given the fact that Flanders, the northern part of Belgium, is already historically a very densely populated area. Habitat fragmentation affects genetic structures of tree populations including population divergence, as was, for instance, shown for Fagus sylvatica L. [33]. A stronger genetic control for bud burst in comparison with autumnal phenophases may be related to the fact that springtime leafing out modulates the timing of carbon assimilation in deciduous forests more than the autumnal phenophases. Early bud burst was associated with a longer period of carbon uptake and with an increase in net annual carbon fluxes [34]. The gain of 1 day extra photosynthesis in spring is more beneficial than an extra day in autumn when days are shorter and trees capture less solar radiation.  Table 1.

Timing of Leaf Bud Burst and Flower Opening
The modelling analysis of the seasonal phenophases observed on clonal replicates of P. padus in the two common garden plantations allowed partitioning of the variance attributable to genetic (between populations and between genotypes) and non-genetic (between ramets within genotypes) components, while accounting for the variability attributable to the two planting sites (location in fixed part of the models).
Bud burst displayed the highest interpopulation variation among all phenophases, indicating a relatively higher population differentiation. Population differentiation for bud burst in woody perennials is often encountered in common garden experiments and may be shaped, at least partly, by local adaptation. Habitat fragmentation may also have contributed to the observed population differentiation for bud burst in our study given the fact that Flanders, the northern part of Belgium, is already historically a very densely populated area. Habitat fragmentation affects genetic structures of tree populations including population divergence, as was, for instance, shown for Fagus sylvatica L. [33]. A stronger genetic control for bud burst in comparison with autumnal phenophases may be related to the fact that springtime leafing out modulates the timing of carbon assimilation in deciduous forests more than the autumnal phenophases. Early bud burst was associated with a longer period of carbon uptake and with an increase in net annual carbon fluxes [34]. The gain of 1 day extra photosynthesis in spring is more beneficial than an extra day in autumn when days are shorter and trees capture less solar radiation.
In trees, timing of bud burst can be serially auto-correlated with the timing of flowering [14,20], and the control of flowering time has been found to share certain regulations at the molecular level with control of the dormancy/growth cycle [35]. Still, our results indicate some level of deviating control between the timing of the two spring phenophases. Although population differentiation is relatively less high in timing of flowering compared with timing of bud burst, the intergenotypic variance within populations is relatively larger. Together with the larger time span between first and last flowering population compared with bud burst, this may point to a stronger influence of assortative mating. Assortative mating occurs when early flowering plants tend to receive pollen from other early flowering plants, and late flowering plants tend to pollinate later flowering plants, which can be suggested to create larger intergenotypic differences within populations. In a simulation study, already large levels of within-population genetic variation in bud burst were predicted [14], which can be hypothesised to be more accentuated in the timing of flowering, as this phenophase is more directly affected by assortative mating than bud burst. In addition, more diversified flowering phenologies in a population (larger intergenotypic variance) supports the hypothesis that phenological spreading in flowering time reduces the competition for generalist pollinators [36]. Likewise, for herbaceous plants, members of a population rarely flower synchronous [19].
In our study, the timing of flower opening displayed relatively less intragenotypic variance (between ramets) than leaf bud burst, suggesting a lesser sensitivity to environmental variation than leaf bud burst. Two main reasons can be suggested. Firstly, as woody plants only start flowering and fructifying several years after establishment, reproductive phenology only comes under selection after several years, at a time when the most vulnerable seedling phase has passed [37]. In addition, the consequences of a poorly timed flowering, which may result in meagre to absent reproduction, may be less detrimental than poorly timed flushing, which may lead to growth retardation, competitive disadvantages and mortality due to frost exposure [37]. As the penalty for a suboptimal flowering may be less detrimental than for a suboptimal bud burst, selection on timing of flowering may be less stringent than on flushing, even when both occur simultaneously. This may also lead to the relatively lesser population differentiation in timing of flowering observed in this study (bud burst displaying a relatively larger interpopulation variance). Secondly, the putative stronger penalty for a suboptimal bud burst, especially in the vulnerable seedling or sapling stage of a tree or shrub, may have evolved a certain level of plasticity on the individual plant level, which is needed to fine-tune the timing of bud burst to (micro-)environmental conditions. This may have caused the less strong genetic control for timing of bud burst in comparison with flowering in our results (a smaller relative intragenotypic variance-between ramets-for timing of flowering). As plasticity comes at a particular cost [3], the stronger genetic control in timing of flowering may in this sense indicate a lesser need for fine tuning of this timing to (micro-)environmental conditions compared with bud burst. The here described deviations between timing of bud burst and flower opening may be related to the phenomenon described for size variables in plants, with vegetative traits being more sensitive to environmental variation than floral traits [38,39].

Timing of Leaf Senescence and Leaf Fall
In the variance analysis, both leaf senescence and fall express relatively less population differentiation than the spring phenophases. As already mentioned, springtime leafing out modulates the timing of carbon assimilation in deciduous forests more than the autumnal phenophases [34]. The correlation between timing of bud burst and flower opening (r = 0.89) is higher than between timing of leaf senescence and fall (r = 0.72), although the first is a correlation between the timing of phenophases that take place in different organs (vegetative buds versus generative buds) whereas the latter are expressed consecutively in the same organ (vegetative leaves). This suggests a lesser genetic control for the autumnal phenophases. Both autumnal phenophases display a larger relative intragenotypic variance (between ramets), being relatively more pronounced in leaf fall. Although a tight control between the two phenophases is suggested [26], the latter observation suggests a partly deviating control. It can be hypothesised that leaf fall is comparatively more sensitive to other abiotic cues than photoperiod, such as temperature (frost) and wind, leading to slightly more variability at the individual plant level (higher relative intragenotypic variance compared with leaf senescence) as a response to deviating micro-climatic conditions. Nutrient resorption during leaf senescence is not total [40]. Falling leaves retain certain concentrations of nutrients. Our data suggest that a certain level of plasticity at the individual plant level, allowing leaves to retain on the shrub for a longer time in more favourable autumnal conditions, may allow further nutrient reallocation, resulting in a higher nutrient resorption efficiency (the proportion of the nutrients in green leaves prior to senescence that are resorbed before abscission [40]).

Variation between the Two Common Gardens
The negative correlation between timing of bud burst and leaf senescence means that genotypes with an advanced flushing have higher chances on a delayed leaf shedding. These timings thus tend to enlarge or diminish the length of the growing season. A negative correlation between bud burst and leaf colouring has been reported in half sib progenies of P. padus in Sweden [41] and in other species [42]. Conversely, in our experiment, all phenophases are delayed in Dentergem compared with Semmerzake and may, at least partly, be attributed to a phenotypic plastic response to the slightly deviating local growth conditions at the two plantation sites. Probably, the observed phenotypic plasticity, which is in itself also genetically controlled but acts on a shorter time scale in comparison to natural selection, delays or advances all phenophases in response to the local growth conditions in P. padus. The timing of one phenophase influences to a certain extent the timing of the following phenophases in the sense that advancement in spring can create advancement in autumn and so on [22]. It can be hypothesised that natural selection, aiming at enlarging the growing season while minimizing frost risk, can surpass the shorter term legacy effects of phenotypic plastic responses.

Conclusions
Studies on genetic variation in common gardens of local populations of woody plants are rare. Still, they may help predicting the way populations will react on a changing climate. We investigated the pattern of genetic control and the putative presence of spatial phenotypic plasticity in the seasonal phenophases of the shrub species P. padus planted in two common gardens consisting of local populations. Clearly, the contribution of genetic and non-genetic components to the total variability varied between the four studied seasonal phenophases. We found indications for natural selection in the four phenophases (differentiation at the level of the populations and expansion or contraction of the growing season's length at the level of genotypes). We also found indications for phenotypic plasticity (advancement or delay of all seasonal phenophases between the two plantation sites). These results, based on the variability of only local populations, corroborate the statement that trees harbor capacities for both types of responses to a changing environment.
In addition, we discovered a relatively large intergenotypic variability (within populations) in the timing of flowering, possibly suggesting that not only natural selection but also assortative mating may play a role. Because of the relatively lower intragenotypic variability in flowering time (in comparison with the other phenophases), a lower sensitivity to environmental variation can be hypothesized, so that disturbances may be less detrimental compared with the vegetative seasonal phenologies.
Supplementary Materials: The following are available online at http://www.mdpi.com/1999-4907/11/10/1070/s1, Figure S1: The difference between the mean temperatures in the closest weather stations to the two planting sites of the common gardens in Dentergem and Semmerzake for the period 1980 till 2014. Figure S2: Modelled timing for bud burst (a), flower opening (b), leaf senescence (c) and leaf fall (d) in the common gardens of P. padus for the different sampled local populations. Population abbreviations are in Table 1. Figure S3: Order of modelled timing of bud burst, flower opening, leaf senescence and leaf fall for the different sampled populations in the plantation of Semmerzake. Standardized days were calculated for which half of the plants of a population had