- Research article
- Open Access
Yield gains and associated changes in an early yellow bi-parental maize population following genomic selection for Striga resistance and drought tolerance
BMC Plant Biology volume 19, Article number: 129 (2019)
Maize yield potential is rarely maximized in sub-Saharan Africa (SSA) due to the devastating effects of drought stress and Striga hermonthica parasitism. This study was conducted to determine the gains in grain yield and associated changes in an early-maturing yellow bi-parental maize population (TZEI 17 x TZEI 11) F3 following genomic selection (GS) for improved grain yield, Striga resistance and drought tolerance. Fifty S1 lines were extracted from each of cycles C0, C1, C2 and C3 of the population and crossed to a tester TZEI 23 to generate 200 testcrosses. The testcrosses were evaluated under drought, artificial Striga-infested and optimal (free from Striga infestation and without limitation of water and nitrogen) environments in Nigeria, 2014-2017.
Gains in grain yield of 498 kg ha− 1 cycle− 1 (16.9% cycle− 1) and 522 kg ha− 1 cycle− 1 (12.6% cycle− 1) were obtained under Striga-infested and optimal environments, respectively. The yield gain under Striga-infested environments was associated with increased plant and ear heights as well as improvement in root lodging resistance, husk cover, ear aspect and Striga tolerance. Under optimal environments, yield gain was accompanied by increase in plant and ear heights along with improvement of husk cover and ear rot resistance. In contrast, genomic selection did not improve grain yield under drought but resulted in delayed flowering, poor pollen-silk synchrony during flowering and increased ear height. Genetic variances and heritabilities for most measured traits were not significant for the selection cycles under the research environments. Ear aspect was a major contributor to grain yield under all research environments and could serve as an indirect selection criterion for simultaneous improvement of grain yield under drought, Striga and optimal environments.
This study demonstrated that genomic selection was effective for yield improvement in the bi-parental maize population under Striga-infested environments and resulted in concomitant yield gains under optimal environments. However, due to low genetic variability of most traits in the population, progress from further genomic selection could only be guaranteed if new sources of genes for Striga resistance and drought tolerance are introgressed into the population.
Maize (Zea mays L.) is the most important cereal crop playing a crucial food and nutrition roles in sub-Saharan Africa (SSA). It also serves as feed and industrial crop in the sub-region [1, 2]. Maize is widely cultivated in all agro-ecologies, however, the savannas of SSA provide the ideal environment for optimum expression of its yield potential due to high incoming solar radiation, low night temperature and minimized disease pressure. The availability of early and extra-early maize varieties has further facilitated the expansion of maize production into new boundaries particularly the marginal areas where annual rainfall is below 500 mm or where the soils are shallow or sandy . Despite the high prospect for maize production, recurrent drought and Striga hermonthica parasitism are major production constraints of maize in the sub-region [1, 4, 5].
Drought has become a regular occurrence in most agro-ecologies of SSA, due to irregular rainfall patterns and climate change in the sub-region. Yield loss due to drought stress could be as high as 90% when the stress occurs from a few days to anthesis to the beginning of grain filling periods in maize . Striga hermonthica parasitism is another major limitation to maize production, threatening the livelihood of over 300 million people in SSA  and accounting for an estimated loss of staple food crop valued at $7 billion yearly. Yield loss due to infestation by Striga vary from 0 to 100% depending on the severity of the infestation, type of variety under cultivation, climatic conditions and fertility status of the soil . Farmers in the sub-region have experienced complete crop failure under severe Striga infestation and have often been forced to abandon their farmlands. The stresses could occur singly, but invariably, they occur simultaneously under field conditions with devastating consequences on maize production [9, 10]. Control measures for drought include planting in hydromorphic soils, application of irrigation and use of drought tolerant maize varieties while methods employed to mitigate the effect of Striga include hand pulling, crop rotation, use of herbicide, application of fertilizer, fallowing, trap and catch crops, seed treatments and use of Striga resistant maize varieties [11, 12]. However, host plant resistance/tolerance is the most effective, economical and sustainable approach to combat the combined effect of drought and Striga in the sub-region [13, 14].
Development of drought tolerant and/or Striga resistant maize populations, and improvement of such populations through phenotypic recurrent selection has proven to be an effective approach of increasing grain yield, while maintaining genetic variability within the populations [1, 15]. Use of molecular markers could fast track the breeding process when pre-flowering genotypic data are used for selection and recombination. This will ultimately lead to increased genetic gain per unit time and cost. However, the conventional marker-assisted selection method has proven ineffective because only major effect QTLs are used for selection, whereas both major and minor effect genes govern the expression of polygenic traits. Furthermore, the QTL effect estimates are not consistent across populations and environments due to epistasis (gene x gene interactions) and genotype x environment interactions. Genomic selection (GS) is an improved marker-based breeding method that could address the limitations of MAS. In the GS method, all available marker information are incorporated into a predictive model to estimate the genetic values of breeding progenies for selection [16, 17]. As such, each marker is considered a putative quantitative trait locus (QTL) for effective marker effect estimation and minimizes risk of missing small-effect QTLs [16,17,18]. The marker estimates are computed from training population, i.e. breeding material with both phenotypic and genome-wide marker data . The marker effects are subsequently used for computation of genomic estimated breeding values (GEBVs) of new breeding lines in a population earmarked for GS. Heffner et al.  compared prediction accuracy of phenotypic selection (PS), conventional marker-assisted selection (MAS), and GS for 13 agronomic traits in a population of 374 winter wheat (Triticum aestivum L.) advanced-cycle breeding lines and found that the average prediction accuracies using GS were 28% greater than with MAS and were 95% as accurate as PS. However, such information comparing the effectiveness of PS and GS in improving maize under the prevalent contrasting biotic and abiotic stresses in SSA is very scarce.
In an effort to mitigate the combined effect of Striga and drought in the sub-region, the International Institute of Tropical Agriculture (IITA) developed a bi-parental maize population, TZEI 17 x TZEI 11 with combined tolerance to drought and resistance to Striga. The population was subjected to GS which involved one cycle of phenotypic selection and two subsequent cycles of marker-only selection for improved grain yield, drought tolerance and Striga resistance. This study was therefore conducted to: (i) determine the gains in grain yield obtained from GS under Striga infested, drought and optimal environments (ii) identify traits associated with yield gains during GS in the population (iii) estimate genetic variances and heritabilities of traits as GS progressed in the population and (iv) investigate interrelationships among measured traits in the population.
Analyses of variance and genotype x environment interaction
Analysis of variance (ANOVA) across drought environments showed significant environment (env), cycle and entry-within-cycle effects for grain yield and most measured traits except cycle mean squares for stalk lodging and entry-within-cycle mean squares for anthesis-silking interval, root lodging, ear aspect and ears per plant (Table 1). The env x cycle interaction mean squares were significant for all measured traits while no trait showed significant effect for the interaction of entry-within-cycle with env. The env and cycle mean squares showed significant variation for grain yield and most measured traits under Striga-infested environments. The few exceptions included env mean square for anthesis-silking interval and ear rot and cycle mean square for ears per plant (Table 1). The entry-within-cycle effect revealed significant effects for days to anthesis and silking, plant height, stalk lodging and Striga damage at 8 WAP. Significant env x cycle interaction was observed for most measured traits except for ear height, stalk lodging, husk cover, ears per plant, Striga damage (8 WAP) and emerged Striga plants (8 WAP). However, no trait displayed significant effect for the interaction of entry-within-cycle with the environment. Under optimal environments, significant env, cycle and entry-within-cycle mean squares were detected for grain yield and other measured traits except cycle mean squares for root lodging and entry-within-cycle mean squares for root lodging, husk cover, ear rot, and ears per plant (Table 1). The env x cycle interaction effects were significant for most measured traits except for days to anthesis, plant height and root lodging while significant entry-within-cycle x env effects were detected for only days to anthesis and silking, anthesis-silking interval and stalk lodging.
Grain yield ranged from 2252 kg ha− 1 for C3 to 2530 kg ha− 1 for C2 under drought, 3214 kg ha− 1 for C0 to 4905 kg ha− 1 for C3 under Striga infestation and 4385 kg ha− 1 for C0 to 6046 kg ha− 1 for C3 under optimal environments (Fig. 1). The mean grain yield of C3 were higher than those of C0 under Striga and optimal environments, whereas, there was no significant difference between the grain yield of C0 and C3 under drought stress.
Gains in grain yield and changes in other traits
Significant gains in grain yield of 498 and 522 kg ha− 1 cycle− 1 corresponding to 16.9 and 12.6% gain cycle− 1 were detected under Striga-infested and optimal environments, respectively (Table 2). The yield gains under Striga were associated with increase in plant and ear height as well as improvement in root lodging resistance, husk cover, ear aspect and Striga tolerance. Gain in grain yield under optimal environments was accompanied by increased plant and ear heights along with improvement of husk cover and ear rot resistance. In contrast, genomic selection did not improve grain yield under drought stress, but resulted in delayed flowering, poor anthesis-silking synchrony during flowering and increased ear height.
Genetic variances and heritability estimates
Under drought stress, genetic variance estimate was significant for one trait each in cycles C1 and C3 (Table 3). However, four to five traits showed significant genetic variances in cycles C0 and C2. The genetic variance estimates for anthesis-silking interval, stalk and root lodging, husk cover, ear aspect, ear rot and stay green characteristic were not significantly different from zero in all the selection cycles. Significant heritability estimates ranged from 0.42 for plant aspect to 0.60 for plant height in C0 and from 0.42 for plant height to 0.65 for days to silking in cycle C2. Only days to anthesis showed significant heritability estimate of 0.49 in C1 while ears per plant had a significant heritability estimate of 0.47 in C3.
Under artificial Striga infestation, significant genetic variances were detected for grain yield, plant height, ear height and emerged Striga plants (8 WAP) in cycle C0, stalk lodging in C1, plant height and Striga damage (8 WAP) in C2 and days to anthesis and silking, plant height and ear height in cycle C3 (Table 4). Significant heritability estimates ranged from 0.43 for grain yield to 0.68 for plant height in C0 and 0.56 for plant height to 0.70 for days to anthesis in C4. Stalk lodging showed significant heritability estimate of 0.42 in C1 while plant height and Striga damage (8 WAP) had significant heritabilities of 0.54 and 0.58 in C2, respectively.
Under optimal environments, significant genetic variances were observed for grain yield and plant height in C0, plant aspect in C1 and plant and ear heights in C2 (Table 5). Significant heritabilities of 0.44 and 0.60 were obtained for grain yield and plant height in C0. Plant aspect had significant heritability of 0.44 in C2 while significant heritabilities of 0.50 and 0.40 were observed for plant and ear heights in C2. It was striking that no trait had significant genetic variances and heritability estimates in cycle C3.
Interrelationships among traits
Under drought stress, the step-wise regression analyses identified ear aspect, stay green character, ear height, husk cover and stalk lodging as first order traits (primary contributors to grain yield), accounting for about 68% of the variation observed in grain yield (Fig. 2). Ear aspect had the highest direct effect (− 0.559) while stalk lodging had the least (− 0.087). Second order traits included plant aspect, ears per plant, plant height, anthesis-silking interval, days to silking, root lodging and ear rot. Plant height contributed indirectly to grain yield through four first order traits while plant aspect and ears per plant contributed indirectly to grain yield through three of the first order traits. Root lodging was an indirect contributor to grain yield through husk cover and stalk lodging, two first order traits. However, days to silking, anthesis-silking interval and ear rot contributed indirectly to grain yield through one first order trait. Days to anthesis was the only trait categorized as third order trait under drought.
Under Striga infestation, ear aspect, husk cover, ears per plant, stalk lodging, plant height and Striga damage (8 and 10 WAP) were the primary contributors to grain yield. These traits accounted for 82% of the variation in grain yield (Fig. 3). Out of the first order traits, ear aspect had the highest direct effect while husk cover had the least. The second order traits consisted of days to silking, ear height, ear rot, root lodging and emerged Striga plants (8 WAP) while the third order traits comprised days to anthesis, anthesis-silking interval and emerged Striga plants. Days to anthesis, plant height, plant and ear aspect and ears per plant were categorized as first order traits, which accounted for 72% of the variation in grain yield under optimal environments (Fig. 4). Plant height had the greatest contribution to grain yield with a direct effect of 0.447 while ears per plant was the least contributor with a direct effect of 0.105. Days to silking, anthesis-silking interval, ear height, husk cover, stalk lodging and ear rot were identified as second order traits. Root lodging was the only trait classified into the third order group under optimal environments.
Analyses of variance and genotype x environment interaction
The significant environment mean square observed for grain yield and other assayed traits under drought, Striga-infested and optimal conditions in the present study indicated that the environments showed uniqueness in discriminating among the genotypes under each research condition . This could be attributed to the varying environmental factors including amount of incoming solar radiation, temperature, soil type, rainfall pattern and disease pressure at the research environments . The significant cycle effects observed for grain yield and most other agronomic traits measured under drought, Striga-infested and optimal growing environments suggested that there were differential responses in the different cycles which could facilitate selection under the contrasting research conditions. Under drought, the observed significant entry-within-cycle effects for grain yield and other measured agronomic traits except for anthesis-silking interval, root lodging, ears per plant and ear aspect implied that the testcrosses within each cycle varied in the expression of the traits. The testcrosses within each selection cycle also showed variability for the expression of traits including days to anthesis, days to silking, plant height, stalk lodging and Striga damage (8 WAP) under Striga infestation and grain yield, days to anthesis, days to silking, anthesis-silking interval, plant height, ear height, stalk lodging, plant aspect and ear aspect under optimal environments. The significant env x cycle interaction effects detected for all agronomic traits assayed under drought, Striga-infested and optimal environments with the exception of ear height, stalk lodging, husk cover, Striga damage (8 WAP), number of emerged Striga plants (8 WAP) and ears per plant under Striga-infested environments and days to anthesis, root lodging and plant height under optimal environments suggested that the traits in the selection cycles were not expressed in a consistent manner from one environment to the other. In contrast, non-significant entry-within-cycle x env interaction effect observed for all traits assayed under drought and Striga as well as traits other than days to anthesis, days to silking, stalk lodging and anthesis-silking interval under optimal environments indicated that the S1 lines within cycles were consistent in the expression of most agronomic traits measured under each of the research conditions.
Yield gains and changes in other agronomic traits
A primary aim of this study was to determine yield gains following genomic selection, involving one PS cycle followed by two selection cycles based on molecular markers for Striga resistance and drought tolerance. The observed yield gains of 498 and 522 kg ha− 1 cycle− 1 which corresponded to 16.9 and 12.6% gain cycle− 1 under Striga-infested and optimal environments implied that genomic selection was effective for improvement of the population for grain yield under Striga-infested environments and resulted in concomitant increase in grain yield of the population under optimal environments. Since there are presently no documented reports on grain yield gains through genomic selection under Striga-infested conditions, the results of this study were compared with those obtained from S1 family selection under artificial Striga infestation. The gain of 498 kg ha− 1 cycle− 1 obtained for yield under Striga infestation is greater than the 52 kg ha− 1 realized gain documented by Badu-Apraku et al.  following four S1 family selection cycles for improved grain yield as well as resistance to Striga hermonthica in TZE-Y Pop DT STR, an early-maturing maize population. Similarly, the gain of 522 kg ha− 1 cycle− 1 obtained for grain yield under optimal environments was higher than that of Badu-Apraku et al.  who reported yield gains of 130 kg ha− 1 cycle− 1 under optimal conditions after subjecting an early maturing maize population to four S1 family selection cycles for improvement of grain yield along with Striga resistance. The 162 kg ha− 1 cycle− 1 grain yield gains obtained by Edmeades et al.  in two maize populations of early maturity, Pool 16 Sequia (after 1 cycle) and Pool 18 Sequia (after 3 cycles) under well-watered conditions is also less than the yield gain obtained under optimal environments in this study. The observed differences between our results and the findings of other authors could be due to the differences in the germplasm and/or the selection procedures. Our findings clearly demonstrated the superiority of genomic selection over conventional phenotypic selection methods for grain yield improvement under Striga infestation and by extension, under optimal environments. However, the lack of yield gain under drought environments implied that genomic selection was ineffective for yield improvement under this stress factor probably due to low marker density, more minor gene loci, low heritability or reliability of phenotyping. Our findings in the present study is contrary to that of Beyene et al.  who obtained an average of 86 kg ha− 1 yield gain in eight bi-parental maize populations following genomic selection under managed drought environments. A plausible explanation for the differences between our findings and the results obtained by Beyene et al.  could be the dissimilarity of the genetic materials and methodology used or the possibility of trade-offs while improving the bi-parental population concurrently for Striga resistance and tolerance to drought. In addition, Bankole et al.  evaluated S1-testcrosses developed from three selection cycles (C0, C1, C2) of a late-intermediate bi-parental population, following marker-assisted recurrent selection. The authors reported higher grain yield from S1-testcrosses involving the C2 population compared with those generated from the C0 of the population, corresponding to a grain yield gain of 7% per cycle. In contrast, Abdulmalik et al.  detected no significant gain in yield following four selection cycles involving molecular markers in a late-intermediate bi-parental maize population under drought. The authors explained that the findings could have been confounded by the combined effects of drought and army worms during evaluation in the field. These contrasting results indicated that GS just like any other selection method could be efficient for improvement of a genetic material or breeding population depending on the genetic variability, breeding objectives and the methodology employed.
Furthermore, the phenotypic selection index used in the present study integrated increased values of grain yield and higher ears per plant along with decreased values of anthesis-silking interval, plant aspect, ear aspect, Striga damage (8 and 10 WAP), number of emerged Striga plants (8 and 10 WAP) and stay green characteristic, which in combination with the genotypic data of the training population were used to compute the GEBV [26,27,28]. However, the association of yield gains under Striga with increase in plant height and ear height as well as improvement in resistance to root lodging, ear aspect, husk cover and Striga tolerance (i.e. decrease in Striga damage scored at 8 and 10 WAP) implied that genomic selection effectively improved grain yield and resistance to Striga and ear aspect, but could not keep the plant and ear heights constant. This is also true for the yield gains obtained under optimal environments which was accompanied by increased plant and ear heights along with improvement of husk cover and ear rot resistance. Under drought, however, genomic selection was ineffective in improving grain yield and resulted in delayed flowering, poor anthesis-silking synchrony during flowering and increased ear height, which were contrary to the expectations.
Genetic variances and heritability estimates
For a recurrent selection program to be considered as effective and successful it should simultaneously increase the average performance of individuals in the population and at the same time maintain or increase the genetic variability of traits in the population [4, 15, 21, 29]. One of the aims of this study was therefore to examine critically the estimates of genetic variances and heritability of traits during GS under Striga infestation, managed drought and optimal environments. The non-significant genetic variances and heritabilities of most agronomic traits assayed in all the selection cycles suggested low gene frequencies of the traits in the population . In the most advanced selection cycle, C3, the lack of significant genetic variances and heritabilities estimated for grain yield and most measured agronomic traits except ears per plant under drought indicated that the residual variability for most traits would not suffice for progress from further GS in the population. Similar trends were observed in C3 for all measured traits under optimal environments and most measured traits under Striga-infested environments except for days to anthesis, days to silking, plant and ear heights. Our findings are not surprising since the population was improved for Striga resistance and drought tolerance using the GS that originated from only two inbred lines and therefore had a narrow genetic base to begin with. Therefore, there is a need to introgress novel sources of favourable Striga resistance and drought tolerance alleles into the population to ensure further progress from GS for yield improvement, Striga resistance, and drought tolerance .
Interrelationships among traits
Index selection has long replaced selection solely for grain yield under stress conditions. This is because grain yield usually has low heritability under stress factors while some reliable secondary traits maintain high heritabilities under these stress conditions [30, 31]. In the present study, step-wise regression and sequential path analyses were used to investigate the cause and effect relationships involving grain yield and other agronomic traits assayed under the different research conditions. The sequential path analysis is a unique approach that identifies traits making direct contributions to grain yield and categorises them as first order traits, followed by traits with indirect contributions to grain yield through traits of the first order, which are categorized as second order traits and so on [32, 33]. This approach facilitates the classification of secondary traits in declining order of their relative contributions to the observed variation in grain yield. Identification of ear aspect, stay green character, ear height, husk cover and stalk lodging as traits of the first order under drought stress was an indication that these secondary traits were important for yield improvement under the stress factor. Contrary to the results of this study, Talabi et al.  identified anthesis-silking interval, stay green characteristic, ear aspect, plant aspect, and ears per plant as the most important traits accounting for the variation in grain yield. Under Striga-infested environments, husk cover, ear aspect, ears per plant, stalk lodging, plant height and Striga damage (8 and 10 WAP) were categorized as traits of the first order, playing key roles in explaining the variation observed in grain yield. This result agrees partially with the results of Badu-Apraku et al.  where ear aspect and ears per plant were identified as the only direct contributors to grain yield. The identification of days to anthesis, plant height, plant aspect, ear aspect and ears per plant as first order traits under optimal environments implied that the traits served as key players for improvement of grain yield under optimal environments even though GS placed emphasis on improvement in grain yield under Striga infestation and water deficit environments. Badu-Apraku et al.  also identified plant and ear aspects as the first order traits responsible for about 94% of the variation in grain yield under optimal or high-N environments. One important information from the present study and those of previous researchers was the identification of ear aspect as a primary contributor to grain yield under Striga-infested, drought, low-N and optimal environments [32, 33, 35]. This indicated that ear aspect could serve as an indirect selection criterion for simultaneous improvement of grain yield under Striga, low-N, drought and optimal environments. However, scoring for this trait is usually done by trained plant breeders. Promotion and wide adoption of ear aspect by plant breeders would improve the efficiency of breeding maize for tolerance to multiple stresses in the sub-region.
This study has led to the conclusion that genomic selection was effective for grain yield improvement under Striga infestation and resulted in concomitant increase in yield performance under optimal environments. However, under drought environments, genomic selection was ineffective for improvement of grain yield. The yield gain under Striga was associated with increase in plant height and ear height as well as improvement in resistance to root lodging, ear aspect, husk cover, and Striga tolerance while gain in grain yield under optimal environments was accompanied by increase in plant height and ear height along with improvement of husk cover and ear rot resistance. The result of genomic selection was delayed flowering, poor anthesis-silking synchrony during flowering and increased ear height under drought environments. In addition, most traits lacked genetic variability in all the selection cycles particularly the cycle C3, necessitating introgression of novel and beneficial drought tolerance and Striga resistance alleles into the population to ensure progress from further GS. Ear aspect is a key trait that could serve as indirect selection criterion for simultaneous improvement of grain yield under Striga, drought and optimal environments.
Development of the genetic material
In 2007, two yellow maize inbred lines of early maturity, TZEI 17 (Striga resistant) and TZEI 11 (drought tolerant) were selected using the available field data, and inter-crossed to form the bi-parental population, TZEI 17 x TZEI 11. Seeds from the parental inbred lines along with leaf samples harvested from the bi-parental cross at the 3–5 leaf stage were genotyped to confirm the presence of parental type alleles in the bi-parental cross. A total of 108 SNPs were initially generated, after screening of over a thousand KASP assays which were developed by LGC Genomics (United Kingdom) following the conversion of 1536 Illumina Golden Gate Array . However, six SNPs were bad and uncallable; thus, they were not used for the GS. The SNP markers used for GS were: (i) uniform and homozygous in the parental lines, (ii) polymorphic in the parental lines and (iii) heterozygous in the bi-parental population. A total of 15, 20, 17, 5, 14, 12, 7, 7, 2 and 3 SNPs were distributed on chromosomes 1 to 10, respectively (Table 6). Although the 102 SNPs used for GS in the present study fell below the required minimum of approximately 200 SNPs recommended by Zhang et al.  for GS in bi-parental populations, they were the only markers that showed polymorphism for the parent type alleles in the bi-parental population and were therefore used for the GS. The bi-parental population, TZEI 17 x TZEI 11 was subjected to two successive cycles of self-pollination to generate S2 lines (F3 progenies) used as cycle C0 (base population) for GS. A total of 382 S2 lines developed from the bi-parental population TZEI 17 x TZEI 11, were inter-mated with TZEI 23, a standard inbred tester of opposing heterotic group, during the dry season of 2009/2010 in Ibadan. These testcrosses were tested at Mokwa and Abuja under artificial Striga infestation and at Ikenne under optimal growing conditions all in 2010. The testcrosses were equally evaluated under random drought at Bagauda in Nigeria and Chiredzi in Zimbabwe, 2010 and under managed drought at Ikenne during the dry season of 2010/2011. Leaf samples of the S2 lines of the bi-parental population were harvested in the field, for genomic DNA extraction and genotyping using the full complement of the 102 SNPs.
In order to advance the population to Cycle 1 of the selection program, a phenotypic selection index (PSI) comprising grain yield, anthesis-silking interval, plant aspect, ear aspect, ears per plant, Striga damage at 8 and 10 WAP, number of emerged Striga plants at 8 and 10 WAP and stay green characteristic was computed using the results of the multi-location trials. The Best Linear Unbiased Prediction (BLUP)  was used to calculate the marker effect of the 382 F3 lines while a genomic relationship matrix following the methods of Habier et al.  and VanRaden , was used to compute the genomic estimated breeding values (GEBV). Based on the presence or absence of parental type alleles, scores were assigned to each of the F3 lines while the marker score of each line was used to multiply the BLUP values per marker of the F3 lines to obtain the GEBVs. The PSI was used to select the top 10% F3 lines involved in the testcrosses developed from (TZEI 17 x TZEI 11) F3 and TZEI 23 (inbred tester), which were inter-mated to form the (TZEI 17 x TZEI 11) C1 in 2011, using the balanced composite approach. The C1 plants were genotyped with the full complement of SNPs used originally for the genotyping of the C0 population. Based on the SNPs data and PSI, the GEBV was computed for each individual plant in C1. The top 10% of the lines in C1 were selected using the GEBV and recombined following the balanced composite approach to constitute (TZEI 17 x TZEI 11) C2 in 2012. Genotyping was done for C2 plants and the protocol used for recombination in C1 was repeated in C2 to develop (TZEI 17 x TZEI 2) C3 during the dry season of 2012/2013 in the breeding nursery of IITA-Ibadan.
Generation of testcross progenies
Fifty S1 families each were extracted from cycles C0, C1, C2 and C3 of TZE 17 x TZEI 11 F3 population during the 2013 growing season in Ibadan. The resulting 200 S1 lines were crossed to a tester of opposing heterotic group, TZEI 23 to generate testcrosses during the 2013/2014 dry season in the IITA-Ibadan breeding nursery. The 200 testcrosses were tested under managed drought, artificial Striga infestation and optimal growing environments in Nigeria, from 2014 to 2017.
Field evaluation of testcross progenies
We conducted three field experiments between 2014 and 2017. The first experiment involved the testing of the 200 testcrosses under managed drought stress at Ikenne (6°53′N, 3°42′E, 60 m altitude, 1200 mm yearly rainfall), during the dry seasons of 2014/2015 and 2016/2017. These testcrosses were also evaluated under combined drought and heat stress at Kadawa (11°45′N, 8°45′E, 468.5 m altitude, 884 mm yearly rainfall) during the dry season of 2015. The drought trials were established at Ikenne during the dry season in November and supplied with 17 mm of water on weekly basis, with the aid of a sprinkler irrigation system. Managed drought stress was imposed at 28 days after planting, when supply of irrigation water was discontinued such that the plants had to depend on residual moisture in the soil for growth, flowering and grain-filling till physiological maturity. The experimental fields in Ikenne sub-station were fairly uniform and flat, with high water-retention capacity, and the soil type is eutric nitrisol . The combined drought and heat trial was established at Kadawa in February, during the dry season. Water was supplied to the trial twice a week for the first 28 days, using furrow irrigation system. At 28 days after planting, irrigation was discontinued and consequently, the plants were subjected to severe combined heat and drought stress for three consecutive weeks in April, with day temperatures varying from 35 to 39 °C and night temperature ranging between 22 and 27 °C. Subsequently, supply of irrigation water continued once a week during grain filling till the crop attained harvest maturity (with day temperature range of 33 to 40 °C and night temperature range of 24 to 28 °C). The soil type at Kadawa is Cambisol . Basal fertilization of the managed drought experiment at Ikenne and the combined heat and drought trial at Kadawa was carried out using 60 kg each of N, P2O5 and K2O ha− 1 at planting, while topdressing with an additional 60 kg of N ha− 1 was done at 4 WAP.
In the second experiment, the 200 testcrosses were tested under artificial infestation of Striga at Abuja (9о16’N, 7о20’E, altitude 300 m, 1500 mm yearly rainfall) and Mokwa (9°18′N, 5° 4′E, altitude 457 m, 1100 mm yearly rainfall) in the Southern Guinea Savanna of Nigeria in 2014. The fields were injected with ethylene gas at about 10 days before planting, to stimulate suicidal germination of residual Striga seeds in the soil. Artificial infestation of Striga was carried out following the method described by Kim  and Kim & Winslow . Striga seeds sourced from neighbouring fields planted to sorghum were stored for about 6 months to break Striga seed dormancy and subsequently used for the infestation. Each hole in the Striga plot received about 500 germinable seeds of Striga mixed with fine sand in the ratio 1:99 following the procedure described by Kim . Fertilizer application was delayed in the Striga experiment fields until about 25 days after planting (DAP) when 30 kg each of N, P and K ha− 1 was applied as NPK (15–15-15). The delay in fertilization and the reduced application rate were necessary precautions undertaken to enhance germination and emergence of Striga seeds as well as to facilitate the attachment of the Striga plants to the roots of the maize plant . With the exception of Striga plants, all weeds were manually eliminated.
In the third experiment, the 200 testcrosses were evaluated under optimal growing environments i.e. environments free from Striga infestation and without limitations of water and nitrogen at Ikenne, Mokwa and Baguada (12°00′N, 8°22′E, 580 m altitude, 800 mm yearly rainfall) during the 2014 growing season. The trials received 60 kg N ha− 1, 60 kg P2O5 ha− 1 and 60 kg K2O ha− 1 at 2 weeks after planting (WAP) and were top-dressed with 60 kg/ha N at 4 WAP.
A 10 × 5 lattice design with two replications was used for the evaluation of the set of 50 testcrosses derived from the crosses involving 50 S1 lines of each selection cycle and the tester TZEI 23. Randomization was restricted within testcrosses of each selection cycle and the four cycles together constituted each trial in the present study. The experimental units were 3 m long single-row plots, with inter- and within-row spacings of 0.75 and 0.40 m, respectively. Three maize seeds were planted per hill but only two seedlings were retained per hill, following thinning at about 2 weeks after emergence. This gave a final population density of 66,666 plants ha− 1. For the Striga experiments, only pre-emergence herbicide was applied to control the weeds and was complemented with manual weeding. In all other experiments, weeds were controlled using Atrazine and Gramoxone as pre- and post-emergence herbicides, respectively at 5 L/ha each of Primextra and Paraquat.
In the managed drought, combined drought and heat, and optimal environments, data were recorded for days to 50% anthesis and silking (DA and DS), anthesis-silking interval (ASI), plant and ear heights (PLHT and EHT), root and stalk lodging (RL and SL), plant and ear aspects (PASP and EASP) and number of ears per plant (EPP). Details on how the traits were measured have been previously described by Badu-Apraku et al. . Stay-green characteristic (STGR) was scored for the drought/combined drought and heat stress experiments at 70 DAP. For trials conducted under managed drought and combined heat and drought, grain yield (kg ha− 1) was computed from the weight of shelled kernels, adjusted to moisture content of 15%. In contrast, grain yield (kg ha− 1) for the optimal and Striga experiments, were estimated from field weight of ears per plot, assuming a shelling percentage of 80, adjusted to moisture content of 15%. Moisture content at harvest was recorded for representative shelled kernels per plot in all experiments using a moisture meter. The data recorded for Striga trials were essentially the same as those assayed under optimal experiment except that plant aspect was not scored. In addition, Striga damage  were scored at 8 and 10 WAP (SDR1 and SDR2) while the number of emerged Striga plants were also counted at 8 and 10 WAP (ESP1 and ESP2) in the Striga-infested plots. Striga damage was scored per plot on a scale of 1 to 9 where 1 = no damage, an indication of normal plant growth and high resistance, and 9 = total collapse or death of the maize plant, i.e., highly susceptible .
The analysis of variance (ANOVA) was conducted across test environments for each experiment on plot mean basis for grain yield and other agronomic traits with PROC GLM in SAS , using a RANDOM statement with TEST option. Location-year combinations were treated as environments. In the combined ANOVA, genotypes were considered as fixed effects, while the test environments, replications, interaction of genotype × environment (G × E), and all other sources of variation were treated as random effects. The location–year combinations, replicate-within-environment, and block-within-replicate of each experiment were random factors, whereas entries were fixed effects. Mean separation was carried out using the LSD. Score and count data were subjected to natural logarithm transformation prior to ANOVA.
Means for grain yield and other agronomic traits of the testcrosses (dependent variables) were individually regressed on the selection cycles (independent variables) to obtain the coefficient of regression or gain cycle− 1 (b-value) under drought, Striga-infested and optimal environments. The b-value divided by the intercept and then multiplied by 100 provided an estimate of the relative percentage genetic gain per cycle.
As described by Hallauer et al. , estimate of broad-sense heritability was computed on a progeny-mean basis as follows:
where r = number of replications within environment; e = number of environments; σg2 is the variance component due to genotypes; σ2ge = variance component due to genotype × environment interactions; σ2 = the experimental error variance estimate. Standard errors associated with the genetic variances and heritabilities estimated were also computed as described by Hallauer et al. . A variance or heritability estimate was considered significant if it had a value greater than two times the standard error while pair-wise comparison of estimates using corresponding standard errors was used to test for differences among the variances and heritability estimates of testcrosses from the different selection cycles.
The software SPSS version 17.0  was used for step-wise regression analyses, to determine the causal relationships among traits of testcrosses in a sequential order under drought, Striga and optimal environments. Firstly, grain yield (primary trait) was regressed on all other agronomic traits to identify traits with direct significant (P ≤ 0.05) contributions to grain yield which were categorized as traits in the first order. Secondly, each of the traits in the first order category was regressed on other measured agronomic traits not in the first order category, to identify those with indirect contributions to grain yield through the first-order traits, which were categorized as second-order traits. This procedure was repeated as applicable, to identify traits in the third, fourth orders. This approach facilitated the categorization of the predictor traits into first, second and third orders with minimized multicolinearity based on the respective contribution of the traits to the total observable variation in grain yield [33, 34]. The sequential path diagrams following the method of Mohammadi et al.  was used to depict the causal relationship among grain yield and other agronomic traits in a pictorial sequential order. The standardized b values obtained from the stepwise regression analyses provided the estimates of the path co-efficient [32,33,34]. The t-test (P ≤ 0.05) was used to test the statistical significance of the path co-efficients. Thus, only traits showing significant path coefficients were captured in the model along with the percentage variation explained in the dependent variable(s).
Analysis of variance
- C0 :
- C1 :
- C2 :
- C3 :
Days to 50% anthesis
Days after planting
Days to 50% silking
Ear number per plant
Genomic estimated breeding value
International Institute of Tropical Agriculture
Stay green characteristic
Weeks after planting
Badu-Apraku B. Genetic variances and correlations in an early tropical white maize population after three cycles of recurrent selection for Striga resistance. Maydica. 2007;52(2):205–17.
Badu-Apraku B, Akinwale RO, Fakorede MAB. Selection of early maturing maize inbred lines for hybrid production using multiple traits under Striga-infested and Striga-free environments. Maydica. 2010a;55:261–74.
Badu-Apraku B, Oyekunle M, Menkir A, Obeng-Antwi K, Yallou CG, Usman IS, Alidu H. Comparative performance of early-maturing maize cultivars developed in three eras under drought stress and well-watered environments in West Africa. Crop Sci. 2013;53(4):1298–311.
Badu-Apraku B. Effects of recurrent selection for grain yield and resistance in an extra-early maize population. Crop Sci. 2010;50(5):1735–43.
Badu-Apraku B, Menkir A, Ajala S, Akinwale R, Oyekunle M, Obeng-Antwi K. Performance of tropical early-maturing maize cultivars in multiple stress environments. Can J Plant Sci. 2010b;90(6):831–52.
NeSmith DS, Ritchie JT. Effects of soil water-deficits during tassel emergence on development and yield component of maize (Zea mays). Field Crops Res. 1992;28(3):251–6.
Ejeta G. Breeding for striga resistance in sorghum: exploitation of an intricate host parasite biology. Crop Sci. 2007;47:216–27.
Kroschel J. Analysis of the striga problem, the first step towards future joint action, Advances in parasitic weed control at on farm level, vol. 1; 1999. p. 3–26.
Cechin I, Press MC. The influence of nitrogen on growth and photosynthesis of sorghum infected with Striga hermonthica from different provenances. Weed Res. 1993;34(4):289–98.
Kim SK, Adetimirin VO. Striga hermonthica seed inoculum rate effects on maize hybrid tolerance and susceptibility expression. Crop Sci. 1997;37(4):1066–71.
Odhiambo G, Ransom JK. Long term strategies for Striga control. In: Maize research for stress environments: proceedings of the fourth eastern and southern africa regional maize conference, Harare, Zimbabwe, vol. 28; 1994. p. 263–6.
Shaxson L, Riches C. Where once there was grain to burn: a farming system in crisis in eastern Malawi. Outlook Agric. 1998;27(2):101–5.
DeVries J. The inheritance of Striga reactions in maize. Weikersheim, Ibadan: Breeding for Striga resistance in cereals. Margraf Verlag; 2000. p. 73–84.
Badu-Apraku B, Fakorede MAB, Menkir A, Kamara AY, Akanvou L, Chaby Y. Response of early maturing maize to multiple-stresses in the Guinea savanna of west and Central Africa [Zea mays L.]. J Gen and Breed. 2004;58:119–30.
Hallauer AR, Carena MJ, Miranda Filho JD. Quantitative genetics in maize breeding. New York: Springer; 2010.
Bernardo R, Yu J. Prospects for genome wide selection for quantitative traits in maize. Crop Sci. 2007;47(3):1082–90.
Heffner EL, Sorrells ME, Jannink JL. Genomic selection for crop improvement. Crop Sci. 2009;49(1):1–12.
Guo Z, Tucker DM, Lu J, Kishore V, Gay G. Evaluation of genome-wide selection efficiency in maize nested association mapping populations. Theor Appl Genet. 2012;124(2):261–75.
Heffner EL, Jannink JL, Sorrells ME. Genomic selection accuracy using multifamily prediction models in a wheat breeding program. Plant Genome. 2011;4(1):65–75.
Badu-Apraku B, Yallou CG, Haruna A, Talabi AO, Akaogu IC, Annor B, Adeoti A. Genetic improvement of extra-early maize cultivars for grain yield and Striga resistance during three breeding eras. Crop Sci. 2016;56(5):2564–78. https://0-doi-org.brum.beds.ac.uk/10.2135/cropsci2016.02.0089.
Badu-Apraku B, Akinwale RO, Fakorede MAB, Oyekunle M, Franco J. Relative changes in genetic variability and correlations in an early-maturing maize population during recurrent selection. Theor Appl Genet. 2012;125(6):1289–301.
Edmeades GO, Banziger M, Chapman SC, Ribaut JM, Bolanos J. Recent advances in breeding for drought tolerance in maize. In: Badu-Apraku B, Akoroda MO, Ouedraogo M, Quin FM, editors. Contributing to food self-sufficiency: maize research and development in west and Central Africa, Proc. of a Regional Maize workshop. Cotonou: IITA; 1995. p. 24–41.
Beyene Y, Semagn K, Mugo S, Tarekegne A, Babu R, Meisel B, Sehabiague P, Makumbi D, Magorokosho C, Oikeh S, Gakunga J, Vargas M, Olsen M, Prasanna BM, Banziger M, Crossa J. Genetic gains in grain yield through genomic selection in eight bi-parental under drought stress. Crop Sci. 2015;55(1):154–63. https://0-doi-org.brum.beds.ac.uk/10.2135/cropsci2014.07.0460.
Bankole F, Menkir A, Olaoye G, Crossa J, Hearne S, Unachukwu N, Gedil M. Genetic gains in yield and yield related traits under drought stress and favorable environments in a maize population improved using marker assisted recurrent selection. Front Plant Sci. 2017;8:808. https://0-doi-org.brum.beds.ac.uk/10.3389/fpls.2017.00808.
Abdulmalik RO, Menkir A, Meseka SK, Unachukwu N, Ado SG, Olarewaju JD, Aba DA, Hearne S, Crossa J, Gedil M. Genetic gains in grain yield of a maize population improved through marker assisted recurrent selection under stress and non-stress conditions in West Africa. Front Plant Sci. 2017;8:841. https://0-doi-org.brum.beds.ac.uk/10.3389/fpls.2017.00841.
Hayes BJ, Bowman PJ, Chamberlain AJ, Goddard ME. Invited review: Genomic selection in dairy cattle: Progress and challenges. J. Dairy Sci. 2009;92(2):433–43.
Habier D, Fernando RL, Dekkers JCM. The impact of genetic relationship information on genome-assisted breeding values. Genetics. 2007;177(4):2389–97.
VanRaden PM. Efficient methods to compute genomic predictions. J Dairy Sci. 2008;91(11):4414–23.
Badu-Apraku B, Oyekunle M, Talabi AO, Annor B, Akaogu IC. Changes in genetic variances and heritabilities in an early white maize population following S1 selection for grain yield, Striga resistance and drought tolerance. J Agric Sci. 2017a;155(4):629–42.
Bolanȏs J, Edmeades GO. The importance of the anthesis-silking interval in breeding for drought tolerance in tropical maize. Field Crops Res. 1996;48:65–80.
Edmeades GO, Bӓnziger M, Elings A, Chapman SC, Ribaut JM. Recent advances in breeding for drought tolerance in maize. In: Kropff MJ, Van Laar HH, Matthews RB, editors. Applications of systems approaches at the field level. Britain: Kluwer Academic Publishers; 1997. p. 63–78.
Mohammadi SA, Prasanna BM, Singh NN. Sequential path model for determining interrelationships among grain yield and related characters in maize. Crop Sci. 2003;43(5):1690–7.
Talabi AO, Badu-Apraku B, Fakorede MAB. Genetic variances and relationship among traits of an early maturing maize population under drought-stress and low nitrogen environments. Crop Sci. 2017;57(2):681–92.
Badu-Apraku B, Akinwale RO, Oyekunle M. Efficiency of secondary traits in selecting for improved grain yield in extra-early maize under Striga-infested and Striga-free environments. Plant Breed. 2014;133(3):373–80. https://0-doi-org.brum.beds.ac.uk/10.1111/pbr.12163.
Badu-Apraku B, Fakorede MAB, Annor B, Talabi AO. Improvement in grain yield and low-nitrogen tolerance in maize cultivars of three eras. Expl Agric. 2017b:1–19. https://0-doi-org.brum.beds.ac.uk/10.1017/S0014479717000394.
Semagn K, Beyene Y, Warburton ML, Tarekegne A, Mugo S, Meisel B, et al. Meta-analyses of QTL for grain yield and anthesis silking interval in 18 maize populations evaluated under water-stressed and well-watered environments. BMC Genomics. 2013;14(1):1–16.
Zhang X, Pérez-Rodríguez P, Semagn K, Beyene Y, Babu R, López-Cruz MA, et al. Genomic prediction in bi-parental tropical maize populations in water-stressed and well-watered environments using low-density and GBS SNPs. Heredity. 2015;114(3):291–9.
Soil Survey Staff. Soil taxonomy: a basic system of soil classification for making and interpreting soil surveys. In: USDA-NRCS Agriculture Handbook No. 436. 2nd ed. Washington, DC: U.S. Gov. Print. Office; 1999. p. 869.
Kim SK. Breeding maize for Striga tolerance and the development of a field technique. In: Combating Striga in Africa. Proceedings of the international workshop; 1991.
Kim SK, Winslow MD. Progress in breeding maize for Striga tolerance/resistance at IITA. In: Ransom et al, editor. Proceedings of the fifth international symposium on parasitic weeds, Jun. 24–30. Nairobi: CIMMYT; 1991. p. 494–9.
SAS Institute Inc. Statistical Analysis Software (SAS) user’s guide. Cary: SAS Inst; 2011.
SPSS Inc. SPSS base 17.0 for windows user’s guide. Chicago: SPSS Inc; 2007.
We appreciate the IITA Maize Program and Bioscience Center staff for technical assistance during the field and laboratory experiments. We are thankful to Jose Crossa for the statistical model and Sarah Hearne for her support during the development of the population.
This research was conducted through funding support from Bill and Melinda Gates foundation under the Drought Tolerant Maize for Africa Project of IITA. The funding body played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
The datasets used and/or analysed in the present study have been deposited at the IITA repository. DOI: https://0-doi-org.brum.beds.ac.uk/10.25502/szwf-he08.
Ethics approval and consent to participate
Consent for publication
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Badu-Apraku, B., Talabi, A.O., Fakorede, M.A.B. et al. Yield gains and associated changes in an early yellow bi-parental maize population following genomic selection for Striga resistance and drought tolerance. BMC Plant Biol 19, 129 (2019). https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-019-1740-z
- Genomic selection
- Striga resistance
- Drought tolerance