Predicting moisture content during maize nixtamalization using machine learning with NIR spectroscopy.

KEY MESSAGE: Moisture content during nixtamalization can be accurately predicted from NIR spectroscopy when coupled with a support vector machine (SVM) model, is strongly modulated by the environment, and has a complex genetic architecture. Lack of high-throughput phenotyping systems for determining moisture content during the maize nixtamalization cooking process has led to difficulty in breeding for this trait. This study provides a high-throughput, quantitative measure of kernel moisture content during nixtamalization based on NIR scanning of uncooked maize kernels. Machine learning was utilized to develop models based on the combination of NIR spectra and moisture content determined from a scaled-down benchtop cook method. A linear support vector machine (SVM) model with a Spearman's rank correlation coefficient of 0.852 between wet laboratory and predicted values was developed from 100 diverse temperate genotypes grown in replicate across two environments. This model was applied to NIR spectra data from 501 diverse temperate genotypes grown in replicate in five environments. Analysis of variance revealed environment explained the highest percent of the variation (51.5%), followed by genotype (15.6%) and genotype-by-environment interaction (11.2%). A genome-wide association study identified 26 significant loci across five environments that explained between 5.04% and 16.01% (average = 10.41%). However, genome-wide markers explained 10.54% to 45.99% (average = 31.68%) of the variation, indicating the genetic architecture of this trait is likely complex and controlled by many loci of small effect. This study provides a high-throughput method to evaluate moisture content during nixtamalization that is feasible at the scale of a breeding program and provides important information about the factors contributing to variation of this trait for breeders and food companies to make future strategies to improve this important processing trait.

ZmCCT and the genetic basis of day-length adaptation underlying the postdomestication spread of maize.

Teosinte, the progenitor of maize, is restricted to tropical environments in Mexico and Central America. The pre-Columbian spread of maize from its center of origin in tropical Southern Mexico to the higher latitudes of the Americas required postdomestication selection for adaptation to longer day lengths. Flowering time of teosinte and tropical maize is delayed under long day lengths, whereas temperate maize evolved a reduced sensitivity to photoperiod. We measured flowering time of the maize nested association and diverse association mapping panels in the field under both short and long day lengths, and of a maize-teosinte mapping population under long day lengths. Flowering time in maize is a complex trait affected by many genes and the environment. Photoperiod response is one component of flowering time involving a subset of flowering time genes whose effects are strongly influenced by day length. Genome-wide association and targeted high-resolution linkage mapping identified ZmCCT, a homologue of the rice photoperiod response regulator Ghd7, as the most important gene affecting photoperiod response in maize. Under long day lengths ZmCCT alleles from diverse teosintes are consistently expressed at higher levels and confer later flowering than temperate maize alleles. Many maize inbred lines, including some adapted to tropical regions, carry ZmCCT alleles with no sensitivity to day length. Indigenous farmers of the Americas were remarkably successful at selecting on genetic variation at key genes affecting the photoperiod response to create maize varieties adapted to vastly diverse environments despite the hindrance of the geographic axis of the Americas and the complex genetic control of flowering time.

Genetic architecture of maize kernel composition in the nested association mapping and inbred association panels.

The maize (Zea mays) kernel plays a critical role in feeding humans and livestock around the world and in a wide array of industrial applications. An understanding of the regulation of kernel starch, protein, and oil is needed in order to manipulate composition to meet future needs. We conducted joint-linkage quantitative trait locus mapping and genome-wide association studies (GWAS) for kernel starch, protein, and oil in the maize nested association mapping population, composed of 25 recombinant inbred line families derived from diverse inbred lines. Joint-linkage mapping revealed that the genetic architecture of kernel composition traits is controlled by 21-26 quantitative trait loci. Numerous GWAS associations were detected, including several oil and starch associations in acyl-CoA:diacylglycerol acyltransferase1-2, a gene that regulates oil composition and quantity. Results from nested association mapping were verified in a 282 inbred association panel using both GWAS and candidate gene association approaches. We identified many beneficial alleles that will be useful for improving kernel starch, protein, and oil content.

Iron bioavailability of maize (Zea mays L.) after removing the germ fraction.

Maize is a staple food for many communities with high levels of iron deficiency anemia. Enhancing the iron concentrations and iron bioavailability of maize with traditional breeding practices, especially after cooking and processing, could help alleviate iron deficiency in many of these regions. Previous studies on a small number of maize genotypes and maize flour products indicated that degermination (germ fraction removed with processing) could improve the iron bioavailability of maize. This study expanded upon this research by evaluating the iron bioavailability, mineral concentrations, and phytate concentrations of 52 diverse maize genotypes before (whole kernels) and after degermination. Whole and degerminated maize samples were cooked, dried, and milled to produce corn flour. Iron bioavailability was evaluated with an in vitro digestion Caco2 cell bioassay. In 30 of the maize genotypes, bioavailable iron increased when degerminated, thus indicating a higher fractional iron uptake because the iron concentrations decreased by more than 70% after the germ fraction was removed. The remaining 22 genotypes showed no change or a decrease in iron bioavailability after degermination. These results confirm previous research showing that the germ fraction is a strong inhibitory component for many maize varieties. Phytate concentrations in maize flours were greatly reduced with degermination. However, the relationship between phytate concentrations and the iron bioavailability of processed maize flour is complex, acting as either inhibitor or promoter of iron uptake depending on the color of the maize kernels and processing method used to produce flour. Other factors in the maize endosperm fractions are likely involved in the effects of degermination on iron bioavailability, such as vitreous or floury endosperm compositions and the polyphenol content of the bran. This study demonstrates that iron nutrition from maize can be enhanced by selecting genotypes where the inhibitory effect of the bran color and endosperm fraction are relatively low, especially after processing via degermination.

Demonstration of local adaptation in maize landraces by reciprocal transplantation.

Populations are locally adapted when they exhibit higher fitness than foreign populations in their native habitat. Maize landrace adaptations to highland and lowland conditions are of interest to researchers and breeders. To determine the prevalence and strength of local adaptation in maize landraces, we performed a reciprocal transplant experiment across an elevational gradient in Mexico. We grew 120 landraces, grouped into four populations (Mexican Highland, Mexican Lowland, South American Highland, South American Lowland), in Mexican highland and lowland common gardens and collected phenotypes relevant to fitness and known highland-adaptive traits such as anthocyanin pigmentation and macrohair density. 67k DArTseq markers were generated from field specimens to allow comparisons between phenotypic patterns and population genetic structure. We found phenotypic patterns consistent with local adaptation, though these patterns differ between the Mexican and South American populations. Quantitative trait differentiation (Q ST) was greater than neutral allele frequency differentiation (F ST) for many traits, signaling directional selection between pairs of populations. All populations exhibited higher fitness metric values when grown at their native elevation, and Mexican landraces had higher fitness than South American landraces when grown in these Mexican sites. As environmental distance between landraces' native collection sites and common garden sites increased, fitness values dropped, suggesting landraces are adapted to environmental conditions at their natal sites. Correlations between fitness and anthocyanin pigmentation and macrohair traits were stronger in the highland site than the lowland site, supporting their status as highland-adaptive. These results give substance to the long-held presumption of local adaptation of New World maize landraces to elevation and other environmental variables across North and South America.

Genetic control of kernel compositional variation in a maize diversity panel.

Maize (Zea mays L.) is a multi-purpose row crop grown worldwide, which, over time, has often been bred for increased yield at the detriment of lower composition grain quality. Some knowledge of the genetic factors that affect quality traits has been discovered through the study of classical maize mutants; however, much of the underlying genetic control of these traits and the interaction between these traits remains unknown. To better understand variation that exists for grain compositional traits in maize, we evaluated 501 diverse temperate maize inbred lines in five unique environments and predicted 16 compositional traits (e.g., carbohydrates, protein, and starch) based on the output of near-infrared (NIR) spectroscopy. Phenotypic analysis found substantial variation for compositional traits and the majority of variation was explained by genetic and environmental factors. Correlations and trade-offs among traits in different maize types (e.g., dent, sweetcorn, and popcorn) were explored, and significant differences and meaningful correlations were detected. In total, 22.9-71.0% of the phenotypic variation across these traits could be explained using 2,386,666 single nucleotide polymorphism (SNP) markers generated from whole-genome resequencing data. A genome-wide association study (GWAS) was conducted using these same markers and found 72 statistically significant SNPs for 11 compositional traits. This study provides valuable insights in the phenotypic variation and genetic control underlying compositional traits that can be used in breeding programs for improving maize grain quality.

Wide variability in kernel composition, seed characteristics, and zein profiles among diverse maize inbreds, landraces, and teosinte.

All crop species have been domesticated from their wild relatives, and geneticists are just now beginning to understand the consequences of artificial (human) selection on agronomic traits that are relevant today. The primary consequence is a basal loss of diversity across the genome, and an additional reduction in diversity for genes underlying traits targeted by selection. An understanding of attributes of the wild relatives may provide insight into target traits and valuable allelic variants for modern agriculture. This is especially true for maize (Zea mays ssp. mays), where its wild ancestor, teosinte (Z. mays ssp. parviglumis), is so strikingly different than modern maize. One obvious target of selection is the size and composition of the kernel. We evaluated kernel characteristics, kernel composition, and zein profiles for a diverse set of modern inbred lines, teosinte accessions, and landraces, the intermediate between inbreds and teosinte. We found that teosinte has very small seeds, but twice the protein content of landraces and inbred lines. Teosinte has a higher average alpha zein content (nearly 89% of total zeins as compared to 72% for inbred lines and 76% for landraces), and there are many novel alcohol-soluble proteins in teosinte relative to the other two germplasm groups. Nearly every zein protein varied in abundance among the germplasm groups, especially the methionine-rich delta zein protein, and the gamma zeins. Teosinte and landraces harbor phenotypic variation that will facilitate genetic dissection of kernel traits and grain quality, ultimately leading to improvement via traditional plant breeding and/or genetic engineering.

Conventional screening overlooks resistance sources: rootworm damage of diverse inbred lines and their B73 hybrids is unrelated.

The western corn rootworm, Diabrotica virgifera virgifera (LeConte), is a major pest of maize, Zea mays L., in the United States and Europe, and it is likely to increase in importance as a trend toward increased nonrotated maize favors larger rootworm populations. Options for rootworm management in nonrotated maize in Europe and in nontransgenic "refuge" areas in countries that permit transgenic maize are limited to insecticides. Development of additional options for growers would be helpful. Screening maize germplasm (e.g., landraces, populations, inbreds) for native resistance to western corn rootworm is labor-intensive and is usually conducted on unfinished germplasm and not on hybrid materials. However, we have recently observed that topcrossed (hybrid) materials tend to have reduced western corn rootworm damage. To formally test whether rootworm damage to inbreds and associated hybrids were correlated, we evaluated 25 diverse inbred lines and their B73 hybrids for western corn rootworm damage in seven environments. Overall, hybrids had significantly less damage than inbreds, but unfortunately, the correlation between inbreds and hybrids was not significant. These findings have important implications regarding screening germplasm for western corn rootworm resistance, namely, that inbred materials and perhaps populations should be topcrossed to form hybrid materials before screening for western corn rootworm damage to ensure that valuable sources of resistance to western corn rootworm are not missed during the screening process.

Structure of linkage disequilibrium in plants.

Future advances in plant genomics will make it possible to scan a genome for polymorphisms associated with qualitative and quantitative traits. Before this potential can be realized, we must understand the nature of linkage disequilibrium (LD) within a genome. LD, the nonrandom association of alleles at different loci, plays an integral role in association mapping, and determines the resolution of an association study. Recently, association mapping has been exploited to dissect quantitative trait loci (QTL). With the exception of maize and Arabidopsis, little research has been conducted on LD in plants. The mating system of the species (selfing versus outcrossing), and phenomena such as population structure and recombination hot spots, can strongly influence patterns of LD. The basic patterns of LD in plants will be better understood as more species are analyzed.

Genetic Analysis of Teosinte Alleles for Kernel Composition Traits in Maize.

Teosinte (Zea mays ssp. parviglumis) is the wild ancestor of modern maize (Zea mays ssp. mays). Teosinte contains greater genetic diversity compared with maize inbreds and landraces, but its use is limited by insufficient genetic resources to evaluate its value. A population of teosinte near isogenic lines (NILs) was previously developed to broaden the resources for genetic diversity of maize, and to discover novel alleles for agronomic and domestication traits. The 961 teosinte NILs were developed by backcrossing 10 geographically diverse parviglumis accessions into the B73 (reference genome inbred) background. The NILs were grown in two replications in 2009 and 2010 in Columbia, MO and Aurora, NY, respectively, and near infrared reflectance spectroscopy and nuclear magnetic resonance calibrations were developed and used to rapidly predict total kernel starch, protein, and oil content on a dry matter basis in bulk whole grains of teosinte NILs. Our joint-linkage quantitative trait locus (QTL) mapping analysis identified two starch, three protein, and six oil QTL, which collectively explained 18, 23, and 45% of the total variation, respectively. A range of strong additive allelic effects for kernel starch, protein, and oil content were identified relative to the B73 allele. Our results support our hypothesis that teosinte harbors stronger alleles for kernel composition traits than maize, and that teosinte can be exploited for the improvement of kernel composition traits in modern maize germplasm.

Genetic Analysis of Kernel Traits in Maize-Teosinte Introgression Populations.

Seed traits have been targeted by human selection during the domestication of crop species as a way to increase the caloric and nutritional content of food during the transition from hunter-gather to early farming societies. The primary seed trait under selection was likely seed size/weight as it is most directly related to overall grain yield. Additional seed traits involved in seed shape may have also contributed to larger grain. Maize (Zea mays ssp. mays) kernel weight has increased more than 10-fold in the 9000 years since domestication from its wild ancestor, teosinte (Z. mays ssp. parviglumis). In order to study how size and shape affect kernel weight, we analyzed kernel morphometric traits in a set of 10 maize-teosinte introgression populations using digital imaging software. We identified quantitative trait loci (QTL) for kernel area and length with moderate allelic effects that colocalize with kernel weight QTL. Several genomic regions with strong effects during maize domestication were detected, and a genetic framework for kernel traits was characterized by complex pleiotropic interactions. Our results both confirm prior reports of kernel domestication loci and identify previously uncharacterized QTL with a range of allelic effects, enabling future research into the genetic basis of these traits.

Expanding Maize Genetic Resources with Predomestication Alleles: Maize-Teosinte Introgression Populations.

Teosinte ( subsp. H. H. Iltis & Doebley) has greater genetic diversity than maize inbreds and landraces ( subsp. ). There are, however, limited genetic resources to efficiently evaluate and tap this diversity. To broaden resources for genetic diversity studies in maize, we developed and evaluated 928 near-isogenic introgression lines (NILs) from 10 teosinte accessions in the B73 background. Joint linkage analysis of the 10 introgression populations identified several large-effect quantitative trait loci (QTL) for days to anthesis (DTA), kernel row number (KRN), and 50-kernel weight (Wt50k). Our results confirm prior reports of kernel domestication loci and identify previously uncharacterized QTL with a range of allelic effects enabling future research into the genetic basis of these traits. Additionally, we used a targeted set of NILs to validate the effects of a KRN QTL located on chromosome 2. These introgression populations offer novel tools for QTL discovery and validation as well as a platform for initiating fine mapping.

The genetic architecture of maize height.

Height is one of the most heritable and easily measured traits in maize (Zea mays L.). Given a pedigree or estimates of the genomic identity-by-state among related plants, height is also accurately predictable. But, mapping alleles explaining natural variation in maize height remains a formidable challenge. To address this challenge, we measured the plant height, ear height, flowering time, and node counts of plants grown in >64,500 plots across 13 environments. These plots contained >7300 inbreds representing most publically available maize inbreds in the United States and families of the maize Nested Association Mapping (NAM) panel. Joint-linkage mapping of quantitative trait loci (QTL), fine mapping in near isogenic lines (NILs), genome-wide association studies (GWAS), and genomic best linear unbiased prediction (GBLUP) were performed. The heritability of maize height was estimated to be >90%. Mapping NAM family-nested QTL revealed the largest explained 2.1 ± 0.9% of height variation. The effects of two tropical alleles at this QTL were independently validated by fine mapping in NIL families. Several significant associations found by GWAS colocalized with established height loci, including brassinosteroid-deficient dwarf1, dwarf plant1, and semi-dwarf2. GBLUP explained >80% of height variation in the panels and outperformed bootstrap aggregation of family-nested QTL models in evaluations of prediction accuracy. These results revealed maize height was under strong genetic control and had a highly polygenic genetic architecture. They also showed that multiple models of genetic architecture differing in polygenicity and effect sizes can plausibly explain a population's variation in maize height, but they may vary in predictive efficacy.

Insights into the effects of long-term artificial selection on seed size in maize.

Grain produced from cereal crops is a primary source of human food and animal feed worldwide. To understand the genetic basis of seed-size variation, a grain yield component, we conducted a genome-wide scan to detect evidence of selection in the maize Krug Yellow Dent long-term divergent seed-size selection experiment. Previous studies have documented significant phenotypic divergence between the populations. Allele frequency estimates for ∼3 million single nucleotide polymorphisms (SNPs) in the base population and selected populations were estimated from pooled whole-genome resequencing of 48 individuals per population. Using FST values across sliding windows, 94 divergent regions with a median of six genes per region were identified. Additionally, 2729 SNPs that reached fixation in both selected populations with opposing fixed alleles were identified, many of which clustered in two regions of the genome. Copy-number variation was highly prevalent between the selected populations, with 532 total regions identified on the basis of read-depth variation and comparative genome hybridization. Regions important for seed weight in natural variation were identified in the maize nested association mapping population. However, the number of regions that overlapped with the long-term selection experiment did not exceed that expected by chance, possibly indicating unique sources of variation between the two populations. The results of this study provide insights into the genetic elements underlying seed-size variation in maize and could also have applications for other cereal crops.

Genetics and consequences of crop domestication.

Phenotypic variation has been manipulated by humans during crop domestication, which occurred primarily between 3000 and 10000 years ago in the various centers of origin around the world. The process of domestication has profound consequences on crops, where the domesticate has moderately reduced genetic diversity relative to the wild ancestor across the genome, and severely reduced diversity for genes targeted by domestication. The question that remains is whether reduction in genetic diversity has affected crop production today. A case study in maize ( Zea mays ) demonstrates the application of understanding relationships between genetic diversity and phenotypic diversity in the wild ancestor and the domesticate. As an outcrossing species, maize has tremendous genetic variation. The complementary combination of genome-wide association mapping (GWAS) approaches, large HapMap data sets, and germplasm resources is leading to important discoveries of the relationship between genetic diversity and phenotypic variation and the impact of domestication on trait variation.

Genetic variability of oxalate oxidase activity and elongation in water-stressed primary roots of diverse maize and rice lines.

A previous study of maize primary roots under water stress showed pronounced increases in oxalate oxidase activity and apoplastic hydrogen peroxide in the apical region of the growth zone where cell elongation is maintained. We examined whether increased oxalate oxidase activity in water-stressed roots is conserved across diverse lines of maize and rice. The maize lines exhibited varied patterns of activity, with some lines lacking activity in the apical region. Moreover, none of the rice lines showed activity in the apical region. Also, although the genotypic response of root elongation to water stress was variable in both maize and rice, this was not correlated with the pattern of oxalate oxidase activity. Implications of these findings for root growth regulation under water stress are discussed.

The genetic architecture of maize stalk strength.

Stalk strength is an important trait in maize (Zea mays L.). Strong stalks reduce lodging and maximize harvestable yield. Studies show rind penetrometer resistance (RPR), or the force required to pierce a stalk rind with a spike, is a valid approximation of strength. We measured RPR across 4,692 recombinant inbreds (RILs) comprising the maize nested association mapping (NAM) panel derived from crosses of diverse inbreds to the inbred, B73. An intermated B73×Mo17 family (IBM) of 196 RILs and a panel of 2,453 diverse inbreds from the North Central Regional Plant Introduction Station (NCRPIS) were also evaluated. We measured RPR in three environments. Family-nested QTL were identified by joint-linkage mapping in the NAM panel. We also performed a genome-wide association study (GWAS) and genomic best linear unbiased prediction (GBLUP) in each panel. Broad sense heritability computed on a line means basis was low for RPR. Only 8 of 26 families had a heritability above 0.20. The NCRPIS diversity panel had a heritability of 0.54. Across NAM and IBM families, 18 family-nested QTL and 141 significant GWAS associations were identified for RPR. Numerous weak associations were also found in the NCRPIS diversity panel. However, few were linked to loci involved in phenylpropanoid and cellulose synthesis or vegetative phase transition. Using an identity-by-state (IBS) relationship matrix estimated from 1.6 million single nucleotide polymorphisms (SNPs) and RPR measures from 20% of the NAM panel, genomic prediction by GBLUP explained 64±2% of variation in the remaining RILs. In the NCRPIS diversity panel, an IBS matrix estimated from 681,257 SNPs and RPR measures from 20% of the panel explained 33±3% of variation in the remaining inbreds. These results indicate the high genetic complexity of stalk strength and the potential for genomic prediction to hasten its improvement.

Comprehensive genotyping of the USA national maize inbred seed bank.

BACKGROUND: Genotyping by sequencing, a new low-cost, high-throughput sequencing technology was used to genotype 2,815 maize inbred accessions, preserved mostly at the National Plant Germplasm System in the USA. The collection includes inbred lines from breeding programs all over the world. RESULTS: The method produced 681,257 single-nucleotide polymorphism (SNP) markers distributed across the entire genome, with the ability to detect rare alleles at high confidence levels. More than half of the SNPs in the collection are rare. Although most rare alleles have been incorporated into public temperate breeding programs, only a modest amount of the available diversity is present in the commercial germplasm. Analysis of genetic distances shows population stratification, including a small number of large clusters centered on key lines. Nevertheless, an average fixation index of 0.06 indicates moderate differentiation between the three major maize subpopulations. Linkage disequilibrium (LD) decays very rapidly, but the extent of LD is highly dependent on the particular group of germplasm and region of the genome. The utility of these data for performing genome-wide association studies was tested with two simply inherited traits and one complex trait. We identified trait associations at SNPs very close to known candidate genes for kernel color, sweet corn, and flowering time; however, results suggest that more SNPs are needed to better explore the genetic architecture of complex traits. CONCLUSIONS: The genotypic information described here allows this publicly available panel to be exploited by researchers facing the challenges of sustainable agriculture through better knowledge of the nature of genetic diversity.

Heterosis is prevalent for multiple traits in diverse maize germplasm.

BACKGROUND: Heterosis describes the superior phenotypes observed in hybrids relative to their inbred parents. Maize is a model system for studying heterosis due to the high levels of yield heterosis and commercial use of hybrids. METHODS: The inbred lines from an association mapping panel were crossed to a common inbred line, B73, to generate nearly 300 hybrid genotypes. Heterosis was evaluated for seventeen phenotypic traits in multiple environments. The majority of hybrids exhibit better-parent heterosis in most of the hybrids measured. Correlations between the levels of heterosis for different traits were generally weak, suggesting that the genetic basis of heterosis is trait-dependent. CONCLUSIONS: The ability to predict heterosis levels using inbred phenotype or genetic distance between the parents varied for the different traits. For some traits it is possible to explain a significant proportion of the heterosis variation using linear modeling while other traits are more difficult to predict.

Maize association population: a high-resolution platform for quantitative trait locus dissection.

Crop improvement and the dissection of complex genetic traits require germplasm diversity. Although this necessary phenotypic variability exists in diverse maize, most research is conducted using a small subset of inbred lines. An association population of 302 lines is now available--a valuable research tool that captures a large proportion of the alleles in cultivated maize. Provided that appropriate statistical models correcting for population structure are included, this tool can be used in association analyses to provide high-resolution evaluation of multiple alleles. This study describes the population structure of the 302 lines, and investigates the relationship between population structure and various measures of phenotypic and breeding value. On average, our estimates of population structure account for 9.3% of phenotypic variation, roughly equivalent to a major quantitative trait locus (QTL), with a high of 35%. Inclusion of population structure in association models is critical to meaningful analyses. This new association population has the potential to identify QTL with small effects, which will aid in dissecting complex traits and in planning future projects to exploit the rich diversity present in maize.