Grain-yield stability among tropical maize hybrids derived from doubled-haploid inbred lines under random drought stress and optimum moisture conditions

Drought is a devastating environmental stress in agriculture and hence a common target of plant breeding. A review of breeding progress on drought tolerance shows that, to a certain extent, selection for high yield in stress-free conditions indirectly improves yield in water-limiting conditions. The objectives of this study were to (i) assess the genotype × environment (GE) interaction for grain yield (GY) and other agronomic traits for maize (Zea mays L.) across East African agro-ecologies; and (ii) evaluate agronomic performance and stability in Uganda and Tanzania under optimum and random drought conditions. Data were recorded for major agronomic traits. Genotype main effect plus GE (GGE) biplot analysis was used to assess the stability of varieties within various environments and across environments. Combined analysis of variance across optimum moisture and random drought environments indicated that locations, mean-squares for genotypes and GE were significant for most measured traits. The best hybrids, CKDHH1097 and CKDHH1090, gave GY advantages of 23% and 43%, respectively, over the commercial hybrid varieties under both optimum-moisture and random-drought conditions. Across environments, genotypic variance was less than the GE variance for GY. The hybrids derived from doubled-haploid inbred lines produced higher GY and possessed acceptable agronomic traits compared with the commercial hybrids. Hybrid CKDHH1098 ranked second-best under optimum-moisture and drought-stress environments and was the most stable with broad adaptation to both environments. Use of the best doubled-haploids lines in testcross hybrids make-up, well targeted to the production environments, could boost maize production among farmers in East Africa.


Introduction
In East Africa, maize (Zea mays L.) is widely used as a major staple. Total maize production for the last 10 years was 236.65 Mt harvested from a total area of 147.17 Mha (FAOSTAT 2014). More than half of the maize produced in East Africa is traded in the commercial maize market, and it makes a major contribution to the economy. Although maize is considered a non-traditional export crop in Africa, its exports from East Africa amounted to 4.71 Mt, valued at US$1.72 million in the last 10 years (FAOSTAT 2014). However, most farmers are unable to attain the potential yield of hybrid maize, owing to biotic and abiotic factors. Among the abiotic factors, drought stress is ranked number one.
In the last decade in Africa, drought stress resulting from climate change has become a major constraint that has gained prominent attention because of its effects on maize productivity (Msowoya et al. 2016). In East Africa alone, 34 drought events were reported, affecting >67.8 million people (EM-DAT 2016). Hence, drought remains a major catastrophe causing much stress and suffering to humanity.
Most places in the maize-producing regions of East Africa experience frequent droughts that often coincide with the flowering period of the maize crop, leading to poor grain yields or total crop failure (Omoyo et al. 2015). Drought has therefore become a major threat to food security for smallholder farmers who grow maize in drought-prone areas of East Africa. The most economically viable and sustainable option for alleviating the situation is breeding and deploying improved drought-tolerant, high-yielding maize cultivars for the farmer to assure profitable yields even in drought years (Cairns et al. 2013). This has called for the identification of different sources of droughttolerance traits and an efficient means of generating new products to cope with the situation. One of the ways is the use of doubled-haploid technology in rapid-development, droughttolerant lines. Application of doubled-haploid technology has made possible the development of inbred lines in one or two generations (Prigge et al. 2011) compared with classical pedigree methods that produce 96.9% homozygous lines after six to 10 generations of selfing heterozygous material (Hallauer et al. 2010). This, in turn, has enabled rapid development of hybrids from doubled-haploid lines that are drought tolerant. These hybrids will improve food security for smallholder farmers in many areas of Africa.
The identification of genotypes with a high yield potential, coupled with wide adaptability and stability, is a key target of any maize-breeding program (Mendes et al. 2012). The major drawback, however, in the selection of genotypes with high yield potential in different environments is the genotype Â environment (GE) interaction, because strong interactions can hamper the selection process. Thus, genotypes that perform well in one environment may not do as well in another (Mendes et al. 2012). GE interactions have been investigated through the use of statistical tools such as additive main effects and multiplicative interaction (AMMI) analysis for grain yield and grain micronutrient concentrations and stability (Oikeh et al. 2004;Gauch 2006); and genotype main effect plus GE interactions (GGE) (Yan 2001) for the analyses of grain yield and stability in tropical maize. Several studies have been undertaken to assess the performance of genotypes under different environmental conditions (Beyene et al. 2013). However, most of these studies were carried out in the same country. Collaborative efforts among research organisations in five countries, Kenya, Mozambique, South Africa, Tanzania and Uganda, have developed quite a diverse germplasm base through the Water Efficient Maize for Africa (WEMA) project. The WEMA project is a public-private partnership established to develop drought-tolerant and insect-protected maize by using conventional breeding, marker-assisted breeding and biotechnology, with a goal to make these varieties available royalty-free to smallholder farmers in Sub-Saharan Africa through African seed companies (Oikeh et al. 2014;Edge et al. 2018). Materials developed are tested in multiple environments that represent drought-stress and nonstress (optimum-moisture) locations to identify high-yielding and adapted varieties for release and cultivation in the respective countries. With the formation of the East African Community, seed companies prefer licensing and marketing maize seed for varieties that have been commercialised in more than one country. The tremendous progress made in WEMA product developments has led to gains in grain yield performance as reported in Maize Regional Trials (Sserumaga et al. 2016;Beyene et al. 2017). With the changing climatic conditions, information and utilisation of knowledge on the genetic performance of different germplasm are paramount in a breeding program. The WEMA project extensively used doubled-haploid technology in inbred line development, which constituted most of the maize hybrids that were to be tested in multiple environments among WEMA countries before commercialisation.
The objectives of this study were to (i) assess GE interaction for grain yield and other agronomic traits in maize across East African agro-ecologies; and (ii) evaluate agronomic performance and stability in Uganda and Tanzania under optimum and random drought conditions.

Genetic materials and performance evaluation
In total, 43 hybrids derived from doubled-haploid inbred lines were selected for this study (see Supplementary materials table 1, available at the journal's website). The 43 hybrids were selected based on grain yield, disease resistance and other agronomic traits from previous preliminary yield testing. Sufficient seed of the 43 hybrids was produced in 2013 at the Kenya Agricultural and Livestock Research Organisation (KARLO) Maize Research Station at Kiboko, Kenya. Each female parent was planted in five rows of 5 m length, and the male donor plants were planted at two different times (-5 and 0 days) to effect nicking at pollination and synchronise flowering. All recommended agronomic practices for maize production for the agroecologies were applied (Sserumaga et al. 2016). Immediately before flowering, all of the ears of the female plants were covered with shoot bags. During hand pollination, pollen was collected and bulked from the male plants when 20% of the males had started to shed pollen. The 43 hybrids along with seven checks were evaluated across eight sites in Uganda and five sites in Tanzania (Supplementary materials table 2). Optimum-moisture sites were selected based on the total amount of precipitation and the trials were established during the main season, whereas at the drought-stress sites, trials were established so that the stress coincided with 2-3 weeks before 50% flowering in order to have sufficient drought stress for evaluation of the materials. The experimental design was a 5 Â 10 a-lattice with two replications at each location. An experimental unit was a two-row plot, 5 m long, spaced 0.75 m between rows and 0.25 m between plants. Two seeds were planted per hill and subsequently thinned to one plant per hill at 4 weeks after emergence, to give a final plant population density of 53 333 plants ha -1 . In all the experiments, standard agronomic and cultural practices including weeding and appropriate fertiliser applications were followed.

Data collection
Under both drought-stress and optimum-moisture conditions, data collected included: days to anthesis (days from planting to when 50% of plants had shed pollen) and days to silking (days from planting to when 50% of plants had extruded silks); anthesis-silking interval (determined as the difference between days to silking and days to anthesis); plant height (measured in cm as the distance from the base of the plant to the height of the first tassel branch); number of ears per plant (determined by dividing the total number of ears produced per plot by the number of plants harvested per plot); husk cover (obtained by dividing the number of ears with poor husk cover by the number of plants harvested per plot); ear aspect (rated on a scale of 1-5 where 1 is uniform ears with the preferred texture and 5 is ears with the undesirable texture); plant aspect (PA) (rated on a scale of 1-5 where 1 is short plant with uniform and short ear placement and 5 is tall plants with high ear placement ear position); and grain moisture. All of the ears harvested from each plot were weighed and representative samples of ears were shelled to determine the percentage moisture of the grain, using a Dickey-John moisture meter at all locations. Grain yield (t ha -1 ) was calculated from ear weight based on a shelling percentage of 80% and grain moisture content of 12.5%.

Statistical analyses Analysis of variance and GE interaction under drought and optimum conditions
Individual analyses of variance (ANOVAs) were carried out on the data for each genotype. For the joint analysis of all datasets, we used a mathematical model that considered sites and progenies as random effects and was equivalent to the following equation described by Cruz and Regazzi (1994): where Y ijk is observed value of the ith progeny of the jth environment in the kth replications; m is general mean; G is effect of the ith genotype (i = 1, 2,. . .i); A is effect of the jth environment (j = 1, 2,. . .j); GA is effects of the interaction of the ith progeny with the jth environment; B/A jk is effect of the kth block within the jth environment; and e ijk is random error. ANOVA for all traits was done separately for each environment and combined across locations by using the PROC MIXED procedure from SAS 9.3 (SAS Institute 2011). For the combined analysis, variances were partitioned into relevant sources of variation to test for differences among genotypes and the presence of GE interactions. In the across-environment ANOVA, genotype effects were tested for significance by using the corresponding interactions with the environment as the error term; the GE interaction was tested by using the pooled error.

Genotype performance under different conditions
Individual locations or across analysis were computed for grain yield and other traits by using the mixed-model analysis in META-R to generate a best linear unbiased estimate (BLUE) for all genotypes (Alvarado et al. 2015). For comparing entries evaluated in different locations, the entry means were expressed as a percentage of the average performance of the best check hybrid at the respective locations.

Genetic variances and heritability under different conditions
Estimates of genotypic (s 2 G ), location (s 2 L ), genotype Â location (s 2 GÂL ) and error (s 2 E ) variance were calculated using PROC MIXED (option = REML) of SAS. Broad-sense heritability (H 2 ) was calculated as the proportion of genetic variance over the total phenotypic variance, and for individual trials was estimated according to Hallauer et al. (2010): where s 2 G is the genotypic variance, s 2 E is the error variance, and r the number of replications.
Broad-sense heritability for traits across environments was estimated by using the variance components according to Hallauer et al. (2010) as: where s 2 G , s 2 GÂL and s 2 E are genotypic, genotype Â location and residual variance components, respectively; E is the number of environments; and R is the number of replications. Genotypic correlations (r g ) among locations were estimated according to Cooper et al. (1996) as: where r p(1,2) is the phenotypic correlation between the traits measured in locations 1 and 2; and H 1 and H 2 are the broad-sense heritabilities for the traits measured in locations 1 and 2, respectively.

Genetic correlations among test locations
The genetic correlation among the pairs of environments for each trait under study was obtained as suggested by Yamada (1962), using the following expression: where r g is coefficient of correlation between the two locations for a certain trait, s g1 is genetic variance at location 1, s g2 is genetic variance at location 2, s 2 g is joint genetic variance of the joint analysis, and s 2 ge is variation of the GE interaction. Cluster analysis using Ward's minimum variance method (Ward 1963) was performed on group environments based on genetic correlations among the environments, using META-R (Alvarado et al. 2015).

GGE biplot analysis of grain yield response and stability
Adjusted data on the grain yield from ANOVA were subjected to GGE biplot analysis (Yan 2001) to determine grain yield stability and the pattern of response of genotypes evaluated across the environments. The analyses were done and biplots generated using the GGEbiplot software version 7 (www. ggebiplot.com/biplot). The GGE biplot Model 3 equation used was: where Y ij is the average yield of genotype i in environment j; Y J is the average yield across all genotypes in environment j; l 1 and l 2 are the singular values for PC1 and PC2, respectively; x i1 and x i2 are the principal component scores PC1 and PC2, respectively, for genotypei; h j1 and h j2 are the PC1 and PC2 scores, respectively for environment j; and e ij is the residual of the model associated with genotype i in environment j.

ANOVA and GE under drought and optimum moisture conditions
Analysis of variance across three locations revealed that environment was highly significant for all traits, whereas genotype was significant for grain yield, husk cover and ear aspect ( Table 1). The GE interaction was significant for ear aspect and husk cover. This implied differences in yield performance among the test materials, although there was no differential response at different locations.
Under conditions of optimum moisture, combined ANOVA across 10 locations revealed that genotype, environment, and GE interaction were highly significant (P < 0.0001) for grain yield (Table 1), suggesting various responses in yield performance among the test materials at different locations. Environment was significant for all the other traits, and genotype was significant for most traits except anthesis-silking interval, ear position and ears per plant (Table 1). The GE interaction was significant only for plant height, grain moisture, ear aspect and plant aspect (Table 1).
Combined ANOVA across all the 13 locations showed that genotype, environment, and GE interaction were all significant for grain yield (Table 1). In addition, environment was highly significant (P < 0.001) for all the traits and genotype was significant for all traits except anthesis-silking interval, ear position and ears per plant (Table 1). The GE interaction was significant for grain yield (P < 0.01), husk cover (P < 0.05), and ear aspect (P < 0.001) ( Table 1). This meant that there were significant differences in yield performance among the test genotypes at different locations.

Hybrid performance under drought-stress and optimummoisture conditions
Genotype performance across drought-stress environments varied with a grand mean of 5.6 t ha -1 . The grain yield of the test hybrids varied from 4.7 t ha -1 for CKDHH1134 (G24) to 7.0 t ha -1 for CKDHH1090 (G14) ( Table 2). The highest yielding testcross hybrid (CKDHH1090) had a 44.2% yield advantage over the best commercial hybrid, Com Check 3 (G48). There were variable responses to phenology, with days to anthesis ranging from 61 to 68.6 days, although the anthesis-silking interval values were comparable (1.8-3.4 days) among all testcrosses and the commercial checks (Table 2). Heritability estimates for the different traits were generally very low (0-0.29), except for PA, which had the highest heritability (Table 3).
Genotype performance varied significantly across optimum moisture environments with a grand mean of 6.7 t ha -1 ( Table 2). The grain yield of the test hybrids varied from 5.8 t ha -1 for CKDHH1076 (G7) to 7.6 t ha -1 for CKDHH1097 (G15). The highest yielding testcross (CKDHH1097) gave 23% yield advantage over the best popular commercial hybrid, Com Check 1 (G46). All genotypes had comparable maturity, with days to anthesis ranging from 61.2 to 64.2 days. Therefore, all testcross hybrids could be categorised as earlymaturing genotypes. Medium to high heritability estimates were found for most traits. The highest heritability of 0.68 was recorded for ear aspect, followed by moisture content (0.60), and the least was for plant aspect (H 2 = 0.13).
Average performance across all test environments showed a grand mean yield of 5.6 t ha -1 , but grain yield of the best performing, top 10 test hybrids varied from 4.7 t ha -1 for CKDHH10960 (G33) to 7.0 t ha -1 for CKDHH1090 (G14) (Supplementary materials table 3). The highest yielding testcross (CKDHH1090) had a 42.9% grain yield advantage over the best popular commercial hybrid, Com Check 3 (G48). The genotypes varied in days to phenological development with days to anthesis ranging from 61.0 to 68.6 days, and entry CZH0616 (G44) with 61 days was considered the earliest maturing genotype. Among the test materials, genotype CKDHH1098 (G16) had the best husk cover of 3.9. Also, heritability estimates were highest for ear aspect (H 2 = 0.76) ( Table 3).

Genetic variances and heritability under different environmental conditions
Estimates of genotypic, location and genotype Â location variances under random drought, optimum moisture conditions and across all test locations are presented in Table 3. Under drought, husk cover was the only factor where the genotypic variance was larger than the location variance. The results showed that environment accounted for 92.2% of the total variation in grain yield. The H 2 for grain yield was only 0.29 under drought conditions.
Similarly, under optimum moisture conditions, genotypic and genotype Â location variances were smaller than location for all traits. The genotypic variance was larger than the genotype Â location variance for days to anthesis alone. The results showed that location accounted for 41% of the total variation for grain yield. The H 2 for grain yield was moderate, 0.57 under optimum moisture condition.
Across locations, genotypic and genotype Â location variances were smaller than the location variance for all traits. The results indicated that genotypic and genotype Â location variances accounted for <5%, whereas location accounted for 69.2% of the total variation in grain yield. The H 2 for grain yield was 0.66.

Genetic correlations among test locations
The genetic correlations among locations were based on grain yield and were used for cluster analysis to classify the environments for their yield potential and stability. Under random drought conditions, genetic correlations among locations were lowest for Ilonga_DT vs Masaka (0.13) and highest for Masaka vs Makutupora_DT (0.6) (Supplementary materials table 4). Clustering based on genetic correlation for grain yield revealed two clusters at 83.64% (Fig. 1a). Cluster I consisted of Ilonga and Cluster II comprised Masaka and Makutupora (Fig. 1a).
Under optimum moisture conditions, genetic correlations among locations were lowest for Abii vs Ikulwe (-0.54) and the highest for Bulindi vs Ngetta (0.39) (Supplementary materials table 4). Clustering based on genetic correlation for grain yield revealed two main clusters at 50.1% (Fig. 1b). Cluster  I consisted of three locations that were separated into Namulonge and a subcluster into Bulindi and Ngetta in Uganda. Cluster II consisted of two major subclusters that were divided into subsubclusters. Across all environments, genetic correlations among locations ranged from -0.54 to +0.58 (Supplementary materials table 4). Clustering based on genetic correlation for grain yield revealed two main clusters at 45.0% (Fig. 1c). Cluster I consisted of two subclusters, the first clustering Masaka and Namulonge and the other Bulindi and Ngetta. Cluster II also had two subclusters, one with Karatu clustered with a sub-subcluster of Kabuku and Makutupora DT, and the other consisting of two subsubclusters that included Ilonga-DT and Ilonga-OPT and then Serere and Ikulwe, respectively.

GGE biplot analysis for grain yield response and stability of maize genotypes evaluated across different environmental conditions
The GGE biplot analysis was used to examine the performance of genotypes in the different environments under different conditions. The view of the polygon of the GGE biplot under drought indicated the best genotypes in each environment (Fig. 2a). The presence of two or more environments within a sector shows that a single genotype has the highest yield in those environments. If the environments fall into different sectors, different genotypes performed well in different environments. PC1 explained 53.1% of total variation, and PC2 explained 30.3%. Thus, these two axes accounted for 83.5% of the GGE variation for grain yield under drought (Fig. 2a). For optimum moisture conditions, PC1 explained 28.14% of total variation, whereas PC2 explained 17.88%. Thus, these two axes accounted for 46.02% of the GGE variation for grain yield (Fig. 2b).
The mean vs stability view of the GGE biplot of yield data under drought was used to assess the stability of the 50 genotypes (Fig. 2c). The genotypes were ranked along the average tester coordinate (ATC) axis (abscissa), with an arrow pointing to a greater value based on their mean performance across all environments. The line with the arrow discriminates or separates entries with below-average means from those with above-average means. The average yield of the genotypes is approximated by the projections of their markers on the average tester axis. In the GGE biplot analysis, the ATC approximates the genotype contribution to GE, which is a measure of their instability. The stability of the genotypes is measured by their projection on the ATC axis; the longer the genotype's projection, the less stable it is. Based on this, the top three stable genotypes under drought were CKDHH1081 (G9), CKDHH1148 (G22) and CKDHH1098 (G16) (Fig. 2c), and under optimum conditions, CKDHH0947 (G29), CKDHH1148 (G22) and CKDHH1097 (G15) (Fig. 3) were the most stable genotypes because they had a near zero projection onto the ATC axis. This implies that their ranking was highly consistent across locations. Conversely, genotype CKDHH1123 (G20) under drought stress, and CKDHH1075 (G6) under optimum moisture conditions, were the most unstable among the genotypes evaluated because they had longer projections away from the ATC axis than the other genotypes. Table 3. Variance decomposition and heritability for grain yield and agronomic traits of 43 testcross hybrids, two internal WEMA hybrid checks, three commercial hybrids and two local hybrid checks across drought, optimum conditions, and environments in Uganda and Tanzania GY, Grain yield; AD, days to anthesis; ASI, anthesis-silking interval; PH, plant height; EH, ear height; EP, ear position; EPP, no. of ears per plant; MOI, grain moisture; HC, husk cover; EA, ear aspect (scale of 1-5, where 1 is uniform cobs with the preferred texture and 5 is cobs with undesirable texture); PA, plant aspect.

Discussion
This study was conducted to examine the potential agronomic performance of drought-tolerant, doubled-haploids inbred lines developed and selected under managed drought stress in their different hybrid combinations across variable growing conditions in East Africa. There were significant differences for most traits, which demonstrated that selection could be made under both optimum conditions and random drought. Beyene et al. (2013), Adebayo and Menkir (2014) and Ertiro et al. (2017) reported differential responses of tropical maize hybrids under different environmental conditions. The non-existence of significant GE interaction effects under random drought suggested that testcross hybrids had consistent performance, which agrees with the work of Adebayo and Menkir (2014  indicating that most of the experimental hybrids were superior to, and had greater stability than, the commercial checks. These results were consistent with the results of other authors who reported better tolerance to drought among new varieties (Beyene et al. 2013;Adebayo and Menkir 2014;Ertiro et al. 2017). However, the average grain-yield reduction for the experimental hybrids under random drought in our study was only 16% compared with the yield under optimum moisture conditions, indicating that the drought stress was low to moderate. The doubled-haploid hybrids exhibited a wide-ranging variation in grain yield and other agronomic traits under both optimum-moisture and drought-stress environments. Similar observations were reported by Sserumaga et al. (2016) while examining performance and GE interactions of doubled-haploid hybrids in Uganda. This might imply that several factors could be affecting the performance of maize in the same or different environments.
Our results showed that the hybrids derived from doubledhaploid inbred lines outperformed the commercial hybrids for grain yield and other agronomic traits. Similar to the present study, Sserumaga et al. (2016) and Beyene et al. (2011) reported superiority in performance of doubled-haploid hybrids over the commercial checks in their studies. Under random drought stress, the highest yielding doubled haploid in the present work (CKDHH1090) had a 44.2% yield advantage over the best popular commercial hybrid, Com Check 3 (G48); however, under optimum environments, the highest yielding doubled haploid (CKDHH1097) had a 23% yield advantage over the best commercial hybrid, Com Check 1 (G46). This implied that the doubled-haploid hybrids were superior in performance to the commercially available hybrids. Therefore, the performance of the doubled-haploid hybrids indicated that the lines used in their development offered potential new sources for accelerating the breeding of high-yielding, drought-tolerant maize hybrids in similar tropical environments.
Under conditions of random drought and optimum moisture, the proportion of the total variance in grain yield attributed to the environment was high, whereas genotype and GE variances were relatively small. The broad-sense heritability for grain yield under drought was very low, which meant that selection for grain yield under this environmental setting led to low genetic gains. However, the broad-sense heritability for grain yield under optimum moisture conditions was higher (0.57) than under drought stress.
Genetic correlations among some locations were positive, indicating that the germplasm selected at any one of these locations could be used at the other, positively correlated locations. However, the reverse is true for those with negative correlations. The genetic correlations of grain yield under drought-stress and optimum-moisture environments were generally weak, which implied that grain yield was mediated mainly by the same set of genes that conditioned similar responses in the hybrids in the two contrasting environments. These results are similar to those reported for grain yield in maize grown under organic and conventional production systems (Lorenzana and Bernardo 2008). Although the correlation in that study was slightly higher than in our study, the results were different from those of Bänziger et al. (1997), who reported that different sets of genes controlled maize grain yield under low and high nitrogen environments. There was a low genetic correlation (r g = 0.04) between some of the pairs of locations, implying that these environments were very different. Burdon (1977) pointed out that locations with low genetic correlations between them should be treated separately.
The stability and performance of different doubled-haploid hybrids were assessed by using GGE biplots to identify the best entries at each location and assess their stability. Emphasis was mainly on high grain yield under optimum moisture and improved yields under random drought, coupled with stable performance across sites with acceptable secondary traits. The most stable varieties in both environmental conditions were all derived from the doubled-haploid homozygous lines. These lines therefore contained favourable genes with additive effects resulting in heterosis (hybrid vigour). The fixation of favourable alleles in parental lines of the hybrids that performed well across stress environments contributes to the superior performance of hybrids (Ertiro et al. 2017). Hybrid CKDHH1098 (G16) ranked second best across all environments (6.82 t ha -1 ) and was the most stable hybrid evaluated. Therefore, it has the potential to be grown in a wide range of drought-prone environments in East Africa because of its broad adaptation. Two of the highest yielding varieties, CKDHH1090 (G14) and CKDHH1075 (G6), were not among the most stable, suggesting that these varieties have specific adaptation to some environments. These results are consistent with those of Badu-Apraku et al. (2012) and Makumbi et al. (2015), who identified high-yielding but unstable varieties in different, contrasting environments. This means that higher average yield indicated a higher response to favourable environments, but could result in lower environmental stability.

Conclusion
This study showed that commercial varieties in the East African market are more vulnerable to drought and therefore less productive in various agro-ecologies in East Africa than the recent, new hybrids developed from the drought-tolerant doubled-haploid inbred lines through the WEMA public-private partnership. In terms of future projections, East Africa is more likely to be more vulnerable to drought if temperatures continue to increase. Results of the present study demonstrate that it is possible to have hybrids bred and released across East Africa and be recommended across different environmental conditions. This study also laid the foundation for exploiting GE interactions, not only to identify stable genotypes but also to classify environments into broader mega-environments, and to identify the most discriminating, high-yielding and stable environment for maize production in East Africa. Therefore, commercialisation of outstanding hybrids, such as CKDHH1098 (G16) identified in the present study, with high mean yield and stable performance across contrasting management conditions would contribute to enhancing maize productivity and yield stability for smallholder farmers.

Conflicts of interest
Authors do not have any conflicts of interest.