Genetic analysis
Initially, was analyzed the data obtained from the experiment, which included various traits such as KY, KW, PH, EH, KR, KIN, KM, and EW. Bartlett's test was employed to assess the uniformity of error across three different environments. The homogeneity of the experimental error for all traits was confirmed by the results of Bartlett's test, allowing to conduct a combined analysis of variance based on the three environments, as detailed in Table 3. The analysis revealed significant environmental influences on all traits except the KR at a 1% probability level, highlighting the impact of environmental conditions on these traits. Notably, there was considerable genetic diversity among the experimental hybrids concerning all traits, also significant at the 1% probability level. The interaction between hybrid and environment was significant for all traits except for KR, underscoring the variability among genotypes across different conditions and the challenges faced by breeders in selection and release of varieties.
Table 3
Combined analysis of variance of studied traits based on the fourth method of Griffing
Source of variation
|
df
|
Mean square
|
|
|
Kernel yield
|
100 kernel weights
|
Plant height
|
Ear height
|
Kernel moisture
|
Ear wood
|
Kernel row
|
Kernel in row
|
Environment
|
2
|
64.56**
|
79248.89**
|
15294.62**
|
23017.71**
|
499.62**
|
493.79**
|
0.52ns
|
290.13**
|
Block (Environment)
|
6
|
11.29
|
690.19
|
1207.00
|
177.98
|
8.87
|
3.36
|
9.05
|
48.28
|
Hybrid
|
35
|
8.82**
|
4099.74**
|
1467.43**
|
745.39**
|
9.04**
|
17.26**
|
13.55**
|
71.87**
|
GCA
|
8
|
12.68**
|
13998.09**
|
4588.20**
|
2427.95**
|
26.42**
|
54.22**
|
43.81**
|
186.53**
|
SCA
|
27
|
7.67**
|
1166.90**
|
542.76**
|
246.95**
|
3.89**
|
6.30**
|
4.59**
|
37.90**
|
Environment × Hybrid
|
70
|
3.68**
|
868.92**
|
249.81**
|
99.10**
|
4.11**
|
3.22**
|
1.49ns
|
24.97**
|
Environment × GCA
|
16
|
6.58**
|
1456.51**
|
274.75**
|
166.22**
|
6.28**
|
6.66**
|
1.52ns
|
41.37**
|
Environment × SCA
|
54
|
2.82**
|
694.83**
|
242.42**
|
79.22**
|
3.47**
|
2.20**
|
1.48ns
|
20.11**
|
Error
|
210
|
1.33
|
326.34
|
136.90
|
34.36
|
2.05
|
1.07
|
1.24
|
11.17
|
CV (%)
|
8.20
|
6.01
|
5.49
|
5.26
|
6.75
|
6.29
|
6.51
|
9.03
|
Baker ratio
|
0.77
|
0.96
|
0.94
|
0.95
|
0.93
|
0.95
|
0.95
|
0.91
|
*, ** and ns: Significant at 5 percent, 1 percent and non- significant. GCA: General Combining Ability, SCA: Specific Combining Ability. |
The genetic diversity observed in the hybrids from crossing was dissected into additive and non-additive variance components using the fourth Griffing method 50. The mean squares for both GCA and SCA for all traits were significant at the 1% probability level. This significance indicates that both additive and non-additive gene effects are crucial in the expression and inheritance of these traits. Previous studies such as those by Vacaro, et al. 53, and Heidari, et al. 54, have emphasized the predominance of the additive effects over non-additive effects in controlling quantitative traits of maize, corroborating present findings.
Further analysis divided the hybrid-environment interaction into two components: GCA × environment and SCA × environment. Both interactions were significant for all traits except the KR at the 1% probability level (Table 3), pointing to environmental impacts on both additive and non-additive gene effects. This consistency with earlier studies by Mostafavi, et al. 55 and Glover, et al. 56 confirms the robustness of these results.
Baker's ratio was used to elucidate the additive and non-additive effects of genes more clearly. For traits where this ratio equals one, control by additive effects is complete. A ratio of 0.5 indicates equal influence of both effects, while a value below 0.5 suggests a dominant role for non-additive effects (dominance, over dominance, and epistasis). In this study, the high values of Baker's ratio for all traits suggest a predominant role of additive effects (Table 3). This finding aligns with other research, such as Saremirad and Mostafavi 27, which utilized Baker's ratio to determine gene effects in controlling various traits.
Figure 1 presents the values of the GCA effects for all studied traits. The GCA for KY varied, with a low of -0.97 in line KE 79017/3211 and a high of 0.42 in line KE 76009/311. Notably, the line KE 76009/311 displayed significant positive GCA for KY at the 5% probability level, identifying it as an optimal general combiner to enhance variance and selection efficiency. Consequently, this line, exhibiting positive GCA, can play a crucial role in boosting KY in selection-based breeding programs. Additionally, this line may increase the presence of genes with additive effects in breeding initiatives. Conversely, lines KE 79017/3211 and KE 77004/1 showed significantly negative impacts on KY at both the 1% and 5% probability levels, thus reducing this trait. For the KW, four lines KE 77008/1, NK79, KE 79017/3211, and KE 76009/311 demonstrated significant positive GCA at both the 1% and 5% probability levels, contributing to increased KW. In terms of PH, lines K 1264/5 − 1, NK79, and KE 75016/321 exhibited positive and significant GCA at the 1% probability level, whereas KE 79017/3211 showed significant negative GCA. For EH, lines KE 75016/321 and K 1264/5 − 1 had positive and significant GCA, while KE 79017/3211 and KE 77004/1 showed negative and significant GCA. A reduced EH can mitigate stem breaking and lodging, thereby preserving yield; using lines with negative and significant GCA can enhance the additive effects of genes and increase selection efficiency. Regarding KM, lines KE 75016/321 and KE 76009/311 showed positive GCA, while line K 1263/1 exhibited negative GCA. In the context of EW, lines KE 76009/311, KE 77004/1, and NK79 demonstrated positive GCA, whereas lines KE 77005/2, K 1264/5 − 1, and K 1263/1 indicated negative GCA at the 1% probability level. Lines K 1264/5 − 1, KE 76009/311, and KE 79017/3211 increased the KR in the ear, while KE 77008/1 and KE 77004/1 decreased them. Conversely, lines K 1263/1 and KE 77004/1 increased, and lines KE 76009/311, NK79, and KE 75016/321 decreased the KIN. Overall, lines K 1264/5 − 1 and K 1263/1 positively impacted the studied traits, while line KE 79017/3211 negatively affected these traits. Given that GCA reflects the additive effects of genes, parents with high GCA also possess substantial additive effects, beneficial for synthesizing new cultivars. Mostafavi, et al. 55 noted that lines K166B, K3615/2, and K3653/2 are among the top general combiners for KR, enhancing the genetic improvement of these traits. Another study by Dehghanpour 57 examined GCA and SCA in early inbred maize lines using the four Griffin method, concluding that additive effects largely influence trait expression, whereas non-additive effects impact traits such as the KR, leaf length, leaves number, PH, and EH.
Figure 2 details the SCA of hybrids for the traits analyzed. In terms of KY, SCA ranged from − 1.83 in the cross KE 77008/1×KE 76009/311 to 1.63 in KE 76009/311×KE 1264/5 − 1. The crosses K 1264/5 − 1×KE 76009/311, KE 77005/2×KE 75016/321, KE 77004/1×KE 77008/1, and KE 77008/1×KE 79017/3211 exhibited positive and significant SCA, suggesting their utility in enhancing non-additive gene effects for hybrid production. Similar findings were reported by Choukan and Mosavat 58 and Mostafavi, et al. 55. Dehghanpour and Ehdaie 59 found that lines KE75039 and K1264/5 − 1 exhibit high GCA stability for KY, minimally interacting with environmental factors; the cross KE75039×K2331 also demonstrated high SCA stability across environments. The KW varied significantly, from − 26.46 in the cross KE 77005/2×KE 76009/311 to 14.86 in KE 77005/2×K 1263/1, yet none of the hybrids showed significant positive SCA for this trait. Regarding PH, hybrids NK79×KE 77005/2 and KE 75016/321×KE 77004/1 showed positive and significant SCA at the 5% probability level, while K 1263/1×KE 75016/321, KE 76009/311×KE 77008/1, and KE 75016/321×KE 79017/3211 had negative SCA. For EH, only the crosses KE 77008/1×KE 77005/2 and KE 75016/321×KE 79017/3211 demonstrated negative and significant SCA. KM changes ranged from − 1.17 in the cross K 1264/5 − 1 NK79 to 1.19 in KE 77005/2×KE 79017/3211, though no significant crosses were noted. For EW, hybrids KE 77008/1×KE 76009/311 and KE 77008/1×KE 77005/2 exhibited positive and significant SCA. Hybrids K 1263/1×KE 77005/2, K 1264/5 − 1×KE 76009/311, and NK79×KE 77008/1 increased the KR, while KE 77005/2×KE 76009/311, KE 77005/2×KE 77008/1, K 1264/5 − 1×KE 79017/3211, and K 1263/1×NK79 increased the KIN. Research on maize has consistently underscored the importance of both GCA and SCA in trait control 55,58. Dehghanpour 60 highlighted that lines OH43/1–42, KE75039, and K1263/1 possess good GCA, which is beneficial for hybrid cultivar preparation. However, for producing single crosses, combinations such as KE75039×K1263/2 − 1, K1263/1×KE72012/1, K1263/1×K2331, and OH43/1/42×K1263/2 − 1, which exhibit considerable SCA, should be considered.
Exploration of correlation and graphical analysis of genotype by yield*trait
Figure 3 presents the results of Pearson's correlation analysis on various traits of experimental maize genotypes. The most notable positive and significant correlation, quantified at 0.73, was observed between PH and EH at a 1% probability level. This strong correlation indicates that an increase in PH is associated with an increase in EH. Other correlations were identified as follows: KM and EH (r = 0.57), KY and KW (r = 0.35), EW and KW (r = 0.33), KY and PH (r = 0.23), KM and KW (r = 0.15), EH and KR (r = 0.12), and KY and KR (r = 0.12). These correlations range from moderate to weak, yet are significant at the 1% and 5% probability levels. A significant negative correlation was most prominently observed between EW and EH (r=-0.57) at a 1% probability level, indicating that an increase in EW is associated with a decrease in EH. Additional correlations included KM and EH (r=-0.49), KW and KIN (r=-0.37), EW and PH (r=-0.30), KM and PH (r=-0.29), EW and KIN (r=-0.24), KM and KIN (r=-0.22), and EH and KW (r=-0.12). These correlations were moderate to weak and significant at the 1% and 5% probability levels. Previous studies, such as those by Akbari, et al. 61 and Nazeran, et al. 62, have also documented similar associations between KY and its components. These complex associations pose challenges in breeding programs. Breeders are thus encouraged to strike a balance while managing the negative correlations among key traits to enhance breeding outcomes, as discussed by Hassani, et al. 2. In this context, graphical methods such as GT and GYT biplots analyses have proven to be comprehensive and effective tools 2,47.
To enhance crop productivity, was conducted a GYT graphical analysis. the findings indicate that the first and second PCs account for 64.54% and 14.38% of the variance, respectively, summing up to 78.92% of the total variance in KY data. This high percentage of variance explained by the first two PCs validates the biplot diagrams' effectiveness in illustrating GYT relationships, as per Cruz, et al. 63. While a complex dataset may not be fully explained by the first two PCs 64, the biplot's utility is not negated 65. Shojaei, et al. 66 demonstrated that approximately 50% of the variation in GYT is captured by the first two PCs. Similarly, Faheem, et al. 48 found that the first two PCs account for nearly 85% of the variance, with the PC1 contributing about 74% and the PC2 about 11%. Hassani, et al. 2 study revealed that the first and second PCs explain 50.53% and 34.96% of the variance, respectively, totaling 85.49% of the yield data variations.
Analyzing the correlation between yield-trait combinations can reveal trait relationships, aiding in the development of new genotypes. Figure 4A illustrates that the acute angles between most vectors indicate a positive correlation among the yield-trait combinations. This positive correlation is attributed to yield being a primary component in these combinations 47. A high correlation between yield-trait combinations implies a high correlation between genotype rankings for those combinations. The graph shows no correlation between the KY*KIN and KY*KW combinations. However, there is a weak positive correlation between KY*KIN and combinations such as KY*PH, KY*KR, and KY*EH. Positive correlations are also observed between KY*KR and KY*EH, KY/EW and KY/KM, and KY*PH height with both KY*KR and KY*EH. These correlations suggest that increasing PH in a genotype may also enhance the KR, EH, and EW, while potentially reducing KM content. However, the correlation between KY*KIN and some traits complicates the breeding process, necessitating strategic planning.
In biplot of Fig. 4B, a polygon formed by connecting genotypes 18, 4, 12, 31, 23, 27, 36, and 20, which are furthest from the origin, is evident. Perpendicular lines drawn from the origin to the polygon sides facilitate the grouping of GYT combinations. Genotypes at the vertices of the polygon exhibit superior performance for the respective combinations, indicating their superiority in terms of those traits. Consequently, genotype 4 was identified as the best for the KY*KIN, while genotypes 12 and 31 excelled in other combinations such as KY/EW, KY/KM, KY*PH, KY*KR, KY*EH, and KY*KW. Genotypes not associated with any combination are considered less desirable and weaker. The polygonal biplot allowed us to categorize the yield-trait combinations into two groups: the first group includes KY*KIN, while the second group encompasses the remaining combinations, such as KY/EW, KY/KM, KY*PH, KY*KR, KY*EH, and KY*KW.
To evaluate the ranking of genotypes based on yield-trait combinations, we utilized the biplot diagram of average tester coordinates (Fig. 4C) and the superiority index (Table 4). In the biplot of tester coordinates, the average axis, represented by an arrow passing through the origin, indicates the mean of all yield-trait combinations. The axis perpendicular to the average axis measures the balance of genotypes across traits. Based on this, genotype 12 and then genotype 31 had the highest average KY and were identified as the best genotypes. In contrast, genotypes 36, 35, 34, and 21 had the lowest average KY. The results from the biplot of average tester coordinates were further confirmed by the superiority index. Various studies 2,47,48,66 and the findings of the present study demonstrate that the biplot diagram of average tester coordinates in the GYT biplot method is a valuable tool that provides insightful information about genotypes. Identifying the hypothetical ideal genotype is based on the concepts of balance and high KY. The desired genotype is one that exhibits the highest KY and maximum balance. Any genotype that is closest to this hypothetical ideal genotype is considered superior, while the genotype furthest from it is deemed the least favorable. Based on this biplot (Fig. 4D), genotypes 12, 11, 19, 9, and 2 were identified as the best genotypes due to their minimal distance from the hypothetical ideal genotype. Conversely, genotypes 36, 27, 35, 34, and 21 were named as undesirable genotypes because they had the greatest distance from the hypothetical ideal genotype (Fig. 4D). Hassani, et al. 2 employed the GYT graphic method to determine trait relationships and select superior genotypes, obtaining valuable results.
Table 4
Standardized genotype by kernel yield*trait data and superiority index for maize genotypes
Genotype
|
KY×KR
|
KY×KIN
|
KY×KW
|
KY×PH
|
KY×EH
|
KY/EW
|
KY/KM
|
Mean superiority index
|
1
|
-0.75
|
0.53
|
-1.28
|
-0.28
|
0.05
|
0.57
|
0.02
|
-0.16
|
2
|
1.24
|
0.91
|
0.37
|
0.60
|
0.46
|
1.55
|
1.41
|
0.93
|
3
|
0.78
|
2.28
|
1.11
|
1.07
|
0.74
|
0.79
|
1.59
|
1.19
|
4
|
-0.55
|
1.67
|
0.73
|
0.94
|
1.37
|
1.78
|
1.86
|
1.12
|
5
|
-0.05
|
-0.22
|
-1.06
|
-0.27
|
0.26
|
-0.54
|
-0.51
|
-0.34
|
6
|
0.25
|
0.77
|
-0.67
|
-0.28
|
0.65
|
-0.27
|
-0.16
|
0.04
|
7
|
-0.13
|
1.55
|
-0.74
|
-0.11
|
-0.64
|
-0.32
|
0.33
|
-0.01
|
8
|
-0.03
|
-0.19
|
-0.77
|
-1.38
|
-1.11
|
-0.05
|
-0.16
|
-0.53
|
9
|
1.36
|
0.18
|
0.24
|
0.70
|
0.68
|
2.18
|
1.47
|
0.97
|
10
|
0.90
|
-0.70
|
0.08
|
0.89
|
0.84
|
0.08
|
1.17
|
0.47
|
11
|
0.05
|
0.38
|
0.82
|
1.24
|
1.34
|
1.57
|
1.28
|
0.95
|
12
|
3.17
|
1.37
|
1.55
|
2.44
|
2.10
|
1.65
|
1.32
|
1.94
|
13
|
-0.34
|
-1.34
|
-0.45
|
0.60
|
1.05
|
-0.15
|
-1.50
|
-0.30
|
14
|
0.10
|
0.01
|
-0.97
|
0.40
|
-0.57
|
-0.54
|
-0.07
|
-0.23
|
15
|
0.65
|
0.05
|
-0.75
|
-0.73
|
-0.77
|
-0.27
|
-0.22
|
-0.29
|
16
|
-0.22
|
-0.93
|
-0.12
|
0.52
|
0.20
|
-0.36
|
0.08
|
-0.12
|
17
|
-0.42
|
0.79
|
0.15
|
-0.60
|
-0.80
|
-0.07
|
-0.19
|
-0.16
|
18
|
0.04
|
1.39
|
-0.56
|
0.07
|
-0.49
|
0.72
|
0.68
|
0.26
|
19
|
0.93
|
0.72
|
0.50
|
0.32
|
1.19
|
1.63
|
0.21
|
0.79
|
20
|
-1.29
|
0.47
|
-2.03
|
-0.68
|
-0.86
|
-0.62
|
-0.52
|
-0.79
|
21
|
-0.87
|
-1.67
|
-1.02
|
-2.04
|
-1.52
|
-0.01
|
-1.72
|
-1.27
|
22
|
-0.34
|
-0.56
|
1.14
|
-0.08
|
0.08
|
-0.28
|
0.46
|
0.06
|
23
|
0.75
|
-0.76
|
1.49
|
0.68
|
0.71
|
-1.24
|
0.00
|
0.23
|
24
|
-0.34
|
-0.54
|
0.47
|
0.58
|
0.64
|
-0.81
|
-0.47
|
-0.07
|
25
|
-0.50
|
0.48
|
-0.04
|
0.22
|
-0.58
|
-0.74
|
0.18
|
-0.14
|
26
|
0.02
|
-0.08
|
1.36
|
-0.37
|
-0.61
|
-0.04
|
-0.27
|
0.00
|
27
|
-0.71
|
-1.82
|
0.38
|
-1.24
|
-1.08
|
-1.87
|
-1.73
|
-1.15
|
28
|
-1.03
|
-1.04
|
0.30
|
0.57
|
0.57
|
0.09
|
-0.65
|
-0.17
|
29
|
-0.80
|
0.72
|
0.74
|
0.65
|
0.08
|
0.30
|
1.10
|
0.40
|
30
|
-0.53
|
-0.71
|
1.38
|
-0.55
|
-0.65
|
-0.16
|
0.12
|
-0.16
|
31
|
2.06
|
0.57
|
1.73
|
1.82
|
2.21
|
0.63
|
0.61
|
1.37
|
32
|
-0.49
|
-0.27
|
-0.27
|
-0.71
|
-1.04
|
-1.37
|
-0.81
|
-0.71
|
33
|
1.08
|
-1.36
|
0.59
|
-0.72
|
-0.51
|
0.08
|
-0.23
|
-0.15
|
34
|
-1.68
|
-0.66
|
-1.63
|
-0.55
|
-0.77
|
-1.67
|
-1.72
|
-1.24
|
35
|
-0.62
|
-1.22
|
-1.44
|
-1.83
|
-1.44
|
-1.06
|
-1.33
|
-1.28
|
36
|
-1.66
|
-0.76
|
-1.33
|
-1.92
|
-1.75
|
-1.18
|
-1.61
|
-1.46
|
Mean
|
0.00
|
0.00
|
0.00
|
0.00
|
0.00
|
0.00
|
0.00
|
-
|
Standard deviation
|
1.00
|
1.00
|
1.00
|
1.00
|
1.00
|
1.00
|
1.00
|
-
|
KY: Kernel yield, KW: 1000 kernel weight, PH: Plant height, EH: ear height, KR: Kernel row number, KIN: Kernels in row, KM: Kernel moisture, EW: Ear wood. |
Multi-trait stability index
The factor analysis was performed based on PCA, and the results were interpreted after varimax rotation. Table 5 presents the outcomes of the factor analysis. Factors with eigenvalues greater than one were selected, and their variance, expressed as a percentage, indicates their importance in explaining the overall data variability. In this analysis, three independent factors accounted for a total of 65.89% of the data variations. The first factor explained 33.11% of the total variance and had an eigenvalue of 2.65. This factor had large negative factor coefficients for KY, KR, PH, and EH. The second factor, with an eigenvalue of 1.55, accounted for 19.34% of the variations and included height positive factor coefficients for the KIN and EW, along with a negative coefficient for the KW. The third factor, with an eigenvalue of 1.08, explained 13.44% of the total data variance and had a height negative factor coefficient for KM. The MTSI of the studied genotypes was calculated based on the factor scores of the three mentioned factors. According to this index, a lower value indicates a genotype that is closer to the stable ideal genotype. Conversely, a higher MTSI value for a genotype suggests a greater distance from the ideal stable genotype, making it less desirable for selection. In Fig. 5A, the experimental genotypes were ranked from the highest to the lowest MTSI value. The genotype with the highest index value is placed at the center, while the genotype with the lowest value is positioned at the outermost circle. Applying a selection pressure of 25%, genotype 10 ranked first, followed by genotypes 9, 13, 11, 1, 2, and 16, which were identified as the most ideal stable genotypes in terms of all traits. Comparing the trait values of the selected genotypes based on the MTSI with those of all experimental genotypes revealed that the average values of KY, KR, PH, and EH increased in the selected genotypes. However, the average values of the KIN, KW, EW, and KM decreased. Increases in KY, KR, PH, and EH were considered favorable, but decreases in the KIN and KW are not among the goals of maize breeding programs. Reducing EW and KM is desirable, and the selected genotypes showed a significant reduction in these traits.
Table 5
Prediction of selection differential for studied traits based on MTSI index
Variable
|
FA1
|
FA2
|
FA3
|
Communality
|
Uniqueness’s
|
Goal
|
h2
|
SD (%)
|
SG (%)
|
KY
|
-0.59
|
0.15
|
0.39
|
0.52
|
0.48
|
increase
|
0.58
|
0.73
|
0.43
|
KR
|
-0.56
|
-0.04
|
-0.13
|
0.33
|
0.67
|
increase
|
0.89
|
2.74
|
2.44
|
KIN
|
0.42
|
0.81
|
0.06
|
0.83
|
0.17
|
increase
|
0.65
|
-2.14
|
-1.40
|
KW
|
-0.11
|
-0.85
|
-0.03
|
0.74
|
0.26
|
increase
|
0.79
|
-1.38
|
-1.09
|
PH
|
-0.84
|
-0.23
|
-0.08
|
0.77
|
0.23
|
increase
|
0.83
|
6.06
|
5.02
|
EH
|
-0.80
|
-0.10
|
0.18
|
0.68
|
0.32
|
increase
|
0.87
|
7.78
|
6.74
|
EW
|
-0.19
|
0.67
|
-0.31
|
0.58
|
0.42
|
decrease
|
0.81
|
-8.34
|
-6.78
|
KM
|
0.01
|
0.12
|
-0.90
|
0.82
|
0.18
|
decrease
|
0.54
|
-4.04
|
-2.20
|
Eig.
|
2.65
|
1.55
|
1.08
|
-
|
-
|
-
|
-
|
-
|
-
|
Var. (%)
|
33.11
|
19.34
|
13.44
|
-
|
-
|
-
|
-
|
-
|
-
|
Cum. Var. (%)
|
33.11
|
52.45
|
65.89
|
-
|
-
|
-
|
-
|
-
|
-
|
KY: Kernel yield, KW: 1000 kernel weight, PH: Plant height, EH: ear height, KR: Kernel row number, KIN: Kernels in row, KM: Kernel moisture, EW: Ear wood, Eig: Eigen value, Var: Variance, Cum. Var: Cumulative variance. |
Overall, the selected genotypes resulted in selection differential and favorable selection progress in all traits except the KIN and the KW (Table 5). In the selected genotypes, all traits except KM and KY had high heritability. Figure 5B illustrates the strengths and weaknesses of the selected genotypes based on the contribution of each factor to the MTSI. According to this diagram, a lower contribution by a factor (closer to the outer edge) indicates that the traits within that factor are closer to the stable ideal state. The dashed line represents the theoretical value if all factors play an equal role. Given that each genotype is closer to the ideal genotype for factors where it has a smaller contribution in terms of traits within that factor, genotypes 13, 10, and 1, which had the lowest value in the first factor, are close to the ideal genotype for KY, KR, PH, and EH, which had the highest factor coefficients in this factor. Genotypes 11, 2, and 9 had the lowest contribution of the second factor, making them very close to the ideal genotype in terms of the KIN, EW, and the KW. Four genotypes- 10, 16, 11, and 9- had the largest share of the third factor, indicating they are very close to the ideal genotype in terms of KM. In other words, these genotypes have low KM. Sharifi, et al. 67 used the MTSI to evaluate yield and other agronomic traits in a set of rice genotypes and showed that this index effectively identifies superior genotypes in terms of yield stability and other agronomic traits. Based on the results of the MTSI index, Rajabi, et al. 15 introduced five genotypes as stable genotypes under conditions infected with rhizomania disease. Taleghani, et al. 68 identified four ideal genotypes simultaneously in terms of different traits using the MTSI index. These results are consistent with the findings obtained from this study regarding the efficiency of the MTSI index in identifying superior genotypes.