Performance and Stability for Grain Yield and Its Components of Some Rice Cultivars under Various Environments

Abstract

Refine current agricultural practices considering environmental changes are crucial for finding tolerant rice varieties that can meet the demands of human consumption. To this end, stability analysis assesses a crop genotype’s ability to adapt to various conditions. Therefore, the objective of this study was to (1) examine the interaction between rice genotypes and environmental conditions; (2) evaluate the stability of twelve rice genotypes using various stability methods; (3) identify representative environments for multi-environment testing; and (4) determine superior genotypes for specific environments. The evaluated rice cultivars were Sakha 101, Sakha 104, Sakha 105, Sakha 106, Sakha 107, Sakha 108, Giza 177, Giza 178, Giza 179, Giza 182, Egyptian Yasmine, and Sakha super 300. The experiment followed a strip-plot design, with three replications. The findings revealed significant differences among the rice varieties across various environments for the majority of the assessed characteristics. The joint regression analysis of variance demonstrated highly significant differences among rice cultivars for all the studied traits in terms of genotype-by-environment interaction (G × E). The statistical significance of the interaction between genetic and environmental factors was evident for all variables demonstrating heritable variation among the rice cultivars, specifically Sakha 108, Sakha 104, Giza 177, and Giza 178, concerning grain yield per feddan. These rice cultivars exhibited stability parameters that were not significantly different from unity for the regression coefficient (b_i) and from zero for the deviations from regression (S²d_i) for those traits. Overall, stability criteria are essential for ensuring reliable rice production, meeting human consumption, advancing genetic improvement, and promoting environmental sustainability in agriculture.

1. Introduction

Rice, a vital cereal that is cultivated as a staple food crop, faces challenges in meeting human consumption demands, which are expected to increase by an additional 112 million metric tons by the year 2035 [1]. To address this growing need, it is imperative to optimize existing farming methods, manage water resources effectively, and explore opportunities for developing new rice varieties [2].

Rice cultivation spans various agroecology and cropping systems, including rain-fed highland and lowland areas, irrigated fields, and mangrove ecosystems [3]. While upland rice constitutes a relatively minor portion (approximately 11%) of global rice output, it holds significant importance in local rice production in both arid and semiarid regions [4,5]. The challenges posed by limited water resources for lowland rice agriculture are exacerbated by the impacts of climate change. In response, upland rice emerges as a potentially sustainable solution to tackle food security concerns, as suggested by [3,5,6].

Before commercial release, plant breeders typically evaluate the performance of the new varieties across various environments to assess their stability and consistency [7]. Genotypes that demonstrate consistent performance across diverse environments are preferred, but the complexity of genotype–environment interactions poses challenges in genotype selection [8,9,10]. Therefore, it is crucial to employ statistical methods to evaluate genotype stability before commercial release [11,12]. To this end, the applied statistical methods for assessing genotype stability can be classified into parametric and nonparametric stability statistics [13,14]. Parametric stability statistics encompass various univariate and multivariate methods. Univariate methods comprise techniques such as Wricke’s equivalence (Wi2) [15], Shukla’s stability variance (S²_Shu.) [16], coefficient of variance (CV) [17], Environmental variance (S²_Env.) [18], Mean-variance component (q) [19], GE variance component (q’) [20], and regression coefficient (b_i) [21], among others. On the other hand, multivariate approaches include the additive main effects and multiplicative interaction (AMMI) model [22] and the GGE biplot method [23], which can accurately predict genotype-by-environment interactions and identify optimal genotypes across diverse testing environments [24].

The objectives of this study are to (1) study the interaction of genotypes and the environment; (2) determine the stability of rice genotypes with various stability methods; (3) determine representative environments for multi-environment testing; (4) determine superior genotypes in specific environments; and (5) identify high-yielding rice genotypes with consistent performance across six distinct locations over two consecutive years (2020–2021) (three sowing dates during two seasons). This is achieved by evaluating the effectiveness of various univariate and multivariate stability criteria. The identification of high-yielding rice genotypes is essential for ensuring food security and enabling adaptation to climate change in rice-producing regions worldwide.

2. Materials and Methods

This study, conducted at the rice research farm Sakha Agricultural station in Kafr Elsheikh, Egypt, aimed to assess 12 rice cultivars over the 2020 and 2021 seasons, which were subjected to three different sowing dates.

2.1. Rice Cultivars and Experimental Design

Twelve rice cultivars, namely Sakha 101, Sakha 104, Sakha 105, Sakha 106, Sakha 107, Sakha 108, Giza 177, Giza 178, Giza 179, Giza 182, Egyptian Yasmine, and Sakha super 300, were evaluated. The experiment was set using a split-plot design with three replications. Sowing dates were allocated in the main plots, and rice varieties were applied in the sub-plots. The plot size was 2 m × 5 m (10 m²).

The chemical and mechanical analyses of the soil and organic matter were conducted following the methods described by [25] at the Agricultural Research Center, Ministry of Agriculture, Egypt. The chemical and physical properties of the soil at the experimental site, measured at a depth of 0–30 cm, are shown in

Table 1. Mechanical and chemical soil analyses of the experimental site during the 2020 and 2021 seasons.

Soil Properties	2020	2021
Mechanical: Clay %	57.22	56.56
Silt %	28.20	27.31
Sand %	14.58	16.13
Texture	Clayey	Clayey
Chemical:
Organic Matter (O.M)%	1.54	1.52
pH (1:2.5 soil suspension)	8.34	8.45
Ec (ds.m−1)	3.34	3.20
Total N (ppm)	467.50	464
Available P (ppm)	14.69	14.32
Available K (ppm)	390	381
Soluble anions, meq.L−1:
${\mathbf{HCO}}_3^{-}$	6.38	6.20
Cl−	9.25	9.10
Soluble Cations, meq.L−1:
Ca2+	10.86	10.72
Mg2+	5.15	5.13
Na+	2.13	2.09
K+	15.71	15.05
Total carbonate %	14.23	14.14

.

Weather data for the two seasons were collected from the Agriculture Research Center, Field Crops Research Institute, and Rice Research and Training Center in Sakha, Kafr El-Sheikh, Egypt (

Table 2. The monthly minimum and maximum air temperature (°C) and relative humidity (%) at Sakha Agricultural research station during the 2020 and 2021 rice growing seasons.

Month	Day	2020					2021
		Air Temp.		RH %		Rain (mm/day)	Air Temp.		RH %		Rain (mm/day)
		Min	Max	7.30	13.30		Min	Max	7.30	13.30
May	1–10	20.98	28.89	76.0	46.9	0	25.4	30.7	77.3	47.7	0
	11–20	24.36	32.93	66.6	41.5	0	26.3	30.3	78.0	44.3	0
	21–31	22.93	29.72	71.5	48.9	0	25.8	30.9	77.9	44.7	0
June	1–10	26.7	33.4	67.3	43.5	0	27.7	32.6	78.7	46.7	0
	11–20	25.5	33.1	76.0	44.3	0	27.8	31.9	81.5	54.2	0
	21–30	26.8	34.2	83.9	52.1	0	28.8	32.9	80.6	53.2	0
July	1–10	26.5	33.9	82.5	57.7	0	29.4	34.7	84.6	57.2	0
	11–20	25.9	33.7	83.0	54.2	0	29.1	35.2	83.4	54.9	0
	21–31	26.0	32.9	81.7	58.6	0	28.4	34.4	85.6	59.1	0
August	1–10	26.2	34.3	84.8	57.6	0	28.8	34.2	86.8	57.2	0
	11–20	26.3	33.0	82.7	56.8	0	29.2	34.9	84.6	58.3	0
	21–31	25.5	36.2	85.0	54.4	0	27.4	35.6	86.2	50.5	0
September	1–10	24.3	33.27	85.4	53.7	0	25.9	33.1	86.0	49.8	0
	11–20	24.93	33.42	82.6	51.3	0	26.0	33.6	87.1	49.5	0
	21–30	23.63	31.09	82.2	50.1	0	26.7	31.5	85.0	51.7	0
October	1–10	20.3	31.5	85.6	50.0	0	24.5	30.0	83.3	53.1	0
	11–20	23.0	29.84	78.4	56.0	0	23.6	27.9	81.5	56.6	0
	21–31	21.90	28.04	83.2	59.8	0	23.9	28.1	78.5	54.6	0

).

2.2. Field Preparation and Planting

The experiment was carried out in plots with seven rows, each 5 m in length, and a plant spacing of 20 m × 20 cm. Bed nursery and permanent field preparation followed the guidelines provided by the Rice Research and Training Center (RRTC). Phosphorus fertilizer in the form of mono-super phosphate (15%) at the rate of 36 kg P₂O₅/ha and potassium in the form of potassium sulfate (48% K₂O) at the rate of 57 kg K₂O/ha were applied during land preparation. Nitrogen in the form of urea (46.5% N) was applied in two splits, i.e., two-thirds as basal application incorporated into the soil immediately before flooding and one-third after 30 days from the first dose.

2.3. Transplantation and Spacing

Seedlings were transplanted from the nursery to the permanent field after 25 days.

Transplantation was performed into sub-sub-plots using plant spacing (20 × 20, cm) for distances between hills and rows. Seedlings were transplanted into hills at a density of 3–4 seedlings per hill with the following schedule: E1 = the first sowing date in the first season 20 April 2020; E2 = the second sowing date in the first season 10 May 2020; E3 = the third sowing date in the first season 30 May 2020; E4 = the first sowing date in the second season 20 April 2021; E5 = the second sowing date in the second season 10 May 2021; and E6 = the third sowing date in the second season 30 May 2021.

2.4. Studied Characteristics

This study proposed the evaluation of additional attributes, including the number of filled spikelets per panicle, seed set percentage, 1000-grain weight, and grain yield per plant. Overall, this study utilized a rigorous experimental design and methodology to evaluate the performance of different rice cultivars under varying sowing dates during the two seasons, providing valuable insights into rice production in the region.

2.5. Statistical Analyses

The stability and adaptability of rice genotypes were determined using parametric and nonparametric stability in addition to multivariate analysis. Parametric stability included linear regression coefficient (b_i) and AMMI stability indexes, whereas nonparametric stability included S_i. Multivariate analysis was performed using AMMI and GGE biplot.

Stability studies were conducted using the Eberhart and Russell model [26] with metan package in R software 2021, with particular attention to the G × E factor. A genotype with a regression coefficient of 1 (b_i = 1) and a deviation not statistically significant from zero (S²d_i = 0) is deemed stable, exhibiting a unity response.

Mean comparisons among genotypes were conducted using the LSD test, employing a significance threshold of 5%. To calculate variance components, ANOVA was performed using a split-plot design. The GenStat statistical tool (12th edition) facilitated this analysis. A Bayesian network is a probabilistic directed acyclic graph, where nodes represent variables such as traits and edges represent causal or conditional dependencies among these variables. The aim of analysis algorithms in Bayesian networks is to infer the relations of conditional or causal dependence among network variables based on the Bayes theorem. So, Bayesian networks belong to the class of structural network models. Bayesian network was built using the package bnlearn in R software.

3. Results

Substantial variations were observed in the mean square values for the number of filled grains per panicle, seed set percentage, 1000-grain weight, and grain yield per panicle among the different Egyptian rice types and sowing dates (

Table 3. Mean square of seed set (%), 1000-grain weight, and grain yield/plant of the 12 rice cultivars.

SOV	Df	MS
		Number of Filled Grains/Panicle	Seed Set (%)	1000-Grain Weight	Grain Yield/Plant
Years (Y)	1	87.53 **	296.83 **	1.38	343.27 **
Date (D)	2	2420.52 **	590.29 **	5.59 **	204.71 **
Y × D	2	53.78 **	132.24 **	0.00	21.80 **
Rep/D/Y = (Ea)	9	25	1.0	0.50	2.0
Genotypes (G)	11	3310.30 **	71.60 **	71.22 **	18.92 **
G × Y	11	105.81 **	15.30 **	0.05	3.84 **
G × D	22	129.95 **	9.21 **	3.01 **	0.66
G × Y×D	22	37.35 **	8.13 **	0.05	1.36
Error/D/years (Pooled Error) = (Eb)	135	14.38	0.71	0.04	0.51

** indicate p < 0.01, respectively.

). Variations in the habitats and genotypes had a significant impact on all the tested traits. Moreover, clear interactions were observed between years and environments, as well as between genotypes and years, across all traits except for the 1000-grain weight, where the disparity was not statistically significant. Furthermore, the analysis revealed a substantial and significant mean square of interaction across genotypes, years, and environments for the number of filled grains per panicle and seed set percentage (Table 3). The results demonstrated that the investigated genotypes had varying responses to environmental variables, highlighting the significance of evaluating genotypes across multiple environments to identify the most well-suited genotypes for specific situations. The environmental mean squares were the primary contributor to the overall mean squares for all traits. Furthermore, the variations attributed to different environments were greater than the variations resulting from the interactions between genotypes and environments for all the observed traits. Hence, the primary factor influencing the variations in the performance of rice genotypes in these trials was the environment, rather than the differences arising from genotype-by-environment interaction.

3.1. Environment Effect

The mean number of filled grains per panicle decreased by 11.6 grains when comparing the sowing dates condition to the third sowing date condition (

Table 4. Means of the 12 rice genotypes over years and environments.

Item	No. of Filled Grains/Panicle	Seed Set %	1000-Grain Weight (g)	Grain Yield/Plant (g)
First year	137.5	93.9	27.3	45.6
Second year	136.2	91.6	27.4	43.1
Env. 1	142.9	95.8	27.7	45.8
Env. 2	136.2	92.4	27.3	44.8
Env. 3	131.3	90.1	27.1	42.5
Mean overall	136.82	92.76	27.36	44.36
LSD at 0.05	1.54	0.33	0.07	0.43

). When comparing the second and third sowing dates to the first sowing date, the seed set percentage, 1000-grain weight (in grams), and grain yield per plant were lower. The results demonstrated that the diverse growing conditions led to a significant difference in average yields, ranging from 45.0 g/plant under favorable initial planting date conditions to 42.5 g/plant under the third stress condition (

Table 5. AMMI analysis for the 12 rice genotypes.

Source	Df	Sum Sq	Mean Sq	F Value	Pr (>F)	Proportion	Accumulated
ENV	5	7.96 × 102	1.59 × 102	111.18	1.27 × 10−9	-	-
REP(ENV)	12	1.72 × 101	1.43 × 100	2.77	2.18 × 10−3	-	-
GEN	11	2.08 × 102	1.89 × 101	36.56	8.63 × 10−35	-	-
GEN:ENV	55	8.67 × 101	1.58 × 100	3.05	1.01 × 10−7	-	-
PC1	15	4.79 × 101	3.20 × 100	6.18	0.00 × 100	55.3	55.3
PC2	13	2.25 × 101	1.73 × 100	3.34	2.00 × 10−4	25.9	81.2
PC3	11	1.19 × 101	1.08 × 100	2.09	2.52 × 10−2	13.7	94.9
PC4	9	4.38 × 100	4.87 × 10−1	0.94	4.93 × 10−1	5.1	100
PC5	7	4.07 × 10−3	5.80 × 10−4	0	1.00 × 100	0	100
Residuals	132	6.83 × 101	5.18 × 10−1	-	-	-	-
Total	270	1.26 × 103	4.68 × 100	-	-	-	-

).

The AMMI 1 biplot (Figure 1a) represents the relationship between grain yield and IPCA 1. In this plot, rice genotypes that fall on the same vertical line indicate similarity in yield, while those on the same horizontal line indicate similarity in interaction pattern. According to the AMMI 1 biplot, the rice genotypes Sakha104, Sakha105, Sakha106, Sakha108, Giza177, Giza178, and SR300 exhibited a high level of stability in grain production, and they are considered extensively adapted lines (Figure 1a). The rice genotypes Sakha101, Sakha107, Giza179, Giza182, and Egyptian Yasmine exhibited relatively low stability in grain yield due to their considerable distance from the origin and their specific adaptability to particular environments. The genotype S101 exhibited superior performance in the E1 and E2 environments, while the genotypes G182 and EY were specifically adapted to environments E3 and E5, respectively. Furthermore, the S106, S108, and G178 genotypes were specifically identified as being well suited to environment E4. In the AMMI 2 biplot (Figure 1b), the scores of IPCA1 plotted against IPCA2 show the stability of the genotypes across different environments. Environments E1, E2, and E4 had relatively small radii and exhibited weak interacting forces, whereas environments E3, E5, E6, and E7 had long radii and exerted substantial interaction. Genotypes S104, S105, and G178 were located in close proximity to the origin, indicating their insensitivity to environmental interaction forces. Conversely, the remaining genotypes were situated further from the zero line, signifying their heightened responsiveness. The genotypes S106 and EY were found to be the most well suited for site E6, whereas the genotypes S108 and SR300 showed a close similarity in environment E4.

3.2. Mean Performance

The average number of filled grains per panicle varied between 121.06 for Giza 179 and 168.0 grains for Giza 178, as shown in

Table 6. The average number of filled grains/panicle and seed sets for the twelve rice cultivars in different environments.

Genotypes	No. of Filled Grains/Panicle							Seed Set %
	E1	E2	E3	E4	E5	E6	Mean	E1	E2	E3	E4	E5	E6	Mean
Sakha 101	144.3	142.0	134.0	145.3	139.0	135.0	139.94	97.7	96.9	94.8	97.9	93.2	90.5	95.2
Sakha 104	141.7	138.0	130.3	140.4	133.3	129.0	135.46	97.3	96.3	94.4	96.7	93.7	90.4	94.8
Sakha105	130.0	124.3	118.3	129.7	124.0	119.0	`124.22	96.7	96.3	93.7	98.6	92.7	88.8	94.5
Sakha 106	127.3	128.3	123.3	131.0	129.0	126.0	127.50	95.1	94.1	91.5	94.6	94.6	93.8	93.9
Sakha 107	126.0	123.3	120.0	127.0	127.3	123.5	124.51	93.3	92.6	90.5	93.8	89.8	87.6	91.3
Sakha 108	152.7	149.0	143.3	151.0	147.0	140.4	147.23	95.5	94.7	93.1	94.7	90.7	86.7	92.6
Giza177	141.7	136.7	129.3	142.7	133.3	128.1	135.30	97.4	96.6	94.3	99.5	90.6	87.1	94.3
Giza178	188.3	163.3	173.0	189.3	150.0	144.0	168.00	95.8	94.2	92.2	98.9	86.2	81.4	91.4
Giza179	129.7	124.3	120.3	130.7	128.3	125.1	126.41	96.7	93.8	93.7	97.3	94.5	91.9	94.7
Giza 182	123.7	117.7	115.3	124.7	124.0	121.0	121.06	91.5	90.7	88.5	92.8	90.4	88.0	90.3
Egyptian Yasmine	149.3	152.7	138.0	150.3	143.7	135.3	144.89	93.2	91.2	89.8	93.8	88.0	83.5	89.9
Sakha super 300	155.7	150.3	142.7	156.7	141.0	137.3	147.3	94.7	92.2	90.8	95.1	84.6	85.1	90.4
Sakha 101	142.53	137.49	132.32	143.23	134.99	130.31	136.82	95.41	94.13	92.28	96.14	90.75	87.90	92.78
LSD 0.05	1.3	15.1	1.7	1.52	3.64	2.01	0.14	0.14	1.32	0.11	0.93	2.40	2.01	0.55
LSD 0.01	1.7	20.6	2.4	2.06	4.95	2.73	0.19	0.19	1.79	0.15	1.27	3.26	2.74	0.72

. The seed set percentages of the twelve rice genotypes exhibited substantial variations, as indicated by the data presented in Table 6. The seed set percentage was highest for Sakha 101 and lowest for E. Yasmin.

The mean 1000-grain weight was highest for Sakha 106, Sakha 105, and Sakha 108, with values of 29.1, 29.0, and 29.0, respectively. The genotype Giza 178 (22.1) had the lowest value among all genotypes (

Table 7. Average 1000-grain weight and grain yield/plant for the twelve rice cultivars over environments.

Genotypes	1000-Grain Weight							Grain Yield/Plant
	E1	E2	E3	E4	E5	E6	Mean	E1	E2	E3	E4	E5	E6	Mean
Sakha 101	28.6	28.4	28.3	28.8	28.6	28.5	28.5	47	47	44.9	46.6	44.3	41.0	45.0
Sakha 104	28.3	28.2	27.7	28.5	28.4	27.9	28.2	47.1	47.1	44.2	46.8	43.1	40.9	44.9
Sakha105	28.1	28.1	28.9	28.3	28.2	29.0	28.4	46.8	46.8	43.8	45.2	43.2	40.8	44.4
Sakha 106	28.3	28.5	28.9	28.5	28.7	29.1	28.6	46.9	46.8	43.2	44.7	43.5	41.6	44.4
Sakha 107	27.5	27.4	27.2	27.7	27.5	27.3	27.4	44.0	44.0	43.5	42.7	41.1	38.4	42.3
Sakha 108	29.1	29.0	28.9	29.3	29.2	29.0	29.1	47.8	47.8	46.9	46.6	44.7	41.5	45.6
Giza177	28.3	28.3	27.9	28.5	28.4	28.1	28.2	46.0	46.0	43.7	43.8	42.3	41.2	43.8
Giza178	22.1	22.0	21.9	22.2	22.2	22.1	22.1	47.3	47.3	44.7	46.1	44.7	42.2	45.4
Giza179	27.5	27.5	27.1	27.7	27.6	27.2	27.4	47.0	47.0	44.1	46.9	45.0	43.3	45.5
Giza 182	27.3	24.2	22.6	27.5	24.4	22.8	24.8	46.6	46.6	44.0	43.7	40.2	39.3	43.4
Egyptian Yasmine	27.7	27.4	27.3	27.9	27.6	27.4	27.5	45.8	45.8	42.6	44.6	42.6	42.1	43.9
Sakha super 300	28.1	27.9	27.8	28.2	28.1	28.0	28.0	45.0	45.0	43.5	44.5	42.8	40.2	43.5
Sakha 101	27.58	27.24	27.04	27.76	27.41	27.20	27.35	46.45	46.44	44.02	44.18	43.13	41.04	44.37
LSD 0.05	0.1	0.5	0.1	0.11	0.55	0.11	0.13	1.1	1.1	0.9	1.28	1.51	1.40	0.47
LSD 0.01	0.2	0.7	0.2	0.16	0.74	0.15	0.17	1.4	1.4	1.2	1.74	2.05	1.91	0.61

). The data presented clearly demonstrate that the genotypes had substantial variations in grain production per plant, which were statistically significant (Table 7). The Sakha 108, Giza 179, and Giza 178 varieties exhibited the highest average weights per plant, measuring 45.6, 45.5, and 45.4 g, respectively. The genotype Sakha 107 had the lowest yield, measuring 42.3 g per plant. The grain production per plant varied between 42.3 and 45.6 g per plant across the different genotypes. The analysis of variations was conducted on the combined data from six different environments, comprising three sowing dates and two seasons. Table 6 and Table 7 demonstrate notable variations in features across different contexts. The genotypic variation (s2g) was greater than the error variation (s2e) for the number of filled grains, seed set percentage, 1000-grain weight, and grain yield per plant. This indicates that environmental effects have little influence on these traits. Conversely, the s2g was significantly larger than s2 for all traits, indicating that genetic effects play a crucial role in the expression of these traits (Table 7).

The acquired results are consistent with those reported by [4]. The phenotypic (PCV) and genotypic (GCV) coefficients of variability were estimated and found to have modest variances for all traits evaluated across the six environments (

Table 8. AMMI-based stability 12 indexes (low numbers indicate stability) and superiority index based on stability and yield (high numbers indicate superiority).

GEN	GY	ASTAB	ASI	ASV	AVAMGE	DA	DZ	EV	FA	MASI	MASV	SIPC	ZA	SI
S101	45.6	0.64	0.36	1.37	2.63	1.48	0.44	0.07	2.18	0.36	1.6	1.26	0.26	77.57
S104	44.9	0.85	0.13	0.48	2.85	1.34	0.64	0.14	1.8	0.17	1.04	1.26	0.17	78.82
S105	44.4	0.11	0.06	0.24	1.2	0.49	0.22	0.02	0.24	0.07	0.43	0.51	0.07	82.82
S106	44.4	0.36	0.23	0.88	1.96	1.06	0.34	0.04	1.13	0.23	1.18	0.87	0.18	69.20
S107	42.3	1.18	0.32	1.23	3.35	1.77	0.69	0.16	3.13	0.34	1.74	1.86	0.31	23.47
S108	45.3	0.35	0.15	0.59	2.16	0.97	0.37	0.04	0.95	0.16	1.07	0.87	0.14	89.28
G177	43.8	0.62	0.12	0.47	2.21	1.18	0.53	0.09	1.39	0.15	1.08	1.13	0.14	68.08
G178	45.4	0.11	0.17	0.66	1.46	0.64	0.17	0.01	0.4	0.17	0.68	0.45	0.1	92.19
G179	45.5	0.95	0.53	2.05	4.4	1.93	0.49	0.08	3.73	0.53	2.05	1.15	0.29	71.91
G182	43.4	2.26	0.73	2.83	6.6	2.87	0.8	0.21	8.24	0.73	3.11	2.14	0.48	17.20
E.Y.	43.9	0.92	0.43	1.66	3.1	1.78	0.53	0.09	3.16	0.43	1.93	1.45	0.31	45.92
SR300	43.5	0.4	0.18	0.7	2.17	1.05	0.39	0.05	1.1	0.19	1.14	1	0.17	57.95

AMMI-based stability parameter (ASTAB), AMMI stability index (ASI), AMMI stability value (ASV), sum across environments of the absolute value of GEIs modeled by AMMI (AVAMGE), Annicchiarico’s D parameter values (DA), Zhang’s D parameter (DZ), sums of the averages of the squared eigenvector values (EV), stability measure based on fitted AMMI model (FA), modified AMMI stability index (MASI), modified AMMI stability value (MASV), sums of the absolute value of the IPC scores (SIPC), absolute value of the relative contribution of IPCs to the interaction (Za), superiority index (SI).

). It is noteworthy that certain rice cultivars consistently exhibited similar characteristics across various conditions, including the number of panicles per plant, panicle weight, the number of filled grains per panicle, seed set percentage, 1000-grain weight, and grain yield per plant. However, certain types of rice varied in performance across different environments. For example, the rice cultivar Sakha 108 showed superiority in terms of the number of filled grains per panicle, 1000-grain weight, and grain weight per panicle in environment E1. On the other hand, Giza 178 performed better in terms of the number of filled grains per panicle, while Sakha 101 had a higher seed set percentage.

3.3. Stability Parameters

The predictability of genotypes for the number of filled grains per panicle varied from 0.04 for Sakha 106 to 1.37 for Sakha super 300 (

Table 9. Mean performance and phenotypic stability measurements of the number of filled grains/panicle and seed set (%) for the twelve genotypes in three environments.

Genotypes	Number of Filled/Grains Panicle		Seed Set (%)		1000-Grain Weight (g)		Grain Yield/Plant (g)
	bi	S2di	bi	S2di	bi	S2di	bi	S2di
Sakha 101	1.08	1.863	0.94	0.106	0.59	0.002	1.21	0.315
Sakha 104	0.85	1.525	0.80	0.343	1.09	0.027	1.20	0.276
Sakha105	1.23	1.601	1.12	0.273	−0.96	0.127	1.11	0.027
Sakha 106	0.04	4.594 **	0.13	1.908	−0.85	0.046	0.96	0.395
Sakha 107	0.96	5.946 **	0.77	0.045	0.64	0.002	0.98	0.788
Sakha 108	1.13	1.971	1.04	1.085	0.48	0.004	1.00	0.244
Giza177	1.33	0.249	1.47	0.480	0.76	0.014	0.88	0.270
Giza178	0.67	144.307 **	2.06	1.439	0.34	0.006	0.92	0.105
Giza179	1.04	8.322 **	0.57	1.256	0.83	0.016	0.72	0.514
Giza 182	1.05	12.466 **	0.50	1.089	7.77	0.509	1.40	1.166
Egyptian Yasmine	1.27	16.268 **	1.22	0.242	0.84	0.000	0.75	0.484
Sakha super 300	1.37	5.268 **	1.39	3.204 *	0.48	0.003	0.86	0.176

Symbols * and ** indicate p < 0.05, p < 0.01, respectively.

). The examined genotypes were classified into three categories based on the regression coefficient (b_i) values. The initial group consisted of the most genetically stable genotypes, specifically Sakha 107, which exhibited the coefficient of regression b_i values equal to 1. The second group consisted of genotypes that were better suited to unfavorable environments, namely Sakha 104, Sakha 106, and Giza 178 (with the lowest b_i values). On the other hand, genotypes Sakha 101, Sakha 105, Sakha 108, Giza 177, Giza 179, Giza 182, E. Yasmin, and Sakha super 300 were sensitive to input and adapted to high-potential environments. Regarding the grain yield, the rice cultivars Sakha 106, Sakha 107, Sakha 108, and Giza 178 had a bi value close to 1 (Table 9), suggesting that these cultivars were well suited to varied environmental conditions. Genotypes Giza 177, E. Yasmin, and Sakha super 300 exhibited b_i values below 1, showing their suitability for unfavorable conditions. Sakha 101, Sakha 104, Sakha 105, and Giza 182 displayed a b_i value greater than 1.

3.4. Environmental Indices

An environment index is a quantitative metric that assesses the appropriateness of a specific environment in a given area. By analyzing these indicators, we may determine the optimal conditions required to fully use the genetic potential. If an environment has a positive index value, it can be inferred that it is a favorable environment for genotypes. The indices of E1 and E4 were 2.08 and 2.07, respectively. The value for E4 was 0.8, while the values for E3, E5, and E6 were −0.37, −1.25, and −3.34, correspondingly.

Figure 2 displays the results of a GGE biplot analysis, which is widely regarded as the most efficient method for summarizing both genotypes and the interaction between genotypes and the environment. This analysis helps select the most suitable genotype for each environment and evaluate the stability of the tested genotypes. Figure 2a demonstrates that genotypes located closer to the origin exhibit greater stability compared to those located further away from the origin. Furthermore, the environments were grouped into distinct clusters, namely (E1, E2, E3, and E4) and (E5 and E6). Furthermore, the genotype characterized by a long-dashed line has a higher GY (grain yield) compared to the genotype with a short-dashed line. The genotypes positioned on the vertices of the polygon in Figure 2b are the most responsive to different settings, meaning that they show greater variability in their performance. On the other hand, the genotypes inside the polygon are more stable across environments, indicating that they are less influenced by changes in environmental conditions. The lines of equality partition the graph (Figure 2b) into five sectors. A total of six habitats were preserved, divided into two sectors, and further categorized into two mega-environments.

In the first mega-environment, which included E1, E2, E3, and E4, the top-performing varieties were Sakha 101, Sakha 104, Sakha 105, and Sakha 108. Conversely, in the second mega-environment encompassing E5 and E6, the leading varieties were Sakha 106, Giza 178, and Giza 179. Figure 2c,d depict the respective rankings of environments and genotypes.

The Bayesian network in Figure 3 depicts the impact of treatments on traits, as well as the influence of traits on other traits, using the Bayesian information criterion (BIC). A greater negative BIC value implies a more pronounced effect. The data showed that genotypes had a significant and direct impact on plant height (PH), days to heading (DTH), grain weight (PW), thousand-kernel weight (TKW), the total grain number per plant (TGNP), and the filled grain number per plant (FGNP). Additionally, the influence of genotypes on FGNP was mediated indirectly through its direct effect on TGNP. Conversely, the environment had a significant direct impact on PL (primary productivity) and PNPP (primary net primary productivity), as well as an indirect impact on PL through PNPP. Additionally, the environment had an indirect impact on FGNP (final gross net primary productivity) and PW (plant water) through PNPP, as well as an indirect impact on GY (growth yield) through PL. PW was influenced by both genotypes and PNPP, but GY was significantly impacted by HI. To summarize, genetics exerted a greater influence on traits compared to the environment.

3.5. Stability and Superiority Indexes

The total of twelve stability indexes, along with one superiority index, is displayed (Table 9). Based on the stability indexes, Sakha 105 exhibited the highest level of stability among the genotypes, whereas Giza 182 displayed the lowest level of stability. The stability indices primarily consider the stability of genotypes in different conditions while disregarding the yield potential of the genotype. In contrast, the superiority index takes both the stability and prospective yield of the genotype into consideration. According to the superiority index, Giza 178 exhibited the highest level of superiority among the genotypes, whilst Giza 182 displayed the lowest level of superiority.

4. Discussion

Multi-environmental trials are crucial in plant breeding projects since they assess crop yield and adaptability across several settings. To enhance the accuracy of genotype selection, one can assess the genotype-by-environment interaction using mega-environment analysis, genotype evaluation, and multi-environment stability analysis [27]. Genetic diversity can be achieved by breeding inter-subspecific cultivars or accessions having distinct genetic lineages [28,29,30]. Environmental variability can be obtained by increasing the magnitudes of location factors such as geographic distance, variances in climatic regions (temperate versus tropical), and differences in elevation [31,32].

Our breeding operations primarily prioritize the preservation of the commercial worth of cultivars. Crop enhancements are primarily achieved by the interbreeding of superior cultivars within restricted genetic pools, resulting in limited genetic diversity among the cultivars [33,34,35]. Furthermore, the test environments were located in geographically proximate agriculturally advantageous flat regions. The stability of a cultivar/genotype pertains to the degree of variation in the phenotypic manifestation of crop performance across diverse environmental conditions. The aim of this study was to assess the idea of “dynamic” stability, which refers to the performance of a stable genotype in different environments. Consequently, not all genotypes exhibited the same response to the environment [36,37].

The results of our study indicate that the genotypes we tested in various locations over two consecutive years are extremely susceptible to climatic zones and environmental factors [38]. The variations seen may be attributed to the disparities in terrain and climatic circumstances among various sites and years where the experiments were carried out. Achieving effective breeding outcomes for yield and related traits in rice necessitates precise measurement and analysis of genotype, environment, years, and their interactions within the breeding protocol [39]. The current discovery of noteworthy sources of variation has been previously documented in rice and other cereal crops [10,27]. Significant effects of both linear and nonlinear components of G × E interaction were seen for all examined traits, as evidenced by highly significant mean squares resulting from G × E (linear) interaction and pooled deviation. The reported results are consistent with those obtained by [40]. The grain production of different rice varieties exhibited significant variations in response to changes in climatic variables, such as planting dates and years [41].

Hence, the cultivar that exhibits consistent and predictable crop output is suitable for implementing risk management strategies to mitigate yield losses in challenging lowland environments [42]. A statistically significant difference has been observed across cultivars and settings in a linear manner, as opposed to a nonlinear manner, in which the deviations are aggregated. This finding highlights the significance of the interplay between genotype and environment in determining grain yield performance and its components across several locations [28,43]. Therefore, the phenotypic stability assessment took into account both the linear (b_i) and nonlinear (S²d_i) components of G × E interactions [23,44]. Moreover, it has been documented that a high average value of the linear regression coefficient, equivalent to the nonlinear coefficient, is advised for achieving a diverse range. Furthermore, it has been suggested that nonlinear regression can be employed to assess stability, whereas linear regression can be utilized to evaluate the impact of different environmental circumstances on varietal response [45]. Hence, it is crucial to take into account the average and the deviation from the regression of each cultivar when assessing the stability and linear regression for evaluating the responsiveness of the variety. Stresses can limit the growth resources available to plants, resulting in a reduction in the size of plant organs such as leaves, tillers, and spikes [41].

The stability regression coefficient (b_i) and deviation from regression (S²d_i) for the genotypes under study are shown in Table 7. A stable genotype refers to a genotype that exhibits a consistently high level of performance, with a regression coefficient (b_i) equal to 1 and a deviation from regression of 0. This indicates that the system is well suited for favorable situations, but its performance will be subpar in unfavorable environments. Therefore, these rice cultivars possess significant value for conducting investigations in harsh environments, as the identified types have been suggested for their suitability in unfavorable environmental conditions, as indicated in Table 9 [2].

In the context of selecting and developing superior rice cultivars, genetic studies on yield and yield quality serve as crucial benchmarks for both farmers and industry professionals. As highlighted in previous research, evaluations incorporating various stability statistics offer a more comprehensive understanding of genotype performance across different environments. For instance, a study on sweet potato genotypes in West Java, Indonesia, demonstrated the importance of genotype-by-environment interactions (GEIs) in identifying high-yield and stable genotypes. This research utilized a randomized block design across three environments and found significant effects of genotypes, environments, and GEIs on yield and quality traits. Specifically, GEIs contributed substantially to the variability observed in yield, sweetness, moisture content, tuber diameter, and tuber length. We identified five genotypes (G4, G6, G7, G31, and G32) as both high-yielding and stable across environments. Furthermore, genotype-by-yield ×trait (GYT) analysis revealed seven superior genotypes based on yield and quality, namely G7, G15, G4, G20, G6, G31, and G14. Notably, G4, G6, G7, and G31 were recognized for their high and stable yields combined with favorable quality traits. These findings underscore the value of integrating multiple stability measures and GYT analysis to identify and release superior varieties [46].

The results of pooled ANOVA for stability according to Eberhart and Russell [29], as well as the AMMI model, revealed that the variance of genotypes was significant for rice grain yield, indicating that the performance of rice genotypes was different; the genotypes also had differential responses to the changes in years and locations (environments).

In this study, the genotype Sakha 105 had regression coefficients (b_i) of 1.11 and was observed to be stable. According to the Eberhart and Russell [29] model, a slope of >1.0 with high mean yield and nonsignificant squared deviation (S²d_i) is suitable for a favorable environment.

The AMMI model defines the genotype–environment interaction. The AMMI model revealed that a major part of the variation in rice grain yield is explained by the environment, which indicates environmental diversity. The AMMI1 biplot analysis showed that the variation was due to the main effect of rice grain yield and the interaction effect. Genotypes with IPCA1 scores close to zero were characterized as having low interaction effects and consequently considered stable.

5. Conclusions

This study highlights that environmental factors are the primary drivers of variations in rice yield and related traits among Egyptian rice genotypes, with genotype–environment interactions playing a secondary role. Certain genotypes, like Sakha 104 and Giza 178, demonstrated stable performance across different environments, while others were better suited to specific conditions. Notably, Sakha 106, Sakha 107, Sakha 108, and Giza 178 showed b_i values close to 1, indicating they are well adapted across diverse environments. In contrast, Giza 177, E. Yasmin, and Sakha Super 300, with b_i values below 1, were more suitable for poor environments, while Sakha 101, Sakha 104, Sakha 105, and Giza 182, with bi values above 1, were better suited to favorable environments but may underperform in less favorable conditions. These findings underscore the importance of multi-environmental trials for selecting stable, high-yielding rice cultivars. The results provide a foundation for targeted breeding programs aimed at optimizing rice productivity by focusing on genotypes that consistently perform well across various growing conditions, considering both stability and adaptability.