Leveraging spatiotemporal genomic breeding value estimates of dry matter yield and herbage quality in ryegrass via random regression models

Joint modeling of correlated multienvironment and multiharvest data of perennial crop species may offer advantages in prediction schemes and a better understanding of the underlying dynamics in space and time. The goal of the present study was to investigate the relevance of incorporating the longitudinal dimension of within-season multiple measurements of forage perennial ryegrass (Lolium perenne L.) traits in a reaction-norm model setup that additionally accounts for genotype × environment (G × E) interactions. Genetic parameters and accuracy of genomic estimated breeding value (gEBV) predictions were investigated by fitting three genomic random regression models (gRRMs) using Legendre polynomial functions to the data. Genomic DNA sequencing of family pools of diploid perennial ryegrass was performed using DNA nanoball-based technology and yielded 56,645 single-nucleotide polymorphisms, which were used to calculate the allele frequency-based genomic relationship matrix. Biomass yield's estimated additive genetic variance and heritability values were higher in later harvests. The additive genetic correlations were moderate to low in early measurements and peaked at intermediates with fairly stable values across the environmental gradient except for the initial harvest data collection. This led to the conclusion that complex (G × E) arises from spatial and temporal dimensions in the early season with lower reranking trends thereafter. In general, modeling the temporal dimension with a second-order orthogonal polynomial improved the accuracy of gEBV prediction for nutritive quality traits, but no gain in prediction accuracy was detected for dry matter yield (DMY). This study leverages the flexibility and usefulness of gRRM models for perennial ryegrass breeding and can be readily extended to other multiharvest crops.