- Key Laboratory of Bio-Resource and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, China
With global climate change, for evaluating warming effect on subalpine forest distribution, the substantial effects of long-term warming on tree growth and soil nutrients need to be explored. In this study, we focused on different responses in the boundaries of trees and soils to warming. Using the open-top chamber (OTC), a 10-year artificial warming experiment was conducted to evaluate the impacts of warming on Abies faxoniana at three different altitudes. We determined metabolites and nutrient concentrations in needles of A. faxoniana and characterized the soil chemistries. Many kinds of sugars, amino acids, and organic acids showed higher contents at high altitude (3,500 m) compared with low altitude (2,600 m), which could have been due to the temperature differences. Warming significantly decreased needle sugar and amino acid concentrations at high altitude but increased them at low altitude. These results indicated contrasting physiological and metabolic responses of A. faxoniana to long-term warming at different altitudes. Furthermore, we found that OTC warming significantly increased the concentrations of soil extractable sodium, aluminum (Al), and manganese (Mn), while decreased potassium (K) and phosphorus (P) concentrations and pH values at low altitude rather than at middle (3,000 m) or high altitude. The soil carbon and nitrogen contents were increased only at the middle altitude. In A. faxoniana at low altitudes, more mineral nutrients iron, K, and P were demand, and a mass of Al, Mn, and zinc was accumulated under warming. Soil P limitation and heavy metals accumulation are disadvantageous for trees at low altitudes with warming. Therefore, compared with high altitudes, A. faxoniana growing at low boundary in alpine regions is expected to be more susceptible to warming.
Introduction
Global warming has attracted much attention around the world. The Tibetan Plateau is experiencing a higher warming rate than elsewhere in the Northern Hemisphere (Duan and Xiao, 2015; Guo et al., 2018), where a 0.32°C increase per decade has occurred over the past 50 years (Luo et al., 2013; Wei and Fang, 2013). Mountain forests are of particular importance in maintaining regional ecosystem stability and succession in China. Temperature acts as a key factor in regulating plant physiology (Rustad et al., 2001; Ziello et al., 2009). Altitudinal gradient is closely related to plant physiology because of a vertical temperature range (Körner, 2007). The rate of warming in mountain forests is amplified with elevation (Pepin et al., 2015). Alpine forests might be especially susceptible to warming, and the sensitivity would vary with elevation (Wang et al., 2017; Panthi et al., 2020). There are many studies on the sink-source balance by exploring alterations in plant growth, photosynthesis, and respiration under experimental warming (Yin et al., 2008; Volder et al., 2013; Harvey et al., 2020; Carter et al., 2021). However, few studies have focused on leaf phenotypic plasticity and nutrient traits combined with soil properties, in long-term field warming experiments.
Many plant physiological processes and traits can be affected by temperature alteration, such as biosynthesis of carbohydrates and leaf nutrient accumulation (Aspinwall et al., 2016). Temperature alteration can damage in the plant cells by disturbing the balance in cellular respiration. Soluble carbohydrates and amino acids play important roles in protecting plant photosynthetic apparatus to maintain photosynthetic capacity, cellular osmotic homeostasis, and membrane stability. A study on Solanum tuberosum revealed that increased temperature could increase the concentration of many primary metabolites including amino acids, fatty acids, fatty alcohols, and organic acids (Hancock et al., 2014). In addition, warming also induced alterations in the nutrient availability in soil and plant nutrient absorption capacity in soils and plants (D’Orangeville et al., 2014; Wu et al., 2020). Previous studies have reported that the effects of warming on soil properties are closely related to tree physiology in subalpine regions (Yin et al., 2013; Zhao et al., 2014). Short-term field warming could change soil microbial carbon (C), nitrogen (N), and the microbial community structure, and there were different physiological responses on the growth of Abies faxoniana seedlings (Wu et al., 2011; Xu et al., 2013, 2015). However, the effect of soil change on tree physiology responding to long-term warming is still unknown.
Previous studies have found that A. faxoniana seedlings from higher altitude have higher leaf N status, net photosynthetic rate, and greater non-structural carbohydrate concentrations than those from lower altitude (Wu et al., 2011; Ran et al., 2013). Phenolic compounds and mineral contents are also varied across altitudes under warming (Wu et al., 2011; Ran et al., 2013; Stromme et al., 2018; Mason et al., 2019). Trees at high elevation are expected to have local adaptation to high frequency and intensity of cold events. The concentrations of sugars, organic acids, and amines would be increased more strongly in high than low-altitude populations (Schrieber et al., 2020). However, tree species in the low treeline ecotones have been found to be more sensitive to climatic factors (Takahashi et al., 2003; Kueppers et al., 2017). Evidence have shown that tree growth declines and mortality increases at lower altitudes as temperature increases (Huang et al., 2017; Wang et al., 2017), whereas at higher altitudes, tree growth and tree density increase synchronously (Liang et al., 2011; Dang et al., 2012; Ponocna et al., 2016; Conlisk et al., 2017).
Temperature is a dominant constraint on tree growth in the high elevations on the eastern Tibetan Plateau (Wang et al., 2017). The growth of firs in subalpine regions in the eastern Tibetan Plateau is sensitive to changed temperature regime rather than the variation of precipitation pattern (Wang et al., 2017). The plant metabolome is the chemical phenotype of individual plants or populations under specific environments (Rivas-Ubach et al., 2012). Additionally, for long-living trees, genetic adaptation is slow (Savolainen et al., 2004). Therefore, they are likely to cope with environmental changes more rapidly through plastic responses than through microevolution (Chevin et al., 2013; Franks et al., 2014). Trees are likely to have unique metabolic profiles due to long-term warming at different altitudes, particularly in perennial forests (Dhuli et al., 2014). However, how much phenotypic plasticity can help them cope with climate change remains largely unknown (Anderson et al., 2012; Duputie et al., 2015). It is important to understand how tree physiology and growth will respond to global changes on a wide scale.
In the present study, subalpine forest systems present unique opportunities to study the warming effects across different altitudes because they are characterized by different temperatures. To explore tree physiological responses to warming, the metabolome and ionome analyses were performed to investigate the responses of trees to a 10 year long-term field warming. We have performed artificial warming experiments using open-top chambers (OTCs) with A. faxoniana in this region since 2007 (Wu et al., 2011). We will address the following hypotheses: (1) long-term warming can induce the different alterations in mineral nutrients and metabolites of A. faxoniana at high and low altitudes; and (2) the low-altitude boundary of A. faxoniana is more susceptible to warming than the high-altitude.
Materials and Methods
Experimental Site Description
The experiments were conducted at the Wanglang National Nature Reserve (32°49′–33°02′N, 103°55′–104°10′E), which is located in the eastern Tibetan Plateau, China. This region is in the monsoon climate zone, and there are distinct dry (November to April) and wet (May to October) seasons. The maximum and minimum temperatures in summer and winter are 26.2 and −17.8°C, and the mean annual temperature is from 5.07 to 1.57°C. The annual mean precipitation ranges from 801 to 825 mm, with most rainfalls happening from May to August (Ran et al., 2013). According to the Chinese soil classification (Chinese Soil Taxonomy Research Group, 1995), the soil type is mountain dark brown forest soil (Tang et al., 2019). A. faxoniana is a dominant tree species and the natural distribution ranges from 2,600 to 3,500 m on shady slopes in this region (Ran et al., 2013). The forest canopy density was 0.80–0.85 with a stand density of 810–900 ha–1 of pure fir trees.
Experimental Design and Sampling
The in situ experiments were designed in early September 2007. To provide an effective and simple method of warming simulation in the field for maintaining precipitation, light, and gas exchange closest to natural levels, OTCs, a passive warming method, were used for above-ground systems warming in this study. OTC warming is a low-maintenance and cost-effective method for field experiments across environmental gradients or in remote areas, especially in alpine area (Marion et al., 1997; Walker et al., 2006). Sites used for OTC are with similar canopy (density of 0.85) and slope angles (35–40°). OTCs are located in open areas where rain and snow could fall into. The OTCs used in this study are 2 m high, with a ground area of 2.25 m2, tapering to 1.5 m2 at the open top. The materials are made of solar-transmitting material (8-mm thick acrylic fiber). The walls of the OTCs have a 35–40° inclination relative to the ground. It ensures that the raining water can fall into the OTC from the open tops and interstitial flow from the inclination. Just like most published papers (Xu et al., 2013; Descombes et al., 2020), each site at any elevation included six OTCs and six control plots, respectively. Control plots were established near the OTC plot at each site. Ten-year-old A. faxoniana saplings (35 ± 5 cm height) were transplanted into 18 OTCs and 18 ambient plots (controls) at altitudes of 2,600, 3,000, and 3,500 m, respectively (Supplementary Figure 1A). In order to minimize the difference in field, saplings were collected from the forest gaps of the corresponding elevations to make them have similar genetic background. Five saplings were included in each OTC or control plot (1.5 m × 1.5 m). Four saplings were planted with a distance of 0.5 m and one sapling was planted in the middle (Supplementary Figure 1B). Finally, a total of 60 saplings (30 warming and 30 control individuals) were included at each altitude. The distance between each subplot is more than 10 m. Six TP-2200 humidity and temperature loggers (Fotel, Shanghai, China) were installed into the experimental plots to measure the temperature and moisture (Xu et al., 2013, 2015). Twice completely interannual temperature and moisture data (from January to December) were collected in 2015 and 2016, respectively. Because the mean annual temperature did not deviate too much from 2015 to 2016, we used the mean temperature of the 2 years as a proxy for the temperatures during the experiment period (Table 1 and Supplementary Figure 2).
After 10 years, on August 26, 2017, only 1-year-old needles from the middle-planted sampling of each plot were harvested at 10:00–13:00, randomly divided into five parts. In order to wash out the inorganic matters from the air and fallen organic matters in the needles, they were washed by deionized water. Then, each part of the needles was wrapped using silver paper and stored into liquid nitrogen immediately for further analyzing. For the soil collection, 15 soil cores (5 cm in diameter) of 5–10 cm depth were taken from each temperature regime at each altitude. The three soil cores from each plot were combined to form one composite sample. Because one control and two OTCs chambers were destroyed by fallen trees at the altitude of 2,600 m, four biological replicates were used for 2,600 m and five biological replicates for 3,000 and 3,500 m, respectively.
Determination of Height and Needle δ13C Values
Plant height and diameter were measured on August 25, 2017. The abundance of stable carbon isotopes in the needle samples was determined with an isotope mass spectrometer (Finnegan MAT Delta V Advantage, Bremen, Germany). Needle samples were oven-dried at 70°C for 48 h. Dried samples (0.20 g) were finely ground, and the relative abundances of 13C and 12C were tested according to Chen et al. (2014). The carbon isotope composition was expressed as δ13C value. The overall precision of the δ values was 0.1‰, as determined by the repetitive samples. The δ13C of needle samples was described relative to the standard Pee Dee Belemnite as follow method: δ13C = (Rsample/Rstandard − 1) × 1,000. Rsample means the 13C/12C ratio of samples, and Rstandard that of standard (Farquhar et al., 1989).
Measurement of Mineral Nutrients in Needles and Soils
Needles and soils were dried at 80°C, and 0.20 g of samples was burned off organic matter using a muffle furnace (450°C). Tissue nutrient concentrations were determined by ash digestion (500°C for 6 h) dissolved in 10 ml of nitric acid and peroxide. Dried soil samples were also digested by nitric acid. The digested plant and soil materials were used to determine mineral nutrient concentrations of calcium (Ca), sodium (Na), aluminum (Al), manganese (Mn), iron (Fe), zinc (Zn), titanium (Ti), chrome (Cr), K, Mg, and P by an Agilent 710 ICP-OES spectrometer (Agilent, Palo Alto, CA, United States) according to instructions (Se et al., 2019; Tang et al., 2019). The total C and N concentrations were analyzed with an elementar vario MACRO cube (Elementar, Germany) (Zuo et al., 2016).
Metabolite Extraction From Needles of Abies faxoniana
The extraction of needle metabolites was performed according to Zhang et al. (2019). Exactly 0.15 g of frozen 1-year-old needles of A. faxoniana was ground in a mortar under liquid nitrogen. The metabolites were extracted using 100% precooled methanol (1,500 μl) and vortexed for 30 s. Then, 60 μl ribitol (0.2 mg ml–1) was added as the internal quantitative standard. Subsequently, they were ultrasonically processed for 15 min, and 750 μl of chloroform (−20°C) and 1,500 μl of dH2O (4°C) were added. After centrifugation at 4,000 × g for 15 min, the supernatants were transferred into new tubes. They were dried in a Termovap sample concentrator (DC150-1A, Yooning Instrument, Hangzhou, China). Then, 60 μl of methoxyamine hydrochloride in pyridine (15 mg ml–1) was added, and the mixture was incubated at 38°C for 16 h. Finally, 60 μl of Bis(trimethylsilyl)trifluoroacetamide (BSTFA), including 1% chlorotrimethylsilane, was added to the mixture and incubated for 60 min at room temperature.
GC–MS and Metabolic Profile Analysis
Metabolites were determined using non-targeted metabolomics with an Agilent 7890A GC/5975C MS system (Agilent, California, Palo Alto, United States). One microliter of extract was injected into a 5% phenyl methyl silox of an HP-5 MS capillary column with a 20:1 split ratio. The columns were 0.25 μm df and 30 m × 0.25 mm ID with 10 m Integra-Guard (Agilent J and W scientific, Folsom, CA, United States). The MS source temperature was 250°C, the MS quad temperature was 150°C and the injection temperature was 280°C. After an initial isothermal heating at 80°C for 5 min, there was a 20°C min–1 ramp-up to 300°C followed by heating for 6 min. Helium (1.0 ml min–1) was used as the carrier gas. The total running time was 22 min. Using a full-scanning method, mass spectrometry was conducted at a range from 35 to 780 (m/z). The corresponding retention times, mass charge ratio and resolution of peaks were resolved by AMDIS software and compared with the NIST1 commercial database according to Behrends et al. (2011). In this study, only one salinization substance (TMS) was used for analysis. The metabolic pathway was constructed by combining the KEGG metabolic database2.
We used XCMS software3 for raw signal extraction, data filtering and peak identification. The setting parameters were as follows: XCMS Set (snthresh = 3, fwhm = 3, mzdiff = 0.5, steps = 2, step = 0.1, max = 300) and retcor [plottype = c (deviation), method = obiwarp, minfrac = 0.3, and bandwidth (bw) = 2]. The final data were imported into Excel tables for multiple linear regression analyses. A quantitative compound (QC) was used to remove misassignment and background noise and to ensure the accuracy of the identified true chemical entities. Before statistical analysis, the total mass of the signal integration area was used to normalize each sample. Then, the normalized data were imported into SIMCA-P 11.0 software (Umetrics AB, Umea, Sweden). Finally, two-way ANOVA was used to identify the effects of altitude, warming, and their interaction. Only metabolites with p < 0.05 and the fold change >1.2 were considered as differentially expressed metabolites (DEMs), according to Yokoi et al. (2015).
Statistical Data Analyses
To detect the difference among treatments, ANOVA with Tukey’s test at the p < 0.05 level was used for identifying significant differences with the SPSS 16.0 software (SPSS, Inc., Chicago, IL, United States). The metabolites were standardized using SIMCA-P 11.04. To better understand the differences between the ambient and OTC warming treatments, t-tests (p < 0.05) were performed within altitudes. In order to reduce differences between groups and make a better separation between groups, partial least squares discrimination analysis (PLS-DA) was performed to display the variation patterns in all the metabolites data instead of principal component analysis (PCA). Orthogonal partial least squares discrimination analysis (OPLS-DA) was conducted to demonstrate the alteration caused by OTC in each attitude. Matrix scatter plot was performed by the package GGally in R. The Pearson’s product–moment correlation coefficients between metabolites and the element contents were calculated and visualized in a heatmap. MetaboAnalyst 4.0 software5 was used to construct a heat map diagram of the metabolite categories (Chong et al., 2018).
Results
Soil Mineral Nutrient Changes in Response to Open-Top Chamber Warming at Three Altitudes
Under ambient conditions, the concentrations of Ca, Na, Al, Mg, and K decreased with a descending elevation gradient (Table 2). OTC warming significantly increased the concentrations of Na, Al, Mn, Zn, and the ratios of C:P and N:P, but decreased the concentrations of Mg, K, and P in soils at 2,600 m. However, it significantly decreased Ca, Na, Zn, Al, and Mg concentrations at 3,500 m. Soil C and N concentrations were significantly increased by OTC warming only at 3,000 m.
Changes in Physiological Traits of Abies faxoniana in Response to Open-Top Chamber Warming at Three Altitudes
The effect of altitude on δ13C values in needles was significant, showing lower values along a descending elevation gradient (Figure 1A). OTC warming did not significantly change needle δ13C values at any altitude. The growth of height and diameter increased with decreasing altitude (Figures 1B,C), and there was a significant difference between altitudes of 2,600 and 3,500 m (Pa < 0.05). However, OTC warming did not cause significant change in the height and diameter of A. faxoniana at any altitude.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766368/fpls-12-766368-HTML/image_m/fpls-12-766368-g001.jpg)
Figure 1. The needle δ13C values (A) and the growth of height (B) and diameter (C) of A. faxoniana at 3,500, 3,000, and 2,600 m after 10-year warming. Values denoted the different letters indicates a significance at p < 0.05 according to Tukey test. Pa, altitude effect; Pt, open-top chamber (OTC) effect; Pa × t, altitude and OTC interaction effect; H1CK, 3,500 m control; H1OTC, 3,500 m OTC warming; H2CK, 3,000 m control; H2OTC, 3,000 m OTC warming; H3CK, 2,600 m control; H3OTC, 2,600 m OTC warming.
Under ambient conditions, the needle Ca and K concentrations were significantly higher at 2,600 m than at 3,500 m (Figure 2). The concentrations of P, Al, and Zn showed increased tendencies with the decline of elevations. OTC warming significantly increased the concentrations of K, Al, and Zn in needles at 2,600 m. The overall effect on element concentration is not significant, but p < 0.05 (Pa) was observed in some elements like P, K, Al, Fe, and Zn. Therefore, we further explored elements correlation, and heatmap cluster analysis for needle ion contents. As shown in Figure 3, the correlations of P with K, Mn, Fe, and Al were changed by warming. Figure 3 and Supplementary Figure 3A showed a clearly changed tendency of the needle nutrients. Among these, P, Al, K, and Zn clustered into the same category and increased with the increased temperature.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766368/fpls-12-766368-HTML/image_m/fpls-12-766368-g002.jpg)
Figure 2. Mineral nutrient contents (μg g–1 Dw) in needles of A. faxoniana at 3,500, 3,000, and 2,600 m under controls and OTCs. Values denoted the different letters indicates a significance at p < 0.05 according to Tukey test. Pa, altitude effect; Pt, OTC effect; Pa × t, altitude and OTC interaction effect; H1, the altitude of 3,500 m; H2, the altitude of 3,000 m; H3, the altitude of 2,600 m; the green bars indicate the control conditions; the blue bars indicate the OTC warming conditions.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766368/fpls-12-766368-HTML/image_m/fpls-12-766368-g003.jpg)
Figure 3. Scatter matrix plot: relationship between needle elements in CK groups at the three attitudes, OTC groups at the three attitudes, and all groups include CK groups and OTC groups at the three attitudes. *0.01 < P ≤ 0.05; **0.001 < P ≤ 0.01; ***P < 0.001.
Metabolic Changes in Responses to Open-Top Chamber Warming at Three Altitudes
In this study, 71 metabolites were identified by GC–MS in the needles of A. faxoniana (Supplementary Table 1), mainly including amino acids, organic acid, and sugars. However, a large number of metabolites cannot be positively identified by using GC–MS, especially secondary metabolites. To visually present the altitude and warming effects, the control and warming groups of the same altitude were divided into the same clusters in the PLS-DA, indicating a clear altitude effect (Supplementary Figure 4). The results of the OPLS-DA indicated that the control and warming samples at each altitude could be clearly separated (Figure 4). Organic acids, sugars, and amino acids were greatly influenced by warming (Figure 4). A heatmap of the 28 samples also showed that the most important factor to causing the obvious separation was altitude but not warming (Supplementary Figure 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766368/fpls-12-766368-HTML/image_m/fpls-12-766368-g004.jpg)
Figure 4. Score plots and loading diagram of orthogonal partial least squares discrimination analysis (OPLS-DA) metabolic profiles in the needles of A. faxoniana at 3,500, 3,000, and 2,600 m under controls and OTCs. H1, the altitude of 3,500 m; H2, the altitude of 3,000 m; H3, the altitude of 2,600 m.
Different Metabolites Between Altitudes Under Open-Top Chamber Warming Conditions
There were 28, 15, and 17 metabolites significantly affected by altitude, warming, and their interaction (Supplementary Table 1), while 23, 12, and 15 metabolites could be considered as DEMs, respectively (Table 3). Under ambient conditions, there were less increased than decreased metabolites compared between low and high altitude. Of these, the contents of most of amino acids and sugars (except sucrose) were higher at high than low altitudes. OTC warming significantly increased the abundances of most of the metabolites at 2,600 m but the trend is the opposite at 3,500 m, e.g., organic acids, ribose, uracil, urea, and threitol.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766368/fpls-12-766368-HTML/image_m/fpls-12-766368-t003.jpg)
Table 3. The significant fold changes of differentially expressed metabolites (DEMs) in needles of A. faxoniana as affected by altitude and warming.
For sugars, the concentrations of fructose-6-phosphate (F6P), glucose-6-phosphate (G6P), and xylose were significantly decreased by warming at 3,500 m but were not decreased at 2,600 m (Figure 5) while sucrose and ribose significantly increased, and xylose decreased at three altitudes. The most significant changes were found in amino acids. Beta-alanine, serine, phenylalanine, leucine, iso-leucine, threonine, valine, and glutamic acid were increased at 2,600 m but most of them were decreased at 3,500 m. Additionally, organic acids also showed different altitude-related variation. Among these metabolites, the most significant increase was found in oxalic acids, with an increase of almost 39-fold at 2,600 m, while the most significant decrease was found in valine, with a decrease of one fourth at 3,500 m by warming (Figure 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766368/fpls-12-766368-HTML/image_m/fpls-12-766368-g005.jpg)
Figure 5. Metabolic pathways in needles of A. faxoniana at 3,500, 3,000, and 2,600 m caused by warming. Values were used for intensities of the detected metabolites. H1, the altitude of 3,500 m; H2, the altitude of 3,000 m; H3, the altitude of 2,600 m. PEP, phosphoenolpyruvic acid; 3-PGA, 3-phosphoglyceric acid.
Correlation Between the Ionome and Metabolome
Comprehensive analysis with the Pearson correlation coefficient were conducted to understand the relationships between soil mineral elements and needle metabolites (Figure 6). Under warming conditions, Na, K, Mg, Al, Zn, Fe, and P were positively correlated with phenylalanine and bata-alanine. The correlations among oxalic acid and P, Fe, Na, and K were negative. Fumaric acid and malic acid were positively correlated with Zn, Fe, K, Na, and Al at ambient conditions, while the correlations were negative after OTC warming. Moreover, the correlations of needle nutrients and metabolites at different altitudes were also explored. Under warming conditions, P was positively correlated with amino acids (threonine, serine, and valine), benzoic acid, G6P, and citric acid at 2,600 m (Figure 7). In addition, oxalic acid, sucrose, and phenylalanine were positively correlated with Al, and the correlation between quinic acid and Mn was increased at 2,600 m under OTC warming (Figure 7). However, P was negatively correlated to the tyrosine, leucine, proine, shikimic acid, oxalic acid, fructose, G6P, and fumaric acid while Al was negatively correlated with malic acid, leucine, glucose, and oxalic acid at 3,500 m. In addition, Mn was negatively correlated with many metabolites involved in carbohydrate metabolism such as malic acid, oxalic acid, fructose, citric acid, fumaric acid, glucose, glycerol, and shikimic acid at 3,500 m (Figure 7).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766368/fpls-12-766368-HTML/image_m/fpls-12-766368-g006.jpg)
Figure 6. The heat maps of the correlation between soil mineral elements and needle metabolites in CK groups (A) and OTC groups (B). Horizontal axis indicates the needle metabolites with significant changes under warming conditions; vertical axis indicates the needle mineral elements with significant changes warming conditions. Abbreviations were listed in Supplementary Table 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766368/fpls-12-766368-HTML/image_m/fpls-12-766368-g007.jpg)
Figure 7. The heat maps of the correlation between needle mineral elements and needle metabolites in CK groups and OTC groups at the three altitudes. (A), H1CK; (B), H1OTC; (C) H2CK; (D) H2OTC; (E) H3CK; and (F) H3OTC. Horizontal axis indicates the needle metabolites with significant changes under warming conditions; vertical axis indicates the needle mineral elements with significant changes warming conditions. CK, control groups; OTC, warming groups. H1, 3,500 m; H2, 3,000 m; H3, 2,600 m. Abbreviations were listed in Supplementary Table 1.
Discussion
The Effects of Altitude and Long-Term Warming on Soil Nutrients
Different altitudes at small geographic scales, such as subalpine ranges, exhibit steep climatic gradients with decreasing temperature (Seastedt and Oldfather, 2021). Forest soil properties differ greatly across altitudes. The soil C:N ratio ranged from 13.49 to 16.77 at three altitudes (Table 2), which was higher than the average value (12.30) in China (Tian et al., 2010). When soil C:N is less than 25:1, soil nitrogen is sufficient, and organic carbon mineralization occurs more easily (Tessier, 2001). Therefore, regardless of growth conditions, the soil C and N supply could be sufficient at these study sites, and microbes utilize organic matter effectively (Xu et al., 2015). Soil C:P reflects the effectiveness of phosphorus utilization, and N:P is often used for the diagnosis of nutrient limitation of soil N or P. Soil and plant N:P ratios are closely related to each other (Fan et al., 2015). In China, the average C:P and N:P of soil are 52.70 and 3.90, respectively (Tian et al., 2010). Additionally, the deficiencies of P, Ca, and Mg are common with warming in the subalpine regions of western China (Ji et al., 2019; Yuan et al., 2019). In this study, C:P and N:P ratios were much higher in the three altitudes than the average values across China, indicating that soil P availability is relatively low, which is not conducive to plant growth. Warming significantly decreased the soil P concentration and increased the needle P concentration at low altitudes (Table 2), suggesting that soil P would be a main limiting factor for A. faxoniana at low altitudes with the increasing global temperature.
The concentrations of soil cations (Ca, Na, Al, Mg, and K) were increased along the altitude. However, OTC warming caused the opposite changes in the soil Na, Al, and Zn between lower and higher altitudes (Table 2). These mineral nutrients come mainly from rock weathering (Uhlig and von Blanckenburg, 2019). Warming can cause more significant changes in the soil microbial community and enzyme activities at low than at high altitude (Zhao et al., 2014; Xu et al., 2015). The microbial activities enhanced by warming at low altitude could accelerate soil nutrient cycling. Compared with the low altitude, the relatively lower annual temperature at high altitude caused a relatively lower rate of nutrient mineralization. At low altitude (2,600 m), the OTC warming effect (1.04°C) was distinctly higher than that at high altitude and would accelerate rock weathering. Therefore, more mineral nutrients, including heavy metals, would be released into the soil under OTC warming at low than high altitude.
Long-Term Warming Effects in Abies faxoniana Saplings at Different Altitudes
Needle δ13C can be a good proxy for evaluating long-term water use efficiency of trees (Gimeno et al., 2021). There was no significant change in needle δ13C values between the control, and OTCs showed there was no influence on water use efficiency in A. faxoniana under long-term warming (Terwilliger et al., 2001), further indicating that water supply was sufficient in the control and warming groups (Figure 1A).
Sugars play important roles in regulating osmotic potential in plants under adverse conditions. Plants have developed an ability to resist cold damage by synthesizing numerous protective substances (e.g., soluble sugars and prolines) and proteins when they are exposed to low temperature (Ding et al., 2019). The metabolic data indicated that the higher content of sugars detected at the high altitude was a sign of greater protection against chilling stress in A. faxoniana cells. Sugars also play important roles in energy acquisition, as they are sources of carbon and nitrogen (Larher et al., 2003). Warming can also induce declines in the concentration of starch and transcripts of genes involved in glycolysis, TCA cycle, and secondary metabolism, accompanied by an increase in glucose (Vicente et al., 2019). In our study, F6P, G6P, and xylose were significantly decreased by OTC warming at high altitude but were not decreased at low altitude. F6P and G6P are important compounds in the first steps of glycolysis. G6P is the starting point for energy metabolism and is also used to generate NADPH in plants (Schluepmann et al., 2012). We think that on the one hand, the decreases in F6P and G6P contents can be attributed to the lower efficiency of hexose phosphorylation by hexokinase rather than the reduction in glucose supply. On the other hand, under warming conditions, A. faxoniana growing at high altitudes possibly consumes less energy than that growing at low altitude (Zhao et al., 2020; Zhang et al., 2021). Therefore, A. faxoniana at low altitude need to obtain more energy to cope with warming, whereas it is most likely to adopt an energy conservation strategy at high altitude.
Among the identified metabolites, the greatest changes caused by OTC warming among altitudes were amino acids. The content of γ-aminobutyric acid (GABA) was significantly increased due to OTC warming. GABA is an essential intermediate in the biosynthesis of many amino acids in plants (Bouché and Fromm, 2004; Batushansky et al., 2014). The increase of GABA, glutamic acid, and alanine in plants is an essential adaption to temperature elevation (Al-Quraan et al., 2012; Li et al., 2020). Many amino acids were significantly increased at 2,600 m (Figure 7). In particular, valine, beta-alanine, glutamic acid, and phenylalanine showed opposite trends due to OTC warming between the low and high altitudes. Moreover, long-term warming also significantly changed the levels of organic acids in the needles. At low altitude, glyceric acid, malonic acid, and oxalic acid were significantly increased in needles, whereas gluconic acid and heptanoic acid were significantly decreased (Figure 7). Notably, the level of oxalic acid was increased greatly at low but not at high altitude. Leaf respiration is commonly positive correlation with the increased temperature (Reich et al., 2016). Therefore, the increase of organic acids related to carbohydrate metabolism possibly indicates the enhanced leaf respiration under warming conditions.
The Interaction Between Nutrients and Physiological Responses
Because soils at high altitudes are exposed to harsher environmental conditions, soil mineral nutrition will be different with that in low-altitude areas. Rock weathering and soil microbial community and enzyme activities can be enhanced greatly under warming at low altitudes (Zhao et al., 2014; Xu et al., 2015; Figure 8). In addition, the availability of soil nutrient and elemental composition can greatly affect plant growth and metabolic processes (Rivas-Ubach et al., 2012). Elevation-driven changes in plant nutrients are associated with the content and quality of soil organic matter and the microbial activity (Mayor et al., 2017). For plant growth, Mn is involved in many physiological processes, but only in small quantity (Foy et al., 1978; Millaleo et al., 2010). In contrast, Al is not necessary but will be toxic to plants at high concentration (Delhaize and Ryan, 1995; Kochian et al., 2005). Tree growth was significantly and negatively related to the concentrations of Mn and Al in soil that would have a toxicity effect due to their solubilization (Zemunik et al., 2018). The soil pH can drastically influence the availability of metals (Lambers et al., 2008), and plant accumulation of Mn and Al was likely due to soil acidification at low altitude (Arduini et al., 1998; Yuan et al., 2019; Figure 8). Therefore, only a small increase in height of A. faxoniana by warming can be due to a decrease in soil P and an increase in Al and Mn at low altitude.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/766368/fpls-12-766368-HTML/image_m/fpls-12-766368-g008.jpg)
Figure 8. Diagram of A. faxoniana responding to warming at low and high altitudes. F6P, fructose-6-phosphate; G6P, glucose-6-phosphate.
According to previous studies, plant amino acid metabolism can be improved by temperature change. In addition to osmotic regulation, the increased amino acids also play important roles in protecting plants against pathogenic bacteria and heavy metals (Xie et al., 2013; Sonawala et al., 2018). However, we did not find visible disease symptoms on A. faxoniana at 2,600 m, and there was a large decrease in the secondary metabolites including catechin and dehydroascorbic acid. Therefore, the disease effects can be excluded. Warming triggered amino acid alteration can significantly up-regulate stress-related proteins (e.g., chalcone isomerase, glutathione S-transferase, and isoflavone reductase) to combat heavy metals (Bilal et al., 2019). Therefore, the increased amino acids by OTC warming at low altitude probably were used to inactivate heavy metals by forming amino acid-cation-polyphosphate complexes to reduce Al and Mn mobility (Figure 8). In addition, oxalic acid can balance the charges of positive ions and amino acids present in the plant cells (Liang et al., 2009; Ma et al., 2011). It has been reported that oxalic acid can form complexes with Al based on the concentration of oxalic acid (Nordstrom and May, 1996; Nakata, 2003). Both the increased internal tolerance and the exclusion of Al from plants are related to oxalate production (Zheng et al., 1998). The increased amount of oxalic acid in needles at low altitudes sequesters Al as non-toxic Al-oxalate salts (Figure 8). Under warming at 2,600 m, amino acid and carbohydrate metabolism related metabolites showed a positive correlation with P, K, Al, and Mn. Importantly, A. faxoniana at low altitude have to increase metabolite contents to cope with Al and Mn toxicity. However, concentrations of nutrients in needles of A. faxoniana at high altitude did not change along with soil nutrients. Many amino acids, organic acids, F6P, and G6P were decreased at high altitude with OTC warming (Figure 8). Succinic acid, glyceric acid, F6P, and G6P are closely related to glycolysis and TCA cycle. Many amino acids and metabolites involved in carbohydrate metabolism showed negative correlation to P, Al, and Mn concentrations in needles at high altitude. This indicates trees could have less energy production and accumulate more sucrose at high altitude (Figure 8).
Overall, OTC warming can directly affect plant physiological processes and indirectly change soil nutrient conditions, especially of P and heavy metals, and sequentially impact plant growth. Contrary to soil nutrient alterations, cations in A. faxoniana needles increased with decreasing altitude. The inconsistent nutrient changes indicated that Ca, Na, Mg, and K in soils were not the limiting nutrients for the growth of A. faxoniana. The warming-facilitated exudation of organic acids from plants to soils can reduce the soil pH value further. In addition, trees experience physiological reprogram with increasing age, including decreased growth rates and shifted carbon distribution pattern (Rossi et al., 2008). Younger trees tend to be more sensitive to climate warming, but older trees show remarkable resistance to warming (Colangelo et al., 2021). In this study, 10-year-old young A. faxoniana seedlings were used. Therefore, the results might be different from the adult trees at certain extent. In addition, there are some limitations in the current warming method. Compared with infrared heaters, it is unavoidably affected the changes of light, air exchange, and soil microbes by the chambers when using OTC warming. In a word, in the future, although the increased temperature will benefit the high-altitude boundary tree population, it is a disadvantage for the low-altitude boundary population. Therefore, the distribution area of A. faxoniana population is unlikely to be enlarged with global warming.
Conclusion
There are significant altitude effects on the soil properties and tree growth. Warming changes the soil cation concentrations differentially at various altitudes, and the concentrations of metals (e.g., Al and Mn) were increased greatly at low altitudes. Needles of A. faxoniana will compete for the mineral nutrients (e.g., Fe, K, and P) at low altitudes in case of warming. The metabolites that differed mostly among the altitudes were the sugars and amino acids. Oxalic acid could be considered an Al and Mn detoxification agent at low altitude. Combined with the growth conditions, we think there are contrasting physiological and metabolic responses of A. faxoniana to long-term warming between low and high altitudes, whereby A. faxoniana growing at low low-altitude boundary will be more susceptible to climate change in alpine regions.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author Contributions
HS performed the experiments and wrote and revised the manuscript. QH did some lab work. SZ designed the research, analyzed the data, and revised the manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by the Second Tibetan Plateau Scientific Expedition and Research Program (2019QZKK0404), the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA20020401), and the Fundamental Research Funds for the Central Universities.
Conflict of Interest
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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.766368/full#supplementary-material
Supplementary Table 1 | The identified metabolites in needles of A. faxoniana at 3500, 3000 and 2600 m caused by warming. H1, the altitude of 3500 m; H2 the altitude of 3000 m; H3, the altitude of 2600 m.
Footnotes
- ^ http://www.nist.gov/index.html
- ^ http://www.genome.jp/kegg/
- ^ www.bioconductor.org/
- ^ http://www.umetrics.com/simca
- ^ http://www.metaboanalyst.ca/faces/ModuleView.xhtml
References
Al-Quraan, N., Locy, B., and Singh, N. (2012). Heat and cold stresses phenotypes of Arabidopsis thaliana calmodulin mutants: regulation of gamma-aminobutyric acid shunt pathway under temperature stress. Int. J. Plant Biol. 3:e2. doi: 10.4081/pb.2012.e2
Anderson, J. T., Inouye, D. W., McKinney, A. M., Colautti, R. I., and Mitchell-Olds, T. (2012). Phenotypic plasticity and adaptive evolution contribute to advancing flowering phenology in response to climate change. Proc. Biol. Sci. 279, 3843–3852. doi: 10.1098/rspb.2012.1051
Arduini, I., Kettner, C., Godbold, D. L., Onnis, A., and Stefani, A. (1998). pH influence on root growth and nutrient uptake of Pinus pinaster seedlings. Chemosphere 36, 733–738. doi: 10.1016/S0045-6535(97)10116-3
Aspinwall, M. J., Drake, J. E., Campany, C., Vårhammar, A., Ghannoum, O., Tissue, D. T., et al. (2016). Convergent acclimation of leaf photosynthesis and respiration to prevailing ambient temperatures under current and warmer climates in Eucalyptus tereticornis. New Phytol. 212, 354–367. doi: 10.1111/nph.14035
Batushansky, A., Kirma, M., Grillich, N., Toubiana, D., Pham, P. A., Balbo, I., et al. (2014). Combined transcriptomics and metabolomics of Arabidopsis thaliana seedlings exposed to exogenous GABA suggest its role in plants is predominantly metabolic. Mol. Plant 7, 1065–1068. doi: 10.1093/mp/ssu017
Behrends, V., Tredwell, G. D., and Bundy, J. G. (2011). A software complement to AMDIS for processing GC-MS metabolomic data. Anal. Biochem. 415, 206–208. doi: 10.1016/j.ab.2011.04.009
Bilal, S., Shahzad, R., Khan, A. L., Al-Harrasi, A., Kim, C. K., and Lee, I. J. (2019). Phytohormones enabled endophytic Penicillium funiculosum LHL06 protects Glycine max L. from synergistic toxicity of heavy metals by hormonal and stress-responsive proteins modulation. J. Hazard Mater. 379:11. doi: 10.1016/j.jhazmat.2019.120824
Bouché, N., and Fromm, H. (2004). GABA in plants: just a metabolite? Trends Plant Sci. 9, 110–115. doi: 10.1016/j.tplants.2004.01.006
Carter, K. R., Wood, T. E., Reed, S. C., Butts, K. M., and Cavaleri, M. A. (2021). Experimental warming across a tropical forest canopy height gradient reveals minimal photosynthetic and respiratory acclimation. Plant Cell Environ. 44, 2879–2897. doi: 10.1111/pce.14134
Chen, J., Duan, B., Wang, M., Korpelainen, H., and Li, C. (2014). Intra- and inter-sexual competition of Populus cathayana under different watering regimes. Funct. Ecol. 28, 124–136. doi: 10.1111/1365-2435.12180
Chevin, L. M., Collins, S., and Lefèvre, F. (2013). Phenotypic plasticity and evolutionary demographic responses to climate change: taking theory out to the field (ed Hoffmann A). Funct. Ecol. 27, 967–979. doi: 10.1111/j.1365-2435.2012.02043.x
Chinese Soil Taxonomy Research Group (1995). Chinese Soil Taxonomy Research Group of Institute of Soil Science, Cooperative Research Group on Chinese Soil Taxonomy (CRGCST) Chinese Soil Taxonomy (Revised Proposal). Beijing: China Agriculture Press.
Chong, J., Soufan, O., Li, C., Caraus, I., Li, S., Bourque, G., et al. (2018). MetaboAnalyst 4.0: towards more transparent and integrative metabolomics analysis. Nucleic Acids Res. 46, 486–494. doi: 10.1093/nar/gky310
Colangelo, M., Camarero, J. J., Gazol, A., Piovesan, G., Borghetti, M., Baliva, M., et al. (2021). Mediterranean old-growth forests exhibit resistance to climate warming. Sci. Total Environ. 801:149684. doi: 10.1016/j.scitotenv.2021.149684
Conlisk, E., Castanha, C., Germino, M. J., Veblen, T. T., Smith, J. M., Kueppers, L. M., et al. (2017). Declines in low-elevation subalpine tree populations outpace growth in high-elevation populations with warming. J. Ecol. 105, 1347–1357. doi: 10.1111/1365-2745.12750
Dang, H., Zhang, Y., Zhang, K., Jiang, M., and Zhang, Q. (2012). Climate-growth relationships of subalpine fir (Abies fargesii) across the altitudinal range in the Shennongjia Mountains, central China. Clim. Change 117, 903–917. doi: 10.1007/s10584-012-0611-5
Delhaize, E., and Ryan, P. R. (1995). Aluminum toxicity and tolerance in plants. Plant Physiol. 107, 315–321. doi: 10.1104/pp.107.2.315
Descombes, P., Kergunteuil, A., Glauser, G., Rasmann, S., and Pellissier, L. (2020). Plant physical and chemical traits associated with herbivory in situ and under a warming treatment. J. Ecol. 108, 733–749. doi: 10.1111/1365-2745.13286
Dhuli, P., Rohloff, J., and Strimbeck, G. R. (2014). Metabolite changes in conifer buds and needles during forced bud break in Norway spruce (Picea abies) and European silver fir (Abies alba). Front. Plant Sci. 5:706. doi: 10.3389/fpls.2014.00706
Ding, Y., Shi, Y., and Yang, S. (2019). Advances and challenges in uncovering cold tolerance regulatory mechanisms in plants. New Phytol. 222, 1690–1704. doi: 10.1111/nph.15696
D’Orangeville, L., Houle, D., Cote, B., and Duchesne, L. (2014). Soil response to a 3-year increase in temperature and nitrogen deposition measured in a mature boreal forest using ion-exchange membranes. Environ. Monit. Assess 186, 8191–8202. doi: 10.1007/s10661-014-3997-x
Duan, A., and Xiao, Z. (2015). Does the climate warming hiatus exist over the Tibetan Plateau. Sci. Rep. 5:13711. doi: 10.1038/srep13711
Duputie, A., Rutschmann, A., Ronce, O., and Chuine, I. (2015). Phenological plasticity will not help all species adapt to climate change. Glob. Change Biol. 21, 3062–3073. doi: 10.1111/gcb.12914
Fan, H., Wu, J., Liu, W., and Yuan, Y. (2015). Linkages of plant and soil C:N:P stoichiometry and their relationships to forest growth in subtropical plantations. Plant Soil. 392, 127–138. doi: 10.1007/s11104-015-2444-2
Farquhar, G. D., Hubick, K. T., Condon, A. G., and Richards, R. A. (1989). “Carbon isotope fractionation and plant water-use efficiency,” in Stable Isotopes in Ecological Research, eds P. W. Rundel, J. R. Ehleringer, and K. A. Nagy (New York, NY: Springer), 21–40. doi: 10.1007/978-1-4612-3498-2_2
Foy, C. D., Chaney, R. L., and White, M. C. (1978). The physiology of metal toxicity in plants. Ann. Rev. Plant Physiol. 29, 511–566. doi: 10.1146/annurev.pp.29.060178.002455
Franks, S. J., Weber, J. J., and Aitken, S. N. (2014). Evolutionary and plastic responses to climate change in terrestrial plant populations. Evol. Appl. 7, 123–139. doi: 10.1111/eva.12112
Gimeno, T. E., Campany, C. E., Drake, J. E., Barton, C. V. M., Tjoelker, M. G., Ubierna, N., et al. (2021). Whole-tree mesophyll conductance reconciles isotopic and gas-exchange estimates of water-use efficiency. New Phytol. 229, 2535–2547. doi: 10.1111/nph.17088
Guo, M. M., Zhang, Y. D., Wang, X. C., Gu, F. X., and Liu, S. R. (2018). The responses of dominant tree species to climate warming at the treeline on the eastern edge of the Tibetan Plateau. For. Ecol. Manage. 425, 21–26. doi: 10.1016/j.foreco.2018.05.021
Hancock, R. D., Morris, W. L., Ducreux, L. J., Morris, J. A., Usman, M., Verrall, S. R., et al. (2014). Physiological, biochemical and molecular responses of the potato (Solanum tuberosum L.) plant to moderately elevated temperature. Plant Cell Environ. 37, 439–450. doi: 10.1111/pce.12168
Harvey, J. E., Smiljanić, M., Scharnweber, T., Buras, A., Cedro, A., Cruz-García, R., et al. (2020). Tree growth influenced by warming winter climate and summer moisture availability in northern temperate forests. Glob. Chang Biol. 26, 2505–2518. doi: 10.1111/gcb.14966
Huang, R., Zhu, H. F., Liu, X. H., Liang, E. Y., Griessinger, J., Wu, G. J., et al. (2017). Does increasing intrinsic water use efficiency (iWUE) stimulate tree growth at natural alpine timberline on the southeastern Tibetan Plateau? Glob. Planet Change 148, 217–226. doi: 10.1016/j.gloplacha.2016.11.017
Ji, H. W., Ossipov, V., Du, B. M., Wen, J. H., and Liu, C. J. (2019). Differences in the relationship between metabolomic and ionomic traits of Quercus variabilis growing at contrasting geologic-phosphorus sites in subtropics. Plant Soil 439, 339–355. doi: 10.1007/s11104-019-04020-1
Kochian, L. V., Pineros, M. A., and Hoekenga, O. A. (2005). The physiology, genetics and molecular biology of plant aluminum resistance and toxicity. Plant Soil 274, 175–195. doi: 10.1093/jxb/erq272
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
Kueppers, L. M., Conlisk, E., Castanha, C., Moyes, A. B., Germino, M. J., de Valpine, P., et al. (2017). Warming and provenance limit tree recruitment across and beyond the elevation range of subalpine forest. Glob. Change Biol. 23, 2383–2395. doi: 10.1111/gcb.13561
Lambers, H., Raven, J. A., Shaver, G. R., and Smith, S. E. (2008). Plant nutrient-acquisition strategies change with soil age. Trends Ecol. Evol. 23, 95–103. doi: 10.1016/j.tree.2007.10.008
Larher, F. R., Aziz, A., Gibon, Y., Trotel-Aziz, P., Sulpice, R., and Bouchereau, A. (2003). An assessment of the physiological properties of the so-called compatible solutes using in vitro experiments with leaf discs. Plant Physiol. Biochem. 41, 657–666. doi: 10.1016/S0981-9428(03)00076-7
Li, Z., Cheng, B., Peng, Y., and Zhang, Y. (2020). Adaptability to abiotic stress regulated by γ-aminobutyric acid in relation to alterations of endogenous polyamines and organic metabolites in creeping bentgrass. Plant Physiol. Biochem. 157, 185–194. doi: 10.1016/j.plaphy.2020.10.025
Liang, E., Wang, Y., Eckstein, D., and Luo, T. (2011). Little change in the fir tree-line position on the southeastern Tibetan Plateau after 200 years of warming. New Phytol. 190, 760–769. doi: 10.1111/j.1469-8137.2010.03623.x
Liang, Y., Strelkov, S. E., and Kav, N. N. V. (2009). Oxalic acid-mediated stress responses in Brassica napus L. Proteomics 9, 3156–3173.
Luo, C., Rodriguezr, L. M., Johnston, E. R., Wu, L., Cheng, L., Xue, K., et al. (2013). Soil microbial community responses to a decade of warming as revealed by comparative metagenomics. Appl. Environ. Microbiol. 80:1777. doi: 10.1128/AEM.03712-13
Ma, Y., Guo, L. Q., Wang, H. X., Bai, B., and Shi, D. C. (2011). Accumulation, distribution, and physiological contribution of oxalic acid and other solutes in an alkali-resistant forage plant, Kochia sieversiana, during adaptation to saline and alkaline conditions. J. Plant Nutr. Soil Sci. 174, 655–663. doi: 10.1002/jpln.201000337
Marion, G. M., Henry, G. H. R., Freckman, D. W., Johnstone, J., Jones, G., Jones, M. H., et al. (1997). Open-top designs for manipulating field temperature in high-latitude ecosystems. Glob. Change Biol. 3, 20–32. doi: 10.1111/j.1365-2486.1997.gcb136.x
Mason, C. J., Keefover-Ring, K., Villari, C., Klutsch, J. G., Cook, S., Bonello, P., et al. (2019). Anatomical defences against bark beetles relate to degree of historical exposure between species and are allocated independently of chemical defences within trees. Plant Cell Environ. 42, 633–646. doi: 10.1111/pce.13449
Mayor, J. R., Sanders, N. J., Classen, A. T., Bardgett, R. D., Clément, J. C., Fajardo, A., et al. (2017). Elevation alters ecosystem properties across temperate treelines globally. Nature 542, 91–95. doi: 10.1038/nature21027
Millaleo, R., Reyes-Diaz, M., Ivanov, A. G., Mora, M. L., and Alberdi, M. (2010). Manganese as essential and toxic element for plants: transport, accumulation and resistance mechanisms. J. Soil Sci. Plant Nutr. 10, 476–494. doi: 10.1016/j.ecoenv.2020.111687
Nakata, P. A. (2003). Advances in our understanding of calcium oxalate crystal formation and function in plants. Plant Sci. 164, 901–909. doi: 10.1016/S0168-9452(03)00120-1
Nordstrom, D. K., and May, H. M. (1996). “Aqueous equilibrium data for mononuclear aluminum species,” in Environment Chemistry of Aluminum, ed. G. Sposito (Boca Raton, FL: CRC Press), 39–80. doi: 10.1201/9780138736781-2
Panthi, S., Fan, Z. X., van der Sleen, P., and Zuidema, P. A. (2020). Long-term physiological and growth responses of Himalayan fir to environmental change are mediated by mean climate. Glob. Change Biol. 26, 1778–1794. doi: 10.1111/gcb.14910
Pepin, N., Bradley, R. S., Diaz, H. F., Baraer, M., Caceres, E. B., Forsythe, N., et al. (2015). Elevation-dependent warming in mountain regions of the world. Nat. Clim. Change. 5, 424–430. doi: 10.1038/nclimate2563
Ponocna, T., Spyt, B., Kaczka, R., Buntgen, U., and Treml, V. (2016). Growth trends and climate responses of Norway spruce along elevational gradients in East-Central Europe. Trees 30, 1633–1646. doi: 10.1007/s00468-016-1396-3
Ran, F., Zhang, X. L., Zhang, Y. B., Korpelainen, H., and Li, C. Y. (2013). Altitudinal variation in growth, photosynthetic capacity and water use efficiency of Abies faxoniana Rehd. et Wils. Seedlings as revealed by reciprocal transplantations. Trees 27, 1405–1416. doi: 10.1007/s00468-013-0888-7
Reich, P. B., Sendall, K. M., Stefanski, A., Wei, X., Rich, R. L., and Montgomery, R. A. (2016). Boreal and temperate trees show strong acclimation of respiration to warming. Nature 531, 633–636. doi: 10.1038/nature17142
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
Rossi, S., Deslauriers, A., Anfodillo, T., and Carrer, M. (2008). Age-dependent xylogenesis in timberline conifers. New Phytol. 177, 199–208. doi: 10.1111/j.1469-8137.2007.02235.x
Rustad, L., Campbell, J., Marion, G., Norby, R., Mitchell, M., Hartley, A., et al. (2001). A meta-analysis of the response of soil respiration, net nitrogen mineralization, and aboveground plant growth to experimental ecosystem warming. Oecologia 126, 543–562. doi: 10.1007/s004420000544
Savolainen, O., Bokma, F., Garcia-Gil, R., Komulainen, P., and Repo, T. (2004). Genetic variation in cessation of growth and frost hardiness and consequences for adaptation of Pinus sylvestris to climatic changes. For. Ecol. Manage. 197, 79–89. doi: 10.1016/j.foreco.2004.05.006
Schluepmann, H., Berke, L., and Sanchez-Perez, G. F. (2012). Metabolism control over growth: a case for trehalose-6-phosphate in plants. J. Exp. Bot. 63, 3379–3390. doi: 10.1093/jxb/err311
Schrieber, K., Caceres, Y., Engelmann, A., Marcora, P., Renison, D., Hensen, I., et al. (2020). Elevational differentiation in metabolic cold stress responses of an endemic mountain tree. Environ. Exp. Bot. 171:103918. doi: 10.1016/j.envexpbot.2019.103918
Se, K. W., Ghoshal, S. K., and Wahab, R. A. (2019). Laser-induced breakdown spectroscopy unified partial least squares regression: an easy and speedy strategy for predicting Ca, Mg and Na content in honey. Measurement 136, 1–10. doi: 10.1016/j.measurement.2018.12.052
Seastedt, T. R., and Oldfather, M. (2021). Climate Change, Ecosystem processes and biological diversity responses in high elevation communities. Climate 9:87. doi: 10.3390/cli9050087
Sonawala, U., Dinkeloo, K., Danna, C. H., McDowell, J. M., and Pilot, G. (2018). Review: functional linkages between amino acid transporters and plant responses to pathogens. Plant Sci. 277, 79–88. doi: 10.1016/j.plantsci.2018.09.009
Stromme, C. B., Julkunen-Tiitto, R., Olsen, J. E., and Nybakken, L. (2018). The dioecious Populus tremula displays interactive effects of temperature and ultraviolet-B along a natural gradient. Environ. Exp. Bot. 146, 13–26. doi: 10.1016/j.envexpbot.2017.09.013
Takahashi, K., Azuma, H., and Yasue, K. (2003). Effects of climate on the radial growth of tree species in the upper and lower distribution limits of an altitudinal ecotone on Mount Norikura, central Japan. Ecol. Res. 18, 549–558. doi: 10.1046/j.1440-1703.2003.00577.x
Tang, D. T., Peng, G. Q., and Zhang, S. (2019). Age-related variations of needles and twigs in nutrient, nonstructural carbon and isotope composition along altitudinal gradients. J. Mt. Sci. 16, 1546–1558. doi: 10.1007/s11629-018-5344-0
Terwilliger, V. J., Kitajima, K., Le Roux-Swarthout, D. J., Mulkey, S., and Wright, S. J. (2001). Intrinsic water-use efficiency and heterotrophic investment in tropical leaf growth of two Neotropical pioneer tree species as estimated from delta C-13 values. New Phytol. 152, 267–281. doi: 10.1046/j.0028-646X.2001.00252.x
Tessier, J. T. (2001). Vernal photosynthesis nutrient retranslocation in Dryopteris intermedia. Am. Fern. J. 91, 187–196. doi: 10.1640/0002-8444(2001)091[0187:VPANRI]2.0.CO;2
Tian, H. Q., Chen, G. S., Zhang, C., Melillo, J. M., and Hall, C. A. (2010). Pattern and variation of C:N:P ratios in China’s soils: a synthesis of observational data. Biogeochemistry 98, 139–151. doi: 10.1007/s10533-009-9382-0
Uhlig, D., and von Blanckenburg, F. (2019). How slow rock weathering balances nutrient loss during fast forest floor turnover in Montane, temperate forest ecosystems. Front. Earth Sci. 7:159. doi: 10.3389/feart.2019.00159
Vicente, R., Bolger, A. M., Martínez-Carrasco, R., Pérez, P., Gutiérrez, E., Usadel, B., et al. (2019). De Novo transcriptome analysis of durum wheat flag leaves provides new insights into the regulatory response to elevated CO2 and high temperature. Front. Plant Sci. 10:1605. doi: 10.3389/fpls.2019.01605
Volder, A., Briske, D. D., and Tjoelker, M. G. (2013). Climate warming and precipitation redistribution modify tree-grass interactions and tree species establishment in a warm-temperate savanna. Glob. Chang Biol. 19, 843–857. doi: 10.1111/gcb.12068
Walker, M. D., Wahren, C. H., and Hollister, R. D. (2006). Plant community responses to experimental warming across the tundra biome. Proc. Natl. Acad. Sci. U.S.A. 103, 1342–1346. doi: 10.1073/pnas.0503198103
Wang, W. Z., Jia, M., Wang, G. X., Zhu, W. Z., and McDowell, N. G. (2017). Rapid warming forces contrasting growth trends of subalpine fir (Abies fabri) at higher- and lower-elevations in the eastern Tibetan Plateau. For. Ecol. Manage. 402, 135–144. doi: 10.1016/j.foreco.2017.07.043
Wei, Y., and Fang, Y. (2013). Spatio-temporal characteristics of global warming in the Tibetan Plateau during the last 50 years based on a generalised temperature zone-elevation model. PLoS One 8:e60044. doi: 10.1371/journal.pone.0060044
Wu, C. C., Peng, G. Q., Zhang, Y. B., Xu, X., Korpelainen, H., Berninger, F., et al. (2011). Physiological responses of Abies faxoniana seedlings to different non-growing-season temperatures as revealed by reciprocal transplantations at two contrasting altitudes. Can. J. For. Res. 41, 599–607. doi: 10.1139/X10-225
Wu, T., Liu, S., Lie, Z., Zheng, M., Duan, H., Chu, G., et al. (2020). Divergent effects of a 6-year warming experiment on the nutrient productivities of subtropical tree species. Forest Ecol. Manag. 461:117952. doi: 10.1016/j.foreco.2020.117952
Xie, X. Y., Weiss, D. J., Weng, B. S., Liu, J. C., Lu, H. L., and Yan, C. L. (2013). The short-term effect of cadmium on low molecular weight organic acid and amino acid exudation from mangrove (Kandelia obovata (S., L.) Yong) roots. Environ. Sci. Pollut. Res. 20, 997–1008. doi: 10.1007/s11356-012-1031-9
Xu, G., Chen, J., Berninger, F., Pumpanen, J., Bai, J. W., Yu, L., et al. (2015). Labile, recalcitrant, microbial carbon and nitrogen and the microbial community composition at two Abies faxoniana forest elevations under elevated temperatures. Soil Biol. Biochem. 91, 1–13. doi: 10.1016/j.soilbio.2015.08.016
Xu, G., Jiang, H., Zhang, Y. B., Korpelainen, H., and Li, C. Y. (2013). Effect of warming on extracted soil carbon pools of Abies faxoniana forest at two elevations. For. Ecol. Manage. 310, 357–365. doi: 10.1016/j.foreco.2013.08.038
Yin, H. J., Li, Y. F., Xiao, J., Xu, Z. F., Cheng, X. Y., and Liu, Q. (2013). Enhanced root exudation stimulates soil nitrogen transformations in a subalpine coniferous forest under experimental warming. Glob. Change Biol. 19, 2158–2167. doi: 10.1111/gcb.12161
Yin, H. J., Liu, Q., and Lai, T. (2008). Warming effects on growth and physiology in the seedlings of the two conifers Picea asperata and Abies faxoniana under two contrasting light conditions. Ecol. Res. 23, 459–469. doi: 10.1007/s11284-007-0404-x
Yokoi, N., Beppu, M., Yoshida, E., Hoshikawa, R., Hidaka, S., Matsubara, T., et al. (2015). Identification of putative biomarkers for prediabetes by metabolome analysis of rat models of type 2 diabetes. Metabolomics 11, 1277–1286. doi: 10.1007/s11306-015-0784-9
Yuan, J., Sun, N. X., Du, H. M., Muhammad, U., Kang, H. Z., Du, B. M., et al. (2019). Correlated metabolic and elemental variations between the leaves and seeds of oak trees at contrasting geologically derived phosphorus sites. Sci. Total. Environ. 691, 178–186. doi: 10.1016/j.scitotenv.2019.07.133
Zemunik, G., Davies, S. J., and Turner, B. L. (2018). Soil drivers of local-scale tree growth in a lowland tropical forest. Ecology 99, 2844–2852. doi: 10.1002/ecy.2532
Zhang, L., Sun, S., Liang, Y., Li, B., Ma, S., Wang, Z., et al. (2021). Nitrogen levels regulate sugar metabolism and transport in the shoot tips of crabapple plants. Front. Plant Sci. 12:626149. doi: 10.3389/fpls.2021.626149
Zhang, S., Tang, D. T., Korpelainen, H., and Li, C. Y. (2019). Metabolic and physiological analyses reveal that Populus cathayana males adopt an energy saving strategy to cope with phosphorus deficiency. Tree Physiol. 39, 1630–1645. doi: 10.1093/treephys/tpz074
Zhao, C. Z., Zhu, L. Y., Liang, J., Yin, H. J., Yin, C. Y., Li, D. D., et al. (2014). Effects of experimental warming and nitrogen fertilization on soil microbial communities and processes of two subalpine coniferous species in Eastern Tibetan Plateau, China. Plant Soil 382, 189–201. doi: 10.1007/s11104-014-2153-2
Zhao, H., Sun, S., Zhang, L., Yang, J., Wang, Z., Ma, F., et al. (2020). Carbohydrate metabolism and transport in apple roots under nitrogen deficiency. Plant Physiol. Biochem. 155, 455–463. doi: 10.1016/j.plaphy.2020.07.037
Zheng, S. J., Ma, J. F., and Matsumoto, H. (1998). High aluminum resistance in buckwheat - I. Al-induced specific secretion of oxalic acid from root tips. Plant Physiol. 117, 745–751. doi: 10.1104/pp.117.3.745
Ziello, C., Estrella, N., Kostova, M., Koch, E., and Menzel, A. (2009). Influence of altitude on phenology of selected plant species in the Alpine region (1971–2000). Clim. Res. 39, 227–234. doi: 10.3354/cr00822
Keywords: warming, OTC, altitude, metabolism, adaption, sub-alpine
Citation: Song H, Han Q and Zhang S (2021) Low-Altitude Boundary of Abies faxoniana Is More Susceptible to Long-Term Open-Top Chamber Warming in the Eastern Tibetan Plateau. Front. Plant Sci. 12:766368. doi: 10.3389/fpls.2021.766368
Received: 02 September 2021; Accepted: 03 November 2021;
Published: 03 December 2021.
Edited by:
Huajun Yin, Chengdu Institute of Biology, Chinese Academy of Sciences (CAS), ChinaReviewed by:
Xiaomin Zhu, Chengdu Institute of Biology, Chinese Academy of Sciences (CAS), ChinaZiliang Zhang, University of Illinois at Urbana-Champaign, United States
Copyright © 2021 Song, Han and Zhang. 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: Sheng Zhang, shengzhang@scu.edu.cn