Next Article in Journal
Changes in Soil Properties, Microbial Quantity and Enzyme Activities in Four Castanopsis hystrix Forest Types in Subtropical China
Next Article in Special Issue
Identification of High-Yielding Genotypes of Barley in the Warm Regions of Iran
Previous Article in Journal
Hydrochorous Seed Transport in the Lower Traisen River before and after Riverbed Restoration
Previous Article in Special Issue
Future Climate Effects on Yield and Mortality of Conventional versus Modified Oil Palm in SE Asia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Selection of High-Yielding and Stable Genotypes of Barley for the Cold Climate in Iran

by
Alireza Pour-Aboughadareh
1,*,
Habibollah Ghazvini
1,*,
Seyed Shahriyar Jasemi
1,
Solaiman Mohammadi
2,
Sayed Alireza Razavi
3,
Mehrdad Chaichi
4,
Marefat Ghasemi Kalkhoran
5,
Hassan Monirifar
6,
Hamid Tajali
3,
Asadollah Fathihafshjani
7 and
Jan Bocianowski
8,*
1
Seed and Plant Improvement Institute, Agricultural Research, Education and Extension Organization (AREEO), Karaj P.O. Box 31587-77871, Iran
2
Field and Horticultural Crops Research Department, Agricultural and Natural Resources Research and Education Center of West-Azarbayjan Province, Agricultural Research, Education and Extension Organization, Urmia P.O. Box 57169-63963, Iran
3
Field and Horticultural Crops Research Department, Agricultural and Natural Resources Research and Education Center of Khorasan Razavi Province, Agricultural Research, Education and Extension Organization, Mashhad P.O. Box 91769-83641, Iran
4
Field and Horticultural Crops Research Department, Agricultural and Natural Resources Research and Education Center of Hamedan Province, Agricultural Research, Education and Extension Organization, Hamedan P.O. Box 65199-91169, Iran
5
Field and Horticultural Crops Research Department, Agricultural and Natural Resources Research and Education Center of Ardabil (Moghan) Province, Agricultural Research, Education and Extension Organization, Ardabil P.O. Box 56951-57451, Iran
6
Crop and Horticultural Science Research Department, East Azarbaijan Agricultural and Natural Resources Research and Education Center, Agricultural Research, Education and Extension Organization (AREEO), Tabriz P.O. Box 51537-15898, Iran
7
Field and Horticultural Crops Research Department, Agricultural and Natural Resources Research and Education Center of Markazi Province, Agricultural Research, Education and Extension Organization, Arak P.O. Box 38135-889, Iran
8
Department of Mathematical and Statistical Methods, Poznań University of Life Sciences, Wojska Polskiego 28, 60-637 Poznań, Poland
*
Authors to whom correspondence should be addressed.
Plants 2023, 12(13), 2410; https://doi.org/10.3390/plants12132410
Submission received: 28 May 2023 / Revised: 19 June 2023 / Accepted: 19 June 2023 / Published: 22 June 2023

Abstract

:
The interaction between genotypes and environments plays an important role in selecting superior genotypes for target locations. The main objectives of the present study were to analyze the effect of the genotype-by-environment interaction (GEI) and identify superior, newly developed, and promising barley genotypes for cold regions in Iran. For these purposes, a set of genotypes obtained from breeding programs for cold climates in Iran, along with two reference genotypes, were investigated at eight research stations (Tabriz, Ardabil, Arak, Miandoab, Mashhad, Jolge Rokh, Karaj, and Hamadan) during two consecutive growing seasons (2019–2020 and 2020–2021). The results of the freezing test (LT50) showed that most of the tested genotypes had significant cold tolerance at the seedling stage. Based on the additive main effect and multiplicative interaction (AMMI) analysis, environment (E) and GEI effects explained 49.44% and 16.55% of the total variation in grain yield, respectively. Using AMMI1 and AMMI2 models, G2 and G20 were found to be superior genotypes in terms of grain yield and stability. Moreover, AMMI-based stability parameters considered the G20 genotype to be the ideal genotype. A two-plot analysis of the genotype-by-environment interaction (GGE) biplot showed that the 16 experimental environments were grouped into 2 mega-environments. Of the test environments, ARK1 and KAJ2 had the highest discriminating power and representativeness ability, and these were identified as ideal environments for testing advanced genotypes for yield and stability performance during early barley breeding practices in cold areas in Iran. In conclusion, both AMMI and GGE biplot models identified several superior genotypes, among which G20, with a high average yield relative to the overall average yield and the lowest IPC1 score, was found to have high yield stability and is recommended for inclusion in breeding programs for cold climates in Iran.

1. Introduction

Barley (Hordeum vulgare L.) is widely adapted to various environmental conditions. This cereal has a higher tolerance to environmental stresses such as salinity, drought, cold, etc. [1]. Barley is the fourth cereal crop in the world, and according to an FAO report, it was estimated that the average harvested area and global grain production would be ~504,000 hectares and 3.50 tons ha−1 in 2019, respectively (Faostat, http://www.fao.org/faostat/en/#home accessed on 1 December 2021). It was proven that barley kernels can be considered a wholesome food commodity, as they provide minerals, calcium and phosphorus, a moderate amount of protein and fiber, and B vitamins [2]. In Iran, and even in some developing countries, barley is known as the second most important food and feed crop after wheat [3]. Iran was one of the ten countries with the highest barley seed production in 2020 [4]. According to the latest statistics, in the 2020–2021 crop season, barley was grown on about 1.44 million hectares in the country, and its grain production was about 2.47 million tons [5]. The mountainous regions in Iran are considered one of the main barley growing centers, where low-temperature or cold stress limits barley production [6]. In these areas, soil temperatures can drop to near zero or even below zero [0 to −4 °C] in winter, which can severely damage seedlings. Therefore, screening barley germplasm for cold tolerance and identifying the most tolerant genotypes at an early stage of plant growth and development is a key breeding strategy for extending barley cultivation in the highlands and cold regions in Iran. Another important criterion for evaluating cold tolerance is to test the stability of grain yield for desirable genotypes in different locations with freezing winter temperatures [7].
As a general rule of thumb, accounting for the differential response of genotypes to environmental factors, especially climatic factors, has been one of the key challenges in plant breeding for decades. One of the major steps toward developing new varieties adapted to a wide range of environments is the evaluation of genotype–environment interactions (GEIs) [8,9]. When crop genotypes are evaluated in different environments (years, locations, and/or their combination), their yield performances can vary significantly [10]. The presence of a large GEI effect often leads to serious challenges for breeders when selecting high-yielding genotypes with overall adaptation [11]. Like other crops, barley is sensitive to environmental changes. According to Khalili and Pour-Aboughadareh [12], Vaezi et al. [13], Ahakpaz et al. [14], Hilmarsson et al. [15], Ghazvini et al. [6], and Nykiel et al. [16], barley yield is strongly affected by environmental changes. This means that GEI analysis plays a key role in commercial purposes and in identifying ideal target production environments to achieve maximum yield. Consequently, having a comprehensive knowledge of the GEI effect, measurements of adaptability and stability, and the influence of the genetic proportion of genotypes and environments on yield performance is essential for breeders when identifying the locations that should be used in their breeding programs [17].
Multi-environment trials (METs) are an important part of any breeding program aiming to identify and release high-yielding and stable genotypes for large-scale production in different environments, according to farmers’ preferences. However, the proper use of these trials requires an understanding of GEI effects [18]. There are many statistical tools for analyzing GEI and identifying stable high-yielding genotypes in different environments. Of all the methods, the additive main effect and multiplicative interaction (AMMI) analysis [19] and the genotype and genotype-by-environment interaction (GGE) biplot methodology [20] are the best models for analyzing such complex data sets. There are two main advantages of using the AMMI model: first, the model partitions the total variation into the main effects of genotypes and environments, as well as their interactions. Indeed, because of this flexibility, complex variation can be studied easily and separately. Second, it provides a pathway to predictive accuracy using the family members model [21]. Overall, the AMMI model allows breeders to use its valuable agricultural recommendation results to leverage both broad and narrow adaptations to increase productivity [22]. On the other hand, the GGE biplot model is a very useful tool and provides graphic images of MET data sets. This methodology is useful for (i) the selection of a high-performance genotype in a production environment, (ii) identifying the target environment for a specific genotype, (iii) determining the representativeness power and discriminating ability of test environments for further genotype assessments, (iv) deciphering the relationships between test environments, and (v) ranking the genotypes tested in term of both stability and yield performance [20]. Most GEI studies on barley have been performed in dry, temperate, and warm climates in Iran [1,12,13,14,23,24,25,26,27,28], and few were performed in cold regions. However, studies conducted in the cold climate in Iran have shown significant GEI in various regions [6,7] and supported the use of the AMMI analysis and the GGE biplot methodology for GEI analysis. With this in mind, the aims of this study were to (i) demonstrate the value of combining AMMI and GGE biplot models for GEI analysis and interpretation, (ii) identify actual mega-environments in cold regions in Iran, and (iii) assess the adaptability and stability of newly promising barley genotypes in different barley production environments.

2. Results

2.1. Response of Investigated Barley Genotypes to the Freezing Test

The results of the analysis of variance for LT50, thousand-grain weight, and grain yield are shown in Table 1. Accordingly, a significant difference was found between G4 and genotypes G12, G15, and G17 in terms of the LT50 indicator. On the contrary, differences between the tested genotypes were significant for thousand-grain weight and grain yield. LT50 ranged from −10 °C to −6 °C with an average of −7.58 °C, and among the genotypes tested, G4, with LT50 −10 °C, was the genotype with the best cold tolerance at the seedling stage. In addition, the thousand-grain weight ranged from 34.50 to 46.60 g with an average of 40.14 g. The highest values of this parameter were observed in genotypes G16, G17, and G19 (reference genotype). Grain yield showed high variability and ranged from 5.27 to 8.71 with an average of 7.13 tons ha−1. Genotype G1 (reference genotype), followed by G20 and G7, showed the highest grain yield compared to the other genotypes tested.

2.2. AMMI Analysis

The result of the AMMI analysis showed that the effects of genotype (G), environment (E), and their interaction were highly significant (Table 2). The main effect of E accounted for 49.44% of the total variation, while G and GEI effects explained 4.99% and 16.55% of the variation, respectively. The GEI effect was divided into five significant IPCAs, each of which explained 22.98%, 18.20%, 13.48%, 11.49%, and 9.19%, respectively. AMMI-based biplots were drawn to dissect the explained variation due to the GEI effect. The AMMI1 biplot represented the main effect of genotypes and environments on the IPCA1 score. Among the genotypes, G2, G6, G11, and G16 showed the highest IPCA1 scores and were therefore considered the least stable genotypes with large distances from the origin of the biplot; hence, these genotypes showed a high interaction with the test environments (Figure 1A).
The first four genotypes recommended for each test environment were identified using the AMM2 model. As shown in Table 3, genotype G2 ranked first in five environments (HAM1, HAM2, JOL2, ARK2, and MIN1) and appeared in eight out of sixteen environments in the top four rankings; genotype G20 ranked in eight out of sixteen environments in the top four rankings and was the first genotype in environment ARD1; genotype G4 ranked first rank in two environments (TAB1 and KAJ2) and appeared in the top four in other two environments (MAS2 and ARK2); genotype G18 was identified in the top four in four environments (TAB2, MAS2, MIN2, and JOL1); and genotype G15 was observed in the top four in six environments (MIN1, MIN2, KAJ1, ARD1, ARD2, and HAM1) and as dominant in the KAJ1 environment. Additional results are shown in Table 3.
The results for the estimated AMMI-based stability parameters along with the average grain yield for each genotype tested are shown in Table 4. Based on grain yield, the reference genotype (cv. Jolge), followed by G20, G2, G18, G3, and G4, showed the best average grain yield among the 16 test environments, while G10, G14, G17, G8, and G11 showed the lowest productivity. Genotype G20, followed by G14, G12, G17, and G6, showed higher stability compared to the other tested genotypes based on their low AMMI stability values (ASVs), average squared eigenvector values (EVs), AMMI-based stability parameter (ASTAB), the average sum between environments of the absolute value of GEI modeled using AMMI (AVAMGE), and IPCA point distance from the origin in space (Dz). According to the fitted AMMI model parameter (FA), genotypes G9, G12, G14, G17, and G20 were the most stable genotypes. The AMMI stability index (ASI) indicated that genotypes G20, G14, G12, G17, and G7 were the most stable genotypes. Based on the distance between the IPCA points and the origin in space (Da), genotypes G5, G9, G14, G17, and G20 showed the greatest stability compared to the other genotypes.
The results of the principal component analysis (PCA) showed that the first two PCAs with eigenvalues greater than 1 explained 90.68% of the total variation in grain yield and stability parameters. The biplot constructed using PCA1 and PCA2 was used to evaluate the correlation between stability parameters and grain yield (Figure 1B). Thus, the large cosine angle between the vectors of the two parameters showed a weak correlation, while the small cosine angle between these vectors showed a strong positive correlation. In contrast, no relationship was displayed at a cosine of 90°, while a strong negative correlation was displayed at a cosine of 180°. Based on the results, positive and strong correlations were observed between AVAMGE, FA, EV, DZ, DA, and ASTAB. On the other hand, grain yield was positively and strongly correlated with the ASV and ASI parameters. Therefore, the latter stability parameters can be expressed as a ‘dynamic’ stability concept.

2.3. GGE Biplots to Visualize GEI

The results of the GGE biplot analysis showed that the two PCAs accounted for 48.06% of the total grain yield variation in the different test environments. The polygon view of the GGE biplot clustered the test environments into two of the six sectors (Figure 2A). The environments of Tabriz (TAB1 and TAB2), Karaj (KAJ1 and KAJ2), and Mashhad (MAS1 and MAS2) along with MIN1 (Miandoab—first year) and ARK1 (Arak—first year) were placed in one sector and formed mega-environment I. The other environments, including Hamadan (HAM1 and HAM2), Joleh-Rokh (JOL1 and JOL2), Ardabil (ARD1 and ARD2), and the second-year data for Arak and Miandoab (ARK2 and MIN2), were considered mega-environment II. On the other hand, the best genotypes were G1 (reference genotype: Jolege), G2, G4, G6, G10, and G11. Among them, G16, G4, and G1 were the highest-yielding and most adapted genotypes for the first mega-environment. In contrast, the G2 genotype was specifically adapted to mega-environment II. Genotypes G6, G10, G11, and G16 showed no specific adaptation to any of the test environments. The “mean vs. stability” biplot viewpoint indicated that G1 (the reference genotype), followed by G20 and G2, had the highest average grain yield in the test environments. Genotypes G3, G14, G18, and G19 showed yields closest to the grand mean value due to their position in the biplot. The G20 genotype, with the high average grain yield, was the most stable, while G2, G4, and G19 showed significant yield variability in different environments. However, some genotypes, such as G10 and G14, with a low average yield, showed high stability (Figure 2B).
Figure 2C shows the discrimination power and representativeness of the test environments. Based on this biplot, test environments TAB1, TAB2, ARD1, MIN1, ARK1, and KAJ2 showed the highest value of discrimination power due to the large length of the environment vectors. Moreover, the representative ability of the test environments was determined using the average environment coordinate (AEC) and the angle between the test environment vectors. Accordingly, the ARK1, KAJ1, KAJ2, and MAS2 environments with smaller angles showed relatively strong representativeness, while MIN2, TAB1, ARD1, and JOL1 showed relatively weak representativeness. Therefore, KAJ2 and ARK1 were considered ideal environments due to their discrimination strength and representativeness potential. A comparative view of the GGE biplot was used to identify ideal genotypes (Figure 2D). According to the theory of Yan and Kang [20], an ideal genotype has the highest yield and stability in all test environments. Of the barley genotypes tested, genotypes G20, G1, G18, G3, G4, G19, and G2 were near the average environment axis (AEA) and were identified as ideal genotypes. Of these, genotype G20 was closest to the central circle in the biplot and showed specific adaptability to ARK1 and KRJ2 environments. However, genotypes G8, G10, G14, and G17 were found to be the most undesirable compared to other genotypes.

3. Discussion

In recent years, major parts of crop fields have been affected by various environmental stresses caused by climate change. These changes can have a negative impact on crop productivity and food security. Of the cereal crops, barley has a high tolerance to abiotic stresses, and this trait has led to its cultivation in highly variable weather environments and even in marginal crop fields. Multi-environment trials (METs) provide an opportunity to dissect GEI effects as well as identify target environments for economical seed production [29]. Moreover, these trials play a key role in identifying superior genotypes and determining their general and specific adaptations to different environments [30]. To obtain superior genotypes with high yield and stability, two strategies have been proposed: (i) dividing mega-environments into homogeneous regions to obtain varieties with specific adaptations and (ii) identifying genotypes with general adaptability in different environments [31]. Based on the literature, both strategies have been widely used to develop stable varieties with high performance [8,15,32,33,34,35,36]. In this study, several new barley genotypes were evaluated in 16 cold environments in Iran. As shown in Figure 3, some test environments are located in regions with cold semi-desert climates (e.g., Hamadan, Arak, Jolge-Rokh, and some parts in Arak); on the other hand, some environments (e.g., Tabriz, Ardabil, and Mashhad) have a wide range of weather conditions (from Mediterranean to cold). In addition, it is worth noting that the largest barley-growing areas in Iran are in cold regions. Therefore, planning breeding projects to identify and develop genotypes with specific and general adaptability is essential to achieve a high level of barley production in Iran.
The results of the AMMI analysis indicated that the effects of the environment (E) and the genotype-by-environment interaction (GEI) accounted for a greater proportion of the total variability than the effect of genotype (Table 2). In line with our results, previous studies have also shown that the environment was the main source of variation in barley samples under different environmental conditions [6,13,15,27,35]. The magnitude of the total sum of squares (TSS) for the GEI effect indicated significant differences in the genotypic response of barley to environmental conditions in cold regions in Iran. Therefore, the adaptation patterns and stability of the genotypes studied were analyzed using different methods. The AMMI1 biplot provided an opportunity to identify high-yielding and stable genotypes due to the distribution of genotypes based on the IPCA1 results and average grain yield. Based on this biplot, the G20, G14, and G10 genotypes, with high average yield relative to the overall average yield and the lowest IPC1 scores, were found to have high yield stability (Figure 1A). Moreover, this result was confirmed using the AMMI2 model, in which the G2 and G20 genotypes were found to be in the top four genotypes in eight of the sixteen test environments (Table 3). The results for the AMMI-based stability parameters showed that genotype G20, with an average grain yield of 7.59 tons ha1, was the most stable genotype compared to the others. Using another AMMI statistic, we found that the ASV and ASI positively and significantly correlated with grain yield, which in turn revealed a ‘dynamic’ concept of stability (Figure 1B). Consequently, barley breeders can use these parameters to select genotypes that have a stable yield similar to the average yield of all genotypes in the tested environments [37].
The biplot of polygons (Figure 2A) classified the test environments into two mega-environments. The first mega-environments included Tabriz (TAB1 and TAB2), Karaj (KAJ1 and KAJ2), and Mashhad (MAS1 and MAS2) along with MIN1 (Miandoab—first year) and ARK1 (Arak—first year), for which genotypes G1 and G4 were identified as the best-adapted genotypes. In addition, genotype G20, with acceptable yield and yield stability, was placed in this mega-environment. The Hamadan (HAM1 and HAM2), Jolgeh-Rokh (JOL1 and JOL2), and Ardabil (ARD1 and ARD2) test environments and the second-year data for Arak and Miandoab (ARK2 and MIN2, respectively) formed the second mega-environment, where the G2 genotype was identified as the best adapted to them. Overall, the clustering of test environments followed the climatic clustering scheme (Supplementary Figure S1). The results of the cluster analysis based on temperature data (minimum, maximum, and average) grouped all test environments into two main clusters. The first cluster included Hamadan, Jolgeh-Rokh, and Ardabil, while the second consisted of Mashhad, Karaj, Arak, and Tabriz. These results were consistent with the results of the mega-environments clustering, so it can be concluded that the tested genotypes can specifically respond to environmental changes. In other words, the obtained mega-environment pattern in this study was related to the classification obtained from the weather data, and this finding was confirmed using the Mantel test (r = 0.548 **) (Figure S2). Moreover, these results were further supported with the AMMI model, in which fourth-order genotypes were identified in target environments. For example, G2 was identified as a high-yielding genotype in the Karaj and Mashhad environments, while G2 was selected for the Hamadan and Jolgeh-Rokh environments (Table 3). Therefore, these results may confirm the possibility of two microclimates in the cold regions in Iran. However, this prediction requires further research using other data sets from different years as well as different field crops.
These results were further supported by the freezing test. As shown in Table 1, genotype G4 appeared to be the most tolerant and produced significant grain yield under field conditions (freezing test at Karaj station) compared to the other genotypes. Therefore, this genotype can be recommended for cultivation in target environments (Tabriz, Mashhad, and Karaj). In addition, the G2 genotype, with acceptable cold tolerance, was suitable for cultivation in the second mega-environment. Since the target environments belonging to this mega-environment are located in semi-arid regions, it is recommended to evaluate the response of this genotype to terminal drought stress.
According to Yan [38], the environment with the most discriminating and representative power can be considered the ideal target environment for evaluating new varieties for their full yield potential. According to this theory, test environments can be divided into three types. Type I includes environments that provide little information on genotypes due to their short vectors; therefore, they should not be used as test environments. Type II includes environments that are suitable for identifying the best high-yielding genotypes due to their small angles with the AEC and long vectors. Type III consists of environments that are useful for eliminating unstable genotypes due to their large angles with AECs and long vectors. In this study, the KAJ2 and ARK1 environments had high discriminating power and representativeness potential. Therefore, these environments can be recommended as desirable environments (Type II) for testing advanced genotypes for yield and yield stability during early breeding programs in cold regions in Iran. On the other hand, among the test locations, the Mashhad location (MAS1 and MAS2) was a Type I environment; therefore, this location can be omitted from stability and adaptability experiments to reduce the cost of field evaluations (Figure 2C). Another result of the GGE biplot was a positive and significant correlation between the Hamadan (HAM1 and HAM2), Jolgeh-Rokh (JOL1 and JOL2), and Arak (ARK1 and ARK2) environments, suggesting that two of the three locations can be removed from barley METs for cold regions in Iran. Another comparison view of the GGE biplot identified the ideal genotypes. Theoretically, genotypes near the beginning of the biplot may have broader adaptation and could be identified as ideal genotypes [39]. Based on this biplot, genotype G20, followed by G1, G18, G3, G4, G19, and G2 were placed closest to the AEA cutoff and selected as ideal genotypes. Of these, the G20 genotype with good correlation in ARK1 and KAJ2 environments showed specific adaptability to these environments (Figure 2D). In summary, our results revealed that both AMMI and GGE models deciphered a similar pattern for determining high-yielding and stable barley genotypes. Therefore, the relative contribution of both models to the identification of superior genotypes in this study was consistent with previous studies on barley and other crops [13,26,27,40,41].
Genotypes that are best suited to specific environmental conditions can be detected using AMMI analyses, which can estimate the effect of the genotype interaction in each environment. For yield, a significant GEI was demonstrated using an AMMI analysis. The high stability of genotypes is related to the AMMI stability value. Determination of the main effect of the genotype, environment, and the most significant GEIs can be completed using AMMI results displayed on the GE biplot. AMMI models are able to measure the importance of environments, genotypes, and their interactions using a value that measures genotype stability across all environments, given the yield [42]. Genotype G20, with the high average yield relative to the overall average yield and the lowest IPC1 score, was shown to have high yield stability and thus is recommended for inclusion in breeding programs.

4. Materials and Methods

4.1. Plant Materials, Field Layout, and Experimental Design

A set of eighteen promising barley genotypes, along with two introduced varieties (Jolge (G1) and Mahtab (G19)) used as reference genotypes, were studied using multi-environment trials (Table 5). Field evaluations were carried out at eight cold-weather test stations (including Ardabil, Miandoab, Arak, Hamadan, Mashhad, Jolge Rokh, Tabriz, and Karaj) in Iran for two consecutive growing seasons (2019–20/2020–21) (Figure 3). Additional information on geographical and climatic data is presented in Supplementary Table S1. In all test locations, field trials were conducted using a randomized complete block design (RCBD) in three replications. The experimental plots consisted of six rows, each 6 m in length with 20 cm spacing between rows. The sowing density in each plot was 400 seeds per m2. Sowing was carried out using an experimental planter (Wintersteiger, Ried, Austria). At each station, basic fertilizers such as N and P2O5 were applied at 32 and 100 kg ha1, respectively, before sowing. During the growing season, five irrigations were applied at growth stages 00, 32, 51, 75, and 85 Zadoks’ [43] at all test locations. In addition, 40 kg ha1 of N was again applied at the stem elongation stage (ZGS 31). At harvest time (based on the time of physiological maturity at each test station), a combine harvester (Wintersteiger, Ried, Austria) was used to harvest the experimental plots. Grain yield was determined for each genotype at the test locations and used for further statistical analysis.

4.2. Evaluation of the Lethal Temperature (LT50) for Genotypes

All genotypes were tested for cold hardiness under field conditions and in a programmable freezer. The field experiment was performed at the Cereal Research Department, Seed and Plant Improvement Institute (SPII) Karaj, Iran (longitude 35° N, latitude 51° N, and altitude of 1132 m above sea level), which is characterized by a moderately cool temperate region. Plant materials were arranged in a randomized complete block design (RCBD) with two replicates. Experimental plots consisted of six rows, each 6 m long, with 20 cm spacing between rows. The seed density in each plot was 450 seeds per m2. Seeds were sown using an experimental planter (Wintersteiger, Ried, Austria) on 3 November 2018. After germination and early growth under field conditions, cold-acclimated seedlings were sampled for freeze testing during the first week of January 2019. Freeze tests to estimate the LT50 (the temperature at which 50% of the plants perish due to cold weather) were performed as described by Limin and Fowler [44]. Plant crowns were planted in pots containing moist sand and transferred to a programmable freezer, which was maintained at −3° C for 12 h. The freezing temperature was lowered in a constant gradient of −2 °C per hour until it reached −21 °C. At 2 °C intervals for each of the ten test temperatures (from −3 °C to −21 °C), five plants from each genotype from both replications were selected and transferred to the outside of the freezer. The crowns of the plants were replanted, and LT50 was recorded based on regrowth after 21 days under controlled glasshouse conditions. At harvest time, the remaining experimental plots for each genotype were harvested under field conditions using an experimental combine harvester (Wintersteiger, Ried, Austria). Finally, thousand-grain weight and grain yield were recorded for each genotype.

4.3. Statistical Methods

Grain yield data obtained from the 16 test environments (a combination of 2 cropping seasons and 8 test locations) were subjected to AMMI analysis. Several AMMI-based stability indices were used to determine the phenotypic adaptability and stability of the barley genotypes tested. To assess the relationships between the estimated stability indices and grain yield, a principal component analysis (PCA) was performed. In addition, grain yield data were graphically analyzed using a GGE biplot model. The biplots were rendered based on the first two principal components (PCAs) obtained using singular value decomposition (SVD). All statistical analyses were calculated using the packages ‘pheatmp’ [45], ‘metan’ [46], ‘FactoMineR’ [47], and ‘factoextra’ [48] in R software [49]. In addition, a cluster analysis based on Ward’s algorithm was conducted to group the test environments using weather data. To examine the relationship between the clustering pattern in the environment and the identified mega-environments, a Mantel test was calculated using XLSTAT software (XLSTAT, Addisonsoft, Paris, France).

5. Conclusions

The collective analysis using the AMMI and GGE biplot models identified two mega-environments for barley-growing cold regions in Iran, which in turn provided important implications for future barley breeding programs. Moreover, the promising new genotypes examined in this study should be recommended for each mega-environment to improve productivity and stability. When considering the results obtained in this study, the high-yielding and stable genotypes in each mega-environment were generally different, which in turn was in accordance with a recommendation from Gauch [22], who noted that genotype stability is a meaningful goal only in each mega-environment and not in multiple mega-environments. In conclusion, our results showed that the G2 and G20 genotypes can be recommended as newly promising barley genotypes for further investigation before commercial release.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12132410/s1, Table S1. Agro-climatic characteristics of the test environments during the 2019–2021 cropping years in the cold regions in Iran; Figure S1. The grouping pattern in the test environments based on weather data during the 2019–2021 cropping years in the cold regions in Iran. ARD, ARK, HAM, MIN, MAS, KAJ, JOL, and TAB indicate Ardabil, Arak, Hamadan, Jolgeh-Rokh, Karaj, Mashhad, Miandoab, and Tabriz locations, respectively. Numbers 1 and 2 indicate the first and second cropping years, respectively; Figure S2. Mantel’s test results between the estimated proximity matrices (Euclidean distance) for the grouping patterns in the test environments based on weather data and grain yield.

Author Contributions

Conceptualization, H.G. and A.P.-A.; methodology, H.G. and S.S.J.; software, A.P.-A.; validation, H.G., S.S.J. and A.P.-A.; formal analysis, A.P.-A.; investigation, S.S.J., S.M., S.A.R., M.C., M.G.K., H.M., H.T. and A.F.; resources, H.G.; data curation, H.G. and A.P.-A.; writing—original draft preparation, A.P.-A.; writing—review and editing, H.G., J.B. and A.P.-A.; visualization, H.G. and A.P.-A.; supervision, H.G.; project administration, H.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data in this manuscript are available from the corresponding authors upon reasonable request.

Acknowledgments

The authors acknowledge the Seed and Plant Improvement Institute (SPII), Agricultural Research, Education and Extension Organization (AREEO), Iran, for providing plant genetic material and supporting the research facilities.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vaezi, B.; Pour-Aboughadareh, A.; Mehraban, A.; Hossein-Pour, T.; Mohammadi, R.; Armion, M.; Dorri, M. The use of parametric and non-parametric measures for selecting stable and adapted barley lines. Arch. Agron. Soil Sci. 2018, 64, 597–611. [Google Scholar] [CrossRef]
  2. Fatemi, F.; Kianersi, F.; Pour-Aboughadareh, A.; Poczai, P.; Jadidi, O. Overview of identified genomic regions associated with various agronomic and physiological traits in barley under abiotic stresses. Appl. Sci. 2022, 12, 5189. [Google Scholar] [CrossRef]
  3. Warzecha, T.; Skrzypek, E.; Bocianowski, J.; Sutkowska, A. Impact of selected PSII parameters on barley DH lines biomass and yield elements. Agronomy 2021, 11, 1705. [Google Scholar] [CrossRef]
  4. FAOSTAT Statistical Database. Available online: http://www.fao.org/faostat/en/#home (accessed on 1 May 2022).
  5. Agricultural Statistics (2020–2021 Cropping Year); Crop Plants; Ministry of Agriculture-Jahad: Jahad, Iran, 2022; pp. 1–100.
  6. Ghazvini, H.; Bagherikia, S.; Pour-Aboughadareh, A.; Sharifalhossaini, M.; Razavi, S.A.; Mohammadi, S.; GhasemiKalkhoran, M.; Fathihafshejani, A.; Khakizade, G. GGE biplot analysis of promising barley lines in the cold regions of Iran. J. Crop Improv. 2022, 36, 461–547. [Google Scholar] [CrossRef]
  7. Ghazvini, H.; Pour-Aboughadareh, A.; Sharifalhosseini, M.; Razavi, S.A.; Mohammadi, S.; Ghasemi Kalkhoran, M.; Fathi Hafshejani, A.; Khakizadeh, G. Phenotypic stability analysis of barley promising lines in the cold regions of Iran. Crop Breed J. 2018, 8, 17–29. [Google Scholar]
  8. Ebem, E.C.; Afuape, S.O.; Chukwu, S.C.; Ubi, B.E. Genotype × environment interaction and stability analysis for root yield in sweet potato [Ipomoea batatas (L.) Lam]. Front. Agron. 2021, 3, 665564. [Google Scholar] [CrossRef]
  9. Jędzura, S.; Bocianowski, J.; Matysik, P. The AMMI model application to analyze the genotype–environmental interaction of spring wheat grain yield for the breeding program purposes. Cereal Res. Commun. 2023, 51, 197–205. [Google Scholar] [CrossRef]
  10. Bocianowski, J.; Księżak, J.; Nowosad, K. Genotype by environment interaction for seeds yield in pea (Pisum sativum L.) using additive main effects and multiplicative interaction model. Euphytica 2019, 215, 191. [Google Scholar] [CrossRef] [Green Version]
  11. Moussa, S.A.M.; Hala, A.A.E.; Nashwa, I.A.E. Stability study of sweet potato yield and its component characters under different environments by joint regression analysis. J. Hortic. Sci. Ornam. Plants 2011, 3, 43–54. [Google Scholar]
  12. Khalili, M.; Pour-Aboughadareh, A. Parametric and nonparametric measures for evaluating yield stability and adaptability in barley doubled haploid lines. J. Agric Sci. Technol. 2016, 18, 789–803. [Google Scholar]
  13. Vaezi, B.; Pour-Aboughadareh, A.; Mohammadi, R.; Mehraban, A.; Hossein-Pour, T.; Koohkan, E.; Ghasemi, S.; Moradkhani, H.; Siddique, K.H.M. Integrating different stability models to investigate genotype × environment interactions and identify stable and highyielding barley genotypes. Euphytica 2019, 215, 63. [Google Scholar] [CrossRef]
  14. Ahakpaz, F.; Abdi, H.; Neyestani, E.; Hesami, A.; Mohammadi, B.; Nader Mahmoudi, K.; Abedi-Asl, G.; Jazayeri Noshabadi, M.R.; 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, 10665. [Google Scholar] [CrossRef]
  15. Hilmarsson, H.S.; Rio, S.; Sanchez, J.I.Y. Genotype by environment interaction analysis of agronomic spring barley traits in Iceland using AMMI, Factorial Regression Model and Linear Mixed Model. Agronomy 2021, 11, 499. [Google Scholar] [CrossRef]
  16. Nykiel, M.; Gietler, M.; Fidler, J.; Graska, J.; Rybarczyk-Płońska, A.; Prabucka, B.; Muszyńska, E.; Bocianowski, J.; Labudda, M. Differential water deficit in leaves is a principal factor modifying barley response to drought stress. Int. J. Mol. Sci. 2022, 23, 15240. [Google Scholar] [CrossRef] [PubMed]
  17. Pour-Aboughadareh, A.; Yousefian, M.; Moradkhani, H.; Poczai, P.; Siddique, K.H.M. STABILITYSOFT: A new online program to calculate parametric and non-parametric stability statistics for crop traits. Appl. Plant Sci. 2019, 7, e01211. [Google Scholar] [CrossRef] [Green Version]
  18. Bocianowski, J.; Radkowski, A.; Nowosad, K.; Radkowska, I.; Zieliński, A. The impact of genotype-by-environment interaction on the dry matter yield and chemical composition in timothy (Phleum pratense L.) examined by using the additive main effects and multiplicative interaction model. Grass Forage Sci. 2021, 76, 463–484. [Google Scholar] [CrossRef]
  19. Gauch, H.G. Statistical Analysis of Regional Trials. AMMI Analysis of Factorial Design, 1st ed.; Elsevier: New York, NY, USA, 1992; p. 278. [Google Scholar]
  20. Yan, W.; Kang, M.S. GGE Biplot Analysis: A Graphical Tool for Breeders, Geneticists, and Agronomists; CRC Press: Boca Raton, FL, USA, 2002; p. 287. [Google Scholar]
  21. Gauch, H.G. Statistical analysis of yield trials by AMMI and GGE. Crop Sci. 2006, 46, 1488–1500. [Google Scholar] [CrossRef]
  22. Gauch, H.G. A simple protocol for AMMI analysis of yield trials. Crop Sci. 2013, 53, 1860–1869. [Google Scholar] [CrossRef]
  23. Dehghani, H.; Ebadi, A.; Yousefi, A. Biplot analysis of genotype by environment interaction for barley yield in Iran. Agron. J. 2006, 98, 388–393. [Google Scholar] [CrossRef]
  24. Ahmadi, J.; Vaezi, B.; Fotokian, M.H. Graphical analysis of multi-environment trials for barley yield using AMMI and GGE-biplot under rain-fed conditions. JPPB 2012, 2, 43–54. [Google Scholar]
  25. Mortazavian, S.M.M.; Nikkhah, H.R.; Hassani, F.A.; Hosseini, S.; Taheri, M.; Mahloohi, M. GGE biplot and AMMI analysis of yield performance of barley genotypes across different environments in Iran. J. Agric. Sci. Technol. 2014, 16, 609–622. [Google Scholar]
  26. Vaezi, B.; Pour-Aboughadareh, A.; Mohammadi, R.; Armion, M.; Mehraban, A.; Hossein-Pour, T.; Dorri, M. GGE biplot and AMMI analysis of barley yield performance in Iran. Cereal Res. Commun. 2017, 45, 500–511. [Google Scholar] [CrossRef] [Green Version]
  27. 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]
  28. Pour-Aboughadareh, A.; Barati, A.; Gholipoor, A.; Zali, H.; Marzooghian, A.; Koohkan, S.A.; Shahbazi-Homonloo, K.; Houseinpour, A. Deciphering genotype-by-environment interaction in barley genotypes using different adaptability and stability methods. J. Crop Sci. Biotechnol. 2023, 1–16. [Google Scholar] [CrossRef]
  29. Pour-Aboughadareh, A.; Khalili, M.; Poczai, P.; Olivoto, T. Stability indices to deciphering the genotype-by-environment interaction (GEI) effect: An applicable review for use in plant breeding programs. Plants 2022, 11, 414. [Google Scholar] [CrossRef]
  30. Ghaffari, M.; Gholizadeh, A.; Andarkhor, S.A.; Zareei Siahbidi, A.; Ahmadi, S.A.; Shariati, F.; Rezaeizad, A. Stability and genotype × environment analysis of oil yield of sunflower single cross hybrids in diverse environments of Iran. Euphytica 2021, 217, 187. [Google Scholar] [CrossRef]
  31. Tai, G.C. Genotypic stability analysis and its application to potato regional trials. Crop Sci. 1971, 11, 184–190. [Google Scholar] [CrossRef]
  32. Nowosad, K.; Tratwal, A.; Bocianowski, J. Genotype by environment interaction for grain yield in spring barley using additive main effects and multiplicative interaction model. Cereal Res. Commun. 2018, 46, 729–738. [Google Scholar] [CrossRef]
  33. Bocianowski, J.; Warzecha, T.; Nowosad, K.; Bathelt, R. Genotype by environment interaction using AMMI model and estimation of additive and epistasis gene effects for 1000-kernel weight in spring barley (Hordeum vulgare L.). J. Appl. Genet. 2019, 60, 127–135. [Google Scholar] [CrossRef] [Green Version]
  34. Bocianowski, J.; Tratwal, A.; Nowosad, K. Genotype by environment interaction for area under the disease-progress curve (AUDPC) value in spring barley using additive main effects and multiplicative interaction model. Australas. Plant Pathol. 2020, 49, 525–529. [Google Scholar] [CrossRef]
  35. Kandel, M.; Rijal, T.R.; Kandel, P. Evaluation and identification of stable and high yielding genotypes for varietal development in amaranthus (Amaranthus hypochondriacus L.) under hilly region of Nepal. J. Agric. Food Res. 2021, 5, 100158. [Google Scholar] [CrossRef]
  36. Wanga, M.A.; Shimelis, H.; Mashilo, J. Genotype by environment interaction of newly developed sorghum lines in Naminia. Euphytica 2021, 218, 147. [Google Scholar] [CrossRef]
  37. Lin, C.S.; Binns, M.R.; Lefkovitch, L.P. Stability analysis: Where do we stand? Crop Sci. 1986, 26, 894–900. [Google Scholar] [CrossRef] [Green Version]
  38. Yan, W. GGE biplot: A windows application for graphical analysis of multi-environment trial data and other types of two-way data. Agron. J. 2001, 93, 1111–1118. [Google Scholar] [CrossRef] [Green Version]
  39. Sharifi, P.; Aminpanah, H.; Erfani, R.; Mohaddesi, A.; Abbasian, A. Evaluation of genotype × environments interaction in rice based on AMMI model in Iran. Rice Sci. 2017, 24, 173–180. [Google Scholar] [CrossRef]
  40. Jamshidmoghaddam, M.; Pourdad, S.S. Genotype × environment interactions for seed yield in rainfed winter safflower (Carthamus tinctorius L.) multi-environment trials in Iran. Euphytica 2013, 190, 357–369. [Google Scholar] [CrossRef]
  41. Dehghani, M.R.; Majidi, M.M.; Mirlohi, A.; Saeidi, G. Integrating parametric and non-parametric measures to investigate genotype ×environment interactions in tall fescue. Euphytica 2016, 208, 583–596. [Google Scholar] [CrossRef]
  42. Bocianowski, J.; Prażak, R. Genotype by year interaction for selected quantitative traits in hybrid lines of Triticum aestivum L. with Aegilops kotschyi Boiss. and Ae. variabilis Eig. using the additive main effects and multiplicative interaction model. Euphytica 2022, 218, 11. [Google Scholar] [CrossRef]
  43. Zadoks, J.C.; Chang, T.T.; Konzak, C.F. A decimal code for the growth stages of cereals. Weed Res. 1974, 14, 415–421. [Google Scholar] [CrossRef]
  44. Limin, A.E.; Fowler, D.B. Cold hardiness expression in interspecific hybrids and amphiploids of the Triticeae. Genome 1988, 30, 261–265. [Google Scholar] [CrossRef]
  45. Kolde, R. Pheatmap: Pretty Heatmaps. R Package Version 1.0.12. Available online: https://cran.rproject.org/web/packages/pheatmap/index.html (accessed on 1 August 2022).
  46. Olivoto, T.; Lucio, A.D. Metan: An R package for multi environment trial analysis. Methods Ecol. Evol. 2020, 11, 783–789. [Google Scholar] [CrossRef] [Green Version]
  47. Husson, F.; Josse, J.; Le, S.; Mazet, J. FactoMineR: Multivariateexploratory Data Analysis and Data Mining. R Package Version 2.4. Available online: http://factominer.free.fr (accessed on 1 August 2022).
  48. Kassambara, A. Factoextra: Extract and Visualize the Results of Multivariate Data Analyses, R Package Version 1:7. Available online: https://CRAN.R-project.org/package=factoextra (accessed on 1 August 2022).
  49. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
Figure 1. The AMMI1 biplot showing the GEI effect for 20 barley genotypes across 16 test environments in the cold climate in Iran (A). Grouping pattern of estimated AMMI-based stability parameters based on the first two components (B). GY, ASTAB, ASI, ASV, AVAMGE, DA, DZ, EV, and FA indicate grain yield, AMMI-based stability parameter, AMMI stability index, AMMI stability value, average sum across environments of the absolute value of the GEI modeled using AMMI, distance between the IPCA points and the origin in space (A and Z), average of the squared eigenvector values, and fitted AMMI model, respectively.
Figure 1. The AMMI1 biplot showing the GEI effect for 20 barley genotypes across 16 test environments in the cold climate in Iran (A). Grouping pattern of estimated AMMI-based stability parameters based on the first two components (B). GY, ASTAB, ASI, ASV, AVAMGE, DA, DZ, EV, and FA indicate grain yield, AMMI-based stability parameter, AMMI stability index, AMMI stability value, average sum across environments of the absolute value of the GEI modeled using AMMI, distance between the IPCA points and the origin in space (A and Z), average of the squared eigenvector values, and fitted AMMI model, respectively.
Plants 12 02410 g001
Figure 2. (A) The ‘which-won-where’ view of the GGE biplot showing mega-environments and their winning genotypes for grain yield. (B) Biplot for the simultaneous selection of grain yield and stability of the investigated barley genotypes. (C) The ‘discriminating power and representativeness’ view of the GGE biplot. (D) Comparison of promising genotypes of barley against the ‘ideal’ genotype for grain yield and stability across the 16 test environments.
Figure 2. (A) The ‘which-won-where’ view of the GGE biplot showing mega-environments and their winning genotypes for grain yield. (B) Biplot for the simultaneous selection of grain yield and stability of the investigated barley genotypes. (C) The ‘discriminating power and representativeness’ view of the GGE biplot. (D) Comparison of promising genotypes of barley against the ‘ideal’ genotype for grain yield and stability across the 16 test environments.
Plants 12 02410 g002
Figure 3. Agro-climatic zone map showing the test locations in this study (https://commons.wikimedia.org/wiki/File:Iran-climate-map-es.svg accessed on 10 January 2023).
Figure 3. Agro-climatic zone map showing the test locations in this study (https://commons.wikimedia.org/wiki/File:Iran-climate-map-es.svg accessed on 10 January 2023).
Plants 12 02410 g003
Table 1. The response of the evaluated barley genotypes to the freezing test under field conditions at the Karaj research station during the 2019–2020 cropping season.
Table 1. The response of the evaluated barley genotypes to the freezing test under field conditions at the Karaj research station during the 2019–2020 cropping season.
Genotype CodeLT50 (°C)Thousand-Grain Weight (g)Grain Yield (Tons ha−1)
G1−840.208.71
G2−837.805.27
G3−839.007.01
G4−1038.807.62
G5−838.807.87
G6−7.541.306.20
G7−734.508.32
G8−842.006.90
G9−7.542.507.02
G10−835.905.94
G11−839.407.33
G12−642.606.56
G13−7.539.306.69
G14−840.407.06
G15−637.906.44
G16−6.545.207.10
G17−646.607.04
G18−8.535.907.21
G19−743.107.86
G20−841.508.51
LSD (0.05)3.686.462.02
Table 2. Result of the AMMI analysis of grain yield data for the evaluated barley genotypes across 16 test environments in cold regions in Iran.
Table 2. Result of the AMMI analysis of grain yield data for the evaluated barley genotypes across 16 test environments in cold regions in Iran.
SourcedfSSMSF-ValueProbability% TSS
Total9591972.32.057
Treatments31914004.3895.61**
Genotype (G)1998.45.186.62**4.99
Environment (E)15975.265.01121.56**49.44
Replication3296.53.0153.85**
GEI285326.41.1451.46**16.55
IPCA133752.2732.9**22.98
IPCA23159.41.9162.45**18.20
IPCA329441.5161.94**13.48
IPCA42737.51.3881.77**11.49
IPCA525301.21.53*9.19
Residuals14080.60.5760.74
Error608475.90.783
* and ** Significant at p < 0.05 and p < 0.01, respectively. TSS indicates the total sum of squares.
Table 3. The first four genotypes of barley selected using the AMMI model in each environment.
Table 3. The first four genotypes of barley selected using the AMMI model in each environment.
EnvironmentMeanIPCA1 Score1234
ARD17.86−0.4858G20G15G2G1
ARD26.212−0.4278G11G15G2G13
ARK17.558−0.0473G7G20G8G2
ARK25.579−0.4088G2G3G1G4
HAM17.289−0.3568G2G6G3G9
HAM210.014−0.5616G2G1G15G11
JOL16.244−0.4681G6G3G10G18
JOL26.838−0.3916G2G6G7G3
KAJ17.3060.0286G15G13G20G12
KAJ26.6140.2851G4G1G20G19
MAS16.850.3332G7G20G16G13
MAS26.3790.1469G18G1G20G4
MIN16.827−0.6028G2G15G20G5
MIN26.5460.8114G16G13G15G18
TAB17.461.2446G4G13G19G1
TAB28.3930.9007G1G18G16G20
ARD, ARK, HAM, MIN, MAS, KAJ, JOL, and TAB indicate Ardabil, Arak, Hamadan, Jolgeh-Rokh, Karaj, Mashhad, Miandoab, and Tabriz locations, respectively. Numbers 1 and 2 indicate the first and second cropping years, respectively.
Table 4. The results for the estimated AMMI-based stability parameter along with grain yield for the investigated barley genotypes.
Table 4. The results for the estimated AMMI-based stability parameter along with grain yield for the investigated barley genotypes.
Genotype CodeGYASTABASIASVAVAMGEDADZEVFA
G17.631.3300.1360.7457.1902.2600.5940.0715.110
G27.490.9900.1931.0607.2502.1400.4670.0434.570
G37.390.5100.1050.5764.5501.4700.3530.0252.160
G47.391.4400.1530.8416.6202.3800.6120.0755.670
G57.110.6200.1360.7455.2201.6700.3770.0282.770
G67.11.3100.2181.2007.7402.4100.5520.0615.810
G77.181.0100.0830.4575.7601.8900.5430.0593.570
G86.770.7900.0840.4645.0901.7000.4730.0442.880
G97.080.4800.0840.4634.5101.3800.3510.0251.910
G106.411.8400.1010.5568.0602.6400.7000.0986.960
G116.931.6600.1941.0708.1202.6100.6460.0836.820
G126.980.4100.0710.3923.9201.2800.3230.0211.630
G137.221.4400.2311.2708.2302.5700.5630.0636.600
G146.640.3400.0470.2583.6001.1200.3040.0181.250
G157.161.7700.1120.6138.8102.6000.6860.0946.760
G166.971.5400.2421.3308.9902.6600.5840.0687.070
G176.690.3800.0720.3953.9101.2200.3160.0201.490
G187.431.3100.1851.0206.9102.3600.5580.0625.580
G197.330.4200.1170.6434.2901.3900.3080.0191.920
G207.590.3800.0390.2143.5801.1900.3210.0211.430
GY, ASTAB, ASI, ASV, AVAMGE, DA, DZ, EV, and FA indicate grain yield, AMMI-based stability parameter, AMMI stability index, AMMI stability value, average sum across environments of the absolute value of GEI modeled using AMMI, distance between IPCA points and the origin in space (A and Z), average of the squared eigenvector value, and fitted AMMI model, respectively.
Table 5. The pedigree of the 18 evaluated, newly promised barley genotypes along with the reference genotype across 8 locations during 2 years in the cold agro-climate in Iran.
Table 5. The pedigree of the 18 evaluated, newly promised barley genotypes along with the reference genotype across 8 locations during 2 years in the cold agro-climate in Iran.
Genotype CodeLine/PedigreeSpike TypeGrowth Type
G1Jolge (Reference 1)Six-rowWinter
G2Bahman/3/Makouee//Zarjow/80-5151Six-rowWinter
G3Alger//CI10117/Choyo/3/Makouee/4/STB-12Six-rowWinter
G4Comp.Cr229//As46/Pro/3/Srs/4/Express/5/Goharan/6/GoharanSix-rowFacultative
G5Zarjow/80-5151//Makouee/3/MakoueeSix-rowWinter
G6Makouee//Zarjow/80-5151/3/BahmanSix-rowWinter
G7Radical/Birgit//Pamir-154/3/GoharanSix-rowFacultative
G8Cali92/Robust//ND16301Two-rowSpring
G9Radical/Birgit//Pamir-154/3/GoharanSix-rowFacultative
G10Yousef/4/82S:510/3/Arinar/Aths//DS 29Six-rowFacultative
G11Courlis/Rhn-03//KaroonSix-rowSpring
G12Mahtab/GoharanSix-rowSpring
G13Comp.Cr229//As46/Pro/3/Srs/4/Express/5/Goharan/6/GoharanSix-rowSpring
G14Pamir-147/Sonata/8/Alpha/Durra/7/P101/5/3896/1-15/3/3896/28//584/28/4/5050/6/TipperTwo-rowWinter
G15Courlis/Rhn-03//KaroonSix-rowWinter
G16Bda/Rhn-03//ICB-107766/3/YousefSix-rowFacultative
G17Sonata/8/Api/CM67//Hma-03/4/Cq/Cm//Apm/3/RM1508/5/Attiki/6/Aths/7/SP(6H)
/Apro//Ca1Mr/3/ROD586/Apm/4/Aths/9/Sararood
Two-rowWinter
G18Nadawa/Rhn-03//BirkaSix-rowSpring
G19Mahtab (Reference 2)Six-rowFacultative
G20Bahman/3/Alger//CI10117/ChoyoSix-rowWinter
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pour-Aboughadareh, A.; Ghazvini, H.; Jasemi, S.S.; Mohammadi, S.; Razavi, S.A.; Chaichi, M.; Ghasemi Kalkhoran, M.; Monirifar, H.; Tajali, H.; Fathihafshjani, A.; et al. Selection of High-Yielding and Stable Genotypes of Barley for the Cold Climate in Iran. Plants 2023, 12, 2410. https://doi.org/10.3390/plants12132410

AMA Style

Pour-Aboughadareh A, Ghazvini H, Jasemi SS, Mohammadi S, Razavi SA, Chaichi M, Ghasemi Kalkhoran M, Monirifar H, Tajali H, Fathihafshjani A, et al. Selection of High-Yielding and Stable Genotypes of Barley for the Cold Climate in Iran. Plants. 2023; 12(13):2410. https://doi.org/10.3390/plants12132410

Chicago/Turabian Style

Pour-Aboughadareh, Alireza, Habibollah Ghazvini, Seyed Shahriyar Jasemi, Solaiman Mohammadi, Sayed Alireza Razavi, Mehrdad Chaichi, Marefat Ghasemi Kalkhoran, Hassan Monirifar, Hamid Tajali, Asadollah Fathihafshjani, and et al. 2023. "Selection of High-Yielding and Stable Genotypes of Barley for the Cold Climate in Iran" Plants 12, no. 13: 2410. https://doi.org/10.3390/plants12132410

APA Style

Pour-Aboughadareh, A., Ghazvini, H., Jasemi, S. S., Mohammadi, S., Razavi, S. A., Chaichi, M., Ghasemi Kalkhoran, M., Monirifar, H., Tajali, H., Fathihafshjani, A., & Bocianowski, J. (2023). Selection of High-Yielding and Stable Genotypes of Barley for the Cold Climate in Iran. Plants, 12(13), 2410. https://doi.org/10.3390/plants12132410

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