Combined Treatment of Salinity Stress and Fruit Thinning Effect on Tomato

This was an experimental investigation of the combined treatments of salinity (SAL) stress and fruit thinning (FT) on the growth, yield, fruit quality, and water use efficiency (WUE) of tomatoes with non-soil cultivation. The experiment was carried out in a plastic tunnel, Japan. Tomato (Solanum lycopersicum) cv. Momotaro seedlings were transplanted in a randomized complete block (RCB) manner with six plants/treatment, and an overall 36 plants in 18 pots (2 plants/pot). The experiment involved varying SAL treatment (no-SAL, moderate SAL, and serious SAL, with electroconductivity of 0.8, 3.0, and 4.5 dS m−1, separately) and FT treatment (NT: no thinning and 3FT: three-fruit treatment). The tomato growth, yield, and WUE were significantly suppressed with increasing SAL. In comparison, FT treatment had less effect on tomato growth and water consumption. Either SAL stress or FT treatment significantly improved fruit quality. The combined treatment proved better than single treatment of either SAL stress or FT, avoided the subsize fruit following SAL stress treatment, reduced fruit cracking found with FT treatment, and greatly improved fruit quality. The SAL thresholds of WUEs in relation to biomass, yield, and marketable yield were approximately 3.0 dS m−1 under these soilless conditions. Path analysis showed that biomass and water consumption were important indexes affecting yield. Logistic equation fitting showed that SAL stress tended to inhibit and delay plant growth; however, FT tended to advance and shorten the period of plant growth.


INTRODUCTION
Global urbanization was accelerating, and it was reported that 60% of the global population will live in urban areas by 2030 (1). People pay more and more attention to urban food supply and food security (2,3). As a result, cities around the world are integrating local food production capacity into the built environment to create food security and environmentally sustainable urban agriculture (4). The home garden (also known as household gardens, kitchen gardens, balcony gardens, or homestead gardens) is one form of urban agriculture that represents critical spaces in the configuration of urban socio-ecological landscapes (5,6). However, the development of home garden was restricted by factors, such as small planting scale, poor professional farming personnel, and low input-output ratio (7).
In recent years, soilless cultivation was used in the home garden to optimize resources and ultimately increase yields (8)(9)(10). Soilless culture is broadly adopted to uplift the regulation of the growth environmental conditions and prevent the indeterminacy of the soil (11). It avoids as well the accumulation of salinity (SAL), pests, and illnesses and realizes minimal environmental pollution from irrigation and fertigation (12,13). Tomato (Solanum lycopersicum) is one of the most important vegetable plants in the world. Global production exceeds 180 million metric tons, with China and the USA as the leading producers in 2019, and production in China accounts for 20% of the world (14). Tomato has moderate tolerance to SAL (1.3-6 dS m −1 ) (15), and greater levels of EC reduce yield and fruit size (16,17) and water and nutrient uptake (18). The negative effects on tomato plant growth are due to SAL stress-causing reductions in root cell growth (19), leaf expansion (20), leaf chlorophyll (21), and plant photosynthesis (22). The decreased output at SAL over the liminal value is due to reductions in the number of fruits produced (19,23) and in fruit size (24). However, many studies have confirmed that tomato fruit's total soluble solids, titratable acidity, and lycopene, fructose, and glucose concentrations are increased with higher SAL (25,26).
The studies also showed that yield, fruit quality, and fruit size were affected by numerous reasons, such as plant community (27,28), fruit thinning (FT) (29,30), pruning, and variety choice (31). FT is adopted for the limitation of the fruit quantity on each girder and decreases competitive activities to elevate fruit weight. The rise of overall numbers of flowers and fruits can bring more competitive activities for photosynthate and hence decreases fruit weight (29). Abdel-Razzak et al. (32) pointed out that controlling the numbers of flowers, fruits, or fruit girder efficiently decreases competitive activities between fruits, so that extra assimilates are diverted to fewer fruits per girder. FT before harvest decreases fruit total yield per plant (33), but increases fruit size and average fruit weight (32,34). However, FT was confirmed to accelerate enlargement and increase the size of the retained fruit, thus increasing the risk of fruit cracking and reducing marketable yield (35). In addition, some researchers reported that lower fruit load could increase total soluble solids, VC, titratable acidity, and total carbohydrate of tomato fruit (32).
The environmental thermal status remarkably affects the tomato physiological process (36,37). Plant growth dynamic models provide the possibility of evaluating the dynamic effects of environmental factors on plant growth. Logistic models are the most commonly used plant dynamic growth simulation models and can be used to depict the association between plant development and thermal unit cumulation (like growing degree days) (38).
According to the above explanation, soilless cultivation is widely used in the home garden because of its convenience and efficiency. SAL stress reduces tomato fruit size and marketable yield, but increases fruit quality. FT increases fruit size and quality, but increases the risk of fruit cracking and so reduces marketable yield. The above two cultivation methods show contradictory advantages and disadvantages for tomato growth. Therefore, it is necessary to determine an appropriate combined treatment of SAL stress and FT for tomatoes. However, the previous research primarily highlighted SAL stress or FT. There has been few considerations of combined treatment, and few relevant studies have used soilless cultivation. Consequently, our research aimed to explore (1) the combined treatment of SAL stress and FT on tomato growth, fruit quality, and WUE under non-soil cultivation and (2) the kinetic features of critical growing indices impacting the plant as per the logistic growing pattern.

Experiment Site
The assay was carried out inside a plastic tunnel (11 m × 5 m, length and width) at the Experiment Farm of Gifu University, Japan (35 • 27 ′ 51 ′′ N, 136 • 44 ′ 14 ′′ E). The experiment duration was for 12 weeks from April to July 2018, with average daytime and nighttime temperatures of 28.1 and 20.4 • C, respectively, and corresponding average humidity of 51.5 and 79.7%.

Experimental Design and Treatments
In the experiment, tomato (Solanum lycopersicum) cv. Momotaro plants were grown via a purchasable hydroponic kit (Minoru Kasei, Co., 388-1 Shimoichi, Japan), which includes a nutrient solution bucket, a plant table, and a plant bowl (Figure 1). The normal cultivation liquid Vegetable Life A (1.3% N, 1.2% P, 1.9% K, 0.32% Mg, 0.008% Mn, and 0.008% B; Otsuka Chemical) was desaturated 200-folds by tap water for the following planting of tomatoes. The desaturated liquor carried an electroconductivity of 0.8 dS m −1 . NaCl was supplemented into the liquor to bring SAL stress to the plants.
Initially, the seedlings were planted at the above-mentioned farm. If they were nearly 20 cm high, then they would be transplanted in an RCB manner as described in the abstract above. All the pots were full of 14 L of desaturated liquor.
The previous study of our team revealed that the SAL liminal value of tomato registered 1.41 dS m −1 via the normal cultivation liquid (21). Therefore, SAL levels of 3.0 dS m −1 (moderate stress) and 4.5 dS m −1 (serious stress) were employed herein. A total of six treatments were utilized, with three electroconductivity  Table 1). In addition, CN, MN, and SN were the treatments without FT; CF, MF, and SF were the treatments in which the truss was pruned to have three remaining fruits. Fruits were pruned when they were marble size (31) or when they were at light-maturity stage (32). Axillary buds were removed at sprout, and the apical bud was removed in the 5th week (28-35 days posterior to transplantation) when the plant was in the 3rd flower cluster stage. Considering the variation in the climate environment in the house, the location of the pots was changed after irrigation. The cultivation liquid was renewed once every 2 weeks. SAL stress began on May 22 and the first fruits were pruned on May 24.

Transpiration
The transpiration (TP, g day −1 ) of each pot was determined via gravimetry, with the following calculation: in which TP denotes the transpiration when the liquor is renewed (g), W a denotes the weight of the fresh liquor (g), W ′ d denotes the weight of the replaced liquor (g), and W ′′ i denotes the weight of watering (g). Watering was carried out supplementing cultivation liquid once every 1 or 2 days. The overall TP of each pot was divided between the two plants in every pot as per their eventual phytomass, and the overall TP of each plant is deemed as the plant's actual TP. The electroconductivity of every pot was surveyed posterior to watering. The solution or NaCl was supplemented when the liquor electroconductivity diverges from the normal values of 3.0 ± 0.2 and 4.5 ± 0.2 dS m −1 .

Plant Growth Parameters
The biomass (g) of every plant was identified posterior to the assay. Leaf chlorophyll was surveyed once a week via a SPAD chlorophyll meter (Minolta). A linear association existed between SPAD results and extracted leaf chlorophyll (39,40); hence, those were employed to reflect the leaf chlorophyll in our study. The dry matter of the superior parts was determined via desiccating the plant material at 105 • C for half an hour and then at 70 • C in the ventilation stove till a steady weight (41).

Yield and WUE
Red fruits were harvested two times a week in harvesting periods. The yield (g plant −1 ), the number of flowers, and the number of fruits on every plant were surveyed in harvesting periods. The harvested fruits were divided into marketable or unmarketable. Marketable ones were divided as per size (≥30 mm) and appearance (no cracking). The sizes were measured by vernier caliper. The marketable fruit number and the marketable yield on every plant were also surveyed. The fruit set (%: number of fruits on every plant as a percentage of the total number of flowers) and marketable fruit set (%: number of marketable fruits on every plant as a percentage of the total number of flowers) on every plant were calculated.
The WUE was calculated according to the biomass, yield, and marketable yield, with the following calculations (42): in which WUE b (g kg −1 ) denotes grams of biomass generated per kg of water; WUE y (g kg −1 ) denotes grams of yield generated per kg of water; and WUE my (g kg −1 ) denotes grams of marketable yield generated per kg of water.

Fruit Quality Parameters
Fully matured tomato fruits were used for the fruit quality measurements. A total of eighteen fruits were selected for every treatment, for an overall 108 fruits for the measuring. Fruit firmness was evaluated with a penetrometer. Fruits that showed symptoms of cracking were separately counted to estimate the fruit cracking ratio. The juice of every fruit was extracted via a squeezer, and the sugar content (%) and acid content (%) were determined with a Pocket Brix-acid meter (PAL-BX/ACID1, ATAGO Co., Ltd., Tokyo). Moreover, taste index was computed via the formula of Hernández-Suárez et al. (43):

Path Analysis
Path coefficient analysis divides the components of direct and indirect traits and provides information on direct and indirect effects of interrelated components on yield. Path analysis of plant-growing factors (e.g., root fresh weight, quantity of clusters per plant, individual fruit weight, number of fruits per plant, and chlorophyll stability index) and the yield not only displays the correlation between these factors, but also presents the direct and indirect effects of growth factors on yield (38,44). In this research, biomass (X1), plant height (X2), leaf chlorophyll (X3), and water consumption (X4) as independent variates and yield (Y) as the dependent variate were selected to calculate path coefficients.

Logistic Equation
A logistic equation was used to quantify the effects of SAL and FT on the growth of tomatoes, which is expressed below (38,45): in which y denotes the dependent growing parameter or plant biomass cumulation (g plant −1 ); k denotes the top asymptotic line reflecting the superior limit of plant growth; a and b denote quotient of the preliminary stage and accretive rate, separately; and t denotes the independent runtime ( • C d).
Feature parameters of the logistic equation were computed below (45): where v m (g plant −1 • C −1 d −1 for plant biomass cumulation) is the maximal slope and is deemed as the maximal growing velocity in the fast-growing stage; t m ( • C d) denotes the relevant time of v m and denotes the time when growing velocity registers its maximal value; t 1 and t 2 ( • C d) denote the inflexion of the logistic equation and represent the nodal points of the start and end time of the fast-growing period, respectively; and t d denotes the duration of the fast-growing period ( • C d). Accumulated thermal unit is adopted in the experiment farm with a small meteorological station, rather than calendar days posterior to cultivation. The formula adopted to calculate thermal time stands thus: (12) in which GDD is the growing degree days; T max and T min denote the daily maximal and minimal temperature ( • C), separately; and T base denotes the basic temperature of 10 • C for greenhouse tomato (38).

Data Analysis
Statistical analysis was carried out via two-way analysis of variance, and the results were afterward in contrast to Duncan's test with a credit level of 5% via the R program language. The logistic equation was fitted using Origin 2019b.

Plant Growth
Plant biomass, dry matter, and chlorophyll significantly reduced with increasing SAL regardless of thinning condition ( Table 2). The plant growth variables exhibited no significant differences with varied FT conditions ( Table 2). The biomass and dry matter were no significant differences with combined treatment (SAL and FT), but leaf chlorophyll was significantly different. Meanwhile, the values of growth variables under MF and SF treatments (combined treatment) were remarkably smaller vs. those under the CN.

Yield
Yield variables significantly decreased with increasing SAL. Compared to NT, the 3FT treatment showed significantly lower fruit number, fruit set, marketable fruit number, and marketable fruit set. For combined treatment, these yield variables showed no evident differences except marketable yield ( Table 3).

WUE
In this study, the daily plant's actual TP of a plant was defined as the daily water consumption. Meanwhile, for the convenience of comparison with soil cultivation, the daily water consumption of a plant was expressed as "mm." The cumulative water consumption (CWC) before and after SAL treatment is shown in Figures 2A,B, respectively. The CWC did not differ among the treatments before applying SAL (Figure 2A) but significantly decreased with increasing SAL (Figure 2B). However, the CWC was not different among the same SAL conditions (i.e., CN and CF; MN and MF; and SN and SF) ( Figure 2B). The WUE for biomass, yield, and marketable yield showed an initial rise and a subsequent reduction with increasing SAL. The WUE variables exhibited no significant differences with varied FT conditions.  For combined treatment, these WUE variables showed significant differences ( Table 4).

Fruit Quality
The parameters for fruit quality (sugar, acid, and taste index) were significantly enhanced with increasing SAL ( Table 5). The values of fruit quality variables under 3FT treatment were significantly higher than for NT treatment. Notably, the sugar, acid, and taste index were higher under 3FT than NT treatment by 7.9, 12.5, and 10.1%, respectively. Fruit firmness decreased with increasing SAL, and FT resulted in significant softening of the fruit. The SAL stress significantly decreased the fruit cracking ratio, and FT significantly increased it. In addition, for combined treatment, despite no evident differences, the fruit quality variables under MF and SF treatments (combined treatment) showed the highest values among the treatments. The fruit firmness showed no evident differences with combined treatment, and the fruit cracking ratio was significantly different.

Path and Logistic Analyses
Yield showed the greatest value of variable coefficient (47.6%), followed by water consumption (35.7%) and biomass (32.4%), with the lowest value for plant height (4.0%) ( Table 6). The Fvalues between plant growth indicators (biomass, plant height, leaf chlorophyll, and water consumption) and yield showed very significant correlations and indicated that regression analysis could be performed ( Table 6). Path analysis showed that biomass had the strongest effect on yield, mainly based on the direct (1.523) and indirect effects of water consumption (−0.609) ( Table 6). Water consumption showed the second greatest effect on yield mainly based on direct (−0.634) and indirect effects of biomass (1.462). This was followed by plant height, mainly based on the direct (−0.103) TABLE 4 | Combined treatment on cultivation water consumption (CWC), water use efficiency for biomass (WUE b ), water use efficiency for yield (WUE y ), and water use efficiency for marketable yield (WUE my ).

Treatment
Water Con. (kg) WUE b (g kg −1 ) WUE y (g kg −1 ) WUE my (g kg −1 ) x * 1,000/~|*1,000/~}*1,000/C     The path analysis showed the significance of biomass on yield. As the logistic equation performed well in the kinetic description of plant growth with the parameters, it was fitted to the association between accumulated temperature and kinetic variations in tomato biomass. Biomass for diverse SAL and FT treatments was fitted via logistic equations (Figure 3). The root mean square error range was 26.47-60.44 g plant −1 with a mean of 49.01 g plant −1 , the index of agreement (d) range was 0.09-0.82 with an average of 0.56, and coefficients of determination (R 2 ) were 0.977-0.994 with an average of 0.987.
The values of t 1 , t 2 , and t m rose with increasing SAL ( Table 7). The start of the rapid growth period (t 1 ) for S 3.0 and S 4.5 was 6 and 12 days later than for S 0.8 , respectively, with a slope of 3.2 days (dS m −1 ) −1 (R 2 = 0.988). The end of the rapid growth period (t 2 ) of S 3.0 and S 4.5 was 4 and 28 days later than S 0.8 , respectively, with a slope of 6.39 days (dS m −1 ) −1 (R 2 = 0.784). The periods of most rapid growth (t m ) of S 3.0 and S 4.5 were 10 and 20 days later than S 0.8 , respectively, with a slope of 5.06 days (dS m −1 ) −1 (R 2 = 0.981). The maximal growing velocity in the fast-growing stage (v m ) dropped with increasing SAL.

DISCUSSION
The tomato plant growth variables decreased with increasing SAL ( Table 2). The decrease of biomass and dry matter is likely because SAL stress restrained plant organ (root, shoot, and leaf) development (20,46). Albacete et al. (47) showed that the fresh weight of tomato root was reduced by 30% posterior to 3 weeks

Treatment
Characteristic parameters under SAL conditions (100 mM NaCl). Kamrani et al. (20) found that SAL up to 20 mM affected tomato shoot development, and that elevated SAL significantly decreased seedling height. Oztekin and Tuzel (48) discovered that the mean seedling height of 21 purchasable tomato cultivars was reduced by 29.03% under 200 mM NaCl treatment in contrast to non-SAL treatment. Azarmi et al. (49) revealed that the overall foliage area dropped with rising SAL (EC: 2.5-6 dS m −1 ). The decrease of leaf chlorophyll content was associated with salt-caused increase in chlorophyllase activity, which casted an unfavorable effect on membranous stability and attenuation of protein-colorant-lipin complex (50). The decrease in leaf chlorophyll was consistent with the results of Azarmi et al. (49), which shows that leaf chlorophyll content was decreased with SAL. Plant growth variables did not significantly differ with varied FT conditions, which indicates that FT did not cast impact casted on the plant growth. Tomato leaves, stems, and roots are the prime source organs for assimilates, so FT had no effect on assimilate products. Although FT may reduce the consumption of assimilating products, the excess may be absorbed by other plant organs (51). In conclusion, FT did not affect the growth of tomato plants.
Although combined treatment (MF and SF) of SAL stress and FT remarkably reduced the biomass, dry matter, and chlorophyll, no difference was evident with varied FT conditions. This indicated that tomato growth variables were reduced due to SAL stress but not FT.
Tomato yield variables significantly decreased with increasing SAL. This conclusion was consistent with many other studies. Qaryouti et al. (52) discovered that the overall yield of tomato cv. Durinta F1 was remarkably decreased at SAL of 5 dS m −1 and above, with a 7.2% yield decrease per unit rise in SAL. Moreover, Magan et al. (53) discovered that tomato overall and marketable yield dropped remarkably with rising SAL. At present, although there are different opinions concerning how SAL stress affects tomato yield, many researchers believe that fruit size is remarkably dropped with rising SAL (54,55).
Most yield variables were significantly lower for 3FT than for NT treatment. This could be explained by lower values of fruit number and marketable fruit number for 3FT treatment, because only three fruits remained for each cluster. Therefore, the fruit set and marketable fruit set were lower than in the NT treatment. However, yield under SF and marketable yield under MF and SF (combined treatment) were slightly higher than CN condition. This result indicated that the FT method (3FT) in this study failed to achieve the threshold for marketable yield. We may attempt in subsequent experiments to retain four or more fruits to maximize marketable yield. This result also confirms the different effects of SAL stress and FT on tomato yield and marketable yield. Salinity stress reduces yield and marketable yield, whereas reasonable FT can increase yield and marketable yield by increasing weight per fruit. The values of yield variables also exhibited slight differences under the same SAL condition (e.g., MN and MF, and SN and SF). This result further elucidates that the loss of tomato yield due to SAL stress can be compensated by reducing fruit load and removing substandard fruit in advance, thus achieving no reduction in yield or even an increase in marketable yield. In other words, proper FT can alleviate the effect of SAL stress on tomato yield.
The discrepancies in CWC before and after applying SAL were caused by SAL stress rather than FT. Zhang et al. (22) showed that SAL would reduce root water absorption via the osmosis effect and afterward induces water stress. Ahanger et al. (18) indicated as well that plant water uptake is reduced with increased SAL. SAL stress disrupts the osmosis balance and causes dropping water absorption and closing of stoma aperture, which leads to the suppression of evapotranspiration (56).
The daily water consumption was 0.92 × 10 −3 mm plant −1 before applying SAL. The daily water consumption after applying SAL ranged from 0.29 × 10 −3 to 1.62 × 10 −3 mm plant −1 with a mean of 0.86 × 10 −3 mm plant −1 . During the whole growth period, daily water consumption ranged from 0.55 × 10 −3 to 1.31 × 10 −3 mm plant −1 with a mean of 0.87 × 10 −3 mm plant −1 . Wang et al. (57) showed that daily water consumption ranged from 71.6 × 10 −3 to 91.9 × 10 −3 mm plant −1 under greenhouse soil cultivation. The daily water consumption of tomatoes under soilless cultivation conditions was significantly lower than that found for soil cultivation environments, although comparability of the results might be affected by such differences as cultivation environment and cultivar.
The WUE variables showed an initial rise and a subsequent reduction with increasing SAL ( Table 4). The SAL interval used in this study was not conducive to evaluating the SAL threshold due to its large range; however, it could be approximately determined that the SAL thresholds of WUE variables were nearly 3.0 dS m −1 under soilless conditions using Vegetable Life A nutrient solution.
It is a common practice in commercial tomato production to add salt to nutrient solutions to improve fruit quality. Many experiments have indicated that SAL stress can increase tomato quality fruit (19). The 3FT treatment showed higher fruit quality. It may be because tomato plant growth is highly dependent on source, sink strength, and their equilibrium (58).
Leaves are both source and sink organs. However, fruit is the strongest sink for assimilate. Therefore, FT treatment will not decrease assimilation substance generation rate, whereas it will decrease the competitive activities for assimilate among fruits (59). The highest values of fruit quality variables were for MF and SF treatments because both SAL stress and FT improved tomato fruit quality, and their dual effects could greatly improve fruit quality.
Fruit cracking is affected by the cultivar, fruit size, fruit firmness, fruit shape, fruit growing, fruit cuticle and saccharinity, irrigation water quality, and environmental conditions (35). Among these factors, firmness has a strong influence on cracking (60). The fruit firmness and fruit cracking ratio decreased with increasing SAL in this study, likely because tomato skin becomes tougher and more resistant to tearing under saline conditions (60). Krauss et al. (61) showed that fruit firmness decreased with increasing SAL. However, most studies have reported increased (62) or unchanged (63) fruit firmness under saline conditions. The differences may result from the use of different cultivars and periods of fruit picking. In addition, lower fruit firmness mitigates fruit cracking due to enhanced cross-linking of pectin in the pericarp cells (64).
The results of path and logistic analyses indicated that the SAL stress tended to inhibit and delay plant growth. Zhang et al. (22) showed that SAL inhibits tomato plant growth by negatively affecting tomato roots, shoots, leaves, yield, and WUE.
The t 1 , t 2 , and t m values under three-fruit treatment were earlier than for the no-thinning treatment ( Table 7), by 17, 19, and 18 days, respectively. The v m under the three-fruit treatment was also lower than for the no-thinning treatment.
Thus, FT tended to advance and shorten the period of plant growth, possibly because the lower fruit load decreased the competition for assimilating (e.g., water, nutrients, and carbon) between the fruits and plant growth (65).

CONCLUSIONS
The experimental results showed that tomato growth variables were significantly decreased with increasing SAL, but FT treatment had less effect on tomato growth.
Salinity stress was confirmed to improve fruit quality and reduce fruit cracking, but reduced marketable yield by reducing fruit size. FT was confirmed to improve fruit quality and marketable yield by removing unhealthy fruit and increasing remaining fruit size, but increasing fruit cracking. Combined treatment of SAL stress and FT avoided the problem of subsize fruit under SAL stress treatment, reduced the fruit cracking found with FT treatment, and greatly improved marketable yield and fruit quality.
Salinity stress suppressed water absorption of tomato, but FT had no effect on water consumption. The daily water consumption before SAL treatment, after treatment, and for the whole growth period was 0.92 × 10 −3 , 0.86 × 10 −3 , and 0.87 × 10 −3 mm plant −1 , respectively. The daily water consumption of tomatoes was significantly lower for soilless than for soil cultivation. The SAL thresholds of WUE variables were determined to be approximately 3.0 dS m −1 under soilless conditions using Vegetable Life A nutrient solution.
Path analysis revealed that biomass and water consumption were vital indices positively influencing yield. Based on the logistic equation fitting, the SAL stress tended to inhibit and delay plant growth. The rapid growth period (t 1 ), end of the rapid growth period (t 2 ), and most rapid growth time (t m ) and SAL stress showed that increasing SAL by 1 dS m −1 resulted in 3.20, 6.39, and 5.06 days of delay, respectively. FT tended to advance and shorten the period of plant growth.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/supplementary materials, further inquiries can be directed to the corresponding author/s.

AUTHOR CONTRIBUTIONS
PZ, SJ, YD, and MS: conceptualization. PZ, SJ, ZZ, and YD: writing-original draft. PZ, YD, and MS: writing, reviewing, and editing. All authors contributed to the article and approved the submitted version.