Experimental site. The experimental site is located in the Arboretum of Southwest Forestry University (Kunming, Yunnan province, E102°46′, N25°03′). It is in the subtropical plateau monsoon climate zone, with an altitude of 1964 m, short frost period, mild climate, average annual temperature of 16.5°C, average annual precipitation of 1035 mm, average annual relative humidity of 67%, and red loam soil type. The temperature inside the shelter is 18.5–37.0℃, and the relative humidity is 22.3–48.0%.
Seedling preparation. Seeds for seedlings were obtained from P. yunnanensis seed orchard in Midu County, Dali City, Yunnan Province. The seeds were cultivated in seedling bags and grown in nursery for 2 years. 200 two-year-old P. yunnanensis seedlings from Yiliang Garden Forestry were selected on August 15, 2020, and transplanted into plastic flowerpots with a diameter of 26.5 cm, a base diameter of 21 cm, and a height of 21 cm. Each pot was filled with an equal amount of 6 kg of sieved soil mixture (red loam: humus = 3:2), with a tray was placed at the bottom, 1 seedling per pot. The transplanted seedlings were placed on the experimental plot, which was covered with mulch to prevent underground moisture from affecting the potted plants, and a rain shelter was built to ensure good ventilation and no shade. The soil water content of the potted plants was under completely controlled conditions. The soil moisture content of the seedlings was maintained at field holding capacity after planting to ensure the healthy growth of the seedlings. Field water holding capacity was determined by ring knife method, The field water holding capacity of the soil used in the experiment was 25.94%, and soil bulk density was 1.14g·cm-3. After the seedlings were transplanted, the root collar diameter and height of all the seedlings were measured with vernier calipers and tape measures with an accuracy of 0.01 mm and 0.1 cm, respectively. Before the experiment started, seedling height, ground diameter and biomass were 20.84 ± 1.70 cm, 18.60 ± 0.47 mm.
Application of drought treatments. All the seedlings with consistent growth were selected, numbered and divided into four groups The difference of seedling height and root collar diameter between the four groups were tested to ensure that the grouping is uniform, 40 plants per group. When there is no significant difference in seedling height and ground diameter between the four groups, hang a tag and record the corresponding number. Soil moisture and density were measured by the ring-knife method before the experiment, and then soil mass moisture content and volumetric moisture content were calculated, and real-time soil volumetric moisture content was detected with a soil moisture meter during the experiment. The potted weighing control method was used to simulate natural drought conditions, and four moisture gradients were set, namely, suitable moisture (CK), light drought stress (LD), moderate drought stress (MD) and severe drought stress (SD), with soil moisture gradients were set at 80% ± 5%, 65% ± 5%, 50% ± 5% and 35% ± 5% of the field water holding capacity, respectively; The actual moisture content was maintained at 19.45–22.05%, 15.56–18.16%, 11.67–14.27%, and 7.78–10.38%, respectively23.The height and ground diameter of all viable seedlings were measured and recorded, and the water supply was then stopped to allow the soil water content to fall naturally to a predetermined range, namely, weighed and sampled as the initial value. The actual soil water content was measured using a Procheck handheld multifunctional soil moisture meter (Decagon, USA), controlled by the potted weighing method23,27, and all potted plants were weighed daily at 17:00, the current weight was noted, and water control or watering was performed according to the target weight. The drought stress experiment started on March 14, 2021, and ended on May 12, 2021, lasting for 60 days in total. 40 seedlings in each of the 4 treatments, after calculating the weight of each pot.
Measurements and sampling. Measurements were performed from March 14, 2021, including growth measurements and destructive sampling on the 15th, 30th, 45th, and 60th days after the the start of the drought stress experiment (i.e., March 29, April 13, April 28, and May 12, respectively). The sampling time is fixed at 17:00 on the sampling date, because the NSC concentration in the leaves fluctuates with the photosynthesis activity28. 4–6 seedlings were selected per treatment at each sampling. They were divided into needles, stems, coarse roots (root diameter > 2 mm) and fine roots (root diameter < 2 mm), and washed separately with distilled water. The seedlings were placed in an oven at 105°C for 30 minutes to prevent enzymatic carbohydrate reactions, and then dried at 80°C to a constant dry weight. For the NSC measurements, all dry samples were ground with a grinder until they could pass through a 60 mesh screen. The NSC was divided into soluble sugars and starch measurements, with the sum equal to the NSC content. Weighing 0.300 g of each dry sample for the determination of soluble sugar and starch. Both soluble sugar and starch contents were determined by the phenol colorimetric method29. Firstly, the weighed dry sample was put into a 10 mL centrifuge tube, 4 mL of 80% ethanol was added and centrifuged at 80℃ for 30 min at 3000 r/min for 10 min, the supernatant was poured into a graduated test tube, the residue was added to 2mL of 80% ethanol and the extraction was repeated twice, the supernatant was decolorized by activated charcoal at 80℃, and the soluble sugar was determined by anthrone colorimetry. The absorbance at 625 nm was measured by anthrone colorimetric method, and the sugar content in the extract was obtained according to the standard curve; secondly, the precipitate after the extraction of soluble sugar was pasteurized with distilled water for 15 min, and then extracted with 9.2 mol/L perchloric acid 2 mL for 15 min and centrifuged, and the supernatant was collected. The absorbance at 625 nm was measured by phenol colorimetric method, and the sugar content in the extract was obtained according to the standard curve, and the obtained sugar content was multiplied by 0.9 as the actual starch content after deducting water when calculating the starch content.
Statistical analysis. All statistical analyses were performed with SPSS 19.0 (IBM). Values are expressed as mean ± standard deviation, with 4–6 replicates per treatment. The data were tested for normality and homogeneity of variance prior to analysis of variance. Two factor analysis of variance (ANOVA) was used to analyze the effects of drought intensity (CK, LD, MD, and SD) and drought duration (March 14, March 29, April 13, April 28, and May 12) on the NSC (sugar and starch) content in each tissue. The differences in soluble sugars, starch, NSC content and biomass were examined by one-way ANOVA with different water treatment levels and different treatment times, and Duncan's method was used for multiple comparisons, at a 0.05 level of significance. The graphs were plotted using Origin (2021 Edition, Origin Lab Company,US).