Connecting mathematical models to genomes: joint estimation of model parameters and genome-wide marker effects on these parameters

Modeling soybean growth: A mixed model approach

The evaluation of plant and animal growth, separately for genetic and environmental effects, is necessary for genetic understanding and genetic improvement of environmental responses of plants and animals. We propose to extend an existing approach that combines nonlinear mixed-effects model (NLMEM) and the stochastic approximation of the Expectation-Maximization algorithm (SAEM) to analyze genetic and environmental effects on plant growth. These tools are widely used in many fields but very rarely in plant biology. During model formulation, a nonlinear function describes the shape of growth, and random effects describe genetic and environmental effects and their variability. Genetic relationships among the varieties were also integrated into the model using a genetic relationship matrix. The SAEM algorithm was chosen as an efficient alternative to MCMC methods, which are more commonly used in the domain. It was implemented to infer the expected growth patterns in the analyzed population and the expected curves for each variety through a maximum-likelihood and a maximum-a-posteriori approaches, respectively. The obtained estimates can be used to predict the growth curves for each variety. We illustrate the strengths of the proposed approach using simulated data and soybean plant growth data obtained from a soybean cultivation experiment conducted at the Arid Land Research Center, Tottori University. In this experiment, plant height was measured daily using drones, and the growth was monitored for approximately 200 soybean cultivars for which whole-genome sequence data were available. The NLMEM approach improved our understanding of the determinants of soybean growth and can be successfully used for the genomic prediction of growth pattern characteristics.

Integration of Crop Growth Models and Genomic Prediction

Crop growth models (CGMs) consist of multiple equations that represent physiological processes of plants and simulate crop growth dynamically given environmental inputs. Because parameters of CGMs are often genotype-specific, gene effects can be related to environmental inputs through CGMs. Thus, CGMs are attractive tools for predicting genotype by environment (G×E) interactions. This chapter reviews CGMs, genetic analyses using these models, and the status of studies that integrate genomic prediction with CGMs. Examples of CGM analyses are also provided.

Adoption and Optimization of Genomic Selection To Sustain Breeding for Apricot Fruit Quality

Genomic selection (GS) is a breeding approach which exploits genome-wide information and whose unprecedented success has shaped several animal and plant breeding schemes through delivering their genetic progress. This is the first study assessing the potential of GS in apricot (Prunus armeniaca) to enhance postharvest fruit quality attributes. Genomic predictions were based on a F1 pseudo-testcross population, comprising 153 individuals with contrasting fruit quality traits. They were phenotyped for physical and biochemical fruit metrics in contrasting climatic conditions over two years. Prediction accuracy (PA) varied from 0.31 for glucose content with the Bayesian LASSO (BL) to 0.78 for ethylene production with RR-BLUP, which yielded the most accurate predictions in comparison to Bayesian models and only 10% out of 61,030 SNPs were sufficient to reach accurate predictions. Useful insights were provided on the genetic architecture of apricot fruit quality whose integration in prediction models improved their performance, notably for traits governed by major QTL. Furthermore, multivariate modeling yielded promising outcomes in terms of PA within training partitions partially phenotyped for target traits. This provides a useful framework for the implementation of indirect selection based on easy-to-measure traits. Thus, we highlighted the main levers to take into account for the implementation of GS for fruit quality in apricot, but also to improve the genetic gain in perennial species.

Leveraging genome-enabled growth models to study shoot growth responses to water deficit in rice

Elucidating genotype-by-environment interactions and partitioning its contribution to phenotypic variation remains a challenge for plant scientists. We propose a framework that utilizes genome-wide markers to model genotype-specific shoot growth trajectories as a function of time and soil water availability. A rice diversity panel was phenotyped daily for 21 d using an automated, high-throughput image-based, phenotyping platform that enabled estimation of daily shoot biomass and soil water content. Using these data, we modeled shoot growth as a function of time and soil water content, and were able to determine the time point where an inflection in the growth trajectory occurred. We found that larger, more vigorous plants exhibited an earlier repression in growth compared with smaller, slow-growing plants, indicating a trade-off between early vigor and tolerance to prolonged water deficits. Genomic inference for model parameters and time of inflection (TOI) identified several candidate genes. This study is the first to utilize a genome-enabled growth model to study drought responses in rice, and presents a new approach to jointly model dynamic morpho-physiological responses and environmental covariates.