Mapping and Predicting Non-Linear Brassica rapa Growth Phenotypes Based on Bayesian and Frequentist Complex Trait Estimation

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.

Rapid Chlorophyll a Fluorescence Light Response Curves Mechanistically Inform Photosynthesis Modeling

Crop improvement is crucial to ensuring global food security under climate change, and hence there is a pressing need for phenotypic observations that are both high throughput and improve mechanistic understanding of plant responses to environmental cues and limitations. In this study, chlorophyll a fluorescence light response curves and gas-exchange observations are combined to test the photosynthetic response to moderate drought in four genotypes of Brassica rapa The quantum yield of PSII (ϕ PSII ) is here analyzed as an exponential decline under changing light intensity and soil moisture. Both the maximum ϕ PSII and the rate of ϕ PSII decline across a large range of light intensities (0-1,000 μmol photons m-2 s-1; β PSII ) are negatively affected by drought. We introduce an alternative photosynthesis model (β PSII model) incorporating parameters from rapid fluorescence response curves. Specifically, the model uses β PSII as an input for estimating the photosynthetic electron transport rate, which agrees well with two existing photosynthesis models (Farquhar-von Caemmerer-Berry and Yin). The β PSII model represents a major improvement in photosynthesis modeling through the integration of high-throughput fluorescence phenotyping data, resulting in gained parameters of high mechanistic value.

Variation within laminae: Semi-automated methods for quantifying leaf venation using phenoVein

PREMISE

Physiological processes may vary within leaf laminae; however, the accompanying heterogeneity in leaf venation is rarely investigated because its quantification can be time consuming. Here we introduce accelerated protocols using existing software to increase sample throughput and ask whether laminae venation varies among three crop types and four subspecies of Brassica rapa.

METHODS

FAA (formaldehyde, glacial acetic acid, and ethanol)-fixed samples were stored in ethanol. Without performing any additional clearing or staining, we tested two methods of image acquisition at three locations along the proximal-distal axis of the laminae and estimated the patterns of venation using the program phenoVein. We developed and made available an R script to handle the phenoVein output and then analyzed our data using linear mixed-effects models.

RESULTS

Beyond fixation and storage, staining and clearing are not necessary to estimate leaf venation using phenoVein if the images are acquired using a stereomicroscope. All estimates of venation required some manual adjustment. We found a significant effect of location within the laminae for all aspects of venation.

DISCUSSION

By removing the clearing and staining steps and utilizing the semi-automated program phenoVein, we quickly and cheaply acquired leaf venation data. Venation may be an important target for crop breeding efforts, particularly if intralaminar variation correlates with variation in physiological processes, which remains an open question.

Integrating transcriptomic network reconstruction and eQTL analyses reveals mechanistic connections between genomic architecture and Brassica rapa development

Plant developmental dynamics can be heritable, genetically correlated with fitness and yield, and undergo selection. Therefore, characterizing the mechanistic connections between the genetic architecture governing plant development and the resulting ontogenetic dynamics of plants in field settings is critically important for agricultural production and evolutionary ecology. We use hierarchical Bayesian Function-Valued Trait (FVT) models to estimate Brassica rapa growth curves throughout ontogeny, across two treatments, and in two growing seasons. We find genetic variation for plasticity of growth rates and final sizes, but not the inflection point (transition from accelerating to decelerating growth) of growth curves. There are trade-offs between growth rate and duration, indicating that selection for maximum yields at early harvest dates may come at the expense of late harvest yields and vice versa. We generate eigengene modules and determine which are co-expressed with FVT traits using a Weighted Gene Co-expression Analysis. Independently, we seed a Mutual Rank co-expression network model with FVT traits to identify specific genes and gene networks related to FVT. GO-analyses of eigengene modules indicate roles for actin/cytoskeletal genes, herbivore resistance/wounding responses, and cell division, while MR networks demonstrate a close association between metabolic regulation and plant growth. We determine that combining FVT Quantitative Trait Loci (QTL) and MR genes/WGCNA eigengene expression profiles better characterizes phenotypic variation than any single data type (i.e. QTL, gene, or eigengene alone). Our network analysis allows us to employ a targeted eQTL analysis, which we use to identify regulatory hotspots for FVT. We examine cis vs. trans eQTL that mechanistically link FVT QTL with structural trait variation. Colocalization of FVT, gene, and eigengene eQTL provide strong evidence for candidate genes influencing plant height. The study is the first to explore eQTL for FVT, and specifically do so in agroecologically relevant field settings.

We estimate the developmental dynamics of plant growth using mathematical functions to fit continuous functions to discrete plant height data collected throughout growth, and we use the parameters defining these mathematical functions as data. We identify genomic regions controlling plant growth and filter a novel transcriptomic data set using network reconstruction models to identify the genes and eigengenes associated with plant height. We combine these genomic and transcriptomic data to predict variation in plant height, and we use quantitative genetics to mechanistically connect plant genetics, transcriptomics, and development. Our approach demonstrates two powerful methods for the type of data reduction (FVT modeling and gene expression network reconstruction for targeted eQTL analyses) and data integration that will be necessary for driving forward the field of genetics in the post-genomic era. To the best of our knowledge, we are the first to apply these techniques to continuous models of plant development, and the first to do so in agroecologically relevant field settings.