Genetic Gains in Yield and Yield Related Traits under Drought Stress and Favorable Environments in a Maize Population Improved Using Marker Assisted Recurrent Selection

The objective of marker assisted recurrent selection (MARS) is to increase the frequency of favorable marker alleles in a population before inbred line extraction. This approach was used to improve drought tolerance and grain yield (GY) in a biparental cross of two elite drought tolerant lines. The testcrosses of randomly selected 50 S1 lines from each of the three selection cycles (C0, C1, C2) of the MARS population, parental testcrosses and the cross between the two parents (F1) were evaluated under drought stress (DS) and well watered (WW) well as under rainfed conditions to determine genetic gains in GY and other agronomic traits. Also, the S1 lines derived from each selection types were genotyped with single nucleotide polymorphism (SNP) markers. Testcrosses derived from C2 produced significantly higher grain field under DS than those derived from C0 with a relative genetic gain of 7% per cycle. Also, the testcrosses of S1 lines from C2 showed an average genetic gain of 1% per cycle under WW condition and 3% per cycle under rainfed condition. Molecular analysis revealed that the frequency of favorable marker alleles increased from 0.510 at C0 to 0.515 at C2, while the effective number of alleles (Ne) per locus decreased from C0 (1.93) to C2 (1.87). Our results underscore the effectiveness of MARS for improvement of GY under DS condition.


INTRODUCTION
Maize is a staple crop in sub-Saharan (SSA) Africa consumed by over 300 million people. The projected annual maize demand in SSA had been estimated to be 500 million tons by the year 2020, which will surpass the demand for both wheat and rice (IFPRI, 2000). The rapidly growing population in the SSA has therefore necessitated the need to increase maize production in the region therefore; the development of maize varieties with enhanced performance under biotic and abiotic stresses continues to be an important objective (Boomsma and Vyn, 2008). The average yield of maize are significantly lower in Africa compared with the other parts of the world as a result of a host of abiotic and biotic stress factors that affects the crop (FAO, 2010). Under drought stress (DS) and sub-optimal soil nitrogen conditions, the yield of maize is reduced by up to 80% (Bänziger and Lafitte, 1997;Betran et al., 1997;Bänziger et al., 2006). Several studies showed that ear characteristics diminished and the economic yield of maize decreased significantly under drought (Edmeades et al., 1995;Ti-da et al., 2006;Mohammadai et al., 2012). The change in the climate is expected to increase the occasions of drought in Africa (Williams and Funk, 2011), coupled with the expansion of maize production into drought prone regions.
The common conventional breeding approach for drought tolerance involves exploitation of native genetic variation and uses selection to incorporate better characteristics into adapted genotypes (Xoconostle-Cazares et al., 2011). This breeding method has a proven track record of improving tolerance to abiotic stresses (DTMA, 2015). The relative performance of genotypes in DS and favorable environments has been a common criterion to identify desirable genotypes for unpredictable rainfed conditions (Nouri et al., 2011). Direct selection for grain yield (GY) under drought is often more difficult due to low heritability of this trait under stress conditions (Edmeades et al., 1999;Venuprasad et al., 2007;Ziyomo and Bernardo, 2013). In an effort to get the most ideal drought tolerant genotypes for farmers, research is now emphasizing on combining conventional breeding methods with the use of molecular tools (Mhike et al., 2011) to speed up the development and deployment of improved varieties (Beyene et al., 2016).
Marker assisted recurrent selection (MARS) uses markers for identifying several regions on the genome involved in the expression of complex traits and assembling them within a single cross or across related populations (Ribaut et al., 2010). MARS aims at accumulating a large number of quantitative traits loci (QTL) in a population using a subset of markers that are significantly associated with target traits (Bernardo, 2008). Several QTLs associated with GY and other traits under DS and well watered (WW) conditions have been reported (Veldboom and Lee, 1996;Ribaut et al., 1997;Tuberosa et al., 2002;Almeida et al., 2013;Semagn et al., 2013). MARS can then facilitate rapid accumulation of such favorable alleles linked to QTLs in a breeding population (Eathington et al., 2007). When parents are crossed continually based on their genotypes obtained through the use of molecular markers, a genotype that has accumulated favorable alleles at all it loci can be obtained. When parents are crossed based on their molecular marker genotypes, it might be possible to get an ideal genotype after successive generations of crossings (Ragimekula et al., 2013).
The International Institute of Tropical Agriculture (IITA) under the maize improvement programme has improved a biparental maize population through MARS using single nucleotide polymorphism (SNP) markers to enhance GY and other traits under DS condition. The parents used to develop this bi-parental cross were elite drought tolerant maize inbred lines. Understanding genetic gain from MARS and the associated changes in the frequencies of favorable alleles in this population may facilitate the potential use of this approach in developing parental lines of hybrids for drought affected areas. Eathington (2005) and Crosbie et al. (2006) indicated that the genetic gain achieved through MARS in maize was about twice that of phenotypic selection in some reference populations. In upland cotton, Yi et al. (2004) reported significant improvements in resistance to Helicoverpa armigera (cotton bollworm) using MARS. Mayor and Bernardo (2009) compared the usefulness of one and two stage phenotypic selection over MARS in testcrosses of double haploid lines and found that multiple cycles of MARS led to the same level of genetic gain as two-stage phenotypic selection, but at a much lower cost. Beyene et al. (2014) compared the effectiveness of MARS with conventional pedigree selection using hybrids derived from 10 bi-parental maize populations and found higher gain for MARS than that for pedigree selection under both DS and optimum growing conditions. In a separate study, comparing MARS and phenotypic selection, hybrids formed from the top C 1 S 2 lines of each population improved through MARS produced higher mean GYs than hybrids formed from each of the populations improved using phenotypic selection .
Unlike the previous studies that used bulk samples of plants representative of the original and advanced selection cycles of populations improved with MARS for crossing with a common tester (Beyene et al., 2014 evaluating crosses of randomly selected lines from the original and advanced cycles of MARS populations with a common tester provides an opportunity to examine the potential usefulness of improved populations for developing superior drought tolerant maize inbred lines. There is also a need to assess performance of the MARS approach in improving traits in multiple locations under rainfed conditions which are typical of farmers' growing conditions. This information is useful to determine the potential of improved cycles as sources of productive inbred lines for developing hybrids adapted to rainfed conditions. Furthermore, molecular markers such as SNP can facilitate monitoring shifts in allele frequencies in response to selection (Baskaran et al., 2009). Assessment of the molecular changes that occur in MARS populations may therefore provide information on the genetic components underlying the response to selection. Several studies that examined changes in allele frequencies in maize populations undergoing MARS found increases in the frequency of favorable alleles (Bernardo, 2008;Mayor and Bernardo, 2009;Massman et al., 2013). The objectives of the present study were, therefore, to determine genetic gains in GY (i) under DS and WW conditions, (ii) in multiple locations under rainfed conditions and, (iii) assess the associated shifts in the frequency of marker alleles after two cycles of MARS in a bi-parental population.

Development of the MARS Population
A MARS population was developed by crossing two elite drought tolerant maize inbreds (DTPL-W-C7-S2-7-1-1-1-1-B-5-B * 4 and Babangoyo/MO17LPA/Babangoyo-23-4-3-3-B * 6) selected for desirable agronomic traits, resistance to foliar diseases. The resulting F 1 was selfed to generate F 2 bulk seeds, which was grown in 50 rows of 5 m length with a spacing of 0.75 m and self-pollinated to generate 300 F 2:3 lines. A total of 250 F 2:3 lines from this population were planted each in a row and crossed to an inbred tester of the opposite heterotic group. The testcrosses were evaluated under DS and WW conditions at the International Institute of Tropical Agriculture (IITA) Ikenne substation during the dry season. Irrigation was withdrawn from 6 weeks after planting up to the harvest to elicit DS at flowering and grain filling stages.

Genomic Estimated Breeding Value and Selection
Marker effects of the genotyped 250 F 2:3 lines were calculated using the best linear unbiased prediction (BLUP) model (Meuwissen et al., 2001), which permitted predicting genomic estimated breeding value (GEBV) (Henderson, 1975;Gianola and Fernando, 1986). The GEBV was calculated per marker across all the lines derived from the original population (C 0 ) using the BLUP values of GY. Each line was scored 0 or 1 based on the presence or absence of parental allele. The BLUP value per marker was multiplied by the marker score of each line and the resultant values per line were the GEBVs. Significant markers on each chromosome were identified by backward elimination. A relaxed significance level (P = 0.10), which was found to maximize the response to MARS (Hospital et al., 1997;Johnson, 2001), was used. The selected C 0 lines were ranked according to their GEBV values and 10% of the F 2:3 lines with the highest GEBV values were planted ear-to-row and inter-crossed. Bulk pollen collected from ten plants in each line was used for inter-crossing with other lines. Four ears were harvested in each row and shelled to obtain more than 100 seeds per ear. Equal amounts of seed randomly selected from each of the harvested ears were bulked to form the subsequent cycle (C 1 ) which was planted in 20 rows. Leaf samples were collected from each of the plants for genotyping at LGC Genomics for SNP genotyping. Genotyping the C 1 lines were done using SNP markers that were initially used for genotyping the C 0 population. GEBV were estimated for all C 1 individual plants and C 2 was formed by intermating the selected top 10% of the C 1 individuals. All recombination were conducted at Ibadan in Nigeria.

Formation of Testcrosses for Phenotypic Evaluation
In 2013, the original (C 0 ) and advanced selection cycles (C 1 and C 2 ) of the MARS population were planted at Ibadan each in 60 rows of 5 m length spaced 0.75 m apart. Several plants were self-pollinated and 120 to 150 S 1 lines were harvested from each selection cycle and retained. Amongst these, 60 S 1 lines were randomly selected and were planted along with parental lines of the bi-parental cross and an inbred tester (IWD-SYN-STR-C3-5-1-1-B-B) at Ibadan in 2013 to generate testcrosses. The S 1 lines and the testers were used as female and male parents, respectively. The pollen obtained male parent was bulked and used to pollinate the emerged silks of several plants in each S 1 lines. The ears from each testcross were harvested and dried to 12.5% moisture content and shelled.

Genomic Deoxyribonucleic Acid (DNA) Extraction and Isolation
Genomic DNA was isolated from randomly selected leaf samples from each of the 60 S 1 lines planted to generate testcrosses and the two parents 2 weeks after planting. DNA extraction was done using a cetyltrimethylammonium bromide (CTAB) based DNA extraction protocol optimized in our laboratory by Azmach et al. (2013). Samples of 10 to 20 leaf disks representing each of the genotypes were cut and were placed into 1.2 ml strip tubes. A 2.4 mm steel grinding balls were placed into each tube before adding the leaf samples. The tubes were covered and the leaf samples were grinded for 1 min using a Genogrinder 2000 R instrument (BT and C Inc., Basking Ridge, NJ, USA). A 600 µl 2% CTAB extraction buffer was added into each tube and mixed by shaking in the genogrinder at a speed of 1000 strokes for 1 min. The homogenized samples were then incubated for 30 min in a 65 • C water bath to free up the nucleic acids from the leaf tissues and cleanse the DNA from other plant materials such as proteins.
Each tube was centrifuged at 3500 rpm for 15 min after adding chloroform isoamyl alcohol (24:1) into each tube and centrifuged at 3500 rpm for 15 min. Subsequently about 200 µl of the supernatant was transferred to new tubes and 400 µl iso-propanol was added. The two solutions were mixed by inverting the tubes gently for 2 min to allow precipitation of the DNA. The DNA precipitated was then pelleted by centrifuging at 3500 rpm for 15 min and the iso-propanol supernatant was discarded. A 500 µl of 70% ethanol was added to the DNA pellet and the tubes were centrifuged at 3500 rpm for 15 min to clean the DNA pellet. The pellets were dried, dissolved in 100 µl ultra-pure water and stored at 4 • C. DNA quality and concentration were checked by running 2 µl of the diluted DNA sample on 1% agarose gel. DNA quantification was carried out using a NANODROP R spectrophotometer (Thermo Fisher Scientific Inc., Denver, CO, USA). Genomic DNA samples were lyophilized to dry powder and sent to LGC genomics (UK) for SNP genotyping on Kbiosciences' KASP assay platform (KBioscence-LGC Genomics 1 ). The SNP data obtained from this assay was used to assign genotype score to the population.

Evaluation of Testcrosses under Drought Stress and Well Watered Conditions
A trial consisting of 50 testcrosses of randomly selected S 1 lines from each of the three cycles of selection, along with testcrosses of each of the parental lines of the bi-parental cross to the same tester (P 1 × T, P 2 × T), a cross of the two parental lines (F 1 ) and commercial hybrid checks (Oba super 2, Oba super 7, and Oba super 9) were evaluated under DS and WW conditions at Ikenne during the 2014 and 2015 dry seasons. The testcrosses were arranged in a 26 × 6 alpha-lattice design with two replications. Each testcross was planted in single rows of 5 m with 0.75 m space between the rows while the plants were spaced 0.25 m apart within the rows. In the DS trial, irrigation was withdrawn from 5 weeks after planting through harvest to impose DS, whereas the WW trial received irrigation until physiological maturity. A compound fertilizer was applied at the rates of 60 kg N, 60 kg P, and 60 kg K ha −1 at the time of sowing and an additional 60 kg N ha −1 4 weeks later. In each trial, gramoxone and atrazine were applied as pre-emergence herbicides at 5.0 l ha −1 and subsequent manual weeding was carried out to keep the experiments weed-free. Significant at * P ≤ 0.05, * * P ≤ 0.01, and * * * P ≤ 0.001 levels, respectively. Significant at * P ≤ 0.05, * * P ≤ 0.01, and * * * P ≤ 0.001 levels, respectively.

Evaluation of Testcrosses in Multiple Locations under Rain Fed Condition
An experiment consisting of 50 randomly selected testcrosses of S 1 lines from each of the three cycles of selection, along with testcrosses of each of the parental line of the bi-parental cross to the same tester (P 1 × T, P 2 × T), a cross of the two parental lines (F 1 ) and commercial hybrid checks (Oba super 2, Oba super 7, and Oba super 9) were evaluated under rainfed condition at Ikenne, Saminaka, and Bagauda during the 2014 and 2015 growing seasons. The testcrosses were also arranged in a 26 × 6 α-lattice design with two replications. Each testcross was planted in single rows of 5 m with 0.75 m space between the rows while the plants were spaced 0.25 m apart within the rows. Basal application of a compound fertilizer was applied at the rates of 60 kg N, 60 kg P, and 60 kg K ha −1 at the time of sowing and an additional 60 kg N ha −1 4 weeks later. In each trial, gramoxone and atrazine were applied as pre-emergence herbicides at 5.0 l ha −1 and subsequent manual weeding was carried out to keep the experiments weedfree.

Data Collection
Days to anthesis (DP) and silking (DS) were recorded in each of the plot as the number of days from sowing to when half of the plants were shedding pollen grains and showed emerged silks, respectively. Anthesis-silking interval (ASI) was computed as the interval in days between silking and anthesis. Plant height (PH) and ear height (EH) were measured in centimeters as the distance from the base of the plant to the height of the first tassel branch and the node bearing the upper ear, respectively. Plant aspect (PA) was rated on a scale of 1 to 5, where 1 = excellent overall phenotypic appeal and 5 = poor overall phenotypic appeal. Ear aspect (EA) was scored on a scale of 1 to 5, where 1 = clean, uniform, large, and well-filled ears and 5 = rotten, variable, small, Significant at * * P ≤ 0.01 and * * * P ≤ 0.001 levels, respectively.
Frontiers in Plant Science | www.frontiersin.org and partially filled ears. Visual leaf death (LD) was scored only under drought at 12 WAP on a scale of 1 to 9, where 1 = less than 10% senesced leaf and 9 = more than 80% senesced leaf area below the ear. The number of ears per plant (EPP) was the proportion of total number of ears divided by the number of harvested plants. All ears harvested from each plot were shelled to determine percentage moisture and GY adjusted to 15% moisture.

Statistical Analyses
Individual analysis of variance was conducted for each trial using a mixed model in SAS version 9.3 (SAS Institute, 2009).
Replication, year, environment and incomplete blocks were considered as random effects while cycle, testcrosses (cycle) and checks were considered as fixed effects. Combined analyses of variance (ANOVA) was performed for traits recorded under both DS and WW conditions. The average gain per cycle was calculated by regressing means of testcrosses averaged over years on the number of cycles of selection. The percentage response per cycle was obtained by dividing the linear regression coefficient with the intercept multiplied by 100 (Menkir and Kling, 1999). PowerMarker Version 3.25 (Liu and Muse, 2005) and Gen ALEX 6.5 (Smouse and Peakall, 2012) were used for the analysis of molecular data. A total of 180 S 1 lines derived from the three selection cycles and the parental lines were genotyped. Observed allelic and genotypic frequencies for each SNP locus were used to estimate the frequency of the most common alleles, number of alleles, private alleles, observed heterozygosity, and expected heterozygosity.

Estimation of Favorable Marker Alleles
Favorable marker alleles are those alleles at the quantitative trait loci (QTL) that made positive contribution to GY during initial evaluation of the bi-parental cross. At initial (cycle 0) the frequencies of the bi-allelic loci wa equal to 0.5 and the alleles were coded "1" and "2" (Liu et al., 2014). The favorable and unfavorable alleles were determined by the positive and negative sign of the estimated marker effect, respectively. The marker covariates were coded as 0, 1, or 2 for estimating the effects of the markers. The score allotted to each of the markers depends on the favorable allele frequencies that vary in each generation. The rationale behind allocation of scores is that selection should be performed for a sufficient number of times to ensure that all favorable alleles are fixed (Liu et al., 2015). For the markers in each cycle, the mean, maximum, and minimum allele frequencies of the favorable alleles were estimated using SAS version 9.3 (SAS Institute, 2009).

RESULTS
Testcrosses of the S 1 lines derived from C 0 , C 1 , and C 2 of the MARS population were evaluated at Ikenne for 2 years under both DS and WW conditions. The average yield of the evaluated testcrosses under DS represented 22% of their average yield under WW condition. In the combined analysis of variance under DS condition, year did not significantly affect all the measured traits except EH and LD. No significant difference was observed among

Genetic Gain in Grain Yield under Drought Stress and Well-Watered Conditions
The overall mean GY of the MARS population under DS was 1032 kg ha −1 ( Table 3). The testcrosses involving S 1 lines derived from C 2 produced higher mean GY than those derived from C 0 with a relative genetic gain of 7% per cycle. Mean GY of the testcrosses at C 2 was also slightly higher than the parental testcrosses with the same tester and the commercial hybrid checks. The highest mean GY was recorded for the cross between the two parental lines (F 1 ). The testcrosses of C 2 did not differ significantly from other selection cycles for days to flowering, plant and EHs, but differ significantly (ρ > 0.05) for plant and aspect scores. Under WW condition, GY of the testcrosses of S 1 lines derived from C 2 was higher than those derived from C 0 , C 1 , the commercial hybrid checks, one of the parental testcrosses and the cross of the two parents ( Table 4). The average genetic gain in GY for the testcrosses under WW condition was 43 kg per cycle representing a relative gain of 1% per cycle. The gains for other traits were either zero or small.

Genetic Gain in Grain Yield under Rainfed Condition
A separate study was conducted to evaluate these testcrosses at three locations for 2 years under rainfed condition to determine genetic gains in GY and other agronomic traits. The testcrosses involving S 1 lines derived from C 2 produced significantly higher GY than those derived from C 0 with a relative genetic gain of 3% per cycle ( Table 5). Mean GY of the testcrosses at C 2 was also higher than the cross of each parent to the tester, the hybrid between the two parents and the average of commercial hybrid checks. The top five testcrosses selected from each selection cycle under DS and their corresponding mean GYs under WW and rainfed conditions are presented in Table 6. Testcrosses involving S 1 lines derived from C 2 had higher GY than those derived from C 0 and C 1 particularly under DS conditions. The best testcrosses from all selection cycles out-yielded the best commercial hybrid check under DS condition by 6 to 44% and the best parental testcross by 22 to 67%. The best testcrosses selected from each selection cycle under well-watered and rainfed conditions were competitive with the commercial hybrid checks in terms of GYs.

Changes in the Frequencies and Number of Favorable Alleles
Genotyping of randomly selected 180 S 1 lines derived from the three selection cycles (C 0 , C 1 , and C 2 ) and two parental lines of the MARS population using SNP markers showed that the mean frequencies of the favorable marker alleles for GY increased slightly from C 0 to C 2 ( Table 7). The minimum number of favorable marker alleles in C 0 was higher than those of C 1 and C 2 with a difference of 38% at C 1 and 30% at C 2 . The maximum number of combinations of favorable alleles present in the S 1 lines derived from C 1 and C 2 was higher than those derived from C 0. The changes in the maximum number of combinations of favorable alleles found in S 1 lines derived from the original and advanced cycles did not follow a consistent trend. On the other hand, the average number of combination of alleles present in the best 10 lines increased slightly from C 0 to C 2 ( Table 8).
The changes in the effective number of alleles, level of heterozygosity, and inbreeding coefficient are presented in Table 9. The number and the frequency of different alleles were similar across cycles but the number of effective alleles was higher at C 0 than at C 2 . The level of heterozygosity decreased from 0.591 at C 0 to 0.500 at C 2 whereas the average expected heterozygosity (H e ) slightly increased from 0.458 at C 0 to 0.478 at C 2 . H o was higher than the H e for all the cycles. As shown in Table 10, five genotypes were lost at C 1 and nine were lost at C 2 .

DISCUSSION
In the present study, the intensity and timing of DS imposed was severe to elicit differential responses of the testcrosses of S 1 lines. The observed yield under DS relative to WW condition fell within the range of 20 to 30% which is considered to be severe DS (Banziger et al., 2000). The significant difference in GY among the testcrosses (cycles) under DS condition indicated that MARS was effective in shifting the frequencies of favorable alleles associated with this trait. The observed significant genetic variation among testcrosses may thus provide us with an opportunity to select maize lines with high yield and acceptable agronomic traits. The absence of interaction of year with testcrosses (cycle) for GY and other traits under DS indicated a consistent performance of the testcrosses across years.
The observed higher mean GY of the testcrosses at C 2 over that at C 0 , parental testcrosses and commercial checks under DS and well-watered condition suggest that MARS was effective for the improvement of complex traits such as tolerance to DS consistent with the findings of Beyene et al. (2015Beyene et al. ( , 2016. Testcrosses of C 2 did not differ significantly from C 0 in PH and days to anthesis under DS and WW conditions, indicating that the increase in GY at C 2 was not associated with an increase in flowering and PH. Previous studies across diverse tropical maize populations Semagn et al., 2015) have also reported a similar trend in which the superiority of MARS was demonstrated without significantly affecting PH and maturity of most populations. The higher GY recorded for the parental testcrosses (F 1 ) under DS condition for 2 years may be related to the superiority of the highly homogeneous F 1 relative to the heterogeneous testcrosses involving S 1 lines derived from C 0 , C 1 , and C 2 .
The changes in the frequencies and combinations of favorable alleles with selection though not significant, showed that MARS experienced a slight but consistent increase in accumulating favorable alleles for GY. The gradual increases in mean frequencies and number of favorable alleles for the MARS population are in line with the findings of Mayor and Bernardo (2009) who reported that MARS led to desirable changes in the frequencies of favorable marker alleles. Edwards and Johnson (1994) had earlier reported increase in the frequencies of favorable alleles from 0.5 to 0.8 in a sweet corn F 2 population. Eathington et al. (2007) also reported changes in the frequencies of favorable alleles from 0.56 to 0.96 in 20 different QTL regions, which significantly increased the possibility of obtaining a desirable genotype.
The maximum number of effective alleles (N e ) is 2 with bi-allelic SNP markers. The fact that the average N e decreased from 1.934 in the starting population to 1.869 in C 2 suggests that this population is gradually moving toward allele fixation (Vogel, 2010). The reduction in the level of heterozygosity and an increase in the values of the fixation index from C 0 to C 2 imply loss of heterozygosity possibly because of positive assortative mating and selection of homozygotes (Smouse and Peakall, 2012). The observed decrease in the number of heterozygotes and the increase in the level of inbreeding resulting from response to selection were in line with the findings of Massman et al. (2013). The loss of genotypes suggests that there were genetic changes in the population due to inadvertent selection against those homozygote alleles. Nonetheless, the loss of these genotypes had no adverse effect on improvements in GY and other traits.
In summary, this study examined the genetic gain in a bi-parental cross of elite drought tolerant lines with good agronomic traits, resistance to maize streak virus disease and other foliar diseases that had undergone two cycles of MARS. The population did respond to MARS for improvement in GY and other agronomic traits. These findings highlight the effectiveness of MARS in improving bi-parental populations as potential success of improved drought tolerant inbred lines with superior agronomic traits under drought condition.

AUTHOR CONTRIBUTIONS
FB: Acquisition of data, analysis and interpretation of both phenotypic and molecular aspects of the work, and drafted the work. AM: Designed the work, revised it critically for intellectual content and approval of final version to be published. MG: Designed and supervised the molecular aspect of the work, revised the aspect of the work critically for intellectual content and also approved the final version to be published. GO: Supervision of the field and quantitative genetics aspect of the work, contributed to revising it critically for intellectual content and approval of final version to be published. JC: Analysis of the molecular data during the development of the population, designed the form of selection to be used, revised and approved the final version to be published. SH: Substantial contribution to the design, analysis and interpretation of molecular data while developing the population. NU: Acquisition, analysis, and interpretation of molecular data.

FUNDING
This project was part of the drought tolerant maize for Africa Project funded by Bill and Melinda Gates foundation.