- 1State Key Laboratory of Tree Genetics and Breeding, Key Laboratory of Tree Breeding and Cultivation of State Forestry Administration, Research Institute of Forestry, Chinese Academy of Forestry, Beijing, China
- 2Research Institute of Forest, Ecology and Environment Protection, Chinese Academy of Forestry, Beijing, China
- 3Guangdong Institute of Eco-environmental Science and Technology, Guangzhou, China
- 4Swiss Federal Research Institute WSL, Birmensdorf, Switzerland
- 5School of Geographical Sciences, Northeast Normal University, Changchun, China
Understanding the mechanisms that determine plant distribution range is crucial for predicting climate-driven range shifts. Compared to altitudinal gradients, less attention has been paid to the mechanisms that determine latitudinal range limit. To test whether intrinsic resource limitation contributes to latitudinal range limits of woody species, we investigated the latitudinal variation in non-structural carbohydrates (NSC; i.e., total soluble sugar plus starch) and nutrients (nitrogen and phosphorus) in mature and juvenile Chinese cork oak (Quercus variabilis Blume) along a 1500 km north-south transect in China. During the growing season and dormant season, leaves, branches, and fine roots were collected from both mature and juvenile oaks in seven sites along the transect. Tissue concentration of NSCs, N, and P did not decrease with increasing latitude irrespective of sampling season and ontogenetic stage. Furthermore, higher levels of NSCs and N in tissues of juveniles relative to mature trees were found during the dormant season. Partial correlation analysis also revealed that during the dormant season, soluble sugar, NSC, the ratio of soluble sugar to starch, and tissue nitrogen concentration were correlated positively with latitude but negatively with precipitation and mean temperature of dormant season. Our results suggest that carbon or nutrient availability may not be the driving factors of the latitudinal range limit of the studied species. Further studies should be carried out at the community or ecosystem level with multiple species to additionally test the roles of factors such as regeneration, competition, and disturbance in determining a species’ northern distribution limit.
Introduction
In the context of climate change, altitudinal and latitudinal gradients can be used as natural laboratories to deduce species’ responses to global warming (Körner, 2007; De Frenne et al., 2013). Nowadays, northward and upward shift of species distribution has been observed in a wide range of studies (Kelly and Goulden, 2008; Chen et al., 2011; Boisvert-Marsh et al., 2014; Desprez et al., 2014; Rinas et al., 2017; Sittaro et al., 2017), but the underlying physiological mechanisms are under debate. Therefore, we may use geographical gradients to explore the mechanisms for species’ distribution under current conditions, and thus to understand and forecast the responses of species’ distribution to predicted climate change.
Compared to a few eco-physiological evidences for the formation of a species’ northernmost distribution, altitudinal tree-line formation has attracted lots of studies during the last several decades, and several important hypotheses have been proposed (Li and Krauchi, 2005). Notably, two mutually exclusive hypotheses: the Carbon Limitation Hypothesis (CLH) (Stevens and Fox, 1991) and the Growth Limitation Hypothesis (GLH) (Körner, 1998; Hoch and Körner, 2003), have been extensively tested, as they both have the potential to be applied worldwide (Körner, 1998) due to intrinsic response or adaptation strategies of plants to a variety of environmental factors at the alpine tree-line. Non-structural carbohydrates (NSC, sum of starch and total soluble sugars) and nutrients (nitrogen or phosphorus) of plant tissues were generally used to evaluate the carbon or nutrient status along altitudinal gradients to test the hypotheses mentioned above (Hoch and Körner, 2003, 2012; Li et al., 2008a; Sullivan et al., 2014; Fajardo and Piper, 2017). NSC, formed during photosynthesis, could act as substrate for respiration to provide the energy needed for growth and maintenance processes (Körner, 2003); and the level of its two main components, soluble sugars and starch, could reflect the balance between carbon gain and carbon utilization and loss within a plant, representing a tree’s capital for growth and acting as a buffer during insufficient source activities due to environmental stress (Li et al., 2002; O’Brien et al., 2014). Meanwhile, nitrogen and phosphorus are the two most limiting elements to terrestrial vegetation (Reich and Oleksyn, 2004), which are not only directly related to carbon assimilation and allocation but also indirectly to stress tolerance or resistance (DeHayes et al., 1989; Villagra et al., 2013; Yan et al., 2016). However, at a single species level, whether the factors in determining the latitudinal distribution range are also associated with carbon and nutrient status still remain unclear.
Apart from the vital roles of carbon and nutrients to support plant survival and growth, ontogenetic variations in responses to environment variability or stress have gained more attention (Niinemets, 2010; le Roux et al., 2013; Klockmann et al., 2017). To date, numerous studies have already reported ontogenetic variations in carbon assimilation and allocation (Portsmuth et al., 2005; Steppe et al., 2011), resource use strategies (Gedroc et al., 1996; Rivas-Ubach et al., 2012), and stress tolerances (Cavender-Bares and Bazzaz, 2000; Niinemets, 2010). Compared to adult trees, however, less attention has been paid to early life stages (e.g., seedlings or saplings), which are more sensitive to environmental changes or stresses (Niinemets, 2010). Although adult trees may persist over hundreds of years, if seedlings or saplings fail to keep pace with the rate of rising temperature or climate change (Zhu et al., 2012; Bell et al., 2014; Máliš et al., 2016), in the long-term the distribution range of a species will shrink. Hence, to predict the responses of tree species to climate change, disentangling ontogenetic variations along environmental gradients is indispensable for a more profound understanding of adaptation strategies among different life stages, as well as for a better understanding of the mechanisms in determining the latitudinal range.
The Chinese cork oak (Quercus variabilis Blume) is one of the most widely distributed tree species in eastern Asia, stretching from approximately 24° to 40°N and 97° to 140°E (Chen et al., 2012). The different regions where the species is present differ greatly in terms of climatic and edaphic conditions, providing an ideal situation for studying the spatial patterns of carbon or nutrient allocation within a single widespread species. In the present study, leaves, branches, and fine roots in healthy juvenile and mature trees in seven sites along a designed 1500 km north-south transect in China, were collected and analyzed. Our aims were to answer the questions of (1) how the carbon or nutrient status vary along the south-north transect; and (2) whether there is a consistent difference in NSCs and nutrients between juvenile and mature trees across latitudes. Our results will contribute toward understanding and predicting the impacts of climate change on the range dynamics of woody plants along latitudinal gradients.
Materials and Methods
Sampling Protocol
A south-north transect covering 14 latitudinal degrees (∼26°–40°N) (c.1500 km) was setup in the natural distribution range of Q. variabilis in China, and seven sites were designed along the transect with a latitudinal interval of ∼2° (Figure 1 and Supplementary Table S1). Using the geographical location-based method by Gong and Jian (1983), we calculated the phenological date for each site (Supplementary Table S1), so that we collected samples at the same phenological stage across the seven sites to avoid effects of plant phenological variation or bias on carbon or nutrient status (Bansal and Germino, 2009; Bazot et al., 2013). The sampling tasks were conducted on two age classes (juvenile and mature) of natural stands during the mid-growing season (August, 2014) and the dormant season (Dec, 2014). At each site, three 50 m × 50 m temporal plots (n = 3), with a minimum distance of 10 km from each other, were set up. Within each plot, three to four canopy trees of the same age class, without browsing and other damages, were selected. The sampling processes were identical for the two sampling seasons. From each sample tree, three to four upper and outer-most sun-exposed branches were cut to collect the leaves (only for growing season) and branches (1-2a). Fine roots (<5 mm in diameter) attached to coarse roots of each sample tree were manually excavated using a mini-spade and carefully collected. All samples were stored in a cool box until they were taken to the laboratory. The same tissue from individuals in the same age class collected from the same site was pooled as one mixed sample. All samples were heated in a microwave oven at 600 W for 60 s, and then dried at 65°C for 72 h and ground to pass a 0.20 mm sieve for further analyses.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/355960/fpls-09-01444-HTML/image_m/fpls-09-01444-g001.jpg)
FIGURE 1. Distribution of the sampling sites of Quercus variabilis in China (PG: Pinggu, Beijing; LC: Lincheng, Hebei; JY: Jiyuan, Henan; NX: Neixiang, Henan; ZG: Zigui, Hubei; CL: Cili, Hunan; and CB: Chengbu, Hunan).
Biochemical Analysis
Total Soluble Sugars and Starch
The powdered material (∼0.10 g) was put into a 10 ml centrifuge tube, and 5 ml of 80% ethanol was added. The mixture was incubated at 80°C in a water bath shaker for 30 min, and then centrifuged at 4000 rpm for 5 min. The pellets were extracted two more times with 80% ethanol. Supernatants were retained, combined, and stored at -20°C for soluble sugar determinations. The ethanol-insoluble pellet was used for starch extraction. Glucose was used as a standard. Soluble sugars were determined using the anthrone method (Seifter et al., 1950). The starch concentration was measured spectrophotometrically at 620 nm using anthrone reagent, and was calculated by multiplying the glucose concentrations by the conversion factor of 0.9 (Osaki et al., 1991). The concentration of soluble sugars and starch was described on a dry matter basis (mg g-1 DW).
Total Nitrogen and Phosphorus
For the determination of tissues’ nitrogen (N) and phosphorus (P) concentrations (mg g-1 DW), finely ground material (∼50 mg) was first digested with H2SO4 and then H2O2 for further analysis. The nitrogen concentration was then measured using the Kjeldahl method (Kjeltec 2200, FOSS, Sweden), while the phosphorus concentration was determined with the molybdenum blue spectrophotometric procedure (6505 UV spectrophotometer, United Kingdom) (Page, 1982).
Environmental Data
The monthly climatic data (2014) were interpolated with the kriging method (Matheron, 1963) from 675 national weather stations around China (download from1) using GIS software (ArcGIS v10.0, Esri, United States). The soil data was obtained from the gridded Global Soil dataset (30 arc-second resolution) which was developed by Shangguan et al. (2014)2. The special values of each sampled plot were then extracted according to the geolocation information (latitude and longitude) of the plot. The climatic variables, monthly mean temperature (°C) and monthly mean precipitation (mm), were divided into two parts as one for the growing season (April to September; represented as gT mean and gPPT) and another for the dormant season (October to December; denoted as dT mean and dPPT). The soil variables, on the other hand, included total nitrogen (SN, %) and total phosphorus (SP, %) for each plot.
Statistical Analyses
All statistical analyses were conducted using R statistical software (RStudio version 1.0.1433). Shapiro-Wilk and Bartlett’s tests were used to test for normality and homogeneity of variances, respectively, and it was found that all data met the assumption for further variance analysis. We utilized a linear mixed effects model with tree age (juvenile and mature trees), sampling site (latitudinal gradient), and season (growing season: Aug; and dormant season: December) as fixed effects and individuals as a random effect (R package ‘lme4’). The response variables included soluble sugar, starch, NSC, nitrogen, and phosphorus in the current-year leaves, branches, and fine roots. The mean and standard error were given if necessary. For both the sampling seasons, a partial pairwise correlation analysis was performed to explore the correlations between the response variables with the environmental factors (e.g., climate or soil), where bio-factors (tree age and tree tissue) were treated as covariant (R package ‘psych’). The contributions of all the environmental and biological factors to variations of responsible variables were performed by the redundancy analysis (RDA). The variation partitioning was achieved by means of partial redundancy analysis (pRDA) to extract the pure effect of environment factors and the pure effect of bio-factors (R package ‘vegan’).
Results
Non-structural Carbohydrate and Its Components
The tree age significantly affected leaf soluble sugars and leaf NSC concentration, as well as branch starch, but did not impact root NSCs (Table 1). The sampling season and the interaction between sampling season and tree age significantly influenced the NSCs in the branches and fine roots (Table 1). For instance, root NSC in both juveniles and mature trees was significantly higher in December than in August (p < 0.001) (Supplementary Figure S1). However, during the growing season, root NSC in juveniles (30.34 ± 1.46 mg g-1) was significantly lower than that in mature trees (42.79 ± 3.06 mg g-1, p < 0.05), whereas it showed an opposite result (80.61 ± 4.79 mg g-1 for juveniles and 60.52 ± 2.94 mg g-1 for mature trees, p < 0.001) during the dormant season. On the other hand, the NSCs including the ratio of soluble sugars to starch varied strongly with latitude and tissue type (Figure 2 and Supplementary Figures S2–S4). The results from simple linear models indicated that the NSC of various tissues for both life stages did not show any decreasing trend with increasing latitude (Figure 2). Inversely, the NSC in the branches and roots of juveniles significantly or marginally significantly increased with the increasing latitude for both sampling seasons (Figure 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/355960/fpls-09-01444-HTML/image_m/fpls-09-01444-t001.jpg)
TABLE 1. Effects of age (juvenile and mature), latitude, sampled season (growing season: Aug, 2014; and dormant season: Dec, 2014) and their interactions on NSCs, nutrients over tissues revealed by linear mixed effects models with individuals as random factors.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/355960/fpls-09-01444-HTML/image_m/fpls-09-01444-g002.jpg)
FIGURE 2. The non-structural carbohydrate concentrations (mg/g) (mean ± se, n = 3) across life stages (juveniles vs. mature), tissues (leaf, branch, and fine root), and sampling seasons (growing season: Aug, 2014; and dormant season: Dec, 2014) along the latitudinal gradient. For each subplot, blue color denotes mature individuals and red denotes juveniles; fitted curves, determinant coefficients, and p-values of simple linear regression are given.
Nitrogen, Phosphorus Concentrations and Their Ratios
The tree age had significant impacts on N, P, and N:P ratio of the studied tissues with exception of leaf N and branch N (Table 1 and Supplementary Figure S5). Mature trees generally had higher P concentration across tissues than juveniles, especially during the dormant season (Supplementary Figure S5). During the growing season, the N, P, and N:P ratio in the leaves of mature trees were 18.93 ± 0.32 mg g-1, 1.63 ± 0.18 mg g-1, and 14.27 ± 1.48, respectively, whereas, the corresponding values for juveniles were 18.54 ± 0.36 mg g-1, 1.13 ± 0.07 mg g-1, and 17.71 ± 1.08, respectively (Supplementary Figure S5). The P and N:P ratio in various tissues were interactively affected by tree age and latitude (Table 1), and N:P ratio in the leaves and roots of mature trees during the growing season increased significantly with latitude (Supplementary Figure S6). The tissues’ N varied significantly with latitude (Table 1 and Figure 3), while root N of both life stages increased with increasing latitude for both sampling seasons with the exception of mature trees’ roots in the dormant season (Figure 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/355960/fpls-09-01444-HTML/image_m/fpls-09-01444-g003.jpg)
FIGURE 3. The total nitrogen and phosphorus concentrations (mg/g) (mean ± se, n = 3) across life stages (juveniles vs. mature), tissues (leaf, branch, and fine root), and sampling seasons (growing season: Aug, 2014; and dormant season: Dec, 2014) along the latitudinal gradient. For each subplot, blue color denotes mature individuals and red denotes juveniles; fitted curves, determinant coefficients, and p-values of simple linear regression are given.
Correlation Between Non-carbohydrates, Nutrients, and Environmental Factors
The season transition strongly affected the intensity of correlations between dependent variables with environmental variables (Figures 4, 5). During the growing season, tissues’ NSC and N have no significant correlations with the environmental factors, e.g., geographical, soil, or climatic variables (Figure 4A). The ratios of soluble sugar to starch were negatively correlated with precipitation of growing season. The ratios of N to P were positively correlated with latitude and longitude but negatively correlated with elevation, soil N, precipitation, and mean temperature of growing season (Figure 4A). The partial RDA revealed that only a total of 26.2% variations of dependent variables were jointly explained by environmental factors (17.5%) and bio-factors (8.3%) (Supplementary Figure S7A). Nevertheless, during the dormant season, soluble sugar, NSC, the ratio of soluble sugar to starch and tissue nitrogen were correlated positively with latitude and longitude but negatively with elevation, soil N, soil P, precipitation, and mean temperature of dormant season (Figure 4B). A total of 48.2% variations of dependent variables were jointly explained by environmental (27.5%) and bio-factors (20.6%) (Supplementary Figure S7B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/355960/fpls-09-01444-HTML/image_m/fpls-09-01444-g004.jpg)
FIGURE 4. Partial correlation analysis for dependent variables (soluble sugar, starch, NSC, ratio of sugar to starch, N, P, and ratio of N to P) and environmental factors (latitude, longitude, elevation, soil N, soil P, precipitation, and monthly mean temperature) where tree age and tree tissue were regarded as control factors [(A): growing season; (B): dormant season]. Significantly positive or negative correlations (p < 0.05) are displayed in blue and in red color, respectively; and color intensity and circle size are proportional to the correlation coefficients.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/355960/fpls-09-01444-HTML/image_m/fpls-09-01444-g005.jpg)
FIGURE 5. Results of redundancy analysis [(A): growing season; (B): dormant season] showing samples (mature tree: dark blue points; juvenile tree: red square), responsive variables (black triangle), and affecting factors (blue arrow).
Discussion
The present study showed that NSCs (total soluble sugars, starch, and NSC) in woody tissues (branches and fine roots) of both life stages had higher concentrations in the dormant season than in the growing season along the latitudinal gradient (Supplementary Figure S1). Our results were consistent with Martínez-Vilalta et al. (2016) who complied data from 121 studies including 177 species under natural conditions and found that NSCs varied seasonally, with a general increase during winter months. Compared to evergreen species, deciduous trees generally stored more NSC in the tissues (e.g., stem or root) to withstand low temperature in the coming winter and to support bud-break and shoot growth in the early spring (Klein et al., 2016). On the other hand, higher level of NSC in reserved tissues during the dormant season may attribute to lower level of growth and maintenance respiration due to low temperature in the winter, with more starch broken down into soluble sugars to promote cold tolerance by adjusting the intracellular osmotic concentration (Morin et al., 2007), as shown by an increase in the sugar-starch ratios in juvenile roots in our study (Supplementary Figure S4). Juveniles’ roots had lower NSC concentration in the growing season but higher NSC concentration in the dormant season than that of mature individuals, especially for north populations (Figure 2), indicating that juveniles are more sensitive to seasonal transition than mature trees. This further suggests that juveniles may use a strategy to invest more NSC into fast growth (thus leading to lower level of NSC in juveniles than in mature trees during the growing season) to get a competitive advantage with relatively larger body size in a community. During the dormant season, juveniles generally had a higher level of NSC to cope with low temperature because of their relatively younger tissue and smaller body size (i.e., smaller carbon pool size) (Bansal and Germino, 2009, 2010). On the other hand, Li et al. (2008b) proposed that trees growing at the elevational or the latitudinal climate limit rely not only on the total NSC concentration but also require a sufficiently high sugar-starch ratio to overwinter successfully. Our results support this view, as both life stages maintained comparable high ratio of soluble sugars to starch in the dormant season (Supplementary Figure S1).
Although there exists considerable disparities between altitudinal and latitudinal gradients (Jump et al., 2009), both gradients are mainly dominated by temperature. In this regard, the altitudinal tree-line hypotheses (e.g., GLH and CLH) present important reference values to deduce the driven mechanism for latitudinal range. In the present study, we found that tissue NSC did not exhibit decreasing trends for both life stages along the latitudinal gradient. Nevertheless, NSC in the woody tissues of juveniles increased significantly with increasing latitude during both the sampling seasons, suggesting that carbon limitation cannot act as a determinant driver to north latitudinal range of the species studied. This result fails to support the CLH along latitudinal gradients. Li et al. (2016) also observed weekly increasing latitudinal trends of leaf NSCs from tropical to cold temperate forests at the levels of species and plant functional groups. Alternatively, other factors, such as recruitment limitation, rather than NSC at the north edges may contribute to the oak’s north-latitudinal limit formation. Our previous study found that the relative densities of seedlings of the species studied were significantly lower in the northern edge than in the core populations (Gao et al., 2017).
Apart from the role of carbon status in determining the range limit, nutrient shortage is another important factor to elucidate the issue which not only limited carbon assimilation but also limited tissue development (McNown and Sullivan, 2013; Sullivan et al., 2014). Ontogenetic variations in tissue P but not in tissue N were found in this study (Supplementary Figure S5), where mature trees had higher tissue P concentration than juveniles, especially during the dormant season, consequently leading to relatively lower N:P ratio occurring in mature trees. This result is inconsistent with Noh et al. (2007) who found that N and P concentrations in tissues of Q. acutissima significantly decreased with tree age or size. Indeed, plant nutrient demand and morphological structure (e.g., root morphology), which vary with ontogeny (Álvarez-Yépiz et al., 2014), may contribute to the discrepancies mentioned above. For example, mature trees featured with larger and deeper root system could help to uptake more available soil P which mainly originates from rock weathering, while available soil N mostly comes from atmospheric deposition in nature ecosystems (Verhoeven and Schmitz, 1991; Güsewell, 2004).
Leaf nutrient concentrations were generally closely correlated with soil available nutrients (Koerselman and Meuleman, 1996; Tessier and Raynal, 2003; Ordoñez et al., 2009; Wang et al., 2017). The mean value of the leaf N:P ratio in the present study was 14.3 ± 1.48 for mature trees and 17.7 ± 1.08 for juveniles for the growing season (Supplementary Figure S5), which may indicate that the juveniles are slightly limited by P according to the threshold for P-limitation established by Koerselman and Meuleman (1996), but this case is still not limited by P availability according to the relaxed threshold suggested by Güsewell (2004). Sun et al. (2015) found that the mean N, P concentration and the N:P ratio of Chinese cork oak leaves were 19.00 ± 0.26 mg g-1, 1.03 ± 0.03 mg g-1, and 20.48 ± 0.63, respectively, across the distribution range of that species in China. Hence, leaf N:P ratio in the present study was lower than that of Sun et al. (2015) but similar to that of Wu et al. (2012) (16.56 for Q. variabilis). In line with our results of P concentration (1.13 ± 0.07 mg g-1 for juvenile, 1.63 ± 0.18 mg g-1 for mature), Han et al. (2005) pointed out that low leaf P across China’s flora, compared to the global average revealed by Reich and Oleksyn (2004) (P: 1.77 mg g-1), is a result of low soil P content in China. On the other hand, we found that the N:P ratio in leaves and roots of mature trees increased with latitude (Supplementary Figure S6), which is a result of increased tissue N accompanied by relatively stable P concentration across latitudes. A higher leaf N concentration in plants in colder habitats is generally considered as an adaptation mechanism that enhances the metabolic activity and growth rates under low temperatures or short growing season (Reich and Oleksyn, 2004; Soolanayakanahally et al., 2009; El Zein et al., 2011). Fajardo and Piper (2017) also found that nutrient limitation is not likely to be involved in the carbon limitations and could not be an explanation for altitudinal tree-line formation. Our recent field investigation found that the Chinese cork oak showed increasing annual basal increment with increasing latitude from south to north (Gao et al., 2018). We therefore speculate that higher availability of resources (e.g., NSC, N, and P) found in northern populations supports higher growth rate, which in turn decreases the resource storage and remains the resource availability at a stable level.
The fact that more pronounced negative correlations between resource storage and climatic variables (e.g., mean temperature) were found in the dormant season (Figure 4B) implies that the projected climate warming could alter trajectories in plant resource allocation into growth and regeneration (De Frenne et al., 2012; Lapenis et al., 2013; Carón et al., 2015), especially as the climate warming in China is predicted to be strongest during the winter months (Ge et al., 2013). This in turn may result in distribution range shift. The results presented here are based on one species with two sampling dates, and should be supplemented by further studies carried out at the community or ecosystem level with multiple species and spanning multiple years, to take into account species-specific responses to environmental variations (Drobyshev et al., 2013; Boisvert-Marsh et al., 2014) as well as their inter-annual fluctuations (Scartazza et al., 2013). Meanwhile, the other factors such as recruitment potential (Graignic et al., 2014), competitive interaction (Ettinger and HilleRisLambers, 2013; Liang et al., 2016), and disturbance (Slaton, 2015; Renwick et al., 2016) had been found to affect distribution range, which should be considered in further exploring the mechanism determining tree species’ northern distribution range other than resource availability.
Conclusion
In the present study, from the perspective of carbon and nutrient allocation strategy, a whole-tree approach integrating with two distinct life stages was used to explore the mechanisms that determine tree latitudinal range limit, with a case study of Chinese cork oak. Our results indicate that tissue NSCs (soluble sugar, starch and sum of soluble sugar and starch, NSC) and nutrient concentrations (nitrogen and phosphorus) in both mature and juvenile trees did not decrease with increasing latitude across seasons, and even root nitrogen and root NSC in the juvenile trees increased with latitude. Our results suggest that available carbon, N, and P may not be the determinant factors driving the latitudinal range limit of the species studied. These findings will greatly improve our understanding of the mechanism involved in determining the latitudinal range limit, and help to understand and predict the dynamics of the northern range under global warming.
Author Contributions
J-FL, M-HL, and J-PL conceived the experiment. Y-PD, X-FW, and Y-YN conducted the experiment and analyzed the samples. J-FL, M-HL, QW, W-FX, and Z-PJ wrote the paper. All authors contributed critically to the drafts and gave final approval for publication.
Funding
This work was supported by the National Natural Science Foundation of China (Nos. 41371075 and 41371076), the Fundamental Research Funds for the Central Non-profit Research Institution of CAF (CAFYBB2018ZB001 and CAFYBB2014ZD001), the Natural Science Foundation of Beijing (No. 8152032), and the China Scholarship Council (No. 201303270003).
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.01444/full#supplementary-material
Footnotes
References
Álvarez-Yépiz, J. C., Cueva, A., Dovčiak, M., Teece, M., and Yepez, E. A. (2014). Ontogenetic resource-use strategies in a rare long-lived cycad along environmental gradients. Conserv. Physiol. 2:cou034. doi: 10.1093/conphys/cou034
Bansal, S., and Germino, M. J. (2009). Temporal variation of nonstructural carbohydrates in montane conifers: similarities and differences among developmental stages, species and environmental conditions. Tree Physiol. 29, 559–568. doi: 10.1093/treephys/tpn045
Bansal, S., and Germino, M. J. (2010). Variation in ecophysiological properties among conifers at an ecotonal boundary: comparison of establishing seedlings and established adults at timberline. J. Veg. Sci. 21, 133–142. doi: 10.1111/j.1654-1103.2009.01127.x
Bazot, S., Barthes, L., Blanot, D., and Fresneau, C. (2013). Distribution of non-structural nitrogen and carbohydrate compounds in mature oak trees in a temperate forest at four key phenological stages. Trees 27, 1023–1034. doi: 10.1007/s00468-013-0853-5
Bell, D. M., Bradford, J. B., and Lauenroth, W. K. (2014). Early indicators of change: divergent climate envelopes between tree life stages imply range shifts in the western United States. Glob. Ecol. Biogeogr. 23, 168–180. doi: 10.1111/geb.12109
Boisvert-Marsh, L., Périé, C., and de Blois, S. (2014). Shifting with climate? Evidence for recent changes in tree species distribution at high latitudes. Ecosphere 5:art83. doi: 10.1890/ES14-00111.1
Carón, M. M., De Frenne, P., Brunet, J., Chabrerie, O., Cousins, S. A. O., Decocq, G., et al. (2015). Divergent regeneration responses of two closely related tree species to direct abiotic and indirect biotic effects of climate change. For. Ecol. Manag. 342, 21–29. doi: 10.1016/j.foreco.2015.01.003
Cavender-Bares, J., and Bazzaz, F. A. (2000). Changes in drought response strategies with ontogeny in Quercus rubra: implications for scaling from seedlings to mature trees. Oecologia 124, 8–18. doi: 10.1007/PL00008865
Chen, D., Zhang, X., Kang, H., Sun, X., Yin, S., Du, H., et al. (2012). Phylogeography of Quercus variabilis based on chloroplast DNA sequence in East Asia: multiple glacial refugia and Mainland-migrated island populations. PLoS One 7:e47268. doi: 10.1371/journal.pone.0047268
Chen, I. C., Hill, J. K., Ohlemuller, R., Roy, D. B., and Thomas, C. D. (2011). Rapid range shifts of species associated with high levels of climate warming. Science 333, 1024–1026. doi: 10.1126/science.1206432
De Frenne, P., Graae, B. J., Brunet, J., Shevtsova, A., De Schrijver, A., Chabrerie, O., et al. (2012). The response of forest plant regeneration to temperature variation along a latitudinal gradient. Ann. Bot. 109, 1037–1046. doi: 10.1093/aob/mcs015
De Frenne, P., Graae, B. J., Rodríguez-Sánchez, F., Kolb, A., Chabrerie, O., Decocq, G., et al. (2013). Latitudinal gradients as natural laboratories to infer species’ responses to temperature. J. Ecol. 101, 784–795. doi: 10.1111/1365-2745.12074
DeHayes, D. H., Ingle, M. A., and Waite, C. E. (1989). Nitrogen fertilization enhances cold tolerance of red spruce seedlings. Can. J. For. Res. 19, 1037–1043. doi: 10.1139/x89-158
Desprez, J., Iannone, B. V., Yang, P., Oswalt, C. M., and Fei, S. (2014). Northward migration under a changing climate: a case study of blackgum (Nyssa sylvatica). Clim. Chang. 126, 151–162. doi: 10.1007/s10584-014-1207-z
Drobyshev, I., Gewehr, S., Berninger, F., and Bergeron, Y. (2013). Species specific growth responses of black spruce and trembling aspen may enhance resilience of boreal forest to climate change. J. Ecol. 101, 231–242. doi: 10.1111/1365-2745.12007
El Zein, R., Maillard, P., Bréda, N., Marchand, J., Montpied, P., and Gérant, D. (2011). Seasonal changes of C and N non-structural compounds in the stem sapwood of adult sessile oak and beech trees. Tree Physiol. 31, 843–854. doi: 10.1093/treephys/tpr074
Ettinger, A. K., and HilleRisLambers, J. (2013). Climate isn’t everything: competitive interactions and variation by life stage will also affect range shifts in a warming world. Am. J. Bot. 100, 1344–1355. doi: 10.3732/ajb.1200489
Fajardo, A., and Piper, F. I. (2017). An assessment of carbon and nutrient limitations in the formation of the southern Andes tree line. J. Ecol. 105, 517–527. doi: 10.1111/1365-2745.12697
Gao, W.-Q., Liu, J.-F., Xue, Z.-M., Zhang, Y.-T., Gao, Z.-H., Ni, Y.-Y., et al. (2018). Geographical patterns and drivers of growth dynamics of Quercus variabilis. For. Ecol. Manag. 429, 256–266. doi: 10.1016/j.foreco.2018.07.024
Gao, W.-Q., Ni, Y.-Y., Xue, Z.-M., Wang, X.-F., Kang, F.-F., Hu, J., et al. (2017). Population structure and regeneration dynamics of Quercus variabilis along latitudinal and longitudinal gradients. Ecosphere 8:e01737. doi: 10.1002/ecs2.1737
Ge, Q., Wang, F., and Luterbacher, J. (2013). Improved estimation of average warming trend of China from 1951–2010 based on satellite observed land-use data. Clim. Chang. 121, 365–379. doi: 10.1007/s10584-013-0867-4
Gedroc, J. J., McConnaughay, K. D. M., and Coleman, J. S. (1996). Plasticity in root/shoot partitioning: optimal, ontogenetic, or both? Funct. Ecol. 10, 44–50. doi: 10.2307/2390260
Gong, G., and Jian, W. (1983). On the geographical distribution of phenodate in China. Acta Geogr. Sin. 38, 33–40. doi: 10.1021/es201859j
Graignic, N., Tremblay, F., and Bergeron, Y. (2014). Geographical variation in reproductive capacity of sugar maple (Acer saccharum Marshall) northern peripheral populations. J. Biogeogr. 41, 145–157. doi: 10.1111/jbi.12187
Güsewell, S. (2004). N : p ratios in terrestrial plants: variation and functional significance. New Phytol. 164, 243–266. doi: 10.1111/j.1469-8137.2004.01192.x
Han, W., Fang, J., Guo, D., and Zhang, Y. (2005). Leaf nitrogen and phosphorus stoichiometry across 753 terrestrial plant species in China. New Phytol. 168, 377–385. doi: 10.1111/j.1469-8137.2005.01530.x
Hoch, G., and Körner, C. (2003). The carbon charging of pines at the climatic treeline: a global comparison. Oecologia 135, 10–21. doi: 10.1007/s00442-002-1154-7
Hoch, G., and Körner, C. (2012). Global patterns of mobile carbon stores in trees at the high-elevation tree line. Glob. Ecol. Biogeogr. 21, 861–871. doi: 10.1111/j.1466-8238.2011.00731.x
Jump, A. S., Mátyás, C., and Peñuelas, J. (2009). The altitude-for-latitude disparity in the range retractions of woody species. Trends Ecol. Evol. 24, 694–701. doi: 10.1016/j.tree.2009.06.007
Kelly, A. E., and Goulden, M. L. (2008). Rapid shifts in plant distribution with recent climate change. Proc. Natl. Acad. Sci. U.S.A. 105, 11823–11826. doi: 10.1073/pnas.0802891105
Klein, T., Vitasse, Y., and Hoch, G. (2016). Coordination between growth, phenology and carbon storage in three coexisting deciduous tree species in a temperate forest. Tree Physiol. 36, 847–855. doi: 10.1093/treephys/tpw030
Klockmann, M., Günter, F., and Fischer, K. (2017). Heat resistance throughout ontogeny: body size constrains thermal tolerance. Glob. Chang. Biol. 23, 686–696. doi: 10.1111/gcb.13407
Koerselman, W., and Meuleman, A. F. M. (1996). The vegetation N:P ratio: a new tool to detect the nature of nutrient limitation. J. Appl. Ecol. 33, 1441–1450. doi: 10.2307/2404783
Körner, C. (1998). A re-assessment of high elevation treeline positions and their explanation. Oecologia 115, 445–459. doi: 10.1007/s004420050540
Körner, C. (2003). Carbon limitation in trees. J. Ecol. 91, 4–17. doi: 10.1046/j.1365-2745.2003.00742.x
Körner, C. (2007). The use of ’altitude’ in ecological research. Trends Ecol. Evol. 22, 569–574. doi: 10.1016/j.tree.2007.09.006
Lapenis, A. G., Lawrence, G. B., Heim, A., Zheng, C., and Shortle, W. (2013). Climate warming shifts carbon allocation from stemwood to roots in calcium-depleted spruce forests. Glob. Biogeochem. Cycles 27, 101–107. doi: 10.1029/2011gb004268
le Roux, P. C., Shaw, J. D., and Chown, S. L. (2013). Ontogenetic shifts in plant interactions vary with environmental severity and affect population structure. New Phytol. 200, 241–250. doi: 10.1111/nph.12349
Li, M. H., Hoch, G., and Körner, C. (2002). Source/sink removal affects mobile carbohydrates in Pinus cembra at the Swiss treeline. Trees 16, 331–337. doi: 10.1007/s00468-002-0172-8
Li, M.-H., and Krauchi, N. (2005). The state of knowledge on alpine treeline and suggestions for future research. J. Sichuan For. Sci. Technol. 26, 36–42.
Li, M.-H., Xiao, W.-F., Shi, P., Wang, S.-G., Zhong, Y.-D., Liu, X.-L., et al. (2008a). Nitrogen and carbon source–sink relationships in trees at the Himalayan treelines compared with lower elevations. Plant Cell Environ. 31, 1377–1387. doi: 10.1111/j.1365-3040.2008.01848.x
Li, M.-H., Xiao, W.-F., Wang, S.-G., Cheng, G.-W., Cherubini, P., Cai, X.-H., et al. (2008b). Mobile carbohydrates in Himalayan treeline trees I. Evidence for carbon gain limitation but not for growth limitation. Tree Physiol. 28, 1287–1296. doi: 10.1093/treephys/28.8.1287
Li, N., He, N., Yu, G., Wang, Q., and Sun, J. (2016). Leaf non-structural carbohydrates regulated by plant functional groups and climate: evidences from a tropical to cold-temperate forest transect. Ecol. Indic. 62, 22–31. doi: 10.1016/j.ecolind.2015.11.017
Liang, E., Wang, Y., Piao, S., Lu, X., Camarero, J. J., Zhu, H., et al. (2016). Species interactions slow warming-induced upward shifts of treelines on the Tibetan Plateau. Proc. Natl. Acad. Sci. U.S.A. 113, 4380–4385. doi: 10.1073/pnas.1520582113
Máliš, F., Kopecký, M., Petřík, P., Vladovič, J., Merganič, J., and Vida, T. (2016). Life stage, not climate change, explains observed tree range shifts. Glob. Chang. Biol. 22, 1904–1914. doi: 10.1111/gcb.13210
Martínez-Vilalta, J., Sala, A., Asensio, D., Galiano, L., Hoch, G., Palacio, S., et al. (2016). Dynamics of non-structural carbohydrates in terrestrial plants: a global synthesis. Ecol. Monogr. 86, 495–516. doi: 10.1002/ecm.1231
Matheron, G. (1963). Principles of geostatistics. Econ. Geol. 58, 1246–1266. doi: 10.2113/gsecongeo.58.8.1246
McNown, R. W., and Sullivan, P. F. (2013). Low photosynthesis of treeline white spruce is associated with limited soil nitrogen availability in the Western Brooks Range, Alaska. Funct. Ecol. 27, 672–683. doi: 10.1111/1365-2435.12082
Morin, X., Améglio, T., Ahas, R., Kurz-Besson, C., Lanta, V., Lebourgeois, F., et al. (2007). Variation in cold hardiness and carbohydrate concentration from dormancy induction to bud burst among provenances of three European oak species. Tree Physiol. 27, 817–825. doi: 10.1093/treephys/27.6.817
Niinemets, Ü. (2010). Responses of forest trees to single and multiple environmental stresses from seedlings to mature plants: past stress history, stress interactions, tolerance and acclimation. For. Ecol. Manag. 260, 1623–1639. doi: 10.1016/j.foreco.2010.07.054
Noh, N. J., Son, Y., Kim, R. H., Seo, K. W., Koo, J. W., Park, I. H., et al. (2007). Biomass accumulations and the distribution of nitrogen and phosphorus within three Quercus acutissima stands in central Korea. J. Plant Biol. 50, 461–466. doi: 10.1007/bf03030683
O’Brien, M. J., Leuzinger, S., Philipson, C. D., Tay, J., and Hector, A. (2014). Drought survival of tropical tree seedlings enhanced by non-structural carbohydrate levels. Nat. Clim. Chang. 4, 710–714. doi: 10.1038/nclimate2281
Ordoñez, J. C., Van Bodegom, P. M., Witte, J.-P. M., Wright, I. J., Reich, P. B., and Aerts, R. (2009). A global study of relationships between leaf traits, climate and soil measures of nutrient fertility. Glob. Ecol. Biogeogr. 18, 137–149. doi: 10.1111/j.1466-8238.2008.00441.x
Osaki, M., Shinano, T., and Tadano, T. (1991). Redistribution of carbon and nitrogen compounds from the shoot to the harvesting organs during maturation in field crops. Soil Sci. Plant Nutr. 37, 117–128. doi: 10.1080/00380768.1991.10415017
Page, A. L. (1982). “Methods of soil analysis. Part 2. Chemical and microbiological properties,” in American Society of Agronomy, Soil Science Society of America, eds A. Klute and A. L. Page (Oakland: The University of California).
Portsmuth, A., Niinemets,Ü., Truus, L., and Pensa, M. (2005). Biomass allocation and growth rates in Pinus sylvestris are interactively modified by nitrogen and phosphorus availabilities and by tree size and age. Can. J. For. Res. 35, 2346–2359. doi: 10.1139/x05-155
Reich, P. B., and Oleksyn, J. (2004). Global patterns of plant leaf N and P in relation to temperature and latitude. Proc. Natl. Acad. Sci. U.S.A. 101, 11001–11006. doi: 10.1073/pnas.0403588101
Renwick, K. M., Rocca, M. E., and Stohlgren, T. J. (2016). Biotic disturbance facilitates range shift at the trailing but not the leading edge of lodgepole pine’s altitudinal distribution. J. Veg. Sci. 27, 780–788. doi: 10.1111/jvs.12410
Rinas, C. L., Dial, R. J., Sullivan, P. F., Smeltz, T. S., Tobin, S. C., Loso, M., et al. (2017). Thermal segregation drives patterns of alder and willow expansion in a montane ecosystem subject to climate warming. J. Ecol. 105, 935–946. doi: 10.1111/1365-2745.12737
Rivas-Ubach, A., Sardans, J., Pérez-Trujillo, M., Estiarte, M., and Peñuelas, J. (2012). Strong relationship between elemental stoichiometry and metabolome in plants. Proc. Natl. Acad. Sci. U.S.A. 109, 4181–4186. doi: 10.1073/pnas.1116092109
Scartazza, A., Moscatello, S., Matteucci, G., Battistelli, A., and Brugnoli, E. (2013). Seasonal and inter-annual dynamics of growth, non-structural carbohydrates and C stable isotopes in a Mediterranean beech forest. Tree Physiol. 33, 730–742. doi: 10.1093/treephys/tpt045
Seifter, S., Dayton, S., Novic, B., and Muntwyler, E. (1950). The estimation of glycogen with the anthrone reagent. Arch. Biochem. 25, 191–200.
Shangguan, W., Dai, Y., Duan, Q., Liu, B., and Yuan, H. (2014). A global soil data set for earth system modeling. J. Adv. Model. Earth Syst. 6, 249–263. doi: 10.1002/2013MS000293
Sittaro, F., Paquette, A., Messier, C., and Nock, C. A. (2017). Tree range expansion in eastern North America fails to keep pace with climate warming at northern range limits. Glob. Chang. Biol. 23, 3292–3301. doi: 10.1111/gcb.13622
Slaton, M. R. (2015). The roles of disturbance, topography and climate in determining the leading and rear edges of population range limits. J. Biogeogr. 42, 255–266. doi: 10.1111/jbi.12406
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
Steppe, K., Niinemets,Ü., and Teskey, R. O. (2011). “Tree size- and age-related changes in leaf physiology and their influence on carbon gain,” in Size- and Age-Related Changes in Tree Structure and Function, eds F. C. Meinzer, B. Lachenbruch, and T. E. Dawson. (Dordrecht: Springer), 235–253.
Stevens, G. C., and Fox, J. F. (1991). The causes of treeline. Annu. Rev. Ecol. Syst. 22, 177–191. doi: 10.2307/2097259
Sullivan, P. F., Ellison, S. B. Z., McNown, R. W., Brownlee, A. H., and Sveinbjörnsson, B. (2014). Evidence of soil nutrient availability as the proximate constraint on growth of treeline trees in northwest Alaska. Ecology 96, 716–727. doi: 10.1890/14-0626.1
Sun, X., Kang, H., Kattge, J., Gao, Y., and Liu, C. (2015). Biogeographic patterns of multi-element stoichiometry of Quercus variabilis leaves across China. Can. J. For. Res. 45, 1827–1834. doi: 10.1139/cjfr-2015-0110
Tessier, J. T., and Raynal, D. J. (2003). Use of nitrogen to phosphorus ratios in plant tissue as an indicator of nutrient limitation and nitrogen saturation. J. Appl. Ecol. 40, 523–534. doi: 10.1046/j.1365-2664.2003.00820.x
Verhoeven, J. T. A., and Schmitz, M. B. (1991). Control of plant growth by nitrogen and phosphorus in mesotrophic fens. Biogeochemistry 12, 135–148. doi: 10.1007/bf00001811
Villagra, M., Campanello, P. I., Montti, L., and Goldstein, G. (2013). Removal of nutrient limitations in forest gaps enhances growth rate and resistance to cavitation in subtropical canopy tree species differing in shade tolerance. Tree Physiol. 33, 285–296. doi: 10.1093/treephys/tpt003
Wang, X., Xu, Z., Yan, C., Luo, W., Wang, R., Han, X., et al. (2017). Responses and sensitivity of N, P and mobile carbohydrates of dominant species to increased water, N and P availability in semi-arid grasslands in northern China. J. Plant Ecol. 10, 486–496. doi: 10.1093/jpe/rtw053
Wu, T., Dong, Y., Yu, M., Geoff Wang, G., and Zeng, D.-H. (2012). Leaf nitrogen and phosphorus stoichiometry of Quercus species across China. For. Ecol. Manag. 284, 116–123. doi: 10.1016/j.foreco.2012.07.025
Yan, B., Ji, Z., Fan, B., Wang, X., He, G., Shi, L., et al. (2016). Plants adapted to nutrient limitation allocate less biomass into stems in an arid-hot grassland. New Phytol. 211, 1232–1240. doi: 10.1111/nph.13970
Keywords: Quercus variabilis, latitudinal distribution range, ontogeny, non-structural carbohydrate, nutrient
Citation: Liu J-F, Deng Y-P, Wang X-F, Ni Y-Y, Wang Q, Xiao W-F, Lei J-P, Jiang Z-P and Li M-H (2018) The Concentration of Non-structural Carbohydrates, N, and P in Quercus variabilis Does Not Decline Toward Its Northernmost Distribution Range Along a 1500 km Transect in China. Front. Plant Sci. 9:1444. doi: 10.3389/fpls.2018.01444
Received: 26 January 2018; Accepted: 12 September 2018;
Published: 17 October 2018.
Edited by:
Hans J. De Boeck, University of Antwerp, BelgiumReviewed by:
Ivika Ostonen, University of Tartu, EstoniaEryuan Liang, Institute of Tibetan Plateau Research (CAS), China
Copyright © 2018 Liu, Deng, Wang, Ni, Wang, Xiao, Lei, Jiang and Li. 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) and the copyright owner(s) 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: Jian-Feng Liu, Liujf2000cn@163.com Jing-Pin Lei, leijp@caf.ac.cn Ze-Ping Jian, Jiangzp@caf.ac.cn Mai-He Li, Maihe.li@wsl.ch