Original Research ARTICLE
Genetic Adaptation vs. Ecophysiological Plasticity of Photosynthetic-Related Traits in Young Picea glauca Trees along a Regional Climatic Gradient
- 1Faculté de Foresterie, de Géographie et de Géomatique, Centre D'étude de la Forêt, Université Laval, Quebec, QC, Canada
- 2Direction de la Recherche Forestière, Ministère des Forêts, de la Faune et des Parcs, Quebec, QC, Canada
Assisted population migration (APM) is the intentional movement of populations within a species range to sites where future environmental conditions are projected to be more conducive to growth. APM has been proposed as a proactive adaptation strategy to maintain forest productivity and to reduce the vulnerability of forest ecosystems to projected climate change. The validity of such a strategy will depend on the adaptation capacity of populations, which can partially be evaluated by the ecophysiological response of different genetic sources along a climatic gradient. This adaptation capacity results from the compromise between (i) the degree of genetic adaptation of seed sources to their environment of origin and (ii) the phenotypic plasticity of functional trait which can make it possible for transferred seed sources to positively respond to new growing conditions. We examined phenotypic variation in morphophysiological traits of six seed sources of white spruce (Picea glauca [Moench] Voss) along a regional climatic gradient in Québec, Canada. Seedlings from the seed sources were planted at three forest sites representing a mean annual temperature (MAT) gradient of 2.2°C. During the second growing season, we measured height growth (H2014) and traits related to resources use efficiency and photosynthetic rate (Amax). All functional traits showed an adaptive response to the climatic gradient. Traits such as H2014, Amax, stomatal conductance (gs), the ratio of mesophyll to stomatal conductance, water use efficiency, and photosynthetic nitrogen-use efficiency showed significant variation in both physiological plasticity due to the planting site and seed source variation related to local genetic adaptation. However, the amplitude of seed source variation was much less than that related to plantation sites in the area investigated. The six seed sources showed a similar level of physiological plasticity. H2014, Amax and gs, but not carboxylation capacity (Vcmax), were correlated and decreased with a reduction of the average temperature of the growing season at seed origin. The clinal variation in H2014 and Amax appeared to be driven by CO2 conductance. The presence of locally adapted functional traits suggests that the use of APM may have advantages for optimizing seed source productivity in future local climates.
In boreal regions, climate is the main factor controlling the spatial distribution of tree species (Davis and Shaw, 2001). For most widespread tree species, natural selection processes operated over thousands of years, has led to local adaptation and differentiation of populations along climatic gradients across their natural range (e.g., Andalo et al., 2005; Aitken and Whitlock, 2013; Savolainen et al., 2013). This climate-related adaptive differentiation along geographic gradients has been demonstrated through provenance tests for several boreal tree species in eastern Canada (Li et al., 1997; Beaulieu et al., 2004; Andalo et al., 2005; Lu et al., 2014; Rossi and Bousquet, 2014; Yang et al., 2015). Therefore, the actual pace of climate change is expected to induce a spatial mismatch between locally adapted populations and their optimal climatic conditions to which they have historically adapted, leading potentially to local maladaptation (Gray and Hamann, 2011; Aitken and Whitlock, 2013). Such a trend could result in a reduction in forest health and productivity, as well as in ecosystem services and the resources they provide (Andalo et al., 2005; Aitken et al., 2008; O'Neill et al., 2008). Given that the evidence for global climate change is now unequivocal (Aitken and Bemmels, 2015), the paradigm of “local seed is the best” may no longer be appropriate for forest plantation programs under a changing climate (Gray and Hamann, 2013; Rainville et al., 2014; Aitken and Bemmels, 2015).
Assisted population migration (APM) is the intentional movement of populations to areas that are expected to harbor future environmental conditions similar to those at their geographical location of origin (Aitken and Whitlock, 2013; Ste-Marie, 2014). APM has been proposed (with current implementation planning in Canada) as a proactive adaptation strategy to maintain forest productivity and reduce the vulnerability of forest ecosystems to climate change (Gray and Hamann, 2011; Pedlar et al., 2011; Aitken and Whitlock, 2013). Despite the potential benefits of this proactive approach, APM implementation should nevertheless be preferentially preceded by trial tests supplemented by the acquisition of more exhaustive ecophysiological knowledge (Isaac-Renton et al., 2014). This is particularly important at the juvenile phase for appropriate adaptation and establishment of seedlings under different site conditions.
Climate-based seed transfer models based on the relationship between climatic parameters such as mean annual temperature (MAT) and summary quantitative traits such as height growth have been widely used to predict optimal climate-transfer distances for APM (Andalo et al., 2005; Wang et al., 2010; O'Neill and Nigh, 2011). These empirical models were generally based on previous provenance test results and transfer distances varying from about 1.5–4°C (MAT) in the boreal region (Andalo et al., 2005; Wang et al., 2010; Pedlar et al., 2011; Yang et al., 2015). However, these models did not take into account several considerations which may increase the risks associated with seed source transfer (O'Neill et al., 2008; Isaac-Renton et al., 2014). These include (i) the higher temporal variability of climatic conditions in boreal regions resulting from climate change, which may affect populations fitness depending on whether or not an extreme event occurs during provenance testing (particularly at the juvenile growing stage); (ii) the uncertainties in regional predictions of the magnitude and the direction of climate changes in boreal regions (Roy et al., 2014; Aitken and Bemmels, 2015); (iii) the predicted increase in the occurrence of extreme events and temporal variability in temperature and rainfall under climate change (Bourque and Simonet, 2007); (iv) the possible contribution of non-climatic factors (nutrient availability, soil properties, and mycorrhizal communities) to local adaptation of tree populations (Aitken and Whitlock, 2013; Kranabetter et al., 2015; Pickles et al., 2015); and (v) the possibility that phenotypic plasticity may vary among populations, with populations exhibiting higher levels of plasticity showing likely a greater survival rate but a lower growth rate.
The assessment of phenotypic variation and adaptive value of functional traits involved in carbon assimilation and resource use efficiencies among seed sources in response to variation in climatic conditions should contribute toward more informed ecophysiological-based methods of matching populations to site conditions under climate change (Gray and Hamann, 2013; Isaac-Renton et al., 2014; Prober et al., 2015). The assessment of such traits should enable to improve APM beyond summary parameters such as growth and yield estimates and thereby, help minimize the risk of establishment failure, and maximize the adaptation potential of transferred seed sources.
Variation in functional traits along a climatic gradient is largely influenced by local adaptation and phenotypic plasticity, i.e., the ability of a genotype to produce different phenotypes under different environmental conditions (Valladares et al., 2014; Anderson and Gezon, 2015). Thus, disentangling genetic (local adaptation) and environmental (phenotypic plasticity) responses could improve our understanding and prediction of the physiological and growth responses of populations to climate change (McLean et al., 2014; Merilä and Hendry, 2014; Anderson and Gezon, 2015) and determine the ability of populations to adjust to novel environmental conditions during and after assisted migration (Isaac-Renton et al., 2014; Anderson and Gezon, 2015).
Photosynthetic-related traits are key adaptive and fitness-related functional traits controlling carbon uptake, resources use efficiencies, and relative growth rates. Many studies conducted on different species have recorded among populations phenotypic differentiation for these traits. For instance, leaf nitrogen concentration (Nmass) and respiration rate (Rd) have been found to increase with latitude and elevation while specific leaf area (SLA) decreases (Friend et al., 1989; Oleksyn et al., 1998; Hultine and Marshall, 2000; Qiuhong et al., 2013). Photosynthetic rate (Amax) has also been reported to increase with latitude and elevation (Cordell et al., 1998; Oleksyn et al., 1998; Soolanayakanahally et al., 2009), but has been also reported to decrease or to have no linear relationships with latitude or elevation (Kogami et al., 2001; Dang et al., 2008; Qiuhong et al., 2013; Benomar et al., 2015). Amax is affected by biochemical and biophysical processes and environmental conditions (Sharkey et al., 2007). Indeed, the trend of Amax along latitude or elevation may be species-specific and may result from the relative limitation imposed by each process.
Carbon isotope discrimination occurring during C3 photosynthesis (Δ13C; i.e., the ratio of the amount of 13C to 12C isotopes in a sample relative to a standard) is related to the supply and demand of CO2 within the needle. Thereby, Δ13C is intrinsically determined by biochemical (carboxylation and photorespiration) and biophysical processes (Farquhar et al., 1989; Cernusak et al., 2013). Leaf Δ13C, which scales negatively with WUE (Farquhar et al., 1989; Warren and Adams, 2006), was found to vary with stomatal and mesophyll conductance (Kogami et al., 2001; Warren and Adams, 2006; Flexas et al., 2008), leaf morphology (Körner et al., 1991), leaf age (Ethier et al., 2006), leaf nitrogen concentration (Sparks and Ehleringer, 1997), and as well as with latitude and elevation (Körner et al., 1991; Cernusak et al., 2013). Therefore, Δ13C may represent an integrated measure of water and carbon economics mediated by the modulation of biophysical and biochemical processes in response to environmental conditions. Little research has been conducted on genetic variation and phenotypic plasticity in Δ13C for white spruce.
White spruce (Picea glauca [Moench] Voss) is one of the most commercially important tree species in the boreal forest of North America (Beaulieu et al., 2009). Its natural stands of the commercial forest zone in Québec are located between 45° and 50° N latitude and between 57° and 79°W longitude (de Lafontaine et al., 2010). Evidence is accumulating that white spruce populations within this narrow geographical range are differentiated with latitude, and more weakly with longitude (Khalil, 1985; Li et al., 1997; Jaramillo-Correa et al., 2001; Namroud et al., 2008; Rainville et al., 2014). Our recent results (Benomar et al., 2015; Villeneuve, 2015), under controlled conditions, showed the presence of clinal variation in photosynthetic-related traits. The relative contribution of both phenotypic plasticity and genetic variation to this phenotypic divergence has fundamental implications to improve our understanding and prediction of the physiological and growth responses of white spruce populations to climate change. However, variation in functional trait plasticity among white spruce seed sources remains poorly studied. To our knowledge, this paper is the first to investigate the possible occurrence of clinal patterns of phenotypic plasticity of physiological traits in white spruce. Consequently, the aim of the present study was to disentangle the contributions of local genetic adaptation and physiological plasticity to the observed geographical variation in photosynthetic-related traits among white spruce seed sources during the establishment phase after outplanting under different sites conditions along a regional climatic gradient. We address the following questions: (1) Do photosynthetic-related traits variation in young white spruce reflect local adaptation at least in part? (2) Does plasticity in photosynthetic-related traits in response to climate transfer vary among seed sources? If so, (3) does the degree of plasticity correlate with the climate of seed origin?
Materials and Methods
Genetic Material and Seedling Production
The white spruce seed sources used in this study came from six first-generation local seed orchards commonly used for reforestation in Québec, Canada (Table 1). The first-generation seed orchards were established about 30 years ago separately for each region using phenotypically selected plus-trees from local natural forests (Figure 1). Open-pollinated seed was collected in each seed orchard for two consecutive years (2008 and 2009). After mixing the collection years for each individual seed orchard, seedlings were produced in a commercial forest nursery at St-Modeste Québec, Canada (47.50°N, 69.23°W) under standard nursery cultural practices in Quebec (Lamhamedi et al., 2006). Seeds having similar seed size among seed sources (seed orchards) were stratified and were manually sown in May 2011 into containers filled with peat:vermiculite (3/1, v/v; bulk density of 0.10 g cm−3). Seedlings were cultivated under an unheated polyethylene tunnel during their first growing season and outdoors during their second growing season. Details relating to usual irrigation and fertilization practices during seedlings production in Québec nursery are described in Lamhamedi et al. (2006).
Table 1. Location and climatic conditionsa (over 1981–2010) of the six white spruce seed sources used in this study and the three plantation sites.
Plantation Sites and Experimental Design
The study area was located in the eastern Canadian boreal forest. Three forest sites were selected for this study and the plantations were established in the spring of 2013: Watford, Asselin, and Deville (Table 1, Figure 1). These sites represent the south, center, and north of the white spruce commercial zone in the eastern part of the boreal forest of Québec, Canada (Figure 1). The Watford (southern) site was formerly occupied by a plantation of black spruce (Picea mariana B.S.P. Voss) which was harvested in 2012. It is located in the sugar maple-yellow birch domain (near the locality of Ste-Rose de Watford, Québec, Canada) on a loamy soil. The Asselin (central) site was a former mixed forest stand growing on a loamy soil, which was commercially harvested in 2011. It is located in the balsam fir-yellow birch domain. The Deville (northern) site with clay loamy soil was also dominated by a mixed forest and was commercially harvested in 2012. It is located in the balsam fir-white birch domain. Soil chemical properties before planting at the three sites are given in Table S1.
Two-year-old seedlings were planted at densities of 2000 stems/ha during the last week of May 2013 at Watford, the first week of June at Asselin, and the second week of June at Deville. Before plantation, 15 seedlings per seed source (SO) were used for a quality check, based upon 27 morphological and nutritional standards developed for large containerized white spruce seedlings in Québec (Table S2). Seedling height and root collar diameter (mean ± standard deviation) were 38.3 ± 3.6 cm and 6.9 ± 0.9 mm, respectively. During the first (2013) and the second (2014) growing seasons, weed competition was controlled mechanically at Watford and Asselin, while no treatment was necessary at Deville due to the very limited presence of weedy vegetation.
At each site, we chose a randomized complete block design with four blocks each one was partitioned into six plots within which the six seed sources (SO) were assigned randomly. The size of each plot was about 730 m2 and contained 144 trees (12 × 12 rows of trees) in which the 64 interior trees were considered, leaving 4 × 4 rows of border trees as a buffer zone.
Climate Data during the Sampling Year
Each site was equipped with a meteorological station. The station was supplied with a shielded air temperature and relative humidity sensor (HMP50, Vaisala, Helsinki, Finland), PAR sensors (Li 190 Campbell Scientific, Logan, UT, USA), pluviometer, and three soil temperature sensors placed at three soil depth (10, 20, and 30 cm). Data were continuously recorded hourly throughout the growing season of 2013 and 2014 (from 152 to 273 day of year) using datalogger (CR10X, Campbell Scientific, Logan, UT, USA) (Figure 2).
Figure 2. Mean daily air temperature, mean daily soil temperature at a depth of 10 cm, and mean daily air vapor pressure deficit (VPD) at the three plantation sites during the growing season of 2013 and 2014. DOY, day of year.
Height and survival were measured at the end of the second growing season (mid-October 2014) on the 64 central trees in each plot. A total of 4608 seedlings were measured (64 trees * 4 blocks * 6 seed sources * 3 sites). Tree height was measured by a standard 1 m-ruler with 0.5 cm precision.
Gas Exchange Measurements
Needle-level gas exchange measurements were made with a portable open-path gas-exchange system (Li-6400, Li-Cor Inc., Lincoln NE), equipped with a Lighted Conifer Chamber (model 6400-22L, Li-Cor Inc., Lincoln NE). Measurements were carried out on one randomly selected plant per plot (6 * 4 * 3 = 72 plants). Dark respiration (Rd) and photosynthetic carbon dioxide response curves (A–Ci) were measured on 1-year-old needles within the uppermost lateral shoot. Dark respiration (net CO2 exchange at PAR = 0 μmol·m−2·s−1) was measured before A–Ci curve to avoid any influence of high level of CO2 on respiration. Rd was recorded after at least 12 min of dark inside the cuvette by turning off LED light source. The A–Ci curve measurements were taken after 20 min of steady state conditions at an ambient atmospheric CO2 partial pressure (Ca = 400 μmol mol−1) and at saturated photosynthetic active radiation (PAR = 1000 μmol m−2 s−1). Thereafter, the reference CO2 (Ca) was changed in the following order: Ca: 400, 300, 200, 100, 50, 400, 500, 600, 800, 1000, 1200, 1400, and 1500 μmol mol−1. Values were recorded based on the stability of photosynthesis, stomatal conductance, CO2, and water vapor concentrations. This was generally achieved after 5–7 min at each step. The attachment points of the shoot to cuvette walls were taped (adhesive tape) allowing for, at most, insignificant leaks into and out of the cuvette. Data were not corrected for a diffusion of CO2. During measurement, the leaf chamber conditions were: T = 25°C, RH = 60 ± 5, VPD = 1.2 ± 0.2 kPa, flow = 300 μmol s−1, and PAR = 1000 μmol m−2 s−1. A total of 72 A–Ci curves were measured. The VPD was easily maintained between 1.1 and 1.4 kPa for days with temperature lower than 25°C and problematic to maintain in days with temperature superior to 25°C. Measurements were made from 08h00 to 12h00 and from 14h00 to 17h00 when VPD within the cuvette still within a suitable range (1.0–1.4 kPa).
Estimation of Vcmax, Jmax, and gm by Curve-Fitting Method
Photosynthetic parameters such as maximum rates of carboxylation (Vcmax), maximum rate of electron transport (Jmax), and mesophyll conductance (gm: CO2 transfer conductance from the intercellular space to the site of carboxylation) were derived from the A–Ci response curves. Following Ethier et al. (2006) and Miao et al. (2009), A–Ci curves were fit with the biochemical model of C3 photosynthesis developed by Farquhar et al. (1980) and modified by Sharkey et al. (2007) using non-linear regression techniques (Proc NLIN, SAS) as described below. The biochemical model of C3 leaf photosynthesis assumes that the rate of CO2 uptake (An) is limited by: (i) the Rubisco activity (Ac), (ii) the rate of regeneration of the RuBP (Aj) and (iii) the rate of use of the triose-P (Ap). The final rate, An, is:
In the present study, Ap was not considered since it has been rarely observed in vivo (von Caemmerer, 2000). Rl, respiration occurring in daylight which is assumed to be primarily mitochondrial respiration (von Caemmerer, 2000) and was assumed to approximate dark respiration (Rd).
The Rubisco-limited rate of CO2 assimilation (Ac) is given by:
where Vcmax is the maximum rate of carboxylation; Cc, the chloroplastic concentration of CO2; Γ*, the CO2 compensation point and Kc and Ko, the Michaelis–Menten constants of Rubisco for CO2 and O2, respectively. The values used for Kc and Ko were 272 and 165 μmol mo1−1, respectively (Sharkey et al., 2007).
The RuBP-limited rate of CO2 assimilation (Aj) is given by:
where J is the rate of electron transport. The rate of electron transport is given by:
where Jmax is the maximum rate of electron transport; Q, the incident PAR; and α, the quantum efficiency which represents the initial slope of the photosynthetic light response curve. Therefore, α is a measure of the efficiency of light absorption by the leaf.
The mesophyll conductance (gm) is given by von Caemmerer (2000):
where Ci, the intercellular concentration of CO2.
The mesophyll to stomatal conductance ratio (gm/gs ratio) was calculated as gm divided by gs-value at the ambient atmospheric CO2 partial pressure (Ca = 400 μmol mol−1). Because gm and gs may vary partially independently, gm/gs ratio is increasingly recognized as a useful trait to improve WUE.
Needle Nitrogen Concentrations, SLA, and Carbon Isotope Discrimination
Following gas exchange measurements (the A/Ci curves), the shoots were carefully removed from the cuvette, and were immediately harvested, placed in plastic bags, and refrigerated (–20°C). Afterward, the needles were scanned and projected leaf area was measured using WinSeedle (Version 2007 Pro, Regent Instruments, Québec, Canada). After oven-drying (56°C for 72 h) and weighting, SLA was calculated as the ratio of the projected needles area (cm2) to the needles dry mass (g). Dried needles were ground to a fine powder in a ball mill. Needle nitrogen concentrations (Nmass: mg g−1) were determined by high-temperature combustion using a LECO elemental analyzer (Leco Corporation, St-Joseph, Michigan). Nitrogen concentration (Nmass) was then converted to a projected area basis (Na) using the SLA measurements for each sample. Photosynthetic nitrogen-use efficiency (PNUE) was calculated as Amax divided by Na (Amax/Na).
Carbon isotope discrimination analysis was performed on needles from the adjacent shoot at the same position in the plant used for gas exchange. After projected leaf area, weighing measurement and grinding, samples were sent to the University of California, Davis, Stable Isotope Facility. The samples were combusted in an online continuous flow dual analyser coupled to an isotope ratio mass spectrometer (Europa Scientific Integra, Cheshire, England, UK). The ratio of 13C to 12C was expressed as delta (δ13)-values in parts per million (‰) with respect to the Vienna Peedee Belemnite (VPDB) carbonate standard following Farquhar et al. (1982):
where (13C/12C)sample and (13C/12C)VPDB are the ratios of 13C/12C in the sample and the arbitrary standard (V-PDB, Vienna Peedee belemnite), respectively.
The accuracy of the δ13C analysis was about 0.05‰, as tested against four standards with known isotopic composition.
The discrimination against 13C (Δ13C) was calculated following (Farquhar et al., 1989):
where Δ13C is an index of discrimination between the two stable carbon isotopes which is generally used as an indicator of long-term water use efficiency (WUE). Theoretically, Δ13C is expected to scale negatively with WUE (Farquhar et al., 1989). δ13Cair is the carbon isotope ratio of atmospheric CO2 (−8‰).
Climatic data collected from adjacent meteorological stations for the period 1981–2010 were interpolated for each seed source (seed orchard) using the BioSIM software (Régnière and St-Amant, 2007). Latitude, longitude, and altitude of the centroid of plus-trees was used to simulate climatic data of each seed orchard (Table 1). This data set contains biologically relevant climatic indices derived from monthly temperature and precipitation, such as mean growing season temperature (MGST), total growing season precipitation (TGSP), mean annual temperature (MAT), total annual precipitation (TAP), number of growing degree-days ≥5°C (GDD5; Table 1), mean July temperature (MJT), frost-free days (Tmin > 0), total radiation during growing season, and an aridity index. Climates of seed origins (seed orchards) showed latitudinal and longitudinal gradients. Although temperature and precipitation varied with latitude and longitude, precipitation was strongly correlated to longitude of seed origin, while MGST, mean July temperature (MJT) and number of growing degree days ≥5°C (GDD5) were strongly correlated to latitude (Benomar et al., 2015).
All analyses were conducted with SAS version 9.4 (SAS Institute, Cary, NC, USA). Except for survival, response variables were analyzed separately by the MIXED procedure using a general linear mixed model with the effects of site and seed source (seed orchard) and their interaction on response variables as fixed effects and block as random effect. A significant seed source effect means that there is genetic differentiation among seed orchards, while a significant site effect is synonymous with trait plasticity and a significant site * seed source interaction effect means that plasticity varies among seed sources. Data were transformed when required to satisfy normality of residuals and homoscedasticity. The proportion of variance explained by each factor (site, seed source, and site * seed source) was estimated by PROC VARCOMP. Means comparisons were performed using the Tukey test; differences were considered significant at P = 0.05. Survival rates of seedlings at the end of the second growing season in the three sites were compared using the χ2-test (PROC FREQ). When ANOVAs showed a significant effect of seed source for a trait, the genetic differentiation along the climatic gradient was tested by univariate linear regression of that trait (mean values across sites) and climate of seed origin (MAP, GDDD, and MJT) using PROC REG (SAS Institute, Cary, NC, USA). The linear regressions were considered significant when P < 0.05 and marginally significant when P < 0.1. When ANOVAS showed a significant effect of the interaction between site and seed source, clinal variation of plasticity along the climatic gradient was tested by univariate linear regression of means of seed sources at each site for a given trait on climate of seed origin.
The phenotypic plasticity was assessed quantitatively for each trait using the relative distance plasticity index (RDPI) as described by Valladares et al. (2007):
where meanij is the mean value of trait of seed source i (i = 1, …., 6) growing at site j (j = 1, …, 3) and n is the total number of distances (n = 3). For each trait, one value of RDPI was calculated for each seed source across treatments which were considered as replicates. This allows statistical comparison of RDPI between traits using simple ANOVAs. As RDPI ranges from 0 (no plasticity) to 1 (maximal plasticity), data were transformed using logit transformation prior to the analysis. An analysis of covariance (ANCOVA) was used to assess the effect of site on the slope and intercept of the relationships between height (H2014) and functional traits.
Total height growth after two growing seasons (H2014) was significantly affected by both site and seed source. However, the interaction between site and seed source was not significant, indicating that the plasticity of seedlings from the six seed sources was similar in response to changes in site conditions along the climatic gradient (Table 2). Variance components showed that, on average, the plantation site effect explained the largest portion of variance (75%) while the seed source effect explained only 9% of the observed variance. Height growth was the greatest in Asselin (central) followed by Watford (south), and Deville (north), respectively (Figure 3). The seedlings from SO2 were significantly taller than those from SO5, and seedlings from other seed orchards (SO1, SO3, SO4, and SO6) had intermediate height values (Figure 3). There was a significant clinal variation among seed sources in height growth; in fact, height growth was negatively associated with latitude of seed origin (P < 0.06; R2 = 0.63) and marginally (P = 0.1) and positively associated with MJT (Figure 4A). Survival rate was similar among sites and seed sources, and it reached 97% at the end of the second growing season.
Table 2. Analysis of variance for growth and functional traits for sources of variation (sites, seed sources and the site * seed source interaction), degrees of freedom (df), F-values, and associated probabilities (P).
Figure 3. Mean total height growth at the end of the second growing season (H2014) of white spruce seedlings from six seed orchards at three plantation sites. Means having the same letters are not significantly different at α = 0.05. Values are the means ± SD.
Figure 4. Relationships between functional traits in six white spruce seed sources and the mean July temperature at their origin (MJT) after growing on three forest sites. (A) Total height at the end of the second growing season (H2014), (B) light-saturated photosynthesis on a projected leaf area basis (Amax), (C) stomatal conductance (gs), (D) mesophyll to stomatal conductance ratio (gm/gs), (E) water use efficiency (WUE), and (F) photosynthetic nitrogen-use efficiency (PNUE). Data points are means of 12 seedlings for functional traits and 768 seedlings for growth (H2014). Regression lines and coefficient of determination (R2) values are shown for significant relationships, with solid line for P < 0.05 and dashed lines for P < 0.1 (n = 6).
Functional Traits: Local Adaptation vs. Plasticity
All measured ecophysiological traits were under the influence of both site (plasticity) and seed source (genetic differentiation) effects. However, and similarly as for height growth, the interaction between site and seed source was not significant for any of the measured traits, indicating a similar plasticity of traits between seed sources in response to the climatic gradient (Table 2). The percent of variance explained by the site effect ranged from 43 to 75% among measured traits.
Light-saturated photosynthetic rate (Amax) expressed on both a projected area (Amax∕area) and a mass basis (Amax∕mass) showed plasticity in response to changes in growing conditions along the climatic gradient (Table 3). Amax∕area and Amax∕mass were greater at Asselin (central) than Watford (south) and Deville (north) (Table 3). There was a significant effect of seed source on Amax (Table 2). Variation in light-saturated photosynthetic rate on a projected area (Amax∕area) was negatively associated with latitude of seed origin (P < 0.05; R2 = 0.75) and positively with MJT (Figure 4B). Stomatal and mesophyll conductances (gs and gm) were significantly impacted by site conditions. However, only gs was significantly influenced by seed source (Table 2). Stomatal conductance was greater at Asselin (central) followed by Watford (south), and Deville (north), respectively. Mesophyll conductance was the greatest in Asselin (central; Table 3), while the mesophyll to stomatal conductance ratio (gm/gs ratio) was the greatest in Deville (north; Table 3). As for Amax, gs was negatively associated with latitude of seed origin (P < 0.05; R2 = 0.78) and positively with MJT (Figure 4C). gm/gs ratio was marginally negatively associated with MJT (Figure 4D). Biochemical limitations to photosynthesis (Vcmax and Jmax) were significantly affected by site, but not by the seed source (Table 2). Vcmax and Jmax were the lowest at Watford (south; Table 3). We found a significant effect of site, but not of the seed source, on SLA and needle nitrogen content (Nmass) (Table 3). Nmass was the lowest (30%) in Watford (south) and SLA was the greatest in Deville (north; Table 3).
Site significantly influenced carbon isotope discrimination (Δ13C) with a substantial difference between the three sites but unexpectedly, no difference was detected among seed sources (Table 2). Intrinsic water use efficiency (WUE) and photosynthetic nitrogen-use efficiency (PNUE) were under the influence of site and the seed source (Table 2). Deville (north) showed the greatest WUE and the lowest PNUE-values (Table 3). Despite the presence of significant variation in WUE and PNUE among seed sources, the latter were not associated with latitude of seed origin or associated climatic variables (Figures 4E,F).
Level of Plasticity between Functional Traits
As assessed by the relative distance plasticity index (RDPI) across all seed sources, plasticity varied between traits. Values of RDPI ranged from 0.03 to 0.17. Among needle functional traits, stomatal conductance, mesophyll conductance, and Nmass were the most plastic traits along the climatic gradient with RDPI = 0.14, 0.11, and 0.10, respectively. Needle morphology (SLA), Δ13C and biochemical limitation to photosynthesis (Vcmax and Jmax) were the least plastic traits (Figure 5). All other traits were of intermediate plasticity.
Figure 5. Comparison of phenotypic plasticity among traits across all seed sources as assessed by relative distance plasticity index (RDPI). Means having the same letters are not significantly different at α = 0.05 (Tukey post hoc). Values are the means ± SD. (see Table 2 for abbreviations).
Phenotypic Correlations among Traits
Photosynthetic capacity on a projected area (Amax∕area) was positively correlated with stomatal conductance (gs), mesophyll conductance (gm), gm/gs ratio, needle nitrogen concentration area basis (Narea), and PNUE. On the other hand, Amax∕area was negatively correlated with Rd and WUE (Table 4). Rd was strongly and positively correlated with Vcmax, Jmax and needle nitrogen concentrations on both area and mass basis (Narea and Nmass). Stomatal and mesophyll conductance were positively correlated. Water use efficiency (WUE) was correlated negatively with gs, but not with gm and it was correlated positively with the gm/gs ratio (Table 4). PNUE was correlated to gm/gs ratio, and Nmass. The trade-off between PNUE and WUE was stronger. Δ13C was correlated with Rd, gm/gs ratio, Vcmax, Jmax, SLA, WUE, and PNUE (Table 4).
Table 4. Pearson correlations (r)a among seed source mean values per plantation site for functional traits (n = 18).
Height was positively related to Amax∕area, and gs and negatively to gm/gs ratio and WUE. Height–Amax∕area, height–gs, and Height–WUE relationships were insensitive to environment conditions (site of plantation; Figures 6A–C). The intercept was significantly greater in Asselin (central) for the relationship of H2014 with gm/gs ratio (Figure 6D). H2014 was unrelated to all other measured functional traits (P ≥ 0.05).
Figure 6. Total height growth at the end of the second growing season (H2014; an index for fitness) plotted against (A) light-saturated photosynthesis on a projected leaf area basis (Amax) and (B) stomatal conductance (gs), (C) water use efficiency (WUE), (D) mesophyll to stomatal conductance ratio (gm/gs). Data points are seed source-specific mean value at each site (n = 18). The strength of relationships is indicated by the p-value and coefficient of determination. The significant effect of site on slope and intercept of relationships was determined by ANCOVA. Where slope and/or intercept were significantly different, regression line and coefficient of determination are given for each relationships and the P-value is given for the common regression. NS, not significant.
Growth and Performance of White Spruce Seedlings from Different Seed Sources under Different Site Conditions
Our study revealed a clinal gradient in height growth at a regional scale among the six tested seed sources (Figure 4A). Irrespective of growing conditions (plantation site), southern seed sources grew faster than northern seed sources (Figure 3). In addition, the survival rate was very high (97%) and similar among seed sources. These results are in agreement with previous reports showing negative latitudinal cline of height growth for young white spruce at regional (Li et al., 1997; Jaramillo-Correa et al., 2001; Thomson et al., 2010; Carles et al., 2011) and range-wide scales (Lu et al., 2014). A negative latitudinal cline in growth rate has also been reported for other boreal conifer species (Matyas and Yeatman, 1992; Wei et al., 2004). Stott and Loehle (1998), suggesting that this pattern may result from the pressure exercised by both biotic and abiotic factors. Populations from the southern range of the species may exhibit a competitive growth strategy as a consequence of higher level of biotic pressures (interspecific competition), while northern populations may exhibit a more conservative growth strategy focused on survival growth strategy as a consequence of higher level of abiotic pressure (cold damage, drought). Our results suggest that the observed trend in height growth with latitude may be driven by the trade-off between Amax and WUE mediated by genetic variation in gs and gm/gs ratio (Figure 6) rather than the trade-off between growth and survival.
The growth of seed sources along the climatic gradient followed a parabolic pattern. Height growth was 30% greater for the site with the intermediate climate (Asselin) followed by the warmest (Watford) and the coldest (Deville) sites, respectively. The lower height growth at the warmest site (Watford) may be interpreted as an existing maladaptation of tested seed sources to warming conditions in southern Québec (Rainville et al., 2014). Considering previous findings showing a strong age-age-genetic correlation of height growth in white spruce (Li et al., 1993; Wahid et al., 2013), one may believe that this observed parabolic growth pattern could be stable over the next years during the juvenile stage. Longer-term data are however required to confirm this pattern, although a similar pattern was observed with older white spruce seed sources on different sites (Andalo et al., 2005).
Functional Trait Expression: Local Adaptation vs. Plasticity
We found a significant effect of site and seed source on Amax, gs, gm/gs ratio, WUE, and PNUE. However, the interaction between site and seed source was not significant suggesting a lack of differential plasticity among seed sources. The low variation explained by seed sources (9–18%) for the above-mentioned traits suggests that genetic differentiation in these traits is limited. Nevertheless, the results indicate that white spruce populations in Québec are locally adapted despite the limited latitudinal range sampled in the present study and the large gene flow existing in this widespread boreal tree species (Jaramillo-Correa et al., 2001). The results of the present study are in agreement with previous findings showing a significant genetic differentiation among populations in the same region for growth, phenology, morphology, and physiology traits (Li et al., 1997; Namroud et al., 2008; Carles et al., 2011; Benomar et al., 2015).
Photosynthetic rate (Amax) showed a clinal variation with southern populations exhibiting a higher Amax than northern populations (Figure 4B). This pattern is in agreement with previous findings (Kogami et al., 2001; Dang et al., 2008; Qiuhong et al., 2013; Benomar et al., 2015). However, the results are in conflict with others that demonstrated a positive latitudinal cline in Amax driven by needle nitrogen (Nmass) and carboxylation capacity (Vcmax) (Oleksyn et al., 1998; Soolanayakanahally et al., 2009). Our results showed no genetic variation in respiration rate (Rd), gm, Vcmax, or Nmass. On the other hand, Amax was strongly correlated with gs, gm, and gm/gs ratio, but not with Vcmax and Jmax. Therefore, the latitudinal cline in Amax in this study likely resulted from a parallel cline in CO2 diffusion (gs and gm/gs ratio). Overall, our study revealed that the variation among white spruce seed sources in height growth, Amax and resources use efficiencies was driven by the physics of CO2 diffusion (CO2 transfer conductance).
Our results are in agreement with the fact that stomatal conductance (gs) is under both environmental and genetic controls. The decrease in stomatal conductance and the increase in gm/gs ratio and WUE along with latitude may arise from a water conservation strategy in northern latitude as an adaptation of stomatal characteristics (stomatal density, stomatal length) and hydraulic conductance to a decrease in soil temperature (Wang et al., 2014). In fact, mean soil temperature during the two growing seasons of experiment was 2.4°C lower in the coldest site (Deville). The decrease in soil temperature with latitude would impose an apparent drought stress likely induced by an increase in water viscosity and a decrease in water and nutrient uptake by roots at cold temperatures (Cochard et al., 2000; Luo et al., 2006; Wang et al., 2014). Under boreal forest conditions, low soil temperature has been reported to decrease both stomatal conductance and root growth (Lamhamedi and Bernier, 1994; Grossnickle, 2000; Dang and Cheng, 2004). Thus, the effect of soil temperature on stomatal conductance may arise from combination of poor root growth, low water uptake and plant resistance to water flow (Grossnickle, 1988).
SLA is an adaptive morphological trait determining plant capacity to succeed in specific biophysical environments. Indeed, SLA influences Amax, Rd, Vcmax Nmass, and leaf life-span (Wright et al., 2004) and SLA has been reported to be under both genetic and environmental control. A decrease in SLA with increases in latitude and elevation has been reported for several conifer species (Friend et al., 1989; Oleksyn et al., 1998; Hultine and Marshall, 2000; Qiuhong et al., 2013). In the present study, we observed only a small change (14%) in SLA along the climatic gradient. As a result, this trait was less informative regarding phenotypic variation in the physiological traits assessed (Table 4). Changes in SLA may result from changes in needle density and/or thickness. Also, needle density and thickness may change in opposite directions leading to unchanged SLA. Thus, a change in needle morphology may possibly occur but not be captured by SLA. Otherwise, the climatic gradient in our study may likely be insufficient to induce a significant change in SLA.
In agreement with previous reports, we found an environmental influence on carbon isotope composition as demonstrated by the significant differences between sites (Körner et al., 1991; Cernusak et al., 2013; McLean et al., 2014). Δ13C increased linearly from the warmest site to the coldest site by 1.4‰ (Kogami et al., 2001). The observed variation in Δ13C in this study resulted from a change in gm/gs ratio, Vcmax, and Rd. The lack of variation in Δ13C among seed sources was likely due to the lack of variation between seed sources in Vcmax and Rd. Our results demonstrate that the change in Δ13C is influenced by both biophysical and biochemical processes (Farquhar et al., 1989; Cernusak et al., 2013). The latter may explain the significant correlation between Δ13C and both WUE and PNUE.
Implications for Assisted Population Migration
The southern transfer of northern seed sources led to an enhancement of growth which was still below that of southern seed sources (Figure 3). These results suggest that as temperature warms, growth of northern seed sources of white spruce in Quebec would increase but without reaching that of seed sources already adapted to these new conditions, at least at the young age. In addition, a trade-off between growth and survival did not occur along the tested climatic gradient. Thus, if the objective is to improve forest plantation productivity under climate change, northward assisted migration of southern seed sources (for instance, by 2° of latitude northward) might be possible without undue risk. Furthermore, the faster growth rate of southern seed sources during the juvenile stage may contribute to the establishment success of plantation during assisted migration by avoiding negative impacts of biotic and abiotic factors such as weedy competition, early and late frost damage, and drought stress. The greatest growth rate of southern seed sources seems to have been reached to the detriment of WUE (Figure 6). For this reason and as it was noted by Andalo et al. (2005) that local moisture conditions may also play an important role in white spruce seed source variation, it would be advisable to choose plantation sites based on projected precipitation before any implementation of long distance transfer of white spruce seed sources, as suggested recently by Lu et al. (2014). On the other hand, investigations regarding variation in drought tolerance among seed sources are necessary. In particular, (i) genetic variation and plastic changes in the vulnerability of xylem to cavitation and embolism resistance under water stress, and (ii) the linkage between hydraulic functioning and stomatal conductance should be further studied.
The present study was motivated by the need to examine the environmental and genetic influences on functional traits variation among white spruce seed sources to assist reforestation decisions under a changing climate. Our study revealed that (i) the seed sources differed in growth, photosynthetic rate (Amax) and its biophysical limitations, and resources use efficiencies; (ii) the seed sources demonstrated similar levels of traits plasticity, suggesting that the observed variation in growth among seed sources resulted largely from genetic variation in functional traits rather than from plasticity. Overall, the results suggest that phenotypic plasticity may allow northern white spruce seed sources to cope over the short-term with a moderate climate change. However, assisted migration under the form of northward seed source transfer from southern areas would still be valuable to increase forest productivity in the northern regions of Québec.
The results of this study were obtained on young material at the juvenile stage. While some short-term conclusions can be drawn with regard to APM, but a long-term follow-up remains necessary to confirm these preliminary findings. Moreover, access to data from additional test sites would help strengthen trends observed. In relation to this, we are in the process of establishing six new tests in order to delineate APM guidelines based on more exhaustive scientific evidence for white spruce in eastern Canada.
HM, ML, JBo, JBe, AR and LB conceived the study, obtained the funding and participated in the drafting of the manuscript. AR participated to designing the genetic tests and providing the seed. ML and his team have installed three weather stations for the acquisition of environmental variables and conducted mineral analyses. LB performed the measurements of different variables at different sites, the analysis and the interpretation of data. All authors read and approved the final version of the manuscript.
This research was funded by grants to HM from the program “Partenariat sur L'aménagement et L'ennvironnement Forestiers” of the Fonds de la Recherche du Québec sur la Nature et les Technologies (FRQ-NT), and the Discovery Grant program of the Natural Sciences and Engineering Research Council of Canada. Major additional support was also provided to ML by the Direction de la Recherche Forestière of the Ministère des Forêts, de la Faune et des Parcs of Québec.
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 Guildo Gagnon, Mario Renaud, and Pascal Desjardins (Ministère des Forêts, de la Faune et des Parcs du Québec), Bakry Mustapha and Marie R. Coyea (University Laval) for their technical assistance throughout the project. This study was conducted in collaboration with the Direction de la Recherche Forestière of the Ministère des Forêts, de la Faune et des Parcs du Québec.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.00048
Aitken, S. N., Yeaman, S., Holliday, J. A., Wang, T., and Curtis-McLane, S. (2008). Adaptation, migration or extirpation: climate change outcomes for tree populations. Evol. Appl. 1, 95–111. doi: 10.1111/j.1752-4571.2007.00013.x
Andalo, C., Beaulieu, J., and Bousquet, J. (2005). The impact of climate change on growth of local white spruce populations in Quebec, Canada. For. Ecol. Manag. 205, 169–182. doi: 10.1016/j.foreco.2004.10.045
Anderson, J. T., and Gezon, Z. J. (2015). Plasticity in functional traits in the context of climate change: a case study of the subalpine forb Boechera stricta (Brassicaceae). Glob. Change. Biol. 21, 1689–1703. doi: 10.1111/gcb.12770
Beaulieu, J., Daoust, G., Deshaies, A., Lamhamedi, M. S., Rainville, A., and Tourigny, M. (2009). “Amélioration génétique des arbres, gestion des vergers à graines et de semences, et production de plants forestiers,” in Manuel de foresterie, ed Ordre. des. ingénieurs. forestiers. du. Québec (Québec, QC: Presses de l'Université Laval), 1093–1146.
Benomar, L., Lamhamedi, M. S., Villeneuve, I., Rainville, A., Beaulieu, J., Bousquet, J., et al. (2015). Fine-scale geographic variation in photosynthetic-related traits of Picea glauca seedlings indicates local adaptation to climate. Tree Physiol. 35, 864–878. doi: 10.1093/treephys/tpv054
Bourque, A., and Simonet, G. (2007). “Québec,” in From Impacts to Adaptation: Canada in a Changing Climate, eds D. S. Lemmen, F. J. Warren, J. Lacroix, and E. Bush (Ottawa, ON: Government of Canada), 171–226.
Carles, S., Lamhamedi, M. S., Beaulieu, J., Stowe, D. C., and Margolis, H. A. (2011). Differences in growth and mineral nutrition of seedlings produced from ten white spruce seed orchards. New For. 42, 195–214. doi: 10.1007/s11056-011-9247-y
Cernusak, L. A., Ubierna, N., Winter, K., Holtum, J. A. M., Marshall, J. D., and Farquhar, G. D. (2013). Environmental and physiological determinants of carbon isotope discrimination in terrestrial plants. New Phytol. 200, 950–965. doi: 10.1111/nph.12423
Cochard, H., Martin, R., Gross, P., and Bogeat−Triboulot, M. B. (2000). Temperature effects on hydraulic conductance and water relations of Quercus robur L. J. Exp. Bot. 51, 1255–1259. doi: 10.1093/jexbot/51.348.1255
Cordell, S., Goldstein, G., Mueller-Dombois, D., Webb, D., and Vitousek, P. M. (1998). Physiological and morphological variation in Metrosideros polymorpha, a dominant Hawaiian tree species, along an altitudinal gradient: the role of phenotypic plasticity. Oecologia 113, 188–196. doi: 10.1007/s004420050367
Dang, Q. L., Maepea, J. M., and Parker, W. H. (2008). Genetic variation of ecophysiological responses to CO2 in Picea glauca seedlings. Open. Forest. Sci. J. 1, 68–79. doi: 10.2174/1874398600801010068
de Lafontaine, G., Turgeon, J., and Payette, S. (2010). Phylogeography of white spruce (Picea glauca) in eastern North America reveals contrasting ecological trajectories. J. Biogeogr. 37, 741–751. doi: 10.1111/j.1365-2699.2009.02241.x
Ethier, G. J., Livingston, N. J., Harrison, D. L., Black, T. A., and Moran, J. A. (2006). Low stomatal and internal conductance to CO2 versus Rubisco deactivation as determinants of the photosynthetic decline of ageing evergreen leaves. Plant. Cell. Environ. 29, 2168–2184. doi: 10.1111/j.1365-3040.2006.01590.x
Farquhar, G. D., Ehleringer, J. R., and Hubick, K. T. (1989). Carbon isotope discrimination and photosynthesis. Annu. Rev. Plant Physiol. Plant Mol. Biol. 40, 503–537. doi: 10.1146/annurev.pp.40.060189.002443
Farquhar, G. D., O'Leary, M. H., and Berry, J. A. (1982). On the relationship between carbon isotope discrimination and the intercellular carbon dioxide concentration in leaves. Aust. J. Plant Physiol. 9, 121–137. doi: 10.1071/PP9820121
Flexas, J., Ribas-Carbo, M., Diaz-Espejo, A., Galmes, J., and Medrano, H. (2008). Mesophyll conductance to CO2: current knowledge and future prospects. Plant. Cell. Environ. 31, 602–621. doi: 10.1111/j.1365-3040.2007.01757.x
Friend, A. D., Woodward, F. I., and Switsur, V. R. (1989). Field measurements of photosynthesis, stomatal conductance, leaf nitrogen and δ13C along altitudinal gradients in Scotland. Funct. Ecol. 3, 117–122. doi: 10.2307/2389682
Isaac-Renton, M. G., Roberts, D. R., Hamann, A., and Spiecker, H. (2014). Douglas-fir plantations in Europe: a retrospective test of assisted migration to address climate change. Glob Change Biol. 20, 2607–2617. doi: 10.1111/gcb.12604
Jaramillo-Correa, J. P., Beaulieu, J., and Bousquet, J. (2001). Contrasting evolutionary forces driving population structure at expressed sequence tag polymorphisms, allozymes and quantitative traits in white spruce. Mol. Ecol. 10, 2729–2740. doi: 10.1046/j.0962-1083.2001.01386.x
Kogami, H., Hanba, Y. T., Kibe, T., Terashima, I., and Masuzawa, T. (2001). CO2 transfer conductance, leaf structure and carbon isotope composition of Polygonum cuspidatum leaves from low and high altitudes. Plant. Cell. Environ. 24, 529–537. doi: 10.1046/j.1365-3040.2001.00696.x
Kranabetter, J. M., Stoehr, M., and O'Neill, G. A. (2015). Ectomycorrhizal fungal maladaptation and growth reductions associated with assisted migration of Douglas-fir. New Phytol. 206, 1135–1144. doi: 10.1111/nph.13287
Lamhamedi, M. S., Labbé, L., Margolis, H. A., Stowe, D. C., Blais, L., and Renaud, M. (2006). Spatial variability of substrate water content and growth of white spruce seedlings. Soil Sci. Soc. Am. J. 70, 108–120. doi: 10.2136/sssaj2005.0109
Li, P., Beaulieu, J., and Bousquet, J. (1997). Genetic structure and patterns of genetic variation among populations in eastern white spruce (Picea glauca). Can. J. For. Res. 27, 189–198. doi: 10.1139/x96-159
Lu, P., Parker, W. H., Cherry, M., Colombo, S., Parker, W. C., Man, R., et al. (2014). Survival and growth patterns of white spruce (Picea glauca [Moench] Voss) rangewide provenances and their implications for climate change adaptation. Ecol. Evol. 4, 2360–2374. doi: 10.1002/ece3.1100
Luo, J., Zang, R., and Li, C. (2006). Physiological and morphological variations of Picea asperata populations originating from different altitudes in the mountains of southwestern China. For. Ecol. Manag. 221, 285–290. doi: 10.1016/j.foreco.2005.10.004
McLean, E. H., Prober, S. M., Stock, W. D., Steane, D. A., Potts, B. M., Vaillancourt, R. E., et al. (2014). Plasticity of functional traits varies clinally along a rainfall gradient in Eucalyptus tricarpa. Plant. Cell. Environ. 37, 1440–1451. doi: 10.1111/pce.12251
Miao, Z., Xu, M., Lathrop, R. G., and Wang, Y. (2009). Comparison of the A–Cc curve fitting methods in determining maximum ribulose 1·5-bisphosphate carboxylase/oxygenase carboxylation rate, potential light saturated electron transport rate and leaf dark respiration. Plant. Cell. Environ. 32, 109–122. doi: 10.1111/j.1365-3040.2008.01900.x
Namroud, M.-C., Beaulieu, J., Juge, N., Laroche, J., and Bousquet, J. (2008). Scanning the genome for gene single nucleotide polymorphisms involved in adaptive population differentiation in white spruce. Mol. Ecol. 17, 3599–3613. doi: 10.1111/j.1365-294X.2008.03840.x
Oleksyn, J., Modrzyñski, J., Tjoelker, M. G., Zytkowiak, R., Reich, P. B., and Karolewski, P. (1998). Growth and physiology of Picea abies populations from elevational transects: common garden evidence for altitudinal ecotypes and cold adaptation. Funct. Ecol. 12, 573–590. doi: 10.1046/j.1365-2435.1998.00236.x
O'Neill, G. A., Hamann, A., and Wang, T. (2008). Accounting for population variation improves estimates of the impact of climate change on species' growth and distribution. J. Appl. Ecol. 45, 1040–1049. doi: 10.1111/j.1365-2664.2008.01472.x
O'Neill, G. A., and Nigh, G. (2011). Linking population genetics and tree height growth models to predict impacts of climate change on forest production. Glob. Change Biol. 17, 3208–3217. doi: 10.1111/j.1365-2486.2011.02467.x
Pedlar, J., McKenney, D., Beaulieu, J., Colombo, S., McLachlan, J., and O'Neill, G. (2011). The implementation of assisted migration in Canadian forests. For. Chron. 87, 766–777. doi: 10.5558/tfc2011-093
Pickles, B. J., Twieg, B. D., O'Neill, G. A., Mohn, W. W., and Simard, S. W. (2015). Local adaptation in migrated interior Douglas-fir seedlings is mediated by ectomycorrhizas and other soil factors. New Phytol. 207, 858–871. doi: 10.1111/nph.13360
Prober, S. M., Byrne, M., McLean, E. H., Steane, D. A., Potts, B. M., Vaillancourt, R. E., et al. (2015). Climate-adjusted provenancing: a strategy for climate-resilient ecological restoration. Front. Ecol. Evol. 3:65. doi: 10.3389/fevo.2015.00065
Rainville, A., Beaulieu, J., Langevin, L., Logan, T., and Lambert, M. (2014). Prédire L'effet des Changements Climatiques Sur le Volume Marchand des Principales Espèces Résineuses Plantées au Québec, Grâce à la Génétique Forestière. Québec, QC: Min Forêts Faune Parcs Québec.
Régnière, J., and St-Amant, R. (2007). Stochastic simulation of daily air temperature and precipitation from monthly normals in North America north of Mexico. Int. J. Biometeorol. 51, 415–430. doi: 10.1007/s00484-006-0078-z
Roy, P., Gachon, P., and Laprise, R. (2014). Sensitivity of seasonal precipitation extremes to model configuration of the Canadian Regional Climate Model over eastern Canada using historical simulations. Clim. Dyn. 43, 2431–2453. doi: 10.1007/s00382-014-2066-z
Sharkey, T. D., Bernacchi, C. J., Farquhar, G. D., and Singsaas, E. L. (2007). Fitting photosynthetic carbon dioxide response curves for C3 leaves. Plant. Cell. Environ. 30, 1035–1040. doi: 10.1111/j.1365-3040.2007.01710.x
Soolanayakanahally, R. Y., Guy, R. D., Silim, S. N., Drewes, E. C., and Schroeder, W. R. (2009). Enhanced assimilation rate and water use efficiency with latitude through increased photosynthetic capacity and internal conductance in balsam poplar (Populus balsamifera L.). Plant. Cell. Environ. 32, 1821–1832. doi: 10.1111/j.1365-3040.2009.02042.x
Ste-Marie, C. (2014). Adapting Sustainable Forest Management to Climate Change: A Review of Assisted Tree Migration and Its Potential Role in Adapting Sustainable Forest Management to Climate Change. Ottawa, ON: Canadian Council of Forest Ministers.
Valladares, F., Matesanz, S., Guilhaumon, F., Araújo, M. B., Balaguer, L., Benito-Garzón, M., et al. (2014). The effects of phenotypic plasticity and local adaptation on forecasts of species range shifts under climate change. Ecol. Lett. 17, 1351–1364. doi: 10.1111/ele.12348
Wahid, N., Lamhamedi, M. S., Rainville, A., Beaulieu, J., and Margolis, H. A. (2013). Genetic control and nursery-plantation genotypic correlations for growth characteristics of white spruce somatic clones. J. Sustain. Forest. 32, 576–593. doi: 10.1080/10549811.2013.791231
Wang, R., Yu, G., He, N., Wang, Q., Xia, F., Zhao, N., et al. (2014). Elevation-related variation in leaf stomatal traits as a function of plant functional type: evidence from Changbai mountain, China. PLoS ONE 9:e115395. doi: 10.1371/journal.pone.0115395
Warren, C. R., and Adams, M. A. (2006). Internal conductance does not scale with photosynthetic capacity: implications for carbon isotope discrimination and the economics of water and nitrogen use in photosynthesis. Plant. Cell. Environ. 29, 192–201. doi: 10.1111/j.1365-3040.2005.01412.x
Wei, R.-P., Han, S. D., Dhir, N. K., and Yeh, F. C. (2004). Population variation in growth and 15-year-old shoot elongation along geographic and climatic gradients in black spruce in Alberta. Can. J. For. Res. 34, 1691–1702. doi: 10.1139/x04-050
Yang, J., Pedlar, J. H., McKenney, D. W., and Weersink, A. (2015). The development of universal response functions to facilitate climate-smart regeneration of black spruce and white pine in Ontario, Canada. For. Ecol. Manag. 339, 34–43. doi: 10.1016/j.foreco.2014.12.001
Keywords: white spruce, climate change, assisted migration, local adaptation, acclimation, functional traits, mesophyll conductance, photosynthetic capacity
Citation: Benomar L, Lamhamedi MS, Rainville A, Beaulieu J, Bousquet J and Margolis HA (2016) Genetic Adaptation vs. Ecophysiological Plasticity of Photosynthetic-Related Traits in Young Picea glauca Trees along a Regional Climatic Gradient. Front. Plant Sci. 7:48. doi: 10.3389/fpls.2016.00048
Received: 29 November 2015; Accepted: 12 January 2016;
Published: 03 February 2016.
Edited by:Judy Simon, University of Konstanz, Germany
Reviewed by:Rubén Retuerto, Universidad de Santiago de Compostela, Spain
Qinglai Dang, Lakehead University, Canada
Copyright © 2016 Benomar, Lamhamedi, Rainville, Beaulieu, Bousquet and Margolis. 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: Lahcen Benomar, firstname.lastname@example.org