Order Matters: Simulations Reveal How to Pair Genomic Prediction with Selection Indices for Better Rice Breeding

Genomic selection has reshaped modern plant breeding, allowing scientists to predict the genetic worth of a seedling from its DNA alone and to make decisions long before a plant ever produces grain. Yet a deceptively simple question has lingered beneath the technology’s rapid adoption: when breeders combine genomic prediction with the classical multi-trait selection indices that have guided crop improvement for nearly a century, does the order in which the two tools are applied change the outcome? A new simulation study published in Theoretical and Applied Genetics suggests that it does, and that the choice of index matters just as much as the choice to use genomics at all.

The research team, led by Roberto Fritsche-Neto of North Carolina State University together with colleagues in Brazil, used stochastic computer simulations to model a rice breeding program over 30 years. Rather than relying on historical field data, which reflects past conditions and cannot easily be re-run under alternative strategies, the researchers built a virtual breeding pipeline in which they could systematically vary the selection method while holding everything else constant. The simulation was anchored in the real-world parameters of the Louisiana State University AgCenter rice breeding program and conducted with the AlphaSimR package, giving the results a practical grounding that purely theoretical exercises often lack.

The simulated program targeted three traits simultaneously: increasing grain yield, reducing the chalkiness rate that degrades rice grain quality, and keeping plant height stable to avoid lodging. Grain yield received the heaviest weight, at 70 percent of the selection objective, while chalkiness and plant height each received 15 percent. The team compared seven selection pipelines built around three classical indices: an empirical index with fixed weights, the Smith–Hazel index that uses genetic and phenotypic covariance matrices, and the Pesek–Baker index, which is constructed from desired gains rather than economic weights. Each index was tested under traditional phenotypic selection and under genomic selection, with the genomic component applied either before index construction or after the index had already been calculated from phenotypes.

The technical distinction between those two architectures turned out to be central. In the first approach, univariate genomic prediction models were fitted separately for each trait using ridge regression best linear unbiased prediction, generating genomic estimated breeding values for yield, chalkiness, and height individually before combining them into an index. In the second, the multi-trait index was first computed from phenotypic data and then treated as a single synthetic trait for genomic prediction. The researchers also refreshed the training population each cycle using a grandparents-parents-offspring strategy, one of the most effective approaches for maintaining prediction accuracy across recurrent cycles, and ran 100 independent replicates of every scenario to capture the stochastic uncertainty inherent in real breeding.

When the 30-year trajectories were compared, no single strategy dominated every trait, but a clear pattern emerged. Applying genomic selection directly to a phenotypically constructed Pesek–Baker index produced the highest gains for grain yield, reaching values close to 10.3 on the standardized scale. When genomic prediction was performed trait by trait before index construction, the Smith–Hazel index delivered the best yield response, near 10.0. For chalkiness, however, the trait-specific prediction approach paired with the Pesek–Baker index was the clear winner, driving the population mean down to roughly negative 4.0, the most favorable reduction observed, compared with only about negative 1.6 under the traditional empirical index.

Plant height told a subtler story. Because the breeding objective was stability rather than directional change, strategies that pushed height upward were actually performing poorly against the stated goal. The Smith–Hazel index combined with genomic selection produced the largest increases in plant height, which ran counter to the program’s intent. The Pesek–Baker index, whether used traditionally or with trait-specific genomic prediction, kept height changes moderate, plateauing near 3.5, while the empirical index was the most conservative of all. Taken together, the results underscore a principle that breeders know intuitively but rarely test so rigorously: the best pipeline is the one that matches the full set of objectives, not the one that maximizes a single headline trait.

The dynamics of selection accuracy and additive genetic variance added a deeper layer to the findings. As expected in recurrent selection, accuracy generally declined over cycles as favorable alleles fixed and the population became more uniform. But the rate and pattern of that decline depended sharply on how genomic information was deployed. When genomic prediction was applied directly to an already aggregated index, accuracy for chalkiness and plant height collapsed dramatically, approaching zero within the first breeding cycle before a partial recovery. The authors attribute this to information compression: folding three traits with different heritabilities, correlations, and selection directions into one value before prediction strips away the trait-specific genetic signals the model needs to learn.

In contrast, estimating genomic values for each trait separately before building the index preserved that information and produced responses more consistent with the multi-trait objectives, with gentler accuracy losses and greater stability after the initial cycles. Additive genetic variance told a parallel story, with genomic scenarios generally depleting variability faster than traditional ones, and the greatest losses appearing in the Smith–Hazel and index-first genomic pipelines. The researchers caution that residual variance alone is not proof of future gain, but the joint pattern of eroding accuracy and shrinking variance in the index-first scenarios points to a genuine methodological limitation rather than simple exhaustion of genetic variation.

The practical implications reach well beyond rice. The study’s authors argue that genomic selection should not be viewed as a wholesale replacement for phenotypic selection but as a tool whose value depends on where it is inserted in the pipeline. For programs juggling traits that must move in opposite directions, such as raising yield while lowering chalkiness, the trait-specific prediction route followed by index construction appears the more reliable architecture. The Pesek–Baker index, built explicitly around desired gains, proved the most faithful vehicle for translating a breeder’s ideotype into a selection criterion, delivering strong yield gains, the deepest chalkiness reduction, and minimal unwanted shifts in plant height. Similar logic applies in crops like malting barley, where quality traits must be held within narrow thresholds rather than pushed to extremes.

The authors acknowledge limitations that shape how far these conclusions can be generalized. The simulations used a fixed marker density of 540 SNPs, and prediction accuracy depends on linkage disequilibrium between markers and causal loci, so different genotyping platforms could shift the magnitude of the observed responses. The genetic architecture was purely additive, with no dominance or epistasis, and only three traits were modeled. Future work, the team notes, should expand the trait set and validate the findings in real breeding populations. Still, for a field racing to feed a growing population under a changing climate, the message is concrete and actionable: in multi-trait breeding, the sequence in which genomic prediction and selection indices are combined is not a technical footnote but a decision that can determine both the magnitude and the direction of genetic progress.

Subject of Research: Integration of genomic prediction with multi-trait selection indices in simulated rice breeding programs

Article Title: Combining genomic prediction and multi-trait indices through stochastic simulations: do index type and deployment order affect genetic gain?

Article References: Fritsche-Neto, R., Queiroz, L. G. C., Viana, J., Gupta, K., Grover, K., & DoVale, J. C. (2026). Combining genomic prediction and multi-trait indices through stochastic simulations: do index type and deployment order affect genetic gain?. Theoretical and Applied Genetics, 139(10), Article 284. https://doi.org/10.1007/s00122-026-05398-0

Image Credits: AI Generated

DOI: 10.1007/s00122-026-05398-0

Keywords: genomic selection, selection index, plant breeding, rice, stochastic simulation, Pesek–Baker index, Smith–Hazel index, genetic gain, grain yield, chalkiness, additive genetic variance, RR-BLUP