High-throughput field phenotyping reveals that selection in breeding has affected the phenology and temperature response of wheat in the stem elongation phase

Crop growth and phenology are driven by seasonal changes in environmental variables, with temperature as one important factor. However, knowledge about genotype-specific temperature response and its influence on phenology is limited. Such information is fundamental to improve crop models and adapt selection strategies. We measured the increase in height of 352 European winter wheat varieties in 4 years to quantify phenology, and fitted an asymptotic temperature response model. The model used hourly fluctuations in temperature to parameterize the base temperature (Tmin), the temperature optimum (rmax), and the steepness (lrc) of growth responses. Our results show that higher Tmin and lrc relate to an earlier start and end of stem elongation. A higher rmax relates to an increased final height. Both final height and rmax decreased for varieties originating from the continental east of Europe towards the maritime west. A genome-wide association study (GWAS) indicated a quantitative inheritance and a large degree of independence among loci. Nevertheless, genomic prediction accuracies (GBLUPs) for Tmin and lrc were low (r≤0.32) compared with other traits (r≥0.59). As well as known, major genes related to vernalization, photoperiod, or dwarfing, the GWAS indicated additional, as yet unknown loci that dominate the temperature response.

Frost Damage Index: The Antipode of Growing Degree Days

Abiotic stresses such as heat and frost limit plant growth and productivity. Image-based field phenotyping methods allow quantifying not only plant growth but also plant senescence. Winter crops show senescence caused by cold spells, visible as declines in leaf area. We accurately quantified such declines by monitoring changes in canopy cover based on time-resolved high-resolution imagery in the field. Thirty-six winter wheat genotypes were measured in multiple years. A concept termed "frost damage index" (FDI) was developed that, in analogy to growing degree days, summarizes frost events in a cumulative way. The measured sensitivity of genotypes to the FDI correlated with visual scorings commonly used in breeding to assess winter hardiness. The FDI concept could be adapted to other factors such as drought or heat stress. While commonly not considered in plant growth modeling, integrating such degradation processes may be key to improving the prediction of plant performance for future climate scenarios.

Image-based phenomic prediction can provide valuable decision support in wheat breeding

KEY MESSAGE

Genotype-by-environment interactions of secondary traits based on high-throughput field phenotyping are less complex than those of target traits, allowing for a phenomic selection in unreplicated early generation trials. Traditionally, breeders' selection decisions in early generations are largely based on visual observations in the field. With the advent of affordable genome sequencing and high-throughput phenotyping technologies, enhancing breeders' ratings with such information became attractive. In this research, it is hypothesized that G[Formula: see text]E interactions of secondary traits (i.e., growth dynamics' traits) are less complex than those of related target traits (e.g., yield). Thus, phenomic selection (PS) may allow selecting for genotypes with beneficial response-pattern in a defined population of environments. A set of 45 winter wheat varieties was grown at 5 year-sites and analyzed with linear and factor-analytic (FA) mixed models to estimate G[Formula: see text]E interactions of secondary and target traits. The dynamic development of drone-derived plant height, leaf area and tiller density estimations was used to estimate the timing of key stages, quantities at defined time points and temperature dose-response curve parameters. Most of these secondary traits and grain protein content showed little G[Formula: see text]E interactions. In contrast, the modeling of G[Formula: see text]E for yield required a FA model with two factors. A trained PS model predicted overall yield performance, yield stability and grain protein content with correlations of 0.43, 0.30 and 0.34. While these accuracies are modest and do not outperform well-trained GS models, PS additionally provided insights into the physiological basis of target traits. An ideotype was identified that potentially avoids the negative pleiotropic effects between yield and protein content.

Combining High-Resolution Imaging, Deep Learning, and Dynamic Modeling to Separate Disease and Senescence in Wheat Canopies

Maintenance of sufficiently healthy green leaf area after anthesis is key to ensuring an adequate assimilate supply for grain filling. Tightly regulated age-related physiological senescence and various biotic and abiotic stressors drive overall greenness decay dynamics under field conditions. Besides direct effects on green leaf area in terms of leaf damage, stressors often anticipate or accelerate physiological senescence, which may multiply their negative impact on grain filling. Here, we present an image processing methodology that enables the monitoring of chlorosis and necrosis separately for ears and shoots (stems + leaves) based on deep learning models for semantic segmentation and color properties of vegetation. A vegetation segmentation model was trained using semisynthetic training data generated using image composition and generative adversarial neural networks, which greatly reduced the risk of annotation uncertainties and annotation effort. Application of the models to image time series revealed temporal patterns of greenness decay as well as the relative contributions of chlorosis and necrosis. Image-based estimation of greenness decay dynamics was highly correlated with scoring-based estimations (r ≈ 0.9). Contrasting patterns were observed for plots with different levels of foliar diseases, particularly septoria tritici blotch. Our results suggest that tracking the chlorotic and necrotic fractions separately may enable (a) a separate quantification of the contribution of biotic stress and physiological senescence on overall green leaf area dynamics and (b) investigation of interactions between biotic stress and physiological senescence. The high-throughput nature of our methodology paves the way to conducting genetic studies of disease resistance and tolerance.

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 Zürich.

Outdoor Plant Segmentation With Deep Learning for High-Throughput Field Phenotyping on a Diverse Wheat Dataset

Robust and automated segmentation of leaves and other backgrounds is a core prerequisite of most approaches in high-throughput field phenotyping. So far, the possibilities of deep learning approaches for this purpose have not been explored adequately, partly due to a lack of publicly available, appropriate datasets. This study presents a workflow based on DeepLab v3+ and on a diverse annotated dataset of 190 RGB (350 x 350 pixels) images. Images of winter wheat plants of 76 different genotypes and developmental stages have been acquired throughout multiple years at high resolution in outdoor conditions using nadir view, encompassing a wide range of imaging conditions. Inconsistencies of human annotators in complex images have been quantified, and metadata information of camera settings has been included. The proposed approach achieves an intersection over union (IoU) of 0.77 and 0.90 for plants and soil, respectively. This outperforms the benchmarked machine learning methods which use Support Vector Classifier and/or Random Forrest. The results show that a small but carefully chosen and annotated set of images can provide a good basis for a powerful segmentation pipeline. Compared to earlier methods based on machine learning, the proposed method achieves better performance on the selected dataset in spite of using a deep learning approach with limited data. Increasing the amount of publicly available data with high human agreement on annotations and further development of deep neural network architectures will provide high potential for robust field-based plant segmentation in the near future. This, in turn, will be a cornerstone of data-driven improvement in crop breeding and agricultural practices of global benefit.

Temperature response of wheat affects final height and the timing of stem elongation under field conditions

In wheat, temperature affects the timing and intensity of stem elongation. Genetic variation for this process is therefore important for adaptation. This study investigates the genetic response to temperature fluctuations during stem elongation and its relationship to phenology and height. Canopy height of 315 wheat genotypes (GABI wheat panel) was scanned twice weekly in the field phenotyping platform (FIP) of ETH Zurich using a LIDAR. Temperature response was modelled using linear regressions between stem elongation and mean temperature in each measurement interval. This led to a temperature-responsive (slope) and a temperature-irresponsive (intercept) component. The temperature response was highly heritable (H2=0.81) and positively related to a later start and end of stem elongation as well as final height. Genome-wide association mapping revealed three temperature-responsive and four temperature-irresponsive quantitative trait loci (QTLs). Furthermore, putative candidate genes for temperature-responsive QTLs were frequently related to the flowering pathway in Arabidopsis thaliana, whereas temperature-irresponsive QTLs corresponded to growth and reduced height genes. In combination with Rht and Ppd alleles, these loci, together with the loci for the timing of stem elongation, accounted for 71% of the variability in height. This demonstrates how high-throughput field phenotyping combined with environmental covariates can contribute to a smarter selection of climate-resilient crops.

Temperature response of wheat affects final height and the timing of stem elongation under field conditions

In wheat, temperature affects the timing and intensity of stem elongation. Genetic variation for this process is therefore important for adaptation. This study investigates the genetic response to temperature fluctuations during stem elongation and its relationship to phenology and height. Canopy height of 315 wheat genotypes (GABI wheat panel) was scanned twice weekly in the field phenotyping platform (FIP) of ETH Zurich using a LIDAR. Temperature response was modelled using linear regressions between stem elongation and mean temperature in each measurement interval. This led to a temperature-responsive (slope) and a temperature-irresponsive (intercept) component. The temperature response was highly heritable (H2=0.81) and positively related to a later start and end of stem elongation as well as final height. Genome-wide association mapping revealed three temperature-responsive and four temperature-irresponsive quantitative trait loci (QTLs). Furthermore, putative candidate genes for temperature-responsive QTLs were frequently related to the flowering pathway in Arabidopsis thaliana, whereas temperature-irresponsive QTLs corresponded to growth and reduced height genes. In combination with Rht and Ppd alleles, these loci, together with the loci for the timing of stem elongation, accounted for 71% of the variability in height. This demonstrates how high-throughput field phenotyping combined with environmental covariates can contribute to a smarter selection of climate-resilient crops.

Global Wheat Head Detection 2021: An Improved Dataset for Benchmarking Wheat Head Detection Methods

The Global Wheat Head Detection (GWHD) dataset was created in 2020 and has assembled 193,634 labelled wheat heads from 4700 RGB images acquired from various acquisition platforms and 7 countries/institutions. With an associated competition hosted in Kaggle, GWHD_2020 has successfully attracted attention from both the computer vision and agricultural science communities. From this first experience, a few avenues for improvements have been identified regarding data size, head diversity, and label reliability. To address these issues, the 2020 dataset has been reexamined, relabeled, and complemented by adding 1722 images from 5 additional countries, allowing for 81,553 additional wheat heads. We now release in 2021 a new version of the Global Wheat Head Detection dataset, which is bigger, more diverse, and less noisy than the GWHD_2020 version.

Ultrasonographic diagnosis of closed pedal bone fractures in bovine claws: An ex-vivo study in slaughterhouse specimens

Pedal bone fractures are one of the most common fracture locations in adult cattle and can be diagnosed by radiographs in two planes. Most bovine practitioners do not have access to such X-ray machines, but many use ultrasound units on a daily basis, primarily for reproductive medicine. For this reason, in this double-masked, randomized controlled study, we aimed to investigate the suitability of ultrasonographic examination using a 5 MHz linear transducer for diagnosing closed fractures of the pedal bone in cattle. A total of 54 hindlimb claws from slaughtered cattle were prepared and approximately 50% of the claws were artificially fractured. All claws were ultrasonographically examined twice by two examiners to determine the presence or absence of fractures and their locations. Ultrasound results were confirmed using radiographs of the claws as the reference standard method. All fracture locations as determined by ultrasonography were situated within ±2 mm of the radiographically-determined fracture zone. Ultrasound examination yielded a calculated sensitivity of 93%, a specificity of 91% and an inter-rater reliability of 0.77. The intra-rater reliability for the examiners were 0.96 and 0.88. Examiner experience with ultrasound examination and using ultrasound images for diagnosis could have influenced diagnostic accuracy. We conclude that artificially-created pedal bone fractures in ex-vivo bovine claws can be diagnosed using ultrasonography; similar results are expected in live animals. These results should encourage veterinarians to use ultrasonography for diagnosing pedal bone fractures in cattle.

Repeated Multiview Imaging for Estimating Seedling Tiller Counts of Wheat Genotypes Using Drones

Early generation breeding nurseries with thousands of genotypes in single-row plots are well suited to capitalize on high throughput phenotyping. Nevertheless, methods to monitor the intrinsically hard-to-phenotype early development of wheat are yet rare. We aimed to develop proxy measures for the rate of plant emergence, the number of tillers, and the beginning of stem elongation using drone-based imagery. We used RGB images (ground sampling distance of 3 mm pixel-1) acquired by repeated flights (≥ 2 flights per week) to quantify temporal changes of visible leaf area. To exploit the information contained in the multitude of viewing angles within the RGB images, we processed them to multiview ground cover images showing plant pixel fractions. Based on these images, we trained a support vector machine for the beginning of stem elongation (GS30). Using the GS30 as key point, we subsequently extracted plant and tiller counts using a watershed algorithm and growth modeling, respectively. Our results show that determination coefficients of predictions are moderate for plant count (R 2 = 0.52), but strong for tiller count (R 2 = 0.86) and GS30 (R 2 = 0.77). Heritabilities are superior to manual measurements for plant count and tiller count, but inferior for GS30 measurements. Increasing the selection intensity due to throughput may overcome this limitation. Multiview image traits can replace hand measurements with high efficiency (85-223%). We therefore conclude that multiview images have a high potential to become a standard tool in plant phenomics.

Global Wheat Head Detection (GWHD) Dataset: A Large and Diverse Dataset of High-Resolution RGB-Labelled Images to Develop and Benchmark Wheat Head Detection Methods

The detection of wheat heads in plant images is an important task for estimating pertinent wheat traits including head population density and head characteristics such as health, size, maturity stage, and the presence of awns. Several studies have developed methods for wheat head detection from high-resolution RGB imagery based on machine learning algorithms. However, these methods have generally been calibrated and validated on limited datasets. High variability in observational conditions, genotypic differences, development stages, and head orientation makes wheat head detection a challenge for computer vision. Further, possible blurring due to motion or wind and overlap between heads for dense populations make this task even more complex. Through a joint international collaborative effort, we have built a large, diverse, and well-labelled dataset of wheat images, called the Global Wheat Head Detection (GWHD) dataset. It contains 4700 high-resolution RGB images and 190000 labelled wheat heads collected from several countries around the world at different growth stages with a wide range of genotypes. Guidelines for image acquisition, associating minimum metadata to respect FAIR principles, and consistent head labelling methods are proposed when developing new head detection datasets. The GWHD dataset is publicly available at http://www.global-wheat.com/and aimed at developing and benchmarking methods for wheat head detection.

Exploring genetic dependence of lipase activity to improve the quality of whole-grain wheat

BACKGROUND

Whole-grain wheat flour is facing the quality challenge of lipid rancidity, which decreases its nutritional, sensory, and technological properties. One of the major causes of lipid rancidity is endogenous esterases and lipases. This study evaluated 66 European wheat varieties grown at a single site over three years (2014, 2015, and 2016).

RESULTS

The 66 wheat varieties showed up to threefold variance on esterase and lipase activities. Wheat varieties that are suitable for lipid-stable whole-grain products ('Julius', 'Lona', and 'Banquet') were selected according to their consistently low esterase and lipase activities. The 3-year mean-based broad-sense heritability of esterase and lipase was 0.75 and 0.44 respectively.

CONCLUSIONS

The findings indicate great genetic dependence of both esterase and lipase activities in wheat. The moderate to high heritability brings a new prospect of breeding selection of low-lipase-activity wheat for stable whole-grain products. This result will improve the use of wheat as raw material, benefit cultivation selection, and provide consumers with better quality products. © 2020 Society of Chemical Industry.

Assessment of Multi-Image Unmanned Aerial Vehicle Based High-Throughput Field Phenotyping of Canopy Temperature

Canopy temperature (CT) has been related to water-use and yield formation in crops. However, constantly (e.g., sun illumination angle, ambient temperature) as well as rapidly (e.g., clouds) changing environmental conditions make it difficult to compare measurements taken even at short time intervals. This poses a great challenge for high-throughput field phenotyping (HTFP). The aim of this study was to i) set up a workflow for unmanned aerial vehicles (UAV) based HTFP of CT, ii) investigate different data processing procedures to combine information from multiple images into orthomosaics, iii) investigate the repeatability of the resulting CT by means of heritability, and iv) investigate the optimal timing for thermography measurements. Additionally, the approach was v) compared with other methods for HTFP of CT. The study was carried out in a winter wheat field trial with 354 genotypes planted in two replications in a temperate climate, where a UAV captured CT in a time series of 24 flights during 6 weeks of the grain-filling phase. Custom-made thermal ground control points enabled accurate georeferencing of the data. The generated thermal orthomosaics had a high spatial accuracy (mean ground sampling distance of 5.03 cm/pixel) and position accuracy [mean root-mean-square deviation (RMSE) = 4.79 cm] over all time points. An analysis on the impact of the measurement geometry revealed a gradient of apparent CT in parallel to the principle plane of the sun and a hotspot around nadir. Averaging information from all available images (and all measurement geometries) for an area of interest provided the best results by means of heritability. Correcting for spatial in-field heterogeneity as well as slight environmental changes during the measurements were performed with the R package SpATS. CT heritability ranged from 0.36 to 0.74. Highest heritability values were found in the early afternoon. Since senescence was found to influence the results, it is recommended to measure CT in wheat after flowering and before the onset of senescence. Overall, low-altitude and high-resolution remote sensing proved suitable to assess the CT of crop genotypes in a large number of small field plots as is required in crop breeding and variety testing experiments.

Spectral Vegetation Indices to Track Senescence Dynamics in Diverse Wheat Germplasm

The ability of a genotype to stay green affects the primary target traits grain yield (GY) and grain protein concentration (GPC) in wheat. High throughput methods to assess senescence dynamics in large field trials will allow for (i) indirect selection in early breeding generations, when yield cannot yet be accurately determined and (ii) mapping of the genomic regions controlling the trait. The aim of this study was to develop a robust method to assess senescence based on hyperspectral canopy reflectance. Measurements were taken in three years throughout the grain filling phase on >300 winter wheat varieties in the spectral range from 350 to 2500 nm using a spectroradiometer. We compared the potential of spectral indices (SI) and full-spectrum models to infer visually observed senescence dynamics from repeated reflectance measurements. Parameters describing the dynamics of senescence were used to predict GY and GPC and a feature selection algorithm was used to identify the most predictive features. The three-band plant senescence reflectance index (PSRI) approximated the visually observed senescence dynamics best, whereas full-spectrum models suffered from a strong year-specificity. Feature selection identified visual scorings as most predictive for GY, but also PSRI ranked among the most predictive features while adding additional spectral features had little effect. Visually scored delayed senescence was positively correlated with GY ranging from r = 0.173 in 2018 to r = 0.365 in 2016. It appears that visual scoring remains the gold standard to quantify leaf senescence in moderately large trials. However, using appropriate phenotyping platforms, the proposed index-based parameterization of the canopy reflectance dynamics offers the critical advantage of upscaling to very large breeding trials.

In-Field Detection and Quantification of Septoria Tritici Blotch in Diverse Wheat Germplasm Using Spectral-Temporal Features

Hyperspectral remote sensing holds the potential to detect and quantify crop diseases in a rapid and non-invasive manner. Such tools could greatly benefit resistance breeding, but their adoption is hampered by i) a lack of specificity to disease-related effects and ii) insufficient robustness to variation in reflectance caused by genotypic diversity and varying environmental conditions, which are fundamental elements of resistance breeding. We hypothesized that relying exclusively on temporal changes in canopy reflectance during pathogenesis may allow to specifically detect and quantify crop diseases while minimizing the confounding effects of genotype and environment. To test this hypothesis, we collected time-resolved canopy hyperspectral reflectance data for 18 diverse genotypes on infected and disease-free plots and engineered spectral-temporal features representing this hypothesis. Our results confirm the lack of specificity and robustness of disease assessments based on reflectance spectra at individual time points. We show that changes in spectral reflectance over time are indicative of the presence and severity of Septoria tritici blotch (STB) infections. Furthermore, the proposed time-integrated approach facilitated the delineation of disease from physiological senescence, which is pivotal for efficient selection of STB-resistant material under field conditions. A validation of models based on spectral-temporal features on a diverse panel of 330 wheat genotypes offered evidence for the robustness of the proposed method. This study demonstrates the potential of time-resolved canopy reflectance measurements for robust assessments of foliar diseases in the context of resistance breeding.

Precision Phenotyping Reveals Novel Loci for Quantitative Resistance to Septoria Tritici Blotch

Accurate, high-throughput phenotyping for quantitative traits is a limiting factor for progress in plant breeding. We developed an automated image analysis to measure quantitative resistance to septoria tritici blotch (STB), a globally important wheat disease, enabling identification of small chromosome intervals containing plausible candidate genes for STB resistance. 335 winter wheat cultivars were included in a replicated field experiment that experienced natural epidemic development by a highly diverse but fungicide-resistant pathogen population. More than 5.4 million automatically generated phenotypes were associated with 13,648 SNP markers to perform the GWAS. We identified 26 chromosome intervals explaining 1.9-10.6% of the variance associated with four independent resistance traits. Sixteen of the intervals overlapped with known STB resistance intervals, suggesting that our phenotyping approach can identify simultaneously (i.e., in a single experiment) many previously defined STB resistance intervals. Seventeen of the intervals were less than 5 Mbp in size and encoded only 173 genes, including many genes associated with disease resistance. Five intervals contained four or fewer genes, providing high priority targets for functional validation. Ten chromosome intervals were not previously associated with STB resistance, perhaps representing resistance to pathogen strains that had not been tested in earlier experiments. The SNP markers associated with these chromosome intervals can be used to recombine different forms of quantitative STB resistance that are likely to be more durable than pyramids of major resistance genes. Our experiment illustrates how high-throughput automated phenotyping can accelerate breeding for quantitative disease resistance.

PhenoFly Planning Tool: flight planning for high-resolution optical remote sensing with unmanned areal systems

BACKGROUND

Driven by a huge improvement in automation, unmanned areal systems (UAS) are increasingly used for field observations and high-throughput phenotyping. Today, the bottleneck does not lie in the ability to fly a drone anymore, but rather in the appropriate flight planning to capture images with sufficient quality. Proper flight preparation for photography with digital frame cameras should include relevant concepts such as view, sharpness and exposure calculations. Additionally, if mapping areas with UASs, one has to consider concepts related to ground control points (GCPs), viewing geometry and way-point flights. Unfortunately, non of the available flight planning tools covers all these aspects.

RESULTS

We give an overview of concepts related to flight preparation, present the newly developed open source software PhenoFly Planning Tool, and evaluate other recent flight planning tools. We find that current flight planning and mapping tools strongly focus on vendor-specific solutions and mostly ignore basic photographic properties-our comparison shows, for example, that only two out of thirteen evaluated tools consider motion blur restrictions, and none of them depth of field limits. In contrast, PhenoFly Planning Tool enhances recent sophisticated UAS and autopilot systems with an optical remote sensing workflow that respects photographic concepts. The tool can assist in selecting the right equipment for your needs, experimenting with different flight settings to test the performance of the resulting imagery, preparing the field and GCP setup, and generating a flight path that can be exported as waypoints to be uploaded to an UAS.

CONCLUSION

By considering the introduced concepts, uncertainty in UAS-based remote sensing and high-throughput phenotyping may be considerably reduced. The presented software PhenoFly Planning Tool (https://shiny.usys.ethz.ch/PhenoFlyPlanningTool) helps users to comprehend and apply these concepts.

Hyperspectral Canopy Sensing of Wheat Septoria Tritici Blotch Disease

Producing quantitative and reliable measures of crop disease is essential for resistance breeding, but is challenging and time consuming using traditional phenotyping methods. Hyperspectral remote sensing has shown potential for the detection of plant diseases, but its utility for phenotyping large and diverse populations of plants under field conditions requires further evaluation. In this study, we collected canopy hyperspectral data from 335 wheat varieties using a spectroradiometer, and we investigated the use of canopy reflectance for detecting the Septoria tritici blotch (STB) disease and for quantifying the severity of infection. Canopy- and leaf-level infection metrics of STB based on traditional visual assessments and automated analyses of leaf images were used as ground truth data. Results showed (i) that canopy reflectance and the selected spectral indices show promise for quantifying STB infections, and (ii) that the normalized difference water index (NDWI) showed the best performance in detecting STB compared to other spectral indices. Moreover, partial least squares (PLS) regression models allowed for an improvement in the prediction of STB metrics. The PLS discriminant analysis (PLSDA) model calibrated based on the spectral data of four reference varieties was able to discriminate between the diseased and healthy canopies among the 335 varieties with an accuracy of 93% (Kappa = 0.60). Finally, the PLSDA model predictions allowed for the identification of wheat genotypes that are potentially more susceptible to STB, which was confirmed by the STB visual assessment. This study demonstrates the great potential of using canopy hyperspectral remote sensing to improve foliar disease assessment and to facilitate plant breeding for disease resistance.

Ranking Quantitative Resistance to Septoria tritici Blotch in Elite Wheat Cultivars Using Automated Image Analysis

Quantitative resistance is likely to be more durable than major gene resistance for controlling Septoria tritici blotch (STB) on wheat. Earlier studies hypothesized that resistance affecting the degree of host damage, as measured by the percentage of leaf area covered by STB lesions, is distinct from resistance that affects pathogen reproduction, as measured by the density of pycnidia produced within lesions. We tested this hypothesis using a collection of 335 elite European winter wheat cultivars that was naturally infected by a diverse population of Zymoseptoria tritici in a replicated field experiment. We used automated image analysis of 21,420 scanned wheat leaves to obtain quantitative measures of conditional STB intensity that were precise, objective, and reproducible. These measures allowed us to explicitly separate resistance affecting host damage from resistance affecting pathogen reproduction, enabling us to confirm that these resistance traits are largely independent. The cultivar rankings based on host damage were different from the rankings based on pathogen reproduction, indicating that the two forms of resistance should be considered separately in breeding programs aiming to increase STB resistance. We hypothesize that these different forms of resistance are under separate genetic control, enabling them to be recombined to form new cultivars that are highly resistant to STB. We found a significant correlation between rankings based on automated image analysis and rankings based on traditional visual scoring, suggesting that image analysis can complement conventional measurements of STB resistance, based largely on host damage, while enabling a much more precise measure of pathogen reproduction. We showed that measures of pathogen reproduction early in the growing season were the best predictors of host damage late in the growing season, illustrating the importance of breeding for resistance that reduces pathogen reproduction in order to minimize yield losses caused by STB. These data can already be used by breeding programs to choose wheat cultivars that are broadly resistant to naturally diverse Z. tritici populations according to the different classes of resistance.

An image analysis pipeline for automated classification of imaging light conditions and for quantification of wheat canopy cover time series in field phenotyping

BACKGROUND

Robust segmentation of canopy cover (CC) from large amounts of images taken under different illumination/light conditions in the field is essential for high throughput field phenotyping (HTFP). We attempted to address this challenge by evaluating different vegetation indices and segmentation methods for analyzing images taken at varying illuminations throughout the early growth phase of wheat in the field. 40,000 images taken on 350 wheat genotypes in two consecutive years were assessed for this purpose.

RESULTS

We proposed an image analysis pipeline that allowed for image segmentation using automated thresholding and machine learning based classification methods and for global quality control of the resulting CC time series. This pipeline enabled accurate classification of imaging light conditions into two illumination scenarios, i.e. high light-contrast (HLC) and low light-contrast (LLC), in a series of continuously collected images by employing a support vector machine (SVM) model. Accordingly, the scenario-specific pixel-based classification models employing decision tree and SVM algorithms were able to outperform the automated thresholding methods, as well as improved the segmentation accuracy compared to general models that did not discriminate illumination differences.

CONCLUSIONS

The three-band vegetation difference index (NDI3) was enhanced for segmentation by incorporating the HSV-V and the CIE Lab-a color components, i.e. the product images NDI3*V and NDI3*a. Field illumination scenarios can be successfully identified by the proposed image analysis pipeline, and the illumination-specific image segmentation can improve the quantification of CC development. The integrated image analysis pipeline proposed in this study provides great potential for automatically delivering robust data in HTFP.

RADIX: rhizoslide platform allowing high throughput digital image analysis of root system expansion

BACKGROUND

Phenotyping of genotype-by-environment interactions in the root-zone is of major importance for crop improvement as the spatial distribution of a plant's root system is crucial for a plant to access water and nutrient resources of the soil. However, so far it is unclear to what extent genetic variations in root system responses to spatially varying soil resources can be utilized for breeding applications. Among others, one limiting factor is the absence of phenotyping platforms allowing the analysis of such interactions.

RESULTS

We developed a system that is able to (a) monitor root and shoot growth synchronously, (b) investigate their dynamic responses and (c) analyse the effect of heterogeneous N distribution to parts of the root system in a split-nutrient setup with a throughput (200 individual maize plants at once) sufficient for mapping of quantitative trait loci or for screens of multiple environmental factors. In a test trial, 24 maize genotypes were grown under split nitrogen conditions and the response of shoot and root growth was investigated. An almost double elongation rate of crown and lateral roots was observed under high N for all genotypes. The intensity of genotype-specific responses varied strongly. For example, elongation of crown roots differed almost two times between the fastest and slowest growing genotype. A stronger selective root placement in the high-N compartment was related to an increased shoot development indicating that early vigour might be related to a more intense foraging behaviour.

CONCLUSION

To our knowledge, RADIX is the only system currently existing which allows studying the differential response of crown roots to split-nutrient application to quantify foraging behaviour in genome mapping or selection experiments. In doing so, changes in root and shoot development and the connection to plant performance can be investigated.

Application of the Prunus spp. Cyanide Seed Defense System onto Wheat: Reduced Insect Feeding and Field Growth Tests

Many crops are ill-protected against insect pests during storage. To protect cereal grains from herbivores during storage, pesticides are often applied. While pesticides have an undoubtable functionality, increasing concerns are arising about their application. In the present study, we investigated a bioinspired cyanogenic grain coating with amygdalin as cyanogenic precursor mimicking the feeding-triggered release of hydrogen cyanide (HCN) found for example in bitter almonds. The multilayer coating consisted of biodegradable polylactic acid with individual layers containing amygdalin or β-glucosidase which is capable of degrading amygdalin to HCN. This reaction occurred only when the layers were ruptured, e.g., by a herbivore attack. Upon feeding coated cyanogenic wheat grains to Tenebrio molitor (mealworm beetle), Rhizopertha dominica (lesser grain borer), and Plodia interpunctella (Indianmeal moth), their reproduction as well as consumption rate were significantly reduced, whereas germination ability increased compared to noncoated grains. In field experiments, we observed an initial growth delay compared to uncoated grains which became negligible at later growth stages. The here shown strategy to artificially apply a naturally occurring defense mechanisms could be expanded to other crops than wheat and has the potential to replace certain pesticides with the benefit of complete biodegradability and increased safety during storage.

The ETH field phenotyping platform FIP: a cable-suspended multi-sensor system

Crop phenotyping is a major bottleneck in current plant research. Field-based high-throughput phenotyping platforms are an important prerequisite to advance crop breeding. We developed a cable-suspended field phenotyping platform covering an area of ~1ha. The system operates from 2 to 5m above the canopy, enabling a high image resolution. It can carry payloads of up to 12kg and can be operated under adverse weather conditions. This ensures regular measurements throughout the growing period even during cold, windy and moist conditions. Multiple sensors capture the reflectance spectrum, temperature, height or architecture of the canopy. Monitoring from early development to maturity at high temporal resolution allows the determination of dynamic traits and their correlation to environmental conditions throughout the entire season. We demonstrate the capabilities of the system with respect to monitoring canopy cover, canopy height and traits related to thermal and multi-spectral imaging by selected examples from winter wheat, maize and soybean. The system is discussed in the context of other, recently established field phenotyping approaches; such as ground-operating or aerial vehicles, which impose traffic on the field or require a higher distance to the canopy.

High-resolution quantification of root dynamics in split-nutrient rhizoslides reveals rapid and strong proliferation of maize roots in response to local high nitrogen

The plant's root system is highly plastic, and can respond to environmental stimuli such as high nitrogen (N) in patches. A root may respond to an N patch by selective placement of new lateral roots, and therewith increases root N uptake. This may be a desirable trait in breeding programmes, since it decreases NO3(-) leaching and N2O emission. Roots of maize (Zea mays L.) were grown without N in split-nutrient rhizoslides. One side of the slides was exposed to high N after 15 d of root development, and root elongation was measured for another 15 d, described in a time course model and parameterized. The elongation rates of crown axile roots on the N-treated side of the plant followed a logistic increase to a maximum of 5.3cm d(-1); 95% of the maximum were reached within 4 d. At the same time, on the untreated side, axile root elongation dropped linearly to 1.2cm d(-1) within 6.4 d and stayed constant thereafter. Twice as many lateral roots were formed on the crown axis on the N side compared to the untreated side. Most strikingly, the elongation rates of laterals of the N side increased linearly with most of the roots reaching an asymptote ~8 d after start of the N treatment. By contrast, laterals on the side without N did not show any detectable elongation beyond the first day after their emergence. We conclude that split-nutrient rhizoslides have great potential to improve our knowledge about nitrogen responsiveness and selection for contrasting genotypes.

Image based phenotyping during winter: a powerful tool to assess wheat genetic variation in growth response to temperature

Having a strong effect on plant growth, temperature adaption has become a major breeding aim. Due to a lack of efficient methods, we developed an image-based approach to characterise genotypes for their temperature behaviour in the field. Twenty-nine winter wheat (Triticum aestivum L.) genotypes were continuously monitored at 3-day intervals on a plot basis during early growth from November to March using a modified digital camera. Canopy cover (CC) was determined by segmentation of leaves in calibrated images. Relative growth rates (RGR) of CC were then calculated for each measurement interval and related to the respective temperature. Also, classical traits used in plant breeding were assessed. Measurements of CC at single dates were highly repeatable with respect to genotype. For the tested range of temperatures (0-7°C), a linear relation between RGR and temperature was observed. Genotypes differed for base temperature and increase in RGR with rising temperature, these two traits showing a strong positive correlation with each other but being independent of CC at a single date. Our simple approach is suitable to screen large populations for differences in growth response to environmental stimuli. Furthermore, the derived parameters reveal additional information that cannot be assessed by usual measurements of static size.

Plant phenotyping: from bean weighing to image analysis

Plant phenotyping refers to a quantitative description of the plant's anatomical, ontogenetical, physiological and biochemical properties. Today, rapid developments are taking place in the field of non-destructive, image-analysis -based phenotyping that allow for a characterization of plant traits in high-throughput. During the last decade, 'the field of image-based phenotyping has broadened its focus from the initial characterization of single-plant traits in controlled conditions towards 'real-life' applications of robust field techniques in plant plots and canopies. An important component of successful phenotyping approaches is the holistic characterization of plant performance that can be achieved with several methodologies, ranging from multispectral image analyses via thermographical analyses to growth measurements, also taking root phenotypes into account.

Remote, aerial phenotyping of maize traits with a mobile multi-sensor approach

BACKGROUND

Field-based high throughput phenotyping is a bottleneck for crop breeding research. We present a novel method for repeated remote phenotyping of maize genotypes using the Zeppelin NT aircraft as an experimental sensor platform. The system has the advantage of a low altitude and cruising speed compared to many drones or airplanes, thus enhancing image resolution while reducing blurring effects. Additionally there was no restriction in sensor weight. Using the platform, red, green and blue colour space (RGB), normalized difference vegetation index (NDVI) and thermal images were acquired throughout the growing season and compared with traits measured on the ground. Ground control points were used to co-register the images and to overlay them with a plot map.

RESULTS

NDVI images were better suited than RGB images to segment plants from soil background leading to two separate traits: the canopy cover (CC) and its NDVI value (NDVIPlant). Remotely sensed CC correlated well with plant density, early vigour, leaf size, and radiation interception. NDVIPlant was less well related to ground truth data. However, it related well to the vigour rating, leaf area index (LAI) and leaf biomass around flowering and to very late senescence rating. Unexpectedly, NDVIPlant correlated negatively with chlorophyll meter measurements. This could be explained, at least partially, by methodical differences between the used devices and effects imposed by the population structure. Thermal images revealed information about the combination of radiation interception, early vigour, biomass, plant height and LAI. Based on repeatability values, we consider two row plots as best choice to balance between precision and available field space. However, for thermography, more than two rows improve the precision.

CONCLUSIONS

We made important steps towards automated processing of remotely sensed data, and demonstrated the value of several procedural steps, facilitating the application in plant genetics and breeding. Important developments are: the ability to monitor throughout the season, robust image segmentation and the identification of individual plots in images from different sensor types at different dates. Remaining bottlenecks are: sufficient ground resolution, particularly for thermal imaging, as well as a deeper understanding of the relatedness of remotely sensed data and basic crop characteristics.

Rhizoslides: paper-based growth system for non-destructive, high throughput phenotyping of root development by means of image analysis

BACKGROUND

A quantitative characterization of root system architecture is currently being attempted for various reasons. Non-destructive, rapid analyses of root system architecture are difficult to perform due to the hidden nature of the root. Hence, improved methods to measure root architecture are necessary to support knowledge-based plant breeding and to analyse root growth responses to environmental changes. Here, we report on the development of a novel method to reveal growth and architecture of maize root systems.

RESULTS

The method is based on the cultivation of different root types within several layers of two-dimensional, large (50 × 60 cm) plates (rhizoslides). A central plexiglass screen stabilizes the system and is covered on both sides with germination paper providing water and nutrients for the developing root, followed by a transparent cover foil to prevent the roots from falling dry and to stabilize the system. The embryonic roots grow hidden between a Plexiglas surface and paper, whereas crown roots grow visible between paper and the transparent cover. Long cultivation with good image quality up to 20 days (four fully developed leaves) was enhanced by suppressing fungi with a fungicide. Based on hyperspectral microscopy imaging, the quality of different germination papers was tested and three provided sufficient contrast to distinguish between roots and background (segmentation). Illumination, image acquisition and segmentation were optimised to facilitate efficient root image analysis. Several software packages were evaluated with regard to their precision and the time investment needed to measure root system architecture. The software 'Smart Root' allowed precise evaluation of root development but needed substantial user interference. 'GiaRoots' provided the best segmentation method for batch processing in combination with a good analysis of global root characteristics but overestimated root length due to thinning artefacts. 'WhinRhizo' offered the most rapid and precise evaluation of root lengths in diameter classes, but had weaknesses with respect to image segmentation and analysis of root system architecture.

CONCLUSION

A new technique has been established for non-destructive root growth studies and quantification of architectural traits beyond seedlings stages. However, automation of the scanning process and appropriate software remains the bottleneck for high throughput analysis.

Can we improve heterosis for root growth of maize by selecting parental inbred lines with different temperature behaviour?

Tolerance to high and low temperature is an important breeding aim for Central and Northern Europe, where temperature fluctuations are predicted to increase. However, the extent to which genotypes differ in their response to the whole range of possible temperatures is not well understood. We tested the hypothesis that the combination of maize (Zea mays L.) inbred lines with differing temperature optima for root growth would lead to superior hybrids. This hypothesis is based on the concept of 'marginal overdominance' in which the hybrid expresses higher relative fitness than its parents, summed over all situations. The elongation rates of axile and lateral roots of the reciprocal cross between two flint and two dent inbred lines were assessed at temperatures between 15°C and 40°C. Indeed, the cross between UH005 and UH250 with lateral root growth temperature optima at 34°C and 28°C, respectively, resulted in intermediate hybrids. At temperatures below and above 31°C, the hybrids' root growth was comparable to the better parent, respectively, thereby increasing temperature tolerance of the hybrid compared with its parents. The implications of and reasons for this heterosis effect are discussed in the context of breeding for abiotic stress tolerance and of putatively underlying molecular mechanisms. This finding paves the way for more detailed investigations of this phenomenon in future studies.

QTLs for the elongation of axile and lateral roots of maize in response to low water potential

Changes in root architecture and the maintenance of root growth in drying soil are key traits for the adaptation of maize (Zea mays L.) to drought environments. The goal of this study was to map quantitative trait loci (QTLs) for root growth and its response to dehydration in a population of 208 recombinant inbred lines from the International Maize and Wheat Improvement Center (CIMMYT). The parents, Ac7643 and Ac7729/TZSRW, are known to be drought-tolerant and drought-sensitive, respectively. Roots were grown in pouches under well-watered conditions or at low water potential induced by the osmolyte polyethylene glycol (PEG 8000). Axile root length (L (Ax)) increased linearly, while lateral root length (L (Lat)) increased exponentially over time. Thirteen QTLs were identified for six seedling traits: elongation rates of axile roots (ER(Ax)), the rate constant of lateral root elongation (k (Lat)), the final respective lengths (L (Ax) and L (Lat)), and the ratios k (Lat)/ER(Ax) and L (Lat)/L (Ax.) While QTLs for lateral root traits were constitutively expressed, most QTLs for axile root traits responded to water stress. For axile roots, common QTLs existed for ER(Ax) and L (Ax). Quantitative trait loci for the elongation rates of axile roots responded more clearly to water stress compared to root length. Two major QTLs were detected: a QTL for general vigor in bin 2.02, affecting most of the traits, and a QTL for the constitutive increase in k (Lat) and k (Lat)/ER(Ax) in bins 6.04-6.05. The latter co-located with a major QTL for the anthesis-silking interval (ASI) reported in published field experiments, suggesting an involvement of root morphology in drought tolerance. Rapid seedling tests are feasible for elucidating the genetic response of root growth to low water potential. Some loci may even have pleiotropic effects on yield-related traits under drought stress.

Mapping of QTLs for lateral and axile root growth of tropical maize

Maize genotypes may adapt to dry environments by avoiding desiccation by means of a deeper root system or by maintaining growth and water extraction at low water potentials. The aim of this study was to determine the quantitative genetic control of root growth and root morphology in a population of 236 recombinant inbred lines (RILs) from the cross between CML444 (high-yielding) x SC-Malawi (low-yielding), which segregates for the response to drought stress at flowering. The RILs and the parental lines were grown on blotting paper in growth pouches until the two-leaf stage under non-stressed conditions; the parents were additionally exposed to desiccation stress induced by polyethylene glycol with a molecular weight of 8000 Dalton (PEG-8000). The lengths of axile and lateral roots were measured non-destructively at 2, 5, 7 and 9 days after germination, by scanning with an A4 scanner followed by digital image analysis. CML444 had a lower rate constant of lateral root elongation (k(Lat)) than SC-Malawi, but the two genotypes did not differ in their response to desiccation. QTLs affecting root vigor, as depicted by increments in k(Lat), the elongation rate of axile roots (ER(Ax)) and the number of axile roots (No(Ax)) were identified in bins 2.04 and 2.05. QTLs for No(Ax) and ER(Ax) collocated with QTLs for yield parameters in bins 1.03-1.04 and 7.03-04. The correspondence of QTLs for axile root traits in bins 1.02-1.03 and 1.08 and QTLs for lateral roots traits in bins 2.04-2.07 in several mapping populations suggests the presence of genes controlling root growth in a wide range of genetic backgrounds.

QTL controlling root and shoot traits of maize seedlings under cold stress

The improvement of early vigour is crucial for the adaptation of maize (Zea mays L.) to the climatic conditions of central Europe and the northern Mediterranean, where early sowing is an important strategy for avoiding the effect of summer drought. The objectives of this study were to identify quantitative trait loci (QTL) controlling cold-related traits and to investigate the relationships among them. A set of 168 F2:4 families of the Lo964 x Lo1016 cross was grown in a sand-vermiculite substrate at 15/13 degrees C (day/night) until the one-leaf stage. Twenty QTL were identified for the four shoot and two seed traits examined. Analysis of root weight and digital measurements of the length and diameter of primary and seminal roots led to the identification of 40 QTL. The operating efficiency of photosystem II (PhiPSII) was related to seedling dry weight at both the phenotypic and genetic level (r = 0.46, two matching loci, respectively) but was not related to root traits. Cluster analysis and QTL association revealed that the different root traits were largely independently inherited and that root lengths and diameters were mostly negatively correlated. The major QTL for root traits detected in an earlier study in hydroponics were confirmed in this study. The length of the primary lateral roots was negatively associated with the germination index (r = -0.38, two matching loci). Therefore, we found a large number of independently inherited loci suitable for the improvement of early seedling growth through better seed vigour and/or a higher rate of photosynthesis.