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

Temporally resolved growth patterns reveal novel information about the polygenic nature of complex quantitative traits

Plant height can be an indicator of plant health across environments and used to identify superior genotypes. Typically plant height is measured at a single timepoint when plants reach terminal height. Evaluating plant height using unoccupied aerial vehicles allows for measurements throughout the growing season, facilitating a better understanding of plant-environment interactions and the genetic basis of this complex trait. To assess variation throughout development, plant height data was collected from planting until terminal height at anthesis (14 flights 2018, 27 in 2019, 12 in 2020, and 11 in 2021) for a panel of ~500 diverse maize inbred lines. The percent variance explained in plant height throughout the season was significantly explained by genotype (9-48%), year (4-52%), and genotype-by-year interactions (14-36%) to varying extents throughout development. Genome-wide association studies revealed 717 significant single nucleotide polymorphisms associated with plant height and growth rate at different parts of the growing season specific to certain phases of vegetative growth. When plant height growth curves were compared to growth curves estimated from canopy cover, greater Fréchet distance stability was observed in plant height growth curves than for canopy cover. This indicated canopy cover may be more useful for understanding environmental modulation of overall plant growth and plant height better for understanding genotypic modulation of overall plant growth. This study demonstrated that substantial information can be gained from high temporal resolution data to understand how plants differentially interact with the environment and can enhance our understanding of the genetic basis of complex polygenic traits.

TWAS facilitates gene-scale trait genetic dissection through gene expression, structural variations, and alternative splicing in soybean

A genome-wide association study (GWAS) identifies trait-associated loci, but identifying the causal genes can be a bottleneck, due in part to slow decay of linkage disequilibrium (LD). A transcriptome-wide association study (TWAS) addresses this issue by identifying gene expression-phenotype associations or integrating gene expression quantitative trait loci with GWAS results. Here, we used self-pollinated soybean (Glycine max [L.] Merr.) as a model to evaluate the application of TWAS to the genetic dissection of traits in plant species with slow LD decay. We generated RNA sequencing data for a soybean diversity panel and identified the genetic expression regulation of 29 286 soybean genes. Different TWAS solutions were less affected by LD and were robust to the source of expression, identifing known genes related to traits from different tissues and developmental stages. The novel pod-color gene L2 was identified via TWAS and functionally validated by genome editing. By introducing a new exon proportion feature, we significantly improved the detection of expression variations that resulted from structural variations and alternative splicing. As a result, the genes identified through our TWAS approach exhibited a diverse range of causal variations, including SNPs, insertions or deletions, gene fusion, copy number variations, and alternative splicing. Using this approach, we identified genes associated with flowering time, including both previously known genes and novel genes that had not previously been linked to this trait, providing insights complementary to those from GWAS. In summary, this study supports the application of TWAS for candidate gene identification in species with low rates of LD decay.

Cotton morphological traits tracking through spatiotemporal registration of terrestrial laser scanning time-series data

Understanding the complex interactions between genotype-environment dynamics is fundamental for optimizing crop improvement. However, traditional phenotyping methods limit assessments to the end of the growing season, restricting continuous crop monitoring. To address this limitation, we developed a methodology for spatiotemporal registration of time-series 3D point cloud data, enabling field phenotyping over time for accurate crop growth tracking. Leveraging multi-scan terrestrial laser scanning (TLS), we captured high-resolution 3D LiDAR data in a cotton breeding field across various stages of the growing season to generate four-dimensional (4D) crop models, seamlessly integrating spatial and temporal dimensions. Our registration procedure involved an initial pairwise terrain-based matching for rough alignment, followed by a bird's-eye view adjustment for fine registration. Point clouds collected throughout nine sessions across the growing season were successfully registered both spatially and temporally, with average registration errors of approximately 3 cm. We used the generated 4D models to monitor canopy height (CH) and volume (CV) for eleven cotton genotypes over two months. The consistent height reference established via our spatiotemporal registration process enabled precise estimations of CH (R 2 = 0.95, RMSE = 7.6 cm). Additionally, we analyzed the relationship between CV and the interception of photosynthetically active radiation (IPAR f ), finding that it followed a curve with exponential saturation, consistent with theoretical models, with a standard error of regression (SER) of 11%. In addition, we compared mathematical models from the Richards family of sigmoid curves for crop growth modeling, finding that the logistic model effectively captured CH and CV evolution, aiding in identifying significant genotype differences. Our novel TLS-based digital phenotyping methodology enhances precision and efficiency in field phenotyping over time, advancing plant phenomics and empowering efficient decision-making for crop improvement efforts.

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.

Identification of potential Auxin Response Candidate genes for soybean rapid canopy coverage through comparative evolution and expression analysis

Glycine max, soybean, is an abundantly cultivated crop worldwide. Efforts have been made over the past decades to improve soybean production in traditional and organic agriculture, driven by growing demand for soybean-based products. Rapid canopy cover development (RCC) increases soybean yields and suppresses early-season weeds. Genome-wide association studies have found natural variants associated with RCC, however causal mechanisms are unclear. Auxin modulates plant growth and development and has been implicated in RCC traits. Therefore, modulation of auxin regulatory genes may enhance RCC. Here, we focus on the use of genomic tools and existing datasets to identify auxin signaling pathway RCC candidate genes, using a comparative phylogenetics and expression analysis approach. We identified genes encoding 14 TIR1/AFB auxin receptors, 61 Aux/IAA auxin co-receptors and transcriptional co-repressors, and 55 ARF auxin response factors in the soybean genome. We used Bayesian phylogenetic inference to identify soybean orthologs of Arabidopsis thaliana genes, and defined an ortholog naming system for these genes. To further define potential auxin signaling candidate genes for RCC, we examined tissue-level expression of these genes in existing datasets and identified highly expressed auxin signaling genes in apical tissues early in development. We identified at least 4 TIR1/AFB, 8 Aux/IAA, and 8 ARF genes with highly specific expression in one or more RCC-associated tissues. We hypothesize that modulating the function of these genes through gene editing or traditional breeding will have the highest likelihood of affecting RCC while minimizing pleiotropic effects.

Next-generation bulked segregant analysis for Breeding 4.0

Functional cloning and manipulation of genes controlling various agronomic traits are important for boosting crop production. Although bulked segregant analysis (BSA) is an efficient method for functional cloning, its low throughput cannot satisfy the current need for crop breeding and food security. Here, we review the rationale and development of conventional BSA and discuss its strengths and drawbacks. We then propose next-generation BSA (NG-BSA) integrating multiple cutting-edge technologies, including high-throughput phenotyping, biological big data, and the use of machine learning. NG-BSA increases the resolution of genetic mapping and throughput for cloning quantitative trait genes (QTGs) and optimizes candidate gene selection while providing a means to elucidate the interaction network of QTGs. The ability of NG-BSA to efficiently batch-clone QTGs makes it an important tool for dissecting molecular mechanisms underlying various traits, as well as for the improvement of Breeding 4.0 strategy, especially in targeted improvement and population improvement of crops.

The Importance of Using Realistic 3D Canopy Models to Calculate Light Interception in the Field

Quantifying canopy light interception provides insight into the effects of plant spacing, canopy structure, and leaf orientation on radiation distribution. This is essential for increasing crop yield and improving product quality. Canopy light interception can be quantified using 3-dimensional (3D) plant models and optical simulations. However, virtual 3D canopy models (VCMs) have often been used to quantify canopy light interception because realistic 3D canopy models (RCMs) are difficult to obtain in the field. This study aims to compare the differences in light interception between VCMs and RCM. A realistic 3D maize canopy model (RCM) was reconstructed over a large area of the field using an advanced unmanned aerial vehicle cross-circling oblique (CCO) route and the structure from motion-multi-view stereo method. Three types of VCMs (VCM-1, VCM-4, and VCM-8) were then created by replicating 1, 4, and 8 individual realistic plants constructed by CCO in the center of the corresponding RCM. The daily light interception per unit area (DLI), as computed for the 3 VCMs, exhibited marked deviation from the RCM, as evinced by the relative root mean square error (rRMSE) values of 20.22%, 17.38%, and 15.48%, respectively. Although this difference decreased as the number of plants used to replicate the virtual canopy increased, rRMSE of DLI for VCM-8 and RCM still reached 15.48%. It was also found that the difference in light interception between RCMs and VCMs was substantially smaller in the early stage (48 days after sowing [DAS]) than in the late stage (70 DAS). This study highlights the importance of using RCM when calculating light interception in the field, especially in the later growth stages of plants.

UAV-based time-series phenotyping reveals the genetic basis of plant height in upland cotton

Plant height (PH) is an important agronomic trait affecting crop architecture, biomass, resistance to lodging and mechanical harvesting. Elucidating the genetic governance of plant height is crucial because of the global demand for high crop yields. However, during the rapid growth period of plants the PH changes a lot on a daily basis, which makes it difficult to accurately phenotype the trait by hand on a large scale. In this study, an unmanned aerial vehicle (UAV)-based remote-sensing phenotyping platform was applied to obtain time-series PHs of 320 upland cotton accessions in three different field trials. The results showed that the PHs obtained from UAV images were significantly correlated with ground-based manual measurements, for three trials (R2  = 0.96, 0.95 and 0.96). Two genetic loci on chromosomes A01 and A11 associated with PH were identified by genome-wide association studies (GWAS). GhUBP15 and GhCUL1 were identified to influence PH in further analysis. We obtained a time series of PH values for three field conditions based on remote sensing with UAV. The key genes identified in this study are of great value for the breeding of ideal plant architecture in cotton.

The elite variations in germplasms for soybean breeding

The genetic base of soybean cultivars (Glycine max (L.) Merr.) has been narrowed through selective domestication and specific breeding improvement, similar to other crops. This presents challenges in breeding new cultivars with improved yield and quality, reduced adaptability to climate change, and increased susceptibility to diseases. On the other hand, the vast collection of soybean germplasms offers a potential source of genetic variations to address those challenges, but it has yet to be fully leveraged. In recent decades, rapidly improved high-throughput genotyping technologies have accelerated the harness of elite variations in soybean germplasm and provided the important information for solving the problem of a narrowed genetic base in breeding. In this review, we will overview the situation of maintenance and utilization of soybean germplasms, various solutions provided for different needs in terms of the number of molecular markers, and the omics-based high-throughput strategies that have been used or can be used to identify elite alleles. We will also provide an overall genetic information generated from soybean germplasms in yield, quality traits, and pest resistance for molecular breeding.

Auxin regulation on crop: from mechanisms to opportunities in soybean breeding

Breeding crop varieties with high yield and ideal plant architecture is a desirable goal of agricultural science. The success of "Green Revolution" in cereal crops provides opportunities to incorporate phytohormones in crop breeding. Auxin is a critical phytohormone to determine nearly all the aspects of plant development. Despite the current knowledge regarding auxin biosynthesis, auxin transport and auxin signaling have been well characterized in model Arabidopsis (Arabidopsis thaliana) plants, how auxin regulates crop architecture is far from being understood, and the introduction of auxin biology in crop breeding stays in the theoretical stage. Here, we give an overview on molecular mechanisms of auxin biology in Arabidopsis, and mainly summarize auxin contributions for crop plant development. Furthermore, we propose potential opportunities to integrate auxin biology in soybean (Glycine max) breeding.

Estimation of soybean yield parameters under lodging conditions using RGB information from unmanned aerial vehicles

The estimation of yield parameters based on early data is helpful for agricultural policymakers and food security. Developments in unmanned aerial vehicle (UAV) platforms and sensor technology help to estimate yields efficiency. Previous studies have been based on less cultivars (<10) and ideal experimental environments, it is not available in practical production. Therefore, the objective of this study was to estimate the yield parameters of soybean (Glycine max (L.) Merr.) under lodging conditions using RGB information. In this study, 17 time point data throughout the soybean growing season in Nanchang, Jiangxi Province, China, were collected, and the vegetation index, texture information, canopy cover, and crop height were obtained by UAV-image processing. After that, partial least squares regression (PLSR), logistic regression (Logistic), random forest regression (RFR), support vector machine regression (SVM), and deep learning neural network (DNN) were used to estimate the yield parameters. The results can be summarized as follows: (1) The most suitable time point to estimate the yield was flowering stage (48 days), which was when most of the soybean cultivars flowered. (2) The multiple data fusion improved the accuracy of estimating the yield parameters, and the texture information has a high potential to contribute to the estimation of yields, and (3) The DNN model showed the best accuracy of training (R²=0.66 rRMSE=32.62%) and validation (R²=0.50, rRMSE=43.71%) datasets. In conclusion, these results provide insights into both best estimate period selection and early yield estimation under lodging condition when using remote sensing.