bioRxivGeneticsPreprint
Software update: Moving the R package sommer to multivariate mixed models for genome-assisted prediction
Abstract
In the last decade the use of mixed models has become a pivotal part in the implementation of genome-assisted prediction in plant and animal breeding programs. Exploiting the use genetic correlation among traits through multivariate predictions has been proposed in recent years as a way to boost prediction accuracy and understand pleiotropy and other genetic and ecological phenomena better. Multiple mixed model solvers able to use relationship matrices or deal with marker-based incidence matrices have been released in the last years but multivariate versions are scarse. Such solvers have become quite popular in plant and animal breeding thanks to user-friendly platforms such as R. Among such software one of the most recent and popular is the sommer package. In this short communication we discuss the update of the package that is able to run multivariate mixed models with multiple random effects and different covariance structures at the level of random effects and trait-to-trait covariance along with other functionalities for genetic analysis and field trial analysis to enhance the genome-assisted prediction capabilities of researchers.
- organism: plant (other)
- organism: non-human primate
- assay: whole-genome sequencing
- research focus: variant effect prediction
- research focus: genome assembly
- research focus: GWAS / association mapping
- research focus: fine-mapping
- approach: statistical model
- approach: mixed model
Abstract summary
Claude Opus 5.5 · xhigh effort · 18 Oct 2026 · from the abstract
The sommer R package gains a multivariate solver: several random effects at once, with covariance structures set per effect and between traits, for genome-assisted breeding prediction.
Genome-assisted prediction in plant and animal breeding rests on mixed models, and borrowing strength across genetically correlated traits has been put forward as a route to better prediction accuracy and to studying pleiotropy. Solvers that take relationship matrices or marker-based incidence matrices exist, but the author describes multivariate ones as scarce, which this update aims to fill inside an R package already in use by breeders.
Key results
ClaimedThe updated solver fits multivariate mixed models in which several random effects can be included and covariance structures can be set both within random effects and for the trait-to-trait covariance.
Described as the content of the software update; the abstract gives no prediction accuracies, no comparison with univariate fits or other solvers, and no runtime or memory figures. · Abstract
ClaimedFurther routines for genetic analysis and for field trial analysis are bundled with the multivariate solver.
Listed as added functionality supporting genome-assisted prediction; the individual routines are not named or evaluated in the abstract. · Abstract
Methods, data and reproducibility
A short communication on an update to the sommer R package, extending its mixed model solver from univariate to multivariate fits with multiple random effects, covariance structures specified at the level of random effects and across traits, together with added genetic and field trial analysis functionality.
- Methods & tools
- sommer, R
Limitations (2)
- Nothing is benchmarked: no dataset, no accuracy comparison against univariate models or competing solvers, and no runtime or memory measurements accompany the update. · noted in this summary
- The claim that multivariate solvers are scarce is made without naming the alternatives, so the gap this update fills cannot be assessed from the text. · noted in this summary
Read this if you run genomic prediction for breeding in R and need multivariate mixed models with covariance structures set across traits and random effects.
- Citations
- 52
Reproducibility ledger
Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.
Sign in to file a reproducibility report.
Related papers
- MegaLMM: Mega-scale linear mixed models for genomic predictions with thousands of traits
Daniel E. Runcie, Jiayi Qu, … Lorin A. Crawford·Genome biologyPeer-reviewed·Jul 23, 2021
- lme4breeding: a fast linear mixed model for multi-trait/multi-environment experiments
Giovanny E. Covarrubias-Pazaran·bioRxiv·May 17, 2024
- Pitfalls and Remedies for Cross Validation with Multi-trait Genomic Prediction Methods
Daniel E. Runcie, Hao Cheng·G3 Genes Genomes GeneticsPeer-reviewed·Sep 11, 2019
- Accuracies Of Univariate And Multivariate Genomic Prediction Models In African Cassava.
Uche Godfrey Okeke, Deniz Akdemir, … Jean‐Luc Jannink·Genetics Selection EvolutionPeer-reviewed·Dec 1, 2017
- Improving Multi-Trait Genomic Prediction Efficiency Through The Incorporation Of Synthetic Traits Selected Based on Co-heritability
Ashmita Upadhyay, Ruhana Azam, … Samuel B. Fernandes·bioRxiv·Jun 23, 2025
- Multi-Trait Improvement by Predicting Genetic Correlations in Breeding Crosses
Jeffrey Neyhart, Aaron Joel Lorenz, Kevin P. Smith·G3 Genes Genomes GeneticsPeer-reviewed·Jul 29, 2019