Reveal the kernel dehydration mechanisms in maize based on proteomic and metabolomic analysis.

BACKGROUND: Kernel dehydration is an important factor for the mechanized harvest in maize. Kernel moisture content (KMC) and kernel dehydration rate (KDR) are important indicators for kernel dehydration. Although quantitative trait loci and genes related to KMC have been identified, where most of them only focus on the KMC at harvest, these are still far from sufficient to explain all genetic variations, and the relevant regulatory mechanisms are still unclear. In this study, we tried to reveal the key proteins and metabolites related to kernel dehydration in proteome and metabolome levels. Moreover, we preliminarily explored the relevant metabolic pathways that affect kernel dehydration combined proteome and metabolome. These results could accelerate the development of further mechanized maize technologies. RESULTS: In this study, three maize inbred lines (KB182, KB207, and KB020) with different KMC and KDR were subjected to proteomic analysis 35, 42, and 49 days after pollination (DAP). In total, 8,358 proteins were quantified, and 2,779 of them were differentially expressed proteins in different inbred lines or at different stages. By comparative analysis, K-means cluster, and weighted gene co-expression network analysis based on the proteome data, some important proteins were identified, which are involved in carbohydrate metabolism, stress and defense response, lipid metabolism, and seed development. Through metabolomics analysis of KB182 and KB020 kernels at 42 DAP, 18 significantly different metabolites, including glucose, fructose, proline, and glycerol, were identified. CONCLUSIONS: In sum, we inferred that kernel dehydration could be regulated through carbohydrate metabolism, antioxidant systems, and late embryogenesis abundant protein and heat shock protein expression, all of which were considered as important regulatory factors during kernel dehydration process. These results shed light on kernel dehydration and provide new insights into developing cultivars with low moisture content.

Linkage and association mapping in multi-parental populations reveal the genetic basis of carotenoid variation in maize kernels.

Carotenoids are indispensable to plants and critical components of the human diet. The carotenoid metabolic pathway is conserved across plant species, but our understanding of the genetic basis of carotenoid variation remains limited for the seeds of most cereal crops. To address this issue, we systematically performed linkage and association mapping for eight carotenoid traits using six recombinant inbred line (RIL) populations. Single linkage mapping (SLM) and joint linkage mapping (JLM) identified 77 unique additive QTLs and 104 pairs of epistatic QTLs. Among these QTLs, we identified 22 overlapping hotspots of additive and epistatic loci, highlighting the important contributions of some QTLs to carotenoid levels through additive or epistatic mechanisms. A genome-wide association study based on all RILs detected 244 candidate genes significantly associated with carotenoid traits, 23 of which were annotated as carotenoid pathway genes. Effect comparisons suggested that a small number of loci linked to pathway genes have substantial effects on carotenoid variation in our tested populations, but many loci not associated with pathway genes also make important contributions to carotenoid variation. We identified ZmPTOX as the causal gene for a QTL hotspot (Q10/JLM10/GWAS019); this gene encodes a putative plastid terminal oxidase that produces plastoquinone-9 used by two enzymes in the carotenoid pathway. Natural variants in the promoter and second exon of ZmPTOX were found to alter carotenoid levels. This comprehensive assessment of the genetic mechanisms underlying carotenoid variation establishes a foundation for rewiring carotenoid metabolism and accumulation for efficient carotenoid biofortification.

Genome-Wide Association Study for Maize Hybrid Performance in a Typical Breeder Population.

Maize is one of the major crops that has demonstrated success in the utilization of heterosis. Developing high-yield hybrids is a crucial part of plant breeding to secure global food demand. In this study, we conducted a genome-wide association study (GWAS) for 10 agronomic traits using a typical breeder population comprised 442 single-cross hybrids by evaluating additive, dominance, and epistatic effects. A total of 49 significant single nucleotide polymorphisms (SNPs) and 69 significant pairs of epistasis were identified, explaining 26.2% to 64.3% of the phenotypic variation across the 10 traits. The enrichment of favorable genotypes is significantly correlated to the corresponding phenotype. In the confident region of the associated site, 532 protein-coding genes were discovered. Among these genes, the Zm00001d044211 candidate gene was found to negatively regulate starch synthesis and potentially impact yield. This typical breeding population provided a valuable resource for dissecting the genetic architecture of yield-related traits. We proposed a novel mating strategy to increase the GWAS efficiency without utilizing more resources. Finally, we analyzed the enrichment of favorable alleles in the Shaan A and Shaan B groups, as well as in each inbred line. Our breeding practice led to consistent results. Not only does this study demonstrate the feasibility of GWAS in F1 hybrid populations, it also provides a valuable basis for further molecular biology and breeding research.

Genetic architecture of ear traits based on association mapping and co-expression networks in maize inbred lines and hybrids.

Ear traits are key contributors to grain yield in maize; therefore, exploring their genetic basis facilitates the improvement of grain yield. However, the underlying molecular mechanisms of ear traits remain obscure in both inbred lines and hybrids. Here, two association panels, respectively, comprising 203 inbred lines (IP) and 246 F1 hybrids (HP) were employed to identify candidate genes for six ear traits. The IP showed higher phenotypic variation and lower phenotypic mean than the HP for all traits, except ear tip-barrenness length. By conducting a genome-wide association study (GWAS) across multiple environments, 101 and 228 significant single-nucleotide polymorphisms (SNPs) associated with six ear traits were identified in the IP and HP, respectively. Of these significant SNPs identified in the HP, most showed complete-incomplete dominance and over-dominance effects for each ear trait. Combining a gene co-expression network with GWAS results, 186 and 440 candidate genes were predicted in the IP and HP, respectively, including known ear development genes ids1 and sid1. Of these, nine candidate genes were detected in both populations and expressed in maize ear and spikelet tissues. Furthermore, two key shared genes (GRMZM2G143330 and GRMZM2G171139) in both populations were found to be significantly associated with ear traits in the maize Goodman diversity panel with high-density variations. These findings advance our knowledge of the genetic architecture of ear traits between inbred lines and hybrids and provide a valuable resource for the genetic improvement of ear traits in maize. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-023-01426-9.

Genome-wide association analysis of kernel nutritional quality in two natural maize populations.

As one of the three staple crops, nutritional traits in maize are important for human and animal nutrition. Grain quality-related traits are closely related to grain commercial value. Understanding the genetic basis of quality-related traits in maize would be helpful for breeding high-quality maize varieties. In this study, two association panels (AM122 and AM180) were subjected to genome-wide association analysis of grain quality-related traits, including protein content, oil content, starch content, and fiber content. In total, 98 SNPs (P < 1 × 10-4) were identified to be significantly associated with these four grain quality-related traits. By integrating two sets of public transcriptome data, 31 genes located in 200 kb regions flanking the associated SNP showed high expression during kernel development and were differentially expressed in two maize inbred lines, KA225 and KB035, with significantly different quality. These genes might regulate maize grain quality by participating in plant hormone processes, autophagy processes, and others. All these results could provide important reference information for breeding high­quality maize varieties. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-023-01360-w.

Genome-wide association analysis for grain moisture content and dehydration rate on maize hybrids.

For mechanized maize production, a low grain water content (GWC) at harvest is necessary. However, as a complex quantitative trait, understand the genetic mechanism of GWC remains a large gap, especially in hybrids. In this study, a hybrid population through two environments including 442 F1 was used for genome-wide association analysis of GWC and the grain dehydration rate (GDR), using the area under the dry down curve (AUDDC) as the index. Then, we identified 19 and 17 associated SNPs for GWC and AUDDC, including 10 co-localized SNPs, along with 64 and 77 pairs of epistatic SNPs for GWC and AUDDC, respectively. These loci could explain 11.39-68.2% of the total phenotypic variation for GWC and 41.07-67.02% for AUDDC at different stages, whose major effect was the additive and epistatic effect. By exploring the candidate genes around the significant sites, a total of 398 and 457 possible protein-coding genes were screened, including autophagy pathway and auxin regulation-related genes, and five inbred lines with the potential to reduce GWC in the combined F1 hybrid were identified. Our research not only provides a certain reference for the genetic mechanism analysis of GWC in hybrids but also provides an added reference for breeding low-GWC materials. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-022-01349-x.

A cost-effective tsCUT&Tag method for profiling transcription factor binding landscape.

Knowledge of the transcription factor binding landscape (TFBL) is necessary to analyze gene regulatory networks for important agronomic traits. However, a low-cost and high-throughput in vivo chromatin profiling method is still lacking in plants. Here, we developed a transient and simplified cleavage under targets and tagmentation (tsCUT&Tag) that combines transient expression of transcription factor proteins in protoplasts with a simplified CUT&Tag without nucleus extraction. Our tsCUT&Tag method provided higher data quality and signal resolution with lower sequencing depth compared with traditional ChIP-seq. Furthermore, we developed a strategy combining tsCUT&Tag with machine learning, which has great potential for profiling the TFBL across plant development.

New insights into the response of maize to fluctuations in the light environment.

Light is the most important environmental cue signaling the transition from skotomorphogenesis to photomorphogenesis, thus affecting plant development and metabolic activity. How the light response mechanisms of maize seedlings respond to fluctuations in the light environment has not been well characterized to date. In this study, we built a gene coexpression network from a dynamic transcriptomic map of maize seedlings exposed to different light environments. Coexpression analysis identified ten modules and multiple genes that closely correlate with photosynthesis and characterized hub genes associated with regulatory networks, duplication events, domestication and improvement. In addition, we identified that 38% of hub genes underwent duplication events, 74% of which are related to photosynthesis. Moreover, we captured the dynamic expression atlas of gene sets involved in the chloroplast photosynthetic apparatus and photosynthetic carbon assimilation in different light environments, which should help to elucidate the key mechanisms and regulatory networks that underlie photosynthesis in maize. Insights from this study provide a valuable resource to better understand the genetic mechanisms of the response to fluctuations in the light environment in maize.

Genome-wide evolutionary characterization and expression analysis of SIAMESE-RELATED family genes in maize.

BACKGROUND: The SIAMESE (SIM) locus is a cell-cycle kinase inhibitor (CKI) gene that has to date been identified only in plants; it encodes a protein that promotes transformation from mitosis to endoreplication. Members of the SIAMESE-RELATED (SMR) family have similar functions, and some are related to cell-cycle responses and abiotic stresses. However, the functions of SMRs are poorly understood in maize (Zea mays L.). RESULTS: In the present study, 12 putative SMRs were identified throughout the entire genome of maize, and these were clustered into six groups together with the SMRs from seven other plant species. Members of the ZmSMR family were divided into four groups according to their protein sequences. Various cis-acting elements in the upstream sequences of ZmSMRs responded to abiotic stresses. Expression analyses revealed that all ZmSMRs were upregulated at 5, 20, 25, and 35 days after pollination. In addition, we found that ZmSMR9/11/12 may have regulated the initiation of endoreplication in endosperm central cells. Additionally, ZmSMR2/10 may have been primarily responsible for the endoreplication regulation of outer endosperm or aleurone cells. The relatively high expression levels of almost all ZmSMRs in the ears and tassels also implied that these genes may function in seed development. The effects of treatments with ABA, heat, cold, salt, and drought on maize seedlings and expression of ZmSMR genes suggested that ZmSMRs were strongly associated with response to abiotic stresses. CONCLUSION: The present study is the first to conduct a genome-wide analysis of members of the ZmSMR family by investigating their locations in chromosomes, identifying regulatory elements in their promoter regions, and examining motifs in their protein sequences. Expression analysis of different endosperm developmental periods, tissues, abiotic stresses, and hormonal treatments suggests that ZmSMR genes may function in endoreplication and regulate the development of reproductive organs. These results may provide valuable information for future studies of the functions of the SMR family in maize.

The distribution pattern of endopolyploidy in maize.

KEY MESSAGE: We discovered that endopolyploidization is common in various organs and tissues of maize at different development stages. Endopolyploidy is not specific in maize germplasm populations. Endopolyploidy is caused by DNA endoreplication, a special type of mitosis with normal DNA synthesis and a lack of cell division; it is a common phenomenon and plays an important role in plant development. To systematically study the distribution pattern of endopolyploidy in maize, flow cytometry was used to determine the ploidy by measuring the cycle (C) value in various organs at different developmental stages, in embryos and endosperm during grain development, in roots under stress conditions, and in the roots of 119 inbred lines from two heterotic groups, Shaan A and Shaan B. Endopolyploidy was observed in most organs at various developmental stages except in expanded leaves and filaments. The endosperm showed the highest C value among all organs. During tissue development, the ploidy increased in all organs except the leaves. In addition, the endopolyploidization of the roots was significantly affected by drought stress. Multiple comparisons of the C values of seven subgroups revealed that the distribution of endopolyploidization was not correlated with the population structure. A correlation analysis at the seedling stage showed a positive relationship between the C value and both the length of the whole plant and the length of main root. A genome-wide association study (GWAS) identified a total of 9 significant SNPs associated with endopolyploidy (C value) in maize, and 8 candidate genes that participate in cell cycle regulation and DNA replication were uncovered in 119 maize inbred lines.

Beyond pathways: genetic dissection of tocopherol content in maize kernels by combining linkage and association analyses.

Although tocopherols play an important role in plants and animals, the genetic architecture of tocopherol content in maize kernels has remained largely unknown. In this study, linkage and association analyses were conducted to examine the genetic architecture of tocopherol content in maize kernels. Forty-one unique quantitative trait loci (QTLs) were identified by linkage mapping in six populations of recombinant inbred lines (RILs). In addition, 32 significant loci were detected via genome-wide association study (GWAS), 18 of which colocalized with the QTLs identified by linkage mapping. Fine mapping of a major QTL validated the accuracy of GWAS and QTL mapping results and suggested a role for nontocopherol pathway genes in the modulation of natural tocopherol variation. We provided genome-wide evidence that genes involved in fatty acid metabolism, chlorophyll metabolism and chloroplast function may affect natural variation in tocopherols. These findings were confirmed through mutant analysis of a particular gene from the fatty acid pathway. In addition, the favourable alleles for many of the significant SNPs/QTLs represented rare alleles in natural populations. Together, our results revealed many novel genes that are potentially involved in the variation of tocopherol content in maize kernels. Pyramiding of the favourable alleles of the newly elucidated genes and the well-known tocopherol pathway genes would greatly improve tocopherol content in maize.

Bivariate flow cytometric analysis and sorting of different types of maize starch grains.

Particle-size distribution, granular structure, and composition significantly affect the physicochemical properties, rheological properties, and nutritional function of starch. Flow cytometry and flow sorting are widely considered convenient and efficient ways of classifying and separating natural biological particles or other substances into subpopulations, respectively, based on the differential response of each component to stimulation by a light beam; the results allow for the correlation analysis of parameters. In this study, different types of starches isolated from waxy maize, sweet maize, high-amylose maize, pop maize, and normal maize were initially classified into various subgroups by flow cytometer and then collected through flow sorting to observe their morphology and particle-size distribution. The results showed that a 0.25% Gelzan solution served as an optimal reagent for keeping individual starch particles homogeneously dispersed in suspension for a relatively long time. The bivariate flow cytometric population distributions indicated that the starches of normal maize, sweet maize, and pop maize were divided into two subgroups, whereas high-amylose maize starch had only one subgroup. Waxy maize starch, conversely, showed three subpopulations. The subgroups sorted by flow cytometer were determined and verified in terms of morphology and granule size by scanning electron microscopy and laser particle distribution analyzer. Results showed that flow cytometry can be regarded as a novel method for classifying and sorting starch granules. © 2017 International Society for Advancement of Cytometry.

Genetic characterization of inbred lines from Shaan A and B groups for identifying loci associated with maize grain yield.

BACKGROUND: Increasing grain yield is a primary objective of maize breeding. Dissecting the genetic architecture of grain yield furthers genetic improvements to increase yield. Presented here is an association panel composed of 126 maize inbreds (AM126), which were genotyped by the genotyping-by-sequencing (tGBS) method. We performed genetic characterization and association analysis related to grain yield in the association panel. RESULTS: In total, 46,046 SNPs with a minor allele frequency (MAF) ≥0.01 were used to assess genetic diversity and kinship in AM126. The results showed that the average MAF and polymorphism information content (PIC) were 0.164 and 0.198, respectively. The Shaan B group, with 11,284 unique SNPs, exhibited greater genetic diversity than did the Shaan A group, with 2644 SNPs. The 61.82% kinship coefficient in AM126 was equal to 0, and only 0.15% of that percentage was greater than 0.7. A total of 31,983 SNPs with MAF ≥0.05 were used to characterize population structure, LD decay and association mapping. Population structure analysis suggested that AM126 can be divided into 6 subgroups, which is consistent with breeding experience and pedigree information. The LD decay distance in AM126 was 150 kb. A total of 51 significant SNPs associated with grain yield were identified at P < 1 × 10- 3 across two environments (Yangling and Yulin). Among those SNPs, two loci displayed overlapping regions in the two environments. Finally, 30 candidate genes were found to be associated with grain yield. CONCLUSIONS: These results contribute to the genetic characterization of this breeding population, which serves as a reference for hybrid breeding and population improvement, and demonstrate the genetic architecture of maize grain yield, potentially facilitating genetic improvement.

Genome-wide recombination dynamics are associated with phenotypic variation in maize.

Meiotic recombination is a major driver of genetic diversity, species evolution, and agricultural improvement. Thus, an understanding of the genetic recombination landscape across the maize (Zea mays) genome will provide insight and tools for further study of maize evolution and improvement. Here, we used c. 50 000 single nucleotide polymorphisms to precisely map recombination events in 12 artificial maize segregating populations. We observed substantial variation in the recombination frequency and distribution along the ten maize chromosomes among the 12 populations and identified 143 recombination hot regions. Recombination breakpoints were partitioned into intragenic and intergenic events. Interestingly, an increase in the number of genes containing recombination events was accompanied by a decrease in the number of recombination events per gene. This kept the overall number of intragenic recombination events nearly invariable in a given population, suggesting that the recombination variation observed among populations was largely attributed to intergenic recombination. However, significant associations between intragenic recombination events and variation in gene expression and agronomic traits were observed, suggesting potential roles for intragenic recombination in plant phenotypic diversity. Our results provide a comprehensive view of the maize recombination landscape, and show an association between recombination, gene expression and phenotypic variation, which may enhance crop genetic improvement.

Deciphering the Genetic Basis of Kernel Composition in a Maize Association Panel.

The primary objective in contemporary maize breeding is to pursue high quality alongside high yield. Deciphering the genetic basis of natural variation in starch, protein, oil, and fiber contents is essential for manipulating kernel composition, thereby enhancing the kernel quality and meeting growing demands. Here, we identified 12 to 88 statistically significant loci associated with kernel composition traits through a genome-wide association study (GWAS) using a panel of 212 diverse inbred lines. A regional association study pinpointed numerous causal candidate genes at these loci. Coexpression and protein-protein interaction network analyses of candidate genes revealed several causal genes directly or indirectly involved in the metabolic processes related to kernel composition traits. Subsequent mutant experiment revealed that nonsense mutations in ZmTIFY12 affect starch, protein, and fiber content, whereas nonsense mutations in ZmTT12 affect starch, protein, and oil content. These findings provide valuable guidance for improving kernel quality in maize breeding efforts.

Maize Autophagy-Related Protein ZmATG3 Confers Tolerance to Multiple Abiotic Stresses.

Abiotic stresses pose a major increasing problem for the cultivation of maize. Autophagy plays a vital role in recycling and re-utilizing nutrients and adapting to stress. However, the role of autophagy in the response to abiotic stress in maize has not yet been investigated. Here, ZmATG3, which is essential for ATG8-PE conjugation, was isolated from the maize inbred line B73. The ATG3 sequence was conserved, including the C-terminal domains with HPC and FLKF motifs and the catalytic domain in different species. The promoter of the ZmATG3 gene contained a number of elements involved in responses to environmental stresses or hormones. Heterologous expression of ZmATG3 in yeast promoted the growth of strain under salt, mannitol, and low-nitrogen stress. The expression of ZmATG3 could be altered by various types of abiotic stress (200 mM NaCl, 200 mM mannitol, low N) and exogenous hormones (500 µM ABA). GUS staining analysis of ZmATG3-GUS transgenic Arabidopsis revealed that GUS gene activity increased after abiotic treatment. ZmATG3-overexpressing Arabidopsis plants had higher osmotic and salinity stress tolerance than wild-type plants. Overexpression of ZmATG3 up-regulated the expression of other AtATGs (AtATG3, AtATG5, and AtATG8b) under NaCl, mannitol and LN stress. These findings demonstrate that overexpression of ZmATG3 can improve tolerance to multiple abiotic stresses.

QTL mapping for plant height and ear height using bi-parental immortalized heterozygous populations in maize.

Introduction: Plant height (PH) and ear height (EH) are key plant architectural traits in maize, which will affect the photosynthetic efficiency, high plant density tolerance, suitability for mechanical harvesting. Methods: QTL mapping were conducted for PH and EH using a recombinant inbred line (RIL) population and two corresponding immortalized backcross (IB) populations obtained from crosses between the RIL population and the two parental lines. Results: A total of 17 and 15 QTL were detected in the RIL and IB populations, respectively. Two QTL, qPH1-1 (qEH1-1) and qPH1-2 (qEH1-4) in the RIL, were simultaneously identified for PH and EH. Combing reported genome-wide association and cloned PH-related genes, co-expression network analyses were constructed, then five candidate genes with high confidence in major QTL were identified including Zm00001d011117 and Zm00001d011108, whose homologs have been confirmed to play a role in determining PH in maize and soybean. Discussion: QTL mapping used a immortalized backcross population is a new strategy. These identified genes in this study can provide new insights for improving the plant architecture in maize.

A dynamic regulome of shoot-apical-meristem-related homeobox transcription factors modulates plant architecture in maize.

BACKGROUND: The shoot apical meristem (SAM), from which all above-ground tissues of plants are derived, is critical to plant morphology and development. In maize (Zea mays), loss-of-function mutant studies have identified several SAM-related genes, most encoding homeobox transcription factors (TFs), located upstream of hierarchical networks of hundreds of genes. RESULTS: Here, we collect 46 transcriptome and 16 translatome datasets across 62 different tissues or stages from the maize inbred line B73. We construct a dynamic regulome for 27 members of three SAM-related homeobox subfamilies (KNOX, WOX, and ZF-HD) through machine-learning models for the detection of TF targets across different tissues and stages by combining tsCUT&Tag, ATAC-seq, and expression profiling. This dynamic regulome demonstrates the distinct binding specificity and co-factors for these homeobox subfamilies, indicative of functional divergence between and within them. Furthermore, we assemble a SAM dynamic regulome, illustrating potential functional mechanisms associated with plant architecture. Lastly, we generate a wox13a mutant that provides evidence that WOX13A directly regulates Gn1 expression to modulate plant height, validating the regulome of SAM-related homeobox genes. CONCLUSIONS: The SAM-related homeobox transcription-factor regulome presents an unprecedented opportunity to dissect the molecular mechanisms governing SAM maintenance and development, thereby advancing our understanding of maize growth and shoot architecture.

The light and hypoxia induced gene ZmPORB1 determines tocopherol content in the maize kernel.

Tocopherol is an important lipid-soluble antioxidant beneficial for both human health and plant growth. Here, we fine mapped a major QTL-qVE1 affecting γ-tocopherol content in maize kernel, positionally cloned and confirmed the underlying gene ZmPORB1 (por1), as a protochlorophyllide oxidoreductase. A 13.7 kb insertion reduced the tocopherol and chlorophyll content, and the photosynthetic activity by repressing ZmPORB1 expression in embryos of NIL-K22, but did not affect the levels of the tocopherol precursors HGA (homogentisic acid) and PMP (phytyl monophosphate). Furthermore, ZmPORB1 is inducible by low oxygen and light, thereby involved in the hypoxia response in developing embryos. Concurrent with natural hypoxia in embryos, the redox state has been changed with NO increasing and H2O2 decreasing, which lowered γ-tocopherol content via scavenging reactive nitrogen species. In conclusion, we proposed that the lower light-harvesting chlorophyll content weakened embryo photosynthesis, leading to fewer oxygen supplies and consequently diverse hypoxic responses including an elevated γ-tocopherol consumption. Our findings shed light on the mechanism for fine-tuning endogenous oxygen concentration in the maize embryo through a novel feedback pathway involving the light and low oxygen regulation of ZmPORB1 expression and chlorophyll content.

Natural variation in the sequence of PSY1 and frequency of favorable polymorphisms among tropical and temperate maize germplasm.

Provitamin A (Pro-VA) is necessary for human vision and immune system health, especially in growing children. The first committed step in the maize carotenoid biosynthesis pathway is catalyzed by phytoene synthase 1 (encoded by PSY1) which controls the flux of substrates into the pathway. The flow of these substrates could be directed into production of the ß-branch carotenoids (the step controlled largely by the lycopene epsilon cyclase gene), but terminated after the production of ß-carotene, rather than allowing it to be converted into the next metabolite (the step controlled largely by the ß-carotenoid hydroxylase gene). In this study, PSY1 was subjected to association mapping in two diverse maize populations, quantitative trait loci (QTL) mapping in one segregating population, and expression analysis of lines polymorphic for sites within PSY1. The results indicated that a 378-bp InDel upstream of the transcription start site and a SNP in the fifth exon resulting in a Thr to Asn substitution, explaining 7 and 8 % of the total carotenoid variation, respectively, may be functional sites associated with total carotenoid levels in maize grain. Analysis of the evolution of PSY1 strongly suggests that there was positive selection for these polymorphic sites after the divergence of yellow maize from white maize.