Genotype-by-Environment Interaction and Yield Stability of Maize Single Cross Hybrids Developed from Tropical Inbred Lines

Nitrogen (N) is one of the most important nutrients required for high productivity of the maize plant. In most farmers’ fields in Sub-Saharan Africa (SSA), there is low availability of N in the soil mainly due to continuous cultivation of the land, crop residues removal, little or no application of fertilizers and rapid leaching. There is a need to develop low N tolerant and adapted maize genotypes. Evaluation of maize genotypes under different nitrogen conditions would therefore be useful in identifying genotypes that combine stability with high yield potential for both stress and non-stress environment. Eighty maize hybrids were evaluated at Mbalmayo and Nkolbisson in Cameroon, during 2012 and 2013 minor and major cropping seasons across 11 environments under low and high N conditions. The objectives of the study were: (i) to determine the effect of genotype x environment interaction (G × E) on grain yield and yield stability of single cross maize hybrids across low N and optimum N environments and (ii) to identify genotypes to recommend for further use in the breeding program. Yield data of 80 hybrids were analyzed initially and the analysis of 20 best performing genotypes was further performed for a better visualization and interpretation of the results. Combined analysis of variance showed highly significant G × E effects for grain yield. The GGE biplot analysis divided the study area into three mega environments: one related to the major cropping season while the two others were related to the minor cropping season. The grain yield of the 20 highest yielding hybrids ranged from 4484.7 to 5198.3 kg ha−1. Hybrid 1368 × 87036 was the highest yielding in the minor season while the most outstanding hybrid, TL-11-A-1642-5 × 87036 was the best for the major season. The latter hybrid showed the potential for production across environments and should therefore be further tested in multiple environments to confirm consistency of its high yield performance and stability, and to facilitate its release as a commercial hybrid. High yielding but not stable hybrids across environments could be recommended for the specific environments where they performed well.


Introduction
Maize is one of the most important cereals in Sub-Saharan Africa (SSA) and a staple food for an estimated 50% of the population.It is an important source of carbohydrate, protein, iron, vitamin B, and minerals [1] and accounts for about 15% of the caloric intake of the population [2].In Cameroon, maize is the most consumed cereal, much more than sorghum, rice and wheat [3].Maize is grown in all the five agro ecological zones of the country, namely: Sudano-Sahelian Zone, High Guinea Savanna, Humid Forest Zone and the Western Highlands with a mono-modal rainfall pattern and the Humid Forest Zone with a bimodal rainfall pattern [4].These agro-ecological zones are within an altitude ranging from zero and 4095 m above sea level.The crop is grown both by small and large scale farmers [3,4] under a wide range of conditions such as different soil types, soil fertility levels, moisture levels, different temperatures and cultural practices.
Low soil nitrogen (low N) limits maize yield production in diverse arable land.It is also one of the most important environmental constraints contributing to yield instability of maize.Ajala et al. [5] found large genotypic and phenotypic variances for maize grain yield under low and high N environments with low heritability estimates but with yield gain of more than 25% in all environments.Therefore, selection for grain yield and correlated traits under low N may result in improving maize performance under low N soil environments.Because of the high genotype × environment interactions involved, stressed experiments often produce rankings that differ significantly from one experiment to another, making it difficult to identify the best germplasm [6].In general, maize yields are considerably low under the smallholders farming systems of the tropics than other environments predominantly due to lack of well-adapted and improved cultivars and due to genotype by environment (G × E) interaction [7].Most farmers, especially small scale farmers usually grow varieties based on many criteria, but they usually do not consider the suitability of the variety to the environment which is usually influenced by many biotic and abiotic stresses among which is low soil nitrogen.Consequently, this always results in low yields compared to yields obtained in research stations [3].In the bimodal forest zone of Cameroon, the average maize yield ranged from 0.8 to 1 t ha −1 at farmer level [8].These low yields obtained by farmers are probably due to poor and unstable environmental conditions.Environmental conditions can fluctuate as a result of drought, reduced soil fertility, pressure from insects and diseases [9].It has also been reported [9] that environmental conditions can further be amplified by socio-economic constraints faced by small scale farmers that result in suboptimal input application.The farmers usually have limited access to technology and inputs, especially fertilizer, irrigation facilities and pesticides and have no means to modify or condition the production environment [9].The authors found that including selection under high priority abiotic stresses, such as drought and low N, in a routine breeding program and with adequate weighting can significantly increase maize yields in a highly variable drought-prone environment and particularly at lower yield levels.
Large genotype by environment interactions (G × E) commonly occur under stress conditions; consequently a variety which performs well in one environment during one season or year may not perform well in a different period or in a different site within the same region [10].This is because genotypes exhibit different levels of phenotypic expression under different environmental conditions resulting in crossover performances [11].Genotype × environment interaction is also the result from differences in the sensitivities of genotypes to the conditions in the target environment [11].Genetic × environment interactions (G × E) are of major importance in developing improved genotypes across different environments.When G × E interaction effects are non-significant, means of evaluated varieties across environments are adequate indicators of genotypic performance across the environments.In this situation, the varieties are said to be stable across the environments [12].Significant G × E indicates that selections from one environment may often perform differently in another and the variety is not stable across the environments [12].Therefore, information on G × E may help in determining a breeding strategy.When G × E exists, it is necessary to determine whether there are important crossovers, i.e., rank changes of the genotypes in different environments, such that different winners are picked up in different environments [13].When there is no change in rank of genotypes over environments, there is non-crossover type of interaction effects, and genotypes with superior means can be recommended for all the environments [13].Breeders can also use information on G × E to choose appropriate locations for selection [13].
Maize growers need cultivars that are reliable and consistent across a wide array of stress conditions and have high yield potential that may be expressed when conditions become more favorable [14].Plant breeders should therefore develop cultivars capable of withstanding unpredictable environmental variations [14].In addition, the varieties developed should be stable across environments in order to be widely accepted by farmers throughout a region [14][15][16][17].It is, therefore, important for newly improved maize cultivars to be evaluated at many sites and for a number of years before release [18,19].Unfortunately, in these multi-location trials, varietal selection is often inefficient due to G × E and relative rankings of varieties usually differ across environments [15,17,18,20,21].As a result, it becomes difficult to demonstrate the superiority of any single variety.This can be done through the use of various statistical models [19,20].These statistical analyses give information on adaptability and stability of varieties across target environments.It would then be possible to identify varieties that are appropriate for a specific environment and those with stable performance across environments.
Many stability analysis models exist: joint regression analysis proposed by Eberhart and Russell [22] to estimate the average performance of a genotype in different environments relative to the mean performance of all genotypes in the same environment; multivariate analysis among which are the additive main effects and multiplicative interaction (AMMI) and the genotype and genotype by environment interaction (GGE) biplot analysis [19,20].The AMMI model [23][24][25] and GGE biplot analysis [26] are the most common statistical tools used for the analysis of multi-environment trials (METs) [18,19].The AMMI model combines analysis of variance for the genotype and environment main effects with principal components analysis of the G × E interactions [27].The AMMI method captures a large portion of the G × E interaction sum of squares [28].The AMMI can also help in informing important decisions in breeding programs, such as which genotypes exhibit specific adaptation and in selecting the testing environments [28].This is important for new breeding programs that have not yet optimized their genotype testing networks.The results of an AMMI analysis are often presented in a biplot, which displays both the genotype and environment values and their relationships using the singular vector technique [29].The GGE integrates the genotypic main effect with the G × E interaction effect (Yan et al., 2000).A GGE biplot can help in grouping mega-environments [27].It can also help to identify more representative environments for cultivar evaluation [30] and to compare and rank genotypes using mean yields and stability [31].The objectives of the study were to determine the effect of G × E on grain yield and yield stability of maize single cross hybrids across low N stress and optimal environments.

Germplasm
Eighty single-cross hybrids were used in this study with four of them serving as hybrid checks.The 76 hybrids were selected from 117 F1 developed by crossing thirty nine tropical inbred lines with three testers using a line × tester mating scheme [4].The selection of hybrids was based on seeds availability.The lines included inbred lines from IRAD Cameroon, IITA, CIMMYT and lines from other African maize breeding programs.Of the 39 inbred lines, six were tolerant to low N, four to drought, five to acid soils and four to aluminum toxicity.The testers are parental lines of high yielding hybrids used as checks in this study.The four checks comprised three hybrids (87036 × Exp1 24, 9071 × Exp1 24, 87036 × 9071) from crosses among the 3 testers and hybrid 88069 × Cam inb gp1 17 a promising yellow hybrid of the national breeding program.The hybrid 87036 × Exp1 24 is a high yielding hybrid released in Cameroon and adapted to the Humid Forest Zone of Cameroon.Exp1 24 × 9071 is also a high yielding hybrid, developed from a cross between tropical lowland × temperate converted lines.Genotypes names and codes are presented in Table 1.

Experimental Sites
The 80 hybrids were evaluated at two locations (Mbalmayo and Nkolbisson) of the Humid Forest Zone of Cameroon with a bimodal rainfall pattern.These locations are among the maize growing areas of the Humid Forest Zone of Cameroon where are located the principal experimental sites of the Institute of Agricultural Research of Cameroon.Nkolbisson is located at 11 • 36 E and 3 • 44 N, 5 km from the main capital city 'Yaoundé'.The altitude is 650 m above sea level (a.s.l.).The mean annual rainfall is 1560 mm with bimodal distribution.The average daily temperature is 23.5 • C. The soil is sandy clay [32].Mbalmayo is located at 11 • 30 E and 3 • 31 N.The altitude is 641 m a.s.l.The mean annual rainfall varies from 1017 to 1990 mm with bimodal distribution.The mean monthly temperature varies from 25 • C to 22 • C. The soil is sandy clay [33].Based on the results of soil analysis in 2012 and 2013, the soil in Mbalmayo had a pH of 5.97 which is moderately acidic, while at Nkolbisson, pH was 4.54 and the soil classified as strongly acid.
The main cropping system in Nkolbisson is maize/groundnut/cassava in sole cropping or mixed cropping while in Mbalmayo, other cultivated crops include banana, melon, plantain and vegetables [32].The hybrids were evaluated in a total of 11 environments.Each environment was assigned a code and consisted of a combination of site × year × season × nitrogen level (Table 2).The soil management consisted of two nitrogen levels; Low N (20 kg ha −1 ) and Optimum N (100 kg ha −1 ).The geographical coordinates, climatic conditions of the localities and the 11 environments are described in (Table 2).

Site Preparation and Soil Analysis
The soil was depleted of available Nitrogen by high density maize cropping without fertilizer application, and complete removal of organic matter after harvest [6], in order to establish low N plots in Mbalmayo and Nkolbisson [4].This was done at Mbalmayo thrice between 2010 and 2011 and at Nkolbisson for six growing seasons between 2008 and 2012.To ensure the low N status of the sites [32], composite soil samples were collected before each cropping season and analyzed at the soil laboratory of the International Institute of Tropical Agriculture (IITA) Cameroon [4].Soils were air-dried and ground to pass through a 2 mm sieve.Soil pH in water was determined in a 1:2.5 (w/v) soil: water suspension.Organic C was determined by chromic acid digestion and spectrophotometric analysis [34].Total N was determined from a wet acid digest [35] and analyzed by colorimetric analysis [36].Exchangeable Ca, Mg, K, and Na were extracted using the Mehlich-3 procedure [37] and determined by atomic absorption spectrophotometry.Exchangeable Al extracted using 1N KCl [38] and analyzed using the pyrocatechol violet method described by Mosquera and Mombiela [39].Available P was extracted by Bray-1 procedure and analyzed using the molybdate blue procedure described by Murphy and Riley [40].P expressed in ppm or µg/g; Al, Ca, Mg, K, and Na reported as cmol(+)/kg or me/100 g.Organic C and Total N expressed as % particle size (three fractions) was determined by the hydrometer method.The results of soil analysis are presented in Table 3.

Experimental Design
The experiment was established in two replicates of single row experimental units with an 8 × 10 alpha lattice design using 80 maize hybrids [4].Rows were 5 m long in Mbalmayo and 4 m long at Nkolbisson.Between row spacing was 0.75 m while spacing between hills within the same row was 0.5 m.At planting, each hill received three seeds which were later thinned to two plants for a final density of 53,330 plants per hectare.Split fertilization, weed and pest control were done on each plot as described in Mafouasson et al. [4].

Data Collection
Data for grain yield was obtained as follows: grains were harvested at maturity from each row.The total number of ears and ear weight was recorded for each plot.Five ears were then randomly selected from each plot and their grains were shelled.The "Dickey John" moisture tester was used to measure the percent grain moisture at harvest.Grain yield ((kg ha −1 ) for every entry from the data of fresh ear weight per plot (adjusted to 15% grain moisture) was calculated using the following formula [4]: Fresh ear weight (kg/plot) × (100 − MC) × 0.8 × 10,000 (100 − 15) × Area harvested/plot where: MC = moisture content in grains at harvest (%) 0.8 = Shelling coefficient 10,000 m 2 corresponds to 1 hectare 15% = moisture content required in maize grain at storage

Statistical Analysis
Data obtained was subject to combined analysis of variance (ANOVA) with the PROC GLM procedure in SAS [41] using the RANDOM statement set to the TEST option.Environments were considered as random effects while genotypes were treated as fixed effects.Entry means were adjusted for block effects with reference to lattice design [42].Each environment was defined as year × season × site × nitrogen treatment and the means were separated using Tukey's test at p < 0.05 [4].
The AMMI statistical analysis of yield data was performed with Breeding View in the Integrated Breeding Platform Breeding Management System version 2.1 [43].
GGE biplot analysis was performed using Genstat 15th edition in order to identify genotypes that were suitable for the different environments as well as genotypes stable across the various environments, and to identify the different mega-environments.It was difficult to present the eighty hybrids on the AMMI and GGE biplot.Therefore, for a better visualization and interpretation of AMMI and GGE biplot, the top 20 best performing hybrids across environments and four checks were used for this analysis.

Analysis of Variance for Grain Yield across Environments
The results of the combined ANOVA across environments for the 80 hybrids showed that genotype main effect (G), environment main effect (E) and G × E were all highly significant (p < 0.001) for grain yield (Table 4).The test environments contributed 60.13% of the total variation in the sum of squares for grain yield, while G and G × E sources of variation accounted for 6.81% and 33.05% of the total variation, respectively.The ratio of genotype (G) effect over genotype + genotype × environment (G + G × E) was 0.17.

Yield Performance of the 20 Best Performing Hybrids and Four Checks across Eleven Environments
The 20 best performing hybrids were selected from the 80 hybrids evaluated across environments based on their highest mean yields across the 11 environments.The four checks were added to the 20 hybrids.Yield performance data of these 24 hybrids across eleven environments is presented in Table 5.The overall mean across the 11 environments for the 20 selected hybrids ranged from 4484.7 kg ha −1 to 5198.3 kg ha −1 .The highest yielding hybrid across environments was TL-11-A-1642-5 × 87036 with a yield of 5198 kg ha −1 .All the 20 hybrids selected yielded higher than the four checks.The best check across environments was Exp1 24 × 9071 (3912.4kg ha −1 ) followed by 87036 × Exp1 24 (3908.9kg ha −1 ).The bold and underlined mean yields are for those hybrids that were the highest yielding in each environment.TL-11-A-1642-5 × 87036 was the highest yielding in two optimum environments E4 and E8 with 9531 kg ha −1 and 8874 kg ha −1 .TL-11-A-1642-5 × Exp1 24 was the best performing in E2 (optimum) and E3 (low N) with 6427 and 5402 kg ha −1 respectively.Entrada 29 × Exp1 24 was also the highest yielding in two environments, E9 (low N) and E11 (optimum).One of the hybrid checks (87036 × Exp1 24) was not the best in any environments but was among the five highest yielding hybrids in E6 (optimum), E9 (low N) and E10 (optimum) with grain yield of 4232 kg ha −1 and 6410 kg ha −1 respectively (Table 5).

Additive Main Effect and Multiplicative Interaction (AMMI) Analysis of 24 Hybrids for Grain Yield
The results of AMMI biplot analysis of the 24 hybrids evaluated in 11 environments showed that environment accounted for 59.82% of the total variation in the sum of squares, while genotype and genotype by environment interaction accounted for 7.89% and 32.28 % of variation observed in grain yield respectively (Table 6).In the AMMI biplot (Figure 1) the genotype and environment main effects for grain yield are on the x-axis while the IPCA1 (Interaction Principal Component Axis 1) scores are on the y-axis.The vertical line is the grand mean for grain yield and the horizontal line (y-ordinate) represents the IPCA1 value of zero.2).
This graph identified E1, E3, E5, E7 and E9 as low yielding environments.These were all low N environments in Mbalmayo and Nkolbisson in 2012 and 2013.Environments E4, E6, E8, E11 were identified as high yielding.These were optimum N plots in both locations in 2012 and 2013.The lowest yielding optimum environment was E2.The highest yielding environment was E4 (optimum N, minor season of 2012 at Nkolbisson) while the lowest was E3 (low N, minor season of 2012 at Nkolbisson).

GGE Biplot Analysis of Best 20 Hybrids and Four Checks
The polygon view of the genotypes in the GGE biplot for 24 genotypes is presented in Figure 3.Primary (PC1) and secondary (PC2) scores were significant and explained 29.98% and 21.44% of the variation, respectively.Together they explained 51.42% of the genotype main effect and G × E interaction for the grain yield of maize hybrids evaluated in the 11 environments at Mbalmayo and Nkolbisson in 2012 and 2013.
The polygon view of a GGE biplot displayed the "which-won-where" pattern (Figure 3).The vertices of the polygon were the genotype markers located farthest away from the biplot origin in various directions, such that all genotype markers were contained within the resulting polygon.The biplot was divided into six sectors and three mega-environments and showed five vertex cultivars 1368 × 87036 (G5), TL-11-A-1642-5 × 87036 (G2), Entrada 29 × Exp1 24 (G31), 87036 × 9071 (G79) and 88069 × Cam inb gp1 17 (G80).The first mega-environment comprised E1, E2, E3 and E4 and had 1368 × 87036 as the highest yielding hybrid.These four environments were low N (E1and E3) and optimum N (E2 and E4), minor season of 2012 at Mbalmayo and Nkolbisson.The second mega-environment consisted of E5, E6, E7, E8 and E11 and had TL-11-A-1642-5 × 87036 (G2) as the highest yielding hybrid.These environments were low N (E5 and E7) and optimum N (E8 and E11) of major season in 2013 at Mbalmayo and Nkolbisson plus E11 which is optimum N plot of minor season of 2013 at Nkolbisson.The third comprised E9 and E10 (low N, and optimum N plots of minor season in 2013 at Mbalmayo), with the highest yielding hybrid as 87036 × 9071 (G79).This mega-environment contained 87036 × Exp1 24 (G46).No environment fell within the sector with Entrada 29 × Exp1 24 (G31) and 88069 × Cam inb gp1 17 (G80), indicating that these hybrids were not the best in any of the mega-environments, or they were the poorest cultivars in some or all of the environments.Genotypes within the polygon were less responsive than the vertex genotypes.Ranking of genotypes based on both mean grain yield and stability performance of the 20 best genotypes and four checks is presented in Figure 4 in order to identify the highest yielding and stable genotypes (Figure 4).Genotypes that are located at the center of the concentric circles are the ideal (highest yielding and stable).The GGE biplot identified CLYN246 × 87036 (G1) and TL-11-A-1642-5 × 87036 (G2) as superior since they were located close to the center of the concentric circles.Both were high yielding but TL-11-A-1642-5 × 87036 (G2) was the highest yielding and therefore the most desirable genotype.These hybrids were followed by CLQRCWQ26 × Exp1 24 (G6) and CLWN201 × 87036 (G9) (Figure 4).from the vertical axis at the left and far from the center of the concentric circle, therefore it was the most inferior hybrid in both mean grain yield and stability of performance.

Discussion
The greater variation contributed by environment than those from genotype and genotype × environment interaction indicated that the test environments were highly variable.This result is in agreement with Badu-Apraku et al. [18] who reported that contribution of test environments are much greater than from the other sources of variation in most multi-environmental trials.The highly significant G × E interaction for grain yield justified the use of AMMI and GGE biplots to decompose the G × E interactions and to determine the yield potential and stability of the evaluated single cross hybrids.
The results of the AMMI biplot analysis of the 24 hybrids evaluated in 11 environments also showed that environment effects accounted for 59.82% of the total variation in the sum of squares and was the highest value compared to the other components.The AMMI biplot revealed large variability among the 11 environments, but the yield range among the 24 hybrids was narrow.This is probably because the 20 hybrids were the best selected.ATP S6-20-Y-1 × 87036 and CML395 × Exp1 24 have IPCA1 scores near zero and therefore had small interaction with the environments.This small interaction with environments suggested that these hybrids are stable across environments [13].TL-11-A-1642-5 × 87036 was identified as the highest yielding hybrid.It was followed by CLYN246 × 87036, ATP S6-20-Y-1 × Exp1 24 and CLWN201 × Exp1 24.All these hybrids, except ATP S6-20-Y-1 × Exp1 24 are crosses between CIMMYT and IRAD lines.Acquaah [44] indicated that the development of adapted high yielding hybrids requires that the varieties used as parents are genetically divergent.The high yields obtained between CIMMYT and IRAD lines could therefore imply that they are genetically diverse.The negative interaction of TL-11-A-1642-5 × 87036 with the IPCA1 suggests that this hybrid was less sensitive to environmental changes and was likely to be adapted to unfavorable environments as indicated by Badu-Apraku et al. [18].In contrast, TL-11-A-1642-5 × Exp1 24, 1368 × 87036 and J16-1 × Exp1 24 had large positive interaction with IPCA1 and might be more sensitive to environmental changes, and probably more adapted to favorable environments.
AMMI biplot displayed the distribution of environments from low to high yielding in different quadrants of the graph.This graph placed all low N environments (E1, E3, E5, E7, E9) in the quadrants of lower yielding genotypes and showed the optimum environments (E4, E6, E8, E11) in quadrants of high yielding genotypes as expected.
The GGE biplot analysis of grain yield response and stability of 24 hybrids showed that PC1 explained 29.98% of total variation while PC2 explained 21.44% and together, the two axes accounted for 51.42%.This suggested that the biplot of PC1 and PC2 adequately approximated the environment centered data.The biplot for 24 hybrids was divided into six sectors and three mega-environments in which different cultivars should be selected and deployed to similar environments as suggested by Yan and Tinker [45].According to Yan and Rajcan [46] a mega-environment is defined as the subset of locations that consistently share the best set of genotypes across years and the growing regions are relatively homogeneous with similar biotic and abiotic stresses and cropping system requirements.
Two out of the three mega-environments identified by the GGE biplot included both low and optimum N plots of the two locations, but they were related to different years and different growing seasons.The third mega-environment was related to one specific season of one specific year, but included two nitrogen treatments of one site.This could imply that the mega-environments constructed are based on growing seasons (minor or major) and not on different sites, or different nitrogen treatments.This suggests that seasons and years may have accounted more for significant environmental differences and to different genotypic responses to environments as indicated by Sibiya et al. [10].It might probably be due to similar variation in rainfall amount and distribution as well as biotic stresses within seasons of each year which might have caused the 24 genotypes to have similar relative performance from one environment to another in the mega-environments.In the Bimodal Humid Forest Zone of Cameroon, there are two growing seasons, the major season and the minor season.During the minor season, the total rainfall was lower, the duration of the rainy period is usually shorter than in the major season.Moreover, during the minor season there is prevalence of many diseases such as fungal diseases (e.g., Maize leaf blight caused by Exserohilum turcicum) and maize stem borers among which the main species is Busseola fusca Fuller [47,48].
The results obtained suggest that highest yielding hybrids identified for each mega-environment should be proposed for environments similar to those of these mega-environments.Therefore, hybrids 1368 × 87036 could be proposed for the minor season and TL-11-A-1642-5 × 87036 for the major season.However, this should be done after further evaluation of hybrids in more environments including more locations, years and seasons as recommended by Yan and Tinker [45] who indicated the need for crossover interactions to be repeatable across the years so that target environments can be divided into mega-environments and genotypes be recommended based on METs (multi-environment trials).Yan and Tinker [13] indicated that an ideal genotype should be one that combines both high mean yield performance and high stability across environments; it should be on average environmental coordinate (AEC) on positive direction and have a vector length equal to the longest vector of the genotype as indicated by an arrow pointed to it.Accordingly, the GGE biplot identified TL-11-A-1642-5 × 87036 and CLYN246 × 87036 as closest to the ideal genotype.According to Badu-Apraku et al. [14], in the process of selecting for broad adaptation in maize production, an ideal genotype should have both high mean performance and high stability.TL-11-A-1642-5 × 87036 and CLYN246 × 87036 which were the highest yielding and the most stable hybrids across environments could therefore be selected for broad adaptation (production across environments).These hybrids were followed by CLQRCWQ26 × Exp1 24 and CLWN201 × Exp1 24.The top 20 hybrids performed better than the checks.The poor performance of the check 87036 × Exp1 24, a commercial hybrid, compared to the other hybrids might be due to the fact that it was developed many years ago and might not be adapted to changes (climatic, diseases) that might have occurred in the environments.

Conclusions
This study revealed that genotypes, environments and genotype × environment interaction were significant for grain yield.The genotypes therefore performed differently with respect to yield in each of the eleven test environments and their relative performance varied from one environment to another.AMMI analysis showed that environment effects accounted for a larger proportion of the total variation in the sum of squares for grain yield than genotype effects and genotype × environment effects.The AMMI biplot showed large variability among the environments but a narrow range for yields among hybrids.The GGE biplot classified the study area into three mega-environments.These mega-environments seemed to be related to the two growing seasons of the year (minor and major).High yielding hybrids were identified for each mega-environment and could be proposed for release for production in similar conditions.These hybrids are 1368 × 87036 for mega-environment 1, which is related to the minor season and TL-11-A-1642-5 × 87036 for Mega environment 2, which is related to the major season.The most outstanding hybrid was TL-11-A-1642-5 × 87036.This hybrid has the potential for production across environments and should therefore be tested further in multiple environments to confirm consistency of its high yield performance and stability to facilitate its release as a commercial hybrid.Hybrids which were selected as high yielding, but were not stable across environments could be recommended for the specific environments where they performed well.The results of this study should therefore be confirmed through further evaluation of hybrids at different locations of the Bimodal Humid Forest Zone during both minor and major seasons for several years.

Figure 3 .
Figure 3.A "which won where" biplot based on grain yield of 24 single hybrids evaluated in 11 environments.

Figure 4 .
Figure 4. Comparison view of 24 hybrids with the ideal genotype based on average grain yield and stability for grain yield across 11 environments in 2012 and 2013.

Table 1 .
Names and codes of 80 single cross hybrids evaluated across 11 environments in 2012 and 2013.

Table 2 .
Description of the eleven environments used to evaluate the 80 hybrids.
a.s.l.= above sea level; Low N = low soil nitrogen; Rainfall data were collected at Mbalmayo by IITA and at Nkolbisson by the Rice Project PRODERiP; Major season: From March to June; Minor season: From September to November; * Data for the entire season in this environment were not available.

Table 3 .
Soil characteristics at Mbalmayo and Nkolbisson before the trials in 2012 and nitrogen level in 2013.
nd = Not determined.

Table 4 .
Combined analysis of variance for grain yield of 80 hybrids across eleven environments.

Table 5 .
Mean grain yield (kg ha −1 ) of 20 hybrids and four checks across 11 environments in Mbalmayo and Nkolbisson in 2012 and 2013.
LO = Low N environment; OP = Optimum environment.

Table 6 .
Analysis of variance for additive main effects and multiplicative interaction model for grain yield of 24 hybrids across 11 environments.