A two-stage approach for the spatio-temporal analysis of high-throughput phenotyping data

High throughput phenotyping (HTP) platforms and devices are increasingly used for the characterization of growth and developmental processes for large sets of plant genotypes. Such HTP data require challenging statistical analyses in which longitudinal genetic signals need to be estimated against a background of spatio-temporal noise processes. We propose a two-stage approach for the analysis of such longitudinal HTP data. In a first stage, we correct for design features and spatial trends per time point. In a second stage, we focus on the longitudinal modelling of the spatially corrected data, thereby taking advantage of shared longitudinal features between genotypes and plants within genotypes. We propose a flexible hierarchical three-level P-spline growth curve model, with plants/plots nested in genotypes, and genotypes nested in populations. For selection of genotypes in a plant breeding context, we show how to extract new phenotypes, like growth rates, from the estimated genotypic growth curves and their first-order derivatives. We illustrate our approach on HTP data from the PhenoArch greenhouse platform at INRAE Montpellier and the outdoor Field Phenotyping platform at ETH Zurich.

Automated summary

In this paper, a two-stage approach is proposed for the analysis of longitudinal HTP data. It corrects for design features and spatial trends in the first stage and focuses on the longitudinal modeling of the spatially corrected data in the second stage. A flexible hierarchical three-level P-spline growth curve model is used to capture the growth patterns. The approach is illustrated using HTP data from two different platforms.