Abstract

The availability of suitable native plant species for local animal husbandry development and ecological restoration is limited on the Qinghai-Tibetan Plateau. Therefore, comparisons of the ecological adaptability of native species to alternative habitats and their introduction into new habitats are of high importance. This study is aimed at identifying the alteration in morphological and physiological characteristics by measuring photosynthetic physiology, nutrient content, and growth associated with adaptation of plants to conditions at different altitudes 2450, 2950, 3100, and 3300 m above sea level (a. s. l.) on the plateau. Seeds of the dominant grass, Elymus nutans, were collected from locations at these altitudes and grown at a test location of 2950 m a. s. l. Results indicated that altitude had no significant effect on plant height and root depth. However, the leaf area and total root surface area of plants derived from 2950 and 3300 m a. s. l. showed a parabolic response, being greater than those of plants derived from the lowest (2450 m) and highest (3300 m a. s. l.). Total (root plus shoot) dry matter reduced progressively from 2450 to 3300 m a. s. l, while root : shoot ratio increased progressively with altitude. Seed yield of plants originating from the test altitude (2950 m a. s. l) was significantly higher than at any other altitude, being 20% lower at 2450 m, and 38% and 58% less in populations originating from the higher altitudes (3100 and 3300 m a. s. l.). There was also a parabolic decline in response of Elymus nutans germplasm from 3100, 3300, and 2450 m, compared with plants from 2950 m a. s. l., to photosynthetic rate, total N, soluble sugar, and starch contents. Germplasm from 2450 m a. s. l. had significantly lower shoot and higher root carbon content, lower shoot nitrogen, and lower root carbon-to-nitrogen ratio compared with plants derived from the other three altitudes. It is suggested that the stable, genetically determined morphological and physiological features of ecotypes showed parabolic responses which means these ecotypes have become adapted to local habitats, whereas parameters such as dry matter, total root : shoot ratio, photosynthetic rate, and intercellular CO2 concentration of plants reflected phenotypic linear response to current abiotic conditions. It is postulated that introduced ecotypes from 2450, 3100, and 3300 m could adapt to the environment at 2950 m a. s. l. gradually. We conclude that the increased thermal regime experienced by plants introduced from high altitude to low altitude may facilitate the increased growth of Elymus nutans subtypes. It is important to preserve local strains of native species, or ecotypes, for reintroduction into degraded environments and to maintain the greatest ecosystem stability in the northeastern Tibetan Plateau.

1. Introduction

Plants have the ability to alter their morphological and physiological traits in response to environmental variations [1, 2] and adjust the expression of these traits to accommodate their adaptability across multiple environments [3]. Altitude can significantly influence plant growth, structure, function, and metabolism [4, 5]. The morphological and physiological characteristics of alpine species reflect their environmental adaption under certain altitudinal gradients [6]. As a consequence, modern plant ecologists have paid much attention to these plant characteristics in different biological and ecological zones to better understand their adaptive mechanisms [7, 8].

Alpine plants exhibit specific morphological, physiological, and biochemical responses to an increasing altitude, such as size and increases in the leaf thickness [911]. Decreases in plant growth and a higher proportion of plant biomass in the roots usually result in a higher leaf N content and lower C : N ratio, which lead to a higher photosynthetic capacity in alpine plants [12]. Moreover, the relationship between plant functional traits and environmental factors has been recently gaining increasing attention. Plant species may have an optimum altitude for production of biomass and net photosynthetic rate; their enzymatic activities can increase or decrease upon deviation from this optimum altitude [1315]. Grasses exposed to high altitudinal areas with extreme climatic factors (cold, frost, drought, salinity, low oxygen, high wind velocity, intense UV radiation, and geographical differentiation) provide opportunities to study such adaptive mechanisms [16, 17]. However, the adaptation mechanism in response to altitude variation and climate change is largely unknown [18].

Elymus nutans is a native perennial grass and plays an important ecological role in maintaining regional species diversity and ecosystem stability. It grows extensively across the Qinghai-Tibetan Plateau and Xinjiang of Northwestern China where it is well adapted to a range of environments and is extensively grazed by livestock [19, 20]. It is also one of the important pasture species sown to increase livestock production and groundcover by perennial species. However, research on the species is not equal to the scale of its resource utilization in China. Most reports [2124] have focused on the biology, genetic diversity, nutrition value, and yield of the species, but there has been little information on the morphological and physiological characteristics of E. nutans and the species’ local adaptation to variations in altitude [25, 26]. An understanding of these characteristics, particularly of the plant morphological attributes and variations in photosynthesis, is essential to help identify the mechanism of plant adaptability to local-scale environments. Our basic question was whether the physiological and ecological traits of plants of the same species grown in different habitats differed in their responses when grown in alternative environments and how they adapted to their new environment. We hypothesized that there was a set of significant differences in plant morphological, physiological, and biochemical responses of E. nutans to increasing altitude that could be taken as a surrogate for differences in thermal regime, solar radiation, and precipitation. Basing on previous work, we attempted to answer the following questions: (1)How do the plant functional traits, photosynthesis, and dry matter respond to altitudinal variation(2)What are the differences in physiological plasticity of the different physiological parameters of E. nutans in the alpine area

We selected changes in shoot and root length, leaf area and root area, dry matter and seed yield, soluble matter concentration, the ratio of carbon and nitrogen, and photosynthetic parameters to elucidate the ecophysiological processes. Because these processes may enable this plant to adapt to varying environments, we could use them to evaluate the potential of this species expanding into higher altitudes. It is expected that findings from this study will help to identify traits for plant breeding strains of E. nutans that may be increasingly adapted to climatically changing habitats and for the development of the management strategies to rehabilitate degraded pastures both for pastoral grazing production and environmental sustainability.

2. Materials and Methods

2.1. Experimental Material

Healthy E. nutans seeds were obtained from wild plants growing in Tianzhu county of Eastern Qilian Mountains (102°23~102°78 E, 37°11~37°18 N), Qilian Mountains, China, which lies in the northeast margin of the Tibetan Plateau, a key site to understanding the species’ adaptation in the northeastern of Tibetan Plateau. At least 15 plants of Elymus nutans per site were cut from four sites (Dachaigou, Jinqianghe, Honggeda, Daiqian) selected randomly from habitats at 2450, 2950, 3100, and 3300 m a. s. l., in an area of 100-1000 m2 that was ecologically homogeneous for the target ecotypes (see Figure 1).

The vegetation type in the sampling areas was composed of alpine meadow grasses and forbs. The soils on which each ecotype grew were of the same soil taxonomic category being shallow, postglacial subalpine meadow soils of high organic matter content. Livestock grazing is the only land use. Elymus nutans is the dominant species in this area and occurs naturally and abundantly at altitudes between 2000 and 5000 m a. s. l. in the Qinghai-Tibetan Plateau. After collection, the samples were bulked, dried to 30°C, and seeds threshed out. Seeds were stored at a constant 4°C temperature at the Seed Centre of the College of Grassland Science, Gansu Agricultural University, Lanzhou, China.

2.2. Site, Climate, and Design

Seeds of the 4 ecotypes were sown at the Tianzhu Alpine Grassland Experimental Station of Gansu Agricultural University at an altitude of 2950 m a. s. l., the northeastern margin of the Qinghai-Tibetan Plateau, China. The climate is alpine, cold, and semiarid, with thin air and low temperature, and strong solar radiation, where there is no frost-free season, only cold and cool seasons. The annual average temperature is 0.1°C, with the highest average temperature in July (mean 12.7°C) and the lowest in January (mean -18.3°C). The average annual precipitation is 328.7 mm, concentrated from June to September. The 4 ecotypes were sown as a randomized complete block design with 3 replicates and the plot size was with a 1 m gap between plots. Each plot has row spacing of 40 cm. Sowing depth was 3-4 cm and sowing rate was 22.5 kg/ha. The soil was a shallow, postglacial, subalpine scrub grassland soil with a higher organic matter content (9.45%) in the 0~20 cm soil layer.

2.3. Measurements
2.3.1. Plant Height and Root Depth, Root Surface Area, and Volume

Establishment and vegetative growth of E. nutans is slow at altitudes above 2000 m a. s. l., and plants only reach flowering and maturity by the second year. The first year’s shoots die back and the growing points lie below the soil surface during the first winter, to reemerge in the second spring [27]. Therefore, all the data presented are from the second year of planting.

Selected measurements were made at different development stages. The plant morphological parameters (shoot height, root length, surface area, and volume) were measured during the period of rapid shoot extension in the second year (the end of July), flag leaf area was measured at full emergence of the flag leaf before flowering (the early of August), and the physiological parameters (soluble sugars, starch, total carbon and nitrogen, and gas exchange properties) and total aboveground (shoot) and belowground (root) dry matter were measured at flowering (August 10th).

Plant height was measured by recording the height of 20 randomly selected plants of E. nutans within a plot before each harvest with a 2 m ruler from the soil surface to the highest point of the plant and averaged to give the mean for each treatment plot. Sections of rows 33 cm long were then randomly selected, whole plants were dug out, and root depth was measured. Whole plants were stored in the self-sealing, labeled bags and refrigerated. Cleaned roots were put into the root plate submerging the root system, and root surface areas and root volume were measured using Photo Scanner and the data analyzed by Root Analysis software.

2.3.2. Leaf Areas

The areas of 10 flag leaves from randomly selected plants taken when flag leaves had fully emerged, and at flowering were determined using an AM-300 Area Meter (Analytical Development Company, Herts, UK).

2.3.3. Herbage Yield

Individual treatment plots were harvested for herbage yield at fully flowering period, with subsamples then taken from each plot, freshly weighed and oven dried at 65°C for 48 hours before being ground.

2.3.4. Seed Yield

Seed yield was made at 80% seed maturity (seeds drop easily) in mid-September. Seed collection was carried out from randomly selected plots of each ecotype, and the seed yield of each experimental plot (10 m2) was determined and converted into kg/ha.

2.3.5. Total Soluble Sugar

The physiological traits were all measured when plants were flowering. Soluble sugar was extracted from a 0.05 g sample with 1-4 ml 80% () aqueous ethanol in a water bath at 80°C for 40 min. This sample was centrifuged for 30 minutes at 12000 r/min the supernatant collected, the residue reextracted with 80% () aqueous ethanol, and finally, the sample was centrifuged and merged with the supernatant. The extracted solution was diluted to 50 ml after extracting 2 times. 1 ml of the extract obtained was mixed with 5 ml of anthrone-sulphuric acid reagent. After cooling to room temperature, the absorbance was measured at 627 nm using a UV-2550 UV-visible spectrophotometer.

2.3.6. Starch Contents

The residue from the soluble sugars was boiled with 6 mol/l HCl in a water bath at 80°C for 5 min, the sample centrifuged for 10 minutes at 4000 r/min the supernatant collected, and the above step repeated twice. Absorbance was measured at 580 nm using a UV-2550 UV-visible spectrophotometer. Starch content was calculated by multiplying by the coefficient of hydrolysis (0.9).

2.3.7. Determination of Total Organic Carbon Content

The samples were air-dried and ground to pass a 0.25 mm mesh for total organic carbon content. Approximately 0.015-0.020 g of each leaf and root sample was put into clear flask, mixed with 10 ml of K2Cr2O7, and then heated 5 min. The solution was cooled and then titrated with FeSO4 [28].

2.3.8. Nitrogen Analysis

N was determined by a small-scale Kjeldahl digestion [29]. 0.5 g of ground herbage from each sample was added to a flask, with 10 ml of H2SO4 and 3 g of CuSO4·5H2O and K2SO4 (CuSO4·5H2O : ), and shaken until thoroughly mixed, heated from a low to high flame until the solution became clear, topped up to 100 ml with ammonium-free distilled water, and 10 ml of the solution mixed with 4 ml 40% NaOH and distilled using a Kjeldahl apparatus for 4 min. NH3 released through distillation was absorbed by 1% boric acid and then titrated against 0.02 mol l-1 HCl to determine the N content.

2.3.9. Gas Exchange

The Li-6400 portable photosynthetic analyzer (Li-COR, USA) was used to determine the daily change (8:30 am-17:30 pm) of photosynthetic parameters of E. nutans in an open gas analysis system. All measurements were carried out on healthy leaves in the field on sunny days at flowering period. Three leaves from each plot were randomly selected. The photosynthetic parameters included a net photosynthesis rate (Pn, mmol CO2 m-2 s-1), transpiration rate (Tr, mmol H2O m-2 s-1), stomatal conductance (Cond, mmol m-2 s-1), and intercellular CO2 concentration (Ci, μmol CO2 mol air-1). The parameters were calculated as the mean value for the day.

2.3.10. Statistical Analyses

Normality of distribution and homogeneity of variance were checked before any further statistical analysis. An investigation into normality of the data was carried out by performing the Kolmogorov-Smirnov test using SPSS 19 software. The data were normally distributed and better meet the assumptions of the statistical tests used. Analysis of variance and calculation of means were conducted using Excel and SPSS19 software. Morphological data and photosynthetic responses of plant were analyzed using a one-way analysis of variance (ANOVA), and the mean values were compared using the least significant difference (LSD) at .

3. Results

3.1. Changes in the Plant Height and Root Depth of Four E. nutans Populations from Different Altitudes

The plants originating from different altitudes showed no statistically significant differences in height and root depth when grown at 2950 m a. s. l. (Figure 2), which indicated that the altitude variations had less effect on the growth of plant height and root depth in alpine zone.

3.2. Changes in Leaf and Root Surface Area

Leaf and root area responses varied significantly along the altitudinal gradient, increasing from 2450 to 2950 m and then decreasing from 2950 to 3300 m a. s. l. Both of leaf areas and root surface areas were significantly different and smaller in value for plants collected from both the highest and lowest altitudes compared with those from 2950 and 3100 m a. s. l. (Figure 3). The leaf area and root surface area showed a significant negative correlation with altitude, and the correlation coefficient was 0.9938 and 0.9393, respectively.

3.3. Changes in Total Root Length and Root Volume

Original altitude significantly influenced total root length and root volume but to different extents (Figure 4). There was no significant difference in total root length from 2450 to 2950 m but root length decreased as altitude increased from 2950 to 3300 m a. s. l. The root length from 3300 m was lower than that from 2450 m and 2950 m a. s. l. 41% and 46%, respectively. Root volume increased as altitude increased from 2450 m to 2950 m a. s. l. and then decreased. The root volume from 3300 m was lower than that from 2950 m a. s. l. 47%.

Figures 3 and 4 show that the smaller leaf areas and root characteristics of plants from the lowest altitude (2450 m a. s. l.) when grown at higher altitude, and the inability of the plants from higher altitude (3310 and 3330 m a. s. l.) to increase their leaf and root size when grown at a lower altitude, provide good evidence that local adaptation has occurred, inhibiting the plants from taking advantage of the conditions prevailing at 2950 m a. s. l.

3.4. Changes in Aboveground and Belowground Dry Matter and Seed Yield

The aboveground dry matter decreased with altitude increased while the belowground dry matter showed a reverse trend. Ecotypes from high altitude allocated significantly () greater dry matter yield to roots than those adapted to low and mid altitudes. Aboveground dry matter yield allocation averaged 875 kg/ha at the lowest altitude, nearly twice the value of the plants from the highest altitude. It was observed that ecotypes from lower altitudes allocated more dry matter to the stems of plants than those from higher altitudes. Total dry matter (combined root and shoot dry matter) values for the four ecotypes showed that there was a reduction with each increase in original elevation (Table 1), while root : shoot ratios also increased, providing a clear indication that the total abiotic stress increased with altitude, as root : shoot ratio increases with nearly every abiotic stress. The plants selected from, and grown at 2450 m a. s. l., were the largest in total weight, while the plants from the 2950 m a. s. l. were only slightly smaller, but they partitioned more of their assimilates into roots or seeds and less into dry matter. Plants from 3000 m a. s. l. were consistently smaller in every respect and were not able to take advantage of the milder conditions when grown at 2950 m a. s. l.

Seed yield of plants from 2450 to 2950 m a. s. l. increased 22.4 kg/ha and then decreased markedly and reached a minimum at the 3300 m (46.2 kg/ha), with the seed yields from 3100 and 3300 m a. s. l. being significantly lower than either of the lower altitudes. While seed production of the 2450 m a. s. l. ecotype was somewhat restricted by being grown at the higher elevation, the lack of capacity of the 3100 and 3300 m ecotypes to exploit the more favourable growing conditions at 2950 m a. s. l. again indicates that substantial local adaptation has taken place.

3.5. Soluble Sugar and Starch Content

Soluble sugar and starch content all showed significant differences in the leaves and roots of plants derived from different altitudes. The soluble sugars followed the similar parabolic curve to the root and shoot morphological traits increasing first with altitude from 2450 to 2950 m a. s. l. and then decreasing (Table 2). Soluble sugars of belowground roots reached the highest value at 3100 m a. s. l. The aboveground starch contents were only marginally different from plants originating from 2950 m a. s. l., but the belowground starch contents of plants from the highest and lowest altitudes were significantly less than those from the intermediate altitudes.

3.6. Total C and Total N Concentration and the Ratios of C : N

Table 3 shows that aboveground C contents were not significantly different except in plants derived from 2450 m, although plants from 3100 and 3300 m a. s. l. had numerically greater C contents. The root C was similarly higher in the plants from the lowest altitude. There were no significant differences in root C content among ecotypes from high altitudes. The plants derived from lowest altitude had a significantly lower aboveground nitrogen (N) content and higher belowground N content than the other ecotypes. The aboveground C : N increased as the altitude increased, but the root C : N ratio was only significantly different in the 2450 m a. s. l. ecotype.

3.7. Photosynthetic Parameters of Leaves of Four E. nutans Populations from Different Altitudes

Photosynthetic characteristics in leaves of the four ecotypes from different altitudes were measured during the period of maximum growth (the end of July, second year). The results (Figure 5) are the means of data collected between 8:30 am-17:30 pm during 25 July. They showed the photosynthetic rate (Pn) and intercellular CO2 concentration (Ci) responded in a similar manner to the shoot dry matter. Pn decreased as the altitude increased. The Pn of plants from 3300 m was 28% lower than those from 2450 m a. s. l., but there was no significant difference among values of plants from 2950, 3100, and 3300 m a. s. l. However, Tr increased as the altitude from which the plants were derived increased. The Tr from 3300 m was 43% higher than that of 2450 m a. s. l. (Figure 5(a)). By calculating the water use efficiency (), it was found the plant WUE gradually decreased from lower altitude to higher altitude. The Cond and Ci responses were highly correlated with the altitude of the original populations and decreased significantly as altitude increased (Figure 5(b)). The Ci and Cond from 3300 m were lower than that of from 2450 m a. s. l. 8% and 48%, respectively.

4. Discussion

4.1. Variation in Morphological Traits of E. nutans with Altitude Gradient Migration

Altitude is an important ecological factor that affects population differentiation and genetic variation. The high altitude especially in the alpine region has a pronounced influence on the plants. Plants have obvious adaptation strategies by modifying their morphological and physiological attributes [30]. The long-term effects of environment on plants can alter their morphology, physiological function, and ecology [31, 32]. Only those that adapt to their environment can survive well, while those that do not will gradually be eliminated from the local ecosystem. Therefore, morphological characterization of ecotypes provides useful information for evaluation of agronomic aptitude and development of new varieties [33].

In most alpine species, small stature is a genetic adaptation to low temperatures, and hence, they retain this feature when grown in a warm environment as low temperature affects the growth process and photosynthetic carbon gain [34]. In this experimental location, plant height, root depth, and root volume of E. nutans from different altitudes showed no significant change. These were the relatively conservative traits among the morphological traits measured along altitudinal gradients in an alpine region, similar to those described by Zhu’s research [35]. However, ecotypes from 3300 m had larger leaf areas and root surface areas than those from 2450 m a. s. l. In the cold temperate zone, an increase in temperature reduces abiotic stress and improves the overall plant growth rate [36], especially affecting the photosynthetic and nutrient absorption processes. Our results are in line with other research findings for trees [37]. Furthermore, the plants derived from 2950 m a. s. l. and grown at the same altitude had the greatest morphological size in attributes such as plant height, root length, leaf area, total root length, and root volume. This indicated that local ecotypes performed better than the introduced plants from other altitudes because they were better adapted to local environmental conditions [38]. In addition, different ecotypes of the same species may adopt different survival strategies to adapt to the environmental stress [39]. For plants growing in the cold environments, for example, the leaf area and leaf size are normally small but the corresponding plant dry matter accumulation and net photosynthesis rate may be more profoundly affected [40]. The physiological measurements relating to photosynthesis taken at the period of fastest daily growth (prior to flowering) showed that the plants derived from 2950 and 3100 m a. s. l. were growing faster, which supports this interpretation. The plants from the highest altitude (3300 m a. s. l.) were not able to adjust to the warmer growing conditions within a single generation, and the plants from the lowest altitude (2450 m a. s. l.) were probably restricted by the lower temperatures they encountered at 2950 m a. s. l. [41].

4.2. Effect of Altitude Migration on Material Accumulation and Distribution

Although the four ecotypes of the species E. nutans were generally similar in appearance and morphological characteristics, differences were observed in the dry matter characteristics, carbon, and nitrogen content. Plants from the lowest altitude maintained higher dry matter accumulation and seed yield than those adapted to 2950 m a. s. l. This is good evidence that there had been local adaptation to the conditions occurring at different altitudes [42]. We could expect that plants selected from the lowest altitude and grown at a higher altitude would not grow so fast because they were exposed to more cold, and that the plants from the higher altitudes grown at lower altitude would grow faster (and perhaps accumulate more belowground dry matter) and also through regulating physiological activity to adapt to the environment changes. When we compared (Table 1) the way the plants from different altitudes partitioned their assimilates, it was clear that the total amount of plant production (aboveground and belowground dry matter and seed yield) was significantly less in the plants from the higher altitudes, despite being grown at a more favourable (warmer) altitude, whereas the plants from the two lower altitudes had a similar total biomass but allocated their assimilates in different ways. The total amount of aboveground herbage was greatest in the plants from the lowest altitude (and probably if all ecotypes had been grown at 2450 m a. s. l., these plants would have been largest overall). The plants originating from the low altitude may have had some differences in morphological and physiological characteristics that could possibly have been due to the results of natural selection not just the altitude factor [43]. However, overall, each ecotype largely kept its original characteristic when being transplanted to another altitude [44]. As we found a linear decrease in seed set with altitudinal origin, it is likely that the higher growth rates of plants from lower altitudes would be maintained in other growing sites.

4.3. Effect of Altitude Migration on Physiological Characteristics

Nonstructural carbohydrate (NSC) which mainly refers to the sugar and starch has a great influence on plant growth and adaptability [45]. NSC is the product of photosynthesis. Therefore, factors affecting photosynthesis will also affect NSC of plants [46]. When photosynthesis changes, the NSC may change a variety of ways that can be induced directly or indirectly by various environmental factors [47, 48]. Previous studies have shown that photosynthetic capacity and plant N content may have a positive correlation while higher photosynthetic activity can increase CO2 assimilation so that NSC is also enhanced [49, 50]. In our study, plants from the middle altitudes (2950 and 3100 m a. s. l.) had higher soluble sugars, and root starch contents, in the typical parabolic response form, whereas the carbon and nitrogen contents and C : N ratios were significantly different only in the plants derived from the lowest (2450 m a. s. l.) altitude suggesting they were partitioning their assimilates into structural components in different ways. The seedlings sourced from higher altitudes had better nitrogen utilization ability than those from lower altitudes. This suggests that ecotypes of Elymus nutans sourced from different altitudes may have beneficial attributes of adaptive value under the condition of future climate change.

We may postulate that plants from the lower elevations (below 3000 m a. s. l.) are therefore most suited to being sown in situations where the amount of grass for yak feed is of most importance, whereas ecotypes that have developed at higher specific altitudes are most suited for use in rehabilitation programs aimed at repairing degraded natural grasslands that have suffered from erosion or overgrazing.

5. Conclusions

We can conclude that local adaptation has occurred in Elymus nutans growing naturally at different altitudes: these adaptations include changes in the size of the plants, the ratio of roots to shoots, the ratio of C : N, the total herbage and seed produced over the growing season, and the extent of the genetic by environment interaction (responsiveness) when grown at intermediate elevations. We suggest that ecotypes from the lower altitudes will be most productive for increasing grazing by yaks below 3000 m a. s. l., but that locally adapted ecotypes should be used within their specific local habitats for ecological restoration of degraded grassland. Adaptation is the result of a long-term evolutionary process in concert with stress conditions. Further studies are needed to fully understand how altitude migration affects Elymus nutans growth on the Qinghai-Tibetan Plateau.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare there are no competing interests.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (grant number 31660684), the Key Laboratory of Superior Forage Germplasm in the Qinghai-Tibetan Plateau (grant number 2020-ZJ-Y09), and the Earmarked Fund for China Agriculture Research System (grant number CARS-34). The authors are grateful to Changlin Xu for his help and valuable technical assistance and those students who were involved in part of the field and lab work.