ZmPTOX1, a plastid terminal oxidase, contributes to redox homeostasis during seed development and germination.

Maize plastid terminal oxidase1 (ZmPTOX1) plays a pivotal role in seed development by upholding redox balance within seed plastids. This study focuses on characterizing the white kernel mutant 3735 (wk3735) mutant, which yields pale-yellow seeds characterized by heightened protein but reduced carotenoid levels, along with delayed germination compared to wild-type (WT) seeds. We successfully cloned and identified the target gene ZmPTOX1, responsible for encoding maize PTOX-a versatile plastoquinol oxidase and redox sensor located in plastid membranes. While PTOX's established role involves regulating redox states and participating in carotenoid metabolism in Arabidopsis leaves and tomato fruits, our investigation marks the first exploration of its function in storage organs lacking a photosynthetic system. Through our research, we validated the existence of plastid-localized ZmPTOX1, existing as a homomultimer, and established its interaction with ferredoxin-NADP+ oxidoreductase 1 (ZmFNR1), a crucial component of the electron transport chain (ETC). This interaction contributes to the maintenance of redox equilibrium within plastids. Our findings indicate a propensity for excessive accumulation of reactive oxygen species (ROS) in wk3735 seeds. Beyond its known role in carotenoids' antioxidant properties, ZmPTOX1 also impacts ROS homeostasis owing to its oxidizing function. Altogether, our results underscore the critical involvement of ZmPTOX1 in governing seed development and germination by preserving redox balance within the seed plastids.

Genetic regulation of self-organizing azimuthal canopy orientations and their impacts on light interception in maize.

The efficiency of solar radiation interception contributes to the photosynthetic efficiency of crop plants. Light interception is a function of canopy architecture, including plant density; leaf number, length, width, and angle; and azimuthal canopy orientation. We report on the ability of some maize (Zea mays) genotypes to alter the orientations of their leaves during development in coordination with adjacent plants. Although the upper canopies of these genotypes retain the typical alternate-distichous phyllotaxy of maize, their leaves grow parallel to those of adjacent plants. A genome-wide association study (GWAS) on this parallel canopy trait identified candidate genes, many of which are associated with shade avoidance syndrome, including phytochromeC2. GWAS conducted on the fraction of photosynthetically active radiation (PAR) intercepted by canopies also identified multiple candidate genes, including liguleless1 (lg1), previously defined by its role in ligule development. Under high plant densities, mutants of shade avoidance syndrome and liguleless genes (lg1, lg2, and Lg3) exhibit altered canopy patterns, viz, the numbers of interrow leaves are greatly reduced as compared to those of nonmutant controls, resulting in dramatically decreased PAR interception. In at least the case of lg2, this phenotype is not a consequence of abnormal ligule development. Instead, liguleless gene functions are required for normal light responses, including azimuth canopy re-orientation.

Current challenges and future of agricultural genomes to phenomes in the USA.

Dramatic improvements in measuring genetic variation across agriculturally relevant populations (genomics) must be matched by improvements in identifying and measuring relevant trait variation in such populations across many environments (phenomics). Identifying the most critical opportunities and challenges in genome to phenome (G2P) research is the focus of this paper. Previously (Genome Biol, 23(1):1-11, 2022), we laid out how Agricultural Genome to Phenome Initiative (AG2PI) will coordinate activities with USA federal government agencies expand public-private partnerships, and engage with external stakeholders to achieve a shared vision of future the AG2PI. Acting on this latter step, AG2PI organized the "Thinking Big: Visualizing the Future of AG2PI" two-day workshop held September 9-10, 2022, in Ames, Iowa, co-hosted with the United State Department of Agriculture's National Institute of Food and Agriculture (USDA NIFA). During the meeting, attendees were asked to use their experience and curiosity to review the current status of agricultural genome to phenome (AG2P) work and envision the future of the AG2P field. The topic summaries composing this paper are distilled from two 1.5-h small group discussions. Challenges and solutions identified across multiple topics at the workshop were explored. We end our discussion with a vision for the future of agricultural progress, identifying two areas of innovation needed: (1) innovate in genetic improvement methods development and evaluation and (2) innovate in agricultural research processes to solve societal problems. To address these needs, we then provide six specific goals that we recommend be implemented immediately in support of advancing AG2P research.

A genetic tradeoff for tolerance to moderate and severe heat stress in US hybrid maize.

Global climate change is increasing both average temperatures and the frequencies of extreme high temperatures. Past studies have documented a strong negative effect of exposures to temperatures >30°C on hybrid maize yields. However, these studies could not disentangle genetic adaptation via artificial selection from changes in agronomic practices. Because most of the earliest maize hybrids are no longer available, side-by-side comparisons with modern hybrids under current field conditions are generally impossible. Here, we report on the collection and curation of 81 years of public yield trial records covering 4,730 maize hybrids, which enabled us to model genetic variation for temperature responses among maize hybrids. We show that selection may have indirectly and inconsistently contributed to the genetic adaptation of maize to moderate heat stress over this time period while preserving genetic variance for continued adaptation. However, our results reveal the existence of a genetic tradeoff for tolerance to moderate and severe heat stress, leading to a decrease in tolerance to severe heat stress over the same time period. Both trends are particularly conspicuous since the mid-1970s. Such a tradeoff poses challenges to the continued adaptation of maize to warming climates due to a projected increase in the frequency of extreme heat events. Nevertheless, given recent advances in phenomics, enviromics, and physiological modeling, our results offer a degree of optimism for the capacity of plant breeders to adapt maize to warming climates, assuming appropriate levels of R&D investment.

High-Throughput Field Plant Phenotyping: A Self-Supervised Sequential CNN Method to Segment Overlapping Plants.

High-throughput plant phenotyping-the use of imaging and remote sensing to record plant growth dynamics-is becoming more widely used. The first step in this process is typically plant segmentation, which requires a well-labeled training dataset to enable accurate segmentation of overlapping plants. However, preparing such training data is both time and labor intensive. To solve this problem, we propose a plant image processing pipeline using a self-supervised sequential convolutional neural network method for in-field phenotyping systems. This first step uses plant pixels from greenhouse images to segment nonoverlapping in-field plants in an early growth stage and then applies the segmentation results from those early-stage images as training data for the separation of plants at later growth stages. The proposed pipeline is efficient and self-supervising in the sense that no human-labeled data are needed. We then combine this approach with functional principal components analysis to reveal the relationship between the growth dynamics of plants and genotypes. We show that the proposed pipeline can accurately separate the pixels of foreground plants and estimate their heights when foreground and background plants overlap and can thus be used to efficiently assess the impact of treatments and genotypes on plant growth in a field environment by computer vision techniques. This approach should be useful for answering important scientific questions in the area of high-throughput phenotyping.

The 2020 derecho revealed limited overlap between maize genes associated with root lodging and root system architecture.

Roots anchor plants in soil, and the failure of anchorage (i.e. root lodging) is a major cause of crop yield loss. Anchorage is often assumed to be driven by root system architecture (RSA). We made use of a natural experiment to measure the overlap between the genetic regulation of RSA and anchorage. After one of the most devastating derechos ever recorded in August 2020, we phenotyped root lodging in a maize (Zea mays) diversity panel consisting of 369 genotypes grown in 6 environments affected by the derecho. Genome-wide and transcriptome-wide association studies identified 118 candidate genes associated with root lodging. Thirty-four percent (40/118) of these were homologs of genes from Arabidopsis (Arabidopsis thaliana) that affect traits such as root morphology and lignin content, expected to affect root lodging. Finally, gene ontology enrichment analysis of the candidate genes and their predicted interaction partners at the transcriptional and translational levels revealed the complex regulatory networks of physiological and biochemical pathways underlying root lodging in maize. Limited overlap between genes associated with lodging resistance and RSA in this diversity panel suggests that anchorage depends in part on factors other than the gross characteristics of RSA.

New insights into trait introgression with the look-ahead intercrossing strategy.

Trait introgression (TI) can be a time-consuming and costly task that typically requires multiple generations of backcrossing (BC). Usually, the aim is to introduce one or more alleles (e.g. QTLs) from a single donor into an elite recipient, both of which are fully inbred. This article studies the potential advantages of incorporating intercrossing (IC) into TI programs when compared with relying solely on the traditional BC framework. We simulate a TI breeding pipeline using 3 previously proposed selection strategies for the traditional BC scheme and 3 modified strategies that allow IC. Our proposed look-ahead intercrossing method (LAS-IC) combines look-ahead Monte Carlo simulations, intercrossing, and additional selection criteria to improve computational efficiency. We compared the efficiency of the 6 strategies across 5 levels of resource availability considering the generation when the major QTLs have been successfully introduced into the recipient and a desired background recovery rate reached. Simulations demonstrate that the inclusion of intercrossing in a TI program can substantially increase efficiency and the probability of success. The proposed LAS-IC provides the highest probability of success across the different scenarios using fewer resources compared with BC-only strategies.

The Anthocyanin Accumulation Related ZmBZ1, Facilitates Seedling Salinity Stress Tolerance via ROS Scavenging.

Anthocyanins are a class of antioxidants that scavenge free radicals in cells and play an important role in promoting human health and preventing many diseases. Here, we characterized a maize Bronze gene (BZ1) from the purple colored W22 introgression line, which encodes an anthocyanin 3-O-glucosyltransferase, a key enzyme in the anthocyanin synthesis pathway. Mutation of ZmBZ1 showed bronze-colored seeds and reduced anthocyanins in seeds aleurone layer, seedlings coleoptile, and stem of mature plants by comparison with purple colored W22 (WT). Furthermore, we proved that maize BZ1 is an aleurone layer-specific expressed protein and sub-located in cell nucleus. Real-time tracing of the anthocyanins in developing seeds demonstrated that the pigment was visible from 16 DAP (day after pollination) in field condition, and first deposited in the crown part then spread all over the seed. Additionally, it was transferred along with the embryo cell activity during seed germination, from aleurone layer to cotyledon and coleoptile, as confirmed by microscopy and real-time qRT-PCR. Finally, we demonstrated that the ZmBZ1 contributes to stress tolerance, especially salinity. Further study proved that ZmBZ1 participates in reactive oxygen scavenging (ROS) by accumulating anthocyanins, thereby enhancing the tolerance to abiotic stress.

Maize Leaf Appearance Rates: A Synthesis From the United States Corn Belt.

The relationship between collared leaf number and growing degree days (GDD) is crucial for predicting maize phenology. Biophysical crop models convert GDD accumulation to leaf numbers by using a constant parameter termed phyllochron (°C-day leaf-1) or leaf appearance rate (LAR; leaf oC-day-1). However, such important parameter values are rarely estimated for modern maize hybrids. To fill this gap, we sourced and analyzed experimental datasets from the United States Corn Belt with the objective to (i) determine phyllochron values for two types of models: linear (1-parameter) and bilinear (3-parameters; phase I and II phyllochron, and transition point) and (ii) explore whether environmental factors such as photoperiod and radiation, and physiological variables such as plant growth rate can explain variability in phyllochron and improve predictability of maize phenology. The datasets included different locations (latitudes between 48° N and 41° N), years (2009-2019), hybrids, and management settings. Results indicated that the bilinear model represented the leaf number vs. GDD relationship more accurately than the linear model (R 2 = 0.99 vs. 0.95, n = 4,694). Across datasets, first phase phyllochron, transition leaf number, and second phase phyllochron averaged 57.9 ± 7.5°C-day, 9.8 ± 1.2 leaves, and 30.9 ± 5.7°C-day, respectively. Correlation analysis revealed that radiation from the V3 to the V9 developmental stages had a positive relationship with phyllochron (r = 0.69), while photoperiod was positively related to days to flowering or total leaf number (r = 0.89). Additionally, a positive nonlinear relationship between maize LAR and plant growth rate was found. Present findings provide important parameter values for calibration and optimization of maize crop models in the United States Corn Belt, as well as new insights to enhance mechanisms in crop models.

KAT4IA: K-Means Assisted Training for Image Analysis of Field-Grown Plant Phenotypes.

High-throughput phenotyping enables the efficient collection of plant trait data at scale. One example involves using imaging systems over key phases of a crop growing season. Although the resulting images provide rich data for statistical analyses of plant phenotypes, image processing for trait extraction is required as a prerequisite. Current methods for trait extraction are mainly based on supervised learning with human labeled data or semisupervised learning with a mixture of human labeled data and unsupervised data. Unfortunately, preparing a sufficiently large training data is both time and labor-intensive. We describe a self-supervised pipeline (KAT4IA) that uses K-means clustering on greenhouse images to construct training data for extracting and analyzing plant traits from an image-based field phenotyping system. The KAT4IA pipeline includes these main steps: self-supervised training set construction, plant segmentation from images of field-grown plants, automatic separation of target plants, calculation of plant traits, and functional curve fitting of the extracted traits. To deal with the challenge of separating target plants from noisy backgrounds in field images, we describe a novel approach using row-cuts and column-cuts on images segmented by transform domain neural network learning, which utilizes plant pixels identified from greenhouse images to train a segmentation model for field images. This approach is efficient and does not require human intervention. Our results show that KAT4IA is able to accurately extract plant pixels and estimate plant heights.

Trajectories of Homoeolog-Specific Expression in Allotetraploid Tragopogon castellanus Populations of Independent Origins.

Polyploidization can have a significant ecological and evolutionary impact by providing substantially more genetic material that may result in novel phenotypes upon which selection may act. While the effects of polyploidization are broadly reviewed across the plant tree of life, the reproducibility of these effects within naturally occurring, independently formed polyploids is poorly characterized. The flowering plant genus Tragopogon (Asteraceae) offers a rare glimpse into the intricacies of repeated allopolyploid formation with both nascent (< 90 years old) and more ancient (mesopolyploids) formations. Neo- and mesopolyploids in Tragopogon have formed repeatedly and have extant diploid progenitors that facilitate the comparison of genome evolution after polyploidization across a broad span of evolutionary time. Here, we examine four independently formed lineages of the mesopolyploid Tragopogon castellanus for homoeolog expression changes and fractionation after polyploidization. We show that expression changes are remarkably similar among these independently formed polyploid populations with large convergence among expressed loci, moderate convergence among loci lost, and stochastic silencing. We further compare and contrast these results for T. castellanus with two nascent Tragopogon allopolyploids. While homoeolog expression bias was balanced in both nascent polyploids and T. castellanus, the degree of additive expression was significantly different, with the mesopolyploid populations demonstrating more non-additive expression. We suggest that gene dosage and expression noise minimization may play a prominent role in regulating gene expression patterns immediately after allopolyploidization as well as deeper into time, and these patterns are conserved across independent polyploid lineages.

Meta-analysis identifies pleiotropic loci controlling phenotypic trade-offs in sorghum.

Community association populations are composed of phenotypically and genetically diverse accessions. Once these populations are genotyped, the resulting marker data can be reused by different groups investigating the genetic basis of different traits. Because the same genotypes are observed and scored for a wide range of traits in different environments, these populations represent a unique resource to investigate pleiotropy. Here, we assembled a set of 234 separate trait datasets for the Sorghum Association Panel, a group of 406 sorghum genotypes widely employed by the sorghum genetics community. Comparison of genome-wide association studies (GWAS) conducted with two independently generated marker sets for this population demonstrate that existing genetic marker sets do not saturate the genome and likely capture only 35-43% of potentially detectable loci controlling variation for traits scored in this population. While limited evidence for pleiotropy was apparent in cross-GWAS comparisons, a multivariate adaptive shrinkage approach recovered both known pleiotropic effects of existing loci and new pleiotropic effects, particularly significant impacts of known dwarfing genes on root architecture. In addition, we identified new loci with pleiotropic effects consistent with known trade-offs in sorghum development. These results demonstrate the potential for mining existing trait datasets from widely used community association populations to enable new discoveries from existing trait datasets as new, denser genetic marker datasets are generated for existing community association populations.

Chromosome-level genome assembly of a regenerable maize inbred line A188.

BACKGROUND: The maize inbred line A188 is an attractive model for elucidation of gene function and improvement due to its high embryogenic capacity and many contrasting traits to the first maize reference genome, B73, and other elite lines. The lack of a genome assembly of A188 limits its use as a model for functional studies. RESULTS: Here, we present a chromosome-level genome assembly of A188 using long reads and optical maps. Comparison of A188 with B73 using both whole-genome alignments and read depths from sequencing reads identify approximately 1.1 Gb of syntenic sequences as well as extensive structural variation, including a 1.8-Mb duplication containing the Gametophyte factor1 locus for unilateral cross-incompatibility, and six inversions of 0.7 Mb or greater. Increased copy number of carotenoid cleavage dioxygenase 1 (ccd1) in A188 is associated with elevated expression during seed development. High ccd1 expression in seeds together with low expression of yellow endosperm 1 (y1) reduces carotenoid accumulation, accounting for the white seed phenotype of A188. Furthermore, transcriptome and epigenome analyses reveal enhanced expression of defense pathways and altered DNA methylation patterns of the embryonic callus. CONCLUSIONS: The A188 genome assembly provides a high-resolution sequence for a complex genome species and a foundational resource for analyses of genome variation and gene function in maize. The genome, in comparison to B73, contains extensive intra-species structural variations and other genetic differences. Expression and network analyses identify discrete profiles for embryonic callus and other tissues.

Identification and utilization of genetic determinants of trait measurement errors in image-based, high-throughput phenotyping.

The accuracy of trait measurements greatly affects the quality of genetic analyses. During automated phenotyping, trait measurement errors, i.e. differences between automatically extracted trait values and ground truth, are often treated as random effects that can be controlled by increasing population sizes and/or replication number. In contrast, there is some evidence that trait measurement errors may be partially under genetic control. Consistent with this hypothesis, we observed substantial nonrandom, genetic contributions to trait measurement errors for five maize (Zea mays) tassel traits collected using an image-based phenotyping platform. The phenotyping accuracy varied according to whether a tassel exhibited "open" versus. "closed" branching architecture, which is itself under genetic control. Trait-associated SNPs (TASs) identified via genome-wide association studies (GWASs) conducted on five tassel traits that had been phenotyped both manually (i.e. ground truth) and via feature extraction from images exhibit little overlap. Furthermore, identification of TASs from GWASs conducted on the differences between the two values indicated that a fraction of measurement error is under genetic control. Similar results were obtained in a sorghum (Sorghum bicolor) plant height dataset, demonstrating that trait measurement error is genetically determined in multiple species and traits. Trait measurement bias cannot be controlled by increasing population size and/or replication number.

TWAS results are complementary to and less affected by linkage disequilibrium than GWAS.

A genome-wide association study (GWAS) is used to identify genetic markers associated with phenotypic variation. In contrast, a transcriptome-wide association study (TWAS) detects associations between gene expression levels and phenotypic variation. It has previously been shown that in the cross-pollinated species, maize (Zea mays), GWAS, and TWAS identify complementary sets of trait-associated genes, many of which exhibit characteristics of true positives. Here, we extend this conclusion to the self-pollinated species, Arabidopsis thaliana and soybean (Glycine max). Linkage disequilibrium (LD) can result in the identification, via GWAS, of false-positive associations. In all three analyzed plant species, most trait-associated genes identified via TWAS are well separated physically from other candidate genes. Hence, TWAS is less affected by LD than is GWAS, demonstrating that TWAS is particularly well suited for association studies in genomes with slow rates of LD decay, such as soybean. TWAS is reasonably robust to the plant organs/tissues used to determine expression levels. In summary, this study confirms that TWAS is a promising approach for accurate gene-level association mapping in plants that is complementary to GWAS, and established that TWAS can exhibit substantial advantages relative to GWAS in species with slow rates of LD decay.

An integrated framework reinstating the environmental dimension for GWAS and genomic selection in crops.

Identifying mechanisms and pathways involved in gene-environment interplay and phenotypic plasticity is a long-standing challenge. It is highly desirable to establish an integrated framework with an environmental dimension for complex trait dissection and prediction. A critical step is to identify an environmental index that is both biologically relevant and estimable for new environments. With extensive field-observed complex traits, environmental profiles, and genome-wide single nucleotide polymorphisms for three major crops (maize, wheat, and oat), we demonstrated that identifying such an environmental index (i.e., a combination of environmental parameter and growth window) enables genome-wide association studies and genomic selection of complex traits to be conducted with an explicit environmental dimension. Interestingly, genes identified for two reaction-norm parameters (i.e., intercept and slope) derived from flowering time values along the environmental index were less colocalized for a diverse maize panel than for wheat and oat breeding panels, agreeing with the different diversity levels and genetic constitutions of the panels. In addition, we showcased the usefulness of this framework for systematically forecasting the performance of diverse germplasm panels in new environments. This general framework and the companion CERIS-JGRA analytical package should facilitate biologically informed dissection of complex traits, enhanced performance prediction in breeding for future climates, and coordinated efforts to enrich our understanding of mechanisms underlying phenotypic variation.

Detection of the Progression of Anthesis in Field-Grown Maize Tassels: A Case Study.

The tassel of the maize plant is responsible for the production and dispersal of pollen for subsequent capture by the silk (stigma) and fertilization of the ovules. Both the amount and timing of pollen shed are physiological traits that impact the production of a hybrid seed. This study describes an automated end-to-end pipeline that combines deep learning and image processing approaches to extract tassel flowering patterns from time-lapse camera images of plants grown under field conditions. Inbred lines from the SAM and NAM diversity panels were grown at the Curtiss farm at Iowa State University, Ames, IA, during the summer of 2016. Using a set of around 500 pole-mounted cameras installed in the field, images of plants were captured every 10 minutes of daylight hours over a three-week period. Extracting data from imaging performed under field conditions is challenging due to variabilities in weather, illumination, and the morphological diversity of tassels. To address these issues, deep learning algorithms were used for tassel detection, classification, and segmentation. Image processing approaches were then used to crop the main spike of the tassel to track reproductive development. The results demonstrated that deep learning with well-labeled data is a powerful tool for detecting, classifying, and segmenting tassels. Our sequential workflow exhibited the following metrics: mAP for tassel detection was 0.91, F1 score obtained for tassel classification was 0.93, and accuracy of semantic segmentation in creating a binary image from the RGB tassel images was 0.95. This workflow was used to determine spatiotemporal variations in the thickness of the main spike-which serves as a proxy for anthesis progression.