Genotype × Environment Interaction of Yield and Grain Quality Traits of Maize Hybrids in Greece

: The interaction of genotype by the environment is very common in multi-environment trials of maize hybrids. This study evaluates the quantity and the quality of grain production and the stability of four maize genotypes in a ﬁeld experiment that was conducted in ﬁve different locations for two years. In order to make a reliable evaluation of the performance of genotypes in the environments, principal components analysis (PCA) was used to investigate the correlation of the yield, soil properties and quality characteristics, while the additive main effects and multiplicative interaction (AMMI) analysis detected the narrow adaptations of genotypes at speciﬁc mega-environments. For the yield, AMMI analysis indicated that a group of ﬁve environments (ENV1, ENV8, ENV6 ENV10 and ENV9) gave higher yields than the mean value and at the same time had low ﬁrst interaction principal components axis (IPC1) scores, indicating small interactions. Regarding protein and ﬁber contents, ENV1 and ENV2, gave the highest values and this could be attributed to the high concentration rates of nutrients like Mg, Ca and the soil texture (C). Speciﬁcally for the protein, the results of the analysis indicated that certain environment would provide more protein content, so in order to obtain higher grain protein, growers should grow in certain locations in order to improve the content of this quality characteristic, certain genotypes should be used in certain environments.


Introduction
Maize (Zea mays L.) is a crop of major importance for the nutrition of the Earth's population. Thus, there has been an urgent need to increase its yield and its quality. There are two main factors that have approximately the same influence on yield increase; improved management practices along with plant breeding have made an impact on this cause. It is notable that their interaction is the one that made such a huge progress to this matter that neither could do alone to this extent [1]. Maize yields have been increasing over the years globally and according to recent data have been doubled from 1961 to 2002 [1][2][3]. Even before the use of hybrids, farmers used to breed plants that seemed to fit their needs, with good adaptation at their specific environmental conditions, while maintaining their quality and morphological traits. Wherever hybrids have been adopted there has been an increase in the maize yield. Even though there was a tendency to select the high yield hybrids, the need for overall stability and dependability favors the selection of hybrids with stress tolerance. The main focus of the new hybrids now is to aim for a high and stable yield in both favorable and unfavorable growing conditions [1]. Indeed, as Rosegrant et al. [4] mention maize production may suffer a huge input reduction since the available water

Experimental Site and Design
The effect of GE interaction in relation to the final grain yield and the various quality characteristics were assessed. The field trials were conducted at 2019 and 2020 and a randomized complete block design was used with three replications. The experimental design had two main factors (genotypes and environments). The four maize genotypes used in the experimental fields were GEN1 (P0937, Pioneer Hi-Bred Hellas S.A.), GEN2 (DKC6050, K. and N. Efthymiadis S.A.), GEN3 (DKC6980, K. and N. Efthymiadis S.A.) and GEN4 (P2105, Pioneer Hi-Bred Hellas S.A.). P0937 hybrid belongs to the 500 FAO group with 120-125 days until maturity, it provides well-balanced plants, resistant to the northern leaf blight of maize, suitable especially for environments that traditionally produce good yield results. DKC6050 belongs to the 600 FAO group with 116-123 days until maturity. It produces average height plants with a very strong stem and root system and provides well-balanced plants. It has solid, filled until the top ears with 16-18 rows of kernels. It is considered to have great agronomic characteristics with high yield and good adaptation in many environments and with great quality grains suitable for human consumption. DKC6980 belongs to the 700 FAO group with 130-136 days until maturity. It produces high plants with a very strong stem and root system. It has big ears with 18-20 rows of kernels. It has stable yield even in high temperatures. P2105 is a hybrid that belongs to the 700 FAO group with 135-140 days until maturity, it provides well-balanced plants with a strong root system. It has a fast, lively growth and high dynamic production especially in early cultivated and well-drained fields. These four hybrids were selected because they are considered to have great characteristics in Greece's environmental variability.
The five locations were Giannouli, in the Prefecture of Larissa, Thessaly; Nea Tyroloi and Kamila, in the Prefecture of Serres, Central Macedonia; Kalamonas, in the Prefecture of Drama, Eastern Macedonia and Thrace and Koutso, in the Prefecture of Xanthi, Eastern Macedonia and Thrace, Greece. Each location-year combination created a different environment; thus 10 environments were used to evaluate the four genotypes.
The environments had different soil texture and variable microclimate conditions (Tables S1 and S2). Planting dates were between 22 and 25 March for the first year and 18 and 26 April for the second year and harvesting dates between 9 and 11 September for the first year and 19 and 21 September for the second year (Table S1). The plant density for all cultivars was 85,000 plants per hectare and the planting depth was 3.5 cm. The field plots (60 m 2 ) consisted of 4 rows with spacing between rows 0.75 m. The measurement of yield was made at a length of 12 m in the 2 middle rows of each plot to avoid the border effect. All five locations are traditionally used for maize production in Greece since they appear to be suitable for this crop.
All field management procedures were standard to ensure the avoidance of deficiencies and the balancing of soil nutrients in all environments. After grain harvesting, the measurement of moisture content was conducted with a portable humidity meter in order to calculate the production per hectare adjusted to 15% humidity. Afterward, grain samples were directed for the determination and measurement of the quality characteristics.

Soil Sample Analyses
At the soil samples, the elements Ca, Mg, K and Na were determined by atomic absorption spectrometry after extraction using BaCl 2 [26]. The measurements of Zn, Mn, Cu and Fe were conducted by atomic absorption spectrometry after extraction using DTPA [27]. The available B was determined using a spectrophotometer, using azomethine-H as the color (yellow) development reagent [28]. Total nitrogen was determined by the Kjeldahl method [29]. Organic carbon was determined with oxidization by K 2 Cr 2 O 7 [30]. Available phosphorus was determined after extraction with NaHCO 3 [31]. Cation exchange capacity was determined according to ISO 11260 [26]. Soil texture was determined using the method of Bouyoucos [32] and the soil taxonomy of USDA (1999). The moisture content was determined in a furnace at 105 • C for 24 h. The value of pH was measured with a pH-meter equipped with a glass electrode in the saturated paste extract.

Quality Characteristics of Harvested Corn Grains
After their harvest, the corn grains were dried in the shade following the farming practices. The moisture content across all cultivars varied from 10.85 to 16.04% with an average of 12.74%.
The color of the corn grains was measured using Minolta Colorimeter (CR-300, Minolta Company, Chuo-Ku, Osaka, Japan). The lightness or brightness of the samples was indicated by the L value where 0-100 represents darkness to lightness color. The index a indicates the redness or greenness of the corn grains, with a positive a value representing more red color. The index b value represents the degree of the yellow-blue color, with a positive b value illustrating more yellow color.
The texture analysis was carried out by HD-Plus texture analyzer (Stable Micro Systems Ltd., Godalming, UK) and the Texture Expert Exceed Software for the data analysis. The determination of the textural characteristics of corn grains was performed by a puncture probe of 5 mm diameter. Probe speeds of 1 mm/s during the test, 2 mm/s for the pretest and 10 mm/s for the post-test were used throughout the study. All the measurements were performed at 25 ± 1 • C and the hardness of the corn seeds was determined and expressed at N.
The corn grains were grinded by using a grinding mill for the determination of the moisture, the ash, the total protein and the total crude fiber content. Ash and crude fiber content of corn flours were determined according to AOAC Official Method 923.03 and 984.04 (Weende Method), respectively and recorded manually. Total protein content analysis of corn flours was conducted by applying the Kjeldahl method (IDF 2008), using a Kjeldahl rapid distillation unit (Protein Nitrogen Distiller DNP-1500-MP, Raypa Spain).

Statistical Analysis
A two-way analysis of variance (ANOVA) was used to evaluate the effect of genotype, environments and their GE interaction on the cultivation and quality characteristics of maize hybrids. The experimental data were analyzed using IBM SPSS software ver. 24 (IBM Corp., Armonk, N.Y., USA). The comparisons of means were calculated using the Duncan test at the 5% level of significance (p < 0.05). Multivariate analysis was conducted by means of principal component analysis (PCA) by using STATISTICA 7 (Statsoft Inc., Tulsa, OK, USA). Additive main effects and multiplicative interaction (AMMI) analysis was conducted by using AMMISOFT version 1.0 (Soil and Crop Sciences, Cornell University, Ithaca, NY, USA).

Grain Yield
Yield (kg/ha) of maize was statistically significantly influenced by genotype, environment and their interaction (Table S3). Genotype and environment effects on maize yield are presented in Table 1. GEN1 had mean yield that was significantly different in all the tested environments. Its highest yield value was recorded in ENV10 and ENV9 (19,288 ± 289 and 19,087 ± 471 kg/ha respectively) and the lowest in ENV2 (12,805 ± 1361 kg/ha). GEN2 presented similar results; the highest mean yield was recorded in ENV10 and ENV9 (18,244 ± 182 and 18,113 ± 86 kg/ha respectively) but the lowest mean yield in this case was recorded in ENV4 and ENV5 (13,617 ± 370 and 13,657 ± 51 kg/ha respectively). On the other hand, GEN3 presented its highest mean yield in ENV10 where it recorded 20,032 ± 179 kg/ha, and significantly lower results in ENV4 (13,220 ± 569 kg/ha). GEN4 had high yields in ENV6 (20,070 ± 346 kg/ha) and the significant lower value was presented in ENV4 (12,077 ± 166 kg/ha). As for the environments, in ENV4, ENV8 and ENV9, GEN1 presented significant higher mean yields. ENV4 had also a statistically significant performance in GEN2 along with ENV2 and ENV3. Environments ENV1, ENV6 and ENV7 showed significant high performance with GEN4. All the environments presented great results when GEN3 was used except ENV6 and ENV7, which had their lowest performance with this genotype.

Quality Characteristics of Harvested Corn Grains
The performance of genotypes on environments for quality characteristics (color parameters, texture, ash, protein and fiber content) of the harvested corn grains is presented in Table 1. The performance (color parameters) for corn grains varied across different environments. The lightness (L) across all cultivars ranged from 34.28 to 78.23, the yellow index b ranged from 34.47 to 58.90 and the red index a ranged from 3.61 to 8.80. For L, a and b color parameters, there was a significant difference between genotype (p < 0.001), environments (p < 0.001) and for their GE interaction, except for the interaction on the red index. The highest values of L color parameter of corn grains were obtained from the environments ENV2, ENV4, ENV6, ENV8 and ENV10 and the lowest from the environments ENV1, ENV3, ENV5 ENV7 and ENV9. According to the average of the tested environments, the highest values of L, a and b color parameters were achieved from genotype GEN2, and the lowest corresponding values were obtained from the genotypes GEN1 and GEN4 (Table 1).
The texture of the harvested corn grains was influenced by the environmental and genotype effects (p < 0.001) ( Table 1). Texture (hardness) of corn grains varied between 12.53 and 26.28 N across the cultivars and environments. The hardness of corn grains at ENV1 (21.38 N), ENV4 (18.75 N), ENV2 (18.64 N), ENV6 (18.63 N), ENV5 (18.33 N) and ENV10 (18.10 N) was relatively higher than the one at ENV3 (13.47 N) and ENV8 (13.38 N). According to the average of the tested environments, the highest values of hardness were achieved from genotypes GEN2 (20.72 N) and GEN1 (20.23 N), and the lowest corresponding values were obtained from the genotypes GEN3 (18.50 N) and GEN4 (17.20 N) ( Table 1).
The ash content of corn grains was significantly affected by genotype (p < 0.001), environments (p < 0.001) and their GE interaction (p < 0.001) ( Table 1). The ash content across all cultivars ranged from 0.73 to 2.40%. The highest ash content was found in corn grains obtained at the environments of ENV1 (1.51%) and ENV4 (1.50%) and the lowest at ENV8 (1.07%). The ash content of corn grains from GEN2 (1.39%) and GEN4 (1.38%) was significantly higher than the one from GEN1 (1.30%) and GEN3 (1.23%).
The genotype (p < 0.001), the environments (p < 0.001) and their GE interaction (p < 0.001) highly influenced the protein content of corn grains ( Table 1). The protein content across all genotypes varied from 6.13 to 8.96%. The highest protein content was found in corn grains obtained at the ENV1 (8.63%). The protein content of corn grains from GEN3 (7.59%) and GEN4 (7.45%) was significantly higher than the ones from the other genotypes. Protein content is a primary quality indicator for corn grains. Mut et al. [33] and Peterson et al. [34] reported that the grain protein content changed from 3 to 4% and 10.0 to 18.0%, respectively within different oat genotypes cultivated in different environments. The protein content was mainly affected by the environment rather than the genotype. This finding is in accordance with other scientific studies [18,35,36]. The protein content of corn grain is illustrated the quality of corn flour and is a desirable trait for the food industry.
The fiber content of corn grains was significantly influenced by genotype (p < 0.01), environments (p < 0.001) and their GE interaction (p < 0.001) ( Table 1). The fiber content across all cultivars ranged from 1.37 to 4.08%. The highest fiber content was found in corn grains obtained at the environments of ENV2 (3.65%) and ENV1 (3.54%) while the lowest at ENV7 (2.01%). The fiber content of corn grains from GEN2 (2.97%) and GEN1 (2.85%) was significantly higher than others. It was observed that the effect of environment on the fiber content of corn grains was stronger as compared to the genotype. This finding is in accordance with other scientific studies [18,37].

Correlation and Evaluation of the Yield, Soil Properties and Quality Characteristics vs. the Genotype and Environment on Maize Cultivation
To investigate the correlation of the yield, soil properties and quality characteristics by using four different genotypes of maize hybrids at ten different environmental conditions, principal components analysis (PCA) was used (Figure 1). Each point on the loading plot represented the contribution of a variable (yield, soil properties: clay, silt, sand, pH, organic matter, total nitrogen, CaCO 3 , K, Ca, Mg, P, Fe, Cu, Zn, Mn, B and quality characteristics: color, texture, moisture, ash, protein and fiber content) to the score, while each point on the score plot represented a tested sample. The first principal component (PC1) described 41.38% of the variation of extraction experiments, whereas the second principal component (PC2) 25.63% respectively, so that they contributed 67.01% of the total variation of extraction experiments. on PC2, while the sand, P and Zn had a positive effect on PC1 and the Mn, silt and yield had a positive effect on PC2. Furthermore, there are correlations between the Mg and fiber, clay and protein, pH and texture, and between L, B, K, Cu and CaCO3. Based on PCA score plot of the tested samples, four main groups of samples were noted. The groups are (a) 1a, 1b, 1c, 1d, 2a, 2b, 2c, 2d (b) 5a, 5b, 5c, 5d, 6a, 6b, 6c, 6d (c) 10a, 10b, 10c, 10d (d) 3a,  3b, 3c, 3d, 8a, 8b, 8c, 8d and (e) 7a, 7b, 7c, 7d, 9a, 9b, 9c, 9d. All five testing locations are suitable for maize production in Greece. However, they have different soil conditions and rainfall, influencing the yield and the quality characteristics of corn grains. The samples of group (a) confirmed that ENV1 and ENV2 were the most effective environments for all the tested hybrids, giving corn grains with the highest According to the PCA plot, the texture, fiber, protein, clay, Mg, Ca, pH and OM had a negative effect on PC1 and the total nitrogen, silt, CaCO 3 , Fe and K had a negative effect on PC2, while the sand, P and Zn had a positive effect on PC1 and the Mn, silt and yield had a positive effect on PC2. Furthermore, there are correlations between the Mg and fiber, clay and protein, pH and texture, and between L, B, K, Cu and CaCO 3 . Based on PCA score plot of the tested samples, four main groups of samples were noted. The groups are (a) 1a, 1b, 1c, 1d, 2a, 2b, 2c, 2d (b) 5a, 5b, 5c, 5d, 6a, 6b, 6c, 6d (c) 10a, 10b, 10c, 10d (d) 3a, 3b, 3c, 3d, 8a, 8b, 8c, 8d and (e) 7a, 7b, 7c, 7d, 9a, 9b, 9c, 9d.
All five testing locations are suitable for maize production in Greece. However, they have different soil conditions and rainfall, influencing the yield and the quality characteristics of corn grains. The samples of group (a) confirmed that ENV1 and ENV2 were the most effective environments for all the tested hybrids, giving corn grains with the highest protein and fiber content. These findings could be attributed to the enhanced soil fertility of ENV1 and ENV2 having the highest concentration of nutrients including Mg, Ca, clay and silt. Many studies demonstrated that the protein content was mostly affected by the environment, indicating its sensitivity to the environment [18,35,36,38] and that the soil nutrient supply affected positively the yield and the quality characteristics of the crop products [39].
The samples of group (b) indicated that ENV5 and ENV6 showed good soil conditions in terms of nutrients resulting in high yields. This finding is in accordance with other studies indicating that any higher nutrient uptake by the plant can result in higher yields [18,40]. The samples of group (c) showed that ENV10 resulted in grains with the maximum hardness concerning their texture essential soil indices such as pH and OM compared with the other environments. Moreover, the samples of group (d) showed that ENV3 and ENV8 had similar characteristics in terms of protein and fiber content growing on a Fe-, P-and N-and sand-rich environment. The samples of group (e) indicated that ENV7 and ENV9 had similar characteristics in terms of yield and quality properties growing on a Zn-rich environment. Concluding, the results depicted by principal components analysis are in agreement with those discussed above.

GE Interaction for Yield, Protein and Fiber Content
According to the ANOVA, genotypes (GEN), environments (ENV) and their interaction (G×E) gave statistically significant differences (p < 0.001) concerning the yield measurement. Moreover, the highest percentage of variation explained by ENV (68.89%), followed by the G×E (27.15%) effect, while GEN explained (3.95%) of the variation (Table S3). Figure 2 shows that GEN3 had the highest mean yield, followed by GEN1, GEN4 and GEN2. Among these, GEN4 was the hybrid with the lowest score of the first interaction principal components axis (IPC1). The great score values of IPC1 mean that these genotypes are adapted to certain environments [24]. As for the environments, ENV10 presented the highest mean yield (18,862 kg/ha) with an IPC1 score close to zero, indicating small interactions and ENV4 the lowest yield (13,227 kg/ha). ENV1, ENV6, ENV8 and ENV9 had IPC1 values close to zero and yield higher than the mean value (17,941,18,839,17,691 and 18,446 kg/ha respectively).   Based on AMMI1 model, GEN3 and GEN1 resulted in the highest narrow adaptations and these were the best adapted genotypes of the two mega-environments delineated ( Table 2). The first one consisted of ENV5, ENV2, ENV6, ENV1, ENV10 and ENV3 in which GEN3 was the best adapted genotype. The other one consisted of ENV8, ENV9, ENV4 and ENV7 with GEN1 presenting better results. Genotypes and environments had been listed according to their IPC1 order ( Table 3), resulting that top and bottom performances have the opposite GE pattern [33]. For instance, GEN3 had positive GE with the environments such as ENV5 and ENV2 and a negative GE with environments like ENV7 and ENV4.  The first two principal components in AMMI analysis were significant (p < 0.001), explaining 82.13% of GE (44.49% IPC1 and 37.64% IPC2) of interaction variation (Table S3). According to the biplot of the first (IPC1) and the second (IPC2) interaction principal components (Figure 3), GEN2 had a positive interaction with five out of ten environments (ENV3, ENV4, ENV8, ENV9 and ENV10). Generally, all genotypes were located far from the biplot origin and contribute to the G×E interaction for yield. AMMI analysis is widely used for the evaluation of maize hybrids yield in multienvironment field trials [20,24,25,[41][42][43], and grain yield of wheat varieties [22,36,[44][45][46], seed yield of oilseed rape [23], nutritional composition of sweet potato [38], yield of sugarcane [47] and yield of chickpea [48].
As for the protein content, according to the ANOVA, genotypes (G), environments (E) and their interaction (GxE) gave statistically significant values (p < 0.001). Moreover, the highest percentage of variation explained by ENV (66.90%) and G×E (29.32%) effects, while G explained the rest of variation (3.79%) (Table S3). Figure 4 shows that all genotypes had protein content percentages close to the mean value, with IPC1 values far from zero. ENV1 was the environment that presented the best results for the protein content (8.63%) of corn grains, while ENV9 presented the lowest one (6.54%). ENV10, ENV8 and ENV9 were the environments that presented the most stable results in terms of protein content. According to AMMI1 model, GEN3 and GEN2 resulted in the highest narrow adaptations delineating two mega-environments ( Table 2). The first one consisted of ENV7, ENV3, ENV9, ENV8 and ENV1, in which GEN3 was the best adapted genotype and the other one consisted of ENV10, ENV5, ENV4, ENV2 and ENV6, where GEN2 was the better suited genotype.
The first two principal components in AMMI analysis were significant (p < 0.001), explaining 83.93% of GE (48.61% IPC1 and 35.32% IPC2) of interaction variation (Table S3). According to the biplot of the first (IPC1) and the second (IPC2) interaction principal components ( Figure 5), GEN1 had a large positive interaction with ENV4 and ENV2, GEN2 had a large positive interaction with ENV6, while GEN4 and GEN3 had a positive interaction with ENV8 and ENV9 and ENV3 respectively. Likewise the yield, all genotypes are located far from the biplot origin and contribute to the G×E interaction for yield.
The first two principal components in AMMI analysis were significant (p < 0.001), explaining 82.13% of GE (44.49% IPC1 and 37.64% IPC2) of interaction variation (Table  S3). According to the biplot of the first (IPC1) and the second (IPC2) interaction principal components (Figure 3), GEN2 had a positive interaction with five out of ten environments (ENV3, ENV4, ENV8, ENV9 and ENV10). Generally, all genotypes were located far from the biplot origin and contribute to the G×E interaction for yield. AMMI analysis is widely used for the evaluation of maize hybrids yield in multienvironment field trials [20,24,25,[41][42][43], and grain yield of wheat varieties [22,36,[44][45][46], seed yield of oilseed rape [23], nutritional composition of sweet potato [38], yield of sugarcane [47] and yield of chickpea [48]. while G explained the rest of variation (3.79%) (Table S3). Figure 4 shows that all genotypes had protein content percentages close to the mean value, with IPC1 values far from zero. ENV1 was the environment that presented the best results for the protein content (8.63%) of corn grains, while ENV9 presented the lowest one (6.54%). ENV10, ENV8 and ENV9 were the environments that presented the most stable results in terms of protein content. According to AMMI1 model, GEN3 and GEN2 resulted in the highest narrow adaptations delineating two mega-environments ( Table 2). The first one consisted of ENV7, ENV3, ENV9, ENV8 and ENV1, in which GEN3 was the best adapted genotype and the other one consisted of ENV10, ENV5, ENV4, ENV2 and ENV6, where GEN2 was the better suited genotype.
The first two principal components in AMMI analysis were significant (p < 0.001), explaining 83.93% of GE (48.61% IPC1 and 35.32% IPC2) of interaction variation (Table  S3). According to the biplot of the first (IPC1) and the second (IPC2) interaction principal components ( Figure 5), GEN1 had a large positive interaction with ENV4 and ENV2, GEN2 had a large positive interaction with ENV6, while GEN4 and GEN3 had a positive interaction with ENV8 and ENV9 and ENV3 respectively. Likewise the yield, all genotypes are located far from the biplot origin and contribute to the G×E interaction for yield.
Regarding fiber content, according to the ANOVA, genotypes (G), environments (E) and their interaction (GxE) gave statistically significant values (p < 0.001). Moreover, the highest percentage of variation explained by ENV (56.43%) and G×E (39.15%) effects, while G explained the rest of variation (4.42%) (Table S3). Figure 6 presents that GEN2 had the highest value, followed by GEN1, GEN3 and GEN4. GEN1 had an IPC1 score close to zero, indicating small interactions. ENV6 scored an IPC1 value near 0, and at the same time had fiber content slightly higher than the mean value. ENV2 was the environment that presented the best results for the fiber content (3.65%) of corn grains, while ENV7 presented the lowest one (2.01%).  Regarding fiber content, according to the ANOVA, genotypes (G), environments (E) and their interaction (GxE) gave statistically significant values (p < 0.001). Moreover, the highest percentage of variation explained by ENV (56.43%) and G×E (39.15%) effects, while G explained the rest of variation (4.42%) (Table S3). Figure 6 presents that GEN2 had the highest value, followed by GEN1, GEN3 and GEN4. GEN1 had an IPC1 score close to zero, indicating small interactions. ENV6 scored an IPC1 value near 0, and at the same time had fiber content slightly higher than the mean value. ENV2 was the environment that presented the best results for the fiber content (3.65%) of corn grains, while ENV7 presented the lowest one (2.01%). According to AMMI1 model, GEN2 and GEN4 resulted in the highest narrow adaptations defining two mega-environments ( Table 2). The first one consisted of ENV10, ENV7, ENV5, ENV2, ENV6, ENV8, ENV9 and ENV4 in which GEN2 was the genotype that presented better results and the other one consisted of ENV1 and ENV3 that had GEN4 as the better suited genotype. The AMMI statistical model has been used to evaluate the effects of genotypes, environments and their interaction for quality characteristics, like iron and zinc concentrations in the grain of maize [49], protein and tryptophan in maize [50], nutritional composition (protein, β-carotene, iron, zinc, starch and sucrose) in sweet potato [38] and vitreousness, SDS sedimentation test, yellow pigment index, protein content and test weight in durum wheat [46].
The first two principal components in AMMI analysis were significant (p < 0.001), explaining 90.79% of GxE (55.20% IPC1 and 35.59% IPC2) of interaction variation (Table S3). According to the biplot of the first (IPC1) and the second (IPC2) interaction principal components (Figure 7), GEN2 had a large positive interaction with ENV7, GEN4 had a large positive interaction with ENV3, while GEN3 had a positive interaction with ENV10. According to AMMI1 model, GEN2 and GEN4 resulted in the highest narrow adaptations defining two mega-environments ( Table 2). The first one consisted of ENV10, ENV7, ENV5, ENV2, ENV6, ENV8, ENV9 and ENV4 in which GEN2 was the genotype that presented better results and the other one consisted of ENV1 and ENV3 that had GEN4 as the better suited genotype. The AMMI statistical model has been used to evaluate the effects of genotypes, environments and their interaction for quality characteristics, like iron and zinc concentrations in the grain of maize [49], protein and tryptophan in maize [50], nutritional composition (protein, β-carotene, iron, zinc, starch and sucrose) in

Conclusions
The evaluation of yield results of different genotypes, under different environmental conditions is a complicated issue, as a lot of parameters have to be considered, in order to lead to reliable results. In such experiments, often the AMMI analysis was used, which provides valuable information that contributes to the understanding of the G×E interaction. In this study, the GEN effect explained a low percentage of the variation and could not lead to the selection of a stable genotype for all environments. However, the results of the AMMI analysis contributed to dividing the region into mega-environments and introduce the most suitable genotype for every environment.
Concerning the yield, based on the AMMI1 model, GEN3 and GEN1 resulted in the highest narrow adaptations and these were the best adapted genotypes of the two megaenvironments delineated (the first consisted of ENV5, ENV2, ENV6, ENV1, ENV10 and ENV3 and the second consisted of ENV8, ENV9, ENV4 and ENV7, respectively). A group of five environments (ENV1, ENV8, ENV6 ENV10 and ENV9) gave higher yields than the mean value and at the same time had low IPC1 scores, which indicated that they gave high yield with all the genotypes used. Regarding the grain quality, GEN3 and GEN2 for the protein content and GEN2 and GEN4 for the fiber content resulted in the highest narrow adaptations delineating two mega-environments. Specifically for the protein, the results of the analysis indicated that in order to obtain higher protein content, certain genotypes should be used in certain environments. It is important to note that ENV1 and ENV2 (location Giannouli for years 2019 and 2020 respectively, gave the highest values concerning protein and fiber content. These findings could be attributed to the high concentration rates of nutrients like Mg, Ca and the soil texture (C).
The results of this study suggest that the target to increase the quantity and quality of grain yield of maize hybrids is a very challenging issue, due to the high G×E interaction, which can be implemented by exploiting positive GE interactions, by dividing the environment into mega-environments.

Conclusions
The evaluation of yield results of different genotypes, under different environmental conditions is a complicated issue, as a lot of parameters have to be considered, in order to lead to reliable results. In such experiments, often the AMMI analysis was used, which provides valuable information that contributes to the understanding of the G×E interaction. In this study, the GEN effect explained a low percentage of the variation and could not lead to the selection of a stable genotype for all environments. However, the results of the AMMI analysis contributed to dividing the region into mega-environments and introduce the most suitable genotype for every environment.
Concerning the yield, based on the AMMI1 model, GEN3 and GEN1 resulted in the highest narrow adaptations and these were the best adapted genotypes of the two mega-environments delineated (the first consisted of ENV5, ENV2, ENV6, ENV1, ENV10 and ENV3 and the second consisted of ENV8, ENV9, ENV4 and ENV7, respectively). A group of five environments (ENV1, ENV8, ENV6 ENV10 and ENV9) gave higher yields than the mean value and at the same time had low IPC1 scores, which indicated that they gave high yield with all the genotypes used. Regarding the grain quality, GEN3 and GEN2 for the protein content and GEN2 and GEN4 for the fiber content resulted in the highest narrow adaptations delineating two mega-environments. Specifically for the protein, the results of the analysis indicated that in order to obtain higher protein content, certain genotypes should be used in certain environments. It is important to note that ENV1 and ENV2 (location Giannouli for years 2019 and 2020 respectively, gave the highest values concerning protein and fiber content. These findings could be attributed to the high concentration rates of nutrients like Mg, Ca and the soil texture (C).
The results of this study suggest that the target to increase the quantity and quality of grain yield of maize hybrids is a very challenging issue, due to the high G×E interac-