1. Introduction
The sugar industry, an integral part of the food and agriculture sector, faces critical global challenges, including increasing world population, hunger, and extreme climate change. These challenges can only be mitigated by increasing agricultural production. This goal can be achieved by expanding both the cultivation area and yield per unit area. Given the limited arable land, increasing yield per unit area is the most viable option [
1]. The sugar beet (
Beta vulgaris L.) is the second most vital and strategic crop for sugar production worldwide. Globally, sugar beets are cultivated on 4.3 million hectares of land, producing approximately 261 million tons with an average yield of 60.8 MG ha
-1. In Iran, the production stands at 5.0 million tons with an average yield of 55.6 MG ha
-1 [
2]. Due to its climatic diversity and latitude, Iran is among the six countries with sugar production infrastructure from sugarcane and sugar beets. Sugar beet cultivation in Iran has a history of more than 100 years. Over the past five years, the average cultivation area has been 109,000 ha, with an average yield of 56 MG ha
-1. Owing to the growing population, the need to supply sugar through sugar beet cultivation is increasing, leading to a cultivation area of 113,000 ha (92,000 and 21,000 ha of spring and autumn cultivation, respectively) in 2022. Notably, 50% of the sugar produced in Iran is obtained from sugar beet [
3].
In Iran, sugar beets account for approximately 7% of the total agricultural crop production and 49% of the total industrial production. West Azerbaijan Province ranks first (approximately 33% of production) [
3]. The high importance of sugar in Iranian households, the need to fulfill domestic demand through imports, and yield reduction due to abiotic stresses, such as drought, underscores the need to increase sugar beets' quantitative and qualitative yield.
The sugar beet germplasm in Iran exhibits favorable diversity in terms of quantitative and qualitative traits. Evaluating agronomic traits, specifically those with the highest impact on the final sugar beet yield, can enhance breeding programs. Understanding the relationships between traits, particularly those that are not easily measurable or have low heritability, is crucial for indirect selection. Identifying and understanding the relationships between traits with complex inheritance and low heritability and those with simple and higher heritability increases the genetic efficiency of breeding complex traits, such as sugar yield [
4]. Selecting correlated traits also impacts the primary trait, but the degree of dependence among traits is limited when the number of independent variables affecting the dependent trait increases. In such cases, correlations alone cannot explain the relationships between variables [
5]. The fundamental relationships between traits need to be expressed through path analysis, which provides additional information that is not generally observed in simple correlations [
6].
Path coefficient analysis is widely used in crop breeding to determine the relationship between yield (response or dependent variable) and its contributing components (predictor or independent variables). This analysis identified components with significant effects on yield for potential use as selection criteria, offering an acceptable explanation of the correlation between traits based on a cause-and-effect model. Many researchers have examined the cause-and-effect relationships between yield and other traits in various crops, including rice [
7], wheat [
8], barley [
9], rapeseed [
10], and sugar beet [
11,
12,
13,
14,
15]. The direct effects of leaf area [
14], root dry weight, total fresh weight, leaf fresh weight, crown fresh weight [
15], root diameter, and molasses sugar content [
12] on root yield, sugar content, white sugar content, alkalinity coefficient, and root potassium content have been previously reported [
15,
16].
Plant breeders generally select only a few traits, which makes it essential to understand their effects on other important characteristics. The current study was conducted to investigate the relationship between sugar yield and its components and to determine the most critical traits affecting sugar beet yield.
2. Materials and Methods
2.1. Experiment Layout
Thirteen different sugar beet monogerm cultivars, including 11 Iranian (SBSI035, SBSI037, SBSI048, Sina, Shokoofa, Paya, Motahar, Ekbatan, Aria, Sharif, and Arta), and two foreign cultivars (Perfecta and Palma) were evaluated during the 2018-2019 crop season. The study was conducted as a split-factorial experiment based on a randomized complete block design (RCBD) with four replications. The main plot consisted of two planting dates (April 21 and June 22, as early and late planting, respectively). The subplots included a factorial combination of three harvest dates (7 and 22 October, and 6 November) and 13 cultivars. The experiment was conducted at the Agricultural and Natural Resources Research and Education Center of West Azerbaijan Province, Miandoab Agricultural Research Station (Longitude, 46º 90' E; latitude, 36º 58' N; altitude 1314 m). This region features a mesic temperature regime (average annual soil temperature from 8 to 15 °C) and a xeric moisture regime (semi-arid). The soil was silty loam with a pH of 7.5-8.0 and an electrical conductivity of approximately 2.3 ds/m. The land was conventionally prepared, including plowing, disking, leveling, and planting. Each plot consisted of three rows, each 8 m long, with a line spacing of 0.5m and a plant distance of 18 cm. Agronomic management followed conventional practices, including irrigation, fertilization, cultivation operations, weeding, and pest control, which were conducted following local customs.
2.2. Measurements
The measured traits included root yield (RY) (MG ha-1), sugar yield (SY), white sugar yield (WSY) (MG ha-1), sugar content (SC), white sugar content (WSC) (%), amount of root impurities (including root sodium, potassium, and α-amino nitrogen content) (meq.100g beet-1), molasses sugar content (MS) (%), extraction coefficient of sugar (ECS) (%), and alkalinity coefficient (Alc), which were determined as follows:
2.3. Root Yield (RY)
After harvesting, weighing, and determining root yield, 30-35 roots were randomly selected. The roots were washed and prepared into pulp in the laboratory. The prepared pulp was then frozen and sent to the Sugar Technology Laboratory of the Sugar Beet Seed Breeding Research Institute in Karaj for further laboratory analysis to determine the desired traits. The methods used for measuring and determining the technological quality and quantity of sugar beets, as assessed in this study, are summarized below.
2.4. Sugar Content (SC)
To determine the sugar content, 177.8 cc of alkaline lead acetate solution was added to 26g of the prepared pulp. The mixture was thoroughly blended and filtered using a magnetic stirrer. A portion of the filtered solution was poured into a polarimeter using a funnel, and the sugar content was measured by observing the degree of deflection of polarized light [
17].
2.5. White Sugar Content (WSC)
This trait, which represents the amount of white sugar extracted in the factory, was calculated by subtracting the sum of molasses sugar and sugar waste in the factory (usually 0.6) from sugar content (SC). The closer the value of this trait to the sugar content, the less waste there is, resulting in a higher extraction coefficient.
2.6. Sugar Yield (SY)
This trait, equivalent to the amount of sugar produced per unit area of the field or the amount of sucrose stored in the roots of sugar beets, is considered the final yield in sugar beet cultivation. It was calculated by multiplying the root yield (RY) by the sugar content (SC) in MG ha
-1 [
18].
2.7. White Sugar Yield (WSY)
The white sugar yield, which is the same amount of sugar that can be extracted from sugar beets per unit area of the field, was obtained by multiplying the root yield (RY) by the white sugar yield (WSC) in MG ha
-1 [
18].
Pure sugar yield (t ha-1) = root yield (t ha-1) × [Sugar Content (%) − sugar efficiency loss (%)]
Sugar efficiency loss (%) = standard loss of molasses (%) + 0.6 (%)
Standard loss of molasses (%) = 0.012 × (K + Na) + 0.024(α-amino N) + 0.48
2.8. Root Impurities (Na, K, and α-Amino N)
Root impurities, including sodium, potassium, and α-amino nitrogen, were measured in milliequivalents per 100 g of root pulp. Another part of the filtered solution was analyzed using a Beet Analyzer to measure impurities. Sodium and potassium levels were determined using flame photometry, which compares the emission spectrum of lithium with previously calibrated standards [
19]. A small portion of the filtered solution was used to measure the α-amino nitrogen content, which included all organic nitrogen compounds such as amino acids. The Abbe number method measures these nitrogen compounds [
20].
2.9. Extraction Coefficient of Sugar (ECS)
The extraction coefficient of sugar (ECS) quantifies the efficiency with which white sugar can be extracted from the sucrose in sugar beet roots. It was calculated by dividing the white sugar content (WSC) by the sugar content (SC) × 100. A value closer to 100 indicates higher sugar production efficiency, implying that a more significant proportion of the sugar in the root is extractable and the root contains a lower percentage of impurities [
21].
2.10. Molasses Sugar Content (MS)
The content or percentage of non-extractable sugar from sugar beet root, known as molasses sugar (MS), was calculated based on the amount and content of root impurities such as sodium (Na), potassium (K), and nitrogen (N), according to the following equation [
22]:MS = 0.343 (Na + K) + 0.094 N − 0.29
2.11. Alkalinity Coefficient (Alc)
This trait was calculated as the total sodium and root potassium content ratio to root α-amino nitrogen content [
23].
2.12. Statistical Analysis
The collected data were subjected to statistical analyses, including variance analysis, genetic correlation analysis, multiple regression analysis, and path analysis. Stepwise multiple regression analysis was employed to investigate the effect of each trait on sugar yield (dependent variable), reduce the number of independent variables, and fit the best regression model. Path analysis based on genetic correlation coefficients was used to understand the relationships between traits and identify those that play the most important roles in sugar yield. Data were analyzed using Excel, SAS 9, and Path2 for statistical calculations.
3. Results
Before data analysis, normality tests were performed using the Kolmogorov-Smirnov and Shapiro-Wilk methods, which confirmed that the error distribution of the data for the traits was normal. Additionally, an F test indicated that the error variance between the two years was insignificant, demonstrating the uniformity of the error variance across the years. Consequently, the effect of the year was merged with the error, and a combined ANOVA was conducted as a factorial split test based on a randomized complete block design (RCBD).
Table 1 presents the results.
The year factor was highly significant for root yield, sugar yield, white sugar yield, sugar content, and white sugar content (p≤0.01) and significant for root α-amino nitrogen and root sodium content (p≤0.05). The planting date, as the main factor (A), significantly affected most of the measured traits. Specifically, it significantly affected the root yield, white sugar content, root α-amino nitrogen, and root potassium content (p≤0.01), as well as the sugar yield, white sugar yield, and sugar content (p≤0.05). The differences between the evaluated cultivars for all measured traits, including root yield, sugar yield, sugar content, root α-amino nitrogen, root sodium content (p≤0.01), white sugar yield, white sugar content, sugar extraction coefficient, root potassium content, alkalinity coefficient, and root molasses sugar content (p≤0.05).
The harvest date factor significantly affected all measured traits except the sugar extraction coefficient, root α-amino nitrogen, root potassium content, alkalinity coefficient, and root molasses sugar content. The interaction between planting date and harvest date for all measured traits was not statistically different, indicating that the influence of these factors was independent. Thus, the main effects on the measured traits were evaluated separately.
The interaction between planting date and cultivar was significant (p≤0.01) for root yield, sugar yield, white sugar yield, sugar content, white sugar content, and potassium content (p≤0.05). Some measured traits, including root yield, sugar yield, white sugar yield, sugar content, white sugar content, root α-amino nitrogen, and root potassium content, were significantly influenced by the interaction between the harvest date and cultivar (p≤0.05). The results suggest that the differences between cultivars vary depending on planting and harvest dates, demonstrating the interaction effect between these factors.
The 3-way interaction effect of cultivar × planting date × harvesting date was significant for root yield, sugar yield, white sugar yield, sugar content, and white sugar content (p≤0.05) but not for other traits. Because of the complexity of interpreting triple interactions, they are less frequently investigated, as they add complexity to the results.
The coefficients of variation ranged from a maximum of 38.85 for root sodium content to a minimum of 5.15 for sugar content. The average coefficients were related to sugar yield (SY), root α-amino nitrogen content (N), and molasses sugar content (MS) (
Table 1).
Table 2 presents the genetic correlation coefficients between different traits that showed significant negative correlations between root yield and sugar content (-0.60**), white sugar content (-0.54**), alkalinity coefficient (-0.52**), root potassium content (-0.47**), molasses sugar content (-0.45**), and root sodium content (-0.33). The results indicated significant positive correlations between root yield and sugar yield (0.86**), white sugar yield (0.79**), and sugar extraction coefficient (0.58**).
The white sugar yield was significantly positively correlated with sugar content (0.71**), white sugar content (0.65**), sugar extraction coefficient (0.59**), and root sodium content (0.50**). Furthermore, white sugar yield was significantly positively correlated with sugar content (0.77**), sugar extraction coefficient (0.76**), and white sugar content (0.71**), but negatively correlated with the alkalinity coefficient (-0.68**), root sodium content (-0.64**), root potassium content (-0.60**), and molasses sugar content (-0.59**). Root α-amino nitrogen content was also positively correlated with root sodium and potassium content (0.24* and 0.29*, respectively).
White sugar, sugar content, root α-amino nitrogen content, and root yield were included in the model (
Table 3). Other traits did not significantly affect the model, indicating that the difference in sugar yield between genotypes could be attributed to variations in these traits. The results showed that the independent variables—white sugar content, sugar content, root α-amino nitrogen content, and alkalinity coefficient—explained 24%, 28%, 11%, and 21% of the variance in sugar yield, respectively. The fitted model's cumulative coefficient of determination (R²) indicated that these independent variables explained 84% of the variance in the sugar yield.
Owing to the relatively high correlation between some variables in the model, the variance inflation factor (VIF) was used to check for the possible existence of collinearity. Only traits with a VIF value less than ten were included in the model. Measurement of the partial coefficients of standardized regression revealed that all independent variables in the model, except for root α-amino nitrogen content, increased sugar yield.
Path analysis was employed to interpret the results of genetic correlation coefficients and stepwise regression more comprehensively and to determine cause-and-effect relationships by specifying the components' direct and indirect effects. Path analysis showed that the correlation between traits and sugar yield could be attributed either to their direct effects on yield or to their indirect effects through other traits.
When the correlation between sugar yield and a trait is due to its direct effect, this indicates that the trait can be selected for improved yield. Conversely, if the correlation is due to an indirect effect through other traits, selection should focus on the trait that has caused the indirect effect.
Considering the important relationships between different traits in developing breeding programs and identifying suitable traits for selection, it is crucial to account for both direct and indirect effects. This approach prevents undesirable results of one-sided selection for crop traits without considering other traits. Thus, correlations between traits and their direct and indirect effects should be considered in breeding programs and the identification of superior genotypes, making path analysis necessary.
Path analysis was performed in this study with sugar yield as the dependent variable (effect) and traits in the stepwise regression model as the independent variables (cause). The results of the path analysis (
Table 4) indicated that the maximum and minimum direct effects were related to the sugar content (1.12) and root α-amino nitrogen content (-0.57), respectively.
Sugar content, with the highest direct and positive effect (1.12), had negative indirect effects through root α-amino nitrogen content (-0.35) and root yield (-0.43) and a positive indirect effect through white sugar content (0.37) on root yield (
Table 4,
Figure 1). Following sugar content, root yield had a direct and highly positive impact (1.03) on sugar yield. The indirect effect of root yield was decreased through white sugar content (-0.14) and sugar content (-0.21) but increased through root α-amino nitrogen content (0.18) (
Table 4,
Figure 1).
The residual effect in this study was 0.26 (
Table 4,
Figure 1), indicating how causal factors justified the variation of the dependent factor (sugar yield). The variables (white and sugar content, root α-amino nitrogen content, and root yield) explained approximately 74% of the variance in sugar yield. This is likely due to the relatively low correlations between traits, especially root α-amino nitrogen content and sugar yield. Additionally, other factors not considered in this study must be included in the path analysis to fully explain the variance in sugar yield.
4. Discussion
The results of the present study revealed substantial diversity among the cultivars. This diversity can be used in breeding programs to select superior cultivars.
Understanding the relationships between traits is crucial for crop breeding programs, as they indirectly influence selection success. Correlation analysis was conducted to evaluate the underlying relationships between different traits to identify the component characters that can be used for yield improvement.
In this study, the significant relationship between sugar yield and sugar content (0.71**) was likely due to the mathematical relationship between these two traits, where sugar content is a component of the formula used to derive the sugar yield. Because sodium is present in the root in an impure form, increasing sodium levels can increase sugar content but decrease sugar purity [
13]. Additionally, the significant positive correlation between root α-amino nitrogen content and roots' sodium and potassium content (0.24* and 0.29*, respectively) indicated that potassium does not inhibit nitrogen uptake. Potassium consumption in saline soils may reduce sodium damage to some extent [
24].
Furthermore, because of the significant negative correlation between root sodium content and white sugar content (-0.86**) and sugar content (-0.74**), it can be inferred that in high-sodium soils (e.g., saline soils), sugar content significantly decreases in the presence of potassium deficiency. The positive correlation between the extraction coefficient of sugar and white sugar yield (0.76**) can be attributed to the common components of these two traits, namely white sugar content. An increase in the white sugar content improved both traits.
The positive and significant correlation between the white sugar content and the extraction coefficient of sugar and white sugar yield supports this claim. Additionally, sugar content shares a common component with white sugar content and white sugar yield, so increasing white sugar content enhances both traits [
15]. The results of this study are consistent with those of previous studies [
13,
23,
25]. However, differences in the study conditions and genotypes affect the results of the correlations in different studies. Therefore, more studies are required to fill this information gap [
7,
23].
From stepwise regression analysis (
Table 3), it can be concluded that white sugar content, sugar content, root α-amino nitrogen content, and root yield are the most important traits affecting sugar yield (Y). Thus, selecting these traits in breeding programs will likely increase sugar yield. However, because the sugar content is part of the final yield, cultivars with higher root yield should be selected, considering the balance between this component and root yield to achieve a higher yield (sugar yield). In contrast with the results of the current study, Bashiri et al. (2015) concluded that sugar and root α-amino nitrogen content explained 61% of the total variance [
5]. The disagreement in the results could be due, at least in part, to differences in the genotypes and cultivars studied.
Path coefficient analysis provides insights into interrelationships by separating the correlation coefficients into the direct and indirect effects of the assessed traits. Plant breeders have used this method to identify useful traits as selection criteria for crop yield improvements. Sugar beet breeders use this method to estimate desirable features to enhance quantitative and qualitative yields. The relatively high and negative direct effect of root α-amino nitrogen content (-0.57) and its significant negative correlation with sugar yield (-0.33) at the 0.05 probability level (p≤0.05) indicate that this trait can be used as a selection criterion to improve sugar yield. However, the reverse relationship should also be considered. The negative direct effect of root α-amino nitrogen content (-0.57) (
Table 4, Fig. 1) suggests that the indirect effects of this trait through other traits are the main factors of its correlation with sugar yield. Therefore, indirect causal factors must be considered in the selection process. The results also indicate that the direct selection of traits with a positive direct effect (sugar content, white sugar content, and root yield) would enhance the breeding efficiency of sugar beet sugar yield.