http://orcid.org/0000-0003-4970-8539
Figures
Abstract
The coniferous forest tree Douglas-fir (Pseudotsuga menziesii) is native to the pacific North America, and is increasingly planted in temperate regions worldwide. Nitrogen (N) metabolism is of great importance for growth, resistance and resilience of trees. In the present study, foliar N metabolism of adult trees of three coastal and one interior provenance of Douglas-fir grown at two common gardens in southwestern Germany (Wiesloch, W; Schluchsee, S) were characterized in two subsequent years. Both the native North American habitats of the seed sources and the common garden sites in Germany differ in climate conditions. Total and mineral soil N as well as soil water content were higher in S compared to W. We hypothesized that i) provenances differ constitutively in N pool sizes and composition, ii) N pools are affected by environmental conditions, and iii) that effects of environmental factors on N pools differ among interior and coastal provenances. Soil water content strongly affected the concentrations of total N, soluble protein, total amino acids (TAA), arginine and glutamate. Foliar concentrations of total N, soluble protein, structural N and TAA of trees grown at W were much higher than in trees at S. Provenance effects were small but significant for total N and soluble protein content (interior provenance showed lowest concentrations), as well as arginine, asparagine and glutamate. Our data suggest that needle N status of adult Douglas-fir is independent from soil N availability and that low soil water availability induces a re-allocation of N from structural N to metabolic N pools. Small provenance effects on N pools suggest that local adaptation of Douglas-fir is not dominated by N conditions at the native habitats.
Citation: Du B, Kreuzwieser J, Dannenmann M, Junker LV, Kleiber A, Hess M, et al. (2018) Foliar nitrogen metabolism of adult Douglas-fir trees is affected by soil water availability and varies little among provenances. PLoS ONE 13(3): e0194684. https://doi.org/10.1371/journal.pone.0194684
Editor: Berthold Heinze, Austrian Federal Research Centre for Forests BFW, AUSTRIA
Received: September 8, 2017; Accepted: March 7, 2018; Published: March 22, 2018
Copyright: © 2018 Du et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: This work was financially supported by the Deutsche Forschungsgemeinschaft (www.dfg.de; grants KR 2010/4-1 to JK, EN829/5-1 to IE and GE1090/7-1 to AG), the Forest Research Institute of the German Federal State Baden-Württemberg (FVA, to IE, UK, HW), the Natural Sciences and Engineering Research Council of Canada (NSERC) and the Swiss National Science Foundation SNF (31003A_159866 to AG). BD obtained funding from the China Scholarship Council, the Ministerium für Wissenschaft, Forschung und Kunst Baden-Württemberg of Germany and the University of Freiburg. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Douglas-fir (Pseudotsuga menziesii) is native to America’s Pacific Northwest (USA, Canada), where it grows in more than 20 million ha of natural forests [1]. Reflecting the great natural range of distribution in longitude, latitude and altitude (from California to British Columbia, and from sea level up to an altitude of >2,000 m), a variety of different ecotypes evolved. These ecotypes are locally adapted, i.e. phenotypes of adaptive traits vary in a way providing highest fitness in the native environments [2–4].
Based on knowledge on intraspecific variation in stress responses and growth, suitable Douglas-fir provenances have been increasingly planted outside North America, e.g. in Central Europe [5,6], New Zealand, Australia, and South America [6]. Available studies on intraspecific variation of Douglas-fir focused on growth [7–11], phenology [12], and photosynthetic gas exchange [13,14], but little is known about natural genetic variation in nitrogen (N) metabolism in field-grown Douglas-fir. N is an important factor for the growth and vitality of trees [15–19]. Differences in climatic, edaphic and biotic factors across large natural ranges of conifer tree species presumably translate into variation in N cycling and, thus, N availability at the native habitats [20–22]. In turn, variation in N availability may cause local adaptation at the level of N uptake [22,23] or metabolism, which is reflected in constitutive variation in N pools among provenances as reported for Norway spruce [20] and Scots pine [24,25]. In a provenance trial with Scots pine, highest needle arginine concentrations were found in the northernmost provenances [25], possibly as a consequence of nutrient imbalances or as an adaptation enabling a C-efficient assimilation of plant available N not required for growth [26–28]. Higher needle total N contents of Norway spruce [20] and Scots pine [24] at higher altitudes or longitudes, or at sites with a cooler climate were concluded to be adaptive [20,24], since they often go along with larger quantities of Rubisco [29,30]. This, in turn, might enable higher net C assimilation rates, which might compensate for environmental limitations for carbon gain.
Adaptation to local environments is very common in trees, and adaptive phenotypic variation is often reflected in genetic differentiation among populations [31]. If genetic differentiation in loci underlying phenotypic variation is constrained by high gene flow, or in temporally highly variable environments, phenotypic plasticity in adaptive traits can confer fitness advantages [32]. Phenotypic plasticity, i.e. the ability of a genotype to express different phenotypes in different environments, may -or may- not vary among genotypes (or populations) [33]. Benomar et al. [34] did not observe variation in plasticity of adaptive traits such as height growth, photosynthetic CO2 assimilation and leaf N concentration among seed sources of Picea glauca. In Douglas-fir, plasticity in height growth varied among provenances (i.e. the provenance-environment interaction effect was significant) [11], but we lack knowledge on variation in plasticity of N metabolism in Douglas-fir.
Besides variation in N pools as a consequence of genotypic differentiation or genotypic variation in plasticity, N metabolism of trees is also affected by acclimatization to variation in environmental conditions [35,36]. Water shortage and elevated temperature exert strong effects on tree N metabolism [37,38]. Under drought and / or heat, foliar total amino acid contents, particularly amino acids from the glutamate family (e.g. glutamate, arginine), increased in deciduous [39,40] and conifer species [35,41–43]. Accumulation of arginine in needles was regarded as a general response of conifers to stresses causing a reduction in growth [26]. For Pinus halepensis, a decreased content of soluble proteins under drought was reported [44], a response reflecting either a shift in N metabolism towards accumulation of osmolytic amino compounds, and / or a reduction of maintenance and synthesis of proteins [45].
In the present study, we characterized the N metabolism of 50-year-old Douglas-fir trees of four provenances grown at two field sites of a common garden experiment located in Germany [46]. The field sites were established in 1961 and differ in annual temperature and precipitation. Trees of one interior and three coastal provenances were grown from seeds collected at North American native habitats differing in precipitation and temperature conditions. The assumption of local adaptation to conditions at the native habitat is supported by significant variation in photosynthetic gas exchange [13] and growth [11] among the provenances studied here. Genetic differentiation of the studied provenances was reported earlier based on both microsatellite [11] and SNP marker analysis [47]. This experimental setting and repeated sampling over the course of two years allowed to study the extent of phenotypic variation in N pools and to attribute this variation to i) differences in environmental conditions among sites and within sites over time, ii) provenance effects reflecting genetic differentiation, or iii) provenance-environment interaction effects (i.e. variation in phenotypic plasticity among provenances). We hypothesized that provenances show constitutive differences in sizes and / or composition of N pools; specifically, we expected that as a consequence of a cooler and dryer climate at their native habitat, the interior provenance shows a higher content of metabolic N including total soluble proteins and highly abundant N-rich amino acids such as arginine and asparagine. We further hypothesized that sizes and composition of N pools are affected by variation in environmental conditions at the common garden sites; in particular, we expected a reduction of structural N and total soluble proteins, and an increase of total amino acids and arginine under reduced soil water availability. Lastly, we hypothesized that effects of environmental factors on N pools differ among provenances; in this context, we expected that the interior provenance shows a smaller extent of re-allocation of structural N towards amino acids due to constitutively higher levels of metabolic N. To analyze the trees’ N metabolism, we determined the total N content as well as the partitioning of N into amino acids, soluble protein and structural N in needles of Douglas-fir provenances by N metabolite profiling.
Materials and methods
Plant material and experimental design
Douglas-fir (Pseudotsuga menziesii (Mirbel) Franco) trees of three coastal provenances (var. menziesii; ‘WA Darrington 3 Conrad Creek’ (CR), ‘BC Cameron Lake’ (LA), ‘OR Santiam River’ (RI)) and one provenance of the interior variety (var. glauca; ‘BC Salmon Arm’ (AR)) were studied. For geographical maps of the origin of these four provenances and the two field sites see Müller et al. [47]. Selection criteria for the provenances were 1) presence on both study sites and 2) high variation in climatic conditions of their natural sites (Table 1) [46]. Long-term averages (period 1961–1990) of yearly climatic moisture deficit (CMD) reveal that the interior provenance AR must acquire more water from other sources than rain (i.e. soil moisture) to avoid the impact of drought stress compared to coastal provenances [48]. CMD data for sites of origin were retrieved from ClimateNA version 5.50 for CR and RI, and from ClimateBC version 5.50 for AR and LA [49].
The trees studied were grown at two field sites in Germany, i.e. “Wiesloch” (“W”) and “Schluchsee” (“S”), that had been established in 1961 as part of an international provenance trial [46]. Seeds collected for this trial originated from open-pollinated trees in naturally regenerated stands in British Columbia (Canada), Washington, and Oregon (USA) [46] as previously described by Neophytou et al. [11]. At each field site, each of the selected provenances had been planted at a density of 3,300 trees ha-1 on two to four plots (each 0.1 ha). The plots had been arranged irregularly, so that identical provenances did not adjoin [11,50]. A height-driven thinning regime has been applied (see Kohnle et al. [50] for details). At the end of the growing season 2011, trees reached heights between 26.6±0.3 m (AR at S) and 33.7±0.3 m (CR at W) [11]. Annual precipitation, precipitation during the vegetation period, mean annual temperature and CMD differ between both sites (Table 2). Long-term averages (period 1961–1990) of CMD are within the range of the native habitats (site W) (Table 1) or lower (site S). CMD data for the field sites in Germany were retrieved with the ClimateEU version 4.63 software package (based on the methodology described by [51]), available at http://tinyurl.com/ClimateEU.
Weather data in one-hour-resolution for the growing seasons of 2010 and 2011 were available from a private meteorological station at Schluchsee (N47°49’16”, E8°11’08”) and from Germany’s National Meteorological Service (Deutscher Wetterdienst) station Kirrlach-Waghäusel (N49°15’0”, E8°32’24”) which is located in close proximity to the site Wiesloch. Soil water availability was modelled in daily resolution using the forest hydrological water budget model WBS3 that estimates daily total available soil water (TAW) using temperature, precipitation, latitude, soil type, plant cover, slope, and slope aspect [54].
Soil sampling and chemical analyses
Soil samples were collected in two layers, 0–10 cm and 10–20 cm, using a hand auger at the studied plots in W and S in August 2013. The litter was manually removed before sampling. Soil samples were stored in plastic bags at 4°C before passed through a sieve (mesh width 2 mm) for further analyses. Half of the samples were dried at 60°C in a drying oven until reaching constant weight. Dried soil was powdered in a ball mill (MM400, RETSCH GmbH, Haan, Germany) and used for analysis of total N and carbon (C) concentrations as well as for analysis of δ15N. The samples were weighed and placed in tin capsules (5 × 9 mm; IVA Analysentechnik e.K., Meerbusch, Germany). Analyses were conducted using a Costech Elemental Analyzer (Costech International S.p.A., Milano, Italy) equipped with a zero-blank auto-sampler and coupled via a ConFloIII to a Thermo Finnigan Delta V Plus isotope ratio mass spectrometer (Thermo Scientific, Waltham, MA, USA) as described in detail by Liu et al. [55].
Ammonium and nitrate were extracted using field-fresh moist soil according to the method described by Schlotter et al. [56]. Fresh soil samples were extracted with distilled water (1:2, w/w) and shaken at 4°C for 24 hours. Extracted solutions were then centrifuged and the supernatants were filtered. Nitrate contents were analyzed with an ion chromatograph (DX 120, Dionex, Idstein, Germany) connected to an autosampler (AS 3500, Thermo Separation Products, Piscataway, NJ, USA). Concentrations were calculated using external nitrate standards. Ammonium was analyzed photometrically (UV-DU 650 spectrophotometer, Beckman Coulter Inc., Fullerton, CA, USA) at 650 nm as based on the color reaction between ammonium and a weakly alkaline mixture of sodium (Na) salicylate and Na hypochlorite as a chlorine source in the presence of Na nitroprusside [57]. Oven-dried ammonium sulfate was used as standard for quantification.
Sampling and processing of needles
Needles of five to six individual trees per provenance were sampled. All trees sampled were part of the original establishment of the site in 1961. Previous year needles were sampled in late spring (W: May 12th 2010, May 12th 2011; S: May 27th 2010, May 26th 2011) and in summer (W: July 14th 2010, July 14th 2011; S: July 28th 2010, July 28th 2011) at each site. Campaigns at W were carried out prior to campaigns at S to account for the later start of the growing season at S [13]. Needles were sampled from the south-facing upper canopy of the trees using hydraulic platforms, pre-installed ropes or a shotgun in the early afternoon between 1 pm and 3 pm. Needles were rapidly excised from the twigs, immediately shock-frozen in liquid nitrogen, homogenized to a fine powder under liquid nitrogen and subsequently stored at -80 °C until analysis. For determination of needle water content, aliquots of 100 mg frozen powder were transferred to 1.5 ml reaction tubes and were dried at 60°C until reaching weight constancy. Needle water content was calculated as the difference between fresh and dry weight.
Determination of N pools and compounds in needles
Total N, total soluble protein (TSP), amino compounds and ammonium in needles were determined as previously described [58]. Total N concentration in oven-dried needle powder was determined after Kjeldahl digestion and analyzed on a AT200 spectrophotometer (Beckman Coulter / Olympus, USA) at 630 nm. TSP content was analyzed as previously described by Dannenmann et al. [59] with some adaptations for Douglas-fir needles [58]. For extraction, 35 mg frozen needle powder was added to a mixture of 1.5 ml extraction buffer (50 mM Tris–HCl, 1 mM ethylenediaminetetraacetic acid, 15% glycerol (v/v), 1 mM phenylmethylsulfonyl fluoride, 5 mM dithiothreitol, 0.1% Triton X-100; pH 8.0) and 150 mg polyvinylpyrrolidone. After 30 min incubation with shaking at 4°C followed by centrifugation at 14,000 g for 15 min at 4°C, 500 μl aliquots of the supernatant were transferred to new tubes. The pellet was extracted again with 0.5 ml extraction buffer, and the supernatants were combined thereafter. To precipitate dissolved proteins, 10% (w/v) trichloroacetic acid was added to the mixed supernatant in the same volume as the combined supernatants. The mixture was incubated for 10 min at 4°C. After centrifugation, the pellets were dissolved in 0.5 ml 1 M KOH. Five μl aliquots of the protein extract were mixed with 200 μl Bradford reagent (Amresco Inc., Solon, Ohio, USA). Absorbance at 595 nm was measured by a Sunrise Microplate Reader (Tecan Austria GmbH, Grödig, Austria) and calibrated with bovine serum albumin (BSA; Sigma-Aldrich Chemie GmbH, Germany) standards.
Amino acids and ammonium were analyzed using an ultra-performance-liquid-chromatography (UPLC) system (Waters Corp., Milford, MA) as described by Winter et al. [60] and Luo et al. [61]. Approximately 50 mg ground plant material was extracted in 1 ml chloroform/methanol (3/7, v/v) and 0.2 ml extraction buffer (20 mM Hepes, 5 mM EGTA, 10 mM NaF, pH 7.0) on ice. After 30 min incubation, 600 μl ddH2O were added to the samples, mixed, and centrifuged for 10 min, 14,000 g at 4°C. Aliquots of supernatants were transferred to new tubes, and 600 μl ddH2O were again added to the first extraction tube followed by 10 min centrifugation. Supernatants from the combined aqueous phases were freeze-dried (Alpha 2–4, Christ, Osterode, Germany) for 3 d. Amino acids and ammonium were separated on an AccQ-Tag Ultra column (2.1×100 mm; Waters) using a gradient of AccQ-Tag Ultra eluents A and B (Waters) at a flow rate of 0.7 ml min-1. The absorbance at 260 nm was detected and quantified using amino acid and ammonium standards.
Structural N and total amino acids calculation in needles
Structural N content in needles was calculated by subtracting the N fractions of soluble proteins, amino acids and ammonium from total N. According to Junker et al. [13], N amounts of chlorophyll in Douglas-fir needles were as low as 1.8% and, therefore, were neglected in the calculation. Total amino acids (TAA) contents were calculated as the sum of the concentrations of individual amino acids determined as described above.
Statistical analysis of data from soil and needle samples
For an initial exploratory analysis of needle N compounds, the publicly available platform “MetaboAnalyst 3.0” (http://www.metaboanalyst.ca/MetaboAnalyst/) [62] was used for Partial Least Squares Discriminant Analysis (PLS_DA). Raw data were cube root transformed and auto-scaled (mean-centered and divided by the standard deviation) because of different ranges and variances of the variables. Missing values were processed with an inbuilt function of MetaboAnalyst which excluded features with >50% missing values and replaced the remaining missing values by the minimum value of the same feature.
For total N, TSP, structural N, TAA, and the three most abundant amino acids arginine, glutamate, and asparagine potential effects of provenance and environmental factors were analysed by linear modelling using the statistical computing environment R v3.3.1 [63]. Based on graphical inspection of the distribution of the data, modelling was done with untransformed data of total N, structural N and glutamate, log2-transformed data of TAA, arginine and asparagine, and square-root transformed data of total soluble protein. As suggested by Zuur et al. [64], we started modelling with a possibly beyond-optimal model including ‘provenance’, ‘site’, ‘TAW’, ‘year’ and ‘season’ as fixed main effects and the interaction effects ‘provenance:site’, ‘provenance:TAW’, ‘site:year’. We considered ‘TAW’ to represent variation in soil water availability within and between sites. To represent variation between sites not directly related to water availability, but to e.g. nutrient concentrations or long-term acclimatization to site-specific environments, we included ‘site’ as a predictor in addition to ‘TAW’. We included ‘season’ (encoded as daylength of sampling day centered to summer solstice as described in Hess et al. [65]) as a predictor to account for effects brought about e.g. by photoperiod, which are independent from ‘TAW’ and ‘site’. ‘Year’ was included as variable integrating variation in environmental conditions over the year of study.
With the above fixed effects model structure, we compared models with a random intercept for tree (function ‘lme’, R package ‘nlme’; [66]) against models without random intercept (function ‘gls’) by BIC and LRT-test. For all dependent variables, models without random effect provided a better fit to the data (dBIC >2). Consequently, the selection of the fixed effects structure was performed on the models without random effect (function ‘gls’) by iteratively dropping non-significant fixed effects (P > 0.01) with highest P-values as computed with the function ‘anova’. This procedure was iterated until all explanatory variables were significant (P < 0.01).
To calculate variance proportions of each significant explanatory variable, the final models were re-fitted using the function ‘lm’, and the proportion of explained variance of a given predictor was calculated by dividing the sum of squares of that predictor by the total sum of squares. Normal distribution and homogeneity of variance of the residuals of the final models were confirmed graphically.
To avoid confounding effects of ‘site’, ‘season’ or ‘year’ on effects of ‘TAW’, we analysed their degree of collinearity. Collinearity of ‘year’ and ‘season’ with ‘TAW’ was small (r2 = 0.1/r = 0.31 and r2 = 0.16/r = 0.4, respectively), while collinearity of ‘TAW’ and ‘site’ was higher (r2 = 0.54/r = 0.73). This prompted us to calculate variance inflation factors (VIFs; calculated with function ‘VIF’, package ‘fmsb’ [67]), which are a commonly used indicator for inaccurate modelling results as a consequence of collinearity (VIFs > 10 indicate distorted model estimation; [68]). VIFs for the starting models and final models ranged from 1.2 to 4.0 (see Results section), indicating reliable model estimation. To further dissect effects of ‘TAW’ and effects of ‘site’ independent from water availability, we used sequential regression. This method is recommended to compute predictors that are completely independent [68]. To achieve this, explanatory variables are regressed against each other [68]. We considered ‘TAW’ as the most direct and obvious indicator for water availability and thus regressed site on ‘TAW’ (‘site’ ~ ‘TAW’) using a generalized linear model (function ‘glm’, binomial family, logit link). The residuals of this regression are orthogonal to ‘TAW’ are thus considered to represent site effects not directly related to water availability. We then modelled (function ‘lm’) N compounds as functions of ‘TAW’ and the residuals (on response scale) of ‘site’ ~ ‘TAW’ (referred to as resid (‘site’)): N compound ~ ‘TAW’ + resid (‘site’). N compounds were transformed as described above for standard multiple regression.
Statistical analysis of differences in soil N and C contents between the two field sites was done by Student’s T test using the software package SigmaPlot 11.0 (Systat Software GmbH, Erkrath, Germany), and raw data were transformed by denary logarithm to match normal distribution if necessary.
Results
Soil nitrogen and environmental conditions during sampling campaigns varied considerably between field sites
Contents of total C, total N, nitrate and ammonium in soil from S were considerably higher than from W in both 0–10 cm and 10–20 cm layers, with the exception of ammonium contents in the upper layer (Table 3). The C / N ratios of soil at W were significantly higher than at S in both layers. In the upper layer, soil δ15N at S was significantly higher than at W, but no significant difference was found between the two sites in the deeper layers (Table 3).
The weather conditions during and before the field campaigns varied considerably between the two sites and the years (Table 4). In 2010, yearly CMD of both sites (i.e., 50 mm at S and 207 mm at W) was within the range of the long-term average ± 1 standard deviation, while 2011 (i.e., 97 mm at S and 325 mm at W) was drier at both sites (compare Table 2). In 2010, the conditions were exceptionally warm and dry in July at W, but rather normal (compared to the meteorological reference period) at S and at W in May. Conditions at both sites were warmer and drier in May 2011 than usually. As a consequence of these dry and warm conditions in July 2011 (at W) and May 2011 (at S and W), total available soil water (TAW) was considerably decreased (Table 4). TAW was always higher in S compared to W (Table 4). Data on pre-dawn and midday twig water potential presented in [13] mirror site-related differences in TAW.
Exploratory analysis suggests effects of site and time-of-the-year on Douglas-fir nitrogen profiles
In a first approach to analyze effects of ‘site’, ‘provenance’ and environmental conditions on Douglas-fir N metabolism, we applied a PLS_DA to data of individual amino acids, total and structural N, TSP and TAA. This analysis clearly indicated a clustering of trees according to site in both years (Fig 1A and 1B). As ranked by Variable Importance in Projection (VIP) scores, TAA, soluble protein and arginine were among the five most important compounds determining the site-related clustering in both years (S1A and S1B Fig). Trees at W generally showed higher contents of these components as compared to site S (S1A and S1B Fig). In July 2010, all four provenances at W showed more than 3 times higher proline contents than at S, especially AR with up to 7 times higher contents (S1 Table).
Sampling was done in May and July 2010 (A) and 2011 (B). Figures represent score plots of a partial least squares–discriminant analysis. Douglas-fir provenances (AR, Salmon Arm; CR, Conrad Creek; LA, Cameron Lake; and RI, Santiam River) are represented by different symbols, open and solid symbols indicate individual trees at field sites Schluchsee (S) and Wiesloch (W) irrespective of month of sampling. Semi-transparent shadings indicate 95% confidence regions.
To investigate any possible clustering according to the month of sampling (May vs. July), PLS_DA were computed for data subsets split per site and year (Fig 2A–2D). A clear separation of samples collected in May and July was evident in all four data subsets. Important N compounds determining this clustering differed between years and sites (S2 Fig). Generally, contents of glutamine and proline were higher in May than July in both years irrespective of sites (S2 Fig and S1 Table).
Sampling was done in May and July 2010 (5/10 and 7/10) and 2011 (5/11 and 7/11). Figures represent score plots of a partial least squares–discriminant analysis computed per year and site. A) and C) show plots for Wiesloch (W), B) and D) for Schluchsee (S). Douglas-fir provenances (AR, Salmon Arm; CR, Conrad Creek; LA, Cameron Lake; and RI, Santiam River) are represented by different symbols, open and solid symbols indicate individual trees. Semi-transparent shadings indicate 95% confidence regions.
A strong clustering of provenances was not indicated by Figs 1 or 2. PLS_DA computed per site, year and month of sampling supported this finding (S3 and S4 Figs). Clustering of provenances was weak with no consistent patterning (S3 and S4 Figs). At W, there was a tendency that the interior provenance AR differed from the other–coastal–provenances in N profiles. AR showed higher contents of tyrosine and proline than the other three provenances in 2010 and May 2011 at site W (S1 Table).
Soil water availability and year of study explain most of the variation in needle N pool sizes and concentrations of major amino acids
In order to dissect effects of ‘provenance’, ‘year’, ‘season’, ‘TAW’ (soil water availability integrating environmental factors including temperature, precipitation, and latitude) and ‘site’, we fitted linear models to data on N pool sizes and the most abundant amino acids.
The total N content of Douglas-fir needles ranged between 13.0 and 19.7 mg g-1 DW (Fig 3A and 3B). Linear modelling revealed that there was a small but significant difference in total N content among the provenances, which was caused mainly by higher total N contents in CR and LA as compared to AR (Table 5 and S5A Fig). Most of the variation in total N content was explained by ‘TAW’ (23% of total variation), which is reflected in a significant and provenance-independent trend towards lower total N contents at higher TAW. In addition, total N contents were significantly higher at site W compared to site S, and there was a significant but very small effect of ‘season’, towards higher total N contents in May compared to July.
Data represent means ± SE of five to six replicates of four Douglas-fir provenances at the two field sites, Wiesloch (left panel) and Schluchsee (right panel), in May and July of 2010 (5/10, 7/10) and 2011 (5/11, 7/11). Stripes above panels indicate modeled total available soil water (TAW in %) at the sampling day.
Total soluble protein content of needles of the different Douglas-fir provenances showed a relatively large range from 13 to 63 mg g-1 DW (Fig 3C and 3D). Like total N, there was a small but significant effect of ‘provenance’ (Table 5 and S5B Fig), due to a lower TSP content in AR as compared to the other provenances. We detected a significantly higher TSP content when TAW was low. There was a small but significant effect towards lower TSP contents in S compared to W. Across all provenances and sites, TSP contents were significantly higher in 2011 compared to 2010, and this effect accounted for 44% of the variation in TSP content.
Structural N content did not differ significantly among the provenances. With increasing TAW, structural N showed a small, significant increase (Table 5 and S5C Fig). Structural N content was significantly higher at W compared to S, and higher in 2010 compared to 2011 on both sites. This variability between the two years under study accounted for 43% of the variation in structural N content.
TAA contents significantly increased from 20 to 85 μmol g-1 DW with decreasing TAW (Fig 3E and 3F, Table 5 and S5D Fig). This effect accounted for 36% of the variation in TAA. The effect of ‘TAW’ was to a small extent dependent on provenance with less impacts of TAW on LA compared to the other three provenances. TAA were significantly lower in trees at S as compared to W.
The three most abundant amino acids were asparagine, arginine and glutamate (Fig 4) which contributed with up to 93% to the TAA content as observed in provenance CR at W in July 2011. Asparagine contents ranged from 12.05 μmol g-1 DW in AR to 27.56 μmol g-1 DW in RI (Fig 4A and 4B). Asparagine contents were significantly affected by ‘TAW’ and the interaction (P<0.005) of ‘TAW’ with ‘provenance’ (Table 5 and S5E Fig). While provenances LA and CR were highly fluctuating and did not display a clear pattern to changes in TAW, asparagine contents in AR constantly decreased with increasing TAW, while asparagine contents of RI were unaffected by increasing TAW until it reached 80%.
Data represent means ± SE of five to six replicates of trees growing at field sites Wiesloch (left panel) and Schluchsee (right panel) in May and July 2010 (5/10, 7/10) and 2011 (5/11, 7/11). Stripes above panels indicate modeled total available soil water (TAW in %) at the sampling day.
Arginine contents in trees of the different provenances varied between concentrations below the limit of detection (about 0.13 μmol g-1 DW, half of the minimal detected value) and approximately 60 μmol g-1 DW (Fig 4C and 4D). This variation among provenances accounted for 9% of the total variance (Table 5 and S5F Fig), and was reflected in significantly lower arginine contents in provenance RI compared to the other provenances. Most of the variation (41%) in arginine contents was explained by the significant trend towards lower arginine concentration with increasing TAW. There was a small but significant interaction effect of ‘TAW’ and ‘provenance’, reflected in a less pronounced decrease of arginine concentrations with increasing TAW in provenance RI. Arginine concentrations were higher in W compared to S: At W, arginine contents ranged from 2.41 to 49.36 μmol g-1 DW. By contrast, at field site S, arginine contents showed a maximum of 3.24 μmol g-1 DW (CR, May 2011); in July 2010, arginine contents of all provenances were below the limit of detection. Similar to TSP and TAA contents, arginine contents were significantly higher in 2011 as compared to 2010, irrespective of site.
Glutamate contents varied only slightly among provenances, but showed a strong and significant trend towards lower glutamate contents with increasing TAW (Table 5 and S5G Fig). This effect accounted for 27% of the variation in glutamate contents. In W, glutamate contents were significantly higher as compared to S. Moreover, glutamate contents were significantly higher in 2011 as compared to 2010. However, the increase in glutamate contents in 2011 was smaller in S as compared to W, reflected in a small but significant ‘site:year’ interaction effect.
For all modelled N compounds, variance inflation factors were substantially smaller than 10 for both the starting models and the final models, indicating reliable model estimation (Table 5). We confirmed that effects of ‘TAW’ were not confounded by effects of ‘site’ by applying sequential regression (S2 Table). Significance and variance proportions of ‘TAW’ and site effects independent from ‘TAW’ estimated by sequential regressions (S2 Table) were very similar to the results obtained by classical multiple regression (Table 5), corroborating that the effects of ‘TAW’ and ‘site’ reported above were reliably estimated.
Discussion
Under low soil water availability, a largely conserved acclimatization of foliar N metabolism involves re-allocation from structural N to metabolic N
Our hypothesis that needle N profiles are affected by variation in environmental conditions at the field sites was supported by both, exploratory and inferential statistical analyses. Linear modelling revealed that total available soil water (TAW) was an important factor (Table 5) underlying variation in N pools between sites indicated by the exploratory analysis (Fig 1 and S1 Fig). The field experiment reported here was complemented by a recently published study of Junker et al. [13] on C metabolism, which revealed a pronounced reduction of photosynthetic CO2 assimilation rate, stomatal conductance and midday twig water potential under low TAW. These results show that responses to low TAW were not limited to N metabolism. In agreement with our hypothesis, we observed a trend to higher total amino acids (TAA) and arginine contents in the needles at lower TAW (Table 5 and S5B, S5D and S5F Fig). Likewise, other studies showed increased TAA contents as an indicator of water stress [69,70] and temperature stress [70]. Furthermore, increased arginine concentrations were reported in earlier studies with Douglas-fir seedlings exposed to controlled drought treatments [43]. Arginine is the proteinogenic amino acid with the lowest C/N ratio. Its accumulation under adverse environmental conditions restricting the N demand for growth might be a response to assimilate excess ammonium in a C-efficient way [26–28]. Moreover, arginine is a precursor for proline, nitric oxide and polyamines, which are compounds playing crucial roles in regulating responses to biotic and abiotic stress [71]. The observed increased arginine contents might thus indicate a physiological acclimatization of the trees to the pronounced reduction of photosynthetic CO2 assimilation rate and stomatal conductance under low TAW in our experiment [13].
In contrast to our initial hypothesis and to results for pine seedlings [44], total soluble proteins (TSP) did not decrease with decreasing TAW. Our data instead indicate that in Douglas-fir, decreasing contents of structural N rather than TSP serve as a source to maintain a high foliar amino acid content, as also observed in oak leaves [40].
Statistically significant although with respect to effect size small ‘provenance:TAW’ interaction effects were detected for TAA, arginine, and asparagine. In our initial hypothesis, we expected a smaller increase in amino acids in the interior provenance AR as a consequence of constitutively higher levels of N-rich amino acids such as arginine or asparagine. In contrast to this hypothesis, a smaller increase in arginine or asparagine under low TAW was not evident for the provenance AR. Moreover, ‘provenance:TAW’ interaction effects for arginine, asparagine and TAA were quite small with respect to effect size, and ‘provenance:TAW’ interaction effects were missing for glutamate and the other N pools (Table 5), suggesting that the adjustments of N metabolism in response to ‘TAW’ observed here is largely independent from constitutive differences in N pool sizes and composition. These findings illustrate that phenotypic plasticity of N metabolism in response to ‘TAW’ is largely conserved among the Douglas-fir provenances studied here.
Foliar N status of adult Douglas-fir trees is independent from site-related differences in soil N content
In our statistical analysis we used TAW to estimate soil water availability integrating important site-related factors including temperature, precipitation, latitude, slope and soil type. TAW was generally higher at site S compared to site W and thus was assumed to capture a large fraction of environmental variation between the two sites. In agreement with this assumption, the effect sizes of the generic factor ‘site’ were low for N pools and major amino acids (Table 5). It is intuitive to consider variation in soil N content as the underlying environmental factor accounting for effects of ‘site’ independent from effects of ‘TAW’. However, needle N pools were significantly higher at W compared to S (Fig 3, Table 5 and S5A–S5D Fig), despite significantly higher soil total N and dissolved mineral N contents at S compared to W (Table 3). This finding contrasts previous reports indicating higher leaf N contents in trees grown at sites of higher N availability or in N fertilization experiments [17–19,26,72,73]. Several possibilities might explain these unexpected observations:
- In forest soils, a high proportion of total N is incorporated in organic macromolecules such as proteins or chitin. Such compounds cannot be taken up directly by the roots of the trees and may lead to a discrepancy between total soil N content and available dissolved soil N content could, therefore, explain the differences between total N in the soil and N abundance in plant tissues [74]. However, both the higher natural 15N-abundance and higher levels of inorganic N at S compared to W, do not support this assumption, because they suggest a faster soil N-turnover at S [75] and, consequently, a better N supply with plant available N forms. Furthermore, the lower C / N ratios at S compared to W (Table 3) suggest reduced soil microbial competition for N against plants due to C limitation of heterotrophic soil microbes, thereby leaving more dissolved N for uptake by trees [76].
- The prevailing environmental conditions, such as the considerably cooler air temperatures at S, could have impaired root N uptake of the trees due to lower soil temperatures and a shorter growing season at S compared to W. This idea is supported by temperature dependencies of N uptake [77] and the observed slower growth of the experimental trees at S compared to W [11]. Moreover, it has been shown that elevated temperatures can cause higher needle N contents in Douglas-fir [78–80]. This is due to increased N uptake as a consequence of improved soil N availability [81,82], which most likely is the result of higher microbial activity in litter and soil organic matter [83,84]. Furthermore, the effect of temperature on root N uptake also depends on soil N status. For instance, N uptake by roots of trees was increased with increasing temperatures at low soil N availability (comparable to site W in our study), but was decreased in N rich soils [85].
- The different air temperatures at the two sites could have influenced the general N allocation patterns in the trees. Hu et al. [40] found enhanced biosynthesis of proteins and amino acids instead of incorporation into structural components in oak under drought and / or air warming. Effects of N availability on N allocation patterns were also reported by Gezelius and Näsholm [86]who observed high TAA concentrations in shoot, needles, and stems of Scots pine seedlings grown under low N supply, but higher concentrations of free amino acids in the roots of seedlings grown under higher N supply.
Consistent with findings on Scots pine [24,28,87] and Norway spruce [88], in the present study up to 19.8% of the total N in Douglas-fir needles was stored in free amino acids, predominantly as arginine, asparagine and glutamate. At W, arginine dominated the TAA profiles whereas this amino acid was virtually absent in Douglas-fir needles at S. Because arginine is also considered a N storage compound, our results suggest enhanced allocation into the N storage pool at W compared to S.
Variation of needle N pools among provenances is small compared to variation related to environmental conditions
Based on the assumption of local adaptation to N availability and environmental restrictions to carbon assimilation at the provenances’ native habitats, we hypothesized that Douglas-fir provenances differ in N pool sizes and -composition. Bond et al. [89] reported needle N contents of 12.5±1.3 mg g-1 for trees from a stand in the Willamette Valley (Oregon, USA). In our study, total foliar N content varied from 13.0 (AR at S) to 19.7 mg g-1 (CR at W). Significant provenance effects on N pools (except structural N), arginine, and glutamate revealed by our study suggest that signatures of local adaptation on N metabolism are detectable in the provenances studied here (Table 5). In contrast to our initial hypothesis, the interior provenance AR did not exhibit the highest levels of TSP, arginine or asparagine. Moreover, the proportion of variance explained by the factor ‘provenance’ is quite low compared to factors capturing environmental variation at the field sites, such as ‘TAW’ and ‘site’ (Table 5). This indicates that genetic differentiation in N metabolism is probably rather limited among the Douglas-fir provenances studied here. A high gene flow among Douglas-fir populations, especially among coastal provenances [89] might contribute to this finding.
Low genetic differentiation in N metabolism was also concluded from a common garden experiment with young Picea glauca trees from six seed sources, which did not differ with respect to their needle N content [34]. In contrast, when trees originating from the same six Picea glauca seed orchards were studied in a greenhouse experiment, small but significant variation in needle N content was detected, which was negatively correlated with mean growing season air temperature (MGST) at the seed orchard [90].
Raitio and Sarjala [25] reported that free amino acids in needles of Scots pine were more significantly affected by provenances as compared to experimental sites. Two common garden studies with Norway spruce [20] and Scots pine [24] reported increasing needle N concentrations with increasing altitude and decreasing mean annual temperature (MAT) [20] or with increasing latitude [24] of seed source origin. It was concluded that higher needle N concentrations were adaptive since they often contain large quantities of Rubisco [29,30] and are correlated with higher maximum photosynthetic capacity (Amax) [7,20,24,91]. A higher Amax, in turn, might at least partially compensate for environmental restrictions in carbon gain in habitats with low temperatures in higher altitudes or latitudes [20,24]. In our study, no obvious trend of total needle N with latitude, altitude, MAT or MGST was evident (Table 1 and Fig 3). However, we found a trend towards increasing TSP content with increasing MAT: provenance AR originating from the coolest habitat had lowest TSP content, provenances CR and RI had intermediate TSP content, and provenance LA, originating from the warmest native habitat, had highest TSP contents (Table 1 and Fig 3).
Does needle N content correlate with photosynthetic CO2 assimilation rates or growth?
To further assess whether intraspecific variation in needle total N was linked to intraspecific variation in adaptive traits related to carbon gain, we analyzed the correlation of total N content with net photosynthetic CO2 assimilation rate (A) collected at the same trees by Junker et al. [13]. For this analysis we selected data obtained in July at Schluchsee, when TAW was least restrictive for A. In contrast to results for Norway spruce, but in agreement with results for Picea glauca [34] no significant correlation between A and total N content was detected (r = -0.11, P = 0.683 for 2010; r = -0.09, P = 0.782 for 2011, S6 Fig). These findings do not support the concept that variation in needle N content among Douglas-fir provenances is related to adaptive variation in maximum C assimilation rates. It is known that in conifer needles Rubisco is not only present in a photosynthetically active form, but also a as storage form [92]. Along this line, Ripullone et al. [91] reported decreasing proportions of active Rubisco with increasing leaf N content in Douglas-fir. Thus, rather than being related to concentrations of Rubisco, provenance effects on CO2 assimilation rates might be a consequence of variation in CO2 diffusion related to stomatal and mesophyll conductance as also concluded from Benomar et al. [24] for field-grown Picea glauca.
A complementary study based on the same experiment revealed significant differences in height growth among the provenances analyzed here [11]. Provenance CR showed vigorous height growth across common garden sites, and this was associated with highest total N content (Fig 3). The interior provenance AR which exhibited lowest total N and TSP contents in our study (Fig 3), showed lowest height growth across common garden sites [11]. Both, total N content and TSP, were significantly correlated with tree height (S7 Fig). These findings are in agreement with previous research suggesting a role of needle N content for height growth [93].
Conclusions
This study reveals that across our study set of four genetically differentiated Douglas-fir provenances, needle N content is independent from soil N content. Small provenance effects suggest that genetic differentiation in N metabolism among the studied provenances is limited. Provenances displayed phenotypic plasticity to decreasing soil water availability by recruiting a largely conserved acclimatization of N metabolism involving a re-allocation of N from structural N to metabolic N. Small or missing interaction effects of ‘provenance’ with soil water availability or ‘site’ indicate that the variation in phenotypic plasticity among provenances is very limited. In the context of our earlier studies on the same experiment [13], our results further reveal that net photosynthetic CO2 assimilation rates do not increase with increasing needle N content. Further studies with a set of provenances spanning strong climatic gradients at low geographical and thus phylogenetic distance are needed to test how far these conclusions can be generalized.
Supporting information
S1 Table. Amino acids and ammonia content (μmol g-1 DW) in previous year needles of the four provenances (AR, CR, RI and LA) in May and July 2010 (5/10 and 7/10) and 2011 (5/11 and 7/11) at two sites (Wiesloch and Schluchsee).
https://doi.org/10.1371/journal.pone.0194684.s001
(PDF)
S2 Table. Results of sequential regression of N compounds in needles of adult Douglas-fir trees of four provenances on soil water availability and site effects independent from soil water availability.
https://doi.org/10.1371/journal.pone.0194684.s002
(PDF)
S1 Fig. Variable Importance in Projection (VIP) scores of component 1 from PLS_DA shown the importance of factors determining the site related (W, Wiesloch; S, Schluchsee) nitrogen partitioning patterns in 2010 (A) and 2011 (B).
https://doi.org/10.1371/journal.pone.0194684.s003
(PDF)
S2 Fig. Variable Importance in Projection (VIP) scores of component 1 from PLS_DA shown the importance determining the seasons related nitrogen partitioning patterns in May (5/10) and July (7/10) 2010 (A, B) and 2011(May, 5/11 and July, 7/11) (C, D) at Wiesloch (W) and Schluchsee (S) respectively.
https://doi.org/10.1371/journal.pone.0194684.s004
(PDF)
S3 Fig. Provenance-related patterns of nitrogen partitioning in previous year needles of the four Douglas-fir provenances sampled in 2010 (AR, Salmon Arm; CR, Conrad Creek; LA, Cameron Lake; RI, Santiam River).
Clustering was revealed by score plots of partial least squares–discriminant analyses.
https://doi.org/10.1371/journal.pone.0194684.s005
(PDF)
S4 Fig. Provenance-related patterns of nitrogen partitioning in previous year needles of the four Douglas-fir provenances sampled in 2011 (AR, Salmon Arm; CR, Conrad Creek; LA, Cameron Lake; RI, Santiam River).
Clustering was revealed by score plots of partial least squares–discriminant analyses.
https://doi.org/10.1371/journal.pone.0194684.s006
(PDF)
S5 Fig. Nitrogen (N) pools and compounds in needles of adult Douglas-fir trees of four provenances grown on two field sites in south-western Germany (Wiesloch and Schluchsee).
Data are plotted over experimental and environmental factors identified as significant predictors (for details of the statistical analysis see Material and methods; data are plotted on the scale used in statistical analysis).
https://doi.org/10.1371/journal.pone.0194684.s007
(PDF)
S6 Fig. Correlation between needle total N pools in needles of adult Douglas-fir trees of four provenances grown on field site ‘Schluchsee’ in south-western Germany with net photosynthetic assimilation rate (A) measured on the same trees in July 2010 and 2011.
https://doi.org/10.1371/journal.pone.0194684.s008
(PDF)
S7 Fig. Correlation between needle total nitrogen (N) pools in needles of adult Douglas-fir trees of four provenances with tree height.
Trees were grown on two field sites in south-western Germany.
https://doi.org/10.1371/journal.pone.0194684.s009
(PDF)
Acknowledgments
We thank Allison Kolar and Gustavo Saiz for their help during the analyses of soil samples at the Center of Stable Isotopes of the Institute of Meteorology and Climate Research Atmospheric Environmental Research (IMK-IFU), Karlsruhe Institute of Technology (KIT).
References
- 1.
Smith WB, Miles PD, Perry CH, Pugh SA. Forest resources of the United States, 2007: a technical document supporting the forest service 2010 RPA Assessment. Gen Tech Rep-USDA For Serv. 2009;Washington, DC: U.S. Department of Agriculture, Forest Service, Washington Office. 336 p.
- 2. Eckert AJ, Bower AD, Wegrzyn JL, Pande B, Jermstad KD, Krutovsky KV, et al. Association genetics of coastal Douglas fir (Pseudotsuga menziesii var. menziesii, Pinaceae). I. Cold-hardiness related traits. Genetics. 2009;182: 1289–1302. pmid:19487566
- 3. Savolainen O, Lascoux M, Merilä J. Ecological genomics of local adaptation. Nat Rev Genet. 2013;14: 807–820. pmid:24136507
- 4. St. Clair JB, Mandel NL, Vance-Borland KW. Genecology of Douglas fir in western Oregon and Washington. Ann Bot. 2005;96: 1199–1214. pmid:16246849
- 5. Brandl H. Ergänzende Untersuchungen zur Ertragslage der Baumarten Fichte, Kiefer, Buche und Eiche in Baden-Württemberg. Allg Forst Jagdztg. 1989;160: 91–98.
- 6. Hermann RK, Lavender DP. Douglas-fir planted forests. New For. 1999;17: 53–70.
- 7. Leites LP, Robinson AP, Rehfeldt GE, Marshall JD, Crookston NL. Height-growth response to climatic changes differs among populations of Douglas-fir: a novel analysis of historic data. Ecol Appl. 2012;22: 154–165. pmid:22471081
- 8. Eilmann B, de Vries SM, den Ouden J, Mohren GM, Sauren P, Sass-Klaassen U. Origin matters! Difference in drought tolerance and productivity of coastal Douglas-fir (Pseudotsuga menziesii (Mirb.)) provenances. For Ecol Manag. 2013;302: 133–143.
- 9. Jansen K, Sohrt J, Kohnle U, Ensminger I, Gessler A. Tree ring isotopic composition, radial increment and height growth reveal provenance-specific reactions of Douglas-fir towards environmental parameters. Trees. 2013;27: 37–52.
- 10. Sergent A-S, Bréda N, Sanchez L, Bastein J-C, Rozenberg P. Coastal and interior Douglas-fir provenances differ in growth performance and response to drought episodes at adult age. Ann For Sci. 2014;71: 709–720.
- 11. Neophytou C, Weisser A-M, Landwehr D, Šeho M, Kohnle U, Ensminger I, et al. Assessing the relationship between height growth and molecular genetic variation in Douglas-fir (Pseudotsuga menziesii) provenances. Eur J For Res. 2016;135: 465–481.
- 12. Gould PJ, Harrington CA, Clair JBS. Growth phenology of coast Douglas-fir seed sources planted in diverse environments. Tree Physiol. 2012;32: 1482–1496. pmid:23135739
- 13. Junker LV, Kleiber A, Jansen K, Wildhagen H, Hess M, Kayler Z, et al. Variation in short-term and long-term responses of photosynthesis and isoprenoid-mediated photoprotection to soil water availability in four Douglas-fir provenances. Sci Rep. 2017;7: 40145. pmid:28071755
- 14. Zhang JW, Marshall JD. Variation in carbon isotope discrimination and photosynthetic gas exchange among populations of Pseudotsuga menziesii and Pinus ponderosa in different environments. Funct Ecol. 1995; 402–412.
- 15. Reich PB, Walters MB, Kloeppel BD, Ellsworth DS. Different photosynthesis-nitrogen relations in deciduous hardwood and evergreen coniferous tree species. Oecologia. 1995;104: 24–30. pmid:28306909
- 16. Rennenberg H, Kreutzer K, Papen H, Weber P. Consequences of high loads of nitrogen for spruce (Picea abies) and beech (Fagus sylvatica) forests. New Phytol. 1998;139: 71–86.
- 17. Luo J, Li H, Liu T, Polle A, Peng C, Luo Z-B. Nitrogen metabolism of two contrasting poplar species during acclimation to limiting nitrogen availability. J Exp Bot. 2013;64: 4207–4224. pmid:23963674
- 18. Luo J, Zhou J, Li H, Shi W, Polle A, Lu M, et al. Global poplar root and leaf transcriptomes reveal links between growth and stress responses under nitrogen starvation and excess. Tree Physiol. 2015;35: 1283–1302. pmid:26420789
- 19. Li H, Li M, Luo J, Cao X, Qu L, Gai Y, et al. N-fertilization has different effects on the growth, carbon and nitrogen physiology, and wood properties of slow-and fast-growing Populus species. J Exp Bot. 2012;63: 6173–6185. pmid:23028021
- 20. Oleksyn J, Modrzýnski J, Tjoelker MG, Reich PB, Karolewski P. Growth and physiology of Picea abies populations from elevational transects: common garden evidence for altitudinal ecotypes and cold adaptation. Funct Ecol. 1998;12: 573–590.
- 21. Kranabetter JM. Ectomycorrhizal fungi and the nitrogen economy of conifers—implications for genecology and climate change mitigation. Botany. 2014;92: 417–423.
- 22. Shay P-E, Winder RS, Trofymow JA. Nutrient-cycling microbes in coastal Douglas-fir forests: regional-scale correlation between communities, in situ climate, and other factors. Front Microbiol. 2015;6: 1097. pmid:26500636
- 23. Gessler A, Schaub M, McDowell NG. The role of nutrients in drought-induced tree mortality and recovery. New Phytol. 2017;214: 513–520. pmid:27891619
- 24. Oleksyn J, Reich PB, Zytkowiak R, Karolewski P, Tjoelker MG. Nutrient conservation increases with latitude of origin in European Pinus sylvestris populations. Oecologia. 2003;136: 220–235. pmid:12756524
- 25. Raitio H, Sarjala T. Effect of provenance on free amino acid and chemical composition of Scots pine needles. Plant Soil. 2000;221: 231–238.
- 26. Näsholm T, Ericsson A. Seasonal changes in amino acids, protein and total nitrogen in needles of fertilized Scots pine trees. Tree Physiol. 1990;6: 267–281. pmid:14972938
- 27. Rabe E, Lovatt CJ. Increased arginine biosynthesis during phosphorus deficiency a response to the increased ammonia content of leaves. Plant Physiol. 1986;81: 774–779. pmid:16664901
- 28. Dijk HF, Roelofs JG. Effects of excessive ammonium deposition on the nutritional status and condition of pine needles. Physiol Plant. 1988;73: 494–501.
- 29. Field C. Allocating leaf nitrogen for the maximization of carbon gain: leaf age as a control on the allocation program. Oecologia. 1983;56: 341–347. pmid:28310214
- 30.
Field C, Mooney HA. Adaptative patterns of energy capture in plants. Givnish TJ, editor in: On the Economy of Plant Form and Function: Proceedings of the Sixth Maria Moors Cabot Symposium. New York, USA: Cambrigde University Press; 1986. pp. 25–55.
- 31. Alberto FJ, Aitken SN, Alía R, González-Martínez SC, Hänninen H, Kremer A, et al. Potential for evolutionary responses to climate change—evidence from tree populations. Glob Change Biol. 2013;19: 1645–1661.
- 32. Anderson JT, Willis JH, Mitchell-Olds T. Evolutionary genetics of plant adaptation. Trends Genet. 2011;27: 258–266. pmid:21550682
- 33. Nicotra AB, Atkin OK, Bonser SP, Davidson AM, Finnegan E, Mathesius U, et al. Plant phenotypic plasticity in a changing climate. Trends Plant Sci. 2010;15: 684–692. pmid:20970368
- 34. Benomar L, Lamhamedi MS, Rainville A, Beaulieu J, Bousquet J, Margolis HA. Genetic Adaptation vs. Ecophysiological Plasticity of Photosynthetic-Related Traits in Young Picea glauca Trees along a Regional Climatic Gradient. Front Plant Sci. 2016;7: 48. pmid:26870067
- 35. Du B, Jansen K, Kleiber A, Eiblmeier M, Kammerer B, Ensminger I, et al. A coastal and an interior Douglas fir provenance exhibit different metabolic strategies to deal with drought stress. Tree Physiol. 2016;36: 148–163. pmid:26491053
- 36. Schneider S, Gessler A, Weber P, Sengbusch D, Hanemann U, Rennenberg H. Soluble N compounds in trees exposed to high loads of N: a comparison of spruce (Picea abies) and beech (Fagus sylvatica) grown under field conditions. New Phytol. 1996;134: 103–114.
- 37. Rennenberg H, Loreto F, Polle A, Brilli F, Fares S, Beniwal RS, et al. Physiological responses of forest trees to heat and drought. Plant Biol. 2006;8: 556–571. pmid:16773557
- 38. Gessler A, Keitel C, Nahm M, Rennenberg H. Water shortage affects the water and nitrogen balance in central European beech forests. Plant Biol. 2004;6: 289–298. pmid:15143437
- 39. Fotelli MN, Nahm M, Heidenfelder A, Papen H, Rennenberg H, Geßler A. Soluble nonprotein nitrogen compounds indicate changes in the nitrogen status of beech seedlings due to climate and thinning. New Phytol. 2002;154: 85–97.
- 40. Hu B, Simon J, Kuster TM, Arend M, Siegwolf R, Rennenberg H. Nitrogen partitioning in oak leaves depends on species, provenance, climate conditions and soil type. Plant Biol. 2013;15: 198–209. pmid:22934888
- 41. Guo J, Yang Y, Wang G, Yang L, Sun X. Ecophysiological responses of Abies fabri seedlings to drought stress and nitrogen supply. Physiol Plant. 2010;139: 335–347. pmid:20230480
- 42. De Diego N, Saiz-Fernández I, Rodríguez JL, Pérez-Alfocea P, Sampedro MC, Barrio RJ, et al. Metabolites and hormones are involved in the intraspecific variability of drought hardening in radiata pine. J Plant Physiol. 2015;188: 64–71. pmid:26433462
- 43. Driessche RVD, Webber JE. Total and soluble nitrogen in Douglas fir in relation to plant nitrogen status. Can J For Res. 1975;5: 580–585.
- 44. Wellburn FA, Lau K-K, Milling PM, Wellburn PM. Drought and air pollution affect nitrogen cycling and free radical scavenging in Pinus halepensis (Mill.). J Exp Bot. 1996;47: 1361–1367.
- 45. Good AG, Zaplachinski ST. The effects of drought stress on free amino acid accumulation and protein synthesis in Brassica napus. Physiol Plant. 1994;90: 9–14.
- 46. Kenk G, Thren M. Ergebnisse verschiedener Douglasienprovenienzversuche in Baden-Wurttemberg. I. Der Internationale Douglasien-Provenienzversuch 1958. Allg Forst Jagdztg. 1984;155: 165–184.
- 47. Müller T, Freund F, Wildhagen H, Schmid KJ. Targeted re-sequencing of five Douglas-fir provenances reveals population structure and putative target genes of positive selection. Tree Genet Genomes. 2015;11: 1–17.
- 48. Wang T, Hamann A, Spittlehouse DL, Murdock TQ. ClimateWNA—high-resolution spatial climate data for western North America. J Appl Meteorol Climatol. 2012;51: 16–29.
- 49. Wang T, Hamann A, Spittlehouse D, Carroll C. Locally downscaled and spatially customizable climate data for historical and future periods for North America. PLoS One. 2016;11: e0156720. pmid:27275583
- 50. Kohnle U, Hein S, Sorensen FC, Weiskittel AR. Effects of seed source origin on bark thickness of Douglas-fir (Pseudotsuga menziesii) growing in southwestern Germany. Can J For Res. 2012;42: 382–399.
- 51. Hamann A, Wang T, Spittlehouse DL, Murdock TQ. A comprehensive, high-resolution database of historical and projected climate surfaces for western North America. Bull Am Meteorol Soc. 2013;94: 1307–1309.
- 52. Russell DJ, Alberti G. Effects of long-term, geogenic heavy metal contamination on soil organic matter and microarthropod communities, in particular Collembola. Appl Soil Ecol. 1998;9: 483–488.
- 53. Feger KH. Nitrogen cycling in two Norway spruce (Picea abies) ecosystems and effects of a (NH4)2SO4 addition. Water Air Soil Pollut. 1992;61: 295–307.
- 54. Keitel C, Matzarakis A, Rennenberg H, Gessler A. Carbon isotopic composition and oxygen isotopic enrichment in phloem and total leaf organic matter of European beech (Fagus sylvatica L.) along a climate gradient. Plant Cell Environ. 2006;29: 1492–1507. pmid:16898013
- 55. Liu M, Dannenmann M, Lin S, Saiz G, Yan G, Yao Z, et al. Ground cover rice production systems increase soil carbon and nitrogen stocks at regional scale. Biogeosciences. 2015;12: 4831–4840.
- 56. Schlotter D, Schack-Kirchner H, Hildebrand EE, von Wilpert K. Equivalence or complementarity of soil-solution extraction methods. J Plant Nutr Soil Sci. 2012;175: 236–244.
- 57. Baethgen WE, Alley MM. A manual colorimetric procedure for measuring ammonium nitrogen in soil and plant Kjeldahl digests. Commun Soil Sci Plant Anal. 1989;20: 961–969.
- 58. Du B, Jansen K, Junker LV, Eiblmeier M, Kreuzwieser J, Gessler A, et al. Elevated temperature differently affects foliar nitrogen partitioning in seedlings of diverse Douglas fir provenances. Tree Physiol. 2014;34: 1090–1101. pmid:25240727
- 59. Dannenmann M, Simon J, Gasche R, Holst J, Naumann PS, Kögel-Knabner I, et al. Tree girdling provides insight on the role of labile carbon in nitrogen partitioning between soil microorganisms and adult European beech. Soil Biol Biochem. 2009;41: 1622–1631.
- 60. Winter H, Lohaus G, Heldt HW. Phloem transport of amino acids in relation to their cytosolic levels in barley leaves. Plant Physiol. 1992;99: 996–1004. pmid:16669030
- 61. Luo Z-B, Janz D, Jiang X, Göbel C, Wildhagen H, Tan Y, et al. Upgrading root physiology for stress tolerance by ectomycorrhizas: insights from metabolite and transcriptional profiling into reprogramming for stress anticipation. Plant Physiol. 2009;151: 1902–1917. pmid:19812185
- 62. Xia J, Sinelnikov IV, Han B, Wishart DS. MetaboAnalyst 3.0—making metabolomics more meaningful. Nucleic Acids Res. 2015;43: W251–W257. pmid:25897128
- 63.
R Core Team. R: A Language and Environment for Statistical Computing [Internet]. Vienna, Austria: R Foundation for Statistical Computing; 2016. https://www.R-project.org/
- 64.
Zuur A, Ieno EN, Walker N, Saveliev AA, Smith GM. Mixed Effects Models and Extensions in Ecology with R. Springer-Verlag New York; 2009.
- 65. Hess M, Wildhagen H, Junker LV, Ensminger I. Transcriptome responses to temperature, water availability and photoperiod are conserved among mature trees of two divergent Douglas-fir provenances from a coastal and an interior habitat. BMC Genomics. 2016;17: 682. https://doi.org/10.1186/s12864-016-3022-6 pmid:27565139
- 66.
Pinheiro J, Bates D, DebRoy S, Sarkar D. nlme: Linear and Nonlinear Mixed Effects Models. In: https://cran.r-project.org/web/packages/nlme/index.html. 2017.
- 67.
Nakazawa M. fmsb: Functions for Medical Statistics Book with some Demographic Data [Internet]. 2017. https://CRAN.R-project.org/package=fmsb
- 68. Dormann CF, Elith J, Bacher S, Buchmann C, Carl G, Carré G, et al. Collinearity: a review of methods to deal with it and a simulation study evaluating their performance. Ecography. 2013;36: 27–46.
- 69. Rai VK. Role of amino acids in plant responses to stresses. Biol Plant. 2002;45: 481–487.
- 70.
Stewart GR, Larher F, Miflin BJ. Accumulation of amino acids and related compounds in relation to environmental stress. Biochem Plants Compr Treatise Vol 5 Amino Acids Deriv. 1980; 609–635.
- 71. Winter G, Todd CD, Trovato M, Forlani G, Funck D. Physiological implications of arginine metabolism in plants. Front Plant Sci. 2015;6: 534. pmid:26284079
- 72. Margolis HA, Waring RH. Carbon and nitrogen allocation patterns of Douglas-fir seedlings fertilized with nitrogen in autumn. I. Overwinter metabolism. Can J For Res. 1986;16: 897–902.
- 73. Turner J, Olson PR. Nitrogen relations in a Douglas-fir plantation. Ann Bot. 1976;40: 1185–1193.
- 74. Rennenberg H, Dannenmann M. Nitrogen nutrition of trees in temperate forests—the significance of nitrogen availability in the pedosphere and atmosphere. Forests. 2015;6: 2820–2835.
- 75. Hogberg P. 15N natural abundance in soil–plant systems. Tansley Review No. 95. New Phytol. 1997;137: 179–203.
- 76. Dannenmann M, Gasche R, Ledebuhr A, Papen H. Effects of forest management on soil N cycling in beech forests stocking on calcareous soils. Plant Soil. 2006;287: 279–300.
- 77. Geßler A, Jung K, Gasche R, Papen H, Heidenfelder A, Börner E, et al. Climate and forest management influence nitrogen balance of European beech forests: microbial N transformations and inorganic N net uptake capacity of mycorrhizal roots. Eur J For Res. 2005;124: 95–111.
- 78. Lewis JD, Lucash M, Olszyk DM, Tingey DT. Relationships between needle nitrogen concentration and photosynthetic responses of Douglas-fir seedlings to elevated CO2 and temperature. New Phytol. 2004;162: 355–364.
- 79. Hobbie EA, Olszyk DM, Rygiewicz PT, Tingey DT, Johnson MG. Foliar nitrogen concentrations and natural abundance of 15N suggest nitrogen allocation patterns of Douglas-fir and mycorrhizal fungi during development in elevated carbon dioxide concentration and temperature. Tree Physiol. 2001;21: 1113–1122. pmid:11581018
- 80. Tingey DT, Mckane RB, Olszyk DM, Johnson MG, Rygiewicz PT, Henry Lee E. Elevated CO2 and temperature alter nitrogen allocation in Douglas-fir. Glob Change Biol. 2003;9: 1038–1050.
- 81. Cleve KV, Oechel WC, Hom JL. Response of black spruce (Picea mariana) ecosystems to soil temperature modification in interior Alaska. Can J For Res. 1990;20: 1530–1535.
- 82. van Breemen N, Jenkins A, Wright RF, Beerling DJ, Arp WJ, Berendse F, et al. Impacts of elevated carbon dioxide and temperature on a boreal forest ecosystem (CLIMEX project). Ecosystems. 1998;1: 345–351.
- 83. Lin G, Ehleringer JR, Rygiewicz P, Johnson MG, Tingey DT. Elevated CO2 and temperature impacts on different components of soil CO2 efflux in Douglas-fir terracosms. Glob Change Biol. 1999;5: 157–168.
- 84. Lin G, Rygiewicz PT, Ehleringer JR, Johnson MG, Tingey DT. Time-dependent responses of soil CO2 efflux components to elevated atmospheric [CO2] and temperature in experimental forest mesocosms. Plant Soil. 2001;229: 259–270.
- 85. Bassirirad H. Kinetics of nutrient uptake by roots: responses to global change. New Phytol. 2000;147: 155–169.
- 86. Gezelius K, Näsholm T. Free amino acids and protein in Scots pine seedlings cultivated at different nutrient availabilities. Tree Physiol. 1993;13: 71–86. pmid:14969902
- 87. Pérez-Soba M, de Visser PH. Nitrogen metabolism of Douglas fir and Scots pine as affected by optimal nutrition and water supply under conditions of relatively high atmospheric nitrogen deposition. Trees. 1994;9: 19–25.
- 88. Weber P, Stoermer H, Gessler A, Schneider S, Von Sengbusch D, Hanemann U, et al. Metabolic responses of Norway spruce (Picea abies) trees to long-term forest management practices and acute (NH4)2SO4 fertilization: transport of soluble non-protein nitrogen compounds in xylem and phloem. New Phytol. 1998;140: 461–475.
- 89. Krutovsky KV, St. Clair JB, Saich R, Hipkins VD, Neale DB. Estimation of population structure in coastal Douglas-fir [Pseudotsuga menziesii (Mirb.) Franco var. menziesii] using allozyme and microsatellite markers. Tree Genet Genomes. 2009;5: 641–658.
- 90. Benomar L, Lamhamedi MS, Villeneuve I, Rainville A, Beaulieu J, Bousquet J, et al. Fine-scale geographic variation in photosynthetic-related traits of Picea glauca seedlings indicates local adaptation to climate. Tree Physiol. 2015;35: 864–878. pmid:26116923
- 91. Ripullone F, Grassi G, Lauteri M, Borghetti M. Photosynthesis–nitrogen relationships: interpretation of different patterns between Pseudotsuga menziesii and Populus × euroamericana in a mini-stand experiment. Tree Physiol. 2003;23: 137–144. pmid:12533308
- 92. Millard P, Grelet G. Nitrogen storage and remobilization by trees: ecophysiological relevance in a changing world. Tree Physiol. 2010;30: 1083–1095. pmid:20551251
- 93. Walters MB, Gerlach JP. Intraspecific growth and functional leaf trait responses to natural soil resource gradients for conifer species with contrasting leaf habit. Tree Physiol. 2013;33: 297–310. pmid:23370548