High-throughput computer vision introduces the time axis to a quantitative trait map of a plant growth response

QTL for the Kinematic Traits That Define the Arabidopsis Root Elongation Zone and Their Relationship to Gravitropism

Cell expansion in a discrete region called the elongation zone drives root elongation. Analyzing time lapse images can quantify the expansion in kinematic terms as if it were fluid flow. We used horizontal microscopes to collect images from which custom software extracted the length of the elongation zone, the peak relative elemental growth rate (REGR) within it, the axial position of the REGR peak, and the root elongation rate. Automation enabled these kinematic traits to be measured in 1575 Arabidopsis seedlings representing 162 recombinant inbred lines (RILs) derived from a cross of Cvi and Ler ecotypes. We mapped ten quantitative trait loci (QTL), affecting the four kinematic traits. Three QTL affected two or more traits in these vertically oriented seedlings. We compared this genetic architecture with that previously determined for gravitropism using the same RIL population. The major QTL peaks for the kinematic traits did not overlap with the gravitropism QTL. Furthermore, no single kinematic trait correlated with quantitative descriptors of the gravitropism response curve across this population. In addition to mapping QTL for growth zone traits, this study showed that the size and shape of the elongation zone may vary widely without affecting the differential growth induced by gravity.

OPEN leaf: an open-source cloud-based phenotyping system for tracking dynamic changes at leaf-specific resolution in Arabidopsis

The first draft of the Arabidopsis genome was released more than 20 years ago and despite intensive molecular research, more than 30% of Arabidopsis genes remained uncharacterized or without an assigned function. This is in part due to gene redundancy within gene families or the essential nature of genes, where their deletion results in lethality (i.e., the dark genome). High-throughput plant phenotyping (HTPP) offers an automated and unbiased approach to characterize subtle or transient phenotypes resulting from gene redundancy or inducible gene silencing; however, access to commercial HTPP platforms remains limited. Here we describe the design and implementation of OPEN leaf, an open-source phenotyping system with cloud connectivity and remote bilateral communication to facilitate data collection, sharing and processing. OPEN leaf, coupled with our SMART imaging processing pipeline was able to consistently document and quantify dynamic changes at the whole rosette level and leaf-specific resolution when plants experienced changes in nutrient availability. Our data also demonstrate that VIS sensors remain underutilized and can be used in high-throughput screens to identify and characterize previously unidentified phenotypes in a leaf-specific time-dependent manner. Moreover, the modular and open-source design of OPEN leaf allows seamless integration of additional sensors based on users and experimental needs.

Dynamic growth QTL action in diverse light environments: characterization of light regime-specific and stable QTL in Arabidopsis

Plant growth is a complex process affected by a multitude of genetic and environmental factors and their interactions. To identify genetic factors influencing plant performance under different environmental conditions, vegetative growth was assessed in Arabidopsis thaliana cultivated under constant or fluctuating light intensities, using high-throughput phenotyping and genome-wide association studies. Daily automated non-invasive phenotyping of a collection of 382 Arabidopsis accessions provided growth data during developmental progression under different light regimes at high temporal resolution. Quantitative trait loci (QTL) for projected leaf area, relative growth rate, and PSII operating efficiency detected under the two light regimes were predominantly condition-specific and displayed distinct temporal activity patterns, with active phases ranging from 2 d to 9 d. Eighteen protein-coding genes and one miRNA gene were identified as potential candidate genes at 10 QTL regions consistently found under both light regimes. Expression patterns of three candidate genes affecting projected leaf area were analysed in time-series experiments in accessions with contrasting vegetative leaf growth. These observations highlight the importance of considering both environmental and temporal patterns of QTL/allele actions and emphasize the need for detailed time-resolved analyses under diverse well-defined environmental conditions to effectively unravel the complex and stage-specific contributions of genes affecting plant growth processes.

Leveraging orthology within maize and Arabidopsis QTL to identify genes affecting natural variation in gravitropism

Plants typically orient their organs with respect to the Earth's gravity field by a dynamic process called gravitropism. To discover conserved genetic elements affecting seedling root gravitropism, we measured the process in a set of Zea mays (maize) recombinant inbred lines with machine vision and compared the results with those obtained in a similar study of Arabidopsis thaliana. Each of the several quantitative trait loci that we mapped in both species spanned many hundreds of genes, too many to test individually for causality. We reasoned that orthologous genes may be responsible for natural variation in monocot and dicot root gravitropism. If so, pairs of orthologous genes affecting gravitropism may be present within the maize and Arabidopsis QTL intervals. A reciprocal comparison of sequences within the QTL intervals identified seven pairs of such one-to-one orthologs. Analysis of knockout mutants demonstrated a role in gravitropism for four of the seven: CCT2 functions in phosphatidylcholine biosynthesis, ATG5 functions in membrane remodeling during autophagy, UGP2 produces the substrate for cellulose and callose polymer extension, and FAMA is a transcription factor. Automated phenotyping enabled this discovery of four naturally varying components of a conserved process (gravitropism) by making it feasible to conduct the same large-scale experiment in two species.

High-throughput measurement of plant fitness traits with an object detection method using Faster R-CNN

Revealing the contributions of genes to plant phenotype is frequently challenging because loss-of-function effects may be subtle or masked by varying degrees of genetic redundancy. Such effects can potentially be detected by measuring plant fitness, which reflects the cumulative effects of genetic changes over the lifetime of a plant. However, fitness is challenging to measure accurately, particularly in species with high fecundity and relatively small propagule sizes such as Arabidopsis thaliana. An image segmentation-based method using the software ImageJ and an object detection-based method using the Faster Region-based Convolutional Neural Network (R-CNN) algorithm were used for measuring two Arabidopsis fitness traits: seed and fruit counts. The segmentation-based method was error-prone (correlation between true and predicted seed counts, r²  = 0.849) because seeds touching each other were undercounted. By contrast, the object detection-based algorithm yielded near perfect seed counts (r²  = 0.9996) and highly accurate fruit counts (r²  = 0.980). Comparing seed counts for wild-type and 12 mutant lines revealed fitness effects for three genes; fruit counts revealed the same effects for two genes. Our study provides analysis pipelines and models to facilitate the investigation of Arabidopsis fitness traits and demonstrates the importance of examining fitness traits when studying gene functions.

Phenotypic plasticity in plant height shaped by interaction between genetic loci and diurnal temperature range

Phenotypic plasticity is observed widely in plants and often studied with reaction norms for adult plant or end-of-season traits. Uncovering genetic, environmental and developmental patterns behind the observed phenotypic variation under natural field conditions is needed. Using a sorghum (Sorghum bicolor) genetic population evaluated for plant height in seven natural field conditions, we investigated the major pattern that differentiated these environments. We then examined the physiological relevance of the identified environmental index by investigating the developmental trajectory of the population with multistage height measurements in four additional environments and conducting crop growth modelling. We found that diurnal temperature range (DTR) during the rapid growth period of sorghum development was an effective environmental index. Three genetic loci (Dw1, Dw3 and qHT7.1) were consistently detected for individual environments, reaction-norm parameters across environments and growth-curve parameters through the season. Their genetic effects changed dynamically along the environmental gradient and the developmental stage. A conceptual model with three-dimensional reaction norms was proposed to showcase the interconnecting components: genotype, environment and development. Beyond genomic and environmental analyses, further integration of development and physiology at the whole-plant and molecular levels into complex trait dissection would enhance our understanding of mechanisms underlying phenotypic variation.

Statistical Methods for the Quantitative Genetic Analysis of High-Throughput Phenotyping Data

The advent of plant phenomics, coupled with the wealth of genotypic data generated by next-generation sequencing technologies, provides exciting new resources for investigations into and improvement of complex traits. However, these new technologies also bring new challenges in quantitative genetics, namely, a need for the development of robust frameworks that can accommodate these high-dimensional data. In this chapter, we describe methods for the statistical analysis of high-throughput phenotyping (HTP) data with the goal of enhancing the prediction accuracy of genomic selection (GS). Following the Introduction in Sec. 1, Sec. 2 discusses field-based HTP, including the use of unoccupied aerial vehicles and light detection and ranging, as well as how we can achieve increased genetic gain by utilizing image data derived from HTP. Section 3 considers extending commonly used GS models to integrate HTP data as covariates associated with the principal trait response, such as yield. Particular focus is placed on single-trait, multi-trait, and genotype by environment interaction models. One unique aspect of HTP data is that phenomics platforms often produce large-scale data with high spatial and temporal resolution for capturing dynamic growth, development, and stress responses. Section 4 discusses the utility of a random regression model for performing longitudinal modeling. The chapter concludes with a discussion of some standing issues.

High-resolution imaging as a tool for identifying quantitative trait loci that regulate photomorphogenesis in Arabidopsis thaliana

A primary component of seedling establishment is the photomorphogenic response as seedlings emerge from the soil. This process is characterized by a reduced growth rate in the hypocotyl, increased root growth, opening of the apical hook and expansion of the cotyledons as photosynthetic organs. While fundamental to plant success, the photomorphogenic response can be highly variable. Additionally, studies of Arabidopsis thaliana are made difficult by subtle differences in growth rate between individuals. High-resolution imaging and computational processing have emerged as useful tools for quantification of such phenotypes. This study sought to: (i) develop an imaging methodology which could capture changes in growth rate as seedlings transition from darkness to blue light in real time, and (ii) apply this methodology to single-quantitative trait locus (QTL) analysis using the Cvi × Ler recombinant inbred line (RIL) mapping population. Significant differences in the photomorphogenic response were observed between the parent lines and analysis of 158 RILs revealed a wide range of growth rate phenotypes. Quantitative trait locus analysis detected significant loci associated with dark growth rate on chromosome 5 and significant loci associated with light growth rate on chromosome 2. Candidate genes associated with these loci, such as the previously characterized ER locus, highlight the application of this approach for QTL analysis. Genetic analysis of Landsberg lines without the erecta mutation also supports a role for ER in modulating the photomorphogenic response, consistent with previous QTL analyses of this population. Strengths and limitations of this methodology are presented, as well as means of improvement.

Estimation of dynamic SNP-heritability with Bayesian Gaussian process models

Motivation

Improved DNA technology has made it practical to estimate single-nucleotide polymorphism (SNP)-heritability among distantly related individuals with unknown relationships. For growth- and development-related traits, it is meaningful to base SNP-heritability estimation on longitudinal data due to the time-dependency of the process. However, only few statistical methods have been developed so far for estimating dynamic SNP-heritability and quantifying its full uncertainty.

Results

We introduce a completely tuning-free Bayesian Gaussian process (GP)-based approach for estimating dynamic variance components and heritability as their function. For parameter estimation, we use a modern Markov Chain Monte Carlo method which allows full uncertainty quantification. Several datasets are analysed and our results clearly illustrate that the 95% credible intervals of the proposed joint estimation method (which ‘borrows strength’ from adjacent time points) are significantly narrower than of a two-stage baseline method that first estimates the variance components at each time point independently and then performs smoothing. We compare the method with a random regression model using MTG2 and BLUPF90 software and quantitative measures indicate superior performance of our method. Results are presented for simulated and real data with up to 1000 time points. Finally, we demonstrate scalability of the proposed method for simulated data with tens of thousands of individuals.

Availability and implementation

The C++ implementation dynBGP and simulated data are available in GitHub: https://github.com/aarjas/dynBGP. The programmes can be run in R. Real datasets are available in QTL archive: https://phenome.jax.org/centers/QTLA.

Supplementary information

Supplementary data are available at Bioinformatics online.

Identification of quantitative trait loci associated with canopy temperature in soybean

A consistent risk for soybean (Glycine max L.) production is the impact of drought on growth and yield. Canopy temperature (CT) is an indirect measure of transpiration rate and stomatal conductance and may be valuable in distinguishing differences among genotypes in response to drought. The objective of this study was to map quantitative trait loci (QTLs) associated with CT using thermal infrared imaging in a population of recombinant inbred lines developed from a cross between KS4895 and Jackson. Heritability of CT was 35% when estimated across environments. QTL analysis identified 11 loci for CT distributed on eight chromosomes that individually explained between 4.6 and 12.3% of the phenotypic variation. The locus on Gm11 was identified in two individual environments and across environments and explained the highest proportion of phenotypic variation (9.3% to 11.5%) in CT. Several of these CT loci coincided with the genomic regions from previous studies associated with canopy wilting, canopy temperature, water use efficiency, and other morpho-physiological traits related with drought tolerance. Candidate genes with biological function related to transpiration, root development, and signal transduction underlie these putative CT loci. These genomic regions may be important resources in soybean breeding programs to improve tolerance to drought.

Increased Power and Accuracy of Causal Locus Identification in Time Series Genome-wide Association in Sorghum

The phenotypes of plants develop over time and change in response to the environment. New engineering and computer vision technologies track these phenotypic changes. Identifying the genetic loci regulating differences in the pattern of phenotypic change remains challenging. This study used functional principal component analysis (FPCA) to achieve this aim. Time series phenotype data were collected from a sorghum (Sorghum bicolor) diversity panel using a number of technologies including conventional color photography and hyperspectral imaging. This imaging lasted for 37 d and centered on reproductive transition. A new higher density marker set was generated for the same population. Several genes known to control trait variation in sorghum have been previously cloned and characterized. These genes were not confidently identified in genome-wide association analyses at single time points. However, FPCA successfully identified the same known and characterized genes. FPCA analyses partitioned the role these genes play in controlling phenotypes. Partitioning was consistent with the known molecular function of the individual cloned genes. These data demonstrate that FPCA-based genome-wide association studies can enable robust time series mapping analyses in a wide range of contexts. Moreover, time series analysis can increase the accuracy and power of quantitative genetic analyses.

The pulse of the tree is under genetic control: eucalyptus as a case study

The pulse of the tree (diurnal cycle of stem radius fluctuations) has been widely studied as a way of analyzing tree responses to the environment, including the phenotypic plasticity of tree-water relationships in particular. However, the genetic basis of this daily phenotype and its interplay with the environment remain largely unexplored. We characterized the genetic and environmental determinants of this response, by monitoring daily stem radius fluctuation (dSRF) on 210 trees from a Eucalyptus urophylla × E. grandis full-sib family over 2 years. The dSRF signal was broken down into hydraulic capacitance, assessed as the daily amplitude of shrinkage (DA), and net growth, estimated as the change in maximum radius between two consecutive days (ΔR). The environmental determinants of these two traits were clearly different: DA was positively correlated with atmospheric variables relating to water demand, while ΔR was associated with soil water content. The heritability for these two traits ranged from low to moderate over time, revealing a time-dependent or environment-dependent complex genetic determinism. We identified 686 and 384 daily quantitative trait loci (QTL) representing 32 and 31 QTL regions for DA and ΔR, respectively. The identification of gene networks underlying the 27 major genomics regions for both traits generated additional hypotheses concerning the biological mechanisms involved in response to water demand and supply. This study highlights that environmentally induced changes in daily stem radius fluctuation are genetically controlled in trees and suggests that these daily responses integrated over time shape the genetic architecture of mature traits.

A Gaussian process model and Bayesian variable selection for mapping function-valued quantitative traits with incomplete phenotypic data

Motivation

Recent advances in high dimensional phenotyping bring time as an extra dimension into the phenotypes. This promotes the quantitative trait locus (QTL) studies of function-valued traits such as those related to growth and development. Existing approaches for analyzing functional traits utilize either parametric methods or semi-parametric approaches based on splines and wavelets. However, very limited choices of software tools are currently available for practical implementation of functional QTL mapping and variable selection.

Results

We propose a Bayesian Gaussian process (GP) approach for functional QTL mapping. We use GPs to model the continuously varying coefficients which describe how the effects of molecular markers on the quantitative trait are changing over time. We use an efficient gradient based algorithm to estimate the tuning parameters of GPs. Notably, the GP approach is directly applicable to the incomplete datasets having even larger than 50% missing data rate (among phenotypes). We further develop a stepwise algorithm to search through the model space in terms of genetic variants, and use a minimal increase of Bayesian posterior probability as a stopping rule to focus on only a small set of putative QTL. We also discuss the connection between GP and penalized B-splines and wavelets. On two simulated and three real datasets, our GP approach demonstrates great flexibility for modeling different types of phenotypic trajectories with low computational cost. The proposed model selection approach finds the most likely QTL reliably in tested datasets.

Availability and implementation

Software and simulated data are available as a MATLAB package ‘GPQTLmapping’, and they can be downloaded from GitHub (https://github.com/jpvanhat/GPQTLmapping). Real datasets used in case studies are publicly available at QTL Archive.

Supplementary information

Supplementary data are available at Bioinformatics online.

Functional QTL mapping and genomic prediction of canopy height in wheat measured using a robotic field phenotyping platform

Genetic studies increasingly rely on high-throughput phenotyping, but the resulting longitudinal data pose analytical challenges. We used canopy height data from an automated field phenotyping platform to compare several approaches to scanning for quantitative trait loci (QTLs) and performing genomic prediction in a wheat recombinant inbred line mapping population based on up to 26 sampled time points (TPs). We detected four persistent QTLs (i.e. expressed for most of the growing season), with both empirical and simulation analyses demonstrating superior statistical power of detecting such QTLs through functional mapping approaches compared with conventional individual TP analyses. In contrast, even very simple individual TP approaches (e.g. interval mapping) had superior detection power for transient QTLs (i.e. expressed during very short periods). Using spline-smoothed phenotypic data resulted in improved genomic predictive abilities (5-8% higher than individual TP prediction), while the effect of including significant QTLs in prediction models was relatively minor (<1-4% improvement). Finally, although QTL detection power and predictive ability generally increased with the number of TPs analysed, gains beyond five or 10 TPs chosen based on phenological information had little practical significance. These results will inform the development of an integrated, semi-automated analytical pipeline, which will be more broadly applicable to similar data sets in wheat and other crops.

Functional QTL mapping and genomic prediction of canopy height in wheat measured using a robotic field phenotyping platform

Genetic studies increasingly rely on high-throughput phenotyping, but the resulting longitudinal data pose analytical challenges. We used canopy height data from an automated field phenotyping platform to compare several approaches to scanning for quantitative trait loci (QTLs) and performing genomic prediction in a wheat recombinant inbred line mapping population based on up to 26 sampled time points (TPs). We detected four persistent QTLs (i.e. expressed for most of the growing season), with both empirical and simulation analyses demonstrating superior statistical power of detecting such QTLs through functional mapping approaches compared with conventional individual TP analyses. In contrast, even very simple individual TP approaches (e.g. interval mapping) had superior detection power for transient QTLs (i.e. expressed during very short periods). Using spline-smoothed phenotypic data resulted in improved genomic predictive abilities (5-8% higher than individual TP prediction), while the effect of including significant QTLs in prediction models was relatively minor (<1-4% improvement). Finally, although QTL detection power and predictive ability generally increased with the number of TPs analysed, gains beyond five or 10 TPs chosen based on phenological information had little practical significance. These results will inform the development of an integrated, semi-automated analytical pipeline, which will be more broadly applicable to similar data sets in wheat and other crops.

Crop Phenomics and High-Throughput Phenotyping: Past Decades, Current Challenges, and Future Perspectives

Since whole-genome sequencing of many crops has been achieved, crop functional genomics studies have stepped into the big-data and high-throughput era. However, acquisition of large-scale phenotypic data has become one of the major bottlenecks hindering crop breeding and functional genomics studies. Nevertheless, recent technological advances provide us potential solutions to relieve this bottleneck and to explore advanced methods for large-scale phenotyping data acquisition and processing in the coming years. In this article, we review the major progress on high-throughput phenotyping in controlled environments and field conditions as well as its use for post-harvest yield and quality assessment in the past decades. We then discuss the latest multi-omics research combining high-throughput phenotyping with genetic studies. Finally, we propose some conceptual challenges and provide our perspectives on how to bridge the phenotype-genotype gap. It is no doubt that accurate high-throughput phenotyping will accelerate plant genetic improvements and promote the next green revolution in crop breeding.

Three-Dimensional Time-Lapse Analysis Reveals Multiscale Relationships in Maize Root Systems with Contrasting Architectures

Understanding how an organism's phenotypic traits are conditioned by genetic and environmental variation is a central goal of biology. Root systems are one of the most important but poorly understood aspects of plants, largely due to the three-dimensional (3D), dynamic, and multiscale phenotyping challenge they pose. A critical gap in our knowledge is how root systems build in complexity from a single primary root to a network of thousands of roots that collectively compete for ephemeral, heterogeneous soil resources. We used time-lapse 3D imaging and mathematical modeling to assess root system architectures (RSAs) of two maize (Zea mays) inbred genotypes and their hybrid as they grew in complexity from a few to many roots. Genetically driven differences in root branching zone size and lateral branching densities along a single root, combined with differences in peak growth rate and the relative allocation of carbon resources to new versus existing roots, manifest as sharply distinct global RSAs over time. The 3D imaging of mature field-grown root crowns showed that several genetic differences in seedling architectures could persist throughout development and across environments. This approach connects individual and system-wide scales of root growth dynamics, which could eventually be used to predict genetic variation for complex RSAs and their functions.

Three-Dimensional Time-Lapse Analysis Reveals Multiscale Relationships in Maize Root Systems with Contrasting Architectures

Understanding how an organism's phenotypic traits are conditioned by genetic and environmental variation is a central goal of biology. Root systems are one of the most important but poorly understood aspects of plants, largely due to the three-dimensional (3D), dynamic, and multiscale phenotyping challenge they pose. A critical gap in our knowledge is how root systems build in complexity from a single primary root to a network of thousands of roots that collectively compete for ephemeral, heterogeneous soil resources. We used time-lapse 3D imaging and mathematical modeling to assess root system architectures (RSAs) of two maize (Zea mays) inbred genotypes and their hybrid as they grew in complexity from a few to many roots. Genetically driven differences in root branching zone size and lateral branching densities along a single root, combined with differences in peak growth rate and the relative allocation of carbon resources to new versus existing roots, manifest as sharply distinct global RSAs over time. The 3D imaging of mature field-grown root crowns showed that several genetic differences in seedling architectures could persist throughout development and across environments. This approach connects individual and system-wide scales of root growth dynamics, which could eventually be used to predict genetic variation for complex RSAs and their functions.

Leveraging Breeding Values Obtained from Random Regression Models for Genetic Inference of Longitudinal Traits

Understanding the genetic basis of dynamic plant phenotypes has largely been limited because of a lack of space and labor resources needed to record dynamic traits, often destructively, for a large number of genotypes. However, the recent advent of image-based phenotyping platforms has provided the plant science community with an effective means to nondestructively evaluate morphological, developmental, and physiological processes at regular, frequent intervals for a large number of plants throughout development. The statistical frameworks typically used for genetic analyses (e.g., genome-wide association mapping, linkage mapping, and genomic prediction) in plant breeding and genetics are not particularly amenable for repeated measurements. Random regression (RR) models are routinely used in animal breeding for the genetic analysis of longitudinal traits and provide a robust framework for modeling trait trajectories and performing genetic analysis simultaneously. We recently used a RR approach for genomic prediction of shoot growth trajectories in rice ( L.) from 33,674 single nucleotide polymorphisms. In this study, we have extended this approach for genetic inference by leveraging genomic breeding values derived from RR models for rice shoot growth during early vegetative development. This approach provides improvements over conventional single time point analyses for discovering loci associated with shoot growth trajectories. The RR approach uncovers persistent as well as time-specific transient quantitative trait loci. This methodology can be widely applied to understand the genetic architecture of other complex polygenic traits with repeated measurements.

The complex genetic architecture of shoot growth natural variation in Arabidopsis thaliana

One of the main outcomes of quantitative genetics approaches to natural variation is to reveal the genetic architecture underlying the phenotypic space. Complex genetic architectures are described as including numerous loci (or alleles) with small-effect and/or low-frequency in the populations, interactions with the genetic background, environment or age. Linkage or association mapping strategies will be more or less sensitive to this complexity, so that we still have an unclear picture of its extent. By combining high-throughput phenotyping under two environmental conditions with classical QTL mapping approaches in multiple Arabidopsis thaliana segregating populations as well as advanced near isogenic lines construction and survey, we have attempted to improve our understanding of quantitative phenotypic variation. Integrative traits such as those related to vegetative growth used in this work (highlighting either cumulative growth, growth rate or morphology) all showed complex and dynamic genetic architecture with respect to the segregating population and condition. The more resolutive our mapping approach, the more complexity we uncover, with several instances of QTLs visible in near isogenic lines but not detected with the initial QTL mapping, indicating that our phenotyping accuracy was less limiting than the mapping resolution with respect to the underlying genetic architecture. In an ultimate approach to resolve this complexity, we intensified our phenotyping effort to target specifically a 3Mb-region known to segregate for a major quantitative trait gene, using a series of selected lines recombined every 100kb. We discovered that at least 3 other independent QTLs had remained hidden in this region, some with trait- or condition-specific effects, or opposite allelic effects. If we were to extrapolate the figures obtained on this specific region in this particular cross to the genome- and species-scale, we would predict hundreds of causative loci of detectable phenotypic effect controlling these growth-related phenotypes.

Genome-wide association of barley plant growth under drought stress using a nested association mapping population

BACKGROUND

Barley (Hordeum vulgare L.) is the fourth most important cereal crop worldwide. Barley production is compromised by many abiotic stresses including drought. Wild barley is a valuable source of alleles that can improve adaptation of cultivated barley to drought stress.

RESULTS

In the present study, a nested association mapping population named HEB-25, consisting of 1420 BC₁S₃ lines that were developed by crossing 25 different wild barley accessions to the elite barley cultivar 'Barke', was evaluated under both control and drought-stressed conditions in the Australian Plant Phenomics Facility, University of Adelaide. Overall, 14 traits reflecting the performance of individual plants in each treatment were calculated from non-destructive imaging over time and destructive end-of-experiment measurements. For each trait, best linear unbiased estimators (BLUEs) were calculated and used for genome-wide association study (GWAS) analysis. Among the quantitative trait loci (QTL) identified for the 14 traits, many co-localise with known inflorescence and developmental genes. We identified a QTL on chromosome 4H where, under drought and control conditions, wild barley alleles increased biomass by 10 and 17% respectively compared to the Barke allele.

CONCLUSIONS

Across all traits, QTL which increased phenotypic values were identified, providing a wider range of genetic diversity for the improvement of drought tolerance in barley.

Time-specific and pleiotropic quantitative trait loci coordinately modulate stem growth in Populus

In perennial woody plants, the coordinated increase of stem height and diameter during juvenile growth improves competitiveness (i.e. access to light); however, the factors underlying variation in stem growth remain unknown in trees. Here, we used linkage-linkage disequilibrium (linkage-LD) mapping to decipher the genetic architecture underlying three growth traits during juvenile stem growth. We used two Populus populations: a linkage mapping population comprising a full-sib family of 1,200 progeny and an association mapping panel comprising 435 unrelated individuals from nearly the entire natural range of Populus tomentosa. We mapped 311 quantitative trait loci (QTL) for three growth traits at 12 timepoints to 42 regions in 17 linkage groups. Of these, 28 regions encompassing 233 QTL were annotated as 27 segmental homology regions (SHRs). Using SNPs identified by whole-genome re-sequencing of the 435-member association mapping panel, we identified significant SNPs (P ≤ 9.4 × 10^-7 ) within 27 SHRs that affect stem growth at nine timepoints with diverse additive and dominance patterns, and these SNPs exhibited complex allelic epistasis over the juvenile growth period. Nineteen genes linked to potential causative alleles that have time-specific or pleiotropic effects, and mostly overlapped with significant signatures of selection within SHRs between climatic regions represented by the association mapping panel. Five genes with potential time-specific effects showed species-specific temporal expression profiles during the juvenile stages of stem growth in five representative Populus species. Our observations revealed the importance of considering temporal genetic basis of complex traits, which will facilitate the molecular design of tree ideotypes.

Genome-Wide Association Studies to Improve Wood Properties: Challenges and Prospects

Wood formation is an excellent model system for quantitative trait analysis due to the strong associations between the transcriptional and metabolic traits that contribute to this complex process. Investigating the genetic architecture and regulatory mechanisms underlying wood formation will enhance our understanding of the quantitative genetics and genomics of complex phenotypic variation. Genome-wide association studies (GWASs) represent an ideal statistical strategy for dissecting the genetic basis of complex quantitative traits. However, elucidating the molecular mechanisms underlying many favorable loci that contribute to wood formation and optimizing GWAS design remain challenging in this omics era. In this review, we summarize the recent progress in GWAS-based functional genomics of wood property traits in major timber species such as Eucalyptus, Populus, and various coniferous species. We discuss several appropriate experimental designs for extensive GWAS in a given undomesticated tree population, such as omics-wide association studies and high-throughput phenotyping technologies. We also explain why more attention should be paid to rare allelic and major structural variation. Finally, we explore the potential use of GWAS for the molecular breeding of trees. Such studies will help provide an integrated understanding of complex quantitative traits and should enable the molecular design of new cultivars.

MuSeeQ, a novel supervised image analysis tool for the simultaneous phenotyping of the soluble mucilage and seed morphometric parameters

BACKGROUND

The mucilage is a model to study the polysaccharide biosynthesis since it is produced in large amounts and composed of complex polymers. In addition, it is of great economic interest for its technical and nutritional value. A fast method for phenotyping the released mucilage and the seed morphometric parameters will be useful for fundamental, food, pharmaceutical and breeding researches. Current strategies to phenotype soluble mucilage are restricted to visual evaluations or are highly time-consuming.

RESULTS

Here, we developed a high-throughput phenotyping method for the simultaneous measurement of the soluble mucilage content released on a gel and the seed morphometric parameters. Within this context, we combined a biochemical assay and an open-source computer-aided image analysis tool, MuSeeQ. The biochemical assay consists in sowing seeds on an agarose medium containing the dye toluidine blue O, which specifically stains the mucilage once it is released on the gel. The second part of MuSeeQ is a macro developed in ImageJ allowing to quickly extract and analyse 11 morphometric data of seeds and their respective released mucilages. As an example, MuSeeQ was applied on a flax recombinant inbred lines population (previously screened for fatty acids content.) and revealed significant correlations between the soluble mucilage shape and the concentration of some fatty acids, e.g. C16:0 and C18:2. Other fatty acids were also found to correlate with the seed shape parameters, e.g. C18:0 and C18:2. MuSeeQ was then showed to be used for the analysis of other myxospermous species, including Arabidopsis thaliana and Camelina sativa.

CONCLUSIONS

MuSeeQ is a low-cost and user-friendly method which may be used by breeders and researchers for phenotyping simultaneously seeds of specific cultivars, natural variants or mutants and their respective soluble mucilage area released on a gel. The script of MuSeeQ and video tutorials are freely available at http://MuSeeQ.free.fr.

Regulation of Root Angle and Gravitropism

Regulation of plant root angle is critical for obtaining nutrients and water and is an important trait for plant breeding. A plant’s final, long-term root angle is the net result of a complex series of decisions made by a root tip in response to changes in nutrient availability, impediments, the gravity vector and other stimuli. When a root tip is displaced from the gravity vector, the short-term process of gravitropism results in rapid reorientation of the root toward the vertical. Here, we explore both short- and long-term regulation of root growth angle, using natural variation in tomato to identify shared and separate genetic features of the two responses. Mapping of expression quantitative trait loci mapping and leveraging natural variation between and within species including Arabidopsis suggest a role for PURPLE ACID PHOSPHATASE 27 and CELL DIVISION CYCLE 73 in determining root angle.

ShinyAIM: Shiny-based application of interactive Manhattan plots for longitudinal genome-wide association studies

Owning to advancements in sensor-based, non-destructive phenotyping platforms, researchers are increasingly collecting data with higher temporal resolution. These phenotypes collected over several time points are cataloged as longitudinal traits and used for genome-wide association studies (GWAS). Longitudinal GWAS typically yield a large number of output files, posing a significant challenge to data interpretation and visualization. Efficient, dynamic, and integrative data visualization tools are essential for the interpretation of longitudinal GWAS results for biologists; however, these tools are not widely available to the community. We have developed a flexible and user-friendly Shiny-based online application, ShinyAIM, to dynamically view and interpret temporal GWAS results. The main features of the application include (a) interactive Manhattan plots for single time points, (b) a grid plot to view Manhattan plots for all time points simultaneously, (c) dynamic scatter plots for p-value-filtered selected markers to investigate co-localized genomic regions across time points, (d) and interactive phenotypic data visualization to capture variation and trends in phenotypes. The application is written entirely in the R language and can be used with limited programming experience. ShinyAIM is deployed online as a Shiny web server application at https://chikudaisei.shinyapps.io/shinyaim/, enabling easy access for users without installation. The application can also be launched on a local machine in RStudio.

The Quantitative Genetic Control of Root Architecture in Maize

Roots remain an underexplored frontier in plant genetics despite their well-known influence on plant development, agricultural performance and competition in the wild. Visualizing and measuring root structures and their growth is vastly more difficult than characterizing aboveground parts of the plant and is often simply avoided. The majority of research on maize root systems has focused on their anatomy, physiology, development and soil interaction, but much less is known about the genetics that control quantitative traits. In maize, seven root development genes have been cloned using mutagenesis, but no genes underlying the many root-related quantitative trait loci (QTLs) have been identified. In this review, we discuss whether the maize mutants known to control root development may also influence quantitative aspects of root architecture, including the extent to which they overlap with the most recent maize root trait QTLs. We highlight specific challenges and anticipate the impacts that emerging technologies, especially computational approaches, may have toward the identification of genes controlling root quantitative traits.

Natural variation of YELLOW SEEDLING1 affects photosynthetic acclimation of Arabidopsis thaliana

Exploiting genetic variation for more efficient photosynthesis is an underexplored route towards new crop varieties. This study demonstrates the genetic dissection of higher plant photosynthesis efficiency down to the genomic DNA level, by confirming that allelic sequence variation at the Arabidopsis thaliana YELLOW SEEDLING1 (YS1) gene explains natural diversity in photosynthesis acclimation to high irradiance. We use a genome-wide association study to identify quantitative trait loci (QTLs) involved in the Arabidopsis photosynthetic acclimation response. Candidate genes underlying the QTLs are prioritized according to functional clues regarding gene ontology, expression and function. Reverse genetics and quantitative complementation confirm the candidacy of YS1, which encodes a pentatrico-peptide-repeat (PPR) protein involved in RNA editing of plastid-encoded genes (anterograde signalling). Gene expression analysis and allele sequence comparisons reveal polymorphisms in a light-responsive element in the YS1 promoter that affect its expression, and that of its downstream targets, resulting in the variation in photosynthetic acclimation.

Natural genetic variation of photosynthesis is an underexplored resource for plant genetic improvement. Here, the authors find allelic variations of YS1 affect Arabidopsis photosynthesis acclimation using genome-wide association study, reverse genetics, and quantitative complementation approaches.

Digital Biomass Accumulation Using High-Throughput Plant Phenotype Data Analysis

Biomass is an important phenotypic trait in functional ecology and growth analysis. The typical methods for measuring biomass are destructive, and they require numerous individuals to be cultivated for repeated measurements. With the advent of image-based high-throughput plant phenotyping facilities, non-destructive biomass measuring methods have attempted to overcome this problem. Thus, the estimation of plant biomass of individual plants from their digital images is becoming more important. In this paper, we propose an approach to biomass estimation based on image derived phenotypic traits. Several image-based biomass studies state that the estimation of plant biomass is only a linear function of the projected plant area in images. However, we modeled the plant volume as a function of plant area, plant compactness, and plant age to generalize the linear biomass model. The obtained results confirm the proposed model and can explain most of the observed variance during image-derived biomass estimation. Moreover, a small difference was observed between actual and estimated digital biomass, which indicates that our proposed approach can be used to estimate digital biomass accurately.

A Comprehensive Image-based Phenomic Analysis Reveals the Complex Genetic Architecture of Shoot Growth Dynamics in Rice (Oryza sativa)

Early vigor is an important trait for many rice ( L.)-growing environments. However, genetic characterization and improvement for early vigor is hindered by the temporal nature of the trait and strong genotype × environment effects. We explored the genetic architecture of shoot growth dynamics during the early and active tillering stages by applying a functional modeling and genomewide association (GWAS) mapping approach on a diversity panel of ∼360 rice accessions. Multiple loci with small effects on shoot growth trajectory were identified, indicating a complex polygenic architecture. Natural variation for shoot growth dynamics was assessed in a subset of 31 accessions using RNA sequencing and hormone quantification. These analyses yielded a gibberellic acid (GA) catabolic gene, , which could influence GA levels to regulate vigor in the early tillering stage. Given the complex genetic architecture of shoot growth dynamics, the potential of genomic selection (GS) for improving early vigor was explored using all 36,901 single-nucleotide polymorphisms (SNPs) as well as several subsets of the most significant SNPs from GWAS. Shoot growth trajectories could be predicted with reasonable accuracy using the 50 most significant SNPs from GWAS (0.37-0.53); however, the accuracy of prediction was improved by including more markers, which indicates that GS may be an effective strategy for improving shoot growth dynamics during the vegetative growth stage. This study provides insights into the complex genetic architecture and molecular mechanisms underlying early shoot growth dynamics and provides a foundation for improving this complex trait in rice.

Time dependent genetic analysis links field and controlled environment phenotypes in the model C4 grass Setaria

Vertical growth of plants is a dynamic process that is influenced by genetic and environmental factors and has a pronounced effect on overall plant architecture and biomass composition. We have performed six controlled growth trials of an interspecific Setaria italica x Setaria viridis recombinant inbred line population to assess how the genetic architecture of plant height is influenced by developmental queues, water availability and planting density. The non-destructive nature of plant height measurements has enabled us to monitor height throughout the plant life cycle in both field and controlled environments. We find that plant height is reduced under water limitation and high density planting and affected by growth environment (field vs. growth chamber). The results support a model where plant height is a heritable, polygenic trait and that the major genetic loci that influence plant height function independent of growth environment. The identity and contribution of loci that influence height changes dynamically throughout development and the reduction of growth observed in water limited environments is a consequence of delayed progression through the genetic program which establishes plant height in Setaria. In this population, alleles inherited from the weedy S. viridis parent act to increase plant height early, whereas a larger number of small effect alleles inherited from the domesticated S. italica parent collectively act to increase plant height later in development.

Growth is a dynamic process that responds to a changing environment. Most of the methods that we have for measuring are static and collecting information throughout an organisms lifecycle is labor and cost prohibitive. Advances in imaging and robotics technology have enabled novel approaches to understanding how plants adapt to the environment. Using the model grass Setaria and new methods for measuring parameters from images, we investigate the genetic architecture of plant height in response to water availability and planting density. Height is one of the most influential components of plant architecture, determining tradeoffs between competition and resource allocation and is an important trait for boosting yields. The non-destructive nature of plant height measurements has enabled us to monitor growth throughout the plant life cycle in both field and controlled environments. We identified several loci controlling height in a population derived from a wild strain of Setaria viridis and its domesticated relative Setaria italica, as well as the developmental time in which these loci act. In this population, alleles inherited from the wild parent act to increase plant height early, whereas a larger number of small effect alleles inherited from the domesticated parent collectively act to increase plant height later in development.

Applying hyperspectral imaging to explore natural plant diversity towards improving salt stress tolerance

Salinity represents an abiotic stress constraint affecting growth and productivity of plants in many regions of the world. One of the possible solutions is to improve the level of salt resistance using natural genetic variability within crop species. In the context of recent knowledge on salt stress effects and mechanisms of salt tolerance, this review present useful phenomic approach employing different non-invasive imaging systems for detection of quantitative and qualitative changes caused by salt stress at the plant and canopy level. The focus is put on hyperspectral imaging technique, which provides unique opportunities for fast and reliable estimate of numerous characteristics associated both with various structural, biochemical and physiological traits. The method also provides possibilities to combine plant and canopy analyses with a direct determination of salinity in soil. The future perspectives in salt stress applications as well as some limits of the method are also identified.

Genetic variation of growth dynamics in maize (Zea mays L.) revealed through automated non-invasive phenotyping

Hitherto, most quantitative trait loci of maize growth and biomass yield have been identified for a single time point, usually the final harvest stage. Through this approach cumulative effects are detected, without considering genetic factors causing phase-specific differences in growth rates. To assess the genetics of growth dynamics, we employed automated non-invasive phenotyping to monitor the plant sizes of 252 diverse maize inbred lines at 11 different developmental time points; 50 k SNP array genotype data were used for genome-wide association mapping and genomic selection. The heritability of biomass was estimated to be over 71%, and the average prediction accuracy amounted to 0.39. Using the individual time point data, 12 main effect marker-trait associations (MTAs) and six pairs of epistatic interactions were detected that displayed different patterns of expression at various developmental time points. A subset of them also showed significant effects on relative growth rates in different intervals. The detected MTAs jointly explained up to 12% of the total phenotypic variation, decreasing with developmental progression. Using non-parametric functional mapping and multivariate mapping approaches, four additional marker loci affecting growth dynamics were detected. Our results demonstrate that plant biomass accumulation is a complex trait governed by many small effect loci, most of which act at certain restricted developmental phases. This highlights the need for investigation of stage-specific growth affecting genes to elucidate important processes operating at different developmental phases.

How can we harness quantitative genetic variation in crop root systems for agricultural improvement?

Root systems are a black box obscuring a comprehensive understanding of plant function, from the ecosystem scale down to the individual. In particular, a lack of knowledge about the genetic mechanisms and environmental effects that condition root system growth hinders our ability to develop the next generation of crop plants for improved agricultural productivity and sustainability. We discuss how the methods and metrics we use to quantify root systems can affect our ability to understand them, how we can bridge knowledge gaps and accelerate the derivation of structure-function relationships for roots, and why a detailed mechanistic understanding of root growth and function will be important for future agricultural gains.

Phenomics for photosynthesis, growth and reflectance in Arabidopsis thaliana reveals circadian and long-term fluctuations in heritability

BACKGROUND

Recent advances in genome sequencing technologies have shifted the research bottleneck in plant sciences from genotyping to phenotyping. This shift has driven the development of phenomics, high-throughput non-invasive phenotyping technologies.

RESULTS

We describe an automated high-throughput phenotyping platform, the Phenovator, capable of screening 1440 Arabidopsis plants multiple times per day for photosynthesis, growth and spectral reflectance at eight wavelengths. Using this unprecedented phenotyping capacity, we have been able to detect significant genetic differences between Arabidopsis accessions for all traits measured, across both temporal and environmental scales. The high frequency of measurement allowed us to observe that heritability was not only trait specific, but for some traits was also time specific.

CONCLUSIONS

Such continuous real-time non-destructive phenotyping will allow detailed genetic and physiological investigations of the kinetics of plant homeostasis and development. The success and ultimate outcome of a breeding program will depend greatly on the genetic variance which is sampled. Our observation of temporal fluctuations in trait heritability shows that the moment of measurement can have lasting consequences. Ultimately such phenomic level technologies will provide more dynamic insights into plant physiology, and the necessary data for the omics revolution to reach its full potential.

Dynamic Quantitative Trait Locus Analysis of Plant Phenomic Data

Advanced platforms have recently become available for automatic and systematic quantification of plant growth and development. These new techniques can efficiently produce multiple measurements of phenotypes over time, and introduce time as an extra dimension to quantitative trait locus (QTL) studies. Functional mapping utilizes a class of statistical models for identifying QTLs associated with the growth characteristics of interest. A major benefit of functional mapping is that it integrates information over multiple timepoints, and therefore could increase the statistical power for QTL detection. We review the current development of computationally efficient functional mapping methods which provide invaluable tools for analyzing large-scale timecourse data that are readily available in our post-genome era.

Mapping Quantitative Trait Loci Underlying Function-Valued Traits Using Functional Principal Component Analysis and Multi-Trait Mapping

We previously proposed a simple regression-based method to map quantitative trait loci underlying function-valued phenotypes. In order to better handle the case of noisy phenotype measurements and accommodate the correlation structure among time points, we propose an alternative approach that maintains much of the simplicity and speed of the regression-based method. We overcome noisy measurements by replacing the observed data with a smooth approximation. We then apply functional principal component analysis, replacing the smoothed phenotype data with a small number of principal components. Quantitative trait locus mapping is applied to these dimension-reduced data, either with a multi-trait method or by considering the traits individually and then taking the average or maximum LOD score across traits. We apply these approaches to root gravitropism data on Arabidopsis recombinant inbred lines and further investigate their performance in computer simulations. Our methods have been implemented in the R package, funqtl.

Genome-wide association mapping of growth dynamics detects time-specific and general quantitative trait loci

Growth is a complex trait determined by the interplay between many genes, some of which play a role at a specific moment during development whereas others play a more general role. To identify the genetic basis of growth, natural variation in Arabidopsis rosette growth was followed in 324 accessions by a combination of top-view imaging, high-throughput image analysis, modelling of growth dynamics, and end-point fresh weight determination. Genome-wide association (GWA) mapping of the temporal growth data resulted in the detection of time-specific quantitative trait loci (QTLs), whereas mapping of model parameters resulted in another set of QTLs related to the whole growth curve. The positive correlation between projected leaf area (PLA) at different time points during the course of the experiment suggested the existence of general growth factors with a function in multiple developmental stages or with prolonged downstream effects. Many QTLs could not be identified when growth was evaluated only at a single time point. Eleven candidate genes were identified, which were annotated to be involved in the determination of cell number and size, seed germination, embryo development, developmental phase transition, or senescence. For eight of these, a mutant or overexpression phenotype related to growth has been reported, supporting the identification of true positives. In addition, the detection of QTLs without obvious candidate genes implies the annotation of novel functions for underlying genes.

Integrating Image-Based Phenomics and Association Analysis to Dissect the Genetic Architecture of Temporal Salinity Responses in Rice

Salinity affects a significant portion of arable land and is particularly detrimental for irrigated agriculture, which provides one-third of the global food supply. Rice (Oryza sativa), the most important food crop, is salt sensitive. The genetic resources for salt tolerance in rice germplasm exist but are underutilized due to the difficulty in capturing the dynamic nature of physiological responses to salt stress. The genetic basis of these physiological responses is predicted to be polygenic. In an effort to address this challenge, we generated temporal imaging data from 378 diverse rice genotypes across 14 d of 90 mm NaCl stress and developed a statistical model to assess the genetic architecture of dynamic salinity-induced growth responses in rice germplasm. A genomic region on chromosome 3 was strongly associated with the early growth response and was captured using visible range imaging. Fluorescence imaging identified four genomic regions linked to salinity-induced fluorescence responses. A region on chromosome 1 regulates both the fluorescence shift indicative of the longer term ionic stress and the early growth rate decline during salinity stress. We present, to our knowledge, a new approach to capture the dynamic plant responses to its environment and elucidate the genetic basis of these responses using a longitudinal genome-wide association model.

Novel imaging-based phenotyping strategies for dissecting crosstalk in plant development

In an era of genomics, proteomics, and metabolomics a large number of mutants are available. The discovery of their phenotypes is fast becoming the bottleneck of molecular plant physiology. This crisis can be overcome by imaging-based phenotyping, an emerging, rapidly developing and innovative approach integrating plant and computer science. A tremendous amount of digital image data are automatically analysed using techniques of 'machine vision'. This minireview will shed light on the available imaging strategies and discuss standard methods for the automated analysis of images to give the non-bioinformatic reader an idea how the new technology works. A number of successful platforms will be described and the prospects that image-based phenomics may offer for elucidating hormonal cross-talk and molecular growth physiology will be discussed.

Integrating Image-Based Phenomics and Association Analysis to Dissect the Genetic Architecture of Temporal Salinity Responses in Rice

Salinity affects a significant portion of arable land and is particularly detrimental for irrigated agriculture, which provides one-third of the global food supply. Rice (Oryza sativa), the most important food crop, is salt sensitive. The genetic resources for salt tolerance in rice germplasm exist but are underutilized due to the difficulty in capturing the dynamic nature of physiological responses to salt stress. The genetic basis of these physiological responses is predicted to be polygenic. In an effort to address this challenge, we generated temporal imaging data from 378 diverse rice genotypes across 14 d of 90 mm NaCl stress and developed a statistical model to assess the genetic architecture of dynamic salinity-induced growth responses in rice germplasm. A genomic region on chromosome 3 was strongly associated with the early growth response and was captured using visible range imaging. Fluorescence imaging identified four genomic regions linked to salinity-induced fluorescence responses. A region on chromosome 1 regulates both the fluorescence shift indicative of the longer term ionic stress and the early growth rate decline during salinity stress. We present, to our knowledge, a new approach to capture the dynamic plant responses to its environment and elucidate the genetic basis of these responses using a longitudinal genome-wide association model.

Breeding-assisted genomics

The revolution of inexpensive sequencing has ushered in an unprecedented age of genomics. The promise of using this technology to accelerate plant breeding is being realized with a vision of genomics-assisted breeding that will lead to rapid genetic gain for expensive and difficult traits. The reality is now that robust phenotypic data is an increasing limiting resource to complement the current wealth of genomic information. While genomics has been hailed as the discipline to fundamentally change the scope of plant breeding, a more symbiotic relationship is likely to emerge. In the context of developing and evaluating large populations needed for functional genomics, none excel in this area more than plant breeders. While genetic studies have long relied on dedicated, well-structured populations, the resources dedicated to these populations in the context of readily available, inexpensive genotyping is making this philosophy less tractable relative to directly focusing functional genomics on material in breeding programs. Through shifting effort for basic genomic studies from dedicated structured populations, to capturing the entire scope of genetic determinants in breeding lines, we can move towards not only furthering our understanding of functional genomics in plants, but also rapidly improving crops for increased food security, availability and nutrition.

Revealing plant cryptotypes: defining meaningful phenotypes among infinite traits

The plant phenotype is infinite. Plants vary morphologically and molecularly over developmental time, in response to the environment, and genetically. Exhaustive phenotyping remains not only out of reach, but is also the limiting factor to interpreting the wealth of genetic information currently available. Although phenotyping methods are always improving, an impasse remains: even if we could measure the entirety of phenotype, how would we interpret it? We propose the concept of cryptotype to describe latent, multivariate phenotypes that maximize the separation of a priori classes. Whether the infinite points comprising a leaf outline or shape descriptors defining root architecture, statistical methods to discern the quantitative essence of an organism will be required as we approach measuring the totality of phenotype.

From association to prediction: statistical methods for the dissection and selection of complex traits in plants

Quantification of genotype-to-phenotype associations is central to many scientific investigations, yet the ability to obtain consistent results may be thwarted without appropriate statistical analyses. Models for association can consider confounding effects in the materials and complex genetic interactions. Selecting optimal models enables accurate evaluation of associations between marker loci and numerous phenotypes including gene expression. Significant improvements in QTL discovery via association mapping and acceleration of breeding cycles through genomic selection are two successful applications of models using genome-wide markers. Given recent advances in genotyping and phenotyping technologies, further refinement of these approaches is needed to model genetic architecture more accurately and run analyses in a computationally efficient manner, all while accounting for false positives and maximizing statistical power.

Lights, camera, action: high-throughput plant phenotyping is ready for a close-up

Anticipated population growth, shifting demographics, and environmental variability over the next century are expected to threaten global food security. In the face of these challenges, crop yield for food and fuel must be maintained and improved using fewer input resources. In recent years, genetic tools for profiling crop germplasm has benefited from rapid advances in DNA sequencing, and now similar advances are needed to improve the throughput of plant phenotyping. We highlight recent developments in high-throughput plant phenotyping using robotic-assisted imaging platforms and computer vision-assisted analysis tools.

R/qtlcharts: interactive graphics for quantitative trait locus mapping

Every data visualization can be improved with some level of interactivity. Interactive graphics hold particular promise for the exploration of high-dimensional data. R/qtlcharts is an R package to create interactive graphics for experiments to map quantitative trait loci (QTL) (genetic loci that influence quantitative traits). R/qtlcharts serves as a companion to the R/qtl package, providing interactive versions of R/qtl’s static graphs, as well as additional interactive graphs for the exploration of high-dimensional genotype and phenotype data.

Dissecting the phenotypic components of crop plant growth and drought responses based on high-throughput image analysis

Significantly improved crop varieties are urgently needed to feed the rapidly growing human population under changing climates. While genome sequence information and excellent genomic tools are in place for major crop species, the systematic quantification of phenotypic traits or components thereof in a high-throughput fashion remains an enormous challenge. In order to help bridge the genotype to phenotype gap, we developed a comprehensive framework for high-throughput phenotype data analysis in plants, which enables the extraction of an extensive list of phenotypic traits from nondestructive plant imaging over time. As a proof of concept, we investigated the phenotypic components of the drought responses of 18 different barley (Hordeum vulgare) cultivars during vegetative growth. We analyzed dynamic properties of trait expression over growth time based on 54 representative phenotypic features. The data are highly valuable to understand plant development and to further quantify growth and crop performance features. We tested various growth models to predict plant biomass accumulation and identified several relevant parameters that support biological interpretation of plant growth and stress tolerance. These image-based traits and model-derived parameters are promising for subsequent genetic mapping to uncover the genetic basis of complex agronomic traits. Taken together, we anticipate that the analytical framework and analysis results presented here will be useful to advance our views of phenotypic trait components underlying plant development and their responses to environmental cues.

Image-based high-throughput field phenotyping of crop roots

Current plant phenotyping technologies to characterize agriculturally relevant traits have been primarily developed for use in laboratory and/or greenhouse conditions. In the case of root architectural traits, this limits phenotyping efforts, largely, to young plants grown in specialized containers and growth media. Hence, novel approaches are required to characterize mature root systems of older plants grown under actual soil conditions in the field. Imaging methods able to address the challenges associated with characterizing mature root systems are rare due, in part, to the greater complexity of mature root systems, including the larger size, overlap, and diversity of root components. Our imaging solution combines a field-imaging protocol and algorithmic approach to analyze mature root systems grown in the field. Via two case studies, we demonstrate how image analysis can be utilized to estimate localized root traits that reliably capture heritable architectural diversity as well as environmentally induced architectural variation of both monocot and dicot plants. In the first study, we show that our algorithms and traits (including 13 novel traits inaccessible to manual estimation) can differentiate nine maize (Zea mays) genotypes 8 weeks after planting. The second study focuses on a diversity panel of 188 cowpea (Vigna unguiculata) genotypes to identify which traits are sufficient to differentiate genotypes even when comparing plants whose harvesting date differs up to 14 d. Overall, we find that automatically derived traits can increase both the speed and reproducibility of the trait estimation pipeline under field conditions.

Image-based high-throughput field phenotyping of crop roots

Current plant phenotyping technologies to characterize agriculturally relevant traits have been primarily developed for use in laboratory and/or greenhouse conditions. In the case of root architectural traits, this limits phenotyping efforts, largely, to young plants grown in specialized containers and growth media. Hence, novel approaches are required to characterize mature root systems of older plants grown under actual soil conditions in the field. Imaging methods able to address the challenges associated with characterizing mature root systems are rare due, in part, to the greater complexity of mature root systems, including the larger size, overlap, and diversity of root components. Our imaging solution combines a field-imaging protocol and algorithmic approach to analyze mature root systems grown in the field. Via two case studies, we demonstrate how image analysis can be utilized to estimate localized root traits that reliably capture heritable architectural diversity as well as environmentally induced architectural variation of both monocot and dicot plants. In the first study, we show that our algorithms and traits (including 13 novel traits inaccessible to manual estimation) can differentiate nine maize (Zea mays) genotypes 8 weeks after planting. The second study focuses on a diversity panel of 188 cowpea (Vigna unguiculata) genotypes to identify which traits are sufficient to differentiate genotypes even when comparing plants whose harvesting date differs up to 14 d. Overall, we find that automatically derived traits can increase both the speed and reproducibility of the trait estimation pipeline under field conditions.