Functional mapping of N deficiency‐induced response in wheat yield‐component traits by implementing high‐throughput phenotyping

From Images to Loci: Applying 3D Deep Learning to Enable Multivariate and Multitemporal Digital Phenotyping and Mapping the Genetics Underlying Nitrogen Use Efficiency in Wheat

The selection and promotion of high-yielding and nitrogen-efficient wheat varieties can reduce nitrogen fertilizer application while ensuring wheat yield and quality and contribute to the sustainable development of agriculture; thus, the mining and localization of nitrogen use efficiency (NUE) genes is particularly important, but the localization of NUE genes requires a large amount of phenotypic data support. In view of this, we propose the use of low-altitude aerial photography to acquire field images at a large scale, generate 3-dimensional (3D) point clouds and multispectral images of wheat plots, propose a wheat 3D plot segmentation dataset, quantify the plot canopy height via combination with PointNet++, and generate 4 nitrogen utilization-related vegetation indices via index calculations. Six height-related and 24 vegetation-index-related dynamic digital phenotypes were extracted from the digital phenotypes collected at different time points and fitted to generate dynamic curves. We applied height-derived dynamic numerical phenotypes to genome-wide association studies of 160 wheat cultivars (660,000 single-nucleotide polymorphisms) and found that we were able to locate reliable loci associated with height and NUE, some of which were consistent with published studies. Finally, dynamic phenotypes derived from plant indices can also be applied to genome-wide association studies and ultimately locate NUE- and growth-related loci. In conclusion, we believe that our work demonstrates valuable advances in 3D digital dynamic phenotyping for locating genes for NUE in wheat and provides breeders with accurate phenotypic data for the selection and breeding of nitrogen-efficient wheat varieties.

Converging functional phenotyping with systems mapping to illuminate the genotype-phenotype associations

Illuminating the phenotype-genotype black box under complex traits is an ambitious goal for researchers. The generation of temporally or spatially phenotypic data today has far outpaced its interpretation, due to their highly dynamic nature depending on the environment and developmental stages. Here, we propose an integrated enviro-pheno-geno functional approach to pinpoint the major challenges of decomposing physiological traits. The strategy first features high-throughput functional physiological phenotyping (FPP) to efficiently acquire phenotypic and environmental data. It then features functional mapping (FM) and the extended systems mapping (SM) to tackle trait dynamics. FM, by modeling traits as continuous functions, can increase the power and efficiency in dissecting the spatiotemporal effects of QTLs. SM could enable reconstruction of a genotype-phenotype map from developmental pathways. We present a recent case study that combines FPP and SM to dissect complex physiological traits. This integrated approach will be an important engine to drive the translation of phenomic big data into genetic gain.

High-throughput phenotyping using hyperspectral indicators supports the genetic dissection of yield in durum wheat grown under heat and drought stress

High-throughput phenotyping (HTP) provides new opportunities for efficiently dissecting the genetic basis of drought-adaptive traits, which is essential in current wheat breeding programs. The combined use of HTP and genome-wide association (GWAS) approaches has been useful in the assessment of complex traits such as yield, under field stress conditions including heat and drought. The aim of this study was to identify molecular markers associated with yield (YLD) in elite durum wheat that could be explained using hyperspectral indices (HSIs) under drought field conditions in Mediterranean environments in Southern Spain. The HSIs were obtained from hyperspectral imagery collected during the pre-anthesis and anthesis crop stages using an airborne platform. A panel of 536 durum wheat lines were genotyped by sequencing (GBS, DArTseq) to determine population structure, revealing a lack of genetic structure in the breeding germplasm. The material was phenotyped for YLD and 19 HSIs for six growing seasons under drought field conditions at two locations in Andalusia, in southern Spain. GWAS analysis identified 740 significant marker-trait associations (MTAs) across all the durum wheat chromosomes, several of which were common for YLD and the HSIs, and can potentially be integrated into breeding programs. Candidate gene (CG) analysis uncovered genes related to important plant processes such as photosynthesis, regulatory biological processes, and plant abiotic stress tolerance. These results are novel in that they combine high-resolution hyperspectral imaging at the field scale with GWAS analysis in wheat. They also support the use of HSIs as useful tools for identifying chromosomal regions related to the heat and drought stress response in wheat, and pave the way for the integration of field HTP in wheat breeding programs.

Genomic regions influencing the hyperspectral phenome of deoxynivalenol infected wheat

The quantitative nature of fusarium head blight (FHB) resistance requires further exploration of the wheat genome to identify regions conferring resistance. In this study, we explored the application of hyperspectral imaging of Fusarium-infected wheat kernels and identified regions of the wheat genome contributing significantly to the accumulation of Deoxynivalenol (DON) mycotoxin. Strong correlations were identified between hyperspectral reflectance values for 204 wavebands in the 397-673 nm range and DON mycotoxin. Dimensionality reduction using principal components was performed for all 204 wavebands and 38 sliding windows across the range of wavebands. The first principal component (PC1) of all 204 wavebands explained 70% of the total variation in waveband reflectance values and was highly correlated with DON mycotoxin. PC1 was used as a phenotype in a genome wide association study and a large effect QTL on chromosome 2D was identified for PC1 of all wavebands as well as nearly all 38 sliding windows. The allele contributing variation in PC1 values also led to a substantial reduction in DON. The 2D polymorphism affecting DON levels localized to the exon of TraesCS2D02G524600 which is upregulated in wheat spike and rachis tissues during FHB infection. This work demonstrates the value of hyperspectral imaging as a correlated trait for investigating the genetic basis of resistance and developing wheat varieties with enhanced resistance to FHB.

Genome-wide association study of agronomic traits related to nitrogen use efficiency in Henan wheat

BACKGROUND

Nitrogen use efficiency (NUE) is closely related to crop yield and nitrogen fertilizer application rate. Although NUE is susceptible to environments, quantitative trait nucleotides (QTNs) for NUE in wheat germplasm populations have been rarely reported in genome-wide associated study.

RESULTS

In this study, 244 wheat accessions were phenotyped by three NUE-related traits in three environments and genotyped by 203,224 SNPs. All the phenotypes for each trait were used to associate with all the genotypes of these SNP markers for identifying QTNs and QTN-by-environment interactions via 3VmrMLM. Among 279 QTNs and one QTN-by-environment interaction for low nitrogen tolerance, 33 were stably identified, especially, one large QTN (r2 > 10%), qPHR3A.2, was newly identified for plant height ratio in one environment and multi-environment joint analysis. Among 52 genes around qPHR3A.2, four genes (TraesCS3A01G101900, TraesCS3A01G102200, TraesCS3A01G104100, and TraesCS3A01G105400) were found to be differentially expressed in low-nitrogen-tolerant wheat genotypes, while TaCLH2 (TraesCS3A01G101900) was putatively involved in porphyrin metabolism in KEGG enrichment analyses.

CONCLUSIONS

This study identified valuable candidate gene for low-N-tolerant wheat breeding and provides new insights into the genetic basis of low N tolerance in wheat.

Genetic regions determine tolerance to nitrogen deficiency in European elite bread wheats grown under contrasting nitrogen stress scenarios

KEY MESSAGE

Clustering 24 environments in four contrasting nitrogen stress scenarios enabled the detection of genetic regions determining tolerance to nitrogen deficiency in European elite bread wheats. Increasing the nitrogen use efficiency of wheat varieties is an important goal for breeding. However, most genetic studies of wheat grown at different nitrogen levels in the field report significant interactions with the genotype. The chromosomal regions possibly involved in these interactions are largely unknown. The objective of this study was to quantify the response of elite bread wheat cultivars to different nitrogen field stress scenarios and identify genomic regions involved in this response. For this purpose, 212 elite bread wheat varieties were grown in a multi-environment trial at different nitrogen levels. Genomic regions associated with grain yield, protein concentration and grain protein deviation responses to nitrogen deficiency were identified. Environments were clustered according to adjusted means for grain yield, yield components and grain protein concentration. Four nitrogen availability scenarios were identified: optimal condition, moderate early deficiency, severe late deficiency, and severe continuous deficiency. A large range of tolerance to nitrogen deficiency was observed among varieties, which were ranked differently in different nitrogen deficiency scenarios. The well-known negative correlation between grain yield and grain protein concentration also existed between their respective tolerance indices. Interestingly, the tolerance indices for grain yield and grain protein deviation were either null or weakly positive meaning that breeding for the two traits should be less difficult than expected. Twenty-two QTL regions were identified for the tolerance indices. By selecting associated markers, these regions may be selected separately or combined to improve the tolerance to N deficiency within a breeding programme.

Phenomics based prediction of plant biomass and leaf area in wheat using machine learning approaches

Introduction

Phenomics has emerged as important tool to bridge the genotype-phenotype gap. To dissect complex traits such as highly dynamic plant growth, and quantification of its component traits over a different growth phase of plant will immensely help dissect genetic basis of biomass production. Based on RGB images, models have been developed to predict biomass recently. However, it is very challenging to find a model performing stable across experiments. In this study, we recorded RGB and NIR images of wheat germplasm and Recombinant Inbred Lines (RILs) of Raj3765xHD2329, and examined the use of multimodal images from RGB, NIR sensors and machine learning models to predict biomass and leaf area non-invasively.

Results

The image-based traits (i-Traits) containing geometric features, RGB based indices, RGB colour classes and NIR features were categorized into architectural traits and physiological traits. Total 77 i-Traits were selected for prediction of biomass and leaf area consisting of 35 architectural and 42 physiological traits. We have shown that different biomass related traits such as fresh weight, dry weight and shoot area can be predicted accurately from RGB and NIR images using 16 machine learning models. We applied the models on two consecutive years of experiments and found that measurement accuracies were similar suggesting the generalized nature of models. Results showed that all biomass-related traits could be estimated with about 90% accuracy but the performance of model BLASSO was relatively stable and high in all the traits and experiments. The R² of BLASSO for fresh weight prediction was 0.96 (both year experiments), for dry weight prediction was 0.90 (Experiment 1) and 0.93 (Experiment 2) and for shoot area prediction 0.96 (Experiment 1) and 0.93 (Experiment 2). Also, the RMSRE of BLASSO for fresh weight prediction was 0.53 (Experiment 1) and 0.24 (Experiment 2), for dry weight prediction was 0.85 (Experiment 1) and 0.25 (Experiment 2) and for shoot area prediction 0.59 (Experiment 1) and 0.53 (Experiment 2).

Discussion

Based on the quantification power analysis of i-Traits, the determinants of biomass accumulation were found which contains both architectural and physiological traits. The best predictor i-Trait for fresh weight and dry weight prediction was Area_SV and for shoot area prediction was projected shoot area. These results will be helpful for identification and genetic basis dissection of major determinants of biomass accumulation and also non-invasive high throughput estimation of plant growth during different phenological stages can identify hitherto uncovered genes for biomass production and its deployment in crop improvement for breaking the yield plateau.

Identifying genomic regions determining shoot and root traits related to nitrogen uptake efficiency in a multiparent advanced generation intercross (MAGIC) winter wheat population in a high-throughput phenotyping facility

In the context of a continuously increasing human population that needs to be fed, with environmental protection in mind, nitrogen use efficiency (NUE) improvement is becoming very important. To understand the natural variation of traits linked to nitrogen uptake efficiency (UPE), one component of NUE, the multiparent advanced generation intercross (MAGIC) winter wheat population WM-800 was phenotyped under two contrasting nitrogen (N) levels in a high-throughput phenotyping facility for six weeks. Three biomass-related, three root-related, and two reflectance-related traits were measured weekly under each treatment. Subsequently, the population was genetically analysed using a total of 13,060 polymorphic haplotypes and singular SNPs for a genome-wide association study (GWAS). In total, we detected 543 quantitative trait loci (QTL) across all time points and traits, which were pooled into 42 stable QTL (sQTL; present in at least three of the six weeks). Besides Rht-B1 and Rht-D1, candidate genes playing a role in gibberellic acid-regulated growth and nitrate transporter genes from the NPF gene family, like NRT 1.1, were linked to sQTL. Two novel sQTL on chromosomes 5 A and 6D showed pleiotropic effects on several traits. The high number of N-specific sQTL indicates that selection for UPE is useful specifically under N-limited conditions.

Applications of hyperspectral imaging in plant phenotyping

Our ability to interrogate and manipulate the genome far exceeds our capacity to measure the effects of genetic changes on plant traits. Much effort has been made recently by the plant science research community to address this imbalance. The responses of plants to environmental conditions can now be defined using a variety of imaging approaches. Hyperspectral imaging (HSI) has emerged as a promising approach to measure traits using a wide range of wavebands simultaneously in 3D to capture information in lab, glasshouse, or field settings. HSI has been applied to define abiotic, biotic, and quality traits for optimisation of crop management.

Advanced high-throughput plant phenotyping techniques for genome-wide association studies: A review

Linking phenotypes and genotypes to identify genetic architectures that regulate important traits is crucial for plant breeding and the development of plant genomics. In recent years, genome-wide association studies (GWASs) have been applied extensively to interpret relationships between genes and traits. Successful GWAS application requires comprehensive genomic and phenotypic data from large populations. Although multiple high-throughput DNA sequencing approaches are available for the generation of genomics data, the capacity to generate high-quality phenotypic data is lagging far behind. Traditional methods for plant phenotyping mostly rely on manual measurements, which are laborious, inaccurate, and time-consuming, greatly impairing the acquisition of phenotypic data from large populations. In contrast, high-throughput phenotyping has unique advantages, facilitating rapid, non-destructive, and high-throughput detection, and, in turn, addressing the shortcomings of traditional methods. Aim of Review: This review summarizes the current status with regard to the integration of high-throughput phenotyping and GWAS in plants, in addition to discussing the inherent challenges and future prospects. Key Scientific Concepts of Review: High-throughput phenotyping, which facilitates non-contact and dynamic measurements, has the potential to offer high-quality trait data for GWAS and, in turn, to enhance the unraveling of genetic structures of complex plant traits. In conclusion, high-throughput phenotyping integration with GWAS could facilitate the revealing of coding information in plant genomes.

Rejuvenation remodels transcriptional network to improve rhizogenesis in mature Juglans tree

Adventitious rooting of walnut species (Juglans L.) is known to be rather difficult, especially for mature trees. The adventitious root formation (ARF) capacities of mature trees can be significantly improved by rejuvenation. However, the underlying gene regulatory networks (GRNs) of rejuvenation remain largely unknown. To characterize such regulatory networks, we carried out the transcriptomic study using RNA samples of the cambia and peripheral tissues on the bottom of rejuvenated and mature walnut (Juglans hindsii × J. regia) cuttings during the ARF. The RNA sequencing data suggested that zeatin biosynthesis, energy metabolism and substance metabolism were activated by rejuvenation, whereas photosynthesis, fatty acid biosynthesis and the synthesis pathways for secondary metabolites were inhibited. The inter- and intra-module GRNs were constructed using differentially expressed genes. We identified 35 hub genes involved in five modules associated with ARF. Among these hub genes, particularly, beta-glucosidase-like (BGLs) family members involved in auxin metabolism were overexpressed at the early stage of the ARF. Furthermore, BGL12 from the cuttings of Juglans was overexpressed in Populus alba × P. glandulosa. Accelerated ARF and increased number of ARs were observed in the transgenic poplars. These results provide a high-resolution atlas of gene activity during ARF and help to uncover the regulatory modules associated with the ARF promoted by rejuvenation.

Functional Mapping of Phenotypic Plasticity of Staphylococcus aureus Under Vancomycin Pressure

Phenotypic plasticity is the exhibition of various phenotypic traits produced by a single genotype in response to environmental changes, enabling organisms to adapt to environmental changes by maintaining growth and reproduction. Despite its significance in evolutionary studies, we still know little about the genetic control of phenotypic plasticity. In this study, we designed and conducted a genome-wide association study (GWAS) to reveal genetic architecture of how Staphylococcus aureus strains respond to increasing concentrations of vancomycin (0, 2, 4, and 6 μg/mL) in a time course. We implemented functional mapping, a dynamic model for genetic mapping using longitudinal data, to map specific loci that mediate the growth trajectories of abundance of vancomycin-exposed S. aureus strains. 78 significant single nucleotide polymorphisms were identified following analysis of the whole growth and development process, and seven genes might play a pivotal role in governing phenotypic plasticity to the pressure of vancomycin. These seven genes, SAOUHSC_00020 (walR), SAOUHSC_00176, SAOUHSC_00544 (sdrC), SAOUHSC_02998, SAOUHSC_00025, SAOUHSC_00169, and SAOUHSC_02023, were found to help S. aureus regulate antibiotic pressure. Our dynamic gene mapping technique provides a tool for dissecting the phenotypic plasticity mechanisms of S. aureus under vancomycin pressure, emphasizing the feasibility and potential of functional mapping in the study of bacterial phenotypic plasticity.

Functional physiological phenotyping with functional mapping: A general framework to bridge the phenotype-genotype gap in plant physiology

The recent years have witnessed the emergence of high-throughput phenotyping techniques. In particular, these techniques can characterize a comprehensive landscape of physiological traits of plants responding to dynamic changes in the environment. These innovations, along with the next-generation genomic technologies, have brought plant science into the big-data era. However, a general framework that links multifaceted physiological traits to DNA variants is still lacking. Here, we developed a general framework that integrates functional physiological phenotyping (FPP) with functional mapping (FM). This integration, implemented with high-dimensional statistical reasoning, can aid in our understanding of how genotype is translated toward phenotype. As a demonstration of method, we implemented the transpiration and soil-plant-atmosphere measurements of a tomato introgression line population into the FPP-FM framework, facilitating the identification of quantitative trait loci (QTLs) that mediate the spatiotemporal change of transpiration rate and the test of how these QTLs control, through their interaction networks, phenotypic plasticity under drought stress.

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.

Management and Characterization of Abiotic Stress via PhénoField®, a High-Throughput Field Phenotyping Platform

In order to evaluate the impact of water deficit in field conditions, researchers or breeders must set up large experiment networks in very restrictive field environments. Experience shows that half of the field trials are not relevant because of climatic conditions that do not allow the stress scenario to be tested. The PhénoField® platform is the first field based infrastructure in the European Union to ensure protection against rainfall for a large number of plots, coupled with the non-invasive acquisition of crops' phenotype. In this paper, we will highlight the PhénoField® production capability using data from 2017-wheat trial. The innovative approach of the PhénoField® platform consists in the use of automatic irrigating rainout shelters coupled with high throughput field phenotyping to complete conventional phenotyping and micrometeorological densified measurements. Firstly, to test various abiotic stresses, automatic mobile rainout shelters allow fine management of fertilization or irrigation by driving daily the intensity and period of the application of the desired limiting factor on the evaluated crop. This management is based on micro-meteorological measurements coupled with a simulation of a carbon, water and nitrogen crop budget. Furthermore, as high-throughput plant-phenotyping under controlled conditions is well advanced, comparable evaluation in field conditions is enabled through phenotyping gantries equipped with various optical sensors. This approach, giving access to either similar or innovative variables compared manual measurements, is moreover distinguished by its capacity for dynamic analysis. Thus, the interactions between genotypes and the environment can be deciphered and better detailed since this gives access not only to the environmental data but also to plant responses to limiting hydric and nitrogen conditions. Further data analyses provide access to the curve parameters of various indicator kinetics, all the more integrative and relevant of plant behavior under stressful conditions. All these specificities of the PhénoField® platform open the way to the improvement of various categories of crop models, the fine characterization of variety behavior throughout the growth cycle and the evaluation of particular sensors better suited to a specific research question.