Variation in shoot architecture traits and their relationship to canopy coverage and light interception in soybean (Glycine max).

BACKGROUND: In soybeans, faster canopy coverage (CC) is a highly desirable trait but a fully covered canopy is unfavorable to light interception at lower levels in the canopy with most of the incident radiation intercepted at the top of the canopy. Shoot architecture that influences CC is well studied in crops such as maize and wheat, and altering architectural traits has resulted in enhanced yield. However, in soybeans the study of shoot architecture has not been as extensive. RESULTS: This study revealed significant differences in CC among the selected soybean accessions. The rate of CC was found to decrease at the beginning of the reproductive stage (R1) followed by an increase during the R2-R3 stages. Most of the accessions in the study achieved maximum rate of CC between R2-R3 stages. We measured Light interception (LI), defined here as the ratio of Photosynthetically Active Radiation (PAR) transmitted through the canopy to the incoming PAR or the radiation above the canopy. LI was found to be significantly correlated with CC parameters, highlighting the relationship between canopy structure and light interception. The study also explored the impact of plant shape on LI and CO2 assimilation. Plant shape was characterized into distinct quantifiable parameters and by modeling the impact of plant shape on LI and CO2 assimilation, we found that plants with broad and flat shapes at the top maybe more photosynthetically efficient at low light levels, while conical shapes were likely more advantageous when light was abundant. Shoot architecture of plants in this study was described in terms of whole plant, branching and leaf-related traits. There was significant variation for the shoot architecture traits between different accessions, displaying high reliability. We found that that several shoot architecture traits such as plant height, and leaf and internode-related traits strongly influenced CC and LI. CONCLUSION: In conclusion, this study provides insight into the relationship between soybean shoot architecture, canopy coverage, and light interception. It demonstrates that novel shoot architecture traits we have defined here are genetically variable, impact CC and LI and contribute to our understanding of soybean morphology. Correlations between different architecture traits, CC and LI suggest that it is possible to optimize soybean growth without compromising on light transmission within the soybean canopy. In addition, the study underscores the utility of integrating low-cost 2D phenotyping as a practical and cost-effective alternative to more time-intensive 3D or high-tech low-throughput methods. This approach offers a feasible means of studying basic shoot architecture traits at the field level, facilitating a broader and efficient assessment of plant morphology.

Cataloging SCN resistance loci in North American public soybean breeding programs.

Soybean cyst nematode (SCN) is a destructive pathogen of soybeans responsible for annual yield loss exceeding $1.5 billion in the United States. Here, we conducted a series of genome-wide association studies (GWASs) to understand the genetic landscape of SCN resistance in the University of Missouri soybean breeding programs (Missouri panel), as well as germplasm and cultivars within the United States Department of Agriculture (USDA) Uniform Soybean Tests-Northern Region (NUST). For the Missouri panel, we evaluated the resistance of breeding lines to SCN populations HG 2.5.7 (Race 1), HG 1.2.5.7 (Race 2), HG 0 (Race 3), HG 2.5.7 (Race 5), and HG 1.3.6.7 (Race 14) and identified seven quantitative trait nucleotides (QTNs) associated with SCN resistance on chromosomes 2, 8, 11, 14, 17, and 18. Additionally, we evaluated breeding lines in the NUST panel for resistance to SCN populations HG 2.5.7 (Race 1) and HG 0 (Race 3), and we found three SCN resistance-associated QTNs on chromosomes 7 and 18. Through these analyses, we were able to decipher the impact of seven major genetic loci, including three novel loci, on resistance to several SCN populations and identified candidate genes within each locus. Further, we identified favorable allelic combinations for resistance to individual SCN HG types and provided a list of available germplasm for integration of these unique alleles into soybean breeding programs. Overall, this study offers valuable insight into the landscape of SCN resistance loci in U.S. public soybean breeding programs and provides a framework to develop new and improved soybean cultivars with diverse plant genetic modes of SCN resistance.

A genomic analysis of the University of Nebraska Replicated Recurrent Selection program.

The inbred-hybrid system of maize breeding closely resembles a reciprocal full-sib (RFS) selection program. Studying changes in genetic variation as a result of RFS selection can help illuminate long-standing questions regarding the relative roles of selection and genetic drift and help understand the nature of adaptation occurring in selection programs. The University of Nebraska-Lincoln Replicated Recurrent Selection (UNL-RpRS) program underwent eight cycles of replicated RFS and S1-progeny selection, making it a powerful system to study genomic changes accompanying selection for inter-population performance. The objectives of this study were to identify regions of the genome under selection after eight cycles of selection and evaluate the effect eight cycles of selection for inter-population full-sib performance had in expanding genome-wide and localized population structure. We address these questions with a large set of individuals sampled from the UNL-RpRS program with dense genotyping-by-sequence data. We found evidence of parallel selection signatures in the UNL-RpRS program, with a region on chromosome 7 being implicated in three of the four selection systems studied. Regions that displayed selection signatures across independently run selection programs represent regions likely to be capitalizing on standing genetic variation and support a soft sweep model of adaptation. We did not find selection to be a strong force in diverging populations undergoing RFS. This could be due to the nature of adaptation occurring in these populations, underlying gene action, or a result of unstable genetic topographies.

A genome-wide analysis of the USDA Soybean Isoline Collection.

The USDA Soybean Isoline Collection has been an invaluable resource for the soybean genetics and breeding community. This collection, established in 1972, consists of 611 near-isogenic lines (NILs) carrying one or multiple genes conferring traits that had been determined to exhibit Mendelian inheritance. It has been used in multiple studies on the genetic basis, physiology, and agronomy of these qualitative traits. Here, we used publicly available genotype (SoySNP50K), phenotype, and pedigree data on this collection to characterize the isogenicity of the NILs and identify chromosomal positions of unmapped genes. A total of 368 NILs had at least 80% identity to their recurrent parent and, thus, were useful for what can be called introgression mapping. Both on-target and off-target introgressions were evaluated. The size of on-target introgressions into individual NILs ranged from 61 kb to 8.4 Mb, whereas off-target introgressions ranged from 2.6 kb to 54.8 Mb. The observed large off-target introgressions indicated that some NILs carry introgressions nearly the size of an entire chromosome. By applying introgression mapping to genes that had never been mapped, we identified the likely chromosomal positions of six such genes: ab, im, lo, Np, pc, and Rpm. The size of mapping intervals was large in some cases (10.28 Mb for im) but small in others (0.21 Mb for Np). The results reported herein will provide future researchers with a resource to help select informative NILs for future studies, and provide a starting point to further fine map, and ultimately clone and functionally characterize these six soybean genes.

Branch angle and leaflet shape are associated with canopy coverage in soybean.

Early canopy coverage is a desirable trait that is a major determinant of yield in soybean (Glycine max). Variation in traits comprising shoot architecture can influence canopy coverage, canopy light interception, canopy-level photosynthesis, and source-sink partitioning efficiency. However, little is known about the extent of phenotypic diversity of shoot architecture traits and their genetic control in soybean. Thus, we sought to understand the contribution of shoot architecture traits to canopy coverage and to determine the genetic control of these traits. We examined the natural variation for shoot architecture traits in a set of 399 diverse maturity group I soybean (SoyMGI) accessions to identify relationships between traits, and to identify loci that are associated with canopy coverage and shoot architecture traits. Canopy coverage was correlated with branch angle, number of branches, plant height, and leaf shape. Using previously collected 50K single nucleotide polymorphism data, we identified quantitative trait locus (QTL) associated with branch angle, number of branches, branch density, leaflet shape, days to flowering, maturity, plant height, number of nodes, and stem termination. In many cases, QTL intervals overlapped with previously described genes or QTL. We also found QTL associated with branch angle and leaflet shape located on chromosomes 19 and 4, respectively, and these QTL overlapped with QTL associated with canopy coverage, suggesting the importance of branch angle and leaflet shape in determining canopy coverage. Our results highlight the role individual architecture traits play in canopy coverage and contribute information on their genetic control that could help facilitate future efforts in their genetic manipulation.

Genomic selection performs as effectively as phenotypic selection for increasing seed yield in soybean.

Increasing the rate of genetic gain for seed yield remains the primary breeding objective in both public and private soybean [Glycine max (L.) Merr.] breeding programs. Genomic selection (GS) has the potential to accelerate the rate of genetic gain for soybean seed yield. Limited studies to date have validated GS accuracy and directly compared GS with phenotypic selection (PS), and none have been reported in soybean. This study conducted the first empirical validation of GS for increasing seed yield using over 1,500 lines and over 7 yr (2010-2016) of replicated experiments in the University of Nebraska-Lincoln soybean breeding program. The study was designed to capture the varying genetic relatedness of the training population to three validation sets: two large biparental populations (TBP-1 and TBP-2) and a large validation set comprised of 457 preselected advanced lines derived from 45 biparental populations (TMP). We found that prediction accuracy (.54) realized in our validation experiments was comparable with what we obtained from a series of cross-validation experiments (.64). Both GS and PS were more effective for increasing the population mean performance compared with random selection (RS). We found a selection advantage of GS over PS, where higher genetic gain and identification of top-performing lines was maximized at 10-20% selected proportion. Genomic selection led to small increases in genetic similarity when compared with PS and RS presumably because of a significant shift on allelic frequencies toward the extremes, suggesting that it could erode genetic diversity more quickly. Overall, we found that GS can perform as effectively as PS but that measures should be considered to protect against loss of genetic variance when using GS.

Plant Breeding for Intercropping in Temperate Field Crop Systems: A Review.

Monoculture cropping systems currently dominate temperate agroecosystems. However, intercropping can provide valuable benefits, including greater yield stability, increased total productivity, and resilience in the face of pest and disease outbreaks. Plant breeding efforts in temperate field crops are largely focused on monoculture production, but as intercropping becomes more widespread, there is a need for cultivars adapted to these cropping systems. Cultivar development for intercropping systems requires a systems approach, from the decision to breed for intercropping systems through the final stages of variety testing and release. Design of a breeding scheme should include information about species variation for performance in intercropping, presence of genotype × management interaction, observation of key traits conferring success in intercropping systems, and the specificity of intercropping performance. Together this information can help to identify an optimal selection scheme. Agronomic and ecological knowledge are critical in the design of selection schemes in cropping systems with greater complexity, and interaction with other researchers and key stakeholders inform breeding decisions throughout the process. This review explores the above considerations through three case studies: (1) forage mixtures, (2) perennial groundcover systems (PGC), and (3) soybean-pennycress intercropping. We provide an overview of each cropping system, identify relevant considerations for plant breeding efforts, describe previous breeding focused on the cropping system, examine the extent to which proposed theoretical approaches have been implemented in breeding programs, and identify areas for future development.

Mining Fiskeby III and Mandarin (Ottawa) Expression Profiles to Understand Iron Stress Tolerant Responses in Soybean.

The soybean (Glycine max L. merr) genotype Fiskeby III is highly resistant to a multitude of abiotic stresses, including iron deficiency, incurring only mild yield loss during stress conditions. Conversely, Mandarin (Ottawa) is highly susceptible to disease and suffers severe phenotypic damage and yield loss when exposed to abiotic stresses such as iron deficiency, a major challenge to soybean production in the northern Midwestern United States. Using RNA-seq, we characterize the transcriptional response to iron deficiency in both Fiskeby III and Mandarin (Ottawa) to better understand abiotic stress tolerance. Previous work by our group identified a quantitative trait locus (QTL) on chromosome 5 associated with Fiskeby III iron efficiency, indicating Fiskeby III utilizes iron deficiency stress mechanisms not previously characterized in soybean. We targeted 10 of the potential candidate genes in the Williams 82 genome sequence associated with the QTL using virus-induced gene silencing. Coupling virus-induced gene silencing with RNA-seq, we identified a single high priority candidate gene with a significant impact on iron deficiency response pathways. Characterization of the Fiskeby III responses to iron stress and the genes underlying the chromosome 5 QTL provides novel targets for improved abiotic stress tolerance in soybean.

The utility of genomic prediction models in evolutionary genetics.

Variation in complex traits is the result of contributions from many loci of small effect. Based on this principle, genomic prediction methods are used to make predictions of breeding value for an individual using genome-wide molecular markers. In breeding, genomic prediction models have been used in plant and animal breeding for almost two decades to increase rates of genetic improvement and reduce the length of artificial selection experiments. However, evolutionary genomics studies have been slow to incorporate this technique to select individuals for breeding in a conservation context or to learn more about the genetic architecture of traits, the genetic value of missing individuals or microevolution of breeding values. Here, we outline the utility of genomic prediction and provide an overview of the methodology. We highlight opportunities to apply genomic prediction in evolutionary genetics of wild populations and the best practices when using these methods on field-collected phenotypes.

Optimization of training sets for genomic prediction of early-stage single crosses in maize.

KEY MESSAGE: Training population optimization algorithms are useful for efficiently training genomic prediction models for single-cross performance, especially if the population is extended beyond only realized crosses to all possible single crosses. Genomic prediction of single-cross performance could allow effective evaluation of all possible single crosses between all inbreds developed in a hybrid breeding program. The objectives of the present study were to investigate the effect of different levels of relatedness on genomic predictive ability of single crosses, evaluate the usefulness of deterministic formula to forecast prediction accuracy in advance, and determine the potential for TRS optimization based on prediction error variance (PEVmean) and coefficient of determination (CDmean) criteria. We used 481 single crosses made by crossing 89 random recombinant inbred lines (RILs) belonging to the Iowa stiff stalk synthetic group with 103 random RILs belonging to the non-stiff stalk synthetic heterotic group. As expected, predictive ability was enhanced by ensuring close relationships between TRSs and target sets, even when TRS sizes were smaller. We found that designing a TRS based on PEVmean or CDmean criteria is useful for increasing the efficiency of genomic prediction of maize single crosses. We went further and extended the sampling space from that of all observed single crosses to all possible single crosses, providing a much larger genetic space within which to design a training population. Using all possible single crosses increased the advantage of the PEVmean and CDmean methods based on expected prediction accuracy. This finding suggests that it may be worthwhile using an optimization algorithm to select a training population from all possible single crosses to maximize efficiency in training accurate models for hybrid genomic prediction.

Characterizing introgression-by-environment interactions using maize near isogenic lines.

KEY MESSAGE: Significant introgression-by-environment interactions are observed for traits throughout development from small introgressed segments of the genome. Relatively small genomic introgressions containing quantitative trait loci can have significant impacts on the phenotype of an individual plant. However, the magnitude of phenotypic effects for the same introgression can vary quite substantially in different environments due to introgression-by-environment interactions. To study potential patterns of introgression-by-environment interactions, fifteen near-isogenic lines (NILs) with > 90% B73 genetic background and multiple Mo17 introgressions were grown in 16 different environments. These environments included five geographical locations with multiple planting dates and multiple planting densities. The phenotypic impact of the introgressions was evaluated for up to 26 traits that span different growth stages in each environment to assess introgression-by-environment interactions. Results from this study showed that small portions of the genome can drive significant genotype-by-environment interaction across a wide range of vegetative and reproductive traits, and the magnitude of the introgression-by-environment interaction varies across traits. Some introgressed segments were more prone to introgression-by-environment interaction than others when evaluating the interaction on a whole plant basis throughout developmental time, indicating variation in phenotypic plasticity throughout the genome. Understanding the profile of introgression-by-environment interaction in NILs is useful in consideration of how small introgressions of QTL or transgene containing regions might be expected to impact traits in diverse environments.

Dominance Effects and Functional Enrichments Improve Prediction of Agronomic Traits in Hybrid Maize.

Single-cross hybrids have been critical to the improvement of maize (Zea mays L.), but the characterization of their genetic architectures remains challenging. Previous studies of hybrid maize have shown the contribution of within-locus complementation effects (dominance) and their differential importance across functional classes of loci. However, they have generally considered panels of limited genetic diversity, and have shown little benefit from genomic prediction based on dominance or functional enrichments. This study investigates the relevance of dominance and functional classes of variants in genomic models for agronomic traits in diverse populations of hybrid maize. We based our analyses on a diverse panel of inbred lines crossed with two testers representative of the major heterotic groups in the U.S. (1106 hybrids), as well as a collection of 24 biparental populations crossed with a single tester (1640 hybrids). We investigated three agronomic traits: days to silking (DTS), plant height (PH), and grain yield (GY). Our results point to the presence of dominance for all traits, but also among-locus complementation (epistasis) for DTS and genotype-by-environment interactions for GY. Consistently, dominance improved genomic prediction for PH only. In addition, we assessed enrichment of genetic effects in classes defined by genic regions (gene annotation), structural features (recombination rate and chromatin openness), and evolutionary features (minor allele frequency and evolutionary constraint). We found support for enrichment in genic regions and subsequent improvement of genomic prediction for all traits. Our results suggest that dominance and gene annotations improve genomic prediction across diverse populations in hybrid maize.

Correction to: Soybean iron deficiency chlorosis high throughput phenotyping using an unmanned aircraft system.

[This corrects the article DOI: 10.1186/s13007-019-0478-9.].

Soybean iron deficiency chlorosis high throughput phenotyping using an unmanned aircraft system.

BACKGROUND: Iron deficiency chlorosis (IDC) is an abiotic stress in soybean [Glycine max (L.) Merr.] that causes significant yield reductions. Symptoms of IDC include interveinal chlorosis and stunting of the plant. While there are management practices that can overcome these drastic yield losses, the preferred way to manage IDC is growing tolerant soybean varieties. To develop varieties tolerant to IDC, breeders may easily phenotype up to thousands of candidate soybean lines every year for severity of symptoms related to IDC, a task traditionally done with a 1-5 visual rating scale. The visual rating scale is subjective and, because it is time consuming and laborious, can typically only be accomplished once or twice during a growing season. RESULTS: The goal of this study was to use an unmanned aircraft system (UAS) to improve field screening for tolerance to soybean IDC. During the summer of 2017, 3386 plots were visually scored for IDC stress on two different dates. In addition, images were captured with a DJI Inspire 1 platform equipped with a modified dual camera system which simultaneously captures digital red, green, blue images as well as red, green, near infrared (NIR) images. A pipeline was created for image capture, orthomosaic generation, processing, and analysis. Plant and soil classification was achieved using unsupervised classification resulting in 95% overall classification accuracy. Within the plant classified canopy, the green, yellow, and brown plant pixels were classified and used as features for random forest and neural network models. Overall, the random forest and neural network models achieved similar misclassification rates and classification accuracy, which ranged from 68 to 77% across rating dates. All 36 trials in the field were analyzed using a linear model for both visual score and UAS predicted values on both dates. In 32 of the 36 tests on date 1 and 33 of 36 trials on date 2, the LSD associated with UAS image-based IDC scores was lower than the LSD associated with visual scores, indicating the image-based scores provided more precise measurements of IDC severity. CONCLUSIONS: Overall, the UAS was able to capture differences in IDC stress and may be used for evaluations of candidate breeding lines in a soybean breeding program. This system was both more efficient and precise than traditional scoring methods.

Genome-Wide Association and Gene Co-expression Network Analyses Reveal Complex Genetics of Resistance to Goss's Wilt of Maize.

Goss's bacterial wilt and leaf blight is a disease of maize caused by the gram positive bacterium Clavibacter michiganensis subsp. nebraskensis (Cmn). First discovered in Nebraska, Goss's wilt has now spread to major maize growing states in the United States and three provinces in Canada. Previous studies conducted using elite maize inbred lines and their hybrids have shown that resistance to Goss's wilt is a quantitative trait. The objective of this study was to further our understanding of the genetic basis of resistance to Goss's wilt by using a combined approach of genome-wide association mapping and gene co-expression network analysis. Genome-wide association analysis was accomplished using a diversity panel consisting of 555 maize inbred lines and a set of 450 recombinant inbred lines (RILs) from three bi-parental mapping populations, providing the most comprehensive screening of Goss's wilt resistance to date. Three SNPs in the diversity panel and 10 SNPs in the combined dataset, including the diversity panel and RILs, were found to be significantly associated with Goss's wilt resistance. Each significant SNP explained 1-5% of the phenotypic variation for Goss's wilt (total of 8-11%). To augment the results of genome-wide association mapping and help identify candidate genes, a time course RNA sequencing experiment was conducted using resistant (N551) and susceptible (B14A) maize inbred lines. Gene co-expression network analysis of this time course experiment identified one module of 141 correlated genes that showed differential regulation in response to Cmn inoculations in both resistant and susceptible lines. SNPs inside and flanking these genes explained 13.3% of the phenotypic variation. Among 1,000 random samples of genes, only 8% of samples explained more phenotypic variance for Goss's wilt resistance than those implicated by the co-expression network analysis. While a statistically significant enrichment was not observed (P < 0.05), these results suggest a possible role for these genes in quantitative resistance at the field level and warrant more research on combining gene co-expression network analysis with quantitative genetic analyses to dissect complex disease resistance traits. The results of the GWAS and co-expression analysis both support the complex nature of resistance to this important disease of maize.

Multi-trait Improvement by Predicting Genetic Correlations in Breeding Crosses.

The many quantitative traits of interest to plant breeders are often genetically correlated, which can complicate progress from selection. Improving multiple traits may be enhanced by identifying parent combinations - an important breeding step - that will deliver more favorable genetic correlations (rG ). Modeling the segregation of genomewide markers with estimated effects may be one method of predicting rG in a cross, but this approach remains untested. Our objectives were to: (i) use simulations to assess the accuracy of genomewide predictions of rG and the long-term response to selection when selecting crosses on the basis of such predictions; and (ii) empirically measure the ability to predict genetic correlations using data from a barley (Hordeum vulgare L.) breeding program. Using simulations, we found that the accuracy to predict rG was generally moderate and influenced by trait heritability, population size, and genetic correlation architecture (i.e., pleiotropy or linkage disequilibrium). Among 26 barley breeding populations, the empirical prediction accuracy of rG was low (-0.012) to moderate (0.42), depending on trait complexity. Within a simulated plant breeding program employing indirect selection, choosing crosses based on predicted rG increased multi-trait genetic gain by 11-27% compared to selection on the predicted cross mean. Importantly, when the starting genetic correlation was negative, such cross selection mitigated or prevented an unfavorable response in the trait under indirect selection. Prioritizing crosses based on predicted genetic correlation can be a feasible and effective method of improving unfavorably correlated traits in breeding programs.

Identification and Fine-Mapping of a Soybean Quantitative Trait Locus on Chromosome 5 Conferring Tolerance to Iron Deficiency Chlorosis.

CORE IDEAS: 'Fiskeby III' harbors a combination of abiotic stress traits, including iron deficiency chlorosis (IDC) tolerance. An IDC quantitative trait locus on chromosome Gm05 was identified in genome-wide association studies and biparental populations. Fine-mapping resolved a 137-kb interval containing strong candidate genes. Iron deficiency chlorosis (IDC) is an important nutrient stress for soybean [Glycine max (L.) Merr.] grown in high-pH soils. Despite numerous agronomic attempts to alleviate IDC, genetic tolerance remains the most effective preventative measure against symptoms. In this study, two association mapping populations and a biparental mapping population were used for genetic mapping of IDC tolerance. Quantitative trait loci (QTLs) were identified on chromosomes Gm03, Gm05, and Gm06. Heterogenous inbred families were developed to fine-map the Gm05 QTL, which was uniquely supported in all three mapping populations. Fine-mapping resulted in a QTL with an interval size of 137 kb on the end of the short arm of Gm05, which produced up to a 1.5-point reduction in IDC severity on a 1 to 9 scale in near isogenic lines.

Genome-Wide Association Mapping of Host-Plant Resistance to Soybean Aphid.

Soybean aphid [ Matsumura (Hemiptera: Aphididae)] is the most damaging insect pest of soybean [ (L.) Merr.] in the Upper Midwest of the United States and is primarily controlled by insecticides. Soybean aphid resistance (i.e., genes) has been documented in some soybean accessions but more sources of resistance are needed. Incorporation of the resistance into marketed varieties has also been slow. Genome-wide association mapping can aid in identifying resistant accessions by correlating phenotypic data with single nucleotide polymorphisms (SNPs) across a genome. Aphid population measures from 2366 soybean accessions were collected from published studies screening cultivated soybean () and wild soybean ( Siebold & Zucc.) with aphids exhibiting Biotype 1, 2, or 3 characteristics. Genotypic data were obtained from the SoySNP50K high-density genotyping array previously used to genotype the USDA Soybean Germplasm Collection. Significant associations between SNPs and soybean aphid counts were found on 18 of the 20 soybean chromosomes. Significant SNPs were found on chromosomes 7, 8, 13, and 16 with known genes. SNPs were also significant on chromosomes 1, 2, 4 to 6, 9 to 12, 14, and 17 to 20 where genes have not yet been mapped, suggesting that many genes remain to be discovered. These SNPs can be used to determine accessions that are likely to have novel aphid resistance traits of value for breeding programs.

Genomic Selection in Preliminary Yield Trials in a Winter Wheat Breeding Program.

Genomic prediction (GP) is now routinely performed in crop plants to predict unobserved phenotypes. The use of predicted phenotypes to make selections is an active area of research. Here, we evaluate GP for predicting grain yield and compare genomic and phenotypic selection by tracking lines advanced. We examined four independent nurseries of F3:6 and F3:7 lines trialed at 6 to 10 locations each year. Yield was analyzed using mixed models that accounted for experimental design and spatial variations. Genotype-by-sequencing provided nearly 27,000 high-quality SNPs. Average genomic predictive ability, estimated for each year by randomly masking lines as missing in steps of 10% from 10 to 90%, and using the remaining lines from the same year as well as lines from other years in a training set, ranged from 0.23 to 0.55. The predictive ability estimated for a new year using the other years ranged from 0.17 to 0.28. Further, we tracked lines advanced based on phenotype from each of the four F3:6 nurseries. Lines with both above average genomic estimated breeding value (GEBV) and phenotypic value (BLUP) were retained for more years compared to lines with either above average GEBV or BLUP alone. The number of lines selected for advancement was substantially greater when predictions were made with 50% of the lines from the testing year added to the training set. Hence, evaluation of only 50% of the lines yearly seems possible. This study provides insights to assess and integrate genomic selection in breeding programs of autogamous crops.