- 1Mision Biologica de Galicia (CSIC), Pontevedra, Spain
- 2Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias (INIFAP), Chihuahua, Mexico
- 3Department of Crop Science, North Carolina State University, Raleigh, NC, United States
Introduction: The study of yield and resistance/tolerance to pest are related traits fundamental for maize breeding programs. Genomic selection (GS), which uses all marker information to calculate genomic breeding values, is presented as an emerging alternative to phenotypic and marker-assisted selections for improving complex traits controlled by many genes with small effects. Therefore, although phenotypic selection (PS) has been effective for increasing resistance and yield under high infestation with maize stem borers, higher genetic gains are expected to be obtained through GS based on the complex architecture of both traits. Our objective was to test whether GS is more effective than PS for improving resistance and/or tolerance to maize stem borers and grain yield.
Methods: For this, we compared different selection programs based on phenotype and genotypic value for a single trait, resistance or yield, and for both traits together.
Results and discussion: We obtained that GS achieved the highest genetic gain for yield, meanwhile phenotypic selection for yield was the program that achieved the highest reduction of tunnel length, but was ineffective for increasing yield. However, phenotypic or genomic selection for increased resistance may be more effective in improving both traits together; although the gains per cycle would be small for both traits.
1 Introduction
The Mediterranean Corn Borer (MCB), Sesamia nonagrioides Lefèbvre (Noctuidae), is a major pest of maize (Zea mays L.) in the Mediterranean region (Velasco et al., 2007; Jika et al., 2020). In Spain, MCB second and subsequent generations larvae cause the main damages inside maize stems by excavating longitudinal tunnels that destroy the pith and negatively affect plant development. Tunnels increase the susceptibility to plant lodging and interferes with assimilate movement to developing cobs, which reduces maize yields (Malvar et al., 2008). This has led to the search for varieties with increased natural resistance to MCB damage, which can be the basis for breeding programs aimed at releasing resistant and/or tolerant varieties (Butrón et al., 1999; Butrón et al., 2005; Revilla et al., 2007). Most studies have evaluated insect resistance in a wide sense, using the pith tunnel length produced by MCB under artificial infestation as the classifying criteria (Malvar et al., 1993; Butrón et al., 1999; Velasco et al., 1999a; Velasco et al., 1999b; Malvar et al., 2004; Velasco et al., 2004a; Butrón et al., 2005; Butrón et al., 2006; Malvar et al., 2007; Revilla et al., 2007). Once sources of resistance were identified (Butrón et al., 1999; Malvar et al., 2004), research was oriented to study trait inheritance (Cartea et al., 1999; Velasco et al., 2002; Velasco et al., 2004b; Butrón et al., 2009) and to develop the most suitable phenotypic selection programs to increase maize resistance and/or tolerance to MCB (Sandoya et al., 2008; Butrón et al., 2012; Ordás et al., 2012; Samayoa et al., 2012; Butrón et al., 2014).
Nevertheless, valuable genetic gains through phenotypic selection have been hindered by the unfavorable genetic correlation between tunnel length by MCB and yield (Sandoya et al., 2010; Butrón et al., 2012). Population yield was maintained in the first three cycles of phenotypic selection (PS) for reducing stem tunneling by MCB maintaining yield (Sandoya et al., 2008), but the yield decreased in the subsequent cycles (Butrón et al., 2012). In addition, the search for QTL for resistance to MCB damage and yield in segregating populations for both traits have rendered co-localized QTL for yield and resistance, indicating the presence of pleiotropism or linkage between genes contributing to the negative relationship between resistance and yield (Ordás et al., 2010; Jiménez-Galindo et al., 2017).
An alternative to PS could be Marker Assisted Selection (MAS) using markers linked to QTL with significant effects on resistance to stem tunneling by MCB but without adverse effects on yield (Ordás et al., 2010; Samayoa et al., 2019). However, commonly employed MAS strategies are not suitable for improving complex agronomic traits, in part, due to the lack of power to capture small effect loci that are often important in the inheritance of these traits since many minor genes are involved in MCB resistance, measured as reduced pith gallery length, and yield (Butrón et al., 2009; Butrón et al., 2012). For this reason, genomic selection (GS), which uses all marker information to calculate genomic estimated breeding values (GEBV), is presented as an emerging alternative to MAS for improving complex traits; genome-wide molecular marker coverage is able to capture both large and small effect QTL. It is also worth noting that genotyping and sequencing technology is advancing at an extremely rapid pace, with reducing cost. This coupled with continued advances in computational techniques to predict GEBV has great potential to increase accuracy at little or no additional cost; whereas phenotyping costs are stagnant or increasing. Therefore, GS would allow the design of selection programs with a higher genetic gain per year and a budget equivalent to that of a MAS program (Heffner et al., 2010; Zhang et al., 2017; Rice and Lipka, 2021). However, these predictions are based on simulation studies and there are few GS programs to validate such predictions (Heffner et al., 2010; Rice and Lipka, 2021). The assessment of genetic gain using GS has been exemplified in maize by rapid cycling recombination of biparental populations for drought tolerance (Beyene et al., 2015; Vivek et al., 2017),multiparent populations for grain yield under different moisture regimes (Zhang et al., 2017; Das et al., 2020) and resistance to F. verticillioides using a diverse set of tropical germplasm (Kuki et al., 2020), but no GS programs have been carried out to improve resistance to pests.
In this study, two selection cycles of six selection programs have been carried out to: 1) reduce the length of the tunnels produced by borers using the phenotypic value as selection criteria; 2) reduce the length of the tunnels using the GEBV; 3) improve yield using the phenotypic value; 4) yield improvement based on GEBV; 5) simultaneously improve both traits (tunnels and yield) using the phenotypic values of the index; and 6) improve yield and tunnels using GEBV of the index. The objective was to test whether SG, based on GEBV, is more effective than SP, based on phenotypic values, to improve resistance (short tunnels) using program 1 versus 2; tolerance (high yield under high pest pressure conditions) comparing programs 3 and 4, and/or tolerance and resistance simultaneously based on the comparison of programs 5 and 6.
2 Materials and methods
A population of F6 recombinant lines (RILs) was constructed from an F2 derived from the cross of two inbred lines (A637 × EP42) with contrasting tolerance to MCB. A637 is tolerant to MCB on the basis that it shows low yield losses caused by MCB damage, while EP42 is sensitive to MCB based on high yield losses caused by MCB damage. However, both parental lines are susceptible because they show a lot of MCB damage measured as tunnel length produced by MCB.
However, both parental lines are susceptible because they show large tunnel lengths by MCB attack. Each F6 was derived from a different F2 plant by using the single-seed descendant method (Butrón et al., 1998; Samayoa et al., 2014).
The phenotypic data for tunnel length produced by MCB and grain yield at 140g of H2O · kg-1, obtained by Samayoa et al. (2014), were used to carry out the phenotypic selection in the (A637 x EP42) RIL populations (original cycle, C0) and, along with genotyping data, to calculate the GEBV that were the basis of genomic selection.
The hundred forty-four RILs were previously genotyped with 130 polymorphic SSRs (Samayoa et al., 2014). These markers were relevant and informative and served to create a linkage map using a LOD threshold of two to declare significantly linked two markers, and a maximum distance of 50 cM among adjacent markers(Samayoa et al., 2014). GEBV for tunnel length and grain yield were estimated via ridge regression-BLUP (RR-BLUP) following Bernardo and Yu (2007). The performance of the RILs was modeled as:
where y is an 144 × 1 vector of RIL phenotypic means (tunnel length, grain yield); μ is the overall mean of the 144 RILs; 1 is an 144 × 1 vector with all elements equal to 1; X is a 144×130 matrix with elements equal to 1 or -1 if RIL is homozygous for SSR allele from A637 (first parental inbred) or from EP42 (second parental inbred), respectively; g is an 130 × 1 breeding value vector associated with the SSR allele from A637 at each of the SSR loci and e is a vector (144 × 1) of residual effects.
The RILs evaluation trials carried out by Samayoa et al. (2014) were used to estimate the GV (genetic variance) and EV (environmental variance) in the original cycle C0. The variance of breeding values at each of the SSR loci was assumed equal to GV/SSR Number. GEBV (g) were obtained by solving the mixed-model equations being μ a fixed effect and g a random effect. SAS/IML software was used for the calculations (SAS 9.4, SAS Institute 2016).
Phenotypic data were used for developing three independent selection programs with three different selection criteria: 1) highest resistance (lower tunnels lengths produced by MCB) (PS-tunlen); 2) highest grain yield (PS-yield); and 3) the most resistant lines among those with a yield above the average (PS-index). For each selection program 20 RILs were chosen (selection intensity of 14%) and recombined to constitute the cycle 1 of selection. On the other hand, the estimated GEBVs were used to carry out three independent genomic selection programs using the same selection criteria and selection intensity than in phenotypic selection: GS-tunlen, GS-yield and GS-index. Three groups of 20 RILs were randomly selected for checking for genetic drift effects. Randomly selected RILs were recombined to constitute the RandomC1.
The recombination in each of the populations were carried out in 2013 by sowing a balanced number of grains from each of the constituent RILs and making plant-to-plant crosses, where each plant is involved in only one cross, either as a male or as female. In addition, the original cycle (C0) was reconstructed by recombining the 144 RILs using a 500 kernel-weighted mixture. In 2014, a second recombination of the C0 and the nine C1 populations was carried out.
In 2015, 150 plants from each C1 population were sown and self-fertilized, obtaining about 100 S1 families from each population. For phenotypic selection self -fertilized plants were infested with approx. 40 MCB eggs per plant deposited between the stem and the sheath of one basal leaf. Eggs reared at the MBG were used as described by Eizaguirre and Albajes (1992) and Khan and Saxena (1997). At harvest, the length of the tunnels produced by MCB were measured, the ears were collected and the yield per plant was evaluated using a subjective scale from 1=poor yield to 9=very good yield. The best 20 families from each population were selected and recombined to constitute the second PS cycles. On the other hand, for genomic selection, S1 families from GS-tunlenC1, GS-yieldC1 and GS-indexC1 were genotyped for the 130 SSRs used in the first cycle of genomic selection. The 20 best S1 families from each population families were selected based on their predicted breeding values, that is, Zg, where Z was an N × M matrix that corresponded to the genotypes for M SSR of the N families evaluated in C1s. The elements in Z were zero if the plants in C1 were heterozygous at the corresponding SSR loci. Selected families were recombined to constitute the second GS cycles. Families were randomly selected and recombined to form the RandomC2 cycle. The recombinations were performed in the same way as those that led to the first selection cycles (Figure S1).
In 2017, the first and second cycles of the 9 selections (3 PS, 3 GS and 3 random) were multiplied to obtain homogeneous seeds for the evaluations. In 2018 and 2019, the selection cycles were evaluated together with the original population (C0) at Pontecaldelas (42°23’31.4 “N 8°30’14.3 “W) and Barrantes (42°30’16.1 “N 8°46’16.2 “W) under natural infestation. Both locations during the two years showed a high intensity of MCB attack estimated as the averaged tunnel length across genotypes in each environment (44.3, 46.5, 21.8 and 10.3 cm in Barrantes 2018, Pontecaldelas 2018, Barrantes 2019 and Pontecaldelas 2019, respectively). Randomized complete block designs with three replications were used. Each experimental plot consisted of two rows of 25 plants · row-1 with a distance between rows of 0.8 m and between plants of 0.18 m. Data were recorded for grain moisture (percentage of water in total grain weight at harvest), yield (plot kernel weight at 14 kernel moisture referred in Mg · ha-1) and length of pith stem tunnels produced by MCB. In 2020, evaluations were also carried out at the Misión Biológica de Galicia (MBG) in Pontevedra (42°24’20.0 “N 8°38’41.9 “W) under artificial infestation and under insecticide protection. Eggs reared at the MBG as described by Eizaguirre and Albajes (1992) and Khan and Saxena (1997) were used to infest all plants in one row of the plot using about 40 eggs that were placed between the stem and the sheath of a basal leaf.
Data were analyzed with the Proc Mixed procedure of SAS software (SAS 9.4, SAS Institute 2016). Genotype was considered a fixed factor while replicates, environments (year-location-infestation condition combinations) and interactions between factors were considered random factors. Means were calculated for each genotype by the least squares procedure and compared by Fisher’s protected Least Significant Difference (LSD) at the 0.05 probability level. Finally, data for each trait were adjusted to linear regression models using a repeated measurement approach. Cycles was the independent variable; selection programs [each method (GS, PS or Random) x selection criteria (yield, tunnel length or index) combination] and the interaction cycle x selection programs were stablished as fixed effects and environments, repetitions and their interaction as random effects.
3 Results
A combined analysis over five environments was performed for tunnel length (the trial protected with insecticide was eliminated because no damage caused by MCB was observed). The genotype x environment interaction was not significant either for grain yield or tunnel length (Table 1). There were significant differences between the cycles within selection methods for the tunnel length produced by MCB. PS-yieldC2 and GS-indexC2 were the only ones that differed significantly (p< 0.05) from the original population. PS-yieldC2 was significantly more resistant than the original cycle but neither both cycles of genomic selection for tunnel length (GS-tunlenC1 and GS-tunlenC2) nor cycle 2 of phenotypic selection for tunnel length (PS-tunlenC2) differed from it. On the other hand, GS-indexC2 was the most susceptible selection cycle as it showed the largest tunnel. It should also be noted that there was no significant difference between the initial cycle and random selection; therefore, the effect of genetic random drift seems negligible for length of tunnels produced by MCB (Figure 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1162440/fpls-14-1162440-HTML/image_m/fpls-14-1162440-t001.jpg)
Table 1 Analysis of variance of the cycles of selection obtained by phenotypic and genomic selections for yield, tunnel length or an index that includes both traits.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1162440/fpls-14-1162440-HTML/image_m/fpls-14-1162440-g001.jpg)
Figure 1 Comparison of means for tunnel length (cm) produced by MCB of the original cycle (C0);cycles of PS and GS to reduce tunnel length (tunlen); to improve grain yield (yield) and to reduce tunnel length while maintaining yield (index) and random selection cycles. Different letters indicate that means differed significantly at the 0.05 probability level.
A combined analysis on six environments was also performed for grain yield. There were significant differences between genotypes while the differences were not significant for the interaction (Table 1). The cycle with the highest yield was the second cycle of the GS for yield (GS-yieldC2) but it did not differ significantly from the second cycle of the PS to reduce tunnels (PS-tunlenC2) nor from the original cycle.
Significant difference for yield between the original cycle and the random selection indicate that the effect of genetic drift was important during the first cycle of selection, and therefore the different selection cycles will be compared to the first cycle of random selection to eliminate as much as possible genetic drift effects. Both GS-yieldC2 and PS-tunlenC2 differed significantly from randomC1 for yield (Figure 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1162440/fpls-14-1162440-HTML/image_m/fpls-14-1162440-g002.jpg)
Figure 2 Comparison of grain yield (Mg·ha-1) of the original cycle (C0); cycles of PS and GS to reduce tunnel length (tunlen); to improve grain yield (yield) and to reduce tunnel length while maintaining yield (index) and random selection cycles. Different letters indicate that means differed significantly at the 0.05 probability level.
Linear regressions allow evaluating how a trait varies throughout the cycles of a selection program. Figure 3 shows the linear regressions of the length of tunnels produced by MCB on cycles of selection in seven selection programs. The regression coefficient was significantly different from zero in two selection programs: 1) phenotypic selection to improve yield (PS-yield) and GS to improve resistance while maintaining yield (GS-index). In the first case, the tunnels are reduced as the selection program progresses (b= -3.98) while in GS-index the tunnel lenght increase as the selection progresses (b=3.72). The best linear response to reduce tunnel length was obtained with PS-yield, but it did not differ significantly from those achieved with PS and GS for tunnel length (p = -0.34 for PS-yield vs. PS-tunlen b coefficients and p = 0.28 for PS-yield vs. GS-tunlen b coefficients) and PS-index (p = 0.15 for b coefficient comparison between PS-yield and PS-index).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1162440/fpls-14-1162440-HTML/image_m/fpls-14-1162440-g003.jpg)
Figure 3 Linear regressions analysis for length of tunnels produced by MCB n selection cycles in six selection programs and one random selection. The linear response equation is shown for each of the selections, indicating whether the linear coefficient was significant at 0.05 (*).
The selection programs GS-yield and PS-tunlen showed linear coefficients significantly different from zero (Figure 4). The best linear response was obtained with GS-yield (0.63), but it did not differ significantly from those achieved with PS-tunlen (p = 0.42) and GS-tunlen (p = 0.42).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1162440/fpls-14-1162440-HTML/image_m/fpls-14-1162440-g004.jpg)
Figure 4 Linear regression analysis for yield (Mg·ha-1) of the selection cycles in six selection programs. The linear response equation is shown for each of the selections, indicating whether the linear coefficient was significant at 0.05 (*). RandomC1 was considered as original cycle (C0).
4 Discussion
Two decades of research on maize resistance to MCB has demonstrated the high negative genetic correlation between yield potential and resistance. Such that pursuing higher resistance, measured as reduced tunnel length made by MCB in the pith, has significantly contributed to reduced kernel yield (Ordás et al., 2010; Butrón et al., 2012; Jiménez-Galindo et al., 2017; Samayoa et al., 2019). Although this negative relationship is not found in all genetic backgrounds (Samayoa et al., 2014; Jiménez-Galindo et al., 2019).
Some studies show that the inheritance of resistance to stem tunneling is basically under additive control by observing that improvement achieved through inbreeding is generally transmitted to hybrids (Butrón et al., 1999; Cartea et al., 1999). However, genotype differences for tunnel length can be biased by MCB genotype and MCB genotype × environment differences among plots, so the experimental error associated to tunnel length is high and makes that heritability for tunnel length is often low (Samayoa et al., 2014; Jiménez-Galindo et al., 2017). Moreover, resistance to MCB is a highly polygenic trait, where many genes are involved with small effects that are often undetected in QTL studies (Ordás et al., 2010; Samayoa et al., 2014; Samayoa et al., 2015; Jiménez-Galindo et al., 2017; Samayoa et al., 2019). The GS strategy based on genome-wide polymorphic markers that can capture minor gene effects on complex traits is presented as a possible alternative to phenotype and marker-assisted selections (Heffner et al., 2010; Liu et al., 2018; Rice and Lipka, 2021). Our results indicate that the GS-tunlen C1 and C2 and PS-tunlenC2 cycles did not differ from the cycle with the shortest tunnels (PS-yieldC2), confirming that genetic variability for resistance to stem tunneling can be managed by phenotypic and genomic selections, although the gains did not differ significantly from zero. Although the additive component was the most important for tunnel length among generations derived from A637 × EP42 (Butrón et al., 2009), heritability in the (A637 × EP42) F2 population was low, h2 = 0.34 (Samayoa et al., 2014) which could explain the lack of a significant improvement by PS. GS was also ineffective to significantly reduce tunnel length likely due to the fact that genotype predictions are often not highly accurate when non-genetic factors are a major source of trait variation as it is the case for tunnel length due to a high experimental error (Hickey et al., 2017). However, genomic prediction accuracy could be improved by increasing the number of markers because less than 200 markers could be enough to achieve the highest prediction accuracy for traits with moderate heritability in F2 and backcrossing populations (Lorenzana and Bernardo, 2009; Lian et al., 2014; Zhang et al., 2015), but Hao et al. (2019) estimated that the highest accuracy was achieved with 800 SNPs when genomic selection is performed in a RIL population because linkage disequilibrium is lower among RILs compared to F2 and backcrossing mapping populations.
PS-yieldC2 was the only selection scheme that significantly reduced tunnel length compared to the original cycle. Phenotypic selection for yield could improve resistance because, in this material, the genetic correlation between resistance and yield seems important and positive, contrary to other genetic backgrounds where this correlation is negative (Butrón et al., 2012). Another aspect that could contribute to the success of PS-yield for reducing tunnel length is the method used to select the S1 families that formed the PS-yieldC2. The ears were not weighed because they come from self-pollinated plants and the amount of grain depends on the moment of pollination; instead, ear yield was assessed by a subjective scale and good looking ears would probably come from healthier plants, tunnels by MCB being determinant in plant health.
Inheritance of yield has an important dominant component, which could lead to a clear inbreeding depression effect when few individuals are selected in a highly variable population. The yield of RandomC1 was significantly lower than that of the original population (4.32 vs. 5.37 ± 0.38 Mg·ha-1) confirming that the selection scheme used generated genetic drift that, in turn, originated inbreeding depression on yield. Therefore, in order to correct for the effect of genetic drift on yield, the RandomC1 cycle was used as the starting point to calculate yield response to selection. GS-yield was the most effective method for increasing yield; meanwhile PS was ineffective to improve the yield. Given that the inheritance of the yield has an important dominant component, the best phenotypes could be those with more favorable allelic interactions for yield. While GS was based on additive variation, which is the heritable genetic component, and was much more effective (Liu et al., 2018). In this work, it has been verified, as stated before by other authors (Zhang et al., 2017; Rice and Lipka, 2021), that GS is a powerful tool for improving the yield of maize.
Selections using indices, both GS and PS, have not been successful in improving resistance to MCB while maintaining yield. Indeed, GS for index values attained a significant increase of tunnels after two cycles of selection and no changes were observed in yield after GS and PS selections. Phenotypic selection for index values did not get a significant gain for tunnel length contrarily to index selection made previously in a maize composite named EPS12 (Sandoya et al., 2008; Butrón et al., 2012). These authors reported that the length of the tunnels and yield had been significantly reduced using the same index selection approach. The lower response to index phenotypic selection observed in the current study would likely be related to the narrowed variability for resistance to MCB in the (A637 × EP42) F2 population compared to the EPS12 composite since A637 and EP42 were both susceptible to MCB attack (Butrón et al., 1998).
In summary, the linear response of yield to GS-yield did not differ significantly from responses to GS and PS selections for tunnel length. Similarly, PS-yield selection, which showed the highest tunnel length reduction, presented a linear response on tunnel length that did not differ significantly from those obtained by GS and PS selections. Suggesting that phenotypic or genotypic selection for tunnel length may be more effective in improving both traits, but the gains per cycle would be small for both traits. However, from a practical point of view, high yield, regardless of the damage caused by MCB, would be the most valuable trait by farmers and genomic selection for yield has been proven as the most successful methodology to increase yield.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
RM and BO conceived the study and performed selection programs with the assistance of AB, JJ-G and LS. MF coordinated and carried out the genotyping work to carry out the genomic selection with the collaboration of AC. RM, JJ-G, LS, and AB took care of field experiments, data recording, and sample collection. RM performed statistical analyses of data with the assistance of NG. NG drafted and edited the manuscript. All authors have read and approved the final version of the manuscript. All authors contributed to the article.
Funding
This research was funded by subsequent coordinated projects financed by MCIU/AEI/FEDER, UE (RTI2018-096776-B-C21, RTI2018-096776-B-C22, PID2021-122196OB-C21 and PID2021-122196OB-C22). We also used funds of the CSIC project 2021AEP132. NG acknowledges an FPI contract (“Contrato predoctoral para la formación de doctores”) funded by project RTI2018-096776-BC21 (financed by MCIU/AEI/FEDER, UE).
Conflict of interest
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1162440/full#supplementary-material
Supplementary Figure 1 | Summary figure of breeding programs based on phenotypic or genotypic values for a single trait (yield or MCB susceptibility measured as tunnel length) or both traits (index).
References
Bernardo, R., Yu, J. M. (2007). Prospects for genome wide selection for quantitative traits in maize. Crop Sci. 47 (3), 1082–1090. doi: 10.2135/cropsci2006.11.0690
Beyene, Y., Semagn, K., Mugo, S., Tarekegne, A., Babu, R., Meisel, B., et al. (2015). Genetic gains in grain yield through genomic selection in eight bi-parental maize populations under drought stress. Crop Sci. 55 (1), 154–163. doi: 10.2135/cropsci2014.07.0460
Butrón, A., Malvar, R. A., Cartea, M. E., Ordás, A., Velasco, P. (1999). Resistance of maize inbreds to pink stem borer. Crop Sci. 39 (1), 102–107. doi: 10.2135/cropsci1999.0011183X003900010016x
Butrón, A., Malvar, R. A., Revilla, P., Ordás, R. A. (1998). Defense mechanisms of maize against pink stem borer. Crop Sci. 38 (5), 1159–1163. doi: 10.2135/cropsci1998.0011183X003800050008x
Butrón, A., Romay, M. C., Peña-Asín, J., Alvarez, A., Malvar, R. A. (2012). Genetic relationship between maize resistance to corn borer attack and yield. Crop Sci. 52 (3), 1176–1180. doi: 10.2135/cropsci2011.11.0584
Butrón, A., Samayoa, F., Santiago, R., Malvar, R. A. (2014). Selection efficiency of tunnel length and stalk breakage to obtain maize inbred lines resistant to stem borer attack. Euphytica 197 (2), 295–302. doi: 10.1007/s10681-014-1067-7
Butrón, A., Sandoya, G., Revilla, P., Malvar, R. A. (2009). Genetics of resistance to the pink stem borer (Sesamia nonagrioides) in maize (Zea mays). Ann. Appl. Biol. 154 (2), 205–217. doi: 10.1111/j.1744-7348.2008.00284.x
Butrón, A., Sandoya, G., Revilla, P., Ordás, A., Malvar, R. A. (2005). Evaluation of popcorn germplasm for resistance to Sesamia nonagrioides attack. J. Econ. Entomol. 98 (5), 1694–1697. doi: 10.1093/jee/98.5.1694
Butrón, A., Sandoya, G., Santiago, R., Ordás, A., Rial, A., Malvar, R. A. (2006). Searching for new sources of pink stem borer resistance in maize (Zea mays l.). Genet. Resour. Crop Evol. 53 (7), 1455–1462. doi: 10.1007/s10722-005-6842-9
Cartea, M. E., Malvar, R. A., Butrón, A., Vales, M. I., Ordás, A. (1999). Inheritance of antibiosis to Sesamia nonagrioides (Lepidoptera: noctuidae) in maize. J. Econ. Entomol. 92 (4), 994–998.
Das, R. R., Vinayan, M. T., Patel, M. B., Phagna, R. K., Singh, S. B., Shahi, J. P., et al. (2020). Genetic gains with rapid-cycle genomic selection for combined drought and waterlogging tolerance in tropical maize (Zea mays l.). Plant Genome 13 (3), 1–15. doi: 10.1002/tpg2.20035
Eizaguirre, M., Albajes, R. (1992). Diapause induction in the stem corn borer, Sesamia nonagrioides (Lepidoptera: noctuidae). Entomol. Gen. 17 (4), 277–283. doi: 10.1127/entom.gen/17/1992/277
Hao, Y., Wang, H., Yang, X., Zhang, H., He, C., Li, D., et al. (2019). Genomic prediction using existing historical data contributing to selection in biparental populations: a study of kernel oil in maize. Plant Genome 12, 180025. doi: 10.2135/cropsci2009.11.0662
Heffner, E. L., Lorenz, A. J., Jannink, J. -L., Sorrells, M. E. (2010). Plant breeding with genomic selection: gain per unit time and cost. Crop Sci. 50, 1681–1690. doi: 10.2135/cropsci2009.11.0662
Hickey, J. M., Chiurugwi, T., Mackay, I., Powell, W. (2017). Genomic prediction unifies animal and plant breeding programs to form platforms for biological discovery. Nat. Genet. 49 (9), 1297–1303. doi: 10.1038/ng.3920
Jika, A. K. N., Ru, B., Capdevielle-Dulac, C., Chardonnet, F., Silvain, J. F., Kaiser, L., et al. (2020). Population genetics of the mediterranean corn borer (Sesamia nonagrioides) differs between wild and cultivated plants. PLoS One 15 (3), 1–17. doi: 10.1371/journal.pone.0230434
Jiménez-Galindo, J. C., Malvar, R. A., Butrón, A., Santiago, R., Samayoa, L. F., Caicedo, M., et al. (2019). Mapping of resistance to corn borers in a MAGIC population of maize. BMC Plant Biol. 19 (1), 431. doi: 10.1186/s12870-019-2052-z
Jiménez-Galindo, J. C., Ordás, B., Butrón, A., Samayoa, L. F., Malvar, R. A. (2017). QTL mapping for yield and resistance against Mediterranean corn borer in maize. Front. Plant Sci. 8 (May). doi: 10.3389/fpls.2017.00698
Khan, Z. R., Saxena, R. C. (1997). Use of a surrogate stem for eliciting ovipositional response of Busseola fusca (Lepidoptera: noctuidae). J. Econ. Entomol. 90 (5), 1426–1429. doi: 10.1093/jee/90.5.1426
Kuki, M. C., Pinto, R. J. B., Bertagna, F. A. B., Tessmann, D. J., Teixeira, A., Scapim, C. A., et al. (2020). Association mapping and genomic prediction for ear rot disease caused by fusarium verticillioides in a tropical maize germplasm. Crop Sci. 60 (6), 2867–2881. doi: 10.1002/csc2.20272
Lian, L., Jacobson, A., Zhong, S., Bernardo, R. (2014). Genomewide prediction accuracy within 969 maize biparental populations. Crop Sci. 54, 1514–1522. doi: 10.2135/cropsci2013.12.0856
Liu, X., Wang, H., Wang, H., Guo, Z., Xu, X., Liu, J., et al. (2018). Factors affecting genomic selection revealed by empirical evidence in maize. Crop J. 6 (4), 341–352. doi: 10.1016/j.cj.2018.03.005
Lorenzana, R. E., Bernardo, R. (2009). Accuracy of genotypic value predictions for marker-based selection in biparental plant populations. Theor. Appl. Genet. 120, 151–161. doi: 10.1007/s00122-009-1166-3
Malvar, R. A., Butrón, A., Alvarez, A., Ordás, B., Soengas, P., Revilla, P., et al. (2004). Evaluation of the European union maize landrace core collection for resistance to Sesamia nonagrioides (Lepidoptera: noctuidae) and Ostrinia nubilalis (Lepidoptera: crambidae). J. Econ. Entomol. 97 (2), 628–634. doi: 10.1603/0022-0493-97.2.628
Malvar, R. A., Butrón, A., Alvarez, A., Padilla, G., Cartea, M. E., Revilla, P., et al. (2007). Yield performance of the European union maize landrace core collection under multiple corn borer infestations. Crop Prot. 26 (5), 775–781. doi: 10.1016/j.cropro.2006.07.004
Malvar, R. A., Butrón, A., Ordás, B., Santiago, R. (2008). “Causes of natural resistance to stem borers in maize,” in Crop protection research advances. Eds. Burton, E. N., Williams, P. V. (Nova Science Publishers, Inc), 57–100.
Malvar, R. A., Cartea, M. E., Revilla, P., Ordás, A., Alvarez, A., Mansilla, J. P. (1993). Sources of resistance to pink stem borer and European corn-borer in maize. Maydica 38 (4), 313–319.
Ordás, B., Butrón, A., Alvarez, A., Revilla, P., Malvar, R. A. (2012). Comparison of two methods of reciprocal recurrent selection in maize (Zea mays l.). Theor. Appl. Genet. 124 (7), 1183–1191. doi: 10.1007/s00122-011-1778-2
Ordás, B., Malvar, R. A., Santiago, R., Butrón, A. (2010). QTL mapping for Mediterranean corn borer resistance in European flint germplasm using recombinant inbred lines. BMC Genomics 11 (1). doi: 10.1186/1471-2164-11-174
Revilla, P., Butrón, A., Ordás, B., Soengas, P., Ordás, A., Malvar, R. A. (2007). Corn borer (Lepidoptera: noctuidae and crambidae) resistance of main races of maize from north America. J. Econ. Entomol. 100 (1), 209–214. doi: 10.1093/jee/100.1.209
Rice, B., Lipka, A. (2021). Diversifying maize genomic selection models. Mol. Breed. 41. doi: 10.1007/s11032-021-01221-4
Samayoa, L. F., Butrón, A., Malvar, R. A. (2014). QTL mapping for maize resistance and yield under infestation with Sesamia nonagrioides. mol. Breed. 34 (3), 1331–1344. doi: 10.1007/s11032-014-0119-y
Samayoa, L. F., Butrón, A., Malvar, R. A. (2019). Usefulness of marker-assisted selection to improve maize for increased resistance to Sesamia nonagrioides attack with no detrimental effect on yield. Ann. Appl. Biol. 174 (2), 219–222. doi: 10.1111/aab.12480
Samayoa, L. F., Butrón, A., Revilla, P., Alvarez, A., Malvar, R. A. (2012). Five cycles of mass selection for earliness and ear appearance under corn borer infestation in the maize synthetic BS17. Crop Sci. 52 (6), 2432–2437. doi: 10.2135/cropsci2011.12.0651
Samayoa, L. F., Malvar, R. A., McMullen, M. D., Butrón, A. (2015). Identification of QTL for resistance to Mediterranean corn borer in a maize tropical line to improve temperate germplasm. BMC Plant Biol. 15 (1). doi: 10.1186/s12870-015-0652-9
Sandoya, G., Butrón, A., Alvarez, A., Ordás, A., Malvar, R. A. (2008). Direct response of a maize synthetic to recurrent selection for resistance to stem borers. Crop Sci. 48 (1), 113–118. doi: 10.2135/cropsci2007.02.0084
Sandoya, G., Malvar, R. A., Santiago, R., Alvarez, A., Revilla, P., Butrón, A. (2010). Effects of selection for resistance to Sesamia nonagrioides on maize yield, performance and stability under infestation with Sesamia nonagrioides and Ostrinia nubilalis in Spain. Ann. Appl. Biol. 156 (3), 377–386. doi: 10.1111/j.1744-7348.2010.00395.x
Velasco, P., Malvar, R. A., Butrón, A., Revilla, P., Ordás, A. (1999b). Ear feeding resistance of sweet corn inbreds to pink stem borer. J. Am. Soc Hortic. Sci. 124 (3), 268–272.
Velasco, P., Malvar, R. A., Revilla, P., Butrón, A., Ordás, A. (1999a). Ear resistance of sweet corn populations to Sesamia nonagrioides (Lepidoptera : noctuidae) and Ostrinia nubilalis (Lepidoptera: pyralidae). J. Econ. Entomol. 92 (3), 732–739.
Velasco, P., Revilla, P., Cartea, M. E., Ordás, A., Malvar, R. A. (2004a). Resistance of early maturing sweet corn varieties to damage caused by Sesamia nonagrioides (Lepidoptera: noctuidae). J. Econ. Entomol. 97 (4), 1432–1437. doi: 10.1603/0022-0493-97.4.1432
Velasco, P., Revilla, P., Malvar, R. A., Butrón, A., Ordás, A. (2002). Resistance to corn borer in crosses between sweet and field corn populations. J. Am. Soc Hortic. Sci. 127 (4), 689–692.
Velasco, P., Revilla, P., Monetti, L., Butron, A., Ordás, A., Malvar, R. A. (2007). Corn borers (Lepidoptera: noctuidae; crambidae) in northwestern Spain: population dynamics and distribution. Maydica 52 (2), 195–203.
Velasco, P., Soengas, P., Revilla, P., Ordás, A., Malvar, R. A. (2004b). Mean generation analysis of the damage caused by Sesamia nonagrioides (Lepidoptera: noctuidae) and Ostrinia nubilalis (Lepidoptera: crambidae) in sweet corn ears. J. Econ. Entomol. 87 (1), 120–126. doi: 10.1603/0022-0493-97.1.120
Vivek, B. S., Krishna, G. K., Vengadessan, V., Babu, R., Zaidi, P. H., Kha, L. Q., et al. (2017). Use of genomic estimated breeding values results in rapid genetic gains for drought tolerance in maize. TPG 10 (1), 1–8. doi: 10.3835/plantgenome2016.07.0070
Zhang, X., Pérez-Rodríguez, P., Burgueño, J., Olsen, M., Buckler, E., Atlin, G., et al. (2017). Rapid cycling genomic selection in a multiparental tropical maize population. G3-Genes Genom.Genet. 7 (7), 2315–2326. doi: 10.1534/g3.117.043141
Keywords: maize, stem borer, resistance, yield, genomic selection, phenotypic selection
Citation: Gesteiro N, Ordás B, Butrón A, de la Fuente M, Jiménez-Galindo JC, Samayoa LF, Cao A and Malvar RA (2023) Genomic versus phenotypic selection to improve corn borer resistance and grain yield in maize. Front. Plant Sci. 14:1162440. doi: 10.3389/fpls.2023.1162440
Received: 09 February 2023; Accepted: 22 June 2023;
Published: 07 July 2023.
Edited by:
Zhenyu Jia, University of California, Riverside, United StatesReviewed by:
Rodrigo Iván Contreras-Soto, Universidad de O’Higgins, ChileQin Yang, Northwest A&F University, China
Copyright © 2023 Gesteiro, Ordás, Butrón, de la Fuente, Jiménez-Galindo, Samayoa, Cao and Malvar. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rosa Ana Malvar, rmalvar@mbg.csic.es