Next Article in Journal
A Review on Epidemiological and Clinical Studies on Buckwheat Allergy
Next Article in Special Issue
Root-Derived Proteases as a Plant Tool to Access Soil Organic Nitrogen; Current Stage of Knowledge and Controversies
Previous Article in Journal
Coordination of Cyclic Electron Flow and Water–Water Cycle Facilitates Photoprotection under Fluctuating Light and Temperature Stress in the Epiphytic Orchid Dendrobium officinale
Previous Article in Special Issue
Nutrients Leaching in Response to Long-Term Fertigation and Broadcast Nitrogen in Blueberry Production
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multivariate Interaction Analysis of Winter Wheat Grown in Environment of Limited Soil Conditions

1
Biosense Institute, University of Novi Sad, 21000 Novi Sad, Serbia
2
Institute of Field and Vegetable Crops, 21000 Novi Sad, Serbia
3
Faculty of Agriculture, University of Novi Sad, 21000 Novi Sad, Serbia
4
Faculty of Economics in Subotica, University of Novi Sad, 21000 Novi Sad, Serbia
5
Institute of Applied Ecology Futura, University of Metropolitan, 21000 Belgrade, Serbia
6
Institute of Applied Sciences in Agriculture, 21000 Belgrade, Serbia
*
Author to whom correspondence should be addressed.
Plants 2021, 10(3), 604; https://doi.org/10.3390/plants10030604
Submission received: 18 February 2021 / Revised: 14 March 2021 / Accepted: 17 March 2021 / Published: 23 March 2021
(This article belongs to the Special Issue Plant-Soil Interactions: From Soil Fertility to Crop Growth)

Abstract

:
The less productive soils present one of the major problems in wheat production. Because of unfavorable conditions, halomorphic soils could be intensively utilized using ameliorative measures and by selecting suitable stress tolerant wheat genotypes. This study examined the responses of ten winter wheat cultivars on stressful conditions of halomorphic soil, solonetz type in Banat, Serbia. The wheat genotypes were grown in field trails of control and treatments with two soil amelioration levels using phosphor gypsum, in amounts of 25 and 50 tha−1. Across two vegetation seasons, phenotypic variability and genotype by environment interaction (GEI) for yield traits of wheat were studied. The additive main effects and multiplicative interaction (AMMI) models were used to study the GEI. AMMI analyses revealed significant genotype and environmental effects, as well as GEI effect. Analysis of GEI using the IPCA (Interaction Principal Components) analysis showed a statistical significance of the first two main components, IPCA1 and IPCA2 for yield, which jointly explained 70% of GEI variation. First source of variation IPCA1 explained 41.15% of the GEI for the grain weight per plant and 78.54% for the harvest index. The results revealed that wheat genotypes responded differently to stressful conditions and ameliorative measures.

1. Introduction

Today, more than ever, great attention has been paid to agro-ecological conditions and issues of climate changes [1,2]. Constant striving to intensify agricultural production, excessive using of chemicals, fertilizers, as well as massive irrigation has led to degradation of arable land, pollution of the environment and significant climate change [3]. At the same time, the less productive soils present one of the major problems in agricultural production. Halomorphic soils, characterized as unfavorable soil for agricultural production, occupy significant areas of the world. The properties of halomorphic soils are controlled by the presence of either soluble salts, exchangeable sodium or both [4,5]. Soil degradation resulting from soil salinity and/or sodicity is recognized as a major problem of soil productivity and quality under arid and semiarid climates [6,7]. According to the current soil national classification, halomorphic soils are divided into a class of saline soils, characterized by high salt content and a class of alkaline soils characterized by a high content of adsorbed sodium in the second (Bt) horizon, which includes solonetz type of soil [8]. Worldwide, solonetz cover about 135 million hectares, with 20 million hectares in Europe. Major solonetz areas are in Ukraine, Russia, Hungary, Bulgaria, Rumania, China, Canada, less humid parts of South America, the south western and north-central part of the United States, South Africa and Australia [9]. In Vojvodina Province (North Serbia), beside of about a million hectares of high productive soils, there are about 80,000 ha less productive soils of halomorphic, solonetz soil type, mainly in Banat region [4]. These soils are most commonly used as extensive pasture and could be intensively utilized using ameliorative measures, which include incorporation of phosphogypsum in soil or by selecting suitable crops genotypes [10]. Long-term results indicated that the application of phosphogypsum in solonetz soil showed positive changes in the content and qualitative composition of cations, increases Ca2+ content in soil and decreases the alkaline reaction, thus improving the content and stability of the structural aggregates, water and air regimens, thermal and nutrient status, as well as the effective fertility of solonetz soils [11,12]. Phosphogypsum, produced by the phosphate fertilizer industry, as a by-product of wet acid production of phosphoric acid from rock phosphate, contains more than 92% calcium sulfate and presents suitable source of calcium, which may be used as a soil conditioner for sodic, solonets and solonetzic soils [13]. Nayak et al. [14] found that with the increasing amounts of phosphogypsum applied on agricultural soil without vegetation, pH was reduced from 7.9 in control to 5.1 in treatment with 20% phosphogypsum. Aside from its benefits, since phosphate ores have increased natural radioactivity, radionuclides are a special type of impurities that do not significantly affect the quality of phosphogypsum, but significantly affect the environment [15]. Hence, there is certain concern that the application of phosphogypsum to agricultural lands may results in plant uptake of radionuclides, fluoride and trace elements [16] Therefore, since phosphogypsum is formed as a by-product in the production of phosphoric acid in such large quantities, its storage is often an environmental problem [17]. However, agricultural production is greatly influenced by soil extremes and variability. The soil degradation process can significantly reduce plant diversity and agricultural yield, land productivity and value in arid and semiarid climate regions [18]. Adverse effects of high Na concentration in the B horizon on management practices and on seedbed preparation has been observed by Cairns [19], while adverse effect to the uptake of plant nutrients by the roots was observed by Peters [20] and Miller and Brierley [11]. Despite wheat being considered as basically grass, durable and more resistant to limiting factors in agricultural production than many other crops, the considerable decreases in yield and yield related traits have been observed in wheat production under stress conditions of halomorphic soils [5,21,22]. At the same time, global demand for wheat is estimated to increase by 2050, as a consequence of growing dietary requirements of an increasing world population [23]. Keeping this in view, the wheat breeders have to share an urgent need to increase wheat grain yield potential by developing new wheat varieties with desirable genetic make-up [24,25]. Considering that much more attention needs to be given not only to the quantity, one of the greatest challenges to wheat breeders is to obtain a wheat genotypes with high mean yield and wide adaptation to the various environments, even in less productive soils, such as a solonetz soil type. Wide adaptation is defined as the genotypes ability to produce relatively high yields consistently across diverse agricultural environments in a growing region [26,27,28]. An assessment of the stability of cultivars yield provides valuable information about their behavior in specific environments [29,30]. The complexity of grain yield and yield related traits is the result of different genotype reactions to changing environmental conditions during plant development [30]. Grain yield is a complex polygenic, quantitative trait, whose expression is the result of the genotype, environment and genotype by environment interaction (GEI). GEI reflects the different responses of the genotypes to environment conditions and distinctly indicates that the best genotypes in a specific environment could not be the best for others. Hence, GEI cannot represent all genetic potential environmental conditions, which makes difficult the recommendation of genotypes by the breeder [31,32,33]. For accessing and better understanding the GEI effects, several statistical methods have been developed. The additive main effects and multiplicative interaction (AMMI) model has been one of the most widely used statistical tools in the analysis of GEI [30,34,35,36,37]. The AMMI model combines firstly the analysis of variance (ANOVA) to partition the variation into genotype main effects (G), environment main effects (E) and genotype by environment interaction (GEI) effects and then it applies principal components analysis (PCA) to GEI in a single analysis [37]. The two main purposes of AMMI analysis of a yield trial’s treatment design are: (i) Understanding complex GEI, which includes delineating mega-environments and selecting genotypes to exploit narrow adaptations; and (ii) increasing accuracy to improve recommendations, repeatability, selections and genetic gains 1 [33]. Therefore, accurate determination of the source of variation, their effects on yield and yield related traits and selection wheat genotypes that have been adapted to limited conditions of solonetz soil is of great importance. In this respect, the aims of the present research were: (i) To determine the influence of a genotype, environment and their interaction on grain yield and yield related traits in different wheat varieties; and (ii) to evaluate stability of the traits using the AMMI model in conditions of different levels of melioration on solonetz soil type. These results could provide the identification of suitable and stable wheat genotypes, which can be successfully used in wheat production on less productive, solonetz soils.

2. Materials and Methods

2.1. Field Trial

The present study was carried out at the experimental trial field on halomorphic soil of solonetz on location Kumane (45.539° N, 20.228° E), in Banat region, in Vojvodina Province (Serbia), during two consecutive vegetation seasons of 2004/2005 and 2005/2006. The experimental material in the study was comprised of ten winter wheat cultivars (Triticum aestivum L.), namely, Mina (G1), Sofija (G2), Tiha (G3), Anastazija (G4), Nevesinjka (G5), Evropa 90 (G6), NSR-5 (G7), Dragana (G8), Ljiljana (G9) and Simonida (G10). The wheat cultivars used in the study were released by the Institute of Field and Vegetable Crops, Novi Sad, in Serbia. All cultivars are agronomically suitable for production in agroecological conditions of Serbia and other surrounding countries. The experimental trial was set up on solonetz type of soil, according to the completely randomized block design (RBCD), with three amelioration treatment and three replications of each treatment. The experimental plots consisted of eight 12.5 cm-spaced rows, 2 m in length. Plant spacing was realized by implementing a plant density of 500 seeds per m−2. Disease, weed and pest control on each plot was performed by using standard cultivation practice, while the NPK (15:15:15) application was split, 50 kg for each treatment. Planting in both growing seasons was completed by the second decade of October, while harvest was ended in the first decade of July. Since the experiment was set on solonetz type of soil, in addition to the results that were analyzed were the results of the soil in two levels maintenance, being 25 and 50 tha−1 phosphogypsum. The objective of applying reclamation measures to solonetz soils is to change their cation content which makes changes in chemical reactions, because increasing the calcium cation level and reducing the alkaline reaction improves the stability of structural aggregates, thermal, water, air properties and nutrient status [38]. The first treatment in the trial was soil without amelioration (control), the second treatment was amelioration using phosphor-gypsum in the amount of 25 tha−1 and the third treatment was amelioration using phosphor-gypsum in the amount of 50 tha−1. In the field trial, each treatment in one growing season was considered as a special agro-environment. This produced six different agro-ecological environments, which were equal in agrotechnical terms, but in different treatments of phosphogypsum melioration of soil repair. The six analyzed environments were labeled as follows: E1 represents control—solonetz soil without melioration in the season 2004/2005, E2 represents solonetz soil with melio-ration of 25 tha−1 of phosphogypsum in the season 2004/2005, E3 represents solonetz soil with melioration of 50 tha−1 of phosphogypsum in the season 2004/2005, E4 represents control—solonetz soil without melioration in the season 2005/2006, E5 represents solonetz soil with melioration of 25 tha−1 of phospho gypsum in the season 2005/2006, E6 represents solonetz soil with melioration of 50 tha−1 of phosphor gypsum in the season 2005/2006, Table 1.
At the stage of full maturity, ten average wheat plants from each replication of each wheat cultivar plot separately were selected and yield traits, such as a plant height (cm), grain weight per plant (g) and harvest index (%), were analyzed. The main sample consisted of 10 plants per replication, while grain yield, expressed in m−2, at 13% moisture, were obtained by harvesting each wheat cultivar plot separately.

2.2. Soil Properties

The experiment was set up at the of halomorphic soil, solonetz type in Banat, Serbia. This type of soil is characterized by unfavorable physical and chemical properties, caused by high content of clay and sodium in the Bt,na horizon. The presence of adsorbed sodium in the Bt horizon causes peptization of colloids and a high alkaline reaction of soil [38]. Table 2 shows the soil properties until a depth of 60 cm of typical solonetz, from which soil samples were taken to establish the experiment.
Solonetz soil is characterized by a lower concentration of salt and possesses a weakly acidic reaction which, with increasing depth, turns into an alkaline reaction. The content of CaCO3 is missing in the soil surface layer and its content increases with depth of soil. The content of humus and nitrogen in the surface layer is high, but drastically decreases with increasing depth. The soil is poor in phosphorus, while the potassium content is optimal. Permeability of solonetz soil is reduced and water penetration into the deeper horizons is possible only through the soil cracks in the dry period [39]. The solonetz type of soil is unsuitable for agricultural production due to the poor water and air regime, the possibility of stagnation and compaction during dry periods.

2.3. Meteorological Conditions

The meteorological conditions of vegetation seasons during the trial were obtained from automatic weather stations of the Republic Hydro-meteorological service of Serbia, situated in Bečej, near to the field trials. The weather conditions in the first vegetation season of 2004/2005 were more favorable for wheat cultivation than climatic conditions of 2005/2006 growing season, Figure 1a,b.
The most significant characteristic of weather conditions in the first vegetation season was unusually high amounts of precipitation. Climatic conditions for wheat growth and development in this period can be characterized as good. Wheat planting was done in the optimal time and plants entered prepared for the conditions of the winter period. Although the winter was cold, the snow cover protected the winter crops from frosts and provided excellent thermal insulation to the crops. At the end of the winter, partial melting of snow lead to an increase in soil moisture reserves. The spring frosts did not cause damage. May and June were characterized by average weather conditions.
Climatic conditions in the second vegetation season were significantly more variable compared to the conditions in the previous season. The second vegetation season had unusually frequent very large and rapid changes in the values of meteorological elements that most affect the plant growth and wheat development stages. On several occasions, there was a change of periods with extremely cold and extremely warm weather, periods with floods and droughts whose length could be measured in weeks or months. Wheat planting was done in the optimal time, but at the time of sowing, due to the deficit of precipitation, there was a strong soil drought and wheat germination was difficult. November saw an extreme precipitation deficit, in relation to the multi-year average. During the winter, there were extremely low temperatures and long-lasting frosts which were unfavorable for wheat due to the lack of snow cover. Late winter and early spring were characterized by heavy rainfall, floods and long-term wetlands. April, May and the middle of June were significantly colder with high precipitation, while in the second half of June, the average daily air temperatures significantly exceeded the multi-year average value. It was characterized by extreme heat and critical wheat phases in the grain development took place in extreme drought conditions. Air temperature values were generally higher than optimal at the time of grain filling, while soil moisture reserves were rapidly depleted due to intensive evapotranspiration, Figure 1a,b.

2.4. Statistical Analyses

Genotype by environment interaction (GEI) was estimated using AMMI (additive main effects and multiplicative interaction) analysis developed by Zobel et al. [33]. The AMMI model Zobel et al. [33] is presented as the following formula:
Yger = μ + αg + βe + Σλnξgnηen + Θge + εger,
where yge is the mean of yield or other observed trait for genotype g in the environment e, μ—grand mean, αg—genotypic mean deviations, βe—environmental mean deviations, n—number of PCA axis retained in the adjusted model, λn—eigenvalue of the PCA axis n, ξgn—genotype score for PCA axis n, ηen—score eigenvector for PCA axis n, Θge—residual and εger—experimental error.
The AMMI model incorporates analysis of variance (ANOVA) and principal components analysis (PCA) into a single statistical model [40]. In the AMMI model, the ANOVA separate additive effects from the interaction, while additional GEI analysis can be carried out by principle component analysis (PCA). The biplot graphic representation shows both the main and interaction effects for genotypes and environments simultaneously and provides analysis of the G × E interaction [22,36,40,41]. The IPCA1 score of a genotype in the AMMI analysis was used as an indicator of the stability of a genotype over environments [12,42]. Zero IPCA value indicates highest stability, while an IPCA value long distance from zero indicates genotype instability [43,44]. The data processing was performed in GenStat 9th Edition statistical software package (trial version) VSN International Ltd. (www.vsn-intl.com, accessed on 11 December 2020). [45].

3. Results and Discussion

Grain yield is a complex super trait consisting of many individual yield components and largely dependent on the genetic potential, which considerably varies primarily as a result of environmental conditions during the growing season [36,46,47]. Changes in certain morphological traits and yield stability have led to an increase in the genetic potential of the genotypes, most probably through the increased tolerance to biotic and abiotic stress factors [48]. Contribution of each individual yield component could vary in the different environmental conditions; testing in one environment could provide only limited information, while testing in different environments provides additional useful information (e.g., a GEI component can be estimated) [49]. Since that interaction represents the heterogeneity of agro-ecological conditions, it is important to identify genotypes that are very efficient, which means that they are productive and stable genotypes, but also genotypes with favorable response to soil amelioration. Estimation of GEI may prove valuable in recommending suitable wheat genotypes, as well as change of agricultural practice.

3.1. Plant Height

The presented results showed that the mean values of plant height of observed wheat genotypes in the trial ranged between 62.1 (E3) and 63.5 cm (E1 and E2) within different environments in the first vegetation season. In the second vegetation season, the mean values of plant height of wheat ranged between 42.8 (E4) and 58.6 cm (E5) within different environments. Between genotypes, the mean values of plant height ranged from 52.4 (G3) to 61.5 cm (G6) within different genotypes (Table 3). In the first vegetation season, at the control variant (E1), solonetz with no amelioration applied, the greatest overall mean values for plant height were denoted for wheat varieties G6 (68.6 cm), G4 (66.9 cm) and G5 (65.6 cm). At the same variant, the lowest mean values of the plant height were observed for genotypes G10 (59.3 cm), G2 (59.9 cm) and G3 (60.2 cm). In the second vegetation season, at the control variant (E4), the greatest overall mean values were denoted for genotypes G6 (47.8 cm), G5 (47.1 cm) and G4 (47.0 cm). The lowest values of plant height were denoted for wheat genotypes G10 (35.6 cm), G2 (38.7 cm), G1 (40.1 cm) and G3 (40.1 cm), Table 3.
In the second environment, with soil amelioration of 25 tha−1 phosphor gypsum applied, during the first vegetation season (E2), the greatest overall mean values for plant height were denoted for wheat genotypes G6 (68.5 cm), G4 (66.4 cm), G8 (65.2 cm) and G7 (64.9 cm). The lowest values of stem height were denoted for wheat genotypes G3 (59.8 cm), G2 (60.0 cm) and G10 (60.9 cm), Table 3.
In the third environment during the first vegetation season at the amelioration of 50 tha−1 phosphor gypsum applied (E3), the greatest overall mean values for plant height were observed for wheat varieties G6 (67.1 cm), G8 (64.7 cm), G4 (64.5 cm) and G7 (63.0 cm), while the smallest values showed for genotypes G3 (57.9 cm), G2 (58.9 cm) and G1 (60.3 cm). In the second vegetation season, with same level of amelioration (E6), the greatest overall mean values of grain yield were denoted for wheat varieties G6 (53.27 cm), G8 (52.63 cm), G10 (51.85 cm) and G4 (49.47 cm), while the smallest values showed for G3 (43.18 cm), G5 (43.99 cm) and G2 (45.69 cm), Table 3.
The results of this study showed that significant plant height variation was noticed due to different environment conditions, but also depended on genotype, as well as on vegetation season. This result is expected since plant height is one of the important yield components and is considered as a quantitative and variable trait whose expression highly depends on the environmental factors [36,50,51]. This was confirmed by high values of variation which ranged from 39.6% to 86.8% (Table 3).
According to the results, a higher sensitivity of plant height to the meteorological conditions was observed in the second vegetation season, since the overall plant height average was the lowest in the control variant. On the other side, more favorable conditions in the first growing season exhibited a more significant impact on reducing the differences between the mean averages per treatment. In terms of the second vegetation season, higher sensitivity of plant height of wheat was mostly caused by drought conditions after planting, winter frost and a less amount of rainfall in the first part of the growing season, and the effects of ameliorative measures were the most evident.
Apart from climatic conditions, tested wheat genotypes on solonetz soil reacted differently on different levels of soil repair. The phenotypic expressions of plant height of certain genotypes were constant in various environments, whereas some other wheat genotypes exposed significant variation over different soil amelioration treatments. In this sense, genotypes G2, G8 and G10 expressed a higher sensitivity to the treatment and reacted favorably on ameliorative measures from 25 (G2 and G8) to 50 tha−1 (G10) in the first season, while in the second season all tested genotypes reacted favorably on both amounts of applied measures. In general, vegetation periods of 2004/2005 provided generally higher mean values of the plant height than vegetation periods of 2005/2006. Varieties responded to ameliorated measures enhancing the plant height in favorable agro-ecological conditions. In addition, more favorable weather conditions of the first vegetation season exhibited significant impact on reducing the differences between the averages per treatment. Hence, the changing weather conditions influenced the effect of amelioration on the soil and the observed variation in the plant height of wheat.
The AMMI analysis of variance showed significant effects for genotypes and environments (Table 4).
The combined ANOVA showed that the phenotypic expression of plant height was significantly influenced by environmental variations, because the significant variance at the 1% level explained 80.14% of the total variation, while genotype contributed 8.64% of the total variation of the experiment (Table 4).
A large sum of squares for the environments indicated that the environments were diverse, with large differences among environmental means causing variation in the plant height. Differences between the growing seasons and diversity of treatments caused a considerable sum of squares for environmental factors in the total variation, indicating that these factors were the most responsible for the variation of the plant height of wheat. GEI contributed 11.23% of the total sum of squares and expressed no significant mean square, indicating that no cross interaction was expressed for the plant height. Large differences among environments and vegetation seasons caused a high sum of square of environmental factors in the overall variation of the trial, and led to the conclusion that they are the most responsible for variations in plant height of wheat. Such results are in agreement with results reported by Popović et al. [36,52] as well as with results by Branković-Radojčić et al. [53], who stated that, in multi-environment trials, often environments explain about 80% of the total variation, while the genotype (G) and GEI share is about 10%. This result was expected since the plant height of wheat is considered as a quantitative and highly variable trait whose expression highly depends on the environmental factors [36,54]. Additional analysis of the GEI using the PCA (interaction principal components) analysis confirmed there were no statistical significance of the first two principal components, IPCA 1 and IPCA 2 (Table 4). Separately, IPCA1 and IPCA2 participated in the GEI variation, with 36.6% and 29.1%, respectively. Both without a statistically significant effect on the GE interaction variation. These two main components explained, jointly, 65.6% of the variation of the genotype by environment interaction (Table 4).
In the AMMI biplot, the genotype and environment main effects for plant height of ten winter wheat genotypes are presented on the x-axis, while the IPCA1 (interaction principal component axis 1) scores are on the y-axis. The vertical line is the grand mean for trait plant height and the horizontal line (y-ordinate) represents the IPCA1 value of zero (Figure 2a). Observing trends in PCA1 scores and based on trends in two-dimensional graphical representation of genotypes and environments, wheat varieties were not differing in relation to main effect, but certain differences could be denoted in multivariate effect.
According to the biplot (Figure 2a) and in terms of average values, it could be noticed that agro-ecological environment E3 was at the level of the experimental overall average. Small distance environmental points from the origin (zero point) were also denoted for environment E2, which indicated that environments E2 and E3 were assessed as the most favorable to achieve a stable reaction to plant height. Environment E1, which correspondent to control variant, was further from the origin, while environments E4, E5 and E6 had the highest interaction values and contributed a lot to the GEI, which determines them the least suitable for stable establishment plant height of wheat. According to the arrangement of environmental points E1, E2, E3 and E5 in relation to the mean value of trait, it can be concluded that the genotypes achieved higher average values of plant height in these environments compared to E4 and E6 points. This result does not favor the E5 environment, as well as E4 and E6, for obtaining higher values of plant height, given the high values of interaction indicate the poor stability of this trait. In general, distribution environment points on the graph revealed that the more stable reaction was achieved in the first vegetation season (Figure 2a).
Based on the graphic presentation of the interaction of genotypes and environments, ten winter wheat genotypes differed in genotype by environment interaction. Wheat genotype G6 (Evropa) expressed the most stable reaction, showing a PCA1 score quite near to 0. Genotypes G1, G2, G3 and G9 showed values close to zero and they had small interaction effects, small contribution to the GE interaction and fitted well into the additive model. Genotypes G4, G7 and G9 showed higher distance from origin, while genotypes G10 and G5 were furthest from the PCA axis, thus they proved to be the most unstable and contributed a lot to the GE interaction captured by the first axis IPCA1. These two genotypes showed the highest interaction at the level of the whole experiment. Genotype G6 appeared to be better adapted to the E3 and E2 environments, which corresponded to soil treatments of 25 and 50 tha−1 phosphogypsum, having a mean value of plant height above the overall mean. Genotypes G9, G7 and G4 expressed also a positive effect in E3 and E2 environments, which correspond to environments with 25 and 50 tha−1 phosphogypsum applied, keeping their average quite above the level of overall mean. Genotype G8 appeared to be better adapted to conditions of E3 and E5, which corresponded to soil treatments of 25 and 50 tha−1 phosphogypsum, while genotypes G2 and G1 were better adapted to conditions of E6, which corresponded to soil treatment of 50 tha−1 phosphogypsum. Genotype G3 appeared to be adapted to conditions of E2, E3 and E1 environments, keeping its average below the overall mean (Figure 2a). According to Gauch [32], the genotypes at the top of the graph table have positive G × E interactions with environments and those genotypes at the bottom have negative G × E interactions with environments. Therefore, genotypes at the top of graph (G3, G4, G5, G6, G7 and G9) showed positive G × E interactions with environments, while genotypes which were at the bottom (G1, G2, G8 and G10) showed negative G × E interactions with environments (Figure 2a).
Obtained results are in accordance with results of AMMI selection of wheat genotypes across the testing environments (Table 5).
Based on the tested genotypes, genotype G6 could be considered as the most promising, since it ranked within the first four AMMI selections in all environments and showed consistent performance, implying it having the best adaptability for all soil conditions. In soil environment E1 (control, soil without amelioration) in the first vegetation season, genotype G6 was first ranked, followed by genotypes G4, G8 and G7. In soil environment E4 (control, soil without amelioration) in the second vegetation season, genotype G6 was first ranked, followed by genotypes G4, G5 and G7. This result indicates that these genotypes have a positive response and stable reaction in less favorable conditions in the control environment. In soil environment E2, which presents soil treatment with 25 tha−1 phosphogypsum in the first season, genotype G6 was first ranked, followed by genotypes G8, G4 and G7. In soil environment E5 (treatment with 25 tha−1 phosphogypsum) in the second season, genotype G6 was first ranked, followed by genotypes G8, G10 and G4. This indicates that these genotypes have positive response on the soil melioration treatment with 25 tha−1 phosphogypsum. In soil environment E3, which corresponds to soil treatment with 50 tha−1 phosphogypsum, in the first season, genotype G6 was first ranked, followed by genotypes G5, G4 and G7. In soil environment E6, which corresponds to soil treatment with 50 tha−1 phosphogypsum, in the second season, genotype G6 was first ranked, followed by genotypes G8, G10 and G4. The results confirm that it is possible to successfully advise the right genotype for all environments, or specific genotypes for specific environments, through AMMI evaluation (Table 5).
However, genotypes which maintain their height in stress environments could be desirable in certain limited conditions, because greater expression of plant height has been frequently attributed to high capacity of stems to accumulate sufficient stem reserves for the partitioning to grain.

3.2. Grain Weight per Plant

Grain weight as a quantitative trait, together with several yield related traits contributes to the formation of total grain yield of wheat. Each yield component is controlled by genes and affected by environmental factors, as well as by the interaction of genetic and environmental factors [55]. Grain weight per plant directly reflects the efficient use of nutrients and their translocation into generative parts of plant [56]. Grain weight per spike, as grain weight per plant, as the last yield components, are the final in the development of many components that occur in the early ontogenic stages [56].
The presented results showed that the mean values of the grain weight per plant of observed wheat genotypes ranged between 6.29 (E3) and 6.51 g (E1) within different environments in the first vegetation season. In the second vegetation season, the mean values of grain yield per plant of wheat ranged between 4.04 (E4) and 6.40 g (E5) within different environments. Between genotypes, the mean values of grain yield per plant ranged from 5.04 (G3) to 6.77 g (G6) within different genotypes (Table 6).
At the control variant (E1), with no amelioration applied, the greatest overall mean value for grain yield per plant was denoted for wheat variety G6 in both years of the study, having values of 7.53 g in the first season and 5.11 g in the second season. In both seasons the same trend was observed between genotypes and the higher values were observed for genotypes G7, G8 and G9. In the second environment (E2), during the first vegetation season at the amelioration of 25 tha−1 phosphor gypsum applied, the greatest overall mean values for grain weight per plant were denoted for wheat varieties G5 (7.09 g), G6 (7.08 g), G9 (6.97 g) and G7 (6.58 g). In the second season, with same level of amelioration (E5), the greatest overall mean values were denoted for wheat varieties G7 (7.74 g), G8 (7.61 g), G6 (7.59 g) and G2 (6.40 g), Table 7. In the third environment, during the first vegetation season at the amelioration of 50 tha−1 phosphor gypsum applied (E3), the greatest overall mean values for grain weight per plant were denoted for wheat varieties G6 (6.78 g), G7 (5.96 g), G8 (5.23 g) and G9 (7.35 g). In the second season, with the same level of amelioration (E6), the greatest overall mean values of grain weight per plant were denoted for wheat varieties G5 (5.56 g), G9 (5.70 g), G6 (6.53 g) and G1 (4.90 g). Presented results revealed that different treatments influenced the differences in grain weight per plant (Table 6).
According to the results, it was observed that the greater response of grain weight per plant were expressed in growing conditions of second vegetation season, due to the meteorological conditions, as the average value was the lowest in the control variant on solonetz without ameliorative measures. In this investigation, the temperature and precipitation regimes were different across years and each cultivar expressed a specific response to environmental conditions. More favorable conditions of the first vegetation season exhibited a more significant impact on minimizing the differences between the average values per treatment. During the second vegetation season, which began with extreme drought conditions after wheat planting, followed with long term frosts during winter without snow cover and unfavorable conditions during critical stages of grain development, the differences were increased and the effects of ameliorative measures were the most manifested (Table 6). According to Khan et al. [57] and Knežević et al. [55], high temperature shortens the grain filling period and induces early maturity, which causes grain shrinkage and low grain weight.
However, aside from climatic conditions, tested wheat genotypes on solonetz soil reacted differently on different levels of soil repair. The phenotypic expressions of grain weight per plant of certain genotypes were constant in various environments, whereas some other wheat genotypes exposed significant variation over different soil amelioration treatments. In this sense, genotypes G1, G4, G5, G9 and G10 expressed a higher sensitivity to the treatment and reacted favorably on both levels of soil ameliorative measures in the first season. While, in the second season, all tested genotypes reacted favorably on both amounts of applied measures. In general, the first vegetation period provided generally higher mean values of the grain weight per plant than the second vegetation period.
However, manifested differences due to years justify the establishment of multi-year trials to evaluate cultivars’ stability. When cultivars are evaluated across several years, the assessment is more reliable and cultivar means are characterized by a smaller error of estimation. These approaches provide support in evaluating cultivar adaptation and selecting cultivars with wide adaptation to environmental conditions [58,59].
The combined ANOVA showed that all three sources of variation (genotypes, treatments and environment) had a significant influence on phenotypic variation of trait grain yield per plant (Table 7). In the combined analysis of variance, the main effects, genotypes and environments were highly significant (52 + 136.1)/284.1 and explain 66.2% of total variation. Participation of the genotype variation in the treatments sum of squares (SS) amounted to 18.3%, while 48.9% of the total sum of squares was attributable to environmental effects (Table 8). A large sum of squares for the environments revealed that the trial environments were diverse, with large differences among environmental means, causing variation on grain yields per plant of wheat. This result is in accordance with results obtained by Rad et al. [51], Mohammadi et al. [42] and Popović et al. [36].
The genotype by environment interaction (GEI) expressed no significant mean square, leading to the conclusion that no cross interaction was expressed for the grain yield per plant of wheat. Furthermore, the additional analysis of the GEI using the PCA analysis confirmed the statistical significance of the first source of variation (i.e., the first main principal component IPCA1), which participated in the GEI variation with 41.15% (Table 7). Therefore, the interaction of the ten wheat genotypes within six environments was best predicted by the first interaction principal component of genotypes and environments. This result is expected, since that grain weight per plant, as well as, per spike, is a quantitative trait controlled by a number of minor genes and under the high influence of environmental conditions. Change in grain weight per plant drastically influences the final grain yield. A high share of the first principal component in the overall GEI, confirmed by AMMI analysis of grain yield, has been reported by Banjac et al. [12] and Popović et al. [36].
According to graphical representation of genotypes and environments in the biplot and in terms of average values (Table 6), a small distance of the environmental points from the origin (zero point) could be noticed, which indicates that the environments E6, E1 and E2 were assessed as the most favorable to achieve a stable response of grain weight per plant of wheat (Figure 2b). According to the arrangement of E1, E2, E3 and E5 points, it can be concluded that the wheat cultivars achieved higher average values of grain yield per plant in these environments compared to E4 and E6 points. This result does not favor the E3 and E5 environments for obtaining higher values of grain weight per plant, as the high values of interaction indicate the poor stability of this trait. Environments E3 and E5, as well as E4, had the highest interaction values, which determines them the least suitable for stable establishment grain weight per plant of wheat. Genotypes G3 and G4 were the most stable over all environments, indicating almost no cross interaction. Small distances of the genotypes points from the origin were also observed for genotypes G10, G6, G9 and G1. Genotypes with above average means, such as G6 and G9, could be selected based on grain weight per plant, while genotypes G8, G7 and G5 had a high distance from the average environment ordinate, thus showing them to be more variable and less stable across environments. Wheat variety G6 appeared to be better adapted to the less favorable conditions of the E1 control environment, showing good stability for this trait with small GEI and responded well, keeping its average above than overall mean. Varieties G2 and G8 responded favorably within treatments of 25 tha−1 phosphogypsum applied in the second growing season, while genotypes G9, G5, G1, G4 and G10 responded favorably within the same treatments applied in the first growing season, (Figure 1b). Genotypes G2 and G3 expressed a positive effect of interaction in E3, within soil treatment of 50 tha−1 phosphogypsum applied (Figure 1b). However, according to Gauch and Moran [32], genotypes at the top of graph (such as G2, G3, G6, G7 and G8) showed positive G × E interactions with environments, while genotypes which were at the bottom (such as a G1, G4, G5, G9 and G10) showed negative G × E interactions with environments (Figure 2b).
These results are in accordance with AMMI selection of wheat genotypes across the testing environments (Table 8). Results of the mean grain weight per plant of the environments, and their respective scores on the first axis of GEI, showed that genotype G6 could be considered as the most promising genotype, since it ranked within the first four AMMI selections in all environments. It showed consistent performance, implying its best adaptability for all soil conditions. In soil environment E1 (control, soil without amelioration) in the first vegetation season, genotype G6 was first ranked, followed by genotypes G7, G8 and G9. In soil environment E4 (control, soil without amelioration) in the second vegetation season, genotype G5 was first ranked, followed by genotypes G9, G6 and G1.
These results indicate that these genotypes have positive response and a stable reaction in less favorable conditions in control environments. In soil environment E2 (treatment with 25 tha−1 phosphogypsum) in the first vegetation season, genotype G5 was first ranked, followed by genotypes G9, G6 and G7. In soil environment E5 (treatment with 25 tha−1 phosphogypsum) in the second vegetation season, genotype G7 was first ranked, followed by genotypes G8, G6 and G2, which indicates that these genotypes have a positive response on soil melioration treatment with 25 tha−1 phosphogypsum. In soil environment E3 (treatment with 50 tha−1 phosphogypsum) in the first vegetation season, genotype G5 was first ranked, followed by genotypes G9, G6 and G1. In soil environment E6 (treatment with 50 tha−1 phosphogypsum) in the second vegetation season, genotype G6 was first ranked, followed by genotypes G7, G8 and G9, which indicates that these genotypes have a positive response on soil melioration treatment with 50 tha−1 phosphogypsum. The results, shown in Table 8, confirm that it is possible to successfully advise the right genotype for all environments, or specific genotypes for specific environments, through AMMI evaluation.

3.3. Grain Yield

The presented results showed that the mean values of grain yield of wheat, expressed in gm−2, ranged between 705.6 (E3) and 728.7 g (E1) within different environments in first vegetation season. In the second vegetation season, the mean values of grain yield of wheat ranged between 466.0 (E4) and 717.2 g (E5) within different environments. Between genotypes, the mean values of grain yield ranged from 564.70 (G3) to 760.70 g (G6) within different genotypes (Table 9). At the control variant (E1), with no amelioration applied, the greatest overall mean values for grain yield in first vegetation season were denoted for wheat varieties G8 (819.90 g), G5 (816.50 g) G4 (785.90 g) and G4 (0.54%). In the second vegetation season, at the control variant (E4), the greatest overall mean values were denoted for genotypes G6 (632.60 g), G7 (631.50 g), G9 (538.90 g) and G8 (535.80 g), Table 9. In the second environment, with soil amelioration of 25 tha−1 phosphor gypsum applied, during the first vegetation season (E2), the greatest overall mean values for grain yield were denoted for wheat varieties G5 (837.20 g), G6 (767.00 g), G4 (740.80 g) and G9 (736.90 g). In the second vegetation season, with the same level of amelioration (E5), the greatest overall mean values were denoted for wheat varieties G7 (889.40 g), G6 (865.30 g) and G8 (848.80 g), Table 9.
In the third environment during the first vegetation season at the amelioration of 50 tha−1 phosphor gypsum applied (E3), the greatest overall mean values for grain yield were denoted for wheat varieties G9 (863.50 g), G5 (822.30 g) and G6 (804.80 g). In the second vegetation season, with same level of amelioration (E6), the greatest overall mean values of grain yield were denoted for wheat varieties G6 (731.90 g), G7 (714.30 g) and G8 (675.40 g), Table 9.
Presented results revealed that different treatments influenced the differences in grain yield of wheat, but also depended on vegetation season. According to the results, a higher sensitivity of grain yield to the meteorological conditions were observed since the overall yield average was the lowest in the control variant in growing conditions of 2005/2006. It can be clearly seen from Table 9 that more favorable conditions in the growing season of 2004/2005 exhibited a more significant impact on reducing the differences between the mean averages per treatment. On the other side, in the second season of 2005/2006, which was characterized by less rainfall in May and stronger precipitation during June, the effects of ameliorative measures were the most evident. Higher sensitivity of grain yield in the second season was caused by high precipitation during flowering and further important developmental stages in wheat. Apart from climatic conditions, genotypes G1, G5, G6, G9 and G10 expressed a higher sensitivity to the treatment and reacted favorably on ameliorative measures in both seasons. Observed variations for grain yield found among these genotypes indicated the high genetic potential for further yield improvement. The similar findings and the influence of climatic conditions on grain yield of wheat has been reported by several authors [22,36].
The combined analysis of variance showed that all three sources of variation, genotypes, treatments and their interactions were highly significant and had significant influence on phenotypic variation of trait grain yield per square meter (Table 10). In the combined analysis of variance, the main effects, genotypes and environments were highly significant (617,250 + 1,552,340)/3,351,528 and explain 64.7% of total variation. The combined ANOVA showed grain yield of wheat was significantly affected by the environment, because significant variance at the 1% level explained 46.3% of the total variation, while the genotype contributed 18.5% of the total variation of the experiment (Table 10).
A large sum of squares for environments indicated that the trial environments were diverse with high differences among environmental means, causing variation on grain yields of wheat. Genotype by environment interaction (GEI) expressed significant mean square, which suggests that grain yield of wheat genotypes varied across environments. The AMMI analysis revealed the complex nature of GEI and contributed 35.2% of the total sum of squares. The additional analysis of the GEI using the PCA analysis confirmed the statistical significance of the first two main components, IPCA 1 and IPCA 2 (Table 10). Separately, IPCA1 and IPCA2 participated in the GE variation with 41.2% and 29.3%, respectively, both with statistically significant effects on the GE interaction variation. These two main components jointly explained more than 70% of the variation of the genotype by environment interaction. This result is expected, since grain yield is a quantitative trait controlled by a number of minor genes and under the high influence of environmental factors.
Significant interactions between environment and wheat cultivars is in accordance with results reported by Banjac et al. [12], Popović et al. [36], Mohammadi et al. [42] and Rad et al. [51]. In addition to the influence of individual plan traits getting more important in terms of grain yield formation per area unit in stressful growing conditions, the combined analysis of variance for grain yield per m−2, compared with grain yield per plant trait, separated the statistical significance of the first two main components, IPCA 1 and IPCA2, which participated in the GE variation.
According to the biplot (Figure 3a) and in terms of the average values (Table 9), it can be observed that environment E6 exposed the lowest distance of the environmental points from the origin (zero point), which indicates that environment E6 can be assessed as the most stable environment to achieve a stable grain yield of wheat (Figure 3a).
According to the arrangement of E1, E2, E3 and E5 environments points, it can be concluded that the wheat cultivars achieved higher average values of grain yield in these environments compared to E4 and E6 points. Considering that high values of interaction indicate the poor stability of a trait, this result does not favor the environments E3 and E5 for obtaining higher values of grain yield. Environments E3 and E5 showed the highest interaction values, which determines them the least suitable for stable establishment grain yield. Regarding genotypes dispersion, genotype G3 was the most stable across all environments, indicating almost no cross interaction. The small distances from the origin were also observed for genotypes G4, G6, G10, G1, G9 and G7. Genotypes G2, G5 and G8 had a high distance from the average environment ordinate, exposed more variables and were less stable across the environments. Genotypes G6 and G7 appeared to be better adapted to the less favorable conditions of the E1 environment, keeping its average above the overall mean. This indicates that these genotypes were better adapted to stressful conditions of abiotic stress solonetz, in comparison to others cultivars. Genotypes G8 and G7 appeared to be better adapted to the conditions of the E5 environment, with its average above the overall mean, while genotypes G9 and G5 reacted favorably and responded well to melioration solonetz, since they achieved higher values of grain yield in environments E2 and E3. Genotypes G1, G2, G3, and G10 appeared to be better adapted in the E6 environment, having mean value of grain yield below the overall mean, while genotype G4 was near the overall mean (Figure 3a). According to Gauch and Moran [32], genotypes at the top of graph (such as G2, G3, G6, G7 and G8) showed positive G × E interactions with environments, while genotypes which were at the bottom (such as a G1, G4, G5, G9 and G10) showed negative G × E interactions with environments (Figure 3a).
These results are in accordance with AMMI selection of wheat genotypes across the testing environments (Table 11). Based on the tested genotypes, genotype G6 could be considered as the most promising genotype, since it ranked within first five AMMI selections in all environments, implying its best adaptability for all soil conditions. Genotypes G8, G5, G4, G2 and G6 were better adapted to stressful conditions of abiotic stress solonetz in the first growing season, while in the second season genotypes G6, G7, G9, G8 and G1 were better adapted.
Genotypes G5, G6, G4, G9 and G7 responded well to melioration solonetz in the amount of 25 tha−1 in the first season, while in the second season, at the same level of repair, genotypes G7, G6, G8, G9 and G2 reacted positively. Genotypes G9, G5, G6, G1 and 10 responded well to melioration solonetz in the amount of 50 tha−1. In the second season, at the same level of repair, genotypes G6, G7, G8, G9 and G5 reacted positively to the reparation of soil.

3.4. Harvest Index

Grain weight per spike plays a significant role in the yield formation, because it directly affects harvest index [46]. Increased yields in winter wheat cultivars have been found to be largely attributable to improved partitioning of biomass to the grain (i.e., higher harvest index). However, there is a biological upper limit to harvest index and; therefore, breeders need to exploit increased biomass production as the mechanism by which yields are increased [60]. Harvest index, as a measure of reproductive efficiency, presents the ratio of grain to total shoot dry matter. It has been determined by interactions between genotypes, environment and crop management [61].
The presented results showed that the mean values of harvest index ranged within different environments between 0.49% (E3) and 0.51% (E1) in the first vegetation season. In the second vegetation season, the mean values of harvest index ranged between 0.40% (E4) and 0.50% (E6) within different environments. Between genotypes, the mean values of harvest index ranged from 0.44% (G2) to 0.51% (G6 and G10) within different genotypes (Table 12). At the control variant (E1), in the first vegetation season, soil with no amelioration applied, the greatest overall mean values for grain yield were denoted for wheat varieties G10 (0.54%), G4 (0.54%), G8 (0.53%) and G6 (0.51%). In the second season at the control variant (E4), the greatest overall mean values were denoted for genotypes G6 (0.51%), G9 (0.48%) and G7 (0.54%), Table 12.
In the second environment, with soil amelioration of 25 tha−1 phosphor gypsum applied, during the first vegetation (E2), the greatest overall mean values for harvest index were denoted for wheat varieties G10 (0.54%), G9 (0.51%), G5 (0.51%) and G1(0.51%). In the second season, with same level of amelioration (E5), the greatest overall mean values were denoted for wheat varieties G10 (0.52%), G8 (0.50%), G4 (0.51%) and G6 (0.48%), Table 12.
In the third environment during the first vegetation season at the amelioration of 50 tha−1 phosphor gypsum applied (E3), the greatest overall mean values for harvest index were denoted for wheat varieties G6 (0.52%), G9 (0.51%), G10 (0.51%) and G1 (0.50%). In the second season, with same level of amelioration (E6), the greatest overall mean values of harvest index were denoted for wheat varieties G10 (0.53%), G4 (0.52%), G8 (0.51%), G6 (0.49%) and G3 (0.50). Presented results revealed that different treatments influenced the differences in grain weight per plant (Table 12).
According to the mean values of both seasons, higher sensitivity of harvest index could be observed in the second growing season, as the environmental average was the lowest in the control variant of soil, without ameliorative measures, as a consequence of variable conditions of the 2005/2006 vegetation period. On the other side, more favorable conditions during the first vegetation season showed significant impact on reducing differences between the mean values per treatment and indicated that the effect of a particular soil treatment was missing. Therefore, the changing weather conditions highly influenced the effect of soil melioration on the harvest index and other grain yield traits of wheat. Since each vegetation season could bring certain differences due to even more unpredictable weather conditions, cultivars’ assessment across several years is more reliable and could bring appropriate recommendation of cultivars characterized by a small error.
The combined ANOVA showed that all three sources of variation (genotypes, treatments and environment) were statistically significant, having significant influence on phenotypic variation of the harvest index. In the combined analysis of variance, the main effects, genotypes and environments were highly significant and explained 60.6% of the total variation. Participation of the genotype variation in the treatment sum of squares (SS) amounted to 13.4%, while 47.1% of the total sum of squares was attributable to environmental effects (Table 13). A large sum of squares for environments indicated the differences between growing seasons, while diversity of treatments caused a considerable sum of squares for environmental factors in the total variation, which indicates that these factors were the most responsible for the variation of the harvest index. Genotype by environment interaction expressed a significant mean square which represented 39.4%. Sum of squares of GEI was 2.9 times higher than the sum of the squares of genotype, indicating that there were significant differences in behavior between genotypes in different environments. Additional analysis of the GEI using the PCA analysis showed a statistical significance of the first source of variation, first main component IPCA1, which participated in the GEI variation with 78.54%. In addition to the two main components jointly explaining more than 90% of the GEI variation, the second source of variation, the second main component IPCA2, expressed no statistically significant effect in the GE interaction variation (Table 13). Significant interactions between environment and wheat cultivars in the harvest index, as a high share of the first main component IPCA1 in the GE variation, have been reported by Dimitrijević et al. [10]. However, a higher influence of genetic factor on the harvest index of wheat was observed by Tayar [62].
According to the biplot (Figure 3b) and in terms of average values (Table 13), it could be noticed that, apart from certain exceptions, agro-ecological environments were quite near the level of the experimental overall average. A small distance of the environmental points from the origin (zero point) indicates that the environments E3 and E2 could be assessed as the most favorable to achieve a stable response of harvest index of wheat. Environments E5, E6 and E1 could be assessed as less stable, while environment E4 had the highest negative interaction value, which determines them the least suitable for the stable establishment harvest index of wheat. Wheat cultivars G1, G5, G6, G7 and G9 achieved higher average values of harvest index in E3 environment, while genotype G10 achieved the greatest value of harvest index in E2 environment. In general, wheat variety G10 (Simonida) expressed the most stable reaction for this trait, showing a PCA1 score close to 0, and it responded well when grown in E2 environment, within soil amelioration of 25 t ha−1 phosphogypsum applied, (Figure 3b). Wheat genotypes G3, G4, G8 and at least G2 responded favorably within treatments of 50 tha−1 phosphogypsum applied. According to Gauch and Moran [32], G × E interactions are estimated based on the position of wheat genotypes within the graph. Therefore, genotypes at the top of graph (such as G2, G3, G4, G8 and G10) showed positive G × E interactions with environments, while genotypes which were at the bottom (such as a G1, G5, G6, G7 and G9) showed negative G × E interactions with environments (Figure 3b). However, in addition to showing the most stable reaction, genotype G10 achieved the greatest mean value of harvest index. In this analysis, wheat genotype G6 also expresses the greatest value of harvest index, but followed with high interaction values, which indicates it as the least suitable for the stable establishment harvest index. Wheat varieties with improved harvest index are desirable since they expressed increased physiological efficiency of nutrient reutilization and efficiency in translocation of nutrients from leaves and stems into grains.
These results are in accordance with AMMI selection of wheat genotypes across the testing environments (Table 14).
Results of the harvest index of the environments and their respective scores on the first axis of GEI showed similar results as previous analysis. In the first vegetation season, in environment E1 (control, soil without amelioration), genotype G10 could be considered as the most promising genotype, since it ranked first, followed with genotypes G4, G8 and G6. In soil environment E4 (control without amelioration) in the second vegetation season, genotype G6 was first ranked, followed by genotypes G9, G7 and G1. These results indicate that these genotypes showed positive response and stable reaction in less favorable conditions on control environments. In soil environment E2 (treatment with 25 tha−1 phosphogypsum), in the first vegetation season, genotype G9 was first ranked, implying its best adaptability for this soil conditions, followed by genotypes G5, G1 and G10. In soil environment E5 (treatment with 25 tha−1 phosphogypsum), in the second vegetation season, genotype G10 was first ranked, followed by genotypes G8, G4 and G6, which indicates that these genotypes showed positive response on soil melioration treatment with 25 tha−1 phosphogypsum. In soil environment E3 (treatment with 50 tha−1 phosphogypsum), in the first vegetation season, genotype G9 could be the most promising genotype, since it ranked first, followed by genotypes G5, G1 and G10. In soil environment E6 (treatment with 50 tha−1 phosphogypsum), in the second season, genotype G10 was first ranked, followed by genotypes G4, G8 and G6, which indicates that these genotypes showed positive response on soil amelioration with amount of 50 tha−1 phosphogypsum. These results support the recommendation of appropriate genotype for all environments, or specific genotypes for specific environments, through AMMI evaluation, Table 14.
According to the results, a greater sensitivity of harvest index was observed in less favorable conditions in the second vegetation period. In the second vegetation period, the general mean value of the harvest index was the lowest in the control variant, without ameliorative measures. More favorable conditions in the first vegetation season had a significant impact on reducing the differences between the averages values per treatment. This result is expected, since the harvest index is a highly variable trait and depends on the weather conditions during the season. Therefore, in a more favorable year the harvest index is higher. Since the harvest index represents efficiency of translocation of the nutrients of plants from vegetative to generative parts, and this trait highly depends on two pronounced quantitative traits, grain weight per plant and plant biomass, it is a highly variable trait. The similar results and high impact of year on this trait were observed by Vuković [63] and Kondic [64]. The results suggested that the harvest index is a critical factor for grain yield of wheat under variable weather conditions, terminal high temperatures, as well as in water shortages.
In general, obtained results indicated that even in wheat, which is considered to be durable and well adapted to a wide range of environments, in addition to the influence of different treatments, different climatic conditions can also cause significant differences in some important yield related traits. Furthermore, variable genotype reactions to different environments point out the importance of the evaluation of different wheat cultivars in order to identify and recommend the optimal genotype for specific growing conditions [64,65,66].

4. Conclusions

Based on the presented findings it can be concluded that the application of complex ameliorative measures could increase the production potential of the solonetz soil near to a level of productive soils. The results also indicated that amelioration soil measures, with adequate selection of appropriate wheat genotypes and adequate management practices, could mean that these soils could provide high and stable wheat yields. Differences in growing seasons contributed considerably to treatment effects. The more favorable weather conditions of the first season exhibited a more significant impact on minimizing the differences between mean values per treatment and; therefore, the greatest overall stability was evident in the first vegetation season.
Based on the tested genotypes, genotype G6 (Evropa 90) expressed a stable reaction over all environments for the traits plant height, grain weight per plant and grain yield per square meter. In addition to this stable reaction, this genotype resulted in greater values of observed traits and could be considered as one of the most promising genotypes in all soil environments. The greater expression for both plant height and grain yield of genotype in all soil conditions could be attributed to the high capacity of the stem to accumulate sufficient stem reserves for the partitioning to grain. This indicates that the genotypes that could maintain their height in stress environments could be desirable. However, genotype G10 (Simonida) could be considered as a wheat variety with improved harvest index, since it expressed the most stable reaction for this trait, showing a PCA1 score close to 0, followed with the greatest mean value of the harvest index. On the other side, genotype G6 (Evropa 90), which showed also the greatest mean values for the harvest index, would be recommended if the GEI was disregarded. However, high yielding, but not stable wheat genotypes across environments could be recommended the specific environments where they performed well. Information in this regard would help wheat producers, as well as breeders, to select the appropriate wheat cultivars that can be successfully exploited in conditions of less productive soils, such as solonetz type.

Author Contributions

Conceptualization, V.P. and N.L.; methodology, M.K.; software, D.P.; formal analysis, M.K., S.E.M. and V.Ć.; investigation, N.L.; resources, V.P.; data curation, D.P., M.P. and B.I.; writing—original draft preparation, V.P., D.P. and N.L.; writing—review and editing, N.L., and V.Ć.; visualization, D.P.; supervision, M.K. and S.J.; project administration, V.P. and B.I.; funding acquisition, N.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors express their appreciation for providing financial support to Ministry of Education, Science and Technological Development of the Republic of Serbia (Grant No. 451-03-9/2021-14/200358, 200032 and 200045); FAO Project 2020–2022—Redesigning the exploitation of small grains genetic resources towards increased sustainability of grain-value chain and improved farmers’ livelihoods in Serbia and Bulgaria—GRAINEFIT; and Bilateral Project Montenegro-Serbia, 2019–2021, Alternative cereals and oil crops as a source of healthcare food and an important raw material for the production of biofuel.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tabaković, M.; Sečanski, M.; Stanisavljević, R.; Mladenović-Drinić, S.; Simić, M.; Knežević, J.; Oro, V. The impact of agroecological factors on morphological traits of maize. Genetika 2020, 52, 1203–1213. [Google Scholar] [CrossRef]
  2. Popović, V.; Jovović, Z.; Marjanović-Jeromela, A.; Sikora, V.; Mikić, S.; Šarčević-Todosijević, L. Climatic Change and Agricultural Production. In Proceedings of the GEA Conference 2021, GEA (Geo-Eco-Eco-Agro), Podgorica, Montenegro, 27–31 March 2020; pp. 160–166. [Google Scholar]
  3. Petrovic, S.; Dimitrijevic, M.; Belic, M.; Banjac, B.; Bosković, J.; Zecević, V.; Pejić, B. The variation of yield components in wheat (Triticum aestivum L.) in response to stressful growing conditions of alkaline soil. Genetika 2010, 42, 545–555. [Google Scholar] [CrossRef]
  4. Belić, M.; Nešić, L.; Ćirić, V. Agriculture Contribution Towards Clean Environment and Healthy Food. Handbook; Faculty of Agriculture, University of Novi Sad: Novi Sad, Serbia, 2014. [Google Scholar]
  5. Pessoa, L.G.M.; Freire, M.N.G.; Filho, J.; Santos, P.R.; Miranda, M.F.A.; Freire, F.J. Characteri-zation and classification of halomorphic soils in the semiarid region of Northeastern Brazil. JAS 2019, 11, 405–418. [Google Scholar] [CrossRef]
  6. Choudhury, S.G.; Srivastava, S.; Singh, R.; Chaudhari, S.K.; Sharma, D.K.; Singh, S.K.; Sarkar, D. Tillage and residue management effects on soil aggregation, organic carbon dynamics and yield attribute in rice–wheat cropping system under reclaimed sodic soil. Soil Till. Res. 2014, 136, 76–83. [Google Scholar] [CrossRef]
  7. Bogunovic, I.; Vukadinovic, V.; Kisic, I.; Chiavalon, S.; Vucic, H.; Durdevic, B. Tillage and Soil Amendments Effect on Soil Physical Properties and Yield of Oats (Avena sativa L.) in Organic Farm in Mediterranean Croatia. Agric. Conspec. Sci. 2018, 83, 17–23. [Google Scholar]
  8. Ćirić, V.; Manojlović, M.; Belić, M.; Nešić, L.; Šeremešić, S. Aggregate Stability and Sealing Risk of Solonetz under Different Land Use Regimes. Ratar. Povrt. 2012, 49, 243–249. [Google Scholar]
  9. Bowser, W.E. Genesis and Characteristics of Solonezt Soils. In Soils in Canada; Legget, R., Ed.; University of Toronto Press: Toronto, ON, Canada, 2019; pp. 165–173. [Google Scholar] [CrossRef]
  10. Dimitrijević, M.; Petrović, S.; Belić, M.; Mladenov, N.; Belić, M.; Hristov, N.; Banjac, B.; Vukosavljev, M. Phenotypic reaction of wheat grown on different soil types. Genetika 2009, 41, 169–177. [Google Scholar] [CrossRef]
  11. Miller, J.J.; Brierley, J.A. Solonetzic soils of Canada: Genesis, distribution, and classification. Can. J. Soil Sci. 2011, 91, 889–916. [Google Scholar] [CrossRef]
  12. Banjac, B.; Mladenov, V.; Dimitrijevic, M.; Petrovic, S.; Boćanski, J. Genotype x environment interactions and phenotypic stability for wheat grown in stressful conditions. Genetika 2014, 46, 799–806. [Google Scholar] [CrossRef]
  13. Nayak, A.K.; Mishra, V.K.; Sharma, D.K.; Jha, S.K.; Singh, C.S.; Shahabuddin, M.; Shahid, M. Efficiency of Phosphogypsum and Mined Gypsum in Reclamation and Productivity of Rice–Wheat Cropping System in Sodic Soil. Commun. Soil Sci. Plant. Anal. 2013, 44, 909–921. [Google Scholar] [CrossRef]
  14. Nayak, S.; Mishra, C.S.; Guru, B.C.; Rath, M. Effect of phos-phogypsum amendment on soil physico-chemical proper-ties, microbial load and enzyme activities. J. Environ. Biol. 2011, 32, 613–617. [Google Scholar]
  15. Mesić, M.; Brezinščak, L.; Zgorelec, Ž.; Perčin, A.; Šestak, I.; Bilandžija, D.; Trdenić, M.; Lisac, H. The application of phosphogypsum in agriculture. Agric. Conspec. Sci. 2016, 81, 7–13. [Google Scholar]
  16. Al-Oudat, M.; Al-Rayyes, A.H.; Kanacri, S. Effects of Phosphogypsum on Chemical Properties of Soil, Wheat Growth and Accumulation of Radionuclides, Trace Elements and Fluoride in Soil and Plants (AECS-A/RSS--234); Department of Agriculture: Damascus, Syrian, 1998.
  17. Milošev, D.; Belić, M.; Nešić, L.; Dimitrijević, M.; Petrović, S.; Šeremešić, S. The Application of Phosphogypsum in the Winter Wheat Production on Solonetz Soil; Zbornik Radova Instituta za Ratarstvo i Povrtarstvo: Novi Sad, Serbia, 2007; Volume 43, pp. 47–53. [Google Scholar]
  18. Temel, S.; Keskin, B.; Simsek, U.; Yilmaz, I. Performance of some forage grass species in halomorphic soil. Turk. J. Field Crop. 2015, 20, 131–141. [Google Scholar] [CrossRef] [Green Version]
  19. Cairns, R.R. Solonetzic Soils and Their Management; Minister of Supply and Servies, Canada Department of Agriculture: Ottawa, ON, Canada, 1978.
  20. Peters, T.W. Solonetzic Soils Technology. In Solonetzic Soils Technology and Management; Toogood, J.A., Cairns, R.R., Eds.; University of Alberta: Edmonton, AB, Canada, 1978; p. 821. [Google Scholar]
  21. Ljubičić, N.; Petrovic, S.; Dimitrijevic, M.; Sekulić, P.; Mladenov, N.; Belić, M.; Vukosavljev, M. Kernel number per spike stability of wheat genotypes grown on meliorated soil (solonetz). Contem. Agric. 2010, 59, 332–338. [Google Scholar]
  22. Mladenov, V.; Dimitrijević, M.; Petrović, S.; Boćanski, J.; Banjac, B.; Kondić-Špika, A.; Trkulja, D. Genetic analysis of spike length in wheat. Genetika 2019, 51, 167–178. [Google Scholar] [CrossRef] [Green Version]
  23. Lantican, M.A.; Braun, H.J.; Payne, T.S.; Singh, R.P.; Sonder, K.; Baum, M.; van Ginkel, M.; Erenstein, O. Impacts of International Wheat Improvement Research, 1994–2014; CIMMYT: Mexico City, Mexico, 2016. [Google Scholar]
  24. Erkul, A.; Unay, A.; Konak, C. Inheritance of yield and yield components in bread wheat (Triticum aestivum L.) cross. Turk. J. Field Crop. 2010, 15, 137–140. [Google Scholar]
  25. Ljubičić, N.; Petrović, S.; Dimitrijević, M.; Hristov, N.; Kostić, M.; Kondić-Špika, A.; Jevtić, R. Diallel analysis of some important grain yield traits in bread wheat. Turk. J. Field Crop. 2017, 22, 1–7. [Google Scholar] [CrossRef] [Green Version]
  26. Gauch, H.G.; Piepho, H.P.; Annicchiarico, P. Statistical analysis of yield trials by AMMI and GGE: Further considerations. Crop. Sci. 2008, 48, 866. [Google Scholar] [CrossRef]
  27. Yan, W. Analysis and handling of G × E in a practical breeding program. Crop. Sci. 2016, 56, 2106–2118. [Google Scholar] [CrossRef]
  28. Studnicki, M.; Kang, M.S.; Iwańska, M.; Oleksiak, T.; Wójcik-Gront, E.; Mądry, W. Consistency of Yield Ranking and Adaptability Patterns of Winter Wheat Cultivars between Multi-Environmental Trials and Farmer Surveys. Agronomy 2019, 9, 245. [Google Scholar] [CrossRef] [Green Version]
  29. Júnior, B.L.A.Y.; Da Silva, C.P.; De Oliveira, L.A.; Nuvunga, J.J.; Pires, L.P.M.; Von Pinho, R.G.; Balestre, M. AMMI Bayesian Models to Study Stability and Adaptability in Maize. Agron. J. 2018, 110, 1765–1776. [Google Scholar] [CrossRef]
  30. Bocianowski, J.; Nowosad, K.; Tomkowiak, A. Genotype—environment interaction for seed yield of maize hybrids and lines using the AMMI model. Maydica Electron. Publ. 2019, 64, 1–8. [Google Scholar]
  31. Gauch, H.G.; Rodrigues, P.C.; Munkvold, J.D.; Heffner, E.L.; Sorrells, M. Two New Strategies for Detecting and Understanding QTL × Environment Interactions. Crop. Sci. 2011, 51, 96–113. [Google Scholar] [CrossRef]
  32. Gauch, H.G. A simple protocol for AMMI analysis of yield trials. Crop. Sci. 2013, 53, 1860. [Google Scholar] [CrossRef]
  33. Hongyu, K.; Garcıa-Pena, M.; de Araujo, L.B.; dos Santos Dias, C.T. Statistical analysis of yield trials by AMMI analysis of genotype × environment interaction. Biom. Lett. 2014, 51, 89–102. [Google Scholar] [CrossRef] [Green Version]
  34. Hongyu, K. Adaptability, stability and genotype by environment interaction using the AMMI model for multienvironmental trials. Biodiversity 2018, 17, 10–21. [Google Scholar]
  35. 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] [PubMed] [Green Version]
  36. Popović, V.; Ljubičić, N.; Kostić, M.; Radulović, M.; Blagojević, D.; Ugrenović, V.; Popović, D.; Ivošević, B. Genotype × Environment Interaction for Wheat Yield Traits Suitable for Selection in Different Seed Priming Conditions. Plants 2020, 9, 1804. [Google Scholar] [CrossRef] [PubMed]
  37. Gauch, H.G.; Moran, D.R. AMMISOFT for AMMI Analysis with Best Practices. Soil Crop. Sci. 2019. [Google Scholar] [CrossRef]
  38. Belić, M.; Nešić, L.; Dimitrijevic, M.; Petrovic, S.; Ćirić, V.; Pekeč, S.; Vasin, J. Impact of reclamation practices on the content and qualitative composition of exchangeable base cations of the solonetz soil. Aust. J. Crop. Sci. 2012, 6, 1471–1480. [Google Scholar]
  39. Živković, B.; Nejgebauer, V.; Tanasijević, Đ.; Miljković, N.; Stojković, L.; Drezgić, P. Zemljišta Vojvodine [Soil of Vojvodina]; Institut za Poljoprivredna Istraživanja [Institute of Field and Vegetable Crops]: Novi Sad, Serbia, 1972; pp. 380–403. [Google Scholar]
  40. Zobel, R.W.; Wright, M.J.; Gauch, H.G. Statistical analysis of yield trial. Agron. J. 1998, 80, 388–393. [Google Scholar] [CrossRef]
  41. Oliveira, J.P.; Moreira, W.N.; Duarte, J.B.; Chaves, L.J.; Pinheiro, J.B. Genotype-environment interaction in maize. Crop. Breed. Appl. Biotechnol. 2003, 3, 185–192. [Google Scholar] [CrossRef] [Green Version]
  42. Mohammadi, R.; Aemion, M.; Zadhasan, E.; Ahmadi, M.M.; Amri, A. The use of AMMI model for interpreting genotype x environment interaction in durum wheat. Exp. Agric. 2018, 54, 670–683. [Google Scholar] [CrossRef]
  43. Crossa, J. Statistical analyses of multi-location trials. Adv. Agron. 1990, 44, 55–85. [Google Scholar]
  44. Shaibu, A.S.; Adnan, A.A.; Jibrin, J.M. Stability performance of extra early maturing maize (Zea mays L.) varieties under high and low nitrogen environments in Sudan Savanna. Cogent. Food Agric. 2016, 2, 1231456. [Google Scholar] [CrossRef]
  45. VSN International. Available online: www.vsn-intl.com (accessed on 11 December 2020).
  46. Popović, V.; Glamočlija, Đ.; Malešević, M.; Ikanović, J.; Dražić, G.; Spasić, M.; Stanković, S. Genotype specificity in nitrogen nutrition of malting barley. Genetika 2011, 43, 197–204. [Google Scholar] [CrossRef]
  47. Vasileva, V.; Vasilev, V. Agronomic characterization and the possibility for potential use of subterranean clover (Trifolium subterraneum L.) in the forage production in Bulgaria. Pak. J. Bot. 2020, 52, 1–5. [Google Scholar] [CrossRef]
  48. Mitrović, B.; Stojaković, M.; Zorić, M.; Stanisavljević, D.; Bekavac, G.; Nastasić, A.; Mladenov, V. Genetic gains in grain yield, morphological traits and yield stability of middle-late maize hybrids released in Serbia between 1978 and 2011. Euphytica 2016, 211, 321–330. [Google Scholar] [CrossRef]
  49. Kang, M. Genotype-Environment Interaction and Stability Analyses: An Update. In Quantitative Genetics, Genomics and Plant Breeding; Kansas State University: Manhattan, KS, USA, 2020. [Google Scholar]
  50. Božović, D.; Popović, V.; Rajičić, V.; Kostić, M.; Filipović, V.; Kolarić, L.; Ugrenović, V.; Spalević, V. Stability of the expression of the maize productivity parameters by AMMI models and GGE-biplot analysis. Not. Bot. Horti. Agrobot. Cluj-Napoca. 2020, 48, 1387–1397. [Google Scholar] [CrossRef]
  51. Rad, M.R.N.; Kadir, M.A.; Rafii, M.Z.; Jaafar, H.Z.E.; Naghavi, M.R.; Ahmadi, F. Genotype × environment interaction by AMMI and GGE biplot analysis in three consecutive generations of wheat (Triticum aestivum) under normal and drought stress conditions. Aust. J. Crop. Sci. 2013, 7, 956–961. [Google Scholar]
  52. Popović, V.; Vučković, S.; Jovović, Z.; Ljubičić, N.; Kostić, M.; Rakaščan, N.; Glamočlija-Mladenović, M.; Ikanović, J. Genotype by year interaction effects on soybean morpho-productive traits and biogas production. Genetika 2020, 52, 1055–1073. [Google Scholar] [CrossRef]
  53. Branković-Radojčić, D.; Babić, V.; Girek, Z.; Živanović, T.; Radojĉić, A.; Filipović, M.; Srdić, J. Evaluation of maize grain yield and yield stability by AMMI analysis. Genetika 2018, 50, 1067–1080. [Google Scholar] [CrossRef]
  54. Ugrenović, V.; Popović, V.; Ugrinović, M.; Filipović, V.; Mačkić, K.; Ljubičić, N.; Popović, S.; Lakić, Ž. Black Oat (Avena strigosa Schreb.) Ontogenesis and Agronomic Performance in Organic Cropping System and Pannonian Environments. Agriculture 2021, 11, 55. [Google Scholar] [CrossRef]
  55. Knežević, D.; Radosavac, A.; Zelenika, M. Variability of grain weight per spike in wheat grown in different ecological conditions. Acta. Agric. Serbica 2015, 20, 85–951. [Google Scholar] [CrossRef]
  56. Protich, R.; Todorovich, G.; Protich, N. Grain weight per spike of wheat using different ways of seed protection. Bulg. J. Agric. Sci. 2012, 18, 185–190. [Google Scholar]
  57. Khan, F.U.; Mohammad, F.; Raziuddin, A.; Shah, Z.; Ahmad, M.; Shah, Z. Genotypic differences and genotype x nitrogen interactions for yield traits in bread wheat. J. Anim. Plant. Sci. 2017, 27, 1264–1268. [Google Scholar]
  58. Vlachostergios, D.N.; Roupakias, D.G. Response to conventional and organic environment of thirty-six lentil (Lens culinaris Medik.) varieties. Euphytica 2008, 163, 449–457. [Google Scholar] [CrossRef]
  59. Windhausen, V.S.; Wagener, S.; Magorokosho, C.; Makumbi, D.; Vivek, B.; Piepho, H.P.; Melchinger, A.E.; Atlin, G.N. Strategies to subdivide a target population of environments: Results from the CIMMYT-led maize hybrid testing programs in Africa. Crop. Sci. 2012, 52, 2143–2152. [Google Scholar] [CrossRef] [Green Version]
  60. White, E.M.; Wilson, F.E.A. Responses of grain yield, biomass and harvest index and their rates of genetic progress to nitrogen availability in ten winter wheat varieties. Ir. J. Agric. Food Res. 2006, 45, 85–101. [Google Scholar]
  61. Porker, K.; Straight, M.; Hunt, J.R. Evaluation of G × E × M Interactions to Increase Harvest Index and Yield of Early Sown Wheat. Front. Plant. Sci. 2020, 11, 994. [Google Scholar] [CrossRef]
  62. Tayyar, S. Grain yield and agronomic characteristics of Romanian bread wheat varieties under the conditions of Northwestern Turkey. Afr. J. Biotechnol. 2008, 7, 1479–1486. [Google Scholar]
  63. Vuković, N. Multivariate Analysis of Wheat Genotypes Grown on Meliorated Soil. Master’s Thesis, Faculty of Agriculture, University of Novi Sad, Novi Sad, Serbia, 7 October 2009. [Google Scholar]
  64. Lakić, Ž.; Glamočlija, Đ.; Kondić, D.; Popović, V.; Pavlović, S. Fodder Plants and Cereals in the Function of Soil Protection from Degradation; Monograph: Banja Luka, Bosnia and Herzegovina, 2018; pp. 1–403. ISBN 978-99938-93-47-9. [Google Scholar]
  65. Kondić-Špika, A.; Mladenov, N.; Grahovac, N.; Zorić, M.; Mikic, S.; Trkulja, D.; Marjanović-Jeromela, A.; Miladinović, D.; Hristov, N. Biometric Analyses of Yield, Oil and Protein Contents of Wheat (Triticum aestivum L.) Genotypes in Different Environments. Agronomy 2019, 9, 270. [Google Scholar] [CrossRef] [Green Version]
  66. Reynolds, M.P.; Trethowan, R.; Crossa, J.; Vargas, M.; Sayre, K.D. Erratum to “Physiological factors associated with genotype by environment interaction in wheat”. Field Crop. Res. 2004, 85, 253–274. [Google Scholar] [CrossRef]
Figure 1. Weather characteristics of the vegetative seasons: Average monthly air temperature (°C) and precipitation (mm) and perennial average for locations in the first year (a); average monthly air temperature (°C) and precipitation (mm) and perennial average for locations in the second year (b).
Figure 1. Weather characteristics of the vegetative seasons: Average monthly air temperature (°C) and precipitation (mm) and perennial average for locations in the first year (a); average monthly air temperature (°C) and precipitation (mm) and perennial average for locations in the second year (b).
Plants 10 00604 g001
Figure 2. (a) AMMI 1 biplot of ten wheat cultivars across six environments (two years × three treatments) for the estimation of main and multivariate (GEI) effects for the plant height of wheat (cm); (b) AMMI 1 biplot of ten wheat cultivars across six environments (two years × three treatments) for the estimation of main and multivariate (GEI) effects for the grain weight per plant of wheat (g).
Figure 2. (a) AMMI 1 biplot of ten wheat cultivars across six environments (two years × three treatments) for the estimation of main and multivariate (GEI) effects for the plant height of wheat (cm); (b) AMMI 1 biplot of ten wheat cultivars across six environments (two years × three treatments) for the estimation of main and multivariate (GEI) effects for the grain weight per plant of wheat (g).
Plants 10 00604 g002
Figure 3. (a) AMMI 1 biplot of ten wheat cultivars across six environments (two years × three treatments) for the estimation of main and multivariate (GEI) effects for the grain yield in gm−2 at 13% moisture; (b) AMMI 1 biplot of ten wheat cultivars across six environments (two years × three treatments) for the estimation of main and multivariate (GEI) effects for the harvest index.
Figure 3. (a) AMMI 1 biplot of ten wheat cultivars across six environments (two years × three treatments) for the estimation of main and multivariate (GEI) effects for the grain yield in gm−2 at 13% moisture; (b) AMMI 1 biplot of ten wheat cultivars across six environments (two years × three treatments) for the estimation of main and multivariate (GEI) effects for the harvest index.
Plants 10 00604 g003
Table 1. Description of the six different environments used to evaluate ten winter wheat genotypes.
Table 1. Description of the six different environments used to evaluate ten winter wheat genotypes.
Environments CodeGrowing SeasonsTreatments by Phosphogypsum
E12004/2005Control–non-ameliorated solonetz
E22004/2005Solonetz ameliorated by 25 tha−1 phosphogypsum
E32004/2005Solonetz ameliorated by 50 tha−1 phosphogypsum
E42005/2006Control–non-ameliorated solonetz–control
E52005/2006Solonetz ameliorated by 25 tha−1 phosphogypsum
E62005/2006Solonetz ameliorated by 50 tha−1 phosphogypsum
Table 2. Chemical composition of solonetz soil.
Table 2. Chemical composition of solonetz soil.
ParameterpHCaCO3HumusTotal NP2O5K2OSalt
Depth (cm)KClH2O(%)(%)(%)mg/100 g Soil(%)
0–104.505.800.005.950.407.2085.000.03
11–306.307.850.261.680.111.6017.500.13
31–606.908.200.021.470.075.3023.400.18
Table 3. Mean values for the trait plant height of wheat (cm) and interaction IPCA1 values of the AMMI model of ten winter wheat varieties grown in six environments.
Table 3. Mean values for the trait plant height of wheat (cm) and interaction IPCA1 values of the AMMI model of ten winter wheat varieties grown in six environments.
Plant Height (cm)
Genotypes (G)
EG1G2G3G4G5G6G7G8G9G10EmIPCAe1V
E161.059.960.366.965.668.665.564.663.959.363.51.18473.4
E261.360.059.866.464.168.564.965.263.760.963.50.52842.0
E360.358.958.064.561.267.163.064.762.261.562.1−0.20986.8
E439.738.740.247.047.147.845.542.643.335.742.82.28767.2
E557.556.053.359.654.063.558.163.058.262.458.6−1.95239.6
E647.245.743.249.544.053.348.052.648.051.948.3−1.83857.6
IPCAg1 0.4688−0.34680.61650.74132.11040.03330.7631−1.08340.2071−2.57256.51.184124.7
Gm54.553.252.459.056.061.557.558.856.555.3
AMMI: The additive main effects and multiplicative interaction; Gm: Genotype mean; Em: Environmental means; IPCAe1: The first interaction principal component axes for environment; IPCAg1: Interaction principal component axes for genotype; G: wheat genotypes; E: Environmental labels with control (K) and 25 or 50 tha−1 phosphogypsum applied in both seasons; V: Variance.
Table 4. AMMI analysis of variance for the plant height of ten winter wheat varieties examined across six environments.
Table 4. AMMI analysis of variance for the plant height of ten winter wheat varieties examined across six environments.
Source 1dfSSMSF-ValueF-probThe Share of Total Variation %
Total17922,318124.7--
Treatments5914,588247.34.81 **0.000065.36
Genotypes912601402.72 **0.00678.64
Environments511,6912338.112.89 **0.000080.14
Block122176181.33.53 **0.000214.92
Interactions45163836.40.71 ns0.903911.23
IPCA11360046.20.90 ns0.557936.63
IPCA21147643.30.84 ns0.599229.06
Residuals2156226.70.52 ns0.95673.85
Error108555351.4---
1 All sources were tested in relation to the error; **: Highly significant at p < 0.01 probability level; ns: Not significant; df: Degree of freedom; F: F value calculated; IPCA1: The first interaction principal components axes; IPCA2: The second interaction principal components axes.
Table 5. Table of first AMMI selections per environment.
Table 5. Table of first AMMI selections per environment.
No.EMeanIPCAGenotypic Rank
(g)Score12345678910
4E442.82.287G6G4G5G7G8G9G1G3G2G10
1E163.51.184G6G4G8G7G5G9G1G10G2G3
2E263.50.528G6G8G4G7G9G10G5G1G2G3
3E362.1−0.209G6G5G4G7G9G8G3G1G2G10
6E648.3−1.838G6G8G10G4G9G7G1G2G5G3
5E558.6−1.952G6G8G10G4G9G7G1G2G5G3
E: Environments; IPCA score: Score based on the first interaction principal component axes.
Table 6. Mean values for the grain weight per plant of wheat (g) and interaction PCA1 values of the AMMI model of ten winter wheat varieties grown in six environments.
Table 6. Mean values for the grain weight per plant of wheat (g) and interaction PCA1 values of the AMMI model of ten winter wheat varieties grown in six environments.
Grain Weight Per Plant (g)
Genotypes (G)
EG1G2G3G4G5G6G7G8G9G10EmIPCAe1V
E15.786.135.696.286.417.537.447.156.596.106.510.34632.0
E26.245.355.556.327.097.086.586.056.976.286.35−0.61711.4
E36.674.785.506.417.706.785.965.237.356.496.29−1.39672.7
E43.203.783.223.783.795.115.104.864.023.574.040.52363.3
E55.316.405.576.065.807.597.747.616.155.796.400.92481.4
E64.905.094.735.355.566.536.386.385.705.185.550.21921.4
IPCAg1−0.6360.6512−0.017−0.201−0.8700.2970.7180.975−0.566−0.351 2.7
Gm5.355.255.045.706.066.776.536.166.135.57
Gm: Genotype mean; Em: Environmental means; IPCAe1: The first interaction principal component axes for environment; IPCAg1: Interaction principal component axes for genotype; E: Environmental labels with control (K) and 25 or 50 tha−1 phosphogypsum applied in both seasons; V: variance.
Table 7. AMMI analysis of variance for the grain weight per plant (g) of ten winter wheat varieties examined across six environments.
Table 7. AMMI analysis of variance for the grain weight per plant (g) of ten winter wheat varieties examined across six environments.
Source 1dfSSMSF-ValueF-ProbThe Share of Total Variation %
Total179488.42.728---
Treatments59284.14.8163.04 **0.000058.17
Genotypes9525.7733.64 **0.000518.30
Environments5136.127.239.85 **0.000048.97
Block1233.22.7651.75 *0.066911.69
Interactions45962.1341.340.107433.79
IPCA 11339.53.0391.92 *0.035541.15
IPCA 21127.12.461.550.123428.23
Residuals2129.51.4030.890.609310.38
Error108171.11.584---
1 All sources were tested in relation to the error; *: Significant at p < 0.05 probability level; **: Highly significant at p < 0.01 probability level; df: Degree of freedom; F: F value calculated; IPCA: Interaction principal components axes.
Table 8. Table of first AMMI selections per environment.
Table 8. Table of first AMMI selections per environment.
No.EMeanIPCAGenotypic Rank
(g)Score12345678910
5E56.400.9248G7G8G6G2G9G4G5G10G3G1
4E44.040.5236G6G7G8G9G5G2G4G10G3G1
1E16.510.3463G6G7G8G9G5G4G2G10G1G3
6E65.550.2192G6G7G8G9G5G4G10G2G1G3
2E26.35−0.6171G5G6G9G7G4G10G1G8G3G2
3E36.29−1.3967G5G9G6G1G10G4G7G3G8G2
E: Environments; IPCA score: Score based on the first interaction principal component axes.
Table 9. Mean values for the grain yield of wheat expressed in gm−2 at 13% moisture and interaction PCA1 values of the AMMI model of ten winter wheat varieties grown in six environments.
Table 9. Mean values for the grain yield of wheat expressed in gm−2 at 13% moisture and interaction PCA1 values of the AMMI model of ten winter wheat varieties grown in six environments.
Grain Yield Per m2 (g)
Genotypes (G)
EG1G2G3G4G5G6G7G8G9G10EmIPCAe1V
E1580.90780.70704.00785.90816.50762.50749.20819.90623.50663.70728.74.475520,653
E2671.70631.00647.60740.80837.20767.00701.50679.60736.90696.00710.9−6.573117,808
E3772.10502.10590.10685.60822.30804.80703.50580.30863.50731.70705.6−14.63830,580
E4448.00384.10322.60375.70367.90632.60631.50535.80538.90423.30466.04.532334,238
E5638.10708.10599.50648.60610.40865.30889.40848.80715.10648.40717.210.001015,810
E6566.70573.20524.20592.90616.40731.90714.30675.40640.20578.90621.42.203015,536
IPCAg1−5.88097.8640−0.1107−1.9985−9.10382.01137.101110.4235−6.4573−3.8487658.30.000030,483
Gm612.90596.50564.70638.20678.40760.70731.60690.00686.30623.70
Gm: Genotype mean; Em: Environmental means; IPCAe1: The first interaction principal component axes for environment; IPCAg1: Interaction principal component axes for genotype; E: Environmental labels with control (K) and 25 or 50 tha−1 phosphogypsum applied in both seasons; V: Variance.
Table 10. AMMI analysis of variance for the grain yield of ten winter wheat varieties in gm−2 at 13% moisture examined across six environments.
Table 10. AMMI analysis of variance for the grain yield of ten winter wheat varieties in gm−2 at 13% moisture examined across six environments.
Source 1dfSSMSF-ValueF-probThe Share of Total Variation %
Total1795,456,44130,483**
Treatments593,351,52856,8063.60 **0.000061.40
Genotypes9617,25068,5834.35 **0.000118.45
Environments51,552,340310,4689.26 **0.000046.30
Block12402,21533,5182.13 *0.020812.00
Interactions451,181,93726,2651.67 **0.016835.20
IPCA113487,08837,4682.38 **0.007741.21
IPCA211346,34931,4862.00 *0.035429.30
Residuals21348,50116,5951.050.409810.40
Error1081,702,69815,766**-
1 All sources were tested in relation to the error; *: Significant at p < 0.05 probability level; **: Highly significant at p < 0.01 probability level; df: Degree of freedom; F: F value calculated; IPCA1: The first interaction principal components axes; IPCA2: The second interaction principal components axes.
Table 11. Table of first AMMI selections per environment.
Table 11. Table of first AMMI selections per environment.
No.EMeanIPCAGenotypic Rank
(g)Score12345678910
5E5717.2010.0010G7G6G8G9G2G4G10G1G5G3
4E4466.004.5320G6G7G9G8G1G10G2G4G5G3
1E1728.704.4760G8G5G4G2G6G7G3G10G9G1
6E6621.402.2030G6G7G8G9G5G4G10G2G1G3
2E2710.90−6.5730G5G6G4G9G7G10G8G1G3G2
3E3705.60−14.6390G9G5G6G1G10G7G4G3G8G2
E: Environments; IPCA score: Score based on the first interaction principal component axes.
Table 12. Mean values for harvest index (%) and interaction PCA1 values of the AMMI model of ten wheat varieties grown in six environments.
Table 12. Mean values for harvest index (%) and interaction PCA1 values of the AMMI model of ten wheat varieties grown in six environments.
Harvest Index (%)
Genotypes (G)
EG1G2G3G4G5G6G7G8G9G10EmIPCAe1V
E10.510.500.520.540.510.510.490.530.510.540.510.13030.001
E20.510.460.480.510.510.520.490.500.510.530.500.03420.003
E30.500.440.460.480.500.520.490.470.510.510.49−0.03950.003
E40.460.270.290.320.440.510.450.310.480.420.40−0.42410.021
E50.480.470.490.510.480.480.460.500.480.520.490.15920.002
E60.490.480.500.520.490.490.480.510.490.530.500.13990.002
PCAg1−0.1270.1880.18670.1703−0.0910−0.2039−0.14330.1728−0.1680.01530.48 0.007
Gm0.490.440.460.480.490.510.480.470.490.51
Gm: Genotype mean; Em: Environmental means; IPCAe1: The first interaction principal component axes for environment; IPCAg1: Interaction principal component axes for genotype; E: Environmental labels with control (K) and 25 or 50 tha−1 phosphogypsum applied in both seasons; V: Variance.
Table 13. AMMI analysis of variance for the harvest index (%) of ten winter wheat varieties examined across six environments.
Table 13. AMMI analysis of variance for the harvest index (%) of ten winter wheat varieties examined across six environments.
Source 1dfSSMSF-ValueF-probThe Share of Total Variation %
Total1791.1960.00668**
Treatments590.57940.009822.59 **0.0000148.44
Genotypes90.07790.008662.28 *0.0219913.44
Environments50.27320.054633.17 **0.0104347.15
Block120.20690.017244.55 **0.0000135.71
Interactions450.22830.005071.35 *0.1129339.40
IPCA 1130.17930.013793.64 **0.0000978.54
IPCA 2110.02810.002550.670.7614712.31
Residuals210.02090.001000.260.999533.61
Error1080.40970.00379**
1 All sources were tested in relation to the error; *: Significant at p < 0.05 probability level; **: Highly significant at p < 0.01 probability level; df: Degree of freedom; F: F value calculated; IPCA: Interaction principal components axes.
Table 14. Table of first AMMI selections per environment.
Table 14. Table of first AMMI selections per environment.
NumberEMeanIPCAGenotypic Rank
(%)Score12345678910
5E50.48830.1592G10G8G4G6G7G2G3G1G9G5
6E60.49730.1399G10G4G8G6G3G9G1G5G2G7
1E10.51470.1303G10G4G8G6G3G9G1G5G7G2
2E20.50070.0343G9G5G1G10G4G3G6G8G7G2
3E30.4873−0.0396G9G5G1G10G6G4G3G7G8G2
4E40.396−0.4241G6G9G7G1G5G10G4G8G3G2
E: Environments; IPCA score: Score based on the first interaction principal component axes.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ljubičić, N.; Popović, V.; Ćirić, V.; Kostić, M.; Ivošević, B.; Popović, D.; Pandžić, M.; El Musafah, S.; Janković, S. Multivariate Interaction Analysis of Winter Wheat Grown in Environment of Limited Soil Conditions. Plants 2021, 10, 604. https://doi.org/10.3390/plants10030604

AMA Style

Ljubičić N, Popović V, Ćirić V, Kostić M, Ivošević B, Popović D, Pandžić M, El Musafah S, Janković S. Multivariate Interaction Analysis of Winter Wheat Grown in Environment of Limited Soil Conditions. Plants. 2021; 10(3):604. https://doi.org/10.3390/plants10030604

Chicago/Turabian Style

Ljubičić, Nataša, Vera Popović, Vladimir Ćirić, Marko Kostić, Bojana Ivošević, Dragana Popović, Miloš Pandžić, Seddiq El Musafah, and Snežana Janković. 2021. "Multivariate Interaction Analysis of Winter Wheat Grown in Environment of Limited Soil Conditions" Plants 10, no. 3: 604. https://doi.org/10.3390/plants10030604

APA Style

Ljubičić, N., Popović, V., Ćirić, V., Kostić, M., Ivošević, B., Popović, D., Pandžić, M., El Musafah, S., & Janković, S. (2021). Multivariate Interaction Analysis of Winter Wheat Grown in Environment of Limited Soil Conditions. Plants, 10(3), 604. https://doi.org/10.3390/plants10030604

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