Unoccupied aerial system enabled functional modeling of maize height reveals dynamic expression of loci

Field-based high-throughput phenotyping enhances phenomic and genomic predictions for grain yield and plant height across years in maize

Field-based phenomic prediction employs novel features, like vegetation indices (VIs) from drone images, to predict key agronomic traits in maize, despite challenges in matching biomarker measurement time points across years or environments. This study utilized functional principal component analysis (FPCA) to summarize the variation of temporal VIs, uniquely allowing the integration of this data into phenomic prediction models tested across multiple years (2018-2021) and environments. The models, which included 1 genomic, 2 phenomic, 2 multikernel, and 1 multitrait type, were evaluated in 4 prediction scenarios (CV2, CV1, CV0, and CV00), relevant for plant breeding programs, assessing both tested and untested genotypes in observed and unobserved environments. Two hybrid populations (415 and 220 hybrids) demonstrated the visible atmospherically resistant index's strong temporal correlation with grain yield (up to 0.59) and plant height. The first 2 FPCAs explained 59.3 ± 13.9% and 74.2 ± 9.0% of the temporal variation of temporal data of VIs, respectively, facilitating predictions where flight times varied. Phenomic data, particularly when combined with genomic data, often were comparable to or numerically exceeded the base genomic model in prediction accuracy, particularly for grain yield in untested hybrids, although no significant differences in these models' performance were consistently observed. Overall, this approach underscores the effectiveness of FPCA and combined models in enhancing the prediction of grain yield and plant height across environments and diverse agricultural settings.

Deciphering temporal growth patterns in maize: integrative modeling of phenotype dynamics and underlying genomic variations

Quantifying the temporal or longitudinal growth dynamics of crops in diverse environmental conditions is crucial for understanding plant development, requiring further modeling techniques. In this study, we analyzed the growth patterns of two different maize (Zea mays L.) populations using high-throughput phenotyping with a maize population consisting of 515 recombinant inbred lines (RILs) grown in Texas and a hybrid population containing 1090 hybrids grown in Missouri. Two models, Gaussian peak and functional principal component analysis (FPCA), were employed to study the Normalized Green-Red Difference Index (NGRDI) scores. The Gaussian peak model showed strong correlations (c. 0.94 for RILs and c. 0.97 for hybrids) between modeled and non-modeled temporal trajectories. Functional principal component analysis differentiated NGRDI trajectories in RILs under different conditions, capturing substantial variability (75%, 20%, and 5% for RILs; 88% and 12% for hybrids). By comparing these models with conventional BLUP values, common quantitative trait loci (QTLs) were identified, containing candidate genes of brd1, pin11, zcn8 and rap2. The harmony between these loci's additive effects and growing degree days, as well as the differentiation of RIL haplotypes across growth stages, underscores the significant interplay of these loci in driving plant development. These findings contribute to advancing understanding of plant-environment interactions and have implications for crop improvement strategies.

Phenomic data-driven biological prediction of maize through field-based high-throughput phenotyping integration with genomic data

High-throughput phenotyping (HTP) has expanded the dimensionality of data in plant research; however, HTP has resulted in few novel biological discoveries to date. Field-based HTP (FHTP), using small unoccupied aerial vehicles (UAVs) equipped with imaging sensors, can be deployed routinely to monitor segregating plant population interactions with the environment under biologically meaningful conditions. Here, flowering dates and plant height, important phenological fitness traits, were collected on 520 segregating maize recombinant inbred lines (RILs) in both irrigated and drought stress trials in 2018. Using UAV phenomic, single nucleotide polymorphism (SNP) genomic, as well as combined data, flowering times were predicted using several scenarios. Untested genotypes were predicted with 0.58, 0.59, and 0.41 prediction ability for anthesis, silking, and terminal plant height, respectively, using genomic data, but prediction ability increased to 0.77, 0.76, and 0.58 when phenomic and genomic data were used together. Using the phenomic data in a genome-wide association study, a heat-related candidate gene (GRMZM2G083810; hsp18f) was discovered using temporal reflectance phenotypes belonging to flowering times (both irrigated and drought) trials where heat stress also peaked. Thus, a relationship between plants and abiotic stresses belonging to a specific time of growth was revealed only through use of temporal phenomic data. Overall, this study showed that (i) it is possible to predict complex traits using high dimensional phenomic data between different environments, and (ii) temporal phenomic data can reveal a time-dependent association between genotypes and abiotic stresses, which can help understand mechanisms to develop resilient plants.

Time series canopy phenotyping enables the identification of genetic variants controlling dynamic phenotypes in soybean

Advances in plant phenotyping technologies are dramatically reducing the marginal costs of collecting multiple phenotypic measurements across several time points. Yet, most current approaches and best statistical practices implemented to link genetic and phenotypic variation in plants have been developed in an era of single-time-point data. Here, we used time-series phenotypic data collected with an unmanned aircraft system for a large panel of soybean (Glycine max (L.) Merr.) varieties to identify previously uncharacterized loci. Specifically, we focused on the dissection of canopy coverage (CC) variation from this rich data set. We also inferred the speed of canopy closure, an additional dimension of CC, from the time-series data, as it may represent an important trait for weed control. Genome-wide association studies (GWASs) identified 35 loci exhibiting dynamic associations with CC across developmental stages. The time-series data enabled the identification of 10 known flowering time and plant height quantitative trait loci (QTLs) detected in previous studies of adult plants and the identification of novel QTLs influencing CC. These novel QTLs were disproportionately likely to act earlier in development, which may explain why they were missed in previous single-time-point studies. Moreover, this time-series data set contributed to the high accuracy of the GWASs, which we evaluated by permutation tests, as evidenced by the repeated identification of loci across multiple time points. Two novel loci showed evidence of adaptive selection during domestication, with different genotypes/haplotypes favored in different geographic regions. In summary, the time-series data, with soybean CC as an example, improved the accuracy and statistical power to dissect the genetic basis of traits and offered a promising opportunity for crop breeding with quantitative growth curves.

Cumulative temporal vegetation indices from unoccupied aerial systems allow maize (Zea mays L.) hybrid yield to be estimated across environments with fewer flights

Unoccupied aerial systems (UAS) based high throughput phenotyping studies require further investigation to combine different environments and planting times into one model. Here 100 elite breeding hybrids of maize (Zea mays L.) were evaluated in two environment trials-one with optimal planting and irrigation (IHOT), and one dryland with delayed planting (DHOT). RGB (Red-Green-Blue) based canopy height measurement (CHM) and vegetation indices (VIs) were estimated from a UAS platform. Time series and cumulative VIs, by both summation (ΣVI-SUMs) and area under the curve (ΣVI-AUCs), were fit via machine learning regression modeling (random forest, linear, ridge, lasso, elastic net regressions) to estimate grain yield. VIs were more valuable predictors of yield to combine different environments than CHM. Time series VIs and CHM produced high accuracies (~68-72%), but inconsistent models. A little sacrifice in accuracy (~60-65%) produced consistent models using ΣVI-SUMs and CHM during pre-reproductive vegetative growth. Absence of VIs produced poorer accuracies (by about ~5-10%). Normalized difference type VIs produced maximum accuracies, and flowering times were the best times for UAS data acquisition. This study suggests that the best yielding varieties can be accurately predicted in new environments at or before flowering when combining multiple temporal flights and predictors.

Genome-wide identification of genes enabling accurate prediction of hybrid performance from parents across environments and populations for gene-based breeding in maize

Accurate prediction of hybrid offspring complex trait phenotype from parents is paramount to enhanced plant breeding, animal breeding, and human medicine. Here we report genome-wide identification of genes enabling accurate prediction of hybrid offspring complex traits from parents using maize grain yield as the target trait. We identified 181 ZmF1GY genes enabling prediction of maize (Zea mays L.) F₁ hybrid grain yield from parents and tested their utility and efficiency for predicting F₁ hybrid grain yields from parents using their expressions, genic SNPs, and number of favorable alleles (NFAs), respectively. The ZmF1GY genes predicted hybrid grain yields from parents at an accuracy of 0.86, presented by correlation coefficient between predicted and observed phenotypes, within an environment, 0.74 across environments, and 0.64 across populations, outperforming genomic prediction by 27-406%, 23%, and 40%, respectively. Furthermore, we identified nine of the ZmF1GY genes containing SNPs or InDels in parents that increased or decreased hybrid grain yields by 14-46%. When the NFAs of these nine ZmF1GY genes were used for hybrid grain yield prediction from parents, they predicted hybrid grain yields at an accuracy of 0.79, outperforming genomic prediction by 21% that was based on up to tens of thousands of genome-wide SNPs. These results demonstrate the feasibility of developing a gene toolkit for a species enabling gene-based breeding across environments and populations that is much more powerful and efficient than current breeding, thereby helping secure the world's food production. The methodology is applicable to all crops, livestock, and humans.

AirMeasurer: open-source software to quantify static and dynamic traits derived from multiseason aerial phenotyping to empower genetic mapping studies in rice

Low-altitude aerial imaging, an approach that can collect large-scale plant imagery, has grown in popularity recently. Amongst many phenotyping approaches, unmanned aerial vehicles (UAVs) possess unique advantages as a consequence of their mobility, flexibility and affordability. Nevertheless, how to extract biologically relevant information effectively has remained challenging. Here, we present AirMeasurer, an open-source and expandable platform that combines automated image analysis, machine learning and original algorithms to perform trait analysis using 2D/3D aerial imagery acquired by low-cost UAVs in rice (Oryza sativa) trials. We applied the platform to study hundreds of rice landraces and recombinant inbred lines at two sites, from 2019 to 2021. A range of static and dynamic traits were quantified, including crop height, canopy coverage, vegetative indices and their growth rates. After verifying the reliability of AirMeasurer-derived traits, we identified genetic variants associated with selected growth-related traits using genome-wide association study and quantitative trait loci mapping. We found that the AirMeasurer-derived traits had led to reliable loci, some matched with published work, and others helped us to explore new candidate genes. Hence, we believe that our work demonstrates valuable advances in aerial phenotyping and automated 2D/3D trait analysis, providing high-quality phenotypic information to empower genetic mapping for crop improvement.

Plant Genotype to Phenotype Prediction Using Machine Learning

Genomic prediction tools support crop breeding based on statistical methods, such as the genomic best linear unbiased prediction (GBLUP). However, these tools are not designed to capture non-linear relationships within multi-dimensional datasets, or deal with high dimension datasets such as imagery collected by unmanned aerial vehicles. Machine learning (ML) algorithms have the potential to surpass the prediction accuracy of current tools used for genotype to phenotype prediction, due to their capacity to autonomously extract data features and represent their relationships at multiple levels of abstraction. This review addresses the challenges of applying statistical and machine learning methods for predicting phenotypic traits based on genetic markers, environment data, and imagery for crop breeding. We present the advantages and disadvantages of explainable model structures, discuss the potential of machine learning models for genotype to phenotype prediction in crop breeding, and the challenges, including the scarcity of high-quality datasets, inconsistent metadata annotation and the requirements of ML models.

Phenomic data-facilitated rust and senescence prediction in maize using machine learning algorithms

Current methods in measuring maize (Zea mays L.) southern rust (Puccinia polyspora Underw.) and subsequent crop senescence require expert observation and are resource-intensive and prone to subjectivity. In this study, unoccupied aerial system (UAS) field-based high-throughput phenotyping (HTP) was employed to collect high-resolution aerial imagery of elite maize hybrids planted in the 2020 and 2021 growing seasons, with 13 UAS flights obtained from 2020 and 17 from 2021. In total, 36 vegetation indices (VIs) were extracted from mosaicked aerial images that served as temporal phenomic predictors for southern rust scored in the field and senescence as scored using UAS-acquired mosaic images. Temporal best linear unbiased predictors (TBLUPs) were calculated using a nested model that treated hybrid performance as nested within flights in terms of rust and senescence. All eight machine learning regressions tested (ridge, lasso, elastic net, random forest, support vector machine with radial and linear kernels, partial least squares, and k-nearest neighbors) outperformed a general linear model with both higher prediction accuracies (92-98%) and lower root mean squared error (RMSE) for rust and senescence scores (linear model RMSE ranged from 65.8 to 2396.5 across all traits, machine learning regressions RMSE ranged from 0.3 to 17.0). UAS-acquired VIs enabled the discovery of novel early quantitative phenotypic indicators of maize senescence and southern rust before being detectable by expert annotation and revealed positive correlations between grain filling time and yield (0.22 and 0.44 in 2020 and 2021), with practical implications for precision agricultural practices.

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.

Unoccupied aerial systems discovered overlooked loci capturing the variation of entire growing period in maize

Traditional phenotyping methods, coupled with genetic mapping in segregating populations, have identified loci governing complex traits in many crops. Unoccupied aerial systems (UAS)-based phenotyping has helped to reveal a more novel and dynamic relationship between time-specific associated loci with complex traits previously unable to be evaluated. Over 1,500 maize (Zea mays L.) hybrid row plots containing 280 different replicated maize hybrids from the Genomes to Fields (G2F) project were evaluated agronomically and using UAS in 2017. Weekly UAS flights captured variation in plant heights during the growing season under three different management conditions each year: optimal planting with irrigation (G2FI), optimal dryland planting without irrigation (G2FD), and a stressed late planting (G2LA). Plant height of different flights were ranked based on importance for yield using a random forest (RF) algorithm. Plant heights captured by early flights in G2FI trials had higher importance (based on Gini scores) for predicting maize grain yield (GY) but also higher accuracies in genomic predictions which fluctuated for G2FD (-0.06∼0.73), G2FI (0.33∼0.76), and G2LA (0.26∼0.78) trials. A genome-wide association analysis discovered 52 significant single nucleotide polymorphisms (SNPs), seven were found consistently in more than one flights or trial; 45 were flight or trial specific. Total cumulative marker effects for each chromosome's contributions to plant height also changed depending on flight. Using UAS phenotyping, this study showed that many candidate genes putatively play a role in the regulation of plant architecture even in relatively early stages of maize growth and development.

Temporal Vegetation Indices and Plant Height from Remotely Sensed Imagery Can Predict Grain Yield and Flowering Time Breeding Value in Maize via Machine Learning Regression

Unoccupied aerial system (UAS; i.e., drone equipped with sensors) field-based high-throughput phenotyping (HTP) platforms are used to collect high quality images of plant nurseries to screen genetic materials (e.g., hybrids and inbreds) throughout plant growth at relatively low cost. In this study, a set of 100 advanced breeding maize (Zea mays L.) hybrids were planted at optimal (OHOT trial) and delayed planting dates (DHOT trial). Twelve UAS surveys were conducted over the trials throughout the growing season. Fifteen vegetative indices (VIs) and the 99th percentile canopy height measurement (CHMs) were extracted from processed UAS imagery (orthomosaics and point clouds) which were used to predict plot-level grain yield, days to anthesis (DTA), and silking (DTS). A novel statistical approach utilizing a nested design was fit to predict temporal best linear unbiased predictors (TBLUP) for the combined temporal UAS data. Our results demonstrated machine learning-based regressions (ridge, lasso, and elastic net) had from 4- to 9-fold increases in the prediction accuracies and from 13- to 73-fold reductions in root mean squared error (RMSE) compared to classical linear regression in prediction of grain yield or flowering time. Ridge regression performed best in predicting grain yield (prediction accuracy = ~0.6), while lasso and elastic net regressions performed best in predicting DTA and DTS (prediction accuracy = ~0.8) consistently in both trials. We demonstrated that predictor variable importance descended towards the terminal stages of growth, signifying the importance of phenotype collection beyond classical terminal growth stages. This study is among the first to demonstrate an ability to predict yield in elite hybrid maize breeding trials using temporal UAS image-based phenotypes and supports the potential benefit of phenomic selection approaches in estimating breeding values before harvest.

Validation of Functional Polymorphisms Affecting Maize Plant Height by Unoccupied Aerial Systems (UAS) Discovers Novel Temporal Phenotypes

Plant height (PHT) in maize (Zea mays L.) has been scrutinized genetically and phenotypically due to relationship with other agronomically valuable traits (e.g. yield). Heritable variation of PHT is determined by many discovered quantitative trait loci (QTLs); however, phenotypic effects of such loci often lack validation across environments and genetic backgrounds, especially in the hybrid state grown by farmers rather than the inbred state more often used by geneticists. A previous genome wide association study using a topcrossed hybrid diversity panel identified two novel quantitative trait variants (QTVs) controlling both PHT and grain yield. Here, heterogeneous inbred families demonstrated that these two loci, characterized by two single nucleotide polymorphisms (SNPs), cause phenotypic variation in inbred lines, but that size of these effects were variable across four different genetic backgrounds, ranging from 1 to 10 cm. Weekly unoccupied aerial system flights demonstrated the two SNPs had larger effects, varying from 10 to 25 cm, in early growth while effects decreased towards the end of the season. These results show that allelic effect sizes of economically valuable loci are both dynamic in temporal growth and dynamic across genetic backgrounds, resulting in informative phenotypic variability overlooked following traditional phenotyping methods. Public genotyping data shows recent favorable allele selection in elite temperate germplasm with little change across tropical backgrounds. As these loci remain rarer in tropical germplasm, with effects most visible early in growth, they are useful for breeding and selection to expand the genetic basis of maize.

High-Throughput Plant Phenotyping Platform (HT3P) as a Novel Tool for Estimating Agronomic Traits From the Lab to the Field

Food scarcity, population growth, and global climate change have propelled crop yield growth driven by high-throughput phenotyping into the era of big data. However, access to large-scale phenotypic data has now become a critical barrier that phenomics urgently must overcome. Fortunately, the high-throughput plant phenotyping platform (HT3P), employing advanced sensors and data collection systems, can take full advantage of non-destructive and high-throughput methods to monitor, quantify, and evaluate specific phenotypes for large-scale agricultural experiments, and it can effectively perform phenotypic tasks that traditional phenotyping could not do. In this way, HT3Ps are novel and powerful tools, for which various commercial, customized, and even self-developed ones have been recently introduced in rising numbers. Here, we review these HT3Ps in nearly 7 years from greenhouses and growth chambers to the field, and from ground-based proximal phenotyping to aerial large-scale remote sensing. Platform configurations, novelties, operating modes, current developments, as well the strengths and weaknesses of diverse types of HT3Ps are thoroughly and clearly described. Then, miscellaneous combinations of HT3Ps for comparative validation and comprehensive analysis are systematically present, for the first time. Finally, we consider current phenotypic challenges and provide fresh perspectives on future development trends of HT3Ps. This review aims to provide ideas, thoughts, and insights for the optimal selection, exploitation, and utilization of HT3Ps, and thereby pave the way to break through current phenotyping bottlenecks in botany.

The potential of remote sensing and artificial intelligence as tools to improve the resilience of agriculture production systems

Modern agriculture and food production systems are facing increasing pressures from climate change, land and water availability, and, more recently, a pandemic. These factors are threatening the environmental and economic sustainability of current and future food supply systems. Scientific and technological innovations are needed more than ever to secure enough food for a fast-growing global population. Scientific advances have led to a better understanding of how various components of the agricultural system interact, from the cell to the field level. Despite incredible advances in genetic tools over the past few decades, our ability to accurately assess crop status in the field, at scale, has been severely lacking until recently. Thanks to recent advances in remote sensing and Artificial Intelligence (AI), we can now quantify field scale phenotypic information accurately and integrate the big data into predictive and prescriptive management tools. This review focuses on the use of recent technological advances in remote sensing and AI to improve the resilience of agricultural systems, and we will present a unique opportunity for the development of prescriptive tools needed to address the next decade's agricultural and human nutrition challenges.

Canopy Roughness: A New Phenotypic Trait to Estimate Aboveground Biomass from Unmanned Aerial System

Cost-effective phenotyping methods are urgently needed to advance crop genetics in order to meet the food, fuel, and fiber demands of the coming decades. Concretely, characterizing plot level traits in fields is of particular interest. Recent developments in high-resolution imaging sensors for UAS (unmanned aerial systems) focused on collecting detailed phenotypic measurements are a potential solution. We introduce canopy roughness as a new plant plot-level trait. We tested its usability with soybean by optical data collected from UAS to estimate biomass. We validate canopy roughness on a panel of 108 soybean [Glycine max (L.) Merr.] recombinant inbred lines in a multienvironment trial during the R2 growth stage. A senseFly eBee UAS platform obtained aerial images with a senseFly S.O.D.A. compact digital camera. Using a structure from motion (SfM) technique, we reconstructed 3D point clouds of the soybean experiment. A novel pipeline for feature extraction was developed to compute canopy roughness from point clouds. We used regression analysis to correlate canopy roughness with field-measured aboveground biomass (AGB) with a leave-one-out cross-validation. Overall, our models achieved a coefficient of determination (R ²) greater than 0.5 in all trials. Moreover, we found that canopy roughness has the ability to discern AGB variations among different genotypes. Our test trials demonstrate the potential of canopy roughness as a reliable trait for high-throughput phenotyping to estimate AGB. As such, canopy roughness provides practical information to breeders in order to select phenotypes on the basis of UAS data.