Trade-off between height growth and spring flushing in common beech (Fagus sylvatica L.)

Abstract

• Introduction

The annual development cycle of boreal and temperate trees results from an evolutionary trade-off between two opposing forces. These are namely, the adjustment of leaf phenology to the timing of frost occurrence at the beginning and/or the end of the growth season countered by an effective adjustment to the duration of the growth season to maximise photosynthesis and biomass production during the growing season.

• Objectives

In a provenance trial of common beech (Fagus sylvatica L.) populations, we tested whether avoiding frost events by late flushing is the primary process of frost tolerance, and also, whether longer growth seasons resulted in an increase in biomass production.

• Results

We detected clines with altitude and longitude of provenances for phenological traits, in contrast to growth variables. A highly significant negative genetic correlation was detected between frost damage and budburst date which confirms that the timing of leaf flushing is the principal determinant of tolerance to late frost.

Height and diameter growth were positively correlated with the length of the vegetation season. However, this effect was associated with the timing of growth cessation rather than with budburst.

• Conclusions

Stabilizing selection, selective pressure on early leaf flushing of juveniles in dense natural beech stands, and epigenetic regulation of budburst are proposed as explanations for these observations.

1 Introduction

Vegetative phenology is considered to be one of the determinants of growth and biomass production in trees (Mimura and Aitken 2010). The timings of budburst and leaf senescence limit canopy duration and thus affect the carbon assimilation rates and the energy budget of forest ecosystems at all levels from a single stand to the global scale (Kindermann et al. 1996).

There is wide agreement that the timing of leaf flushing in the spring results from opposing evolutionary pressures in forest trees. An early budburst increases the length of the available growing season (capacity adaptation). On the other hand, later budburst decreases the risk of frost damage to actively growing parts (survival adaptation; Leinonen and Hänninen 2002; Loehle 1998; Vitasse et al. 2009). An optimised trade-off between prolonged vegetation and frost damage avoidance is expected to ensure maximal fitness of an individual under a given climate in the long-term, and may thus be considered an evolutionarily stable strategy sensu Maynard-Smith (1982).

There are controversial reports about the mechanisms controlling budburst in European beech. In temperate deciduous trees, the developmental processes associated with bud dormancy release are generally controlled by temperature in the form of a required chilling period followed by the accumulation of forcing temperatures (Kramer 1994). The role of photoperiod is unclear, particularly in beech. Heide (1993) or Lechowicz (1984) emphasized the importance of day length. In contrast, Kramer (1994) concluded that bud dormancy in beech is determined mainly by winter chilling. There might also be a geographic variation in the signals for breaking bud dormancy in beech; Falusi and Calamassi (1996) reported that a possible compensation between winter chilling and day length may depend on the latitude of origin.

Less is known about the climatic signals controlling the onset of dormancy in the autumn, showing overall less year-to-year variability. Both the photoperiod and the temperature regime seem to play a role in the timing of leaf senescence (Delpierre et al. 2009; Estrella and Menzel 2006).

In spite of the general consent about the evolutionary trade-off between growth and frost damage avoidance reflected in the optimization of the timing of growing season duration, this hypothesis apparently has not been tested explicitly. To address it, we studied relationships between vegetative phenology, height and diameter growth, and frost damage in a provenance experiment with European beech (Fagus sylvatica L.). It has been demonstrated that at least the timing of budburst is a highly heritable trait in beech (von Wuehlisch et al. 1995a). Clear geographic clines and a strong correlation of budburst date with the temperature regime at the site of origin (von Wuehlisch et al. 1995b) indicate an adaptive significance of this trait. Under the trade-off hypothesis, a positive genetic correlation is expected between growth and growing season length, whereas budburst date and frost damage are expected to be negatively correlated. The objective of this study is testing, whether these expectations conform with reality in common beech.

2 Materials and methods

2.1 Experimental material

The study is based on a provenance experiment with European beech established under the coordination of the Federal Forest Research Centre, Institute of Forest Genetics in Grosshansdorf, Germany, comprising 23 trials. The majority of trial plots were planted in 2000 with 2-year-old seedlings. The Slovak trial, located at Tále near the town of Zvolen, central Slovakia (48°38′ N, 19°02′ E, 810 m a.s.l.), comprises 32 provenances covering most of the distribution range of beech in Europe, with the exception of the south western Balkans, south Italy and Spain (Table 1). The test was established under a randomized complete block design with three blocks, each plot containing initially 50 plants planted at 2 × 1-m spacing.

Table 1 List of provenances with geographical coordinates and provenance means (BLUPs) of growth, frost damage and phenological traits

Phenotypic traits were assessed in 2007 and 2008. Spring leaf flushing was scored at 3–5-day intervals between April 1 and June 1, 2007 (13 scoring dates) using a seven-point scale (Table 2). Autumn leaf colouring was scored at 6–10-day intervals between September 5 and November 1, 2007 (six scoring dates) using a six-point scale (Table 2). Scoring was repeated in 2008 between April 10 and May 27 (eight dates) and between August 29 and October 27 (ten dates) using the same methodology.

Table 2 Scoring scales for phenological stages and frost damage

Both budburst and leaf senescence represent processes of irreversible transition between two steady states. Therefore, phenological scores were fitted (least squares) to the sigmoid function to estimate the midpoint and duration of these processes:

$$ p = \frac{{1 + 2\tanh \frac{{t - m}}{d}}}{2} $$

where p is the phenological stage, t is the scoring date (Julian day), m is the midpoint (Julian day) and d is the duration (days) of leaf flushing and/or leaf colouring.

In 2007, two late frost events occurred during the flushing season, leading to extensive damage of beech forests in central Slovakia (April 22/23, May 1/2). Frost damage was scored in the trial on May 3, 2007, using a five-point scale (Table 2).

Height, root-collar diameter (at the 20 cm height) and breast-height diameter (at the 1.3 m height) were measured in 2008 before the beginning of shoot growth.

2.2 Data analysis

Phenotypic data were analysed under the mixed linear model using the procedure MIXED of SAS (2009). Blocks were considered a fixed-effect factor and provenances a random-effect factor. The best linear unbiased predictors (BLUP) of provenance means and variance components associated with random effects were also estimated by the procedure MIXED (Table 1).

To assess spatial independence of data, simple Mantel tests (correlations between the matrix of geographic distances and the matrix of differences in provenance means of phenotypic traits) were calculated and their significances were tested using 100,000 random permutations.

Existence of clinal variation was tested using linear regression of phenotypic data against geographical coordinates. Significances of regression coefficients were tested by the t test against H ₀: b = 0.

Genotypic and phenotypic correlations between phenotypic traits were estimated using the method of moments based on a multivariate extension of the ordinary least squares (multivariate analysis of variance, cf. Mode and Robinson 1959). The GLM procedure (SAS 2009) using the MANOVA option was employed for the analysis. Provenance components of variances and covariances were used to calculate genetic correlations, whereas the sums of all random components were employed for phenotypic correlations.

3 Results

Differences among provenances were significant for almost all traits except vegetation-season length in 2008 (Table 3). There were also differences among blocks in phenology and frost damage, and provenances reacted differently to this micro-environmental heterogeneity, as indicated by significant provenance-by-block interactions. This methodological failure is associated with the layout of the trial plot. It was established within a limited area of a former forest nursery and was surrounded by an adult beech stand, where any spatial arrangement of blocks would produce differences in the behaviour of the planted material due to differences in lateral shading. Nevertheless, similar situations are quite common in breeding and transplant experiments with forest trees, and we consider the existence of unambiguous trends in phenology and frost damage to be proof of the validity of our observations.

Table 3 Mixed-model analysis (F tests for the fixed effects, Wald's Z-statistics for the random effects) of growth, frost damage and phenological traits

Percentage of variance attributable to the effect of provenance varies; it is particularly high in the case of budburst date (34–38%), but also reaches almost 20% for height growth (Table 3). Despite this, the residual within-provenance variation, partly caused by genetic variation among trees within populations, represents by far the largest component in all traits, always exceeding 55%.

Growth traits (height, diameters) did not exhibit any geographical patterns: they were neither spatially autocorrelated (Table 4), nor did they show any range-wide clinal trend (Fig. 1, Table 5). However, the results confirmed the existence of clear geographic clines in the timing of both spring flushing and autumn leaf colouring. Provenances from the eastern part of the range and from lower elevations tend to flush earlier in spring and terminate photosynthesis earlier in autumn (Fig. 1, Table 5). On the other hand, no latitudinal cline was found (see a complete overview of phenology—geography relationships provided in the Eletronic Supplementary Material). There was no clear pattern in the growing season length (Fig. 1).

Table 4 Simple Mantel tests of spatial autocorrelation of growth, frost damage and phenological traits

Fig. 1

Relationships between selected phenotypic traits and geographical coordinates. Superscripted number one:average over 2007 and 2008

Table 5 Linear regressions (R ² and regression coefficients) of growth, frost damage and phenological traits against geographical coordinates

Budburst date is very stable across years; genetic as well as phenotypic correlations between 2007 and 2008 are positive and highly significant for midpoint, but much less for the duration of leaf flushing. Leaf senescence midpoint is less correlated between years, which may indicate a more complicated regulation of this process. On the other hand, the duration of leaf colouring is more stable, compared to budburst. The vegetation period duration as the difference between leaf flushing and leaf discoloration dates shows an intermediate interannual stability (Table 6).

Table 6 Genetic and phenotypic correlations of phenological traits between years (2007 and 2008)

There is a very close relationship between budburst date and duration (R ² = 0.67); later-flushing provenances tend to flush rapidly. Interestingly, just the opposite was observed for leaf discolouration (Fig. 2).

Fig. 2

Relationships between budburst date and duration (left), and leaf discolouring date duration (right). Superscripted number one: average over 2007 and 2008

Frost damage was closely correlated with the leaf-flushing date: both genetic and phenotypic correlations are negative and highly significant (Table 7; since none of the variables are spatially independent, significance tests must be considered approximate).

Table 7 Genetic and phenotypic correlations between growth/frost damage and phenological traits

A temporal shift of budburst by 1 day decreases late frost damage by 2% (Fig. 3).

Fig. 3

Relationships between vegetation season length and height growth (left), and budburst date and late frost damage (right). Superscripted number one: average over 2007 and 2008

Genetic correlations between growth traits (height, diameter) and the vegetation period duration are positive, fairly high and highly significant in both years (Table 7; see also the Electronic supplementary material). This relationship is illustrated by the example of height growth: a prolongation of the growing season by 1 day improves height at the age of 10 years by 7.4 cm (Fig. 3). At the same time, the correlations of growth variables with growing season length are higher than those with spring flushing or autumn leaf colouring dates. The only exception is the correlation with breast-height diameter, where correlations with the timing of leaf senescence were found to be the greatest. However, this trait was measured only on a subset of the highest individuals, exceeding the recording height of 1.3 m.

4 Discussion

Strictly speaking, the definition of genetic correlation implies using clonal replicates of single genotypes. However, in evolutionary studies as well as in the practice of tree breeding, a wider genetic or taxonomic entity such as family, population or taxon is often widely suggested to constitute a genotype (Schlichting 1986; Chambel et al. 2005). We used the terms “genetic” and “phenotypic” correlations as a sort of simplification rather than the more cumbersome terms “provenance-level” and “individual-level” correlations, respectively. In fact, tree-level correlations also include a genetic component, associated with the genetic variation among individuals within a population.

The phenotypic variation in trees cannot be separated from its geographical context. Beech phenology is known to vary along both longitudinal and altitudinal gradients (von Wuehlisch et al. 1995b), which was also confirmed by our observations. Climatic contrasts among beech populations are huge, as illustrated by Fig. 4. At an oceanic site (Aguilar, La Rioja, Spain), the accumulation of temperatures is much more rapid than that in the interior (Kráľová, near the trial site). An arbitrarily chosen temperature sum of 500°C was reached on March 5 in Aguilar, but this was not achieved until as late as April 27 in Kráľová, despite an almost similar altitude for both sites. Using constant sum of temperatures as a signal for vegetative processes would cause either an ineffective use of the vegetation season under the continental climate or a shift of flushing into wintertime under the oceanic climate. Populations of perennial plants have to develop mechanisms to cope with this climate variation.

Fig. 4

Course of average daily temperatures and cumulative degree days during the winter/spring 2007 at two contrasting sites. Sources of data: Technical University in Zvolen and Dirección General de Justicia e Interior de La Rioja (http://ias1.larioja.org/estaciones/estaciones/mapa/portada/index.jsp)

As mentioned in the introduction, the adaptation of the timing of vegetative processes of trees to climate (at least in the boreal and temperate zones) has traditionally been considered to result from stabilizing selection driven by two opposing forces: the avoidance of frost damage during the cold season of the year, as opposed to maximum use of favourable growing conditions for photosynthesis and growth (Leinonen and Hänninen 2002). Such selection mechanisms require that avoiding frost events by late flushing is the primary mechanism of freezing tolerance, and that longer time available for growth is reflected in a higher biomass production.

Plants are known to protect themselves against cold by the accumulation of various compounds in the cytoplast and cell structures, such as oligosaccharides (Stushnoff et al. 1997), free proline (Reyes-Diaz et al. 2005), dehydrins (Kosová et al. 2007), heat-shock proteins (Lopez-Matas et al. 2004), or by changes in lipid composition and the degree of fatty acid saturation of the cell membranes (Čiamporová and Trgiňová 1999). These substances are believed to be involved in the freezing point depression of cell sap, the prevention of plasmolysis during cell dehydration caused by freezing, and also in membrane stabilization (Reyes-Diaz et al. 2005). However, cold acclimation requires a certain time, as the maximum freezing tolerance of plants is not a constitutive trait, but is induced in response to low but non-freezing temperatures (below 10°C; Alberdi and Corcuera 1991). Therefore, such defence mechanisms apply to mature vegetative tissues (such as evergreen leaves or needles) rather than to flushing leaves in the spring. In our case, there was a very close negative provenance-level correlation between frost damage and the leaf-flushing stage, implying that the primary risk-avoidance strategy of beech with respect to late frosts is to delay bud break. Naturally, trees are capable of re-flushing following frost defoliation. In fact, beech woods in central Slovakia (including the damaged part of the provenance plot) recovered from the damage in May 2007. However, the second flush leaves differ in morphology and have lower nutrient and sucrose concentrations, which indicates a depletion of nutrient and carbohydrate leaf resources and may place constraints on growth capacity (St. Clair et al. 2009). Moreover, flowering in beech occurs simultaneously with leaf flushing. Systematically repeated destruction of the flowers by late frosts may hinder the natural regeneration of beech forests and thus endanger the persistence of beech populations over several generations.

Prolongation of the growing season was reflected in a higher growth rate. However, this effect seems to be associated primarily with delayed cessation rather than with an early beginning of vegetation. Genetic correlations of height or diameter with budburst date are consistently negative but non-significant. This is in agreement with the observations in other tree species, where such correlations are generally weak (El-Kassaby and Park 1993; Rehfeldt 1992), and in some cases even positive, meaning that trees may benefit from delayed bud break (Hawkins and Stoehr 2009; Li and Adams 1993). On the other hand, growing season length or vegetation cessation date are generally more closely correlated with height, diameter, dry weight or other indicators of biomass production (Dewald et al. 1992; Hawkins and Stoehr 2009; Rehfeldt 1992).

In temperate deciduous forests, early flushing of seedlings and saplings prior to the canopy closure is critical for their survival and growth (Augspurger 2008), thus an advance in bud break decreases with increasing height and age (Seiwa 1999). The advantage of early-flushing juveniles in dense beech woods may therefore be another mechanism maintaining the average timing of budburst at earlier dates.

The existence of clear clines has traditionally been interpreted as an outcome of adaptation by natural selection. There is, however, growing evidence of epigenetic memory mechanisms regulating the phenological behaviour of tree populations (Johnsen et al. 2005; Rohde and Junttila 2008). Most of this evidence is associated with conifers (namely Norway spruce), having an incomparably larger genome and differing in epigenetic markings associated with heterochromatin (Fuchs et al. 2008). Moreover, beech as a long-lived species and difficult to be artificially crossed or propagated by somatic embryogenesis is among the least appropriate objects for the study of memory effects in angiosperms. Nevertheless, phenological behaviour of beech is quite fine-tuned to the temperature regime of the sites of origin (von Wuehlisch et al. 1995b), and it is questionable as to whether this can be explained by natural selection alone.

The absence of latitudinal trends and the presence of clear longitudinal and elevational clines support the hypothesis of climatic (temperature) rather than photoperiodic signals for the course of phenological processes. However, the genetic backgrounds of spring and autumn vegetative phenology must differ. The beginning of budburst is obviously heritable, but its course is driven by actual temperatures: provenances starting early accomplish budburst under generally colder weather and thus need a longer time to finish flushing. However, the relationship between mid-date and duration is just the opposite in autumn discoloration: western, late provenances require more time compared to local provenances to finish leaf colouring, in spite of decreasing temperatures during the late autumn and the risk of early frost damage. Autumn phenology is thus less affected by the environment of the plantation site; the behaviour of provenances apparently results from accommodation to the climate of the site of origin.

In conclusion, the data on common beech confirmed that an evolutionary trade-off between growth and frost avoidance exists to a certain extent, and is associated mainly with the autumn vegetative phenology. However, the underlying genetic (and potentially epigenetic) mechanisms are largely unknown and deserve further study.