{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,9,29]],"date-time":"2026-09-29T10:00:42Z","timestamp":1790676042317,"version":"4.1.0"},"reference-count":41,"publisher":"Frontiers Media SA","license":[{"start":{"date-parts":[[2026,9,29]],"date-time":"2026-09-29T00:00:00Z","timestamp":1790640000000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":["frontiersin.org"],"crossmark-restriction":true},"short-container-title":["Front. Genet."],"abstract":"<jats:p>Accurate modeling of genotype-by-environment (G \u00d7 E) interaction is a central challenge in plant breeding, as it determines how genetic potential is expressed across heterogeneous and changing environments. In genomic prediction, a widely adopted strategy represents G \u00d7 E interaction through the Hadamard (element-wise) product of genotype and environment covariance matrices, capturing interaction effects via joint similarity. While effective, this formulation may not fully represent the structured complexity of G \u00d7 E relationships observed in multi-environment trials. In this study, we extend a recently proposed matrix-product framework for covariance construction to explicitly model G \u00d7 E interaction. Unlike the Hadamard approach, which captures similarity-by-similarity effects, the proposed method derives symmetric kernels from the matrix multiplication of genotype and environment kernels. These kernels retain complementary subsets of matrix-product information, but the numerical upper\/lower assignment is tied to record ordering and should not be interpreted as an intrinsic biological direction. Because the direct matrix product and its symmetrized subcomponents are not guaranteed to be valid covariance matrices, positive semidefiniteness must be examined explicitly before the kernels are used as random-effect covariance structures. The proposed framework was evaluated using eleven multi-environment datasets in wheat and rice, comprising diverse traits and experimental structures. On average across datasets, models containing matrix-product-derived kernels showed higher predictive performance than the conventional Hadamard-only formulation, with an average COR difference of up to 13.2% for the combined model. These empirical gains should be interpreted as average predictive differences under the validation scheme used here rather than as a claim of universal or statistically significant superiority. These results show that matrix-product-derived terms can carry predictive information complementary to the conventional Hadamard interaction term. At the same time, the present revision makes explicit two structural limitations that must be considered when interpreting these kernels: positional triangular extraction is record-order dependent, and balanced incidence combined with centered genomic relationships can make the matrix product degenerate or nearly zero. These findings define both the practical scope of the current formulation and priorities for further methodological development.<\/jats:p>","DOI":"10.3389\/fgene.2026.1961511","type":"journal-article","created":{"date-parts":[[2026,9,29]],"date-time":"2026-09-29T09:39:11Z","timestamp":1790674751000},"update-policy":"https:\/\/doi.org\/10.3389\/crossmark-policy","source":"Crossref","is-referenced-by-count":0,"title":["A novel matrix multiplication framework for modeling genotype-by-environment interaction in genomic prediction"],"prefix":"10.3389","volume":"17","author":[{"given":"Osval A.","family":"Montesinos-L\u00f3pez","sequence":"first","affiliation":[{"name":"Facultad de Telem\u00e1tica, Universidad de Colima","place":["Colima, Mexico"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Abelardo","family":"Montesinos-L\u00f3pez","sequence":"additional","affiliation":[{"name":"Centro Universitario de Ciencias Exactas e Ingenier\u00edas (CUCEI), Universidad de Guadalajara","place":["Guadalajara, Mexico"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Cricelio","family":"Montesinos-Lopez","sequence":"additional","affiliation":[{"name":"Consultant","place":["Mexico, Mexico"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Jos\u00e9","family":"Crossa","sequence":"additional","affiliation":[{"name":"Centro de Estad\u00edstica y Ciencias de Datos, Colegio de Postgraduados (COLPOS)","place":["Montecillo, Mexico"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Jaime","family":"Cuevas","sequence":"additional","affiliation":[{"name":"Divisi\u00f3n de Ciencias, Ingenier\u00eda y Tecnolog\u00edas (DCIT), Universidad Aut\u00f3noma del Estado de Quintana Roo","place":["Chetumal, Mexico"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Susanne","family":"Dreisigacker","sequence":"additional","affiliation":[{"name":"International Maize and Wheat Improvement Center (CIMMYT)","place":["Texcoco, Mexico"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Carlos M.","family":"Hern\u00e1ndez Su\u00e1rez","sequence":"additional","affiliation":[{"name":"Coordinaci\u00f3n General de Investigaci\u00f3n Cient\u00edfica, Universidad de Colima","place":["Colima, Mexico"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Rodomiro","family":"Ortiz","sequence":"additional","affiliation":[{"name":"Department of Plant Breeding, Swedish University of Agricultural Sciences (SLU)","place":["Lomma, Sweden"]}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"1965","published-online":{"date-parts":[[2026,9,29]]},"reference":[{"key":"B1","doi-asserted-by":"publisher","first-page":"552","DOI":"10.1016\/j.molp.2024.03.007","article-title":"Genomic selection in plant breeding: key factors shaping two decades of progress","volume":"17","author":"Alemu","year":"2024","journal-title":"Mol. Plant"},{"key":"B2","doi-asserted-by":"publisher","first-page":"2633","DOI":"10.1093\/bioinformatics\/btm308","article-title":"TASSEL: software for association mapping of complex traits in diverse samples","volume":"23","author":"Bradbury","year":"2007","journal-title":"Bioinformatics"},{"key":"B3","doi-asserted-by":"publisher","first-page":"707","DOI":"10.2135\/cropsci2011.06.0299","article-title":"Genomic prediction of breeding values when modeling genotype \u00d7 environment interaction using pedigree and dense molecular markers","volume":"52","author":"Burgue\u00f1o","year":"2012","journal-title":"Crop Sci."},{"key":"B4","first-page":"33","article-title":"Modelling plant breeding programs","volume":"2","author":"Cooper","year":"1999","journal-title":"Trends Agron."},{"key":"B5","doi-asserted-by":"publisher","first-page":"2141","DOI":"10.2135\/cropsci2015.08.0512","article-title":"Use of crop growth models with whole\u2010genome prediction: application to a maize multienvironment trial","volume":"56","author":"Cooper","year":"2016","journal-title":"Crop Sci."},{"key":"B6","doi-asserted-by":"publisher","first-page":"998","DOI":"10.2135\/cropsci1999.0011183X003900040012x","article-title":"Prediction assessment of shrinkage estimators of multiplicative models for multi-environment cultivar trials","volume":"39","author":"Cornelius","year":"1999","journal-title":"Crop Sci."},{"key":"B7","doi-asserted-by":"publisher","first-page":"193","DOI":"10.2135\/cropsci1992.0011183X003200010043x","article-title":"Statistical tests and estimators of multiplicative models for genotype-by-environment interaction","volume":"32","author":"Cornelius","year":"1992","journal-title":"Crop Sci."},{"key":"B8","doi-asserted-by":"publisher","first-page":"92","DOI":"10.1038\/s41437-020-00353-1","article-title":"Nonlinear kernels, dominance, and envirotyping data increase the accuracy of genome-based prediction in multi-environment trials","volume":"126","author":"Costa-Neto","year":"2021","journal-title":"Heredity"},{"key":"B9","doi-asserted-by":"publisher","first-page":"55","DOI":"10.1016\/S0065-2113(08)60818-4","article-title":"Statistical analyses of multilocation trials","volume":"44","author":"Crossa","year":"1990","journal-title":"Adv. Agron."},{"key":"B10","doi-asserted-by":"publisher","first-page":"713","DOI":"10.1534\/genetics.110.118521","article-title":"Prediction of genetic values of quantitative traits in plant breeding using pedigree and molecular markers","volume":"186","author":"Crossa","year":"2010","journal-title":"Genetics"},{"key":"B11","doi-asserted-by":"publisher","first-page":"961","DOI":"10.1016\/j.tplants.2017.08.011","article-title":"Genomic selection in plant breeding: methods, models, and perspectives","volume":"22","author":"Crossa","year":"2017","journal-title":"Trends Plant Sci."},{"key":"B12","doi-asserted-by":"publisher","first-page":"651480","DOI":"10.3389\/fpls.2021.651480","article-title":"The modern plant breeding triangle: optimizing the use of genomics, phenomics, and enviromics data","volume":"12","author":"Crossa","year":"2021","journal-title":"Front. Plant Sci."},{"key":"B13","doi-asserted-by":"publisher","first-page":"41","DOI":"10.1534\/g3.116.035584","article-title":"Bayesian genomic prediction with genotype \u00d7 environment interaction kernel models","volume":"7","author":"Cuevas","year":"2017","journal-title":"G3 Genes|Genomes|Genetics"},{"key":"B14","doi-asserted-by":"publisher","first-page":"1605202","DOI":"10.3389\/fpls.2025.1605202","article-title":"Enhancing wheat genomic prediction by a hybrid kernel approach","volume":"16","author":"Cuevas","year":"2025","journal-title":"Front. Plant Sci."},{"key":"B15","doi-asserted-by":"publisher","first-page":"36","DOI":"10.2135\/cropsci1966.0011183X000600010011x","article-title":"Stability parameters for comparing varieties","volume":"6","author":"Eberhart","year":"1966","journal-title":"Crop Sci."},{"key":"B16","doi-asserted-by":"publisher","first-page":"2090","DOI":"10.2135\/cropsci2015.01.0061","article-title":"Half a century of studying genotype \u00d7 environment interactions in plant breeding experiments","volume":"56","author":"Elias","year":"2016","journal-title":"Crop Sci."},{"key":"B17","doi-asserted-by":"publisher","first-page":"e19379","DOI":"10.1371\/journal.pone.0019379","article-title":"A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species","volume":"6","author":"Elshire","year":"2011","journal-title":"PLoS ONE"},{"key":"B18","doi-asserted-by":"publisher","first-page":"250","DOI":"10.3835\/plantgenome2011.08.0024","article-title":"Ridge regression and other kernels for genomic selection with R package rrBLUP","volume":"4","author":"Endelman","year":"2011","journal-title":"Plant Genome"},{"key":"B19","doi-asserted-by":"publisher","first-page":"742","DOI":"10.1071\/AR9630742","article-title":"The analysis of adaptation in a plant-breeding programme","volume":"14","author":"Finlay","year":"1963","journal-title":"Aust. J. Agric. Res."},{"key":"B20","volume-title":"Statistical Methods for Research Workers","author":"Fisher","year":"1925"},{"key":"B21","doi-asserted-by":"publisher","first-page":"705","DOI":"10.2307\/2531585","article-title":"Model selection and validation for yield trials with interaction","volume":"44","author":"Gauch","year":"1988","journal-title":"Biometrics"},{"key":"B22","volume-title":"Applications of Linear Models in Animal Breeding","author":"Henderson","year":"1984"},{"key":"B23","doi-asserted-by":"publisher","first-page":"463","DOI":"10.1007\/s00122-013-2231-5","article-title":"Integrating environmental covariates and crop modeling into the genomic selection framework to predict genotype by environment interactions","volume":"127","author":"Heslot","year":"2014","journal-title":"Theor. Appl. Genet."},{"key":"B24","doi-asserted-by":"publisher","first-page":"e20034","DOI":"10.1002\/tpg2.20034","article-title":"Genome\u2010based prediction of multiple wheat quality traits in multiple years","volume":"13","author":"Ibba","year":"2020","journal-title":"Plant Genome"},{"key":"B25","doi-asserted-by":"publisher","first-page":"595","DOI":"10.1007\/s00122-013-2243-1","article-title":"A reaction norm model for genomic selection using high-dimensional genomic and environmental data","volume":"127","author":"Jarqu\u00edn","year":"2014","journal-title":"Theor. Appl. Genet."},{"key":"B26","doi-asserted-by":"publisher","first-page":"180017","DOI":"10.3835\/plantgenome2018.03.0017","article-title":"Prospects and challenges of applied genomic selection\u2014A new paradigm in breeding for grain yield in bread wheat","volume":"11","author":"Juliana","year":"2018","journal-title":"Plant Genome"},{"key":"B27","doi-asserted-by":"publisher","first-page":"151","DOI":"10.1016\/j.eja.2018.01.007","article-title":"Leveraging biological insight and environmental variation to improve phenotypic prediction: integrating crop growth models (CGM) with whole genome prediction (WGP)","volume":"100","author":"Messina","year":"2018","journal-title":"European J. Agron."},{"key":"B28","doi-asserted-by":"publisher","first-page":"1819","DOI":"10.1093\/genetics\/157.4.1819","article-title":"Prediction of total genetic value using genome-wide dense marker maps","volume":"157","author":"Meuwissen","year":"2001","journal-title":"Genetics"},{"key":"B29","doi-asserted-by":"publisher","first-page":"2383","DOI":"10.1534\/g3.115.021667","article-title":"LinkImpute: fast and accurate genotype imputation for nonmodel organisms","volume":"5","author":"Money","year":"2015","journal-title":"G3 Genes|Genomes|Genetics"},{"key":"B30","doi-asserted-by":"publisher","first-page":"1165","DOI":"10.1534\/g3.116.028118","article-title":"Genomic Bayesian prediction model for count data with genotype \u00d7 environment interaction","volume":"6","author":"Montesinos-L\u00f3pez","year":"2016","journal-title":"G3 Genes|Genomes|Genetics"},{"key":"B31","doi-asserted-by":"publisher","first-page":"131","DOI":"10.1534\/g3.117.300309","article-title":"Prediction of multiple-trait and multiple-environment genomic data using recommender systems","volume":"8","author":"Montesinos-L\u00f3pez","year":"2018","journal-title":"G3 Genes|Genomes|Genetics"},{"key":"B32","doi-asserted-by":"publisher","first-page":"19","DOI":"10.1186\/s12864-020-07319-x","article-title":"A review of deep learning applications for genomic selection","volume":"22","author":"Montesinos-L\u00f3pez","year":"2021","journal-title":"BMC Genomics"},{"key":"B33","doi-asserted-by":"publisher","first-page":"1636438","DOI":"10.3389\/fgene.2025.1636438","article-title":"Genomic prediction powered by multi-omics data","volume":"16","author":"Montesinos-L\u00f3pez","year":"","journal-title":"Front. Genet."},{"key":"B34","doi-asserted-by":"publisher","first-page":"e70171","DOI":"10.1002\/tpg2.70171","article-title":"Hybrid kernels integrating genomic and multispectral data improve wheat genomic prediction accuracy","volume":"19","author":"Montesinos\u2010L\u00f3pez","year":"","journal-title":"Plant Genome"},{"key":"B35","doi-asserted-by":"publisher","first-page":"1519","DOI":"10.1534\/g3.119.400064","article-title":"Integrating molecular markers and environmental covariates to interpret genotype by environment interaction in rice (Oryza sativa L.) grown in subtropical areas","volume":"9","author":"Monteverde","year":"2019","journal-title":"G3 Genes|Genomes|Genetics"},{"key":"B36","doi-asserted-by":"publisher","first-page":"483","DOI":"10.1534\/genetics.114.164442","article-title":"Genome-wide regression and prediction with the BGLR statistical package","volume":"198","author":"P\u00e9rez","year":"2014","journal-title":"Genetics"},{"key":"B37","doi-asserted-by":"publisher","first-page":"e32253","DOI":"10.1371\/journal.pone.0032253","article-title":"Development of high-density genetic maps for barley and wheat using a novel two-enzyme genotyping-by-sequencing approach","volume":"7","author":"Poland","year":"2012","journal-title":"PLoS ONE"},{"key":"B38","doi-asserted-by":"publisher","first-page":"4414","DOI":"10.3168\/jds.2007-0980","article-title":"Efficient methods to compute genomic predictions","volume":"91","author":"VanRaden","year":"2008","journal-title":"J. Dairy Sci."},{"key":"B39","doi-asserted-by":"publisher","first-page":"1783939","DOI":"10.3389\/fgene.2026.1783939","article-title":"Beyond QTL and GWAS: how deep learning, graph models, and multi-omics are reshaping plant genomic prediction analysis","volume":"17","author":"Wang","year":"2026","journal-title":"Front. Genet."},{"key":"B40","doi-asserted-by":"publisher","first-page":"622","DOI":"10.1002\/csc2.20052","article-title":"Predictive breeding for maize: making use of molecular phenotypes, machine learning, and physiological crop models","volume":"60","author":"Washburn","year":"2020","journal-title":"Crop Sci."},{"key":"B41","doi-asserted-by":"publisher","first-page":"597","DOI":"10.2135\/cropsci2000.403597x","article-title":"Cultivar evaluation and mega\u2010environment investigation based on the GGE biplot","volume":"40","author":"Yan","year":"2000","journal-title":"Crop Sci."}],"container-title":["Frontiers in Genetics"],"original-title":[],"link":[{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fgene.2026.1961511\/full","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,9,29]],"date-time":"2026-09-29T09:39:13Z","timestamp":1790674753000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fgene.2026.1961511\/full"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,9,29]]},"references-count":41,"alternative-id":["10.3389\/fgene.2026.1961511"],"URL":"https:\/\/doi.org\/10.3389\/fgene.2026.1961511","relation":{},"ISSN":["1664-8021"],"issn-type":[{"value":"1664-8021","type":"electronic"}],"subject":[],"published":{"date-parts":[[2026,9,29]]},"article-number":"1961511"}}