Genetic analysis of pericarp pigmentation variation in Corn Belt dent maize.

The US standard for maize commercially grown for grain specifies that yellow corn can contain at maximum 5% corn of other colors. Inbred parents of commercial hybrids typically have clear pericarp, but transgressive segregants in breeding populations can display variation in pericarp pigmentation. We identified 10 doubled haploid biparental populations segregating for pigmented pericarp and evaluated qualitative genetic models using chi-square tests of observed and expected frequencies. Pigmentation ranged from light to dark brown color, and pigmentation intensity was quantitatively measured across 1,327 inbred lines using hue calculated from RGB pixel values. Genetic mapping was used to identify loci associated with pigmentation intensity. For 9 populations, pigmentation inheritance best fit a hypothesis of a 2- or 3-gene epistatic model. Significant differences in pigment intensity were observed across populations. W606S-derived inbred lines with the darkest pericarp often had clear glumes, suggesting the presence of a novel P1-rw allele, a hypothesis supported by a significant quantitative trait locus peak at P1. A separate quantitative trait locus region on chromosome 2 between 221.64 and 226.66 Mbp was identified in LH82-derived populations, and the peak near p1 was absent. A genome-wide association study using 416 inbred lines from the Wisconsin Diversity panel with full genome resequencing revealed 4 significant associations including the region near P1. This study supports that pericarp pigmentation among dent maize inbreds can arise by transgressive segregation when pigmentation in the parental generation is absent and is partially explained by functional allelic variation at the P1 locus.

Leveraging data from the Genomes-to-Fields Initiative to investigate genotype-by-environment interactions in maize in North America.

Genotype-by-environment (G×E) interactions can significantly affect crop performance and stability. Investigating G×E requires extensive data sets with diverse cultivars tested over multiple locations and years. The Genomes-to-Fields (G2F) Initiative has tested maize hybrids in more than 130 year-locations in North America since 2014. Here, we curate and expand this data set by generating environmental covariates (using a crop model) for each of the trials. The resulting data set includes DNA genotypes and environmental data linked to more than 70,000 phenotypic records of grain yield and flowering traits for more than 4000 hybrids. We show how this valuable data set can serve as a benchmark in agricultural modeling and prediction, paving the way for countless G×E investigations in maize. We use multivariate analyses to characterize the data set's genetic and environmental structure, study the association of key environmental factors with traits, and provide benchmarks using genomic prediction models.

Genetic architecture of source-sink-regulated senescence in maize.

Source and sink interactions play a critical but mechanistically poorly understood role in the regulation of senescence. To disentangle the genetic and molecular mechanisms underlying source-sink-regulated senescence (SSRS), we performed a phenotypic, transcriptomic, and systems genetics analysis of senescence induced by the lack of a strong sink in maize (Zea mays). Comparative analysis of genotypes with contrasting SSRS phenotypes revealed that feedback inhibition of photosynthesis, a surge in reactive oxygen species, and the resulting endoplasmic reticulum (ER) stress were the earliest outcomes of weakened sink demand. Multienvironmental evaluation of a biparental population and a diversity panel identified 12 quantitative trait loci and 24 candidate genes, respectively, underlying SSRS. Combining the natural diversity and coexpression networks analyses identified 7 high-confidence candidate genes involved in proteolysis, photosynthesis, stress response, and protein folding. The role of a cathepsin B like protease 4 (ccp4), a candidate gene supported by systems genetic analysis, was validated by analysis of natural alleles in maize and heterologous analyses in Arabidopsis (Arabidopsis thaliana). Analysis of natural alleles suggested that a 700-bp polymorphic promoter region harboring multiple ABA-responsive elements is responsible for differential transcriptional regulation of ccp4 by ABA and the resulting variation in SSRS phenotype. We propose a model for SSRS wherein feedback inhibition of photosynthesis, ABA signaling, and oxidative stress converge to induce ER stress manifested as programed cell death and senescence. These findings provide a deeper understanding of signals emerging from loss of sink strength and offer opportunities to modify these signals to alter senescence program and enhance crop productivity.

Effects of endosperm type and storage length of whole-plant corn silage on nitrogen fraction, fermentation products, zein profile, and starch digestibility.

Zeins are commercially important proteins found in corn endosperms. The objective of this study was to evaluate the effect of altering zein levels in corn inbred lines carrying endosperm mutations with differential allelic dosage and analyze the effects on the composition, nutritive value, and starch digestibility of whole-plant corn silage (WPCS) at 5 storage lengths. Three inbred lines carrying 3 different endosperm modifiers (opaque-2 [o2], floury-2 [fl2], and soft endosperm-1 [h1]) were pollinated with 2 pollen sources to form pairs of near-isogenic lines with either 2 or 3 doses of the mutant allele for each endosperm modifier. The experiment was designed as a split-plot design with 3 replications. Pollinated genotype was the main plot factor, and storage length was the subplot-level factor. Agronomic precautions were taken to mimic hybrid WPCS to the extent possible. Samples were collected at approximately 30% dry matter (DM) using a forage harvester and ensiled in heat-sealed plastic bags for 0, 30, 60, 120, and 240 d. Thus, the experiment consisted of 30 treatments (6 genotypes × 5 storage lengths) and 90 ensiling units (3 replications per treatment). Measurements included nutrient analysis, including crude protein, soluble crude protein, amylase-treated neutral detergent fiber, acid detergent fiber, lignin, starch, fermentation end products, zein concentration, and in vitro starch digestibility (ivSD). The nutritional profile of the inbred-based silage samples was similar to hybrid values reported in literature. Significant differences were found in fresh (unfermented) sample kernels for endosperm vitreousness and zein profiles between and within isogenic pairs. The o2 homozygous (3 doses of mutant allele) had the highest reduction in vitreousness level (74.5 to 38%) and zein concentration (6.2 to 4.7% of DM) compared with the heterozygous counterpart (2 doses of mutant allele). All genotypes showed significant reduction of total zeins and α-zeins during progressive storage length. In vitro starch digestibility increased with storage length and had significant effects of genotype and storage length but not for genotype by storage length interaction, which suggests that the storage period did not attenuate the difference in ivSD between near-isogenic pairs caused by zeins in WPCS. Both total zeins and α-zeins showed a strong negative correlation with ivSD, which agrees with the general hypothesis that the degradation of zeins increases ruminal starch degradability. Homozygous o2 was the only mutant with significantly higher ivSD compared with the heterozygous version, which suggests that, if all other conditions remain constant in a WPCS systems, substantial reductions in endosperm α-zeins are required to significantly improve ivSD in the silo.

A key to totipotency: Wuschel-like homeobox 2a unlocks embryogenic culture response in maize (Zea mays L.).

The ability of plant somatic cells to dedifferentiate, form somatic embryos and regenerate whole plants in vitro has been harnessed for both clonal propagation and as a key component of plant genetic engineering systems. Embryogenic culture response is significantly limited, however, by plant genotype in most species. This impedes advancements in both plant transformation-based functional genomics research and crop improvement efforts. We utilized natural variation among maize inbred lines to genetically map somatic embryo generation potential in tissue culture and identify candidate genes underlying totipotency. Using a series of maize lines derived from crosses involving the culturable parent A188 and the non-responsive parent B73, we identified a region on chromosome 3 associated with embryogenic culture response and focused on three candidate genes within the region based on genetic position and expression pattern. Two candidate genes showed no effect when ectopically expressed in B73, but the gene Wox2a was found to induce somatic embryogenesis and embryogenic callus proliferation. Transgenic B73 cells with strong constitutive expression of the B73 and A188 coding sequences of Wox2a were found to produce somatic embryos at similar frequencies, demonstrating that sufficient expression of either allele could rescue the embryogenic culture phenotype. Transgenic B73 plants were regenerated from the somatic embryos without chemical selection and no pleiotropic effects were observed in the Wox2a overexpression lines in the regenerated T0 plants or in the two independent events which produced T1 progeny. In addition to linking natural variation in tissue culture response to Wox2a, our data support the utility of Wox2a in enabling transformation of recalcitrant genotypes.

Transcription factor bHLH121 regulates root cortical aerenchyma formation in maize.

Root anatomical phenotypes present a promising yet underexploited avenue to deliver major improvements in yield and climate resilience of crops by improving water and nutrient uptake. For instance, the formation of root cortical aerenchyma (RCA) significantly increases soil exploration and resource capture by reducing the metabolic costs of root tissue. A key bottleneck in studying such phenotypes has been the lack of robust high-throughput anatomical phenotyping platforms. We exploited a phenotyping approach based on laser ablation tomography, termed Anatomics, to quantify variation in RCA formation of 436 diverse maize lines in the field. Results revealed a significant and heritable variation for RCA formation. Genome-wide association studies identified a single-nucleotide polymorphism mapping to a root cortex-expressed gene-encoding transcription factor bHLH121. Functional studies identified that the bHLH121 Mu transposon mutant line and CRISPR/Cas9 loss-of-function mutant line showed reduced RCA formation, whereas an overexpression line exhibited significantly greater RCA formation when compared to the wild-type line. Characterization of these lines under suboptimal water and nitrogen availability in multiple soil environments revealed that bHLH121 is required for RCA formation developmentally as well as under studied abiotic stress. Overall functional validation of the bHLH121 gene's importance in RCA formation provides a functional marker to select varieties with improved soil exploration and thus yield under suboptimal conditions.

A novel high-throughput hyperspectral scanner and analytical methods for predicting maize kernel composition and physical traits.

Large-scale investigations of maize kernel traits important to researchers, breeders, and processors require high throughput methods, which are presently lacking. To address this bottleneck, we developed a novel flatbed platform that automatically acquires and analyzes multiwavelength near-infrared (NIR hyperspectral) images of maize kernels precisely enough to support robust predictions of protein content, density, and endosperm vitreousness. The upward facing-camera design and the automated ability to analyze the embryo or abgerminal sides of each individual kernel in a sample with the appropriate side-specific model helped to produce a superior combination of throughput and prediction accuracy compared to other single-kernel platforms. Protein was predicted to within 0.85% (root mean square error of prediction), density to within 0.038 g/cm3, and endosperm vitreousness percentage to within 6.3%. Kernel length and width were also accurately measured so that each kernel in a rapidly scanned sample was comprehensively characterized.

Genetic mapping and prediction of flowering time and plant height in a maize Stiff Stalk MAGIC population.

The Stiff Stalk heterotic pool is a foundation of US maize seed parent germplasm and has been heavily utilized by both public and private maize breeders since its inception in the 1930s. Flowering time and plant height are critical characteristics for both inbred parents and their test crossed hybrid progeny. To study these traits, a 6-parent multiparent advanced generation intercross population was developed including maize inbred lines B73, B84, PHB47 (B37 type), LH145 (B14 type), PHJ40 (novel early Stiff Stalk), and NKH8431 (B73/B14 type). A set of 779 doubled haploid lines were evaluated for flowering time and plant height in 2 field replicates in 2016 and 2017, and a subset of 689 and 561 doubled haploid lines were crossed to 2 testers, respectively, and evaluated as hybrids in 2 locations in 2018 and 2019 using an incomplete block design. Markers were derived from a practical haplotype graph built from the founder whole genome assemblies and genotype-by-sequencing and exome capture-based sequencing of the population. Genetic mapping utilizing an update to R/qtl2 revealed differing profiles of significant loci for both traits between 635 of the DH lines and 2 sets of 570 and 471 derived hybrids. Genomic prediction was used to test the feasibility of predicting hybrid phenotypes based on the per se data. Predictive abilities were highest on direct models trained using the data they would predict (0.55-0.63), and indirect models trained using per se data to predict hybrid traits had slightly lower predictive abilities (0.49-0.55). Overall, this finding is consistent with the overlapping and nonoverlapping significant quantitative trait loci found within the per se and hybrid populations and suggests that selections for phenology traits can be made effectively on doubled haploid lines before hybrid data is available.

Root angle in maize influences nitrogen capture and is regulated by calcineurin B-like protein (CBL)-interacting serine/threonine-protein kinase 15 (ZmCIPK15).

Crops with reduced nutrient and water requirements are urgently needed in global agriculture. Root growth angle plays an important role in nutrient and water acquisition. A maize diversity panel of 481 genotypes was screened for variation in root angle employing a high-throughput field phenotyping platform. Genome-wide association mapping identified several single nucleotide polymorphisms (SNPs) associated with root angle, including one located in the root expressed CBL-interacting serine/threonine-protein kinase 15 (ZmCIPK15) gene (LOC100285495). Reverse genetic studies validated the functional importance of ZmCIPK15, causing a approximately 10° change in root angle in specific nodal positions. A steeper root growth angle improved nitrogen capture in silico and in the field. OpenSimRoot simulations predicted at 40 days of growth that this change in angle would improve nitrogen uptake by 11% and plant biomass by 4% in low nitrogen conditions. In field studies under suboptimal N availability, the cipk15 mutant with steeper growth angles had 18% greater shoot biomass and 29% greater shoot nitrogen accumulation compared to the wild type after 70 days of growth. We propose that a steeper root growth angle modulated by ZmCIPK15 will facilitate efforts to develop new crop varieties with optimal root architecture for improved performance under edaphic stress.

Genetic loci associated with winter survivorship in diverse lowland switchgrass populations.

High winter mortality limits biomass yield of lowland switchgrass (Panicum virgatum L.) planted in the northern latitudes of North America. Breeding of cold tolerant switchgrass cultivars requires many years due to its perennial growth habit and the unpredictable winter selection pressure that is required to identify winter-hardy individuals. Identification of causal genetic variants for winter survivorship would accelerate the improvement of switchgrass biomass production. The objective of this study was to identify allelic variation associated with winter survivorship in lowland switchgrass populations using bulk segregant analysis (BSA). Twenty-nine lowland switchgrass populations were evaluated for winter survival at two locations in southern Wisconsin and 21 populations with differential winter survivorship were used for BSA. A maximum of 10% of the individuals (8-20) were bulked to create survivor and nonsurvivor DNA pools from each population and location. The DNA pools were evaluated using exome capture sequencing, and allele frequencies were used to conduct statistical tests. The BSA tests revealed nine quatitative trait loci (QTL) from tetraploid populations and seven QTL from octoploid populations. Many QTL were population-specific, but some were identified in multiple populations that originated across a broad geographic landscape. Four QTL (at positions 88 Mb on chromosome 2N, 115 Mb on chromosome 5K, and 1 and 100 Mb on chromosome 9N) were potentially the most useful QTL. Markers associated with winter survivorship in this study can be used to accelerate breeding cycles of lowland switchgrass populations and should lead to improvements in adaptation within USDA hardiness zones 4 and 5.

Genomic variation within the maize stiff-stalk heterotic germplasm pool.

The stiff-stalk heterotic group in Maize (Zea mays L.) is an important source of inbreds used in U.S. commercial hybrid production. Founder inbreds B14, B37, B73, and, to a lesser extent, B84, are found in the pedigrees of a majority of commercial seed parent inbred lines. We created high-quality genome assemblies of B84 and four expired Plant Variety Protection (ex-PVP) lines LH145 representing B14, NKH8431 of mixed descent, PHB47 representing B37, and PHJ40, which is a Pioneer Hi-Bred International (PHI) early stiff-stalk type. Sequence was generated using long-read sequencing achieving highly contiguous assemblies of 2.13-2.18 Gbp with N50 scaffold lengths >200 Mbp. Inbred-specific gene annotations were generated using a core five-tissue gene expression atlas, whereas transposable element (TE) annotation was conducted using de novo and homology-directed methodologies. Compared with the reference inbred B73, synteny analyses revealed extensive collinearity across the five stiff-stalk genomes, although unique components of the maize pangenome were detected. Comparison of this set of stiff-stalk inbreds with the original Iowa Stiff Stalk Synthetic breeding population revealed that these inbreds represent only a proportion of variation in the original stiff-stalk pool and there are highly conserved haplotypes in released public and ex-Plant Variety Protection inbreds. Despite the reduction in variation from the original stiff-stalk population, substantial genetic and genomic variation was identified supporting the potential for continued breeding success in this pool. The assemblies described here represent stiff-stalk inbreds that have historical and commercial relevance and provide further insight into the emerging maize pangenome.

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.

The importance of dominance and genotype-by-environment interactions on grain yield variation in a large-scale public cooperative maize experiment.

High-dimensional and high-throughput genomic, field performance, and environmental data are becoming increasingly available to crop breeding programs, and their integration can facilitate genomic prediction within and across environments and provide insights into the genetic architecture of complex traits and the nature of genotype-by-environment interactions. To partition trait variation into additive and dominance (main effect) genetic and corresponding genetic-by-environment variances, and to identify specific environmental factors that influence genotype-by-environment interactions, we curated and analyzed genotypic and phenotypic data on 1918 maize (Zea mays L.) hybrids and environmental data from 65 testing environments. For grain yield, dominance variance was similar in magnitude to additive variance, and genetic-by-environment variances were more important than genetic main effect variances. Models involving both additive and dominance relationships best fit the data and modeling unique genetic covariances among all environments provided the best characterization of the genotype-by-environment interaction patterns. Similarity of relative hybrid performance among environments was modeled as a function of underlying weather variables, permitting identification of weather covariates driving correlations of genetic effects across environments. The resulting models can be used for genomic prediction of mean hybrid performance across populations of environments tested or for environment-specific predictions. These results can also guide efforts to incorporate high-throughput environmental data into genomic prediction models and predict values in new environments characterized with the same environmental characteristics.

Single-parent expression drives dynamic gene expression complementation in maize hybrids.

Single-parent expression (SPE) is defined as gene expression in only one of the two parents. SPE can arise from differential expression between parental alleles, termed non-presence/absence (non-PAV) SPE, or from the physical absence of a gene in one parent, termed PAV SPE. We used transcriptome data of diverse Zea mays (maize) inbreds and hybrids, including 401 samples from five different tissues, to test for differences between these types of SPE genes. Although commonly observed, SPE is highly genotype and tissue specific. A positive correlation was observed between the genetic distance of the two inbred parents and the number of SPE genes identified. Regulatory analysis showed that PAV SPE and non-PAV SPE genes are mainly regulated by cis effects, with a small fraction under trans regulation. Polymorphic transposable element insertions in promoter sequences contributed to the high level of cis regulation for PAV SPE and non-PAV SPE genes. PAV SPE genes were more frequently expressed in hybrids than non-PAV SPE genes. The expression of parentally silent alleles in hybrids of non-PAV SPE genes was relatively rare but occurred in most hybrids. Non-PAV SPE genes with expression of the silent allele in hybrids are more likely to exhibit above high parent expression level than hybrids that do not express the silent allele, leading to non-additive expression. This study provides a comprehensive understanding of the nature of non-PAV SPE and PAV SPE genes and their roles in gene expression complementation in maize hybrids.

Characterizing introgression-by-environment interactions using maize near isogenic lines.

KEY MESSAGE: Significant introgression-by-environment interactions are observed for traits throughout development from small introgressed segments of the genome. Relatively small genomic introgressions containing quantitative trait loci can have significant impacts on the phenotype of an individual plant. However, the magnitude of phenotypic effects for the same introgression can vary quite substantially in different environments due to introgression-by-environment interactions. To study potential patterns of introgression-by-environment interactions, fifteen near-isogenic lines (NILs) with > 90% B73 genetic background and multiple Mo17 introgressions were grown in 16 different environments. These environments included five geographical locations with multiple planting dates and multiple planting densities. The phenotypic impact of the introgressions was evaluated for up to 26 traits that span different growth stages in each environment to assess introgression-by-environment interactions. Results from this study showed that small portions of the genome can drive significant genotype-by-environment interaction across a wide range of vegetative and reproductive traits, and the magnitude of the introgression-by-environment interaction varies across traits. Some introgressed segments were more prone to introgression-by-environment interaction than others when evaluating the interaction on a whole plant basis throughout developmental time, indicating variation in phenotypic plasticity throughout the genome. Understanding the profile of introgression-by-environment interaction in NILs is useful in consideration of how small introgressions of QTL or transgene containing regions might be expected to impact traits in diverse environments.

Assessing the Viability of Recovery of Hydroxycinnamic Acids from Lignocellulosic Biorefinery Alkaline Pretreatment Waste Streams.

Invited for this month's cover is the research team from the D.O.E. Great Lake Bioenergy Research Center (GLBRC) at the University of Wisconsin-Madison. The cover image shows how a diverse team with expertise in many different fields works together in an integrated fashion to address complex problems. Only when the whole system, from field to the liquid fuels and co-products, is assessed, can we identify the key parameters needed to design an economically viable biorefinery-based economy. Cover art by Chelsea Mamott. The Full Paper itself is available at 10.1002/cssc.201903345.

Maize genomes to fields (G2F): 2014-2017 field seasons: genotype, phenotype, climatic, soil, and inbred ear image datasets.

OBJECTIVES: Advanced tools and resources are needed to efficiently and sustainably produce food for an increasing world population in the context of variable environmental conditions. The maize genomes to fields (G2F) initiative is a multi-institutional initiative effort that seeks to approach this challenge by developing a flexible and distributed infrastructure addressing emerging problems. G2F has generated large-scale phenotypic, genotypic, and environmental datasets using publicly available inbred lines and hybrids evaluated through a network of collaborators that are part of the G2F's genotype-by-environment (G × E) project. This report covers the public release of datasets for 2014-2017. DATA DESCRIPTION: Datasets include inbred genotypic information; phenotypic, climatic, and soil measurements and metadata information for each testing location across years. For a subset of inbreds in 2014 and 2015, yield component phenotypes were quantified by image analysis. Data released are accompanied by README descriptions. For genotypic and phenotypic data, both raw data and a version without outliers are reported. For climatic data, a version calibrated to the nearest airport weather station and a version without outliers are reported. The 2014 and 2015 datasets are updated versions from the previously released files [1] while 2016 and 2017 datasets are newly available to the public.

Assessing the Viability of Recovery of Hydroxycinnamic Acids from Lignocellulosic Biorefinery Alkaline Pretreatment Waste Streams.

The hydroxycinnamic acids p-coumaric acid (pCA) and ferulic acid (FA) add diversity to the portfolio of products produced by using grass-fed lignocellulosic biorefineries. The level of lignin-bound pCA in Zea mays was modified by the alteration of p-coumaroyl-CoA monolignol transferase expression. The biomass was processed in a lab-scale alkaline-pretreatment biorefinery process and the data were used for a baseline technoeconomic analysis to determine where to direct future research efforts to couple plant design to biomass utilization processes. It is concluded that future plant engineering efforts should focus on strategies that ramp up accumulation of one type of hydroxycinnamate (pCA or FA) predominantly and suppress that of the other. Technoeconomic analysis indicates that target extraction titers of one hydroxycinnamic acid need to be >50 g kg-1 biomass, at least five times higher than observed titers for the impure pCA/FA product mixture from wild-type maize. The technical challenge for process engineers is to develop a viable process that requires more than 80 % reduction of the isolation costs.