Performance of grain sorghum hybrids under drought stress using GGE biplot analyses

The objective of this study was to estimate the adaptability and stability of grain sorghum hybrids grown under postflowering water stress and non-stress conditions. The trials were carried out in Nova Porteirinha-MG during the season of 2014 and 2015, and in Teresina-PI in the 2014 season. Twenty-nine-grain sorghum hybrids were evaluated, in a randomized complete block design, with three replications. Plots consisted of four lines with 3 m long. The grain yield data were submitted to the individual variance analysis, having considered the effects of the hybrids as fixed and the other effects as random. The joint analysis was carried out, and when the interaction genotypes x environments was significant, the grain yield data were submitted to the adaptability and stability analysis by the GGE biplot 2 P.S.C. Batista et al. Genetics and Molecular Research 16 (3): gmr16039761 method. A substantial reduction in the grain yield in environments with water stress was found. The highest yielding hybrids under water stress conditions in Nova Porteirinha-MG were 50A50, AG1080, AG1090, DKB550, DKB590, Jade, and BM737, and the highest yielding hybrids under the water stress in Teresina-PI were 1G282, 1G244, and A9721R. Considering all environments, the highest yielding hybrids were 1G282, DKB540, A9721R, 1G100, and AG1090.


INTRODUCTION
Abiotic stresses are among the main constraints to the development of agriculture in the world.Among the abiotic stresses, water deficiency is the most limiting crop performance (Araus et al., 2002), with a direct effect on yield reduction, causing damages in the physiological and metabolic processes in all phases of plant development (Taiz and Zeiger, 2009).
The use of tolerant crops such as sorghum [Sorghum bicolor (L.) Moench] consists in one of the main strategies that can reduce the yield losses caused by water stress.The good adaptation of sorghum in environments with low water availability is associated with the xerophytic features and to efficient drought tolerance mechanism that this crop owns (Landau and Sans, 2012).Sorghum stands out among other cereals of economic importance like corn by the possibility of cultivation in areas without irrigation, in seasons and locals with the occurrence of erratic rainfall distribution (Xin et al., 2009).In Brazil, sorghum has been grown as a succession crop, i.e., second season after soybean.The sowing occurs after February when rain season is diminishing, and post-flowering water stress is widespread.
Even being one of the most drought-tolerant crops, under marked water stress the sorghum plant can suffer damages in all the development phases, and the reproductive phase is the most affected by the stress, reducing weight and number of grains (Subudhi et al., 2000;Menezes et al., 2015).The impacts caused by water stress in sorghum can be mitigated by genetic improvement for drought tolerance, exploiting the genetic variability that this species presents (Rosenow and Dahlberg, 2000).
For an accurate selection of drought tolerance, it is essential the evaluation of the genotypes in different environments, allowing the precise identification of those most adapted and stable.Thus, according to Menezes et al. (2015), the interaction genotypes x environments (GxE) plays an important role in the phenotypic expression of the drought tolerance.
Several methods, based on different concepts, have been described for evaluation of the GxE interaction and the determination of the phenotypic adaptability and stability.The methods of Lin and Binns (1988), AMMI (Zobel et al., 1988), Annicchiarico (1992), and GGE biplot (Yan et al., 2000) are most frequently used.Among these methods, both AMMI and GGE biplot, which utilized multivariate techniques associated with descriptive analyses, in combination or substitution of the univariate analyses, make easier the observation of complex interactions (Ma et al., 2004).
The GGE biplot method takes into consideration the effect of the genotype and the GxE interaction.In this method, only the main effect of the genotype and the GxE interaction are important.Thus, both must be considered simultaneously.So, it is considered that the main effect of environment does not present relevance in the selection of genotypes.The GGE biplot model does not separate the genotype from the GxE interaction, maintaining them together in two multiplicative terms (Yan et al., 2000).
In this way, the objective of this study was to check the adaptability and stability through the GGE biplot analysis of grain sorghum hybrids in environments with and without water stress.

MATERIAL AND METHODS
Three trials were carried out to evaluate the performance of the hybrids.Two at the Experimental Station of Embrapa Maize and Sorghum, in Nova Porteirinha-MG, in the seasons of 2014 and 2015, and one at the Experimental Station of the Embrapa Mid-North, in Teresina-PI, in the season of 2014.Both sites have a very well-defined rain season, with no rain from May to October.In the case of sorghum drought selection, the main target has been to find stress tolerance during grain-filling period.
Nova Porteirinha is situated in the mesoregion of the North of Minas Gerais, considered as a semi-arid area.The geographical coordinates are 15°48'S latitude and 43°18'W longitude.The climate, according to Köppen, is of the type Aw (tropical with dry winter).The soil of the experimental area is characterized as medium-textured Red-Yellow Latosol.
Teresina presents the geographic coordinates of 05°05'S latitude and 42°48'W longitude.The climate, according to the classification of Thornthwaite and Mather, is C 1 sA'a', characterized as dry sub-humid, mega-thermal, with moderate water surplus in the summer.Teresina is located in a semi-arid area.The soil of the experimental area is a sandy loamtextured Dystrophic Yellow Argisol.
The maximum and minimum temperatures and the rainfall during the grown of the trials are presented in Figures 1 and 2. Twenty-nine commercial grain sorghum hybrids, belonging to the main seed companies of Brazil and among the most planted cultivars during the off-season in the Bioma Cerrado, were subjected to two moisture regimes.In treatment-1, the hybrids were wellwatered throughout the growing period to give full potential of the genotypes under nonstress conditions.In treatment-2, the cultivars received adequate watering from germination to boot stage (just before the flowering stage) after which no more watering was applied.This treatment simulated the post-flowering (terminal) moisture stress situation.Each combination of local and water regimes was considered one environment, totaling six environments.
Trials were irrigated using sprinklers spaced 12.0 x 12.0 m apart, operating pressure of 250 kPa, nozzles of 4.0 x 2.6 mm in diameter, and flow of 1.6 m 3 /h.The irrigation management was performed based on the evapotranspiration of the crop and the irrigation summed to the rainfall in the trials with water stress ranged from 334.3 to 400.0 mm and in the trials without water stress ranged from 462.7 to 600.0 mm.
An experimental design of randomized complete blocks was used, with three replications and plots of four rows of 3 m in length and 0.5 m between rows.The two central rows of each plot were harvested.Sowing was carried out with excess seeds and, after thinning at 30 days after sowing, the plant density was 180,000 plants/ha.
The crop management consisted of two hand hoeings and applications of insecticides for the control of armyworm.Soon after flowering, the panicles were covered with polyethylene bags to prevent birds damaging.Grain yield was evaluated by weighing the grain mass, corrected to 13% moisture, transforming the results to kg/ha.
Data were submitted to the individual variance analysis, considering the effect of the hybrids as fixed and the other effects as random.Since it was detected that the ratio between the greatest and the smallest mean square of the individual variance analysis did not exceed the 7:1 ratio, the joint analysis of the trials was performed (Banzatto and Kronka, 2006).After that, the data were submitted to the adaptability and stability analysis using the GGE biplot method (Yan et al., 2000).
The GGE Biplot model used was: where Y ij is the grain yield average of genotype in environment j; μ is the general mean of the observations; β j is the principal effect of environment; y 1 and y 2 are the scores associated with the first (PC1) and second principal component (PC2), respectively; ε 1 and ε 2 are the values of PC1 and PC2 of the genotype i, respectively; ρ j1 and ρ j2 are the values of PC1 and PC2 for the environment j, respectively; and ε ij is the error associated with the model of the i-th genotype and j-th environment (Yan et al., 2000).The analysis was developed using the GGEGui package implemented in the R software (R Development Core Team, 2016).

RESULTS AND DISCUSSION
There were significant differences between the sources of hybrids (G), environments (E) and interaction GxE (Table 1).The results showed that the hybrids performed differently under the two water regimes, and the significant interaction GxE made more difficult the selection of a hybrid that performed well in both conditions.The coefficient of variation was below 20%, and it is considered satisfactory for field trials when yield is measured (Pimentel-Gomes, 2009).The general mean considering the six trials was 4323 kg/ha (Table 2), superior to the national mean achieved in the 2014 season, which was 2587 kg/ha (CONAB, 2015).The average grain yield of the hybrids was inferior to the national mean only in Nova Porteirinha-MG, in 2015, and in Teresina-PI, in the trials with water stress (Table 2).total variation of the yield performance of the hybrids, plus the interaction with environments (G+GxE).2.
In Figure 3, there is a formation of a polygon by the connection of the vertices farthest from the origin of the biplot, being formed by the hybrids DKB 540 (18), Jade (21), FOX (26), 1G 220 (5), BRS 332 (3), A 9735 R (28), 1G 100 (4), and 1G 282 (7).All the other hybrids are situated inside that polygon.The graphic was divided into eight sectors by the vectors coming from the biplot center (0; 0), perpendicular to the polygon sides.
When different hybrids are adapted to groups of different environments, the variation between groups is greater than within the group, and there is the formation of a mega-environment (Yan and Kang, 2003).In biplot, the mega-environments are the sectors that contain one or more environments.Thus, in the present study, there was the formation of two mega-environments.In mega-environment 1, the environments from Nova Porteirinha were grouped under water stress in 2014 and 2015 (NC14 and NC15).In mega-environment 2, the environments from Teresina-PI were grouped under water stress (TC14) and non-stress conditions (TS14), and from Nova Porteirinha under non-stress conditions in 2014 (NS14) and 2015 (NS15).
The hybrids situated on the vertices of each sector present the best or worst performance in that sector (Yan and Tinker, 2006).Hybrid DKB 540 (18) was the vertex of mega-environment 1.Thus, it was the most adapted to the environments of Nova Porteirinha under water stress in 2014 and 2015, presenting greater yield in those environments.Hybrids 50A50 (9), AG 1080 (13), AG 1090 ( 14), DKB 550 ( 19), DKB 590 (20), Jade (21), and BM 737 (23) presented good adaptation to these environments.Thus, these are the hybrids that presented higher yields under water stress in Nova Porteirinha.
For the other environments, hybrid 1G 282 (7) was the most adapted because it was the vertex of mega-environment 2 and hybrids 1G 244 (6) and A 9721 R ( 27) also showed good adaptation to these environments.As the TC14 and NS15 environments were positioned on the division of sectors seven and eight, the hybrids belonging to the two sectors presented good adaptation to these environments (Figure 3).So, the hybrids with higher yields under water stress conditions in Nova Porteirinha were DKB 540, 50A50, AG 1080, AG 1090, DKB 550, DKB 590, Jade, and BM 737 and the hybrids with higher yields under stress conditions in Teresina-PI were 1G 282, 1G 244, and A 9721 R.
Even the highest yielding hybrids presented a reduction in grain yield when comparing trials under stress to those under non-stress trials.The reduction in grain yield in Nova Porteirinha in 2014 for hybrids 50A50, AG 1080, AG 1090, DKB 550, DKB 590, and BM 737 was of 73.4, 66.6, 59.2, 69.3, 68.5, and 48.7%, respectively.In Nova Porteirinha, in the year of 2015, the reduction in grain yield for hybrids 50A50, AG 1080, AG 1090, DKB 550, DKB 590, and BM 737 was of 73.4, 66.6, 59.2, 69.3, 68.5, and 48.7%, respectively.Hybrid Jade, in Nova Porteirinha, in the year of 2014, presented equivalent yield in the environments with and without water stress.However, in 2015, the reduction in the water stress was 54.1%.For the highest yielding hybrids in Teresina; 1G 282, 1G 244, and A 9721 R, the reduction in grain yield under the water stress condition was of 40.9, 56.9, and 31.4%,respectively (Table 2).
The hybrids grouped in sectors that do not contain environments are not among the highest yielders in any of the environments, being considered unfavorable to the evaluated locals (Karimizadeh et al., 2013).The selection of the best hybrids, based on their performance, allows a safer recommendation to the farmer concerning planting seasons and locals where water stress can occur, like the off-season in Cerrado and the Northwestern of Brazil.Another advantage is the selection of the hybrids most responsive to a possible water increasing.
In the "mean vs stability" biplot (Figure 4), it is possible to obtain hybrids with high yield and high stability using a mega-environment, making it possible, in that way, to indicate the best hybrids for each particular environment (Yan et al., 2007).
In Figure 4, a straight line was drawn with an arrow passing through the origin of the biplot.The direction of the arrow shows the highest grain yield, with the sorghum hybrids being ranked along the line.The greater the projection distance of the hybrid regarding this straight line, the lower its stability is (Yan, 2011).Therefore, the highest yielding hybrids with mean above the general mean of all the environments were 1G 282 (7), DKB 540 (18), A 9721 R (27), 1G 100 (4), and AG 1090 (14), presenting average yield in all environments of 5716, 5514, 5243, 5193, and 4916 kg/ha, respectively (Table 2).
In Figure 5, it is possible to observe the formation of concentric circles around the average grain yield, and the closer to the center of the concentric circles the more useful the hybrid is; what means that it presents the best average performance.When estimating the adaptability and stability, this is very important because it is not worth if a hybrid is stable but has a low yield (Yan, 2011).The most useful hybrid was 1G282 (7), followed by hybrids A9721R (27), 1G100 (4), and DKB540 (18).The high stability associated with a high yielding performance is essential for the recommendation of a hybrid, principally to environments or seasons with high rainfall instability, as in the off-season cultivation in Brazil.Vectors with origin on the center of biplot link the environments.According to Yang et al. (2009), the size of the vectors indicates the stability of the environments.So, environments with smaller vectors are more stable.Thus, NC15 was the environment that less contributed to the GxE interaction, presenting the smallest yield difference among the evaluated hybrids.Furthermore, NC14 was the environment that presented the largest vector, thus, the environment that most contributed to the GxE interaction.2.
The changes in grain yield stability in Nova Porteirinha in the environments under water stress in the seasons 2014 and 2015 can be explained by weather conditions, which varied from a year to another.In 2014, planting was carried out in May, with higher temperatures at the beginning of plant development, but at flowering time and grain filling, the temperatures were lower.In 2015, the planting was carried out in June, with milder temperatures during the initial development of the plants, but in the reproductive period of the plants, temperatures were highest, maximizing the effects of water stress on the plants (Figure 1).A different effect was observed in the same place in the non-stress environments, with smaller variation in the grain yield stability between the two periods.The normal water supplying minimized the effect of the high temperatures observed in 2015.

Figure 3 .
Figure 3. Sectors and mega-environments obtained by the GGE biplot model for grain yield of twenty-nine sorghum hybrids grown under water stress and non-stress conditions in Nova Porteirinha-MG, in 2014 and 2015 seasons, and in Teresina-PI, in the 2014 season.For abbreviations, see legend to Table2.

Figure 4 .
Figure 4. Mean versus stability according to the GGE biplot model for grain yield of twenty-nine sorghum hybrids grown under water stress and non-stress conditions in Nova Porteirinha-MG, in 2014 and 2015 seasons, and in Teresina-PI, in the 2014 season.For abbreviations, see legend to Table2.

Figure 5 .
Figure 5. Classification of twenty-nine-grain sorghum genotypes according to the GGE biplot model for grain yield of 29 sorghum hybrids grown under water stress and non-stress conditions in Nova Porteirinha-MG, in 2014 and 2015 seasons, and in Teresina-PI, in the 2014 season.For abbreviations, see legend to Table2.

Figure 6
Figure6illustrates the relationship between grain yield and environmental stability.Vectors with origin on the center of biplot link the environments.According toYang et al. (2009), the size of the vectors indicates the stability of the environments.So, environments with smaller vectors are more stable.Thus, NC15 was the environment that less contributed to the GxE interaction, presenting the smallest yield difference among the evaluated hybrids.Furthermore, NC14 was the environment that presented the largest vector, thus, the environment that most contributed to the GxE interaction.

Figure 6 .
Figure 6.Relationship between the environments obtained by the GGE biplot model for grain yield of twenty-nine sorghum hybrids grown under water stress and non-stress conditions in Nova Porteirinha-MG, in 2014 and 2015 seasons, and in Teresina-PI, in the 2014 season.For abbreviations, see legend to Table2.

Table 1 .
Summary of the joint analysis of variance concerning grain yield of twenty-nine-grain sorghum hybrids grown under water stress and non-stress conditions in Nova Porteirinha-MG, in 2014 and 2015 seasons, and in Teresina-PI, in the 2014 season.*Significant at 1% of probability.*Significant at 5% of probability.ns Nonsignificant by the F-test.d.f.= degrees of freedom, CV = coefficient of variation. *