- Research Institute of Xinjiang Uygur Autonomous Region Academy of Agricultural Sciences/Northwest Center of National Saline-alkali Tolerant Rice Technology Innovation Center, Urumqi, China
Introduction: This study aimed to determine the optimal nitrogen (N) application rate for mitigating low-temperature (LT) stress damage at different growth stages and to elucidate the underlying mechanisms in rice cultivation in Xinjiang.
Methods: A split-plot field experiment was conducted in 2024 using the rice cultivar ‘Xinjing 8’. The main plots included three temperature environments: control (CK), lowtemperature stress at the tillering (TLS), and low-temperature stress at the jointing stage (JLS), perform a 3-day LT treatment at 16°C. Subplots were assigned five N application rates: 295.5, 267.9, 240.3, 212.7, and 0 kg·ka-1. The effects of N rate on rice panicle traits, yield, and economic benefits under these growing environments were evaluated.
Results: The results indicated that under lowtemperature stress, rice exhibited significant reductions in panicle weight, plant height, number of secondary branches, effective panicles, and grain yield per panicle, particularly during the jointing stage. Nitrogen supplementation under low-temperature stress significantly improved rice tillering and panicle traits, simultaneously elevating grain yield and economic benefits. Specifically, at optimal nitrogen application rate of 267.9 kg ha-1, the crop yields under TLS and JLS treatments were 4424.40 and 2350.82 kg ha-1. At the same N application rate, a positive net income was achieved only in the CK treatment.
Discussion: This study identifies a nitrogen application rate of 267.9 kg ha-1 as an effective measure to mitigate low-temperature stress-induced limitations on rice panicle development, resulting in enhanced yield and economic gains. Consequently, it provides a practical solution for advancing sustainable rice cultivation in coldprone regions such as Xinjiang.
1 Introduction
Rice (Oryza sativa L.) is one of the most important global food crops, serving as a staple food for nearly half of the world’s population. However, with the continuous growth of the global population, it is estimated that by 2050, rice production must double compared to 2005 to meet the increasing food demand (Verma et al., 2021). Recently, global warming has markedly shifted the temperature variability and extremes around the world, and the occurrence of periodic low-temperatures (LT) may increase (Shen et al., 2018). China has the world’s second largest rice cultivation area and the world’s highest rice production, accounting for 29.9% of global production. Abnormal chilling events increase risks and losses to rice production, decrease agricultural income stability, and thereby threaten food supply stability. Therefore, exploring strategies to mitigate the impact of low-temperature on yield loss and improving the efficiency of low-temperature prevention and control for rice in China is of critical research significance for ensuring China’s food security (Shimono et al., 2007).
Rice originated in tropical regions (Yuan et al., 2021), and exhibits a preference for warm environments, with an optimal growth temperature range is 25-32°C, while being highly sensitive to low-temperature stress (Ma et al., 2024). Low-temperature affect rice starting from seed germination until physiological maturity. If rice seedlings encounter low-temperatures after spring sowing, they may exhibit slow growth, yellowing, poor development, wilting, and even death, seriously affecting rice growth, development, and final yield (Jia et al., 2022b). During the vegetative growth period, Low-temperatures dropping below 15°C can inhibit the differentiation of tiller primordia, reducing the number of effective tillers by 40-60% (Guo et al., 2022). Low-temperature usually exacerbates the balance between the source of energy and the metabolic sink, which inhibits the rate of photosynthesis and results in reduced green leaf area and biomass accumulation (Shi et al., 2022a). Low-temperatures during growth stages of booting result in pollen sterility, delayed and inadequate heading, and restricted anther dehiscence (Ali et al., 2021; Gao, et al., 2024; Shi et al., 2022b), which can lead to a 30% to 40% reduction in rice yield, and in extreme cases, even total crop failure. Low-temperatures at the rice grain filling stages delays the rice maturity date and cause frost disasters, which ultimately leads to insufficient grain filling, rice grain weight decrease and yield reduction, along with significant alterations in amylose and protein content in rice grains (Guo et al., 2024). Researchers have been working on screening rice varieties with enhanced cold tolerance (Lee et al., 2025; Nie et al., 2020). However, large-scale adoption of such varieties remains challenging. Most existing studies on low-temperature stress in rice have focus on the Northeast China (Ma et al., 2023; Shi et al., 2023). Xinjiang, characterized by a typical arid and semi-arid continental climate, achieves a rice yield per unit area that is 35.18% higher than the national average (National Statistical Office, 2024). However, research on the effects of low-temperatures on rice in this region remains relatively limited. In production practice, nitrogen fertilizer is usually applied when rice leaves show yellow due to low-temperatures, aiming to promote tillering and improve yield (Wang et al., 2016). Yet, excessive nitrogen application after temperature returns to normal can lead to redundant tillering, prolong the vegetative growth phase, and delay maturity (Peng et al., 2018). It is uncertain whether nitrogen application under low-temperature conditions is effective in supporting the recovery and normal growth of rice.
This study focused on the critical growth stages of tillering and jointing for rice. Based on the experimental data for Nitrogen application treatments, we compared the response of yield traits to chilling treatment at different growth stages for the selected rice variety. This study was conducted with the following objectives: (1) To evaluate the effects of low-temperature (LT) stress at critical vegetative stages on rice panicle morphology, yield components, and economic returns. (2) To identify the optimal nitrogen (N) application rate for alleviating LT stress and maximizing productivity and profitability in Xinjiang rice production systems. (3) To provide a scientific basis for developing defensive strategies against LT stress and rational N management practices in cold-affected seasons.
2 Materials and methods
2.1 Overview of the experimental site
This experiment was conducted from April to October 2024 at the Chabuchar Rice Experiment Station of Xinjiang Academy of Agricultural Sciences (43°50′ N, 81°7′ E). The region is characterized by an average annual sunshine duration of 2810.4-2846.0 hours, an average annual precipitation of 222 mm, and a frost-free period ranging from 169 to 177 days. It experiences cold winters and hot summers, with a mean annual temperature of 7.9 °C. The soil at the experimental site is clayey, with the following physicochemical properties (0–30 cm depth): pH 8.0, electrical conductivity 144 mS·cm-¹, organic matter 20.11 g·kg-¹, total nitrogen 1.65 g·kg-¹, available nitrogen 113.40 mg·kg-¹, available phosphorus 36.25 mg·kg-¹, and available potassium 124.20 mg·kg-¹. Temperature variations during the experimental period (April–October 2024) are presented in Figure 1.
Figure 1. Shows the air temperature during the rice growing period (April to October) in 2024. “MM” is the month and “DD” is the day. The same as below.
2.2 Test materials
The experimental material employed in this study was the rice cultivar ‘Xinjing 8’, which was developed locally by the Xinjiang Academy of Agricultural Sciences. Nitrogen, phosphorus, and potassium were supplied in the form of urea (N ≥ 46%; Xinjiang, China), double superphosphate (46% P2O5; Yunnan, China), and potassium sulfate (K2O ≥ 50%; Qinghai, China), respectively.
2.3 Experimental design
A split-plot design with two factors was employed. The main-plot factor consisted of three temperature environments: low-temperature treatment at the tillering stage (TLS), low-temperature treatment at the jointing stage (JLS), and conventional field conditions as the control (CK). The sub-plot factor involved five nitrogen (N) application rates based on side-deep fertilization technology. Previous studies indicated that an annual application of 295.5 kg·ha-¹ (N) could achieve high yields under local conditions. Based on this rate, the following five N levels were set: 100% (N295.5), 80% (N267.9), 60% (N240.3), 40% (N212.7), and a zero-N control (N0), resulting in a total of 15 treatment combinations.
Cold tolerance test method: At the time of transplanting, 9 plants of each treatment to be tested were transplanted into individual pots (10 cm×10 cm) respectively, and labeled. After transplanting, they were placed in the field under the same water and fertilizer management as the field control. At the tillering or jointing stages, they were transferred to an artificial climate chamber for low-temperature treatment for 3 days, with the constant temperature set at 16°C and the light duration at 16 hours. After 3 days, the pots were removed, and the plants were buried in the field until maturity. The field control plants were labeled but received no low-temperature treatment. Phosphorus (P) and potassium (K) fertilizers were applied uniformly across all treatments. Phosphorus was applied as a one-time basal dose at 138 kg·ha-¹ (P2O5). Potassium was split-applied, with 75 kg·ha-¹ (K2O) as basal fertilizer and 37.5 kg·ha-¹ (K2O) as panicle fertilizer. Detailed nitrogen application schedules for each growth stage are provided in Table 1.
2.4 Field management
Bowl-shaped blanket seedling trays were used for seedling raising. During transplantation, side-deep fertilization was simultaneously applied using a Yanmar YR60D (2ZGQ-60D; Jiangsu, China) rice transplanter. With this technology, both the base and tillering fertilizers were placed 3 cm to the side and 5 cm deep relative to the rice roots. For transplanting cultivation, the seedlings were raised in the nursery for 28 days before being transplanted into the field with a row spacing of 30 cm and a distance of 14 cm between hills. A large-plot experimental design was employed to accommodate mechanized transplanting operations. Each treatment was assigned to a 90 m2 main plot, which was subdivided into three 30 m2 (10 m × 3 m) sub-plots. Ridges measuring 30 cm in height and 50 cm in width were constructed between sub-plots and covered with plastic film to enable separate irrigation and drainage, thereby preventing water and fertilizer interaction among adjacent experimental units. All other field management practices followed local commercial production standards.
2.5 Measured items and methods
2.5.1 Growth stages recording
The dates of each growth stage of rice (sowing period, transplanting period, heading period, full heading period, and maturity period) were accurately recorded based on the standard for each growth period of rice (Agricultural Industry Standard of the People's Republic of China, 2021).
2.5.2 SPAD value
Starting from 30 days after transplanting, the SPAD (Soil Plant Analysis Development) values of rice leaves were measured at 10-day intervals. On each sampling date, between 11:00 and 13:30 under clear and windless conditions, nine uniformly growing hills per plot were selected. The chlorophyll content of the uppermost fully expanded leaf (flag leaf) was measured using a hand-held SPAD-502 Plus chlorophyll meter (produced by Konica Minolta, Japan), and the average value per plot was calculated.
2.5.3 Tiller dynamics
Seven days after transplanting, nine representative hills in each plot were randomly selected and tagged. The number of tillers per hill was counted every 7 days until the full heading stage. Data on basic seedlings, maximum seedling number, and effective panicle number were recorded. The tillering rate and panicle-forming rate (effective tiller percentage) were calculated as follows:
2.5.4 Yield component analysis and grain quality measurement
Before the rice harvest, the number of effective panicles in 9 rice clusters was investigated in each plot. After harvesting and threshing, the samples were air-dried to a constant weight for yield component analysis. The following parameters were determined: panicle length, weight per panicle, number of primary branches, number of secondary branches, number of effective panicles, number of grains per panicle, number of filled grains, number of empty and shriveled grains, and 1000-grain weight, etc. Theoretical yield was calculated based on these components. Grain characteristics were measured using SC-G (Wanshen, Hangzhou, China) seed detector.
2.5.5 Economic analysis
The total economic input (expressed in CNY·ha-¹) was calculated by aggregating the costs associated with all cultivation factors. These included expenditures on seeds, nursery trays, substrate, agricultural film, irrigation, personnel labor, fertilizers and agrochemicals, mechanical operations, and other incidentals. All cost data were acquired through local market surveys. The average market price of rice was taken as 3.00 CNY·kg-¹, and the currency conversion rate used was 1.00 USD = 7.11 CNY. The output income, net income, input-output ratio were calculated as followed:
2.6 Data processing
Data processing was conducted using Microsoft Excel 2010, graphs were made with Origin 2022. The data were subjected to variance analysis using SPSS 26.0 to determine the treatment effects on the agronomic traits. The differences between means were determined using the least significant difference test at P < 0.05. For comparisons under a specific temperature treatment, Duncan’s multiple range test was applied to identify significant differences in nitrogen uptake.
3 Results
3.1 Impact of low-temperature on agronomic traits of rice under different nitrogen application rates
3.1.1 Phenological development
Under consistent environmental conditions, the first heading date, full heading date and maturing date of rice are delayed with the increase of nitrogen application, and the total growth period is extended (Table 2). The total growth period of rice under N295.5 treatment was the longest, lasting 173 days (CK), 175 days (TLS) and 178 days (JLS) respectively, and it was prolonged by 1–13 days compared to other nitrogen application treatments. Under the same nitrogen application treatment, the total growth period of JLS was extended by 1–3 days compared to TLS. Specifically, the N267.9 treatment showed only a 1-day extension, whereas all other treatments exhibited a 3-day extension.
Table 2. Effects of low-temperature on phenological period of rice with different nitrogen application rates.
3.1.2 SPAD
Under the same environmental conditions, the SPAD values of rice leaves under different nitrogen application rates exhibited an initial increase followed by a decreasing trend throughout the growth cycle (Figure 2). The SPAD values of rice leaf blades reached the highest on July 15th, with the maximum under the N267.9 treatment being 43.99 (CK), 45.58 (TLS), and 43.77 (JLS), respectively. The overall trend indications were N267.9 > N295.5 > N240.3 > N212.7 > N0. Under the same nitrogen application treatment, the decline in leaf SPAD values TLS and JLS was less than that of CK, with the smallest decline occurring JLS.
Figure 2. Effect of low-temperature on chlorophyll SPAD value in rice with different nitrogen application rates. CK, conventional field conditions as the control; TLS, low-temperature treatment at the tillering stage; JLS, low-temperature treatment at the jointing stage.
3.1.3 Tillering rate and spike rate
Under consistent nitrogen application rates, the tillering rate of rice followed a similar trend in both CK and JLS, with both being significantly higher than that in TLS. Notably, the tillering rate in TLS was on average 40.19% lower than that in JLS and the CK (Figure 3A). In contrast, the spike rate of rice in TLS increased by an average of 89.89% and 14.14%, respectively, compared to JLS and CK (Figure 3B).
Figure 3. Effects of low-temperature on tillering rate and spike rate of rice with different nitrogen application rates. Different letters in lowercase indicate significant difference of the data in the same temperature stress treatment at P < 0.05. Data are the means ± standard deviation (means ± SD). (A) Tillering rate; (B) Spike rate.
Under the same environmental conditions, the tillering rate of rice decreased with reduced nitrogen application. The highest tillering rate was observed under the N295.5 treatment. However, the difference between N267.9 and N295.5 treatment was not significant, though both differed significantly from other nitrogen treatments (P<0.05). Meanwhile, the spike rate exhibited an initial increase followed by a decrease as nitrogen application was reduced, with the maximum spike rate recorded under the N212.7 treatment.
3.1.4 Plant height and stem diameter
Under consistent nitrogen application rates, rice plant height followed the order: CK > JLS > TLS (Figure 4A). In contrast, the stem diameter of rice was exactly the opposite, with JLS being the thickest and CK the thinnest (Figure 4B). Under the same environmental conditions, as the nitrogen application rate decreases, the plant height of rice plants decreases accordingly, while the stem thickness of rice shows a trend of increasing first and then decreasing. Among them, the stem diameter reached the peak under the treatment of N267.9, which was 4.26 mm (CK), 4.40 mm (TLS), and 4.53 mm (JLS), respectively. Notably, there was no significant difference between the N267.9 treatment and the N295.5 and N240.3 treatments, while significant differences were detected compared to the N212.7 and N0 treatments.
Figure 4. Effects of low-temperature on plant height and stem diameter of rice with different nitrogen application rates. Different letters in lowercase indicate significant difference of the data in the same temperature stress treatment at P <0.05. Data are the means ± standard deviation (means ± SD). (A) Plant height; (B) Stem diameter.
3.2 Effects of low-temperature on panicle morphology of rice with different nitrogen application rates
3.2.1 Panicle length and single spike weight
When the nitrogen application rate was consistent, relative to CK, the panicle lengths of TLS and JLS exhibited average reductions of 17.93% and 15.49%, respectively (Figure 5A). A similar trend was observed for the single spike weight, which average decreased by 67.74% and 51.54%, respectively. These findings indicate that low-temperature adversely affects the morphological development of rice panicles, especially with a significant reduction in the weight per panicle (Figure 5B).
Figure 5. Effects of low-temperature on panicle length and single spike weight of rice with different nitrogen application rates. Different letters in lowercase indicate significant difference of the data in the same temperature stress treatment at P <0.05. Data are the means ± standard deviation (means ± SD). (A) Panicle length; (B) Single spike weight. (C) Rice panicle morphology.
Under the identical environmental conditions, both panicle length and single panicle weight of rice exhibited a trend of “increasing first and then decreasing” with the reduction of nitrogen application rate and were the largest under the N267.9 treatment. These findings suggest that under low-temperature conditions, an appropriate nitrogen application level is conducive to increasing panicle length and single panicle weight, whereas excessive nitrogen may inhibit their growth and development (Figure 5C).
3.2.2 Number of branches
When the nitrogen application rate was consistent, the number of primary branches and secondary branches in CK was significantly higher than that of TLS and JLS (Figure 6). Among them, there was no significant difference in the primary branches between TLS and JLS. The number of secondary branches in TLS was significantly higher than that in JLS (P < 0.05). Under the same environmental conditions, the number of primary branches and secondary branches of rice exhibited a trend of “increasing first and then decreasing” with a reduction in nitrogen application rate, with the largest values recorded in the N267.9 treatment. Among them, the number of secondary branches was 22.05 per panicle (CK), 7.91 per panicle (TLS), and 2.93 per panicle (JLS), respectively. This indicates that CK had 178.76% and 652.56% more secondary branches than TLS and JLS, respectively.
Figure 6. Effects of low-temperature on number of branches with different nitrogen application rates. Different letters in lowercase indicate significant difference of the data in the same temperature stress treatment at P <0.05. Data are the means ± standard deviation (means ± SD). (A) Primary branches; (B) secondary branches.
3.2.3 Grain morphology of rice
Under consistent nitrogen application rates, grain area, circumference, length, width, and diameter all followed the order: CK > JLS > TLS, with JLS showing a significant reduction in these parameters (Table 3). In contrast, the grain length-to-width ratio exhibited an opposite trend. Specifically, the length-to-width ratio of JLS increased significantly by 6.77% and 6.88% compared to CK and TLS, respectively. Additionally, grain roundness in JLS was significantly lower than that in CK and TLS.
Table 3. Effects of low-temperature on grain traits of rice with different nitrogen application rates.
Under identical environmental conditions, all grain dimensional traits (area, circumference, length, width, and diameter) were highest under the N0 treatment. The grain length-to-width ratio peaked under the N295.5 treatment in both CK and JLS, with values of 2.716 and 3.007, respectively. In TLS, however, the maximum length-to-width ratio was observed under the N267.9 treatment (2.726), which showed a significant difference only when compared to the N240.3 treatment (P < 0.05).
3.3 Impact of low-temperature on rice yield and economic benefits under different nitrogen application rates
3.3.1 Yield and yield components
Under consistent nitrogen application rates, grain yield, effective panicle number, total grain number per panicle, seed setting rate, and 1000-grain weight all followed the order: CK > JLS > TLS (Table 4). Among these traits, effective panicle number and grain yield showed significant differences among the three treatments. With the exception of the N0 treatment, no significant difference was observed in total grain number per panicle between TLS and JLS. In addition, no significant difference was detected in seed setting rate or 1000-grain weight between CK and TLS, although both were significantly higher than JLS.
Table 4. Effects of low-temperature on rice yield and yield components at different nitrogen application rates.
Under identical environmental conditions, grain yield, effective panicle number, and total grain number per panicle exhibited an initial increase followed by a decrease as nitrogen application was reduced, with maximum values recorded under the N267.9 treatment. Specifically, grain yields were 10,549.42 kg·ha-¹ (CK), 4,424.40 kg·ha-¹ (TLS), and 2,350.82 kg·ha-¹ (JLS), respectively. In contrast, the 1000-grain weight increased with decreasing nitrogen application, reaching its highest value under the N0 treatment, which differed significantly from all other nitrogen treatments (P < 0.05).
3.3.2 Economic benefits
Under consistent nitrogen application rates, net income was positive only in the CK treatment among all fertilized groups (excluding N0), with an input-output ratio exceeding 1 (Table 5). In contrast, both TLS and JLS resulted in negative net income and input-output ratios below 1, with JLS experiencing more severe economic losses. Under identical environmental conditions, the N0 treatment (no nitrogen application) also yielded negative net income, indicating that nitrogen fertilization generally enhances the economic returns of rice production.
Table 5. Effects of low-temperature on the economic benefits of rice with different nitrogen application rates.
When environmental conditions were maintained constant, total input cost increased with rising nitrogen application rates. However, output income, net income, and input-output ratio all showed an initial increase followed by a decrease as nitrogen application was reduced, with peak values observed under the N267.9 treatment.
4 Discussion
4.1 The impact of low-temperature on yield formation
Chlorophyll, the most critical component for photosynthesis in plants, is responsible for absorbing and converting light energy (Fotovat et al., 2007). Leaf chlorophyll content not only serves as an indicator of plant nutritional status but also reflects the crop’s tolerance to adverse stress conditions (Neufeld et al., 2006). Zhao et al. reported that low-temperature stress disrupts chlorophyll biosynthesis and chloroplast development in rice, leading to leaf yellowing or chlorosis and a significant decrease in chlorophyll content (Zhao et al., 2020). However, the present study observed higher SPAD values in later growth stages under low-temperature stress compared to field controls. This apparent contradiction may be attributed to the suspension of rice growth and development under low-temperature stress, which requires additional time for metabolic recovery and chilling damage mitigation. The extended growth duration likely allowed for compensatory chlorophyll synthesis and chloroplast repair (Lu et al., 2021), resulting in delayed phenological development but ultimately higher chlorophyll content in later stages. This aligns with findings from other studies (Li et al., 2003; Guo et al., 2024), who also reported that low-temperatures can extend the rice growth period.
The present study found that low-temperature at the jointing stage (JLS) significantly delayed both vegetative and reproductive development, extending the total growth duration. The application of appropriate nitrogen rates (N267.9) enhanced recovery capacity after cold damage and reduced the extension of growth days. This suggests that nitrogen nutrition plays a crucial role in modulating plant responses to low-temperature stress, possibly through its involvement in the synthesis of protective compounds, membrane components, and stress-related proteins (Ye et al., 2022). The tillering and jointing stages represent critical phases for rice vegetative growth. Most studies have shown that low-temperatures at the tillering stage reduce both tillering rate and panicle formation rate (Shimono and Okada, 2013). However, some research suggests that although low-temperature decreases tiller numbers, it may improve individual development and increase panicle formation rate (Zhang et al., 2024). The results of this study generally support the latter view. Specifically, low-temperature at the tillering stage (TLS) significantly reduced rice tillering rate and plant height but increased panicle formation rate, whereas low-temperature at the jointing stage (JLS) significantly reduced panicle formation rate. This pattern suggests that plants under stress may employ a resource optimization strategy, sacrificing tiller quantity to maintain the quality of surviving tillers. The increased panicle formation rate under TLS could result from reduced competition among fewer tillers for limited resources.
The panicle represents a critical site for rice yield formation, and its morphological development constitutes a key research focus in ideal plant-type breeding (Song et al., 2022). Panicle morphology is influenced by multiple factors, including grain number per panicle, number and length of branches, grain distribution density, and panicle length (Li et al., 2021). Most studies indicate that low-temperature stress significantly reduces the number of primary branches, secondary branches, and total branches in rice (Jia et al., 2022a). Consistent with these findings, our study observed significant reductions in panicle length, single panicle weight, and numbers of primary and secondary branches in both TLS and JLS treatments compared to CK. The more pronounced reductions in single panicle weight and secondary branch number in JLS than TLS suggest that the jointing stage represents a particularly sensitive period for panicle development. This sensitivity may be attributed to the active differentiation of panicle components during this stage, making them vulnerable to temperature stress through disruption of hormonal balance and carbohydrate partitioning (Sarma et al., 2023).
Numerous studies have demonstrated that low-temperatures reduce rice yield by decreasing effective panicle number, grains per panicle, seed setting rate, and 1000-grain weight (Huang et al., 2013; Shi et al., 2025). In our experiment, while no significant differences were observed in seed setting rate and 1000-grain weight between TLS and CK (both being significantly higher than JLS), the effective panicle number of JLS was significantly lower than that of TLS. The differential impact on yield components across developmental stages reflects the stage-specific nature of low-temperature damage. The jointing stage appears particularly crucial for determining panicle number and grain development, possibly due to its role in panicle initiation and differentiation processes (Ren et al., 2017). Consequently, the average yields of TLS and JLS decreased by 61.60% and 81.99%, respectively, compared to CK. Among treatments with the same nitrogen application rate, only CK demonstrated positive net income, while both TLS and JLS resulted in economic losses. These results indicate that low-temperature exposure during the jointing stage leads to more severe yield reduction and economic loss compared to stress during earlier growth phases. This has significant implications for rice production in temperate regions and highlights the need for targeted management strategies during sensitive growth stages. This has significant implications for rice production in subtropical regions and highlights the need for targeted management strategies during sensitive growth stages.
4.2 Influence of nitrogen application rate on yield formation
In this experiment, under consistent environmental conditions, the first heading date, full heading date, and maturation date of rice were progressively delayed with increasing nitrogen application rates, resulting in a corresponding extension of the total growth period. This observation aligns with previous findings that excessive nitrogen application leads to delayed crop maturity (Zhao et al., 2022). Regarding chlorophyll content, the highest nitrogen treatment (N295.5) resulted in a 2.79% lower SPAD value compared to the N267.9 treatment. This indicates that moderate nitrogen application promotes chlorophyll synthesis in rice leaves (Zhou et al., 2023), whereas excessive nitrogen may disrupt nitrogen metabolism homeostasis, potentially inducing feedback inhibition or reactive oxygen species accumulation that impairs chloroplast structure and function (Kim and Kim, 2024). This nonlinear response highlights the existence of an optimal range for nitrogen’s regulatory effect on the photosynthetic apparatus, beyond which inhibitory effects prevail (Song et al., 2023).
Nitrogen is a crucial nutrient for rice growth and development, playing a central role in tiller formation (He et al., 2024). In this experiment, a moderate increase in nitrogen fertilizer promoted tiller initiation and improved tiller quality, thereby increasing the effective panicle number and grains per panicle (Peng et al., 2018). In contrast, excessive nitrogen application suppressed the growth of tillers at the middle and lower leaf positions, reducing the number of effective tillers (Zhou et al., 2022). This suppression may be attributed to intensified light competition within the canopy and reduced photo assimilate allocation to tillers under high nitrogen conditions. Furthermore, this study found that appropriate nitrogen application under low-temperature conditions significantly promoted rice tillering, increased effective panicle number, and simultaneously enhanced both yield and nitrogen use efficiency, consistent with previous reports (Zhou et al., 2018). These results suggest that optimal nitrogen availability may alleviate the inhibitory effects of low-temperature on tiller bud germination and growth by improving the coordination of carbon and nitrogen metabolism (Liu et al., 2019; Wang et al., 2023).
Multiple studies have demonstrated that the number of primary and secondary branches in rice increases with rising nitrogen application rates (Ding et al., 2014; Zhou et al., 2017), a trend also confirmed in this study. Specifically, panicle length, single panicle weight, and branch number increased with nitrogen application up to the N267.9 treatment but declined under the N295.5 treatment, generally following the order: N267.9 > N295.5 > N240.3 > N212.7 > N0. This threshold response indicates that optimal nitrogen application simultaneously promotes the growth of both vegetative and reproductive organs, effectively mitigating the adverse effects of low-temperature, and enhancing cold tolerance and yield potential (Xiong et al., 2025). Conversely, excessive nitrogen may lead to insufficient carbohydrate supply to the panicles or disrupt hormonal regulation, ultimately limiting the full expression of panicle traits.
Over the past six decades, rice yields and nitrogen application rates in China have increased almost simultaneously (Zhang and Liu, 2022). However, a universally applicable optimal nitrogen rate for achieving maximum yield remains undefined. A comprehensive analysis of 514 field experiments across seven rice-producing regions in Hubei Province demonstrated that yields in nitrogen-fertilized plots were significantly higher than in non-fertilized plots, with medium nitrogen treatment yielding the best results (Wang et al., 2012). Numerous studies have reported that excessive nitrogen application can lead to over-luxuriant growth and delayed maturation, impairing grain filling, reducing the seed-setting rate and 1000-grain weight, and ultimately causing yield loss (Chen et al., 2022; Huang et al., 2020). Under identical environmental conditions in this study, effective panicle number, total grains per panicle, yield, net income, and input-output ratio all exhibited an initial increase followed by a decrease as nitrogen application was reduced, peaking under the N267.9 treatment. Specifically, yields were 10,549.42 kg·ha-1 (CK), 4,424.40 kg·ha-1 (TLS), and 2,350.82 kg·ha-1 (JLS), respectively. Under CK conditions, output income, net income, and input-output ratio reached 31,648.26 CNY·ha-1, 14,018.16 CNY·ha-1, and 1.80, respectively. Notably, net income under the N267.9 treatment increased by 29.60%, 31.11%, and 98.66% compared to the N295.5, N240.3, and N212.7 treatments, respectively. These results clearly demonstrate that excessive nitrogen application under low-temperature conditions at both the tillering (TLS) and jointing (JLS) stages is detrimental to achieving high yield and income in rice production, while also resulting in significant nitrogen wastage (Ai et al., 2024), a conclusion consistent with observations under normal temperature conditions.
Rice yield is primarily determined by effective panicle number, with secondary contributions from grains per panicle, seed-setting rate, and 1000-grain weight (Bahmutsky et al., 2024). This study also found that no significant differences in seed-setting rate or 1000-grain weight across different nitrogen application rates, whereas nitrogen application significantly influenced effective panicle number and total grains per panicle. Therefore, the synergistic optimization of yield components through appropriate nitrogen fertilization is key to stabilizing rice production in low-temperature environments (Liu et al., 2019).
5 Conclusions
The results demonstrated that chilling treatment at the rice tillering and jointing stages affected the maturity date, grain yield, yield components and economic benefits, with more pronounced effects during the jointing stage. Nitrogen supplementation under low-temperature stress significantly improved rice tillering and panicle traits, simultaneously elevating grain yield and economic benefits. The single panicle weight, number of first branches and secondary branches, effective panicle number, total grain number per panicle, net income, and input-output ratio of rice all reached peak values at a nitrogen application rate of 267.90 kg·ha-1. Further increases in nitrogen application beyond this level resulted in a decline across all measured parameters. At this rate, yields under TLS and JLS conditions were 4424.40 kg·ha-1 and 2350.82 kg·ha-1, representing increases of 16.59%–102.04% and 11.09%–273.95%, respectively, compared to other nitrogen treatments. Among all nitrogen treatments, only the field control (CK) yielded positive net income under standard temperature conditions. The N267.9 treatment achieved a net income of 14,018.16 CNY·ha-1, exceeding the N295.5, N240.3, and N212.7 treatments by 29.60%, 31.11%, and 98.66%, respectively. In summary, a nitrogen application rate of 267.90 kg·ha-1 can effectively mitigate low-temperature damage in rice and serves as a recommended strategy for achieving synergistic improvements in both yield and economic return under chilling stress conditions. Future research should focus on elucidating the physiological and molecular mechanisms by which nitrogen regulates rice adaptation to low-temperature, providing a theoretical basis for climate-smart nitrogen management strategies.
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 authors.
Author contributions
XD: Investigation, Data curation, Formal analysis, Writing – original draft, Writing – review & editing. YZ: Project administration, Methodology, Writing – original draft, Writing – review & editing. TH: Validation, Software, Writing – original draft, Writing – review & editing. ZZ: Data curation, Supervision, Investigation, Writing – review & editing. DL: Formal analysis, Investigation, Writing – review & editing. MK: Validation, Investigation, Writing – review & editing. XW: Conceptualization, Methodology, Writing – original draft. FT: Resources, Investigation, Writing – review & editing. BA: Visualization, Investigation, Writing – review & editing. YQ: Data curation, Resources, Methodology, Writing – original draft. QY: Software, Methodology, Writing – original draft, Writing – original draft. YJ: Project administration, Methodology, Writing – review & editing. FW: Funding acquisition, Conceptualization, Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This work was supported by the Key research and development Project of Xinjiang Uygur Autonomous Region (2022B02001-5), the Project of Fund for Stable Support to Agricultural Sci-Tech Renovation (xjnkywdzc-2023001-19, xjnkywdzc-2025001-14), Xinjiang Uygur Autonomous Region’s “Tianshan Talent” Training Program-High-level Talents for Development of System Agriculture and Animal Husbandry Industry (2024SNGGGCC035), the Earmarked Fund for China Agriculture Research System (CARS-01-74), Youth Science Fund Project of the Xinjiang Uyghur Autonomous Region (2022D01B161).
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.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
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.
References
Agricultural Industry Standard of the People's Republic of China (2021). Technical specification for the identification of rice germplasm resources. (Beijing: China Agriculture Press).
Ai, H.-s., Fan, B., Zhou, Z.-Q., and Liu, J. (2024). The impact of nitrogen fertilizer application on air pollution: evidence from China. J. Environ. Manage. 370, 122880. doi: 10.1016/j.jenvman.2024.122880
Ali, I., Tang, L., Dai, J., Kang, M., Mahmood, A., Wang, W., et al. (2021). Responses of grain yield and yield related parameters to post-heading low-temperature stress in japonica rice. Plants. 10, 1425. doi: 10.3390/plants10071425
Bahmutsky, S., Grassauer, F., Arulnathan, V., and Pelletier, N. (2024). A review of life cycle impacts and costs of precision agriculture for cultivation of field crops. Sustainable Production and Consumption. 52, 347–362. doi: 10.1016/j.spc.2024.11.010
Chen, Y., Teng, Z., Yuan, Y., Yi, Z., Zheng, Q., Yu, H., et al. (2022). Excessive nitrogen in field-grown rice suppresses grain filling of inferior spikelets by reducing the accumulation of cytokinin and auxin. Field Crops Res. 283, 108542. doi: 10.1016/j.fcr.2022.108542
Ding, C., You, J., Chen, L., Wang, S., and Ding, Y. (2014). Nitrogen fertilizer increases spikelet number per panicle by enhancing cytokinin synthesis in rice. Plant Cell Rep. 33, 363–371. doi: 10.1007/s00299-013-1536-9
Fotovat, R., Valizadeh, M., and Toorchi, M. (2007). Association between water-use efficiency components and total chlorophyll content (SPAD) in wheat (Triticum aestivum L.) under well-watered and drought stress conditions. J. Food Agric. Environ. 5, 225. doi: 10.3168/jds.2006-793
Gao, Y., Li, Y., Huang, L., Zhao, J., Li, S., Lu, J., et al. (2024). Identification of the effects of low temperature on grain-setting rate of different types of late-season rice (Oryza sativa) during heading. Field Crops Res. 318, 109584. doi: 10.1016/j.fcr.2024.109584
Guo, E., Wang, L., Jiang, S., Xiang, H., Shi, Y., Chen, X., et al. (2022). Impacts of chilling at the tillering phases on rice growth and grain yield in Northeast China. J. Agron. Crop Sci. 208, 510–522. doi: 10.1111/jac.12596
Guo, E., Wang, L., Shi, Y., Jiang, S., Cheng, X., Chen, X., et al. (2024). Response of rice yield, grain quality to chilling at grain filling stages. J. Sci. Food Agric. 104, 9312–9323. doi: 10.1002/jsfa.13754
He, X., Zhu, H., Shi, A., and Wang, X. (2024). Optimizing nitrogen fertilizer management enhances rice yield, dry matter, and nitrogen use efficiency. Agronomy. 14, 919. doi: 10.3390/agronomy14050919
Huang, M., Jiang, L., Zou, Y., and Zhang, W. (2013). On-farm assessment of effect of low temperature at seedling stage on early-season rice quality. Field Crops Res. 141, 63–68. doi: 10.1016/j.fcr.2012.10.019
Huang, M., Lei, T., Cao, F., Chen, J., Shan, S., and Zou, Y. (2020). Grain yield responses to nitrogen rate in two elite double-cropped inbred rice cultivars released 41 years apart. Field Crops Res. 259, 107970. doi: 10.1016/j.fcr.2020.107970
Jia, M., Meng, X., Song, X., Zhang, D., Kou, L., Zhang, J., et al. (2022b). Chilling-induced phosphorylation of IPA1 by OsSAPK6 activates chilling tolerance responses in rice. Cell Discov. 8, 71. doi: 10.1038/s41421-022-00413-2
Jia, Y., Liu, H., Wang, H., Zou, D., Qu, Z., Wang, J., et al. (2022a). Effects of root characteristics on panicle formation in japonica rice under low temperature water stress at the reproductive stage. Field Crops Res. 277, 108395. doi: 10.1016/j.fcr.2021.108395
Kim, T.-H. and Kim, S.-M. (2024). Effects of SPAD value variations according to nitrogen application levels on rice yield and its components. Front. Plant Sci. 15. doi: 10.3389/fpls.2024.1437371
Lee, Y., Jang, S.-G., Jin, B., Jo, S.-M., Kang, J.-W., Lee, S.-M., et al. (2025). Integrating molecular markers and phenotypic analysis to assess cold tolerance in rice germplasm. Mol. Breed. 45, 23. doi: 10.1007/s11032-025-01543-7
Li, X., Qian, Q., Fu, Z., Wang, Y., Xiong, G., Zeng, D., et al. (2003). Control of tillering in rice. Nature. 422, 618. doi: 10.1038/nature01518
Li, G., Zhang, H., Li, J., Zhang, Z., and Li, Z. (2021). Genetic control of panicle architecture in rice. Crop J. 9, 590–597. doi: 10.1016/j.cj.2021.02.004
Liu, Z., Tao, L., Liu, T., Zhang, X., Wang, W., Song, J., et al. (2019). Nitrogen application after low-temperature exposure alleviates tiller decrease in rice. Environ. Exp. Bot. 158, 205–214. doi: 10.1016/j.envexpbot.2018.11.001
Lu, X., Song, S., Xiao, Y., Fan, F., Zhou, Y., Jia, G., et al. (2021). Circadian clock-coordinated response to chilling stress in rice. Environ. Exp. Bot. 185, 104398. doi: 10.1016/j.envexpbot.2021.104398
Ma, H., Shi, Y., Li, T., Wang, L., Zhang, X., Yang, X., et al. (2024). Chilling at grain filling stage reduced rice grain protein content: An experimental and modeling study. Environ. Exp. Bot. 226, 105891. doi: 10.1016/j.envexpbot.2024.105891
Ma, H., Zhu, X., Cao, T., Li, S., Guo, E., Shi, Y., et al. (2023). Responses of rice cultivars with different cold tolerance to chilling in booting and flowering stages: an experiment in Northeast China. J. Agron. Crop Sci. 209, 864–875. doi: 10.1111/jac.12663
National Statistical Office (2024). China Statistical Yearbook. Vol. 12 (Beijing: China Statistics Press).
Neufeld, H. S., Chappelka, A. H., Somers, G. L., Burkey, K. O., Davison, A. W., and Finkelstein, P. L. (2006). Visible foliar injury caused by ozone alters the relationship between SPAD meter readings and chlorophyll concentrations in cutleaf coneflower. Photosynthesis Res. 87, 281–286. doi: 10.1007/s11120-005-9008-x
Nie, L., Liu, H., Zhang, L., and Wang, W. (2020). Enhancement in rice seed germination via improved respiratory metabolism under chilling stress. Food Energy Secur. 9, e234. doi: 10.1002/fes3.234
Peng, X. L., Zhang, X. H., Liu, Z. L., Tao, L. Y., Liu, T. T., and Yu, C. L. (2018). Effect of nitrogen application after low temperature on rice yield and nitrogen use efficiency in cold area. J. Northeast Agric. Univ. 49, 10–20. doi: 10.3969/j.issn.1005-9369.2018.04.002
Ren, C., Fan, S., Liu, J., Zhu, L., Yang, Q., Liu, D., et al. (2017). QTL analysis for heading stage and panicle traits at different years in rice. J. Shenyang Agric. Univ. 48, 78–83. doi: 10.3969/j.issn.1000-1700.2017.01.011
Sarma, B., Kashtoh, H., Lama Tamang, T., Bhattacharyya, P. N., Mohanta, Y. K., and Baek, K. H. (2023). Abiotic stress in rice: visiting the physiological response and its tolerance mechanisms. Plants. 12, 3948. doi: 10.3390/plants12233948
Shen, X., Liu, B., and Lu, X. (2018). Weak cooling of cold extremes versus continued warming of hot extremes in China during the recent global surface warming hiatus. J. Geophysical Research: Atmospheres. 123, 4073–4087. doi: 10.1002/2017JD027819
Shi, Y., Guo, E., Cheng, X., Wang, L., Jiang, S., Yang, X., et al. (2022a). Effects of chilling at different growth stages on rice photosynthesis, plant growth, and yield. Environ. Exp. Bot. 203, 105045. doi: 10.1016/j.envexpbot.2022.105045
Shi, Y., Guo, E., Wang, L., Li, T., Jiang, S., Xiang, H., et al. (2022b). Effects of chilling at the booting and flowering stages on rice phenology and yield: A case study in Northeast China. J. Agron. Crop Sci. 208, 197–208. doi: 10.1111/jac.12549
Shi, Y., Ma, H., Li, T., Guo, E., Zhang, T., Zhang, X., et al. (2025). Innovative modeling on the effects of low-temperature stress on rice yields. J. Exp. Bot. 76, 1230–1243. doi: 10.1093/jxb/erae452
Shi, Y., Wang, L., Jiang, S., Guo, E., Li, T., Zhou, L., et al. (2023). An experimental study on the effects of intermittent chilling at different growth stages on rice yield in Northeast China. J. Agron. Crop Sci. 209, 317–329. doi: 10.1093/jxb/erae452
Shimono, H. and Okada, M. (2013). Plasticity of rice tiller production is related to genotypic variation in the biomass response to elevated atmospheric CO2 concentration and low temperatures during vegetative growth. Environ. Exp. Bot. 87, 227–234. doi: 10.1016/j.envexpbot.2012.11.008
Shimono, H., Okada, M., Kanda, E., and Arakawa, I. (2007). Low temperature-induced sterility in rice: evidence for the effects of temperature before panicle initiation. Field Crops Res. 101, 221–231. doi: 10.1016/j.fcr.2006.11.010
Song, Y., Dong, M., Hu, Y., Shi, L., Gu, J., Wang, Y., et al. (2023). Effects of nitrogen fertilization on nitrogen metabolism and grain protein in premium-tasty japonica rice. Cereal Res. Commun. 52, 1831–1842. doi: 10.1007/s42976-023-00480-5
Song, X., Meng, X., Guo, H., Cheng, Q., Jing, Y., Chen, M., et al. (2022). Targeting a gene regulatory element enhances rice grain yield by decoupling panicle number and size. Nat. Biotechnol. 40, 1403–1411. doi: 10.1038/s41587-022-01281-7
Verma, V., Vishal, B., Kohli, A., and Kumar, P. P. (2021). Systems-based rice improvement approaches for sustainable food and nutritional security. Plant Cell Rep. 40, 2021–2036. doi: 10.1007/s00299-021-02790-6
Wang, W., Lu, J., Ren, T., Li, X., Su, W., and Lu, M. (2012). Evaluating regional mean optimal nitrogen rates in combination with indigenous nitrogen supply for rice production. Field Crops Res. 137, 37–48. doi: 10.1016/j.fcr.2012.08.010
Wang, L. D., Zhao, H. W., Cai, H. L., Jia, Y., and Liu, H. L. (2016). Effects of nitrogen application on tillering and yield of japonica rice under cold water stress at tillering stage. Soil Fertilizer Sci. China. 6, 86–92. doi: 10.11838/sfsc.20160614
Wang, H., Zhong, L., Fu, X., Huang, S., Zhao, D., He, H., et al. (2023). Physiological analysis reveals the mechanism of accelerated growth recovery for rice seedlings by nitrogen application after low temperature stress. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1133592
Xiong, Z., Zheng, F., Wu, C., Tang, H., Xiong, D., Cui, K., et al. (2025). Nitrogen supply mitigates temperature stress effects on rice photosynthetic nitrogen use efficiency and water relations. Plants. 14, 961. doi: 10.3390/plants14060961
Ye, J. Y., Tian, W. H., and Jin, C. W. (2022). Nitrogen in plants: from nutrition to the modulation of abiotic stress adaptation. Stress Biol. 2, 4. doi: 10.1007/s44154-021-00030-1
Yuan, S., Linquist, B. A., Wilson, L. T., Cassman, K. G., Stuart, A. M., Pede, V., et al. (2021). Sustainable intensification for a larger global rice bowl. Nat. Commun. 12, 7163. doi: 10.1038/s41467-021-27424-z
Zhang, W. Y., Lei, Y. B., Zhang, Y. J., Guo, Z. H., Liu, N. S., Ma, W. D., et al. (2024). Research status on cold tolerance in rice and biotechnological breeding strategies for cold-tolerant early geng/japonica in Heilongjiang Province. J. Plant Genet. Resour. 25, 1428–1440. doi: 10.13430/j.cnki.jpgr.20240303001
Zhang, J. and Liu, Y. (2022). Agronomic improvements, not climate, underpin recent rice yield gains in changing environments. Agronomy. 12, 2071. doi: 10.3390/agronomy12092071
Zhao, Y., Han, Q., Ding, C., Huang, Y., Liao, J., Chen, T., et al. (2020). Effect of low temperature on chlorophyll biosynthesis and chloroplast biogenesis of rice seedlings during greening. Int. J. Mol. Sci. 21, 1390. doi: 10.3390/ijms21041390
Zhao, C., Liu, G., Chen, Y., Jiang, Y., Shi, Y., Zhao, L., et al. (2022). Excessive nitrogen application leads to lower rice yield and grain quality by inhibiting the grain filling of inferior grains. Agriculture. 12, 962. doi: 10.3390/agriculture12070962
Zhou, Y., Li, X., Cao, J., Li, Y., Huang, J., and Peng, S. (2018). High nitrogen input reduces yield loss from low temperature during the seedling stage in early-season rice. Field Crops Res. 228, 68–75. doi: 10.1016/j.fcr.2018.08.018
Zhou, W., Lv, T., Chen, Y., Hu, J., Zhang, Q., and Ren, W. (2017). Late nitrogen application enhances spikelet number in indica hybrid rice (Oryza sativa L.). Scientia Agricola. 74, 127–133. doi: 10.1590/1678-992X-2016-0094
Zhou, Z., Struik, P. C., Gu, J., van der Putten, P. E., Wang, Z., Yin, X., et al. (2023). Enhancing leaf photosynthesis from altered chlorophyll content requires optimal partitioning of nitrogen. Crop Environ. 2, 24–36. doi: 10.1016/j.crope.2023.02.001
Keywords: rice, low-temperature stress, nitrogen, panicle traits, grain yield, economic benefits
Citation: Du X, Zhang Y, Hou T, Zhao Z, Li D, Kang M, Wen X, Tang F, Abliz B, Qi Y, Yuan Q, Yuan J and Wang F (2025) Effect of nitrogen fertilizer regulation on rice panicle morphology, yield and economic benefits under low-temperature stress. Front. Plant Sci. 16:1720280. doi: 10.3389/fpls.2025.1720280
Received: 07 October 2025; Accepted: 24 November 2025; Revised: 18 November 2025;
Published: 19 December 2025.
Edited by:
Miguel Angel Prieto Lage, University of Vigo, SpainReviewed by:
Javier E. Alvarez, University of Vigo, SpainPauline Donn, University of Yaounde I, Cameroon
Ana Pérez Vázquez, University of Vigo, Spain
Copyright © 2025 Du, Zhang, Hou, Zhao, Li, Kang, Wen, Tang, Abliz, Qi, Yuan, Yuan and Wang. 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: Jie Yuan, eXVhbmppZTgwMTAyM0AxNjMuY29t; Fengbin Wang, eGpua3l3ZmJAMTYzLmNvbQ==
†These authors have contributed equally to this work and share first authorship
‡ORCID: Tianyu Hou, orcid.org/0009-0001-9551-9491
Zhiqiang Zhao, orcid.org/0009-0004-0013-6992
Dong Li, orcid.org/0009-0001-0319-7967
Mintai Kang, orcid.org/0009-0007-1575-666X
Xiaorong Wen, orcid.org/0009-0009-0597-7885
Fusen Tang, orcid.org/0009-0002-9675-4036
Buhaliqem Abliz, orcid.org/0009-0001-9748-4966
Yuhong Qi, orcid.org/0009-0003-2910-0603
Quan Yuan, orcid.org/0009-0004-5826-6511
Jie Yuan, orcid.org/0009-0009-4035-1626
Fengbin Wang, orcid.org/0000-0001-6906-5425
Tianyu Hou†‡