Original Research ARTICLE
Intraspecific Variation in Wood Anatomical, Hydraulic, and Foliar Traits in Ten European Beech Provenances Differing in Growth Yield
- 1Plant Ecology, Albrecht von Haller Institute for Plant Sciences, University of Göttingen, Göttingen, Germany
- 2Faculty of Forestry, Technical University in Zvolen, Zvolen, Slovakia
- 3Federal Research Institute for Rural Areas, Forestry and Fisheries, Thuenen Institute for Forest Genetics, Großhansdorf, Germany
- 4UMR BIOGECO Institut National de la Recherche Agronomique-UB, University of Bordeaux, Talence, France
In angiosperms, many studies have described the inter-specific variability of hydraulic-related traits and little is known at the intra-specific level. This information is however mandatory to assess the adaptive capacities of tree populations in the context of increasing drought frequency and severity. Ten 20-year old European beech (Fagus sylvatica L.) provenances representing the entire distribution range throughout Europe and differing significantly in aboveground biomass increment (ABI) by a factor of up to four were investigated for branch wood anatomical, hydraulic, and foliar traits in a provenance trial located in Northern Europe. We quantified to which extend xylem hydraulic and leaf traits are under genetic control and tested whether the xylem hydraulic properties (hydraulic efficiency and safety) trades off with yield and wood anatomical and leaf traits. Our results showed that only three out of 22 investigated ecophysiological traits showed significant genetic differentiations between provenances, namely vessel density (VD), the xylem pressure causing 88% loss of hydraulic conductance and mean leaf size. Depending of the ecophysiological traits measured, genetic differentiation between populations explained 0–14% of total phenotypic variation, while intra-population variability was higher than inter-population variability. Most wood anatomical traits and some foliar traits were additionally related to the climate of provenance origin. The lumen to sapwood area ratio, vessel diameter, theoretical specific conductivity and theoretical leaf-specific conductivity as well as the C:N-ratio increased with climatic aridity at the place of origin while the carbon isotope signature (δ13C) decreased. Contrary to our assumption, none of the wood anatomical traits were related to embolism resistance but were strong determinants of hydraulic efficiency. Although ABI was associated with both VD and δ13C, both hydraulic efficiency and embolism resistance were unrelated, disproving the assumed trade-off between hydraulic efficiency and safety. European beech seems to compensate increasing water stress with growing size mainly by adjusting vessel number and not vessel diameter. In conclusion, European beech has a high potential capacity to cope with climate change due to the high degree of intra-population genetic variability.
European beech (Fagus sylvatica L.) dominated the natural vegetation types of forests in Central Europe for centuries, forming large stands of resilient forest ecosystems (Ellenberg and Leuschner, 2010). Despite the competitive superiority of European beech for tree populations in temperate forests, this species is more vulnerable to drought-induced stem growth reductions than other temperate broad-leaved trees (Leuschner et al., 2001; Zimmermann et al., 2015). Consequently, the Central European beech populations are expected to suffer high mortality rates probably altering the distribution range at its dry distributional limit as a consequence of increased physiological stress due to a higher risk of drought exposure associated with recent climate change.
According to its large geographic distributional range, European beech is expected to exhibit substantial genetic diversity (Bolte et al., 2007). The inherent capability of this species to survive and reproduce successfully across such a wide range of habitats is either maintained by long-term adaptation (i.e., genotypic variation) or short-term acclimation (i.e., phenotypic plasticity; Lindner et al., 2010; Kremer et al., 2014). Consequently, evolutionary adaptation may have caused the selection of ecotypes adapted to the regional climatic conditions, which is manifested in phenotypic variation of plant functional traits across various provenances (Hamrick, 2004). Despite an increasing amount of studies on functional traits which are known to be related to drought resistance in beech, e.g., leaf stomatal features (Stojnić et al., 2015), wood structure (Eilmann et al., 2014), resistance to xylem cavitation (Herbette et al., 2010, Wortemann et al., 2011), or hydraulic architecture (Aranda et al., 2015; Schuldt et al., 2016), still very little information is available on the relevance of phenotypic plasticity or genetic variation of traits on the adaptation of populations.
The high trait plasticity and remarkable regeneration potential of beech trees after disturbance leads to the assumption that beech is able to bear water shortage to a certain degree (Kahle, 2006; van der Werf et al., 2007). However, other studies report on drought susceptibility or dramatic yield loss and regeneration periods of decades subsequent to drought events and there is still an ongoing debate on the response of beech populations to drought exposure (Leuschner et al., 2001; Peuke et al., 2002; Leuzinger et al., 2005; Bréda et al., 2006; Jump et al., 2006). The mechanisms behind the adaptation and adaptability of beech to drought is a main challenge for a better understanding of the effects of climate change on this economically and ecologically important tree species. Because productivity is closely coupled with hydraulic efficiency (Hajek et al., 2014; Hoeber et al., 2014; Kotowska et al., 2015) and may only be achieved at the cost of hydraulic safety (Cochard et al., 2007), genetic differentiation in wood anatomical and hydraulic traits should be most pronounced between provenances differing in yield. However, the way how different provenances of the same tree species cope with such trade-offs still remains poorly investigated.
The main objective of this study was to investigate the genetic differentiation in vulnerability to xylem embolism and other related hydraulic and foliar properties across ten European beech provenances differing in growth performance. We further aimed to investigate the relationships and potential trade-offs between hydraulic, wood anatomical and leaf traits within species. We hypothesized that (i) drought-related traits are under genetic control and therefore populations are locally adapted to their environment, (ii) hydraulic efficiency and embolism resistance are related to anatomical traits, and (iii) a high hydraulic efficiency leads to high growth rates at the cost of xylem safety.
Materials and Methods
Experimental Site, Plant Material, and Microclimatic Conditions at the Place of Origin
The field trial studied is part of the International Beech Provenance Trial Series 1993/1995 under EU funding (AIR3-CT94-2091), which investigates the role of the genetic variation of beech for adaptability, productivity and selected ecosystem functions considering risks of global climate change (von Wuehlisch et al., 1998). For the present study conducted in August 2014, the common-garden field trial with 100 different European beech (Fagus sylvatica L.) provenances established in 1995 in Northern Germany (Schleswig-Holstein) with 2-year old saplings near the Schädtbek Experimental Farm (54° 18′N, 10°16′E, 40 m a.s.l.) was used. The climate at the site is oceanic, moderately cold with a mean annual temperature (MAT) of 8.3°C, a mean annual precipitation (MAP) of 742 mm and a mean early growing season precipitation from April to June (MSP) of 149 mm (data obtained from the German Meteorological Service). The common garden trial consists of 10 × 10 m plots planted with 50 beech saplings per provenance in a 3-times replicated randomized block design. The beech trees are arranged in a rectangular grid with a planting distance of 2 m between and 1 m within rows. The trial is surrounded by a single bordering tree row serving as buffer zone to avoid edge effects. As all provenances were grown in a single environment (common garden), we were able to assess the genetic differentiations between provenances for several functional traits. The common garden experienced no management operations since plantation establishment in 1995, resulting in competitive selection among tree individuals leading to heterogeneity in plot characteristics.
From the 100 available provenances, ten provenances differing in aboveground growth increment (Figure S1) were selected in order to cover not only a gradient in growth yield but also a climatic gradient at the place of origin (Figure S2). By this selection, we covered provenances from different geographic regions and climates throughout Europe with a broad range of MAT (3.4–15.3°C) and MAP (575–1080 mm; Table 1). Mean annual climate data from 1950 to 2000 at the place of origin for each of the ten provenances were obtained from the WorldClim database with 30 arc-seconds resolution (Hijmans et al., 2005). In order to estimate climate dryness at the place of origin, we used the WorldClim database to access the global aridity index (GAI, Zomer et al., 2008) and calculated the forest aridity index (FAI) according to Fuehrer et al. (2011) as FAI = 100 × TJul−Aug/(PMai−Jul + PJul−Aug), where T is the temperature and P the precipitation of the associated interval. Since the atmospheric evaporative demand in the growing season is highest in mid-summer (July and August), the July precipitation was weighted by a factor of two in the denominator. We further calculated Ellenberg's climate quotient (EQ, Ellenberg and Leuschner, 2010) as EQ = (TJul/Pannual) × 1000. All measures of water availability and climatic aridity at the places of provenance origin were highly interrelated (data not shown). We therefore decided to use the forest aridity index (FAI) as integrate variable for all subsequent analyses.
Table 1. Climatic data at the place of origin of the ten Fagus sylvatica provenances including elevation, mean annual temperature (MAT), mean annual precipitation (MAP), mean early growing season precipitation from April to June (MSP), Ellenberg's climate quotient (EQ), the global aridity index (GAI) and the forest aridity index (FAI).
Provenances were selected according to differences in aboveground growth increment based on a tree inventory in January 2013 (Figure S1). Within a given provenance, ten tree individuals of comparable size close to the population average (diameter at breast height and tree height) were selected assuming uniform growth conditions (site conditions and intraspecific competition) and consequently representing the average growth increment of a given provenance, yielding 100 processed tree individuals in total. A 2- to 4-year-old branch segment (mean ± SE : 2.25 ± 0.05 year) was collected in August 2014 from the uppermost canopy of each selected tree individual with a long-reaching telescope pruner and recut to approximately 50 cm length on the ground. Selected segments were defoliated and immediately transferred to plastic tubes containing deionized water and Micropur (Katadyn, Wallisellen, Switzerland) to prevent microbial activity and stored at 4°C until further processing within 4 weeks. All leaves of the respective segments were stored separately for foliar analyses. A list of all measured traits, their symbols and units are given in Table 2.
Timber Volume, Above-Ground Biomass, and Basal Area Increment
Aboveground growth performance (i.e., stem increment and height growth) of the respective genotypes under the local environment was evaluated from diameter at breast height (DBH, cm) and tree height (m). Aboveground productivity expressed as aboveground biomass increment (ABI, kg yr−1) was calculated for the entire growth period (1995–2014), and basal area increment (BAI, cm2 yr−1) from two inventories in January 2013 and August 2014, respectively. For calculating the volume of stems and branches of at least 7 cm in diameter (standing volume of timber), we used the allometric equation V = π × ((D/100)/2)2 × H × f, where V is the volume of timber (m3), D the diameter at breast height (cm), H tree height (m), and f an empirically derived form factor for beech trees (Bergel, 1973) with f = 0.4039 + (0.0017335 × H) + (1.1267/H)–(118.188/D3) + (4.2 × 106 × D2). Above-ground biomass was estimated from an empirical equation given by Wutzler et al. (2008) as AGB = 0.00523 × D2.12 × H0.655.
From each branch segment, all leaves were removed from the basipetal segment upwards to determine mean leaf size (Aleaf, cm2) and cumulative leaf area (AL, m2) using a flatbed scanner and the WinFOLIA software (Régent Instruments, Quebec, Canada). Per branch segment, 19 to 194 leaves were scanned, yielding 6008 leaves in total. Specific leaf area (SLA, cm2 g−1) was calculated by dividing the total leaf area by the leaf dry weight (70°C, 48 h). The Huber value, i.e., sapwood to leaf area ratio (AS : AL, 10−4 m2 m−2), was calculated by dividing maximal sapwood area by AL. Subsequently, leaf samples were ground and the leaf dry mass analyzed for foliar C and N concentrations as well as the carbon isotope signature (δ13C, ‰) by using a Delta Plus Isotope mass ratio spectrometer (Finnigan MAT, Bremen, Germany), a Conflo III interface (Thermo Electron Corporation, Bremen, Germany) and a NA2500 elemental analyser (CE-Instruments, Rodano, Milano, Italy) using standard δ notion: δ = (Rsample/Rstandard–1) × 1000 (‰) at the Centre for Stable Isotope Research and Analysis (KOSI), University of Göttingen. The foliar concentrations of Ca, K, Mg, and P were measured by ICP analysis (Optima 5300 DV, PerkinElmer Inc., USA).
Branch Xylem Anatomy, Theoretical Conductivity, and Growth Rate
Measurements of anatomical parameters were carried out on 9–10 branch segments from the basipetal end of the samples used for the hydraulic measurements of each provenance, yielding 98 samples in total. Prior to cutting semi-thin (10–20 μm) transverse sections with a sliding microtome (G.S.L.1, Schenkung Dapples, Zürich, Switzerland), the ethanol-stored (70%) segments were completely stained with a safranin solution (1 in 50% ethanol, Merck, Darmstadt, Germany) and subsequently embedded in Euparal medium. For digitalization, a stereo-microscope equipped with an automatic stage was used (SteREOV20, Carl Zeiss MicroImaging GmbH, Jena, Germany; Software: AxioVision v4.8.2, Carl Zeiss MicroImaging GmbH, Jena, Germany), enabling a time-efficient digitalization of the complete cross-section at 100x magnification. Image analysis was performed using the software Adobe Photoshop CS2 (Version 9.0, Adobe Systems Incorporated, USA) and ImageJ (v1.44p, http://rsb.info.nih.gov/ij) applying the particle analysis function. For all subsequent calculations, the complete xylem cross-section without pith and bark was analyzed, yielding 888,358 analyzed vessels in total. The following parameters were calculated: idealized vessel diameter (D, μm) as obtained from major (a) and minor (b) vessel radii according to Lewis and Boose (1995) as D = ((32 × (a × b)3)/(a2 + b2))¼, vessel density (VD, n mm−2) and cumulative vessel lumen area. The lumen to sapwood area ratio (Alumen : Axylem, %), was obtained by dividing cumulative vessel lumen area by the corresponding sapwood area. The diameter of individual vessels was used to calculate hydraulically-weighted vessel diameter (Dh, μm) according to Sperry et al. (1994) as Dh = ΣD5/ΣD4. Theoretical specific conductivity (, kg m−1 MPa−1 s−1) was calculated according to the Hagen-Poiseuille equation as = (((π × ΣD4)/128 η) × ρ)/Axylem, where η is the viscosity of water (1.002 10−9 MPa s), ρ the density of water (998.2 kg m−3), both at 20°C, and Axylem (m2) the corresponding xylem area without pith and bark. Branch growth rate (Agrowth, mm2 yr−1) was calculated by dividing Axylem by the number of growth rings (i.e., branch age, BA). In addition to , theoretical leaf-specific conductivity was calculated by division by AL as = (((π × ΣD4)/128 η) × ρ)/AL.
Hydraulic Conductivity Measurement
Hydraulic traits were measured in ten branch segments (mean diameter ± SE: 6.83 ± 0.06 mm) per provenance using the Xyl'em apparatus (Bronkhorst, Montigny-les-Cormeilles, France). In the laboratory, all lateral branches were cut off and the scares sealed with quick-drying superglue (Loctite 431, Henkel, Düsseldorf, Germany) applicable to wet surfaces. Subsequently, the segments were shortened to a length of 290.5 ± 0.8 mm (mean ± SE). For the determination of maximal hydraulic conductivity (Kh, kg m MPa−1 s−1) at 6 kPa, demineralized filtered (0.22 μm) and degassed water (10 mM KCl and 1 mM CaCO3) was used, interrupted by three 10-min flushes at 120 kPa to assure removal off all potential emboli. The diameter of each segment was measured twice at the basipetal and distal end, and at four positions along the segment. The following regression coefficients were used to calculate sapwood area without pith and bark for a given beech branch segment diameter according to Schuldt et al. (2016): Axylem = −3.715 + 0.770 Across. Subsequently, empirical specific conductivity (, kg m−1 MPa−1 s−1) was calculated by dividing Kh by the maximal basipetal, and not average, sapwood area (Hajek et al., 2014; Hoeber et al., 2014; Schuldt et al., 2016). Kh was further used to calculate empirical leaf-specific conductivity (, kg m−1 MPa−1 s−1) by division by AL.
Xylem Resistance to Cavitation
Vulnerability to xylem cavitation was determined on 8 to 10 branch samples (replicated trees) per provenance using the Cavitron technique (Cochard et al., 2005), yielding 94 samples in total. Segments with a standardized length of 27.5 cm were mounted in a custom-built honeycomb rotor chamber of the Cavitron, which uses a commercially available centrifuge as basis (Sorvall RC-5C, Thermo Fisher Scientific, Waltham, MA, USA), and spun at defined velocities recorded with the software CaviSoft (version 220.127.116.11, University of Bordeaux, France). Conductivity measurements started at 1.0 MPa and were stepwise repeated at intervals of 0.2 to 0.3 MPa until the percent loss of conductivity (PLC) reached at least 90%. All vulnerability curves of the present study measured on 27.5 cm long segments were s-shaped, indicating that no open vessels were present. For each branch segment, a sigmoid function (Willigen and Pammenter, 1998) was fitted to describe the relationship between PLC and xylem pressure using the expression PLC = 100/(1 + exp(s/25 × (Pi – P50))), where P50 (MPa) is the xylem tension causing 50% loss of hydraulic conductivity and s (% MPa−1) is the slope of the curve at the inflexion point. The xylem pressures causing 12% (P12) and 88% (P88) loss of conductivity were calculated as well.
In order to assess the significance of differentiation between provenances for all functional traits, we used linear mixed effect models (LME) starting with a random effect model without any fixed effect but with “provenance” and “block” added as a random effects using the “lme” function of the R package “nlme” according to the following model : Yijk = μ + Pi + bj + εijk, where Yijk is the observation of tree individual k for one of the analyzed characters from provenance i and block j, μ is the overall mean, and Pi the error term for provenance, bj for block and εijk for the residual variation. During the analysis, normal distribution of the residuals and homogeneity of variance were assessed visually using residual diagnostics and quantile-quantile plots; some data (, , AS : AL) had to be log-transformed in order to achieve normal distribution. In order to test whether the variance component accounting for genetic differentiation is significantly different from zero, a likelihood ratio test (LRT) was performed against a reduced model without a random effect for provenance using the restricted maximum likelihood (REML) method. As under those conditions the LRT is performed on the boundary of the parameter space, the resulting P-values had to be corrected by multiplying them by 0.5 (Verbeke and Molenberghs, 2009).
To quantify the influence of climate at the place of origin, FAI was added to the LME as fixed variable according to the following model : Yijk = α + β FAIi + Pi + bj + εijk, where α is the intercept and β the slope of the linear model with FAIi, which all three together determine μ of model . A significant influence of climate at provenance origin was determined by a LRT using the maximum likelihood (ML) method.
We further calculated coefficients of variation for within-provenance variation (CVintra) and between-provenance variation (CVinter) for each trait in order to allocate total measured trait variation to a genetic component (CVinter) and a predominantly phenotypic component (CVintra); the between-provenance variability (CVinter) was calculated from the between-provenance standard deviation (SD) and the overall mean value. Additionally, the ratio of provenance variance component to total variance was calculated using the R package “varComp” according to a variance component analyses with the program “lme” to calculate the proportion of total variance () explained by the variability between provenances (), replicated randomly distributed “blocks” () and residual variance within provenances (). Inter-population variance component (VCinter) was calculated according to VCinter = (/( + + )) × 100, variance component between replicated plots as VCblock = (/( + + )) × 100, and intra-population variance component (VCintra) as VCintra = (/( + + )) × 100, all in percentage.
Pearson correlation analysis was used to test for inter-relationships between different branch traits of the trees and for detecting relationships between traits based on data pooled across all provenances. Applying correlations to raw measurements for such exploratory analysis may mix several sources of variation and result in inflated degrees of freedom and potentially overestimation of significance. Statistical analyses were performed with the software R (version 3.1.3, R Development Core Team, 2011), and all linear and non-linear regression analyses were carried out with the software Xact 8.03 (SciLab, Hamburg, Germany).
Genetic Differentiation between Provenances
The estimate of aboveground growth potential as inferred from the basal area increment (BAI) from January 2013 to August 2014 and aboveground biomass increment (ABI) for the entire length of the experiment from August 1995 to August 2014 revealed large significant differences between the ten selected beech provenances (Table 3). The provenances with highest ABI from Slovakia and Ukraine (SK: 3.6; UA 4.0 kg yr−1, P < 0.001; Table S1) showed four times higher growth rates than the provenances with lowest growth rates from Slovenia and Sweden (SL: 1.02 and SE: 1.69 kg yr−1; Table S1). Both provenances had the highest radial growth rates (SK: 33.36 and UA: 23.20 cm2 yr−1) and grew more than seven-fold compared to the two provenances with lowest basal area increment (ES: 4.34 and SL: 5.68 cm2 yr−1; Table S1). The other six provenances reached intermediate growth rates (13.05–22.55 cm2 yr−1). ABI and BAI scaled highly positive with tree size (DBH, r2 = 0.98, P < 0.001 and r2 = 0.73, P < 0.001; height, r2 = 0.66, P < 0.01 and r2 = 0.31, P < 0.001; Table 5).
Table 3. Results of a random effects model on the genetic differentiation between provenances and the coefficient of variation for all traits measured for the variability between provenances (CVinter) and within provenances (CVintra).
Significant genetic differentiations between provenances were found in three out of 22 ecophysiological traits only, namely the xylem pressure inducing 88% loss of hydraulic conductance (P88), vessel density (VD) and mean leaf size (Aleaf; Table 3). The P88 values ranged from −3.60 MPa for the most vulnerable provenance from Germany (DE-BB) to −4.21 MPa for the most resistant provenance from Slovenia (SL), i.e., by 15% (Table S1). The differences in embolism resistance between provenances were, however, not mirrored in xylem anatomical adjustments. VD was significantly higher in the branch xylem of the most-productive provenances from Slovakia (SK) and the Ukraine (UA) while the largest leaves were found for the provenance adapted to the local climate of Schleswig Holstein (DE-SH; Table S1). Other provenances from Sweden (SE) and Czech Republic (CZ) produced on average 30% smaller leaves, which was not mirrored in adjustments of specific leaf area (SLA) or the Huber value (AS : AL). In general, the leaf chemistry was uniform across the provenances.
Significant genetic differentiation between provenances (Table 3) was mirrored in high values of inter-population variance component (VCinter) for these three ecophysiological traits (P88, VD, Aleaf) ranging from 8.38 to 14.02% (Table 3). Interestingly, all three measures of embolism resistance showed a higher variance between blocks (VCblock : ~25%) than between populations (VCinter : ~10%). The remaining measured hydraulic traits, empirical specific conductivity () and empirical leaf-specific conductivity (), likewise differed by ~10% between blocks but showed no variance between populations. Contrary to these hydraulic traits, all wood anatomical and calculated theoretical specific conductivity () revealed no variance between blocks. For all traits measured excluding the growth-related variables, variability within provenances (CVintra) was about two times higher than variability between provenances (CVinter; Table 3).
Trait-Relatedness to Climate at Origin
We observed several significant linear relationships between functional traits and the forest aridity index (FAI) as a measure of the climatic conditions at the place of provenance origin (Figure 1). These simple linear regression analyses were supported by linear mixed effect models confirming that FAI had a significant influence on seven of the 27 measured functional traits (Table 4). With increasing FAI, vessel diameter (D), lumen to sapwood area ratio (Alumen : Axylem), theoretical specific conductivity (), theoretical leaf-specific conductivity () and the C:N ratio significantly increased while the carbon isotope signature (δ13C) declined (Figures 1A–F; Table 4). Surprisingly, the three traits showing genetic differentiation between provenances (P88, VD and Aleaf) were not significantly affected by FAI. However, our results show that provenances originating from dry habitats with high FAI values (e.g., BG, ES, DE-BB) form particularly wide vessels and have high lumen to sapwood area ratios (Figures 1A,B) resulting in a high (Figure 1C) and (Figure 1D) compared to the other provenances when grown at this humid site with high precipitation rates received during the entire growing season. The C:N ratio was also significantly positively related to FAI, while a strong significant negative relationship between FAI and δ13C was observed (Figure 1F).
Figure 1. Relationship between the forest aridity index (FAI) at the place of origin and vessel diameter (D, A), lumen to sapwood area ratio (Alumen : Axylem, B), theoretical specific conductivity (, C), theoretical leaf-specific conductivity (, D), the carbon to nitrogen ratio (C:N, E) and the carbon isotope signature (δ13C, F). Given values are means per provenance. Each symbol represents one provenance, for provenance abbreviations see Table 1 (BG; CZ; DE-BB; DE-SH; ES; RO; SE; SK; SL; UA).
Table 4. Results of a linear mixed effect model examining the influence of the forest aridity index (FAI) at the place of origin on 27 measured parameters 19 years after planting of the provenance trial in Northern Germany.
Determinants of Aboveground Growth Performance
Among the xylem anatomical traits, only VD was positively related to aboveground biomass increment (ABI; r2 = 0.66, P < 0.005; Figure 2A) and negatively to branch growth rate (Agrowth; r2 = 0.35, P < 0.05; Figure 2B), which is a simplified measure of the annually produced branch sapwood, at the provenance level (n = 10). At the tree level (n = 100), only the negative relation between VD and Agrowth could be confirmed (r2 = 0.26, P < 0.001; Table 5). Moreover, VD strongly increased with increasing tree height (r2 = 0.83, P < 0.001; Figure 2C). The remaining branch wood anatomical and hydraulic properties, however, varied independently from aboveground growth performance, and no correlations were found between ABI, embolism resistance and hydraulic efficiency, respectively (Table 5). However, a weak though significant negative relation was found for basal area increment (BAI) and P50 at the tree level (r2 = 0.08, P < 0.01; Table 5). Although we only observed weak or no relations between ABI and the wood anatomical and hydraulic traits, several of them were related to Agrowth (Table 5). Hence, the expected trade-off between hydraulic efficiency and growth could be confirmed at branch level, but not at the tree level.
Figure 2. Relationship between vessel density (VD) and aboveground biomass increment (ABI, A), annual branch sapwood area increment (Agrowth, B), and tree height (C). Given values are means per provenance, for symbol definition see Figure 1.
Table 5. Pearson correlation coefficients for linear relationships between 27 functional trait variables at the tree level (n = 100).
Inter-Relationships between Functional Traits
We found that hydraulic efficiency in beech depends on the xylem properties. The vessel lumen to sapwood area ratio (Alumen : Axylem), annual branch growth rate (Agrowth), vessel diameter (D) and hydraulically-weighed vessel diameter (Dh) were closely positively related to empirical specific conductivity () and empirical leaf-specific conductivity (). In contrast, correlations between vessel density (VD) and the hydraulic traits or were not found (Table 5). Furthermore, none of the wood anatomical traits were related to embolism resistance and the expected trade-off between hydraulic conductivity and embolism resistance was absent. In addition, we found significant links between functional leaf traits such as specific leaf area (SLA) and empirical as well as theoretical leaf-specific conductivity (, ). The carbon isotope signature (δ13C) was positively correlated with ABI (fast growth was associated with a frequent stomatal closure) but negatively with the P50 value (Figures 3A,B). The provenances with the highest growth rates (UA, SK, and BG) exhibited highest carbon isotope signature (r2 = 0.29, P < 0.05; Figure 3A). The foliar C:N ratio was positively related to vessel diameter (D) and the theoretical specific conductivity (; Figures 4A,B), indicating that higher growth rates (associated with a high C:N ratio) were indirectly related to xylem hydraulic properties of the branch wood.
Figure 3. Aboveground biomass increment (ABI, A) and the xylem pressure causing 50% loss of hydraulic conductivity (P50, B) in relation to the carbon isotope signature. Given values are means per provenance, for symbol definition see Figure 1.
Figure 4. The carbon to nitrogen ratio (C:N) in relation to vessel diameter (D, A) and theoretical specific conductivity (, B). Given values are means per provenance, for symbol definition see Figure 1.
Only recently, a growing number of studies have investigated the intraspecific genetic differentiation in ecophysiological traits between populations (e.g., Bresson et al., 2011; Lamy et al., 2011; Eilmann et al., 2014; Aranda et al., 2015; Schreiber et al., 2015). Despite this positive trend, common-garden studies that have reported anatomical or hydraulic traits are still scarce (Anderegg and Meinzer, 2015). Our study on ten European beech provenances native to different localities in Europe revealed little intraspecific variation of the ecophysiological traits covered. In agreement with former common-garden experiments we observed no significant genetic differences between populations for the xylem pressure causing 50% loss of conductivity. The same pattern has been described for conifers (Sáenz-Romero et al., 2013; Lamy et al., 2014) as well as for beech (Wortemann et al., 2011, but see Aranda et al., 2015). Contrary to these studies, however, we observed significant genetic differentiation for the xylem pressure inducing 88% loss of hydraulic conductance (P88), which corresponds to the threshold for catastrophic hydraulic failure leading to irreversible drought-induced dysfunction in angiosperms (Urli et al., 2013; Li et al., 2015). In our study, genetic differentiation was further found for the xylem anatomical trait vessel density and the foliar trait mean leaf size. Our first hypothesis postulating that wood hydraulic properties are under genetic control due to local adaptation is further supported by the observation that most wood anatomical and derived hydraulic traits as well as some foliar traits were closely related to the climate at provenance origin. These results may indicate that wood anatomical and hydraulic traits of beech reflect genetic differentiation between provenances although no significant intraspecific differences became evident. This is in line with the observation that several wood properties related to vessel size are predominantly under genetic predisposition (Eilmann et al., 2014).
However, the high intra-population variability in all ecophysiological traits covered, which was on average two to three times higher than inter-population variability, is also a strong evidence for the high genetic diversity within provenances. For example, Aranda et al. (2015) stated that beech exhibits a high degree of intra-population genetic variability for embolism resistance. This indicates that European beech might have a high potential capacity to adapt to climate change. In agreement hereon, both hydraulic efficiency and embolism resistance were unrelated to climatic aridity at the place of provenance origin and showed the highest variance within a given provenance (~65%), followed by the second highest variance between replicated blocks (~25%) due to micro-environmental effects including within-canopy variability, and only a comparatively low variance between populations (~10%). This likewise confirms that beech exhibits a high degree of intra-population genetic variability for embolism resistance. Surprisingly, provenances originating from drought prone habitats (e.g., ES, BG, and DE-BB) grown at our site do not represent the most drought tolerant provenances in terms of functional trait adaptation as expected (see Eilmann et al., 2014). This finding is in agreement with former studies quantifying the high phenotypic plasticity of embolism resistance across environmental gradients (Herbette et al., 2010; Wortemann et al., 2011; Schuldt et al., 2016). Moreover, the vessel density of the branch xylem seems to play an outstanding role in the adaptation of the branch hydraulic system to different climates. Likewise with embolism resistance, vessel density varied independently of climatic conditions at the place of origin but was significantly different among provenances presumably as a consequence of differences in branch height. European beech trees may accordingly hold the capacity to adapt to local climate conditions by modifying their branch hydraulic traits primarily through adjustments in the relative abundance of vessels. This assumption is further supported by the close relation between vessel density and productivity as tree growth performance is hypothesized to scale with hydraulic efficiency (Tyree, 2003), which contrary to our results has been confirmed for both temperate and tropical tree species (Hajek et al., 2014; Hoeber et al., 2014; Kotowska et al., 2015).
It has been suggested that embolism resistance decreases with growth rate due to conflicting carbon allocation either to the construction of thicker cell walls, or to the building of foliar and axial tissues destined to increase canopy carbon gain (Cochard et al., 2007). However, empirical data from different species or genotypes do not unequivocally support this trade-off. For example, Fichot et al. (2010) reported on embolism-resistant genotypes of poplar which grew faster than more vulnerable genotypes. We observed no relation between P50 in branches and aboveground biomass increment across the provenances contradicting the hypothesized relation between embolism resistance and growth, which should trade-off with hydraulic efficiency. Several recent studies also failed to detect this relationship between embolism resistance and growth rate (Sterck et al., 2012; Hajek et al., 2014; Guet et al., 2015), indicating that embolism resistance is partly decoupled from hydraulic efficiency and biomass production (Fichot et al., 2015). In agreement hereon and in line with our results, hydraulic safety seems mostly decoupled from hydraulic efficiency, both at the inter-specific (Gleason et al., 2016) and intra-specific level (Schuldt et al., 2016; but see Hajek et al., 2014). In addition to these anticipated trade-offs it is necessary to include the foliage functionality for a holistic understanding of growth-related processes (Carlquist, 2012). In our sample, the P50 and P88 values decreased in parallel with increasing carbon isotope signature. The most embolism resistant provenances thus presumably had to close their stomata more frequently at the expense of significantly reduced carbon assimilation. This observation is contradicting observations of the functional coordination between vulnerability to xylem cavitation and the regulation of stomatal conductance (Sparks and Black, 1999; Brodribb et al., 2002). Either a sensitive stomatal regulation or simply a lower water demand as a consequence of smaller stature of drought prone provenances may explain this trait interrelation.
Moreover, the P50 and P88 values were not related to xylem anatomical traits like vessel diameter or vessel density as evidenced in several studies among and across species (Domec et al., 2010; Hajek et al., 2014). In contrast, the range of variation was small in our study and might have hampered a significant relationship as likewise described by other authors who also failed to detect a relation between vessel diameter and embolism resistance in closely related genotypes or different hybrids of poplar (Cochard et al., 2007; Fichot et al., 2010). These findings support the growing evidence that variation in P50 is mainly determined by other wood structural properties such as the topology of the xylem network (Loepfe et al., 2007; Martínez-Vilalta et al., 2012) as well as the pit membrane morphology (Plavcova et al., 2013; Li et al., 2016). The relation with vessel size in studies covering a sufficiently steep range in diameters is thus rather indirect as already speculated by Tyree and Sperry (1989).
In our study, provenances originating from drought-prone habitats (e.g., Spain) were most vulnerable to xylem cavitation. This contradicts the general observation of a less vulnerable structure for provenances originating from drought-prone sites (Rose et al., 2009; Robson et al., 2012). Surprisingly, the Spanish provenance developed the largest vessels of all ten provenances under these favorable environmental conditions contrary to our expectation. Due to the genetic control of this trait, beech provenances adapted to drier climates seem to over-develop their vascular system when grown in more humid environments. However, this did not translate into a better growth performance of these provenances.
Our study on ten beech provenances from all over Europe revealed that most wood anatomical and derived hydraulic traits and some foliar traits are under genetic predisposition according to significant differentiation between provenances or the relation with climatic aridity at the place of provenance origin. Although a certain degree of genetic differentiation was observed for embolism resistance, the high ratio of intra- vs. inter-population variance suggests that the genotypic variability of this trait is high between genotypes of a given provenance. In European beech, this high adaptive capacity in xylem function seems predominantly to be a consequence of adjusting vessel number but not necessarily vessel size. Nevertheless, we could not confirm the anticipated trade-off between hydraulic efficiency, xylem safety and growth in agreement with several recent studies. It further remains unclear if the observed high adaptive capacity of young European beech tree individuals can be extrapolated to old-growth forest trees, potentially enabling them to withstand increased drought by a flexible hydraulic response.
BS and GvW designed the study, PH, DK, and BS collected the field samples, DK and PH performed the hydraulic, wood anatomical and leaf morphological measurements and PH, BS, and SD analyzed the data and performed the statistical analyses. PH and BS wrote the first version of the manuscript, which was intensively discussed and revised by all authors.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank Roman Link for statistical advice, Ana Sapoznikova and all other assistants for their invaluable contributions to the experiment and the two reviewers, who provided helpful suggestions for improving the manuscript. The work has been financially supported by the University of Göttingen granted to BS. We further acknowledge the support by the Open Access Publication Funds of the University of Göttingen.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.00791
Anderegg, W. R. L., and Meinzer, F. C. (2015). “Wood anatomy and plant hydraulics in a changing climate,” in Functional and Ecological Xylem Anatomy, ed U. Hacke (Berlin: Springer International Publishing), 235–253.
Aranda, I., Cano, F. J., Gasco, A., Cochard, H., Nardini, A., Mancha, J. A., et al. (2015). Variation in photosynthetic performance and hydraulic architecture across European beech (Fagus sylvatica L.) populations supports the case for local adaptation to water stress. Tree Physiol. 35, 34–46. doi: 10.1093/treephys/tpu101
Bréda, N., Huc, R., Granier, A., and Dreyer, E. (2006). Temperate forest trees and stands under severe drought: a review of ecophysiological responses, adaptation processes and long-term consequences. Ann. For. Sci. 63, 625–644. doi: 10.1051/forest:2006042
Bresson, C. C., Vitasse, Y., Kremer, A., and Delzon, S. (2011). To what extent is altitudinal variation of functional traits driven by genetic adaptation in European oak and beech? Tree Physiol. 31, 1164–1174. doi: 10.1093/treephys/tpr084
Brodribb, T. J., Holbrook, N. M., and Gutiérrez, M. V. (2002). Hydraulic and photosynthetic co-ordination in seasonally dry tropical forest trees. Plant Cell Environ. 25, 1435–1444. doi: 10.1046/j.1365-3040.2002.00919.x
Cochard, H., Casella, E., and Mencuccini, M. (2007). Xylem vulnerability to cavitation varies among poplar and willow clones and correlates with yield. Tree Physiol. 27, 1761–1767. doi: 10.1093/treephys/27.12.1761
Cochard, H., Damour, G., Bodet, C., Tharwat, I., Poirier, M., and Améglio, T. (2005). Evaluation of a new centrifuge technique for rapid generation of xylem vulnerability curves. Physiol. Plant 124, 410–418. doi: 10.1111/j.1399-3054.2005.00526.x
Domec, J.-C., Schäfer, K., Oren, R., Kim, H. S., and McCarthy, H. R. (2010). Variable conductivity and embolism in roots and branches of four contrasting tree species and their impacts on whole-plant hydraulic performance under future atmospheric CO2 concentration. Tree Physiol. 30, 1001–1015. doi: 10.1093/treephys/tpq054
Eilmann, B., Sterck, F., Wegner, L., de Vries, S. M., von Arx, G., Mohren, G. M., et al. (2014). Wood structural differences between northern and southern beech provenances growing at a moderate site. Tree Physiol. 34, 882–893. doi: 10.1093/treephys/tpu069
Fichot, R., Barigah, T. S., Chamaillard, S., Le Thiec, D., Laurans, F., Cochard, H., et al. (2010). Common trade-offs between xylem resistance to cavitation and other physiological traits do not hold among unrelated Populus deltoides × Populus nigra hybrids. Plant Cell Environ. 33, 1553–1568. doi: 10.1111/j.1365-3040.2010.02164.x
Fichot, R., Brignolas, F., Cochard, H., and Ceulemans, R. (2015). Vulnerability to drought-induced cavitation in poplars: synthesis and future opportunities. Plant Cell Environ. 38, 1233–1251. doi: 10.1111/pce.12491
Gleason, S. M., Westoby, M., Jansen, S., Choat, B., Hacke, U. G., Pratt, R. B., et al. (2016). Weak tradeoff between xylem safety and xylem-specific hydraulic efficiency across the world's woody plant species. New Phytol. 209, 123–136. doi: 10.1111/nph.13646
Guet, J., Fichot, R., Lédée, C., Laurans, F., Cochard, H., Delzon, S., et al. (2015). Stem xylem resistance to cavitation is related to xylem structure but not to growth and water-use efficiency at the within-population level in Populus nigra L. J. Exp. Bot. 66, 4643–4652. doi: 10.1093/jxb/erv232
Hajek, P., Leuschner, C., Hertel, D., Delzon, S., and Schuldt, B. (2014). Trade-offs between xylem hydraulic properties, wood anatomy and yield in Populus. Tree Physiol. 34, 744–756. doi: 10.1093/treephys/tpu048
Herbette, S., Wortemann, R., Awad, H., Huc, R., Cochard, H., and Barigah, T. S. (2010). Insights into xylem vulnerability to cavitation in Fagus sylvatica L.: phenotypic and environmental sources of variability. Tree Physiol. 30, 1448–1455. doi: 10.1093/treephys/tpq079
Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G., Jarvis, A., and others (2005). Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 25, 1965–1978. doi: 10.1002/joc.1276
Hoeber, S., Leuschner, C., Köhler, L., Arias-Aguilar, D., and Schuldt, B. (2014). The importance of hydraulic conductivity and wood density to growth performance in eight tree species from a tropical semi-dry climate. For. Ecol. Manag. 330, 126–136. doi: 10.1016/j.foreco.2014.06.039
Jump, A. S., Hunt, J. M., and Peñuelas, J. (2006). Rapid climate change-related growth decline at the southern range edge of Fagus sylvatica. Glob. Change Biol. 12, 2163–2174. doi: 10.1111/j.1365-2486.2006.01250.x
Kahle, H. P. (2006). “Impact of the drought in 2003 on intra-and inter-annual stem radial growth of beech and spruce along an altitudinal gradient in the Black Forest, Germany,” in TRACE—Tree Rings in Archaeology, Climatology and Ecology, Vol. 4: Proceedings of the Dendrosymposium 2005, eds H. I. Gärtner, M. Monbaron, and G. Schleser (Fribourg: Forschungszentrums Jülich), 152–164. Schriften des Forschungszentrums Jülich, Reihe Umwelt Vol. 61.
Kotowska, M. M., Hertel, D., Rajab, Y. A., Barus, H., and Schuldt, B. (2015). Patterns in hydraulic architecture from roots to branches in six tropical tree species from cacao agroforestry and their relation to wood density and stem growth. Front. Plant Sci. 6:191. doi: 10.3389/fpls.2015.00191
Lamy, J.-B., Bouffier, L., Burlett, R., Plomion, C., Cochard, H., and Delzon, S. (2011). Uniform selection as a primary force reducing population genetic differentiation of cavitation resistance across a species range. PLoS ONE 6:e23476. doi: 10.1371/journal.pone.0023476
Lamy, J.-B., Delzon, S., Bouche, P. S., Alia, R., Vendramin, G. G., Cochard, H., et al. (2014). Limited genetic variability and phenotypic plasticity detected for cavitation resistance in a Mediterranean pine. New Phytol. 201, 874–886. doi: 10.1111/nph.12556
Leuschner, C., Backes, K., Hertel, D., Schipka, F., Schmitt, U., Terborg, O., et al. (2001). Drought responses at leaf, stem and fine root levels of competitive Fagus sylvatica L. and Quercus petraea (Matt.) Liebl. trees in dry and wet years. For. Ecol. Manag. 149, 33–46. doi: 10.1016/S0378-1127(00)00543-0
Li, S., Feifel, M., Karimi, Z., Schuldt, B., Choat, B., and Jansen, S. (2015). Leaf gas exchange performance and the lethal water potential of five European species during drought. Tree Physiol. 36, 179–192. doi: 10.1093/treephys/tpv117
Li, S., Lens, F., Espino, S., Karimi, Z., Klepsch, M., Schenk, J., et al. (2016). Intervessel pit membrane thickness as a key determinant of embolism resistance in angiosperm xylem. IAWA J. 37, 152–171. doi: 10.1163/22941932-20160128
Lindner, M., Maroschek, M., Netherer, S., Kremer, A., Barbati, A., Garcia-Gonzalo, J., et al. (2010). Climate change impacts, adaptive capacity, and vulnerability of European forest ecosystems. For. Ecol. Manage. 259, 698–709. doi: 10.1016/j.foreco.2009.09.023
Loepfe, L., Martinez-Vilalta, J., Piñol, J., and Mencuccini, M. (2007). The relevance of xylem network structure for plant hydraulic efficiency and safety. J. Theor. Biol. 247, 788–803. doi: 10.1016/j.jtbi.2007.03.036
Martínez-Vilalta, J., Mencuccini, M., Álvarez, X., Camacho, J., Loepfe, L., and Piñol, J. (2012). Spatial distribution and packing of xylem conduits. Am. J. Bot. 99, 1189–1196. doi: 10.3732/ajb.1100384
Peuke, A. D., Schraml, C., Hartung, W., and Rennenberg, H. (2002). Identification of drought-sensitive beech ecotypes by physiological parameters. New Phytol. 154, 373–387. doi: 10.1046/j.1469-8137.2002.00400.x
Plavcova, L., Jansen, S., Klepsch, M., and Hacke, U. G. (2013). Nobody's perfect: can irregularities in pit structure influence vulnerability to cavitation? Front. Plant Sci. 4:453. doi: 10.3389/fpls.2013.00453
R Development Core Team (2011). R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing. Available onlie at: http://www.R-project.org
Robson, T. M., Sánchez-Gómez, D., Cano, F. J., and Aranda, I. (2012). Variation in functional leaf traits among beech provenances during a Spanish summer reflects the differences in their origin. Tree Genet. Genomes 8, 1111–1121. doi: 10.1007/s11295-012-0496-5
Rose, L., Leuschner, C., Köckemann, B., and Buschmann, H. (2009). Are marginal beech (Fagus sylvatica L.) provenances a source for drought tolerant ecotypes? Eur. J. For. Res. 128, 335–343. doi: 10.1007/s10342-009-0268-4
Sáenz-Romero, C., Lamy, J.-B., Loya-Rebollar, E., Plaza-Aguilar, A., Burlett, R., Lobit, P., et al. (2013). Genetic variation of drought-induced cavitation resistance among Pinus hartwegii populations from an altitudinal gradient. Acta Physiol. Plant 35, 2905–2913. doi: 10.1007/s11738-013-1321-y
Schreiber, S. G., Hacke, U. G., and Hamann, A. (2015). Variation of xylem vessel diameters across a climate gradient: insight from a reciprocal transplant experiment with a widespread boreal tree. Funct. Ecol. 29, 1392–1401. doi: 10.1111/1365-2435.12455
Schuldt, B., Knutzen, F., Delzon, S., Jansen, S., Müller-Haubold, H., Burlett, R., et al. (2016). How adaptable is the hydraulic system of European beech in the face of climate change related precipitation reduction? New Phytol. 210, 443–458. doi: 10.1111/nph.13798
Sparks, J. P., and Black, R. A. (1999). Regulation of water loss in populations of Populus trichocarpa: the role of stomatal control in preventing xylem cavitation. Tree Physiol. 19, 453–459. doi: 10.1093/treephys/19.7.453
Sperry, J. S., Nichols, K. L., Sullivan, J. E. M., and Eastlack, S. E. (1994). Xylem embolism in ring-porous, diffuse-porous, and coniferous trees of Northern Utah and Interior Alaska. Ecology 75, 1736–1752. doi: 10.2307/1939633
Sterck, F. J., Martínez-Vilalta, J., Mencuccini, M., Cochard, H., Gerrits, P., Zweifel, R., et al. (2012). Understanding trait interactions and their impacts on growth in Scots pine branches across Europe. Funct. Ecol. 26, 541–549. doi: 10.1111/j.1365-2435.2012.01963.x
Stojnić, S., Orlović, S., Miljković, D., Galić, Z., Kebert, M., and von Wuehlisch, G. (2015). Provenance plasticity of European beech leaf traits under differing environmental conditions at two Serbian common garden sites. Eur. J. For. Res. 134, 1109–1125. doi: 10.1007/s10342-015-0914-y
Urli, M., Porté, A. J., Cochard, H., Guengant, Y., Burlett, R., and Delzon, S. (2013). Xylem embolism threshold for catastrophic hydraulic failure in angiosperm trees. Tree Physiol. 33, 672–683. doi: 10.1093/treephys/tpt030
van der Werf, G. W., Sass-Klaassen, U. G. W., and Mohren, G. M. J. (2007). The impact of the 2003 summer drought on the intra-annual growth pattern of beech (Fagus sylvatica L.) and oak (Quercus robur L.) on a dry site in the Netherlands. Dendrochronologia 25, 103–112. doi: 10.1016/j.dendro.2007.03.004
von Wuehlisch, G., Liesebach, M., Muhs, H. J., and Stephan, R. (1998). “A network of international beech provenance trials,” in First EUFORGEN Meeting on Social Broadleaves, 23-25 October 1997, Bordeaux, France, eds J. Turok, A. Kremer, and S. De Vries (Rome: International Plant Genetic Resources Institute), 164–172.
Willigen, C. V., and Pammenter, N. W. (1998). Relationship between growth and xylem hydraulic characteristics of clones of Eucalyptus spp. at contrasting sites. Tree Physiol. 18, 595–600. doi: 10.1093/treephys/18.8-9.595
Wortemann, R., Herbette, S., Barigah, T. S., Fumanal, B., Alia, R., Ducousso, A., et al. (2011). Genotypic variability and phenotypic plasticity of cavitation resistance in Fagus sylvatica L. across Europe. Tree Physiol. 31, 1175–1182. doi: 10.1093/treephys/tpr101
Wutzler, T., Wirth, C., and Schumacher, J. (2008). Generic biomass functions for Common beech (Fagus sylvatica) in Central Europe: predictions and components of uncertainty. Can. J. For. Res. 38, 1661–1675. doi: 10.1139/X07-194
Zimmermann, J., Hauck, M., Dulamsuren, C., and Leuschner, C. (2015). Climate warming-related growth decline affects Fagus sylvatica, but not other broad-leaved tree species in Central European mixed forests. Ecosystems 1, 560–572. doi: 10.1007/s10021-015-9849-x
Zomer, R. J., Trabucco, A., Bossio, D. A., and Verchot, L. V. (2008). Climate change mitigation: a spatial analysis of global land suitability for clean development mechanism afforestation and reforestation. Agric. Ecosyst. Environ. 126, 67–80. doi: 10.1016/j.agee.2008.01.014
Keywords: adaptive capacity, Fagus sylvatica L., genetic variability, hydraulic conductivity, leaf morphology, phenotypic plasticity, provenance trial, vulnerability to cavitation
Citation: Hajek P, Kurjak D, von Wühlisch G, Delzon S and Schuldt B (2016) Intraspecific Variation in Wood Anatomical, Hydraulic, and Foliar Traits in Ten European Beech Provenances Differing in Growth Yield. Front. Plant Sci. 7:791. doi: 10.3389/fpls.2016.00791
Received: 28 February 2016; Accepted: 22 May 2016;
Published: 15 June 2016.
Edited by:Sergio Rossi, Université du Québec à Chicoutimi, Canada
Reviewed by:Tommaso Anfodillo, University of Padova, Italy
Jordi Martínez-Vilalta, Autonomous University of Barcelona, Spain
Copyright © 2016 Hajek, Kurjak, von Wühlisch, Delzon and Schuldt. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bernhard Schuldt, firstname.lastname@example.org