Register      Login
Crop and Pasture Science Crop and Pasture Science Society
Plant sciences, sustainable farming systems and food quality
RESEARCH ARTICLE (Open Access)

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

Julius Pyton Sserumaga A G , Yoseph Beyene B , Kiru Pillay C , Alois Kullaya D , Sylvester O. Oikeh E , Stephen Mugo B , Lewis Machida B , Ismail Ngolinda F , Godfrey Asea A , Justin Ringo F , Michael Otim A , Grace Abalo A and Barnabas Kiula F
+ Author Affiliations
- Author Affiliations

A National Agricultural Research Organisation, National Crops Resources Research Institute, Namulonge, PO Box 7084 Kampala, Uganda.

B International Maize and Wheat Improvement Center (CIMMYT), ICRAF House, UN Avenue, Gigiri, Village Market, PO Box 1041-00621, Nairobi, Kenya.

C Monsanto, 2 Vermeulen Straat, Petit, 1512, South Africa.

D Mikocheni Agricultural Research Institute, PO Box 6226, Dar es Salaam, Tanzania.

E African Agricultural Technology Foundation (AATF), PO Box 30709-00100, Nairobi, Kenya.

F Ilonga Agricultural Research Institute, PO Box 33, Kilosa, Morogoro, Tanzania.

G Corresponding author. Email: j.serumaga@gmail.com

Crop and Pasture Science 69(7) 691-702 https://doi.org/10.1071/CP17348
Submitted: 20 September 2017  Accepted: 13 June 2018   Published: 4 July 2018

Journal Compilation © CSIRO 2018 Open Access CC BY

Abstract

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.

Additional keywords: correlation, East Africa, G-E interaction, heritability, management.

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 drought-tolerance 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, drought-tolerant 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 non-stress (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.


Materials and methods

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 agro-ecologies 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 α-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):

E1

where Yijk is observed value of the ith progeny of the jth environment in the kth replications; µ 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/Ajk is effect of the kth block within the jth environment; and ε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 (σ2G), location (σ2L), genotype × location (σ2G×L) and error (σ2E) variance were calculated using PROC MIXED (option = REML) of SAS. Broad-sense heritability (H2) 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):

E2

where σ2G is the genotypic variance, σ2E 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:

E3

where σ2G, σ2G×L and σ2E are genotypic, genotype × location and residual variance components, respectively; E is the number of environments; and R is the number of replications. Genotypic correlations (rg) among locations were estimated according to Cooper et al. (1996) as:

E4

where rp(1,2) is the phenotypic correlation between the traits measured in locations 1 and 2; and H1 and H2 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:

E5

where rg is coefficient of correlation between the two locations for a certain trait, σg1 is genetic variance at location 1, σg2 is genetic variance at location 2, σ2g is joint genetic variance of the joint analysis, and σ2ge 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:

E6

where Yij is the average yield of genotype i in environment j; YJ is the average yield across all genotypes in environment j; λ1 and λ2 are the singular values for PC1 and PC2, respectively; ξi1 and ξi2 are the principal component scores PC1 and PC2, respectively, for genotypei; ηj1 and ηj2 are the PC1 and PC2 scores, respectively for environment j; and εij is the residual of the model associated with genotype i in environment j.


Results

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.


Table 1.  Mean-squares and degrees of freedom from ANOVA 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
GEI, Genotype × environment interaction. *P < 0.05; **P < 0.01; ***P < 0.001; n.s., not significant
Click to zoom

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 optimum-moisture 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).


Table 2.  Mean performance of the top 15 hybrids and bottom five hybrids across different drought and optimum environments in Uganda and Tanzania
GY, Grain yield; AD, days to anthesis; ASI, anthesis–silking interval; 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); EH, ear height; EPP, no. of ears per plant; EP, ear position; MOI, grain moisture; PA, plant aspect; PH, plant height. Ranking based on yield
Click to zoom


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. σ2G, σ2G×L, σ2E: Genotypic, genotype × location, and residual variance components, respectively: H2, heritability
Click to zoom

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 early-maturing 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 (H2 = 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 (H2 = 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 H2 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 H2 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 H2 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).


Fig. 1.  Clustering of (a) three locations with significant grain yield under drought conditions, (b) eight locations with significant grain yield under optimum conditions, and (c) 11 locations based on grain yield, in Uganda and Tanzania.
Click to zoom

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 sub-subclusters.

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 sub-subclusters 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).


Fig. 2.  Which-won-where biplot of grain yield of 50 maize varieties evaluated across (a) three locations under drought conditions, and (b) 10 locations under optimum conditions in Uganda and Tanzania. (c) Mean vs stability view of the genotype main effect plus genotype × environment interaction biplot based on yield data of 50 maize varieties evaluated across three locations under drought conditions in Uganda and Tanzania.
Click to zoom

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.


Fig. 3.  Mean vs stability view of the genotype main effect plus genotype × environment interaction biplot based on yield data of 50 maize varieties evaluated across 10 locations under optimum conditions in Uganda and Tanzania.
F3


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). The average grain yield of the top 15 experimental hybrids was higher than the best check under all management conditions, 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 doubled-haploid 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 (rg = 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.



Acknowledgements

This research was supported by the Bill and Melinda Gates and the Howard G. Buffett Foundations, and the United States Agency for International Development through the Water Efficient Maize for Africa project. The constructive comments of anonymous reviewers, which helped to improve the manuscript, are highly appreciated. The authors also thank all field technicians and assistants for data collection at the various experimental sites.


References

Adebayo MA, Menkir A (2014) Assessment of hybrids of drought tolerant maize (Zea mays L.) inbred lines for grain yield and other traits under stress managed conditions. Nigerian Journal of Genetics 28, 19–23.
Assessment of hybrids of drought tolerant maize (Zea mays L.) inbred lines for grain yield and other traits under stress managed conditions.Crossref | GoogleScholarGoogle Scholar |

Alvarado G, López M, Vargas M, Pacheco Á, Rodríguez F, Burgueño J, Crossa J (2015) META-R (Multi Environment Trail Analysis with R for Windows) Version 5.0. International Maize and Wheat Improvement Center (CIMMYT). Available at: http://hdl.handle.net/11529/10201

Badu-Apraku B, Oyekunle M, Obeng-Antwi K, Osuman AS, Ado SG, Coulibay N, Yallou CG, Abdulai M, Boakyewaa GA, Didjeira A (2012) Performance of extra-early maize cultivars based on GGE biplot and AMMI analysis. The Journal of Agricultural Science 150, 473–483.
Performance of extra-early maize cultivars based on GGE biplot and AMMI analysis.Crossref | GoogleScholarGoogle Scholar |

Bänziger M, Betrán FJ, Lafitte HR (1997) Efficiency of high-nitrogen selection environments for improving maize for low-nitrogen target environments. Crop Science 37, 1103–1109.
Efficiency of high-nitrogen selection environments for improving maize for low-nitrogen target environments.Crossref | GoogleScholarGoogle Scholar |

Beyene Y, Mugo S, Pillay K, Tefera TSA, Njoka S, Karaya H, Gakunga J (2011) Testcross performance of doubled haploid maize lines derived from tropical adapted backcross populations. Maydica 56, 351–358.

Beyene Y, Mugo S, Semagn K, Asea G, Trevisan W, Tarekegne A, Tefera T, Gethi J, Kiula B, Gakunga J, Karaya H, Chavangi A (2013) Genetic distance among doubled haploid maize lines and their testcross performance under drought stress and non-stress conditions. Euphytica 192, 379–392.
Genetic distance among doubled haploid maize lines and their testcross performance under drought stress and non-stress conditions.Crossref | GoogleScholarGoogle Scholar |

Beyene Y, Mugo S, Oikeh SO, Juma C, Olsen M, Prasanna BM (2017) Hybrids performance of doubled haploid lines derived from 10 tropical bi-parental maize populations evaluated in contrasting environments in Kenya. African Journal of Biotechnology 16, 371–379.
Hybrids performance of doubled haploid lines derived from 10 tropical bi-parental maize populations evaluated in contrasting environments in Kenya.Crossref | GoogleScholarGoogle Scholar |

Burdon RD (1977) Genetic correlation as a concept for studying genotype-environment interaction in forest tree breeding. Silvae Genetica Journal 26, 168–175.

Cairns JE, Hellin J, Sonder K, Araus LJ, MacRobert JF, Thierfelder C, Prassana BM (2013) Adapting maize production to climate change in sub-Saharan Africa. Food Security 5, 345–360.
Adapting maize production to climate change in sub-Saharan Africa.Crossref | GoogleScholarGoogle Scholar |

Cooper M, Delacy IH, Basford KE (1996) Relationships among analytical methods used to analyze genotypic adaptation in multi-environment trials. In ‘Plant adaptation and crop improvement’. (Eds M Cooper, GL Hammer) pp. 193–224. (CAB International: Wallingford, UK)

Cruz CD, Regazzi AJ (1994) Modelos biométricos aplicados ao melhoramento genético. Universidade Federal de Viçosa, Viçosa, MG, Brazil.

Edge M, Oikeh SO, Kyetere D, Mugo S, Mashingaidze K (2018) Water efficient maize for Africa: A public–private partnership in technology transfer to smallholder farmers in Sub-Saharan Africa. In ‘From agriscience to agribusiness: theories, policies and practices in technology transfer and commercialization’. (Eds N Kalaitzandonakes, EG Carayannis, E Grigoroudis, S Stelios-Rozakis) (Springer: Dordrecht, The Netherlands)

EM-DAT (2016) The International Disaster Database. Centre for Research on the Epidemiology of Disasters. Available at: http://www.emdat.be/database (accessed 25 January 2017).

Ertiro BT, Beyene Y, Das B, Mugo S, Olsen M, Oikeh S, Juma C, Labuschagne M, Prasanna BM (2017) Combining ability and testcross performance of drought-tolerant maize inbred lines under stress and non-stress environments in Kenya. Plant Breeding 136, 197–205.
Combining ability and testcross performance of drought-tolerant maize inbred lines under stress and non-stress environments in Kenya.Crossref | GoogleScholarGoogle Scholar |

FAOSTAT (2014) FAO Statistical Database. Food and Agricultural Organization of the United Nations. Available at: http://faostat3.fao.org/download/Q/QC/E_FAOSTAT (accessed 25 April 2017).

Gauch HG (2006) Statistical analysis of yield trials by AMMI and GGE. Crop Science 46, 1488–1500.
Statistical analysis of yield trials by AMMI and GGE.Crossref | GoogleScholarGoogle Scholar |

Hallauer AR, Carena M, Miranda Filho JB (2010) ‘Quantitative genetics in maize breeding.’ 3rd edn (Iowa State University Press: Ames, IA, USA)

Lorenzana RE, Bernardo R (2008) Genetic correlation between corn performance in organic and conventional production systems. Crop Science 48, 903–910.
Genetic correlation between corn performance in organic and conventional production systems.Crossref | GoogleScholarGoogle Scholar |

Makumbi D, Diallo A, Kanampiu K, Mugo S, Karaya H (2015) Agronomic performance and genotype × environment interaction of herbicide-resistant maize varieties in eastern Africa. Crop Science 55, 540–555.
Agronomic performance and genotype × environment interaction of herbicide-resistant maize varieties in eastern Africa.Crossref | GoogleScholarGoogle Scholar |

Mendes FF, Guimarães LJM, Souza JC, Guimarães PEO, Pacheco CAP, Machado JRDA, Meirelles WF, Silva ARD, Parentoni SN (2012) Adaptability and stability of maize varieties using mixed model methodology. Crop Breeding and Applied Biotechnology 12, 111–117.
Adaptability and stability of maize varieties using mixed model methodology.Crossref | GoogleScholarGoogle Scholar |

Msowoya K, Madani K, Davtalab R, Mirchi A, Lund RJ (2016) Climate change impacts on maize production in the warmer heart of Africa. Water Resources Management 30, 5299–5312.
Climate change impacts on maize production in the warmer heart of Africa.Crossref | GoogleScholarGoogle Scholar |

Oikeh SO, Menkir A, Maziya-Dixon B, Welch RM, Glahn RP, Gauch HG (2004) Environmental stability of iron and zinc concentrations in grain of elite early-maturing tropical maize genotypes grown under field conditions. Journal of Agricultural Science (Cambridge) 142, 1–9.

Oikeh SO, Nganyamo-Majee D, Mugo SIN, Mashingaidze K, Cook V, Stephens M (2014) Water efficient maize for Africa: an example of public-private partnership. In ‘Biotechnology in agriculture and forestry: convergence of food security, energy security, and sustainable agriculture’. Vol. 67, Ch. 13. (Eds DD Songstad, JL Hatfield, DT Tomes) (Springer: Dordrecht, The Netherlands)

Omoyo NN, Wakhungu J, Oteng’i S (2015) Effects of climate variability on maize yield in the arid and semi arid lands of lower eastern Kenya. Agriculture & Food Security 4, 8
Effects of climate variability on maize yield in the arid and semi arid lands of lower eastern Kenya.Crossref | GoogleScholarGoogle Scholar |

Prigge V, Sanchez BSD, Schipprack W, Araus J, Bänziger M, Melchinger AE (2011) Doubled haploids in tropical maize: I. Effects of inducers and source germplasm on in vivo haploid induction rates. Crop Science 51, 1498–1506.
Doubled haploids in tropical maize: I. Effects of inducers and source germplasm on in vivo haploid induction rates.Crossref | GoogleScholarGoogle Scholar |

SAS Institute (2011) SAS/IML 9.3 user’s guide. (SAS Institute Inc.: Cary, NC, USA)

Sserumaga JP, Oikeh SO, Mugo S, Asea G, Otim M, Beyene Y, Abalo G, Kikafunda J (2016) Genotype by environment interactions and agronomic performance of doubled haploids testcross maize (Zea mays L.) hybrids. Euphytica 207, 353–365.
Genotype by environment interactions and agronomic performance of doubled haploids testcross maize (Zea mays L.) hybrids.Crossref | GoogleScholarGoogle Scholar |

Ward JHJ (1963) Hierarchical grouping to optimize an objective function. Journal of the American Statistical Association 58, 236–244.
Hierarchical grouping to optimize an objective function.Crossref | GoogleScholarGoogle Scholar |

Yamada Y (1962) Genotype by environment interaction and genetic correlation of the same trait under different environments. Japanese Journal of Genetics 37, 498–509.
Genotype by environment interaction and genetic correlation of the same trait under different environments.Crossref | GoogleScholarGoogle Scholar |

Yan W (2001) GGE biplot: A windows application for graphical analysis of multi-environment trial data and other types of two-way data. Agronomy Journal 93, 1111–1118.
GGE biplot: A windows application for graphical analysis of multi-environment trial data and other types of two-way data.Crossref | GoogleScholarGoogle Scholar |