Heterosis and combining abilities in a diverse seven-parent pearl millet population diallel tested in West Africa

Pearl millet [Pennisetum glaucum (L.) R. Br.] is an important food-security crop to smallholder farmers in West Africa (WA). Breeding for high yield and stability is a major challenge in the harsh environments of WA but could be tackled by a more systematic exploitation of hybrid vigor and heterosis in breeding of both open-pollinated varieties (OPVs) and different types of hybrids. Knowledge of combining ability patterns and quantitative-genetic parameters is required for an efficient development of hybrid vigor and heterosis in breeding programs. Hence, our objectives were to complement other existing studies and estimate the combining ability of seven unique, highly diverse Sahelian pearl millet populations from Senegal, Mali, Benin, Burkina Faso, Niger, Sudan and Nigeria and the heterosis and stability of their 42 diallel-derived population hybrids (or hybrid populations) to inform pearl millet OPV and hybrid breeding. The materials were evaluated in six environments in WA in 2007. Grain yield (GY) exhibited an average panmictic mid-parent heterosis of 24%, ranging from − 1.51 to 64.69%. General combining ability (GCA) was significant across test environments as reflected by high heritability estimates and high GCA:SCA variance ratios. Thus, early selection for parental per se performance would be rewarding. The parental population from Sudan (IP8679) had strongly negative GCA for GY. Its lack of adaptation contributed to the predominance of additive effects in the present germplasm set. Parental populations PE02987 (Senegal), PE05344 (Mali) and ICMV IS 92222 (Niger) showed large positive GCA for GY. Their offspring, especially PE02987 × PE05344 and Kapelga × ICMV IS 92222, exhibited a high and stable GY across all test environments. Tapping the regional pearl millet genetic diversity and preselecting the crossing parents seem beneficial for OPV and hybrid breeding to increase pearl millet productivity in WA.


Introduction
Pearl millet [Pennisetum glaucum (L.)] is one of the most important cereal crops in arid and semi-arid regions of Africa and India, due to its ability to tolerate severe abiotic stress such as drought, high temperature, poor soil conditions and salinity (Gemenet et al. 2015(Gemenet et al. , 2016Krishnamurthy et al. 2007;Payne et al. 1998;Pucher et al. 2015;van Staveren and Stoop 1985). It is mainly produced by African and Indian smallholder farmers as their major source of food along with its usage as green biomass (Kanfany et al. 2018) and building material (Mason et al. 2015). As it can serve both human and animal nutrition, pearl millet is often used as ''dual purpose'' crop (Kelley et al. 1996) in those countries. In addition, it is used as forage and cover crop in some regions of the United States, South America, and Australia (Bashir 2015). Pearl millet is a diploid cross-pollinated crop, which can profit from the heterosis effect of hybrids (Pandey et al. 2015;Pucher et al. 2016;Sattler et al. 2019;Sattler and Haussmann 2020). Its protogynous type of flowering where the stigmas emerge earlier than the anthers and remain receptive for pollination for 3-4 days (Brunken 1977) increases the rate of outcrossing to [ 70% (Bashir 2015;Burton 1974). The protogynous flowering together with drought tolerance, high seed numbers, good tillering ability, high heterosis and a low seeding rate make it suitable for population hybrid breeding in droughtprone regions (Bashir 2015).
Pearl millet yields are much lower in West Africa (WA) than in India. Grain yields in Niger and Senegal were on average 471 kg ha -1 and 694 kg ha -1 between 2010 and 2014, respectively. In contrast, India had an average pearl millet yield of 1157 kg ha -1 over the same period (FAO 2017). In addition to the very harsh production conditions, a main reason behind this low grain yield in WA is the dependency on landraces for the production of pearl millet (Pucher et al. 2016). Furthermore, the adoption of improved open-pollinated varieties (OPVs) was very limited (Ndjeunga 1997;Smale et al. 2012) due to the lack of access to improved OPV seed along with their performance that did not always fulfill farmers' preferences (Bidinger et al. 2008). The adoption rate of improved varieties in WA countries, where local seed system development was implemented together with farmer-led seed cooperatives varies from 5 to 37% (Christinck et al. 2014). In contrast, this rate wavers around 65-70% in India (Kumara Charyulu et al. 2014).
Using hybrid cultivars combined with improved management practices has brought a tremendous yield increase in pearl millet in India from 305 kg ha -1 in 1951-1955 to 1157 kg ha -1 in 2010-2014(Dave 1986FAO 2017;Yadav and Rai 2013). However, hybrid seed production is more time consuming and cost intensive compared to OPVs and landraces. Therefore, a sufficiently large heterosis effect is required to cover the higher production costs (Duvick and Cassman 1999).
The heterosis effect of population hybrids (based on a cross between two open-pollinated populations) is termed panmictic mid-parent heterosis (PMPH). PMPH is the difference between the population hybrid and the mean of its two open-pollinated parents (Lamkey and Edwards 1999), while panmictic betterparent heterosis (PBPH) refers to the heterosis effect relative to the superior open-pollinated parent. Ouendeba et al. (1993) observed a 55% higher average GY of crosses compared to their parents in a set of 10 WA population hybrids. In comparison, Pucher et al. (2016) found a mean PMPH of 17% in 100 factorialderived WA population hybrids and Sattler et al. (2019) observed a mean PMPH of 23% in 136 diallelderived WA population hybrids. This indicates good scope for pearl millet hybrid breeding in Africa.
Heterotic groups are genetically distinct germplasms groups, which upon crossing show high combining ability. Therefore, keeping those groups separately and not mixing them during the process of hybrid parent development is necessary (Melchinger and Gumber 1998). Significantly superior population hybrids can be developed from crossing genetically distinct WA OPVs and recurrent population improvement methods can be used to stabilize these hybrid populations and develop stable OPVs out of them. Such population hybrids can be a good short-term solution in regions where commercial seed sectors are scarce and underdeveloped, which is the case for many parts of West Africa (Pucher et al. 2016;Sattler et al. 2019;Sattler and Haussmann 2020).
Compared to single-cross hybrids, alternatives like 3-way, 4-way, top-cross or population hybrids could bring rewarding results in these stress-prone environments, due to their higher population buffering capacity (Haussmann et al. 2012). Single-cross hybrids are genetically uniform as they are developed from homozygous inbreed parents and do not express population buffering capacity (Haussmann et al. 2012). This uniformity makes the use of single cross hybrids risky for WA smallholder farmers as they do not have the means to counteract the multitude of biotic and abiotic stresses (Haussmann et al. 2012;van Oosterom et al. 1996;Weltzien et al. 1998;). Furthermore, producing single cross hybrids in WA requires a great effort on resource management, skilled labor and a viable seed industry (Kumara Charyulu et al. 2014). Therefore, single-cross hybrids are considered less appropriate for the extremely variable climate conditions that are prevalent in the Sahel (Sattler and Haussmann 2020).
Genetic diversity analysis in WA pearl millets revealed absence of genetically clearly distinct germplasm groups, rather a high admixture was found (Bashir et al. 2014;Pucher et al. 2015;Sattler et al. 2018;Stich et al. 2010). WA is the center of pearl millets' origin (Burgarella et al. 2018;Manning et al. 2011). In combination with the high outcrossing rate, variable planting dates, prolonged flowering growth stage of OPVs and landraces, and drought-tolerant wind-borne pollen (Haussmann et al. 2007;Hoekstra et al. 1989), this explains the observed high genetic admixture in WA. Consequently, this genetically admixed germplasm complicates the identification of putative heterotic groups in pearl millet using genomic tools (Pucher et al. 2016). However, Diack et al. (2017) were able to find distinct groups of early and late flowering genotypes in Senegalese pearl millet using twelve simple sequence repeat (SSR) markers. Additionally, grouping of distinct Indian pearl millet inbred seed (B) and pollinator (R) lines was possible in more recent studies by using marker-assisted genetic distances, which may shed some light on the presence of putative heterotic gene pools in B and R lines (Ponnaiah et al. 2019;Ramya et al. 2018). Finally, Varshney et al. (2017) studied 345 cultivated and 31 wild accessions using 450,000 single nucleotide polymorphism (SNP) markers. They found an East, a Central and a West Sahel group, and a group of cultivated accessions that was closely related to the wild Central Sahelian accessions. Furthermore, they identified Indian B and R lines, whose offspring showed an 8% above average yield, with a genomic prediction approach. However, the best hybrids were partly derived from intra-group crosses, indicating a large genetic diversity within the B and the R group that could be explored. Pucher et al. 2016 found some high yielding hybrids from the crosses between Senegal 9 Niger or Senegal 9 Nigeria in an experiment of diallel cross with four pearl millet populations from five different WA countries. Sattler et al. 2019 also found good combining ability patterns between the OPVs from Niger and Senegal while crossing 17 pearl millet OPVs in diallel design. These studies yielded interesting results but were not able to propose distinct initial heterotic groups or patterns. However, the populations having higher GCA effects and per se performance could be allocated into eastern and western West and Central Africa to develop heterotic pools for successful population hybrid breeding in pearl millet (Sattler and Haussmann 2020). Therefore, it is crucial to expand existing results with further factorial or diallel mating design-based yield trials to validate combining ability patterns between groups of genotypes across multiple locations and years.
Thus, the present study aims to complement previous results and has the following specific objectives: (1) to determine heterosis effects in pearl millet population hybrids derived from geographically and phenotypically highly diverse Sahelian parental populations; (2) to validate and augment information about combining ability patterns among Sahelian pearl millets; (3) to study G 9 E interaction and stability of pearl millet population hybrids for GY and other agronomic traits; and (4) to complement existing recommendations for pearl millet population hybrid breeding in WA.

Plant materials and experimental design
Seven diverse early to medium maturity pearl millet parental populations (landraces or OPVs) from seven different Sahelian countries were selected from a pearl millet collection (Table 1). Selection was based on the complementarity of different agronomic traits and geographic origin. The control varieties were contributed from the following countries (Institutions): Nigeria (Lake Chad Research Institute, LCRI); Niger/ Mali (International Crops Research Institute for the Semi-Arid Tropics-Sahelian Center, ICRISATSC / Institute of Rural Economy, IER); Niger (National Institute of Agricultural Research, INRAN); Mali (IER); Senegal (Senegalese Institute of Agricultural Research, ISRA). The control varieties used were released from the respective regions and they are the primary OPVs of the respective regions (Table 1).
The populations were crossed in a 7 9 7 diallel mating design by following Griffing's diallel method 3 (F1's and reciprocals) to develop 42 population hybrids (Griffing 1956). The population crosses were created in the offseason 2006/07 at the ICRISAT station in Niger, using three planting dates and at least 25 panicles per female parent were pollinated with a pollen mixture of 25 panicles of the male parent. The spike of the female parent was covered in bags to avoid contamination. Crossing was done by collecting fresh pollen bulk from the male parent and by dusting them on female spike. All 42 population hybrids, their 7 parental populations, and 7 controls were evaluated in 2007 in six environments (locations) of WA ( Table 2). The environments were Bambey (Senegal), Cinzana (Mali), Gampela (Burkina Faso), Maiduguri (Nigeria), Maradi (Niger) and Sadoré (Niger). There was variation among the test environments in terms of rainfall, soil type, temperature and geographical distribution ( Table 2).
At each experimental site, the seeds were sown by following a 7 9 8 alpha lattice design with three replications. Each plot consisted of two rows of 4.8 m length with 0.75 m distance between them resulting in a plot size of 4.5 m 2 . According to data from Pucher et al. (2015) or ICRISAT unpublished data *Data missing In each environment, except Gampela, seven hills were planted per row at an interval of 0.8 m, having 14 hills per plot. In Gampela, each plot had 18 hills, 9 hills per row, as recommended by the Institute of the Environment and Agricultural Research (INERA). At all locations, thinning was done to two plants per hill at around three weeks after planting. Microdoses of 6 g 15 N-15P-15 K or 2 g DAP fertilizers per hill were applied at the time of seed sowing.
Phenotypic data were recorded on 13 agronomic traits: seedling vigor (SV; 5 = excellent; 4 = good; 3 = medium, 2 = poor; 1 = very poor); percentage of downy mildew infested plants (DM%); days to 50% flowering (FLO), measured from the time of sowing to the time when 50% of plant flowered in a plot; plant height (PH; cm); panicle length (PL; cm); panicle circumference (CIR; cm); panicle compactness (CMP; 3 = excellent, 2 = medium and 1 = loose); number of productive panicles per hill (NPP); panicle yield (PY; g m -2 ); grain yield (GY; g m -2 ); dry stover weight (STW; g m -2 ); thousand seed weight (TSW; g); and volumetric weight of grains (VWG; g cm -3 ). PH, CIR, and PL were measured on three representative plants per plot and then the average was calculated. Unfortunately, it was not possible to measure more than three plants per plot because of budget and time constraints. But the fact that we obtained reasonable heritability estimates for these traits (0.86-0.97), illustrate that the method was sufficient for our purpose, and it wasn't necessary to measure more than three plants per plot for these traits. Besides, we did not intend to characterize intra-population variability for these traits, which would have indeed required to measure 10-15 plants per plot. Due to technical problems, VWG in Maiduguri and Maradi, DM% in Maradi and TSW in Maiduguri were not recorded.

Statistical analysis
Raw data of each location were checked for outliers and tested for normality. The statistical analysis was performed following a one-step approach. Adjusted entry means were calculated based on the performance of the parents, controls, and hybrids using the following model across six environments: where y ijkl is the observed phenotype; l represents the general mean; g i refers to the effect of the ith genotype; e j refers to the effect of the jth environment; (g 9 e) ij is the interaction effect between genotype i and environment j; r jk is the effect of the kth replication in the jth environment; b jkl is the effect of lth block nested in the kth replication of jth environment; and e ijkl is the residual effect. For the estimation of entry means, the genotypic effects were considered fixed and all other effects were random. For the variance component estimation, all effects were regarded as random (Piepho et al. 2003). Broad-sense heritability (H 2 ) for all 13 agronomic traits was calculated across environments, using the following formula, where r 2 g denotes genotypic variance component, r 2 gxe corresponds to genotype 9 environment interaction variance component, E is the number of environments, R is the number of replications and r 2 e residual variance component. For single environment repeatability estimation, environment factors and interaction with the environment were discarded from the formula.
where y ijklm is the observed phenotype; l is the general mean; e k refers to the effect of the kth environment; gca i and gca j are the GCA effects of parents i and j, respectively; sca ij is the SCA effect of parent combination i and j; rec ij is the reciprocal effect of parent combination i and j; m i is the maternal effect of female parent i; gca i 9 e k and gca j 9 e k are the GCA interaction effects from parents i and j, respectively with environment k; sca ij 9 e k is the SCA interaction effect of parent combination i and j with environment k; rec ij 9 e k is the reciprocal interaction effect of parent combination i and j with environment k; m i 9 e k is the maternal interaction effect of female parent i with environment k; r kl is the effect of replication l in environment k; b klm is the effect of incomplete block m nested in replication l in environment k; and e ijklm is the residual error variance. For variance component analysis, all the effects were treated as random except for the maternal effect.
PMPH and PBPH were estimated using the following formulas (Pucher et al. 2016): where HP = hybrid performance, MP = mid-parent ((P1 ? P2)/2) and BP = better parent (better performance according to the maximum value, but for FLO and DM%, the minimum value is better). Standard heterosis (SH) was calculated relative to each control variety for GY using the following formula: where HP = hybrid performance, and CV = performance of the control variety. Pearson's correlation coefficients were estimated among the agronomic traits in each environment separately and across environments using the adjusted entry means. Correlations coefficients were also computed between MP and HP; SCA and HP; per se performance of parental population and GCA, and PMPH and SCA. Geographic distances between parental populations' origins were calculated and correlation coefficients between geographic distance and hybrid performance, and geographic distance and PMPH were estimated. For stability analysis of genotypes across environments, the Additive Main Effects and Multiplicative Interaction Model (AMMI) was applied to calculate AMMI Stability Values (ASV; Purchase et al. 2000). Yield stability index (YSI) for each genotype was calculated using the ranking of overall mean performance and the raking of ASV (Farshadfar 2008;Mahmodi et al. 2011). YSI represent both GY performance and stability of the genotypes, with small scores indicating high stability combined with high GY. All calculations were computed within the R environment version 3.5.0 (R Development Core Team, 2018). Correlation coefficient and AMMI stability were analyzed using Agricolae package (de Mendiburu 2014). ASReml-R 3.0 (Butler et al. 2009) package was used for a single environment ANOVA and combined ANOVA.

Mean performance of genotypes at single environments and across environments
Environmental means for GY (across parents, hybrids and controls) ranged from 252 g m -2 at Gampela, making it the best performing environment, to 43 g m -2 at Sadoré, which was the lowest performing environment (Supplemental Table 1). When combined across environments, mean GY ranged from 110 to 223 g m -2 for hybrids, from 96 to 188 g m -2 for parents and from 154 to 217 g m -2 for control varieties (Fig. 1, Supplemental Table 2). The per se performance of parents 1, 2 and 5 was higher than the marginal mean of all parents and they were considered high yielding out of 7 parents. Hybrid 22 (Kapelga 9 ICMV IS 92222) and hybrid 1 (PE02987 9 PE05344) were the two best hybrids showing a mean GY of 223 g m -2 and 220 g m -2 , respectively, which was greater than the best control variety ICMV IS 89305 (217 g m -2 ). These two hybrids were derived from inter-country crosses, namely Burkina Faso (Kapelga) 9 Niger (ICMV IS 92222) and Senegal (PE02987) 9 Mali (PE05344). Wide ranges were observed within each entry group (parents, hybrids, checks) for most other traits, which can be observed from the boxplots below (Fig. 1).

Estimates of heritabilities and variance components
A wide range of heritability was observed for the13 agronomic traits across the six environments. For single locations, the repeatability for GY ranged from 0.30 in Gampela to 0.82 in Bambey and Maradi, while combined across environments, GY showed a heritability of 0.60 (Table 3). The highest heritabilities were observed for PL (0.97) and FLO (0.94) across environments. DM% showed a very low heritability of 0.03. Estimates of r 2 gca were significant for most traits except for SV and DM% (Table 3). Significant r 2 sca were observed for traits CIR, FLO, NPP, PH, TSW, VWG. r 2 rec were only significant for NPP, PL, SV and TSW. The r 2 gca9E interaction variances were significant for most of the traits except for SV and DM%. Significant r 2 sca9E and r 2 rec9E variances were observed for CIR, NPP, PH, PL, and TSW. The GCA: SCA variance ratio of 13 agronomic traits ranged from 0.21 for SV to 72.83 for PL, with all values larger than 1.0 indicating the preponderance of additive gene effects.

Estimates of combining abilities across environments
The GCA effects of the tested parents differed widely (Table 4). PE02987 was generally a good combiner for GY, PY, VWG, PL and CMP as it exhibited a significant positive GCA for these traits. Negative significant GCA effects for FLO exhibited by PE02987, PE03942 and ICMV IS 92222 are effective for breeding early flowering hybrids. Mixed results were observed for the SCA effects of hybrids across all traits (Supplemental Table 4). None of the hybrids was able to show consistent good SCA effects for all the traits. Hybrid 22 (Kapelga 9 ICMV IS 92222) showed the highest and significant positive SCA effect for GY, CMP, PY, and VWG. No significant SCA effects were observed for DM%. REC effects were mostly non-significant for the observed traits.
The correlation coefficients between GY and yield component traits such as CIR, NPP and PL were ranging from -0.01 to 0.49 across environments (data not shown). PH had a significant positive correlation with GY (r = 0.30). The significant negative correlation between DM% and GY (r = -0.29) indicated that the higher is the disease infestation, the higher is the yield loss. TSW was negatively correlated with CMP (-0.32), while CMP was positively correlated with GY (0.36).

Relationship among quantitative-genetic parameters for different traits
Correlations between MP and HP were significant for most of the traits except for SV and DM% (Table 5). The highest correlation was recorded for PL  (r = 0.97). Correspondingly, a highly significant positive correlation was observed between parental performance per se and GCA for most of the traits except for SV, DM% and VWG. No significant association was observed between HP and SCA. Except for DM% (r = -0.31), no significant association was found between PMPH and SCA. Significant negative correlations were observed between geographic distances of parents' origins and HP for the traits PH, PY, GY, and TSW (Supplemental Fig. 1). Geographic distances and PMPH were also significantly negatively correlated with GY (Table 5 and Fig. 2). The nonadapted parent 6 from Sudan largely contributed to the negative relationship (Fig. 2).
G 9 E interaction pattern and yield stability The AMMI biplots (Figs. 3 and 4) for GY illustrate G 9 E interaction patterns and stability of 42 hybrids, 7 parents and 7 control varieties across six environments. In the AMMI1 biplot (Fig. 3), genotypes which are situated on the right-hand side of the grand mean level and close to the line of PC1 = 0, are the high yielding and stable genotypes. Hybrids 22, 28 and 12 show above-average GY and high stability.
In Fig. 4 (AMMI2), PC1 and PC2 explained 52.5% and 18.3% of the total variation, respectively. The vectors of Gampela and Sadoré and, Maiduguri and Cinzana point at a similar direction and therefore, the G 9 E interaction between these environments was low. The vectors of Bambey and Maiduguri point at opposite directions in the AMMI2 biplot, which indicated the presence of high G 9 E interactions between them. Therefore, the ranking of genotypes based on GY among these environments was different. In the AMMI2 biplot, genotypes which were situated close to the center point (PC1, PC2 = 0) are more stable than those, which were distant from the center point. For example, the position of hybrid 22 (H22) fell close to the center point of the AMMI2 plot, therefore the performance of this genotype is highly stable across all six environments. ASV, showing the ranking of stability of the genotypes across environments, ranged from 1.11 to 14.27, where lower ASV indicate more stable genotype and higher ASV indicates less stable genotypes (Supplemental Table 9). YSI ranged from 3 to 110 and hybrid 22 showed the least YSI followed by hybrids 12, 14, 21, 8 and 1. Traits are: see Fig. 1 Significance levels are *significant at p \ 0.05 and **significant at p \ 0.01, respectively Table 5 Mean and ranges of panmictic mid-parent heterosis (PMPH) and panmictic better parent heterosis (PBPH); correlation between geographic distances of parental populations' origins and hybrid performance r(geodist, HP); geographic distances and PMPH r (geodist., PMPH); mid-parent and hybrid performance r (MP,HP); hybrid performance and SCA r (HP, SCA) of 42 population hybrids and parental populations per se performance and GCA r(per se, GCA) for 13 agronomic traits evaluated across six environments in West Africa

Selection of genotypes based on trait-specific performance
In this study, we found significant genetic variation for all agro-morphologic traits, except SV and DM%, among the evaluated WA pearl millet genotypes. This variation could be used effectively for the trait specific selection of genotypes. Sadoré exhibited the lowest average GY among all the test environments. This could be the result of deficient rainfall and high temperature in the growing season because Sadoré received 28% less rainfall than Gampela, the best yielding environment. The low precipitation in September and October in Sadoré led to post-flowering drought stress, which affects the grain filling and results in lower GY (Yadav 2008). The genotypes that performed well in Sadoré can therefore be used to develop drought tolerant hybrids. Despite receiving the highest amount of rainfall, Maiduguri showed a medium average GY compared to the high yielding environments Gampela and Bambey. High precipitation promotes downy mildew infestation, which explains the occurrence of the higher magnitude of average DM % in Maiduguri (Supplemental Table 3). Parents that produced hybrids with both high GY and downy mildew resistance in Maiduguri can be used in downy mildew resistance breeding. However, downy mildew resistance highly depends on G 9 E interactions (Kountche et al. 2013). Therefore, resistant genotypes found in Maiduguri have to be tested in multiple locations to confirm their resistance. Heavy rainfalls, as recorded in Maiduguri, can hamper pollination by washing away the pollen during flowering time resulting in a reduced seed setting and lower GY (Mar et al. 2018;Sattler et al. 2018). Early flowering genotypes have a higher chance to escape terminal drought stress and perform better in environments with marginal rainfall. Early flowering genotypes like PE03942, PE02987 and ICMV IS 92222 could be subjected to further selection for developing early to extra early maturity varieties. The positive relationship of GY with PL indicates that selection of plants with longer PL will increase yield (Bidinger and Raju 2000;Ezeaku et al. 2015;Drabo et al. 2018). Additionally, PL plays an important role when it comes to variety acceptance by Sahelian farmers (Ouendeba 1993). For example, Nigerien farmers prefer long panicles, while Malian Farmers prefer the short type. GY had a significant positive correlation with PH, which is directly related to STW. It is therefore possible to select for GY and fodder yield (Upadhyaya et al. 2018). In general, the selection of traits should be prioritized depending on the preference of farmers and the condition of the agroecological target zones. The extent of heterosis in pearl millet hybrid breeding in West Africa Successful hybrid breeding mainly relies on the magnitude of heterosis for economically important traits in the target region. Presence of heterosis was found in many of the previous studies on pearl millet. The trials that were based on line 9 tester designs with inbred parents cannot be compared directly with our results, but showed the potential of pearl millet top-cross hybrids for WA (Bidinger et al. 1994(Bidinger et al. , 2005Kanfany et al. 2018;Mahalakshmi et al. 1992;Yadav et al. 2000). In our study, the PMPH and PBPH level of GY was relatively higher than of the other traits which agrees with several studies on heterosis in pearl millet (Bidinger et al. 2003;Pucher et al. 2016;Yadav et al. 2006). This study with 24% PMPH and 9% PBPH on average for GY indicates the possibility of increasing total GY productivity and shows good potential for population hybrid breeding in WA (Table 4). A PMPH with 23% on average was observed by Sattler et al. (2019) in 136 population hybrids. Pucher et al. (2016) observed an average PMPH of 17% for GY in 100 WA population hybrids. Similarly, 17% of PMPH was found in 19 hybrids from a cross between 4 Indian landraces and 3 elite populations in African Indian origin by Yadav (2006). Higher prevalence of PBPH is often found in pearl millet. For example, Ouendeba et al (1993) found an average PBPH of 44% in a study on African landraces carried out in two locations in Niger over two years. In our study, the two best yielding population hybrids (1 and 22) exhibited a PBPH of 39.87% and 18.82% respectively, showing the potential to increase yield by breeding population hybrid varieties. However, farmers can also use the control ICMV IS 89305 as a short-term alternative to the best hybrid population (Kapelga 9 ICMV IS 92222) as it yielded only 3% less than the best population hybrid tested.
Comparing the population hybrids to control varieties is necessary to test the yield superiority of the best population hybrids over the best commercially available varieties. The best population hybrid 22 (Kapelga 9 ICMV IS 92222) had a 3-45% higher GY relative to the seven OPV and landrace control varieties. Moreover, this hybrid showed a 17% yield superiority over the improved control OPV Souna 3, which is recommended for production in the National catalogs of Senegal and Niger. The low standard heterosis estimates in some population hybrids were due to inclusion of the non-adapted parent from Sudan in the data set. The good performance of most control varieties underlines their excellent local adaptation.
We observed the highest amount of PMPH for population hybrids in Sadoré (69%), which was significantly positively correlated with the mean GY of population hybrids in the respective environment. The extent of PMPH was also significantly positively correlated with the environmental mean of Cinzana, Gampela and Maradi. These are the environments which received low rainfall during the growing season and might have faced a certain amount of drought stress, which gives a clear indication of the fact that the population hybrids were highly advantageous under drought environments. A similar trend was observed in a study of 36 pearl millet population hybrids in eight environments in India, where hybrids in droughtprone environments showed the highest average PMPH for GY (Presterl and Weltzien 2003). Haussmann et al. (1998) observed similar results from a study with 24 sorghum (Sorghum bicolor (L.) Moench) parental lines and 12 single-cross hybrids in Kenya. However, the ability of population hybrids to perform well under stress-prone conditions depends on the type and selection procedure of parental genetic material used to produce these hybrids (Pucher et al. 2016).
In summary, several studies including ours generally point at the great prospect of pearl millet hybrid breeding in WA and the potential benefits of growing population hybrids for WA smallholder farmers. Such population hybrids can be developed into superior OPVs using recurrent selection methods.
Multilocation testing for the identification of stable hybrids Farmers in WA often face environmental calamities such as drought stress and irregular rainfall. In our study, the precipitation in Bambey, Gampela, and Maiduguri in 2007 differed noticeably from the average annual rainfall. Therefore, it is necessary to develop stable pearl millet cultivars for GY to cope with such environmental variations. We used AMMI biplots to show the GY stability of the population hybrids and their parents in six different environments. Those environments were clearly separated, which reflects the different levels of precipitation and other environmental stress factors. Population hybrids and their parents showed a similar kind of dispersion in the AMMI biplots, which represents a similar GY stability. The high G 9 E interaction for GY observed in our study shows the importance of selecting genotypes based on their GY stability. Some genotypes performed well only in one specific environment like Hybrid 19 (Kapelga 9 PE02987) in Maradi or Hybrid 27 (ICMV IS 92222 9 PE03942) in Bambey. Such site-specific adaptation needs to be validated across years before it can be built on in breeding programs.
As stated before, the risk-averse smallholder farmers in WA mainly prefer stability with high minimum yield (Haussmann et al. 2012). Such static stability is usually due to population buffering, caused by the heterogeneity within population or landrace varieties (Becker and Lèon 1988). Static stability is often associated with low GY (Fasahat et al. 2016), but population hybrid 22 showed both high stability and high GY across all test environments. However, stable hybrids should be chosen based on multiple year trials due to high interannual weather variation. Monitoring their performance consistency is essential before releasing those population hybrids as commercial varieties in the target environments.

Implications of combining abilities for selection of potential parents
Our study showed significantly positive GCA effects for GY for several parental OPVs that can be used in future breeding programs. An effective population hybrid breeding program mainly depends on the accurate identification and selection of parental populations, which can be achieved precisely by using GCA effects of the parents rather than using midparent values (Gowda et al. 2013;Guo et al. 2013). Therefore, parental populations PE02987 and ICMV IS 92222 with significant positive GCA effects are the best candidates for developing high yielding hybrids. GY loss due to drought stress or downy mildew disease is common in many African countries. It is therefore highly beneficial to select populations with drought or downy mildew resistance. Kapelga and PE05344 with a negative GCA for DM% are useful for further breeding programs to reduce GY losses due to downy mildew. However, these OPVs have to be screened under multiple downy mildew affected environments to validate their disease resistance. PH is an important trait to the farmers in WA, especially in Senegal and Niger, as they use the stover for fuel or building materials (Ouendeba et al. 1993). However, the Sahel zone of WA is often characterized by strong winds, which causes lodging especially of taller plants (Izge et al. 2007). Therefore, very elastic stems and possibly a medium plant height are required to reduce the lodging risk. Parental OPVs with a negative GCA for PH could be used to produce a hybrid with medium height and reduced lodging. As a yield component trait, PL is highly important to the WCA farmer. We observed a comparatively large GCA:SCA ratio for PL and STW among the population hybrids in our study. According to several authors (Melchinger, 1999;Melchinger and Gumber 1998;Reif et al. 2007;Schrag et al. 2006), parents belonging to distinct heterotic groups increase the ratio of GCA:SCA variance in their derived crosses. The use of seven highly distinct parental populations (landraces and OPVs) from seven different Sahelian countries in our diallel study resulted in an increased additive variance which translated into a higher GCA:SCA ratio. Pucher et al (2016) also found relatively higher GCA:SCA ratio variance for PL (32.50) in a study of 100 pearl millet population crosses from different WA countries (Senegal, Mali, Burkina Faso, Niger and Nigeria). These results represent the existence of underlying initial heterotic groups in WCA pearl millet populations which needs to be further experimented.
The magnitude of r 2 gca was higher than r 2 sca for GY and other traits indicating the prevalence of additive gene action in the inheritance of these traits, which agrees with other studies on pearl millet (Kanfany et al. 2018;Ouendeba et al. 1993;Zaveri et al. 1988). This ultimately facilitates the prediction of population hybrids' performance, because it can be predicted from their parental per se performance. Large GCA: SCA ratios (Baker 1978) and high H 2 of GY underline the preponderance of additive gene actions. A high correlation between GCA: SCA ratio and H 2 indicating that additive gene actions largely control highly heritable traits was also observed by Pucher et al. (2016) in their study on WA population hybrids. It must be underlined here that the lack of adaptation of the parental population from Sudan (IP8679) largely contributed to the predominance of additive effects in the present germplasm set.
Identification of heterotic patterns and prediction of hybrid performance We observed significant negative associations between PMPH of GY and geographic distances of parental populations' origins. Crosses between two distant parental populations may cause outbreeding depression due to the loss of co-adapted gene complexes and thus restricts the favorable gene flow from one population to another, ultimately leading to a lower GY and heterosis of hybrids (Raabová et al. 2008). In our study, the negative correlation between PMPH and geographic distance primarily reflects the poor performance of the Sudanese parental OPV and its offspring.
Developing heterotic groups based on combining ability patterns are essential for pearl millet population hybrid breeding, as comprehensively described by Sattler and Haussmann (2020). In our study, crosses between parental materials from Senegal and Mali or Burkina Faso and Niger produced some outstanding population hybrids. Potential parental materials from these countries could be separated and crossed systematically via reciprocal recurrent selection to form heterotic groups (Pucher et al. 2016;Sattler and Haussmann 2020). Our study also showed that nonadapted hybrid parents from Sudan should not be considered in a pearl millet hybrid breeding program targeting WA, as their non-adaptation largely reduced hybrid performance in WA.

Conclusions
The present study was meant to complement other existing combining ability studies in Sahelian pearl millets. It generally confirmed the potential benefits of population hybrid breeding and exploitation of heterosis to improve pearl millet production in WA. Crosses of populations from Senegal and Mali or Burkina Faso and Niger produced high yielding hybrids. Incorporating genetic materials from these countries to establish heterotic pools based on combining ability is a promising approach. The best genotype in the present study, hybrid 22 (Kapelga 9 ICMV IS 92222) with 47% and 19% PMPH and PBPH for GY, and outstanding yield stability across all six test environments, seems potentially very suitable for farmers in the WA region. However, data from multiple years are necessary to underpin our findings and to confirm the observed yield stability. Our study also showed that pre-selection of hybrid parents for adaptation to the target environment is key -and that inclusion of too distant or non-adapted parents (for example from Sudan) in the crossing work might not be effective in a breeding program targeting WA.