Next Article in Journal
Extracts of Fruits and Plants Cultivated In Vitro of Aristotelia chilensis (Mol.) Stuntz Show Inhibitory Activity of Aldose Reductase and Pancreatic Alpha-Amylase Enzymes
Next Article in Special Issue
Magnetic Iron–Improved Growth, Leaf Chemical Content, Yield, and Fruit Quality of Chinese Mandarin Trees Grown under Soil Salinity Stress
Previous Article in Journal
Biogenic CuO and ZnO Nanoparticles as Nanofertilizers for Sustainable Growth of Amaranthus hybridus
Previous Article in Special Issue
Proteomic Changes in Paspalum fasciculatum Leaves Exposed to Cd Stress
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Detection of Superior Rice Genotypes and Yield Stability under Different Nitrogen Levels Using AMMI Model and Stability Statistics

1
Rice Research and Training Center, Field Crops Research Institute, Agricultural Research Center, Kafrelsheikh 33717, Egypt
2
Department of Biology, College of Science, Princess Nourah Bint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia
3
Genetics Department, Faculty of Agriculture, Tanta University, Tanta 31527, Egypt
4
Center of Excellence in Biotechnology Research, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
5
Genetics Department, Faculty of Agriculture, Kafrelsheikh University, Kafrelsheikh 33716, Egypt
*
Author to whom correspondence should be addressed.
Plants 2022, 11(20), 2775; https://doi.org/10.3390/plants11202775
Submission received: 25 September 2022 / Revised: 13 October 2022 / Accepted: 15 October 2022 / Published: 19 October 2022

Abstract

:
Sustainable agriculture is a prerequisite for food and environmental security. Chemical fertilization, especially nitrogenous fertilization, is considered the most consumed for field crops. In rice crops, plants consume much less than half of the applied N-fertilizer. In the current investigation, multiple N environments were generated by applying different N doses of urea fertilizer to a permanent transplanted field for two successive summer growing seasons at the rice research and training center, Kafrelsheikh, Egypt. A set of 55 genotypes consisting of 25 Jabonica, 4 Tropical Japonica, 20 Indica, and 6 Indica/Japonica were transplanted under no N (0N), Low N (LN), medium N (MN), and High N (HN) (i.e., 0, 48, 96, and 165 Kg N ha−1, respectively). Highly significant differences were detected among the tested genotypes. AMMI analysis of variance revealed the existence of the genotype via nitrogen interaction (GNI) on yield performance. The GNI principal components (IPCA); IPCA1 and IPCA2 scores were significant and contributed values of 71.1 and 21.7%, respectively. The highest-ranked genotypes were MTU1010, IR22, SK2046, SK2058, IR66, and Yabani LuLu based on their grain yield means (30.7, 29.9, 29.5, 29.3, 28.8, and 28.5 g plant−1). These genotypes were grouped into the same subcluster (SCL) according to the stability analysis ranking matrix. Based on AMMI analysis and biplots, MTU1010 and Yabani LuLu showed yield stability across environments. Meanwhile, the which-won-where biplot showed that IR22 was superior under unfavorable N-levels and MTU1010 was stable across the different environments. These findings are considered to be of great importance to breeders for initiating low-nitrogen-input breeding programs for sustainable agriculture.

1. Introduction

Rice (Oryza sativa L.) is a staple food in many countries and a major component of diets in many others. In many parts of the world, a lack of rice supplies could lead to starvation. Globally, 755 million tons of paddy rice are produced from 162 million ha of land [1]. To meet the needs of the world’s rapidly growing population, increasing grain yield per unit area of rice is critical to combating poverty. This goal could be met by cultivating rice cultivars with high yield potential and using appropriate management [2].
Nitrogen (N) is the most important nutrient that is needed by plants for growth and completing their life cycle. However, the rice crop consumes much less than half of the applied N fertilizer. The recovery efficiency of nitrogen fertilizer in rice amounted to 30–39%, and the nitrogen use efficiency decreases with the increase in the applied nitrogen [3]. N application losses contribute to soil deterioration, groundwater pollution, and emissions of ammonia and greenhouse gases [4,5]. Accordingly, for rice production sustainability that is less harmful to the environment, it is necessary to identify the superior genotypes under low application of N-fertilizer. Hence, rice breeders are always evaluating a large set of genotypes under different environments to identify the high-yielding genotypes that have adapted to different environments [6,7].
Grain yield (GY) is the economic determinant of the best-performing genotypes [8]. GY is a quantitative trait determined by the additive main effect of environment (E) and genotype (G) in addition to the nonadditive effect of the G X E interaction (GEI) [9]. Breeders focus on the GEI effect to identify the yield stability of genotypes across different conditions and environments, which cannot be revealed by the separate effects of genotype or environment [10,11,12,13,14]. The GY heritability is exposed to variability across different environments [15,16], which hinders the accuracy of superior varietal selection processes [17]. Therefore, widely adapted genotypes with the ability to produce stable high yields across diversified environments constitute a major goal for rice breeders. Hence, it is critical to evaluate the adaption of several genotypes across different N inputs to identify superior genotypes with stable GY under different N environments to detect GNI. Furthermore, the analysis of yield stability under different soil conditions is found to be associated with imbalanced yield stability [18]. Several statistical analysis methods have been reported and developed to manipulate GEI including parametric and non-parametric stability statistical methods. Parametric methods include univariate and multivariate methods. Wricke’s ecovalence (Wi2) [19], Shukla’s stability variance (s2) [20], the coefficient of variance (CV) [21], Environmental variance (S2) [22], the Mean-variance component (q) [23], the GE variance component (q’) [24], the Regression coefficient (bi) [25], and many others constitute univariate methods. Multivariate methods include the additive main effects and multiplicative interaction (AMMI) model [26] and the GGE biplot method [27]. AMMI analysis considers both ANOVA and principal component analysis (PCA). The results of AMMI provide genotypes’ yield stability under different environments and facilitate the precise selection of the best-performing genotypes for the environment under study [28,29,30]. Multivariate methods estimate GEI by following approaches such as the ‘which-won-where’ pattern, identifying mega environments and superior genotypes across different environments under study, and ranking environments [31]. Non-parametric methods include Nassa and Huhn’s and Huhn’s statistics (S) [32], Kang’s rank-sum (KR) [33], TOP-Fox (TOP) [34], Thennarasu’s non-parametric statistics (NP) [35], and the Genotype stability index (GSI) [36]. Recently, a new R package has been reported to manipulate the multi-environment trials (MET), called metan [37]. The package is a work-flow based approach that has a collection of functions to compute the most used parametric and nonparametric stability statistics. Unlike other R software packages that could possibly be used for analyzing MET data, metan is specifically coded for a complete analysis of MET trials (checking, manipulating, analyzing, and visualizing the data).
In recent years, there has been substantial interest in detecting both the yield and yield stability under different fertilizer conditions in different crops [38]. In our recent investigation, we estimated the yield stability of newly developed rice lines across different water treatments [15]. The present study aims to identify superior rice genotypes under low nitrogen input with stable yield performance over different N levels. The study is expected to enrich our knowledge regarding low-nitrogen-input genotypes with high yielding ability.

2. Results

2.1. Genotypic Variability of GY under Different N Environments

The mean performance of the evaluated genotypes under different N environments is presented in Supplementary Table S1. The data indicate that the genotypes exhibited different behavior as a response to the different N treatments, which was confirmed by the analysis of variance as indicated from the genotypic mean square value under different N levels (Supplementary Table S2). The mean performance of the tested genotypes across the different N environments showed that MTU1010 has the highest GY mean value across N environments, with 30.7 g plant−1. Moreover, environments’ mean performances showed an increase in the overall GY mean performance of the genotypes with each incremental dose of N, which ranged between 19.14 and 28.84 g plant−1 for the 0N and HN levels under study.

2.2. Combined AMMI Analysis of Variance of the GY and the Decomposition of GEI Effect

The AMMI model analysis for the genotypes’ GY revealed significant differences for both the main (G and E) and GEI effects (Supplementary Table S3). These findings show that there is considerable variability among the genotypes, environments, and their interactions. The G factor accounted for 49.97% of the total variance in the AMMI analysis, followed by E (33.69%) and GEI (7.40%). The decomposition of the GEI effect by AMMI analysis yielded three significant IPCAs. These three IPCAs explained 100 percent of the total GEI effects explained, which were all significant accounting for 71.1%, 21.7%, and 7.2% for IPCA1, IPCA 2, and IPCA3, respectively.

2.3. Environmental Effect on the Performance of the Genotypes

Based on the output of the AMMI analysis, Table 1 presents the mean genotypic GY performance in each N-environment, EPCA1, and the top five ranked genotypes for each N environment. The environments had different records of the overall genotypes’ GY mean. According to the environmental index the environments were classified into unfavorable and favorable for the grown genotypes under these environments. Environments with negative index estimates were unfavorable while MN and HN were favorable. 0N was the most unfavorable while HN was the most-favorable condition for the genotypes under study. Furthermore, according to the estimated EPC1, the results showed that N-environments contributed differently to the genotype stability for the GY. Among the different N- levels, 0N has the lowest PC1 estimate while HN has the highest PC1 record. These records indicate that the HN was the main contributor to the genotypic GY stability and 0N has the lowest contribution to the GEI.

2.4. Graphical Representation of Genotypes and N-Environments in the AMMI Biplots

The biplots for AMMI1-Means vs. PC1 and AMMI2-PC1 vs. PC2 were generated to identify the mega environment and stable genotypes across N treatments (Figure 1A, Supplementary Table S4). The AMMI1 biplot was generated to plot genotypes’ and environments’ means against their IPCA1. Subsequently, the environments 0N and LN were identified as low-yielding ones (Figure 1A AMMI1 Biplot). The AMMI2 biplot reveals the environmental and genotypical scores corresponding to both the GY IPCA1 and IPCA2 (Figure 1B AMMI2 Biplot). Based on this plot, the environments and/or the genotypes that are located in close proximity to the origin are less influenced by the effect of GEI. On the contrary, those genotypes and/or environments that are distanced from the origin either on the negative or positive side are more influenced by the GEI effect. Among the different N environments, LN was the nearest to the GEI origin compared with all other environments. Meanwhile, on the genotypes side, C22, WAB 450-1-B-P-91-HB, Black Rice, Arabi, Yun Lu No. 48, GZ 6522-15-1-1-3, Giza178, SKC 23808-28-5-2-1-1, Taikeng Yu 1420, IR28, Sakha104, and SK2058 generated a polygon where they were the most interactive genotypes with the changes in N levels.

2.5. Which-Won-Where GGE Biplot Analysis

The sum of the first and second PCA axes explained 96.36% of the total GEI variation, as can be seen in the GGE biplot presented in Figure 2A. The genotypes that were farthest away from the biplot origin include Black Rice, Yun Lu No. 48, Arabi, IR22, MTU1010, Giza178, Reiho, SKC 23808-28-5-2-1-1, and Giza177. These genotypes, called vertex genotypes, have the longest vectors with respect to their direction. Vertex genotypes are the most responsive genotypes to the environment in their direction. The vertex genotypes with no environmental indicators nearby are the poorly performing ones. Accordingly, Black Rice had the poorest performance in all environments under study. The which-won-where GGE biplot of GY divided the four N treatments into two sectors. The low treatments, 0N and LN, were located together in one sector where the IR22 genotype was the best-performing genotype, while MN and HN were located together in the other sector. MTU1010 showed a stable performance under different N treatments, placing it in the top five ranked genotypes across the tested N environments.

2.6. Superior Genotypes Selection Based on GY Means and Stability Parameters

The ranking based on the various stability statistics is presented in Supplementary Table S5. Based on GY, which considers the main selection criteria for genotype selection, the genotypes MTU1010, IR22, SK2046, SK2058, and IR66 were the highest-ranked genotypes in this regard. From another stability statistics analysis, the weighted average of absolute scores (WAASB) was used to better identify the best genotypes based on the mean GY and stability. The biplot shows the distribution of the tested rice genotypes and environments based on the genotypes’ GY mean and WAASB values, as shown in Figure 2B.
The first quadrant, I, contains the genotypes Black Rice, WAB 450-1-B-P-91-HB, Yun Lu No. 48, Sabieny, Pusa Basmati 1, Giza177, SKC 23808-28-5-2-1-1, GZ 6522-15-1-1-3, GZ5830-59-10-2, Taikeng Yu 1420, and IR 67075-2B-5-2. These genotypes showed lower grain yield compared to the mean grain yield. In addition, this quadrant has low GY environments 0N and LN. Accordingly, the genotypes and environments located in this quadrant have the largest response to GEI. The second quadrant contains the environments with GY above average, MN and HN environments, and the genotypes that have a GY above average with a high GEI response, namely, Arabi, Reiho, Giza178, SK2034, WAB 880 SG 73, Giza14, Nahda, and IR22. The genotypes GZ7718-13-3-1-3, SKC23822-304-3-1-1-1, GZ6214-4-1-1-1, GZ7718-13-3-2-2, IR68353-35-3-3-2-2-1-2, IR68373-R-R-B-22-2-2, Yen Geng 135, Agami M.1, IR 73689-31-1, Sakha103, Giza182, IR70, and GZ6903-3-4-2-1 showed low yield but stable performance across the different environments. These genotypes were located in the third quadrant and there were no environments in this one. The fourth quadrant contains the high-yielding stable genotypes. Those genotypes are MTU1010, SK2046, SK2058, GZ7922-B-44-1, Egyptian Yasmin, Yabani LuLu, IET 1444C22, SK2035, BG 304, IR 74, Milyang63, GZ6910-28-1-3-1, IR74, IR64, GZ6903-3-4-2-1, GZ6910-28-1-3-1, E 7034, and Nabatat Asmar.

2.7. Cluster Analysis and Dendrogram Based on the Stability Statistics Values

A hierarchical cluster analysis was conducted based on the squared Euclidean distance via Ward’s method using values of the stability statistics for each genotype. The grouping pattern resulting from this analysis revealed the distribution of the tested rice genotypes into two main clusters (CL) (Figure 3). CL-1 was further divided into two subclusters (SCL). SCL1 contains the genotypes WAB 450-1-B-P-91-HB, Black Rice, Arabi, and Yun Lu No. 48, while SCL2 has the genotypes Sabieny, IR 67075-2B-5-2, GZ7922-B-44-1, Sakha103, Giza175, GZ6214-4-1-1-1, IR7421-35-1-1-2, GZ7718-13-3-2-2, IR74, Agami M.1, and SKC23822-304-3-1-1-1. At the same time, CL2 was divided into two main SCLs, SCL3 and SCL4. SCL3 was divided into two groups, SCL3-1 and SCL3-2. SCL3-1 contains the top six ranking genotypes, MTU1010, IR22, SK2046, SK2058, IR66, and Yabani LuLu, with the highest GY mean across environments. These genotypes are indica, indicating that indica-type genotypes have more adaptability to withstand different N doses and provide high yield ability. Meanwhile, the genotypes Taikeng Yu 1420, Giza178,48, BG 304, IET1444, IR70, SK2034, C22, SK2035, Giza14, and Egyptian Yasmin were clustered in SCL3-2. Similarly, SCL4 has two main clusters, SCL4-1 and SCL4-2. The genotypes Pusa Basmati 1, Milyang63, Nabatat Asmar, GZ6903-3-4-2-1, Giza159, and Nahda clustered in one group together with the other group of IR64, Sakha101, GZ6910-28-1-3-1, E7034, Sakha104, WAB880 SG 73, Yen Geng 135, IR68373-R-R-B-22-2-2, IR73689-31-1, Giza182, and IR28. Both groups were clustered from SCL4-1. Meanwhile, the genotypes GZ5830-59-10-2, IR 68353-35-3-3-2-2-1-2, Giza177, GZ 6522-15-1-1-3, and SKC 23808-28-5-2-1-1 were grouped in SCL4-2.

3. Discussion

GY is the main criterion for a breeder based on which new genotypes are selected. Basically, it is the final product of the genotype performance, which considers the output of GEI. In recent years, great attention has been paid to sustainable cultivation. It mainly focuses on reaching the maximum yield potential of a genotype with a lower input of water, chemical fertilizers, and pesticides [4,15].
In the current investigation, considerable variability was detected among the tested genotypes based on their GY performance under different N levels. The genotypes have different response patterns for each incremental dose of N fertilizer (Supplementary Figure S1 and Supplementary Table S1). The AMMI-ANOVA analysis for the combined two-season data revealed highly significant differences between genotypes, N-environments, and GEI. However, the proportion for the genotypic variability was the highest. This finding explains the existence of selection capacity to find better-performing genotypes under different N levels since all other elements were consistent except for the application of N- chemical fertilizer. This suggestion is further confirmed by our findings of an interaction between the genotypes and the different N levels (GEI) regarding the amount of variability. However, several investigations confirmed the findings of significant variability due to G, E, and GEI, but the environments had the highest source of variation (i.e., [39,40]). These trials were conducted at different locations, which means the environmental conditions are entirely different from each other. In our investigation, the environments 0N and HN were the lowest and highest GY mean records, indicating the possibility of further enhancing the varietal productivity under low-N-input environments. These findings confirm the need for stability analysis for the tested genotypes under multi-N-levels.
Among the widely used stability analysis methods, the AMMI model combines both ANOVA and GEI to measure the magnitude of genotype and environment variability [26,41]. In our investigation, a significant interaction between N fertilizer and the tested genotypes in GY, as a high proportion of the first two IPCAs, was found. Variability in response to GEI identifies that genotypes with different ranks of GY potential correspond to specific environments and their stability across the environments [42]. For the tested environments, high N levels were favored by the genotypes. However, genotypes that efficiently developed high yields while growing under poor N levels have the ability to tolerate N deficiency in the soil. Genotypes IR22, Giza14, C22, MTU1010, and SK2058 were the highest yielding under the 0N environment. Interestingly, those genotypes are the indica type, except Giza14, which is the japonica type. Furthermore, these genotypes were comparable to the newly developed cultivars under this level. The AMMI biplots were used to measure the distribution of G, E, and GEI between them [43], as indicated in the AMMI1 and AMMI 2 biplots shown in Figure 1A,B. The AMMI1 biplot explained the distribution of the tested genotypes based on the GY mean and IPCA1 values, which are useful to identify the genotypes MTU1010, IR66, and Yabani LuLu with the lowest IPCA1 scores and a high GY mean. Furthermore, the tested environments exhibited wide dispersion relative to the GY performance of the genotypes under study. The AMMI2 biplot further described the multiplicative effects of GEI through the first two IPCAs (Figure 1B). Accordingly, the genotypes G25, MTU1010, BG 304, and Yabani LuLu were the best genotypes having the lowest IPCAs values and located in close proximity to the origin. Using AMMI biplots, several investigations were able to identify their corresponding genotypes such as in barley [44], rice [40], Sorghum [45], cassava [46], and maize [47]. MTU1010 and Yabani LuLu are among the two groups most productive and stable. Intriguingly, IR22 showed specific adaptability to the low N environments (0N, LN, and MN), which was further confirmed by the which-won-where GGE biplot, as shown in Figure 2A. Moreover, MTU1010 was located at the boundary between the two sectors where the four environments were scattered, indicating that this genotype has stable performance across the different N-environments. Both IR22 and MTU1010 are indica genotypes, while the study has other types such as japonica and indica/japonica types (Table 2). IR22 is the cultivar derived from the miracle rice, IR8. MTU1010 is an elite, mega Indian cultivar grown extensively in India [48]. These findings may indicate that indica type has better performance than other types under study environments, which may be due to high grain nitrogen use efficiency [49]. This point requires further analysis and investigation.
We further estimated the stability measure of the weighted average of the absolute score (WAASB), which was taken into consideration as the sum of absolute values of the IPCAs [50]. To identify the highly adaptive and best-performing genotypes, we generated a biplot to simultaneously consider the GY and stability measure of WAASB. Among the tested environments, the two MN and HN environments were productive with average GY means higher than the grand mean. The genotypes MTU1010, IR66, Yabani LuLu, and GZ7922-B-44-1, with high GY and the lowest WAASB values, were identified as the most stable genotypes across the different N environments under study. The WAASB index was proven to be an efficient index in the selection of stable, high GY barley genotypes across different environments [51,52].
To identify groups of the genotypes based on the different stability statistical analysis ranking matrices, CLA was conducted. The most notable findings were that the low-yielding genotypes WAB 450-1-B-P-91-HB, Black Rice, and Yun Lu No. 48 were grouped in a single SCL while the high-yielding ones, IR22, SK2046, SK2058, MTU1010, Yabani LuLu, and IR66, were grouped in another SCL. The former group was low yielding under favorable and unfavorable N environments, while the latter group of genotypes was high yielding under all environments.
Considering our results, MTU1010 and Yabani LuLu showed high performance and GY stability across the different N-levels, while IR22 showed superior performance under the unfavorable N-levels. Our results need more investigations regarding the performance of these genotypes with respect to nitrogen use efficiency-related parameters.

4. Materials and Methods

4.1. Plant Materials

To maintain diversity, a set of 55 rice genotypes consisting of 25 Jabonica, 4 Tropical Japonica, 20 Indica, and 6 Indica/Japonica were used in the current experiment (Table 2). Among these genotypes were 15 old Egyptian rice cultivars. These genotypes were collected from rice GeneBank unit of the rice research and training center (RRTC), Kafrelsheikh, Egypt.

4.2. Experimental Location and Soil Properties

The 55 rice genotypes were evaluated for two successive seasons at the Sakha Agricultural Research Station experimental farm (31.09° N and 30.9° E), located in Kafrelsheikh, Egypt. During the two seasons, the climate was typified by high temperatures and no rainfall. The soil analysis was conducted according to Piper [53] and Black et al., [54], where the soil texture of both seasons was clayey with a total N of 759 and 770 mg kg−1 for the two seasons, as shown in Supplementary Table S6.

4.3. Field Experiments and Environmental Conditions

These diverse genotypes were evaluated following a transplanting method under 4 different N-levels (Urea form) to generate the different environments of no N (0N), Low N (LN), medium N (MN), and High N (HN) (i.e., 0, 48, 96, and 165 Kg N ha−1, respectively). The experimental design of each environment was employed according to the randomized complete block design (RCBD) with three replicates. Each plot was five one-meter-long rows, with 20 cm × 20 cm spacing. N fertilizer was applied in two splits: 2/3 as a basal application and 1/3 a month after transplanting according to the N fertilizer amount for each N environment. Other cultural practices (such as field preparation, fertilizers, and weed management) were applied according to the standard to maintain the consistency of all factors except N-levels across the studied environments. At the harvest stage, the plots were harvested, and GY (g plant−1) was estimated for each experimental plot after adjusting the moisture content.

4.4. Statistical Analysis

All statistical analyses were conducted in the statistical software R (R Core Team, [55]) version 4.1.1. The “metan” package [37] was employed to conduct the analysis of variance, AMMI analysis of variance [43], genotype plus genotype by environment (GGE) biplot analysis [56], stability statistical analysis, and weighted average of absolute scores [50]. The “Nbclust” package was used to conduct hierarchical cluster analysis [57,58]. The AMMI analysis was conducted based on the following mathematical formula:
y i j N =   μ +   g i +   e j +   Σ λ k Y i k α j k +   ε i j
where yij is the yield of the ith genotype in the jth environment, N is the number of PCI in the AMMI model, μ is the overall mean of the genotypes, and gi and ej are the genotype and environment diversions from the overall mean. λk is the eigenvalue of the PCA axis k, Yikand αjk are the GE-PCs scores for axis k, and Σij is the remaining value. Meanwhile, the GGE model was considered by the following formula:
y i j N =   μ + e j +   Σ λ k Y i k α j k +   ε i j

5. Conclusions

In the current investigation, 55 rice genotypes were evaluated under 4 different N- levels to identify the genotypes with high yields under low N input and with stable performance under different N-levels. The results showed that the highest-ranked genotypes were IR22, SK2046, SK2058, MTU1010, YABANI LULU, and IR66 based on the GY mean. These genotypes were grouped in the same SCL according to the stability analysis ranking matrix. Based on AMMI analysis and biplots, MTU1010 and YABANI LULU showed stability, while the which-won-where biplot showed that IR22 was superior under unfavorable N-levels and MTU1010 was stable across the different environments. The study recommends incorporating these genotypes into the breeding for low-nitrogen fertilizer input programs for sustainable agriculture.

Supplementary Materials

The following supporting information can be downloaded at: www.mdpi.com/article/10.3390/plants11202775/s1, Figure S1: Heat map of the GY performance of the tested 55 genotypes across the 4 different environments; Table S1: Genotypes GY means (g plant-1) of the two successive seasons under the different N environments under study; Table S2: RCBD ANOVA table for the mean GY values of the two seasons for the tested 55 genotypes under the different 4 N levels grown under 3 replicated; Table S3: AMMI analysis of variance for GY (g plant-1) for the tested 55 genotypes under 4 different N levels; Table S4: GY means and IPCAs for the genotypes and the environments under study; Table S5: Ranking of the genotypes based on the different stability statistical analysis; Table S6: Physical and chemical properties of the experimental soil before rice transplanting at the two growing seasons.

Author Contributions

Conceptualization, M.A., A.-S.E.D., M.H.A., K.A. (Kotb Attia) and S.A.D.; methodology, M.A.; software, M.A., K.A. (Khadiga Alharbi) and N.A.A.-S.; validation, M.A., M.H.A., K.A. (Khadiga Alharbi), N.A.A.-S. and K.A. (Kotb Attia); formal analysis, M.A.; investigation, M.A., S.A.D., T.A.E.-M., S.M. and M.E.E.-D.; resources, T.A.E.-M.; data curation, M.A.; writing—original draft preparation, M.A.; writing—review and editing, M.A, E.-S.N., S.A.D., S.M., K.A. (Khadiga Alharbi) and N.A.A.-S.; visualization, M.A.; supervision, S.A.D. and A.-S.E.D.; project administration, M.A.; funding acquisition, M.A., K.A. (Khadiga Alharbi) and N.A.A.-S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are available in the Supplementary Materials.

Acknowledgments

The authors would like to thank Princess Nourah bint Abdulrahman University Researchers Supporting Project Number (PNURSP2022R188), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia. The authors extend their appreciation for the support and resources provided by the RRTC/FCRI, Egypt.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shaheen, S.M.; Antoniadis, V.; Shahid, M.; Yang, Y.; Abdelrahman, H.; Zhang, T.; Hassan, N.E.; Bibi, I.; Niazi, N.K.; Younis, S.A. Sustainable applications of rice feedstock in agro-environmental and construction sectors: A global perspective. Renew. Sustain. Energy Rev. 2022, 153, 111791. [Google Scholar] [CrossRef]
  2. Naeem, E.; Abd El-Megeed, T.; Emadeldin, Y.; Abushady, A.M.; Abdelrahman, M. Injected Anhydrous Ammonia Is More Effective Than Broadcast Urea as a Source of Nitrogen for Drill Seeded Rice. Agronomy 2022, 12, 942. [Google Scholar] [CrossRef]
  3. Fageria, N. Yield physiology of rice. J. Plant Nutr. 2007, 30, 843–879. [Google Scholar] [CrossRef]
  4. Abdel-Rahman, M.; El-Denary, M.; Ammar, M.; Abdelkhalik, A.; Draz, A.; Dora, S. QTL validation for grain yield and nitrogen use efficiency under different nitrogen levels in rice. Egypt. J. Genet. Cytol. 2015, 44, 235–251. [Google Scholar] [CrossRef] [Green Version]
  5. Ju, X.-T.; Xing, G.-X.; Chen, X.-P.; Zhang, S.-L.; Zhang, L.-J.; Liu, X.-J.; Cui, Z.-L.; Yin, B.; Christie, P.; Zhu, Z.-L. Reducing environmental risk by improving N management in intensive Chinese agricultural systems. Proc. Natl. Acad. Sci. USA 2009, 106, 3041–3046. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Krishnamurthy, S.; Sharma, P.; Sharma, D.; Ravikiran, K.; Singh, Y.; Mishra, V.; Burman, D.; Maji, B.; Mandal, S.; Sarangi, S. Identification of mega-environments and rice genotypes for general and specific adaptation to saline and alkaline stresses in India. Sci. Rep. 2017, 7, 7968. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Oladosu, Y.; Rafii, M.Y.; Abdullah, N.; Magaji, U.; Miah, G.; Hussin, G.; Ramli, A. Genotype× Environment interaction and stability analyses of yield and yield components of established and mutant rice genotypes tested in multiple locations in Malaysia. Acta Agric. Scand. Sect. B—Soil Plant Sci. 2017, 67, 590–606. [Google Scholar] [CrossRef]
  8. Mehana, M.; Abdelrahman, M.; Emadeldin, Y.; Rohila, J.S.; Karthikeyan, R. Impact of genetic improvements of rice on its water use and effects of climate variability in Egypt. Agriculture 2021, 11, 865. [Google Scholar] [CrossRef]
  9. Bocianowski, J.; Tratwal, A.; Nowosad, K. Genotype by environment interaction for main winter triticale varieties characteristics at two levels of technology using additive main effects and multiplicative interaction model. Euphytica 2021, 217, 26. [Google Scholar] [CrossRef]
  10. Chairi, F.; Aparicio, N.; Serret, M.D.; Araus, J.L. Breeding effects on the genotype× environment interaction for yield of durum wheat grown after the Green Revolution: The case of Spain. Crop J. 2020, 8, 623–634. [Google Scholar] [CrossRef]
  11. Khan, M.M.H.; Rafii, M.Y.; Ramlee, S.I.; Jusoh, M.; Al Mamun, M. AMMI and GGE biplot analysis for yield performance and stability assessment of selected Bambara groundnut (Vigna subterranea L. Verdc.) genotypes under the multi-environmental trials (METs). Sci. Rep. 2021, 11, 22791. [Google Scholar] [CrossRef]
  12. Omrani, A.; Omrani, S.; Khodarahmi, M.; Shojaei, S.H.; Illés, Á.; Bojtor, C.; Mousavi, S.M.N.; Nagy, J. Evaluation of Grain Yield Stability in Some Selected Wheat Genotypes Using AMMI and GGE Biplot Methods. Agronomy 2022, 12, 1130. [Google Scholar] [CrossRef]
  13. Yohane, E.N.; Shimelis, H.; Laing, M.; Mathew, I.; Shayanowako, A. Genotype-by-environment interaction and stability analyses of grain yield in pigeonpea [Cajanus cajan (L.) Millspaugh]. Acta Agric. Scand. Sect. B—Soil Plant Sci. 2021, 71, 145–155. [Google Scholar] [CrossRef]
  14. Yue, H.; Gauch, H.G.; Wei, J.; Xie, J.; Chen, S.; Peng, H.; Bu, J.; Jiang, X. Genotype by Environment Interaction Analysis for Grain Yield and Yield Components of Summer Maize Hybrids across the Huanghuaihai Region in China. Agriculture 2022, 12, 602. [Google Scholar] [CrossRef]
  15. Abdelrahman, M.; Selim, M.E.; ElSayed, M.A.; Ammar, M.H.; Hussein, F.A.; ElKholy, N.K.; ElShamey, E.A.; Khan, N.; Attia, K.A. Developing Novel Rice Genotypes Harboring Specific QTL Alleles Associated with High Grain Yield under Water Shortage Stress. Plants 2021, 10, 2219. [Google Scholar] [CrossRef]
  16. Ferreira Coelho, I.; Peixoto, M.A.; Marcal, T.d.S.; Bernardeli, A.; Silva Alves, R.; de Lima, R.O.; Reis, E.F.d.; Bhering, L.L. Accounting for spatial trends in multi-environment diallel analysis in maize breeding. PLoS ONE 2021, 16, e0258473. [Google Scholar] [CrossRef]
  17. LI, C.-x.; Yun, S.; ZHANG, L.-l. Effects of long-term organic fertilization on soil microbiologic characteristics, yield and sustainable production of winter wheat. J. Integr. Agric. 2018, 17, 210–219. [Google Scholar] [CrossRef] [Green Version]
  18. Zhang, L.; Yuan, J.; Zhang, M.; Zhang, Y.; Wang, L.; Li, J. Long term effects of crop rotation and fertilization on crop yield stability in southeast China. Sci. Rep. 2022, 12, 14234. [Google Scholar] [CrossRef]
  19. Wricke, G. Uber eine Methode zur Erfassung der okologischen Streubreite in Feldverzuchen. Z. Pflanz. 1962, 47, 92–96. [Google Scholar]
  20. Shukla, G. Some statistical aspects of partitioning genotype-environmental components of variability. Heredity 1972, 29, 237–245. [Google Scholar] [CrossRef]
  21. Francis, T.; Kannenberg, L. Yield stability studies in short-season maize. I. A descriptive method for grouping genotypes. Can. J. Plant Sci. 1978, 58, 1029–1034. [Google Scholar] [CrossRef]
  22. Roemer, J. Sinde die ertagdreichen Sorten ertagissicherer. Mitt DLG 1917, 32, 87–89. [Google Scholar]
  23. Plaisted, R.; Peterson, L. A technique for evaluating the ability of selections to yield consistently in different locations or seasons. Am. Potato J. 1959, 36, 381–385. [Google Scholar] [CrossRef]
  24. Plaisted, R. A shorter method for evaluating the ability of selections to yield consistently over locations. Am. Potato J. 1960, 37, 166–172. [Google Scholar] [CrossRef]
  25. Finlay, K.W.; Wilkinson, G. The analysis of adaptation in a plant-breeding programme. Aust. J. Agric. Res. 1963, 14, 742–754. [Google Scholar] [CrossRef] [Green Version]
  26. Zobel, R.W.; Wright, M.J.; Gauch, H.G., Jr. Statistical analysis of a yield trial. Agron. J. 1988, 80, 388–393. [Google Scholar] [CrossRef]
  27. Yan, W.; Cornelius, P.L.; Crossa, J.; Hunt, L. Two types of GGE biplots for analyzing multi-environment trial data. Crop Sci. 2001, 41, 656–663. [Google Scholar] [CrossRef] [Green Version]
  28. Agahi, K.; Ahmadi, J.; Oghan, H.A.; Fotokian, M.H.; Orang, S.F. Analysis of genotype× environment interaction for seed yield in spring oilseed rape using the AMMI model. Crop Breed. Appl. Biotechnol. 2020, 20, e26502012. [Google Scholar] [CrossRef]
  29. Katsenios, N.; Sparangis, P.; Leonidakis, D.; Katsaros, G.; Kakabouki, I.; Vlachakis, D.; Efthimiadou, A. Effect of genotype× Environment interaction on yield of maize Hybrids in Greece Using AMMI Analysis. Agronomy 2021, 11, 479. [Google Scholar] [CrossRef]
  30. Mafouasson, H.N.A.; Gracen, V.; Yeboah, M.A.; Ntsomboh-Ntsefong, G.; Tandzi, L.N.; Mutengwa, C.S. Genotype-by-environment interaction and yield stability of maize single cross hybrids developed from tropical inbred lines. Agronomy 2018, 8, 62. [Google Scholar] [CrossRef] [Green Version]
  31. Gauch, H.G., Jr.; Piepho, H.P.; Annicchiarico, P. Statistical analysis of yield trials by AMMI and GGE: Further considerations. Crop Sci. 2008, 48, 866–889. [Google Scholar] [CrossRef]
  32. Nassar, R.; Huehn, M. Studies on estimation of phenotypic stability: Tests of significance for nonparametric measures of phenotypic stability. Biometrics 1987, 43, 45–53. [Google Scholar] [CrossRef]
  33. Kang, M. A rank-sum method for selecting high-yielding, stable corn genotypes. Cereal Res. Commun. 1988, 16, 113–115. [Google Scholar]
  34. Fox, P.; Skovmand, B.; Thompson, B.; Braun, H.-J.; Cormier, R. Yield and adaptation of hexaploid spring triticale. Euphytica 1990, 47, 57–64. [Google Scholar] [CrossRef]
  35. Thennarasu, K. On Certain Non-Parametric Procedures for Studying Genotype-Environment Inertactions and Yield Stability; IARI, Division of Agricultural Statistics: New Delhi, India, 1995. [Google Scholar]
  36. Farshadfar, E. Incorporation of AMMI stability value and grain yield in a single non-parametric index (GSI) in bread wheat. Pak. J. Biol. Sci. 2008, 11, 1791. [Google Scholar] [CrossRef] [Green Version]
  37. Olivoto, T.; Lúcio, A.D.C. metan: An R package for multi-environment trial analysis. Methods Ecol. Evol. 2020, 11, 783–789. [Google Scholar] [CrossRef] [Green Version]
  38. Mastrodomenico, A.T.; Haegele, J.W.; Seebauer, J.R.; Below, F.E. Yield stability differs in commercial maize hybrids in response to changes in plant density, nitrogen fertility, and environment. Crop Sci. 2018, 58, 230–241. [Google Scholar] [CrossRef] [Green Version]
  39. Stomph, T.; Dordas, C.; Baranger, A.; de Rijk, J.; Dong, B.; Evers, J.; Gu, C.; Li, L.; Simon, J.; Jensen, E.S. Designing intercrops for high yield, yield stability and efficient use of resources: Are there principles? Adv. Agron. 2020, 160, 1–50. [Google Scholar]
  40. Zewdu, Z.; Abebe, T.; Mitiku, T.; Worede, F.; Dessie, A.; Berie, A.; Atnaf, M. Performance evaluation and yield stability of upland rice (Oryza sativa L.) varieties in Ethiopia. Cogent Food Agric. 2020, 6, 1842679. [Google Scholar] [CrossRef]
  41. Zaid, I.U.; Zahra, N.; Habib, M.; Naeem, M.K.; Asghar, U.; Uzair, M.; Latif, A.; Rehman, A.; Ali, G.M.; Khan, M.R. Estimation of Genetic Variances and Stability Components of Yield-Related Traits of Green Super Rice at Multi-Environmental Conditions in Pakistan. Agronomy 2022, 12, 1157. [Google Scholar] [CrossRef]
  42. Gauch, H., Jr.; Zobel, R.W. Optimal replication in selection experiments. Crop Sci. 1996, 36, 838–843. [Google Scholar] [CrossRef]
  43. George, N.; Lundy, M. Quantifying genotype× environment effects in long-term common wheat yield trials from an agroecologically diverse production region. Crop Sci. 2019, 59, 1960–1972. [Google Scholar] [CrossRef]
  44. Gauch, H., Jr. A simple protocol for AMMI analysis of yield trials. Crop Sci. 2013, 53, 1860–1869. [Google Scholar] [CrossRef]
  45. Pour-Aboughadareh, A.; Barati, A.; Koohkan, S.A.; Jabari, M.; Marzoghian, A.; Gholipoor, A.; Shahbazi-Homonloo, K.; Zali, H.; Poodineh, O.; Kheirgo, M. Dissection of genotype-by-environment interaction and yield stability analysis in barley using AMMI model and stability statistics. Bull. Natl. Res. Cent. 2022, 46, 19. [Google Scholar] [CrossRef]
  46. Da Silva, K.J.; Teodoro, P.E.; da Silva, M.J.; Teodoro, L.P.R.; Cardoso, M.J.; Godinho, V.d.P.C.; Mota, J.H.; Simon, G.A.; Tardin, F.D.; da Silva, A.R. Identification of mega-environments for grain sorghum in Brazil using GGE biplot methodology. Agron. J. 2021, 113, 3019–3030. [Google Scholar] [CrossRef]
  47. Bakare, M.A.; Kayondo, S.I.; Aghogho, C.I.; Wolfe, M.D.; Parkes, E.Y.; Kulakow, P.; Egesi, C.; Rabbi, I.Y.; Jannink, J.-L. Exploring genotype by environment interaction on cassava yield and yield related traits using classical statistical methods. PLoS ONE 2022, 17, e0268189. [Google Scholar] [CrossRef]
  48. Omar, M.; Rabie, H.A.; Mowafi, S.A.; Othman, H.T.; El-Moneim, D.A.; Alharbi, K.; Mansour, E.; Ali, M. Multivariate Analysis of Agronomic Traits in Newly Developed Maize Hybrids Grown under Different Agro-Environments. Plants 2022, 11, 1187. [Google Scholar] [CrossRef]
  49. Anila, M.; Mahadeva Swamy, H.; Kale, R.; Bhadana, V.; Anantha, M.; Hajira, S.; Balachiranjeevi, C.; Ayyappa Dass, M.; Bhaskar, S.; Dilip, T. Breeding lines of the Indian mega-rice variety, MTU 1010, possessing protein kinase OsPSTOL (Pup1), show better root system architecture and higher yield in soils with low phosphorus. Mol. Breed. 2018, 38, 147. [Google Scholar] [CrossRef]
  50. Koutroubas, S.; Ntanos, D. Genotypic differences for grain yield and nitrogen utilization in Indica and Japonica rice under Mediterranean conditions. Field Crops Res. 2003, 83, 251–260. [Google Scholar] [CrossRef]
  51. Olivoto, T.; Lúcio, A.D.; da Silva, J.A.; Marchioro, V.S.; de Souza, V.Q.; Jost, E. Mean performance and stability in multi-environment trials I: Combining features of AMMI and BLUP techniques. Agron. J. 2019, 111, 2949–2960. [Google Scholar] [CrossRef]
  52. Ahakpaz, F.; Abdi, H.; Neyestani, E.; Hesami, A.; Mohammadi, B.; Mahmoudi, K.N.; Abedi-Asl, G.; Noshabadi, M.R.J.; Ahakpaz, F.; Alipour, H. Genotype-by-environment interaction analysis for grain yield of barley genotypes under dryland conditions and the role of monthly rainfall. Agric. Water Manag. 2021, 245, 106665. [Google Scholar] [CrossRef]
  53. Pour-Aboughadareh, A.; Sanjani, S.; Nikkhah-Chamanabad, H.; Mehrvar, M.R.; Asadi, A.; Amini, A. Identification of salt-tolerant barley genotypes using multiple-traits index and yield performance at the early growth and maturity stages. Bull. Natl. Res. Cent. 2021, 45, 117. [Google Scholar] [CrossRef]
  54. Piper, C. Soil and Plant Analysis; The University of Adelaide: Adelaide, Australia, 1950; Volume 190, p. 4. [Google Scholar]
  55. Black, C.A. Method of soil analysis part 2. Chem. Microbiol. Prop. 1965, 9, 1387–1388. [Google Scholar]
  56. Team, R.C. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2012. [Google Scholar]
  57. Yan, W.; Kang, M.S. GGE Biplot Analysis: A Graphical Tool for Breeders, Geneticists, and Agronomists; CRC Press: Boca Raton, FL, USA, 2002. [Google Scholar]
  58. Charrad, M.; Ghazzali, N.; Boiteau, V.; Niknafs, A. NbClust: An R package for determining the relevant number of clusters in a data set. J. Stat. Softw. 2014, 61, 1–36. [Google Scholar] [CrossRef]
Figure 1. The AMMI1 and AMMI2 (A,B), respectively) biplots indicating the GEI for the 55 rice genotypes across 4 N environments. The genotypes’ and environments’ legends are presented in Table 2.
Figure 1. The AMMI1 and AMMI2 (A,B), respectively) biplots indicating the GEI for the 55 rice genotypes across 4 N environments. The genotypes’ and environments’ legends are presented in Table 2.
Plants 11 02775 g001
Figure 2. Which-won-where pattern based on the GGE biplot polygon to identify the best rice genotypes under the four different N environments (A). The GY × WAAS statistic biplot for selecting the high-yielding and stable rice genotypes (B).
Figure 2. Which-won-where pattern based on the GGE biplot polygon to identify the best rice genotypes under the four different N environments (A). The GY × WAAS statistic biplot for selecting the high-yielding and stable rice genotypes (B).
Plants 11 02775 g002
Figure 3. Hierarchical classification dendrogram of the 55 tested rice genotypes based on their ranks for GY and stability statistics via Ward’s method. SCL 1 (light blue), SCL 2 (Blue), SCL3 (green), and SCL 4 (red).
Figure 3. Hierarchical classification dendrogram of the 55 tested rice genotypes based on their ranks for GY and stability statistics via Ward’s method. SCL 1 (light blue), SCL 2 (Blue), SCL3 (green), and SCL 4 (red).
Plants 11 02775 g003
Table 1. AMMI analysis based on GY means for genotypic performance and N environment EPC1 value and the 5 top ranking genotypes for each N environment.
Table 1. AMMI analysis based on GY means for genotypic performance and N environment EPC1 value and the 5 top ranking genotypes for each N environment.
EMeanEPC1IndexClass12345
0N19.13842−2.35395−4.575242UnfavorableIR22GIZA14C22MTU1010SK2058
LN21.94194−1.39131−1.771727UnfavorableIR22MTU1010Yabani LuluSK2035IR66
MN24.936120.0545211.222455favorableSK2046IR22MTU1010Egyptian YasminGiza178
HN28.838183.6907415.124515favorableSK2058MTU1010ReihoWAB 880 SG 73Giza178
Table 2. Name, type, parentage, and origin of the plant materials.
Table 2. Name, type, parentage, and origin of the plant materials.
No.GenotypeTypeParentageOrigin
1SabienyJSelection from IntroductionsEGYPT
2Nabatat AsmarJSelection from Agami M1EGYPT
3Giza 159JGiza14/Agami M.1EGYPT
4Yabani LuLuJSelection from IntroductionsEGYPT
5GZ 5830-59-10-2JGZ4120/Suweon349EGYPT
6Giza 14JYabani Pearl/Iraki16EGYPT
7GZ 7718-13-3-1-3JSakha101/HR4856-1-1-2EGYPT
8NahdaJSelection from IntroductionsEGYPT
9GZ 6214-4-1-1-1JGZ4122-23-4-2/IRI396EGYPT
10Sakha 101JGiza176/Milyang79EGYPT
11IR 68373-R-R-B-22-2-2T.J.JINMIBYEO/YR14987-91IRRI
12Giza 182IGiza181/IR39422//Giza181EGYPT
13GZ 7922-B-44-1JGiza177/IDSAEGYPT
14Pusa Basmati 1IIndia selectionINDIA
15Sakha 104JGZ4096/GZ4100EGYPT
16IR 28IIR8333-6-2-1///IR1561-149-1//IR24*4/O. NIVARAIRRI
17GZ 6522-15-1-1-3JGZ5581/GZ4316EGYPT
18IR 64IIR5657-33-2-1/IR2061-4665-1-5-5IRRI
19GZ 7718-13-3-2-2JSakha101/HR4856-1-1-2EGYPT
20Giza 177JGiza171/Yamji No.1//PI NO.4EGYPT
21IR 66IIR13240-108-2-2-3/IR9129-209-2-2-2-1IRRI
22GZ 6910-28-1-3-1JSakha101/GZ24316(MUT)EGYPT
23IR 70IIR19660-73-4/IR54//IR9828-36-3IRRI
24Agami M.1JSelection from cultivated varietiesEGYPT
25Sakha 103JGZ4120/Suweon349EGYPT
26ArabiI/JJava3/Yabani Montkhab 3EGYPT
27Milyang 63I/JTONGIL/IR946-33-2-2-2//YR675-131-2KOREA
28Yen Geng 135JChinese selectionCHINA
29IR 73689-31-1T.J.SR18977-TB-4/JINMIBYEOIRRI
30Giza 178I/JGiza175/Milyang49EGYPT
31WAB 450-1-B-P-91-HBI---Africa Rice
32BG 304I---SRILANKA
33MTU 1010I---INDIA
34IR 68353-35-3-3-2-2-1-2T.J.CHEOLWEON49/KYWHA9IRRI
35Giza 175I/J(IR28/IR1541)/(Giza180/Giza14)EGYPT
36WAB 880 SG 73I---Africa Rice
37E 7034JEWAN NO.5/857CHINA
38SKC 23808-28-5-2-1-1I/J98-Y-116/Sakha102IRRI
39IET 1444ITN1/CO.29INDIA
40Black RiceJJingo9601China
41IR 7421-35-1-1-2T.J.IR2035-290-2-1-1/MASINOIRRI
42GZ 6903-3-4-2-1JSakha101/Suweon313EGYPT
43SKC 23822-304-3-1-1-1I/JM202/Giza177IRRI
44Taikeng Yu 1420JC253///J692130/BL6//TAINUNG67/IR4547-2-1-2TAIWAN
45Egyptian YasmineIIR262-43-8-11/KDML105EGYPT
46IR 67075-2B-5-2IIR10198-66-2//GZ2175/CSR1IRRI
47IR 74IIR19661-131-1-2/IR15795-199-3-3IRRI
48ReihoJHOYOKU/AYANISHKIJAPAN
49C 22ITJRERMAS/BPI76//PALAWAN/AZUCENAIRRI
50Yun Lu No. 48JLUYIN NO.7/YUNANJINGDAO-38CHINA
51IR 22IPETA/DEE GEO WOO GEN//TADUKANIRRI
52SK2034IIR69625A/Giza178REGYPT
53SK2046IIR69625A/Giza181REGYPT
54SK2035IIR70368A/Giza178REGYPT
55SK2058IIR69625A/Giza182REGYPT
J, Japonica; I, Indica; T.J., Tropical japonica; I/J, Indica japonica; IRRI, International Rice Research Institute.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abdelrahman, M.; Alharbi, K.; El-Denary, M.E.; Abd El-Megeed, T.; Naeem, E.-S.; Monir, S.; Al-Shaye, N.A.; Ammar, M.H.; Attia, K.; Dora, S.A.; et al. Detection of Superior Rice Genotypes and Yield Stability under Different Nitrogen Levels Using AMMI Model and Stability Statistics. Plants 2022, 11, 2775. https://doi.org/10.3390/plants11202775

AMA Style

Abdelrahman M, Alharbi K, El-Denary ME, Abd El-Megeed T, Naeem E-S, Monir S, Al-Shaye NA, Ammar MH, Attia K, Dora SA, et al. Detection of Superior Rice Genotypes and Yield Stability under Different Nitrogen Levels Using AMMI Model and Stability Statistics. Plants. 2022; 11(20):2775. https://doi.org/10.3390/plants11202775

Chicago/Turabian Style

Abdelrahman, Mohamed, Khadiga Alharbi, Medhat E. El-Denary, Taher Abd El-Megeed, El-Sayed Naeem, Samah Monir, Najla A. Al-Shaye, Megahed H. Ammar, Kotb Attia, Said A. Dora, and et al. 2022. "Detection of Superior Rice Genotypes and Yield Stability under Different Nitrogen Levels Using AMMI Model and Stability Statistics" Plants 11, no. 20: 2775. https://doi.org/10.3390/plants11202775

APA Style

Abdelrahman, M., Alharbi, K., El-Denary, M. E., Abd El-Megeed, T., Naeem, E. -S., Monir, S., Al-Shaye, N. A., Ammar, M. H., Attia, K., Dora, S. A., & Draz, A. -S. E. (2022). Detection of Superior Rice Genotypes and Yield Stability under Different Nitrogen Levels Using AMMI Model and Stability Statistics. Plants, 11(20), 2775. https://doi.org/10.3390/plants11202775

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop