Genetics, not environment, contributed to winter wheat yield gains in Ontario, Canada.

KEY MESSAGE: Changes in entries' market classes and genetic improvements within classes-not environmental changes-enhanced yields over thirty-one years of wheat trials. Correlations between yields and ancestries drove genomic prediction accuracies. Increasing crop yields is important for enhancing farmers' livelihoods, meeting market demands, and reducing the environmental impact of agriculture. We analyzed the yield trends of Ontario winter wheat variety trials between 1988 and 2018. Over this period, wheat yields steadily increased by 38 kg ha-1 yr-1, or 0.68% yr-1 relative to the mean. While fungicide treatment of trials contributed a one-time 670 kg ha-1 yield increase, yields were otherwise unaffected by long-term changes in agronomic practice, climate, or other non-genetic factors. Genetic improvement entirely accounted for yield improvement. Market class changes over the 31 year span accounted for some yield improvement. More importantly, genetic improvement occurred within each market class. Entry yield estimates calculated from genomic prediction models strongly correlated with field estimated yields with a mean r of 0.68. Genomic prediction accuracies were high because yields differed across genetically distinct subpopulations. Despite environmental changes, genetic improvement will likely increase Ontario winter wheat yields into the future.

Zea mays RNA-seq estimated transcript abundances are strongly affected by read mapping bias.

BACKGROUND: Genetic variation for gene expression is a source of phenotypic variation for natural and agricultural species. The common approach to map and to quantify gene expression from genetically distinct individuals is to assign their RNA-seq reads to a single reference genome. However, RNA-seq reads from alleles dissimilar to this reference genome may fail to map correctly, causing transcript levels to be underestimated. Presently, the extent of this mapping problem is not clear, particularly in highly diverse species. We investigated if mapping bias occurred and if chromosomal features associated with mapping bias. Zea mays presents a model species to assess these questions, given it has genotypically distinct and well-studied genetic lines. RESULTS: In Zea mays, the inbred B73 genome is the standard reference genome and template for RNA-seq read assignments. In the absence of mapping bias, B73 and a second inbred line, Mo17, would each have an approximately equal number of regulatory alleles that increase gene expression. Remarkably, Mo17 had 2-4 times fewer such positively acting alleles than did B73 when RNA-seq reads were aligned to the B73 reference genome. Reciprocally, over one-half of the B73 alleles that increased gene expression were not detected when reads were aligned to the Mo17 genome template. Genes at dissimilar chromosomal ends were strongly affected by mapping bias, and genes at more similar pericentromeric regions were less affected. Biased transcript estimates were higher in untranslated regions and lower in splice junctions. Bias occurred across software and alignment parameters. CONCLUSIONS: Mapping bias very strongly affects gene transcript abundance estimates in maize, and bias varies across chromosomal features. Individual genome or transcriptome templates are likely necessary for accurate transcript estimation across genetically variable individuals in maize and other species.

Patterns of stability and change in the maize genome: a case study of small RNA transcriptomes in two recombinant inbred lines and their progenitors.

Small RNAs (sRNAs) are epigenetic regulators of eukaryotic genes and transposable elements (TEs). Diverse sRNA expression patterns exist within a species, but how this diversity arises is not well understood. To provide a window into the dynamics of maize sRNA patterning, sRNA and mRNA transcriptomes were examined in two related Zea mays recombinant inbred lines (RILs) and their inbred parents. Analysis of these RILs revealed that most clusters of sRNA expression retained the parental sRNA expression level. However, expression states that differ from the parental allele were also observed, predominantly reflecting decreases in sRNA expression. When RIL sRNA expression differed from the parental allele, the new state was frequently similar between the two RILs, and similar to the expression state found at the allele in the other parent. Novel sRNA expression patterns, distinct from those of either parent, were rare. Additionally, examination of sRNA expression over TEs revealed one TE family, Gyma, which showed consistent enrichment for RIL sRNA expression differences compared to those found in parental alleles. These findings provide insights into how sRNA silencing might evolve over generations and suggest that further investigation into the molecular nature of sRNA trans regulators is warranted.

The neonicotinoid insecticide thiamethoxam enhances expression of stress-response genes in Zea mays in an environmentally specific pattern.

Recent studies indicate that thiamethoxam (TMX), a neonicotinoid insecticide, can affect plant responses to environmental stressors, such as neighboring weeds. The molecular mechanisms behind both stable and environmentally specific responses to TMX likely involve genes related to defense and stress responses. We investigated the effect of a TMX seed treatment on global gene expression in maize coleoptiles both under normal conditions and under low ratio red to far-red (R:FR) light stress induced by the presence of neighboring plants. The neighboring plant treatment upregulated genes involved in biotic and abiotic stress responses and affected specific photosynthesis and cell-growth related genes. Low R:FR light may enhance maize resistance to herbivores and pathogens. TMX appears to compromise resistance. The TMX treatment stably repressed many genes that encode proteins involved in biotic stress responses, as well as cell-growth genes. Notably, TMX effects on many genes' expression were conditional on the environment. In response to low R:FR, plants treated with TMX engage genes in the JA pathway, as well as other stress-related response pathways. Neighboring weeds may condition TMX-treated plants to become more stress tolerant.

Introgressed DNA within a Zea mays near-isogenic line displays lower levels of 24nt sRNA expression than the homologous region from the recurrent parent.

Near-isogenic lines (NILs) are classical genetic tools used to dissect the actions of an allele when placed in a uniform genetic background. Although the goal of NIL creation is to examine the effects of a single allele in isolation, DNA linked to the allele is invariably retained and can confound any allele-specific effects. In addition to genetic variation, highly polymorphic species such as Zea mays will contain introgressed polymorphisms encompassing transposable elements (TEs) and the cis-acting small RNA (sRNA) that represses them. Through transcriptomics, we described the differences in sRNA and TE transcriptional expression between a W22-derived introgression and its homologous B73 region. As anticipated, many differences in sRNA expression were observed. Unexpectedly, however, 24nt sRNA expression over the introgressed region was low overall compared to both the homologous B73 region and the rest of the genome. Across the introgression, low sRNA expression was accompanied by increased TE transcription. Possible explanations for the observed trends in sRNA and TE expression across the introgression region are discussed. These findings support the notion that any introgressed allele is in an epigenetic environment distinct from that found at the allele from the recurrent parent. Additionally, these results suggest that further study of sRNA expression levels during the introgression process is warranted.

Genomic regions associated with important seed quality traits in food-grade soybeans.

BACKGROUND: The production of soy-based food products requires specific physical and chemical characteristics of the soybean seed. Identification of quantitative trait loci (QTL) associated with value-added traits, such as seed weight, seed protein and sucrose concentration, could accelerate the development of competitive high-protein soybean cultivars for the food-grade market through marker-assisted selection (MAS). The objectives of this study were to identify and validate QTL associated with these value-added traits in two high-protein recombinant inbred line (RIL) populations. RESULTS: The RIL populations were derived from the high-protein cultivar 'AC X790P' (49% protein, dry weight basis), and two high-yielding commercial cultivars, 'S18-R6' (41% protein) and 'S23-T5' (42% protein). Fourteen large-effect QTL (R2 > 10%) were identified associated with seed protein concentration. Of these QTL, seven QTL were detected in both populations, and eight of them were co-localized with QTL associated with either seed sucrose concentration or seed weight. None of the protein-related QTL was found to be associated with seed yield in either population. Sixteen candidate genes with putative roles in protein metabolism were identified within seven of these protein-related regions: qPro_Gm02-3, qPro_Gm04-4, qPro_Gm06-1, qPro_Gm06-3, qPro_Gm06-6, qPro_Gm13-4 and qPro-Gm15-3. CONCLUSION: The use of RIL populations derived from high-protein parents created an opportunity to identify four novel QTL that may have been masked by large-effect QTL segregating in populations developed from diverse parental cultivars. In total, we have identified nine protein QTL that were detected either in both populations in the current study or reported in other studies. These QTL may be useful in the curated selection of new soybean cultivars for optimized soy-based food products.

Shared and genetically distinct Zea mays transcriptome responses to ongoing and past low temperature exposure.

BACKGROUND: Cold temperatures and their alleviation affect many plant traits including the abundance of protein coding gene transcripts. Transcript level changes that occur in response to cold temperatures and their alleviation are shared or vary across genotypes. In this study we identify individual transcripts and groups of functionally related transcripts that consistently respond to cold and its alleviation. Genes that respond differently to temperature changes across genotypes may have limited functional importance. We investigate if these genes share functions, and if their genotype-specific gene expression levels change in magnitude or rank across temperatures. RESULTS: We estimate transcript abundances from over 22,000 genes in two unrelated Zea mays inbred lines during and after cold temperature exposure. Genotype and temperature contribute to many genes' abundances. Past cold exposure affects many fewer genes. Genes up-regulated in cold encode many cytokinin glucoside biosynthesis enzymes, transcription factors, signalling molecules, and proteins involved in diverse environmental responses. After cold exposure, protease inhibitors and cuticular wax genes are newly up-regulated, and environmentally responsive genes continue to be up-regulated. Genes down-regulated in response to cold include many photosynthesis, translation, and DNA replication associated genes. After cold exposure, DNA replication and translation genes are still preferentially downregulated. Lignin and suberin biosynthesis are newly down-regulated. DNA replication, reactive oxygen species response, and anthocyanin biosynthesis genes have strong, genotype-specific temperature responses. The ranks of genotypes' transcript abundances often change across temperatures. CONCLUSIONS: We report a large, core transcriptome response to cold and the alleviation of cold. In cold, many of the core suite of genes are up or downregulated to control plant growth and photosynthesis and limit cellular damage. In recovery, core responses are in part to prepare for future stress. Functionally related genes are consistently and greatly up-regulated in a single genotype in response to cold or its alleviation, suggesting positive selection has driven genotype-specific temperature responses in maize.

Distinct gene networks modulate floral induction of autonomous maize and photoperiod-dependent teosinte.

Temperate maize was domesticated from its tropical ancestor, teosinte. Whereas temperate maize is an autonomous day-neutral plant, teosinte is an obligate short-day plant that requires uninterrupted long nights to induce flowering. Leaf-derived florigenic signals trigger reproductive growth in both teosinte and temperate maize. To study the genetic mechanisms underlying floral inductive pathways in maize and teosinte, mRNA and small RNA genome-wide expression analyses were conducted on leaf tissue from plants that were induced or not induced to flower. Transcriptome profiles reveal common differentially expressed genes during floral induction, but a comparison of candidate flowering time genes indicates that photoperiod and autonomous pathways act independently. Expression differences in teosinte are consistent with the current paradigm for photoperiod-induced flowering, where changes in circadian clock output trigger florigen production. Conversely, differentially expressed genes in temperate maize link carbon partitioning and flowering, but also show altered expression of circadian clock genes that are distinct from those altered upon photoperiodic induction in teosinte. Altered miRNA399 levels in both teosinte and maize suggest a novel common connection between flowering and phosphorus perception. These findings provide insights into the molecular mechanisms underlying a strengthened autonomous pathway that enabled maize growth throughout temperate regions.

Proanthocyanidin accumulation and transcriptional responses in the seed coat of cranberry beans (Phaseolus vulgaris L.) with different susceptibility to postharvest darkening.

BACKGROUND: Edible dry beans (Phaseolus vulgaris L.) that darken during postharvest storage are graded lower and are less marketable than their non-darkened counterparts. Seed coat darkening in susceptible genotypes is dependent upon the availability of proanthocyanidins, and their subsequent oxidation to reactive quinones. Mature cranberry beans lacking this postharvest darkening trait tend to be proanthocyanidin-deficient, although the underlying molecular and biochemical determinants for this metabolic phenomenon are unknown. RESULTS: Seed coat proanthocyanidin levels increased with plant maturation in a darkening-susceptible cranberry bean recombinant inbred line (RIL), whereas these metabolites were absent in seeds of the non-darkening RIL plants. RNA sequencing (RNA-seq) analysis was used to monitor changes in the seed coat transcriptome as a function of bean development, where transcript levels were measured as fragments per kilobase of exon per million fragments mapped. A total of 1336 genes were differentially expressed between darkening and non-darkening cranberry bean RILs. Structural and regulatory genes of the proanthocyanidin biosynthesis pathway were upregulated in seed coats of the darkening RIL. A principal component analysis determined that changes in transcript levels for two genes of unknown function and three proanthocyanidin biosynthesis genes, FLAVANONE 3-HYDROXYLASE 1, DIHYDROFLAVONOL 4-REDUCTASE 1 and ANTHOCYANIDIN REDUCTASE 1 (PvANR1) were highly correlated with proanthocyanidin accumulation in seed coats of the darkening-susceptible cranberry bean RIL. HPLC-DAD analysis revealed that in vitro activity of a recombinant PvANR1 was NADPH-dependent and assays containing cyanidin yielded epicatechin and catechin; high cyanidin substrate levels inhibited the formation of both of these products. CONCLUSION: Proanthocyanidin oxidation is a pre-requisite for postharvest-related seed coat darkening in dicotyledonous seeds. In model plant species, the accumulation of proanthocyanidins is dependent upon upregulation of biosynthetic genes. In this study, proanthocyanidin production in cranberry bean seed coats was strongly associated with an increase in PvANR1 transcripts during seed maturation. In the presence of NADPH, PvANR1 converted the physiologically relevant substrate cyanidin to epicatechin and catechin.

Beyond the single gene: How epistasis and gene-by-environment effects influence crop domestication.

Domestication is a multifaceted evolutionary process, involving changes in individual genes, genetic interactions, and emergent phenotypes. There has been extensive discussion of the phenotypic characteristics of plant domestication, and recent research has started to identify the specific genes and mutational mechanisms that control domestication traits. However, there is an apparent disconnect between the simple genetic architecture described for many crop domestication traits, which should facilitate rapid phenotypic change under selection, and the slow rate of change reported from the archeobotanical record. A possible explanation involves the middle ground between individual genetic changes and their expression during development, where gene-by-gene (epistatic) and gene-by-environment interactions can modify the expression of phenotypes and opportunities for selection. These aspects of genetic architecture have the potential to significantly slow the speed of phenotypic evolution during crop domestication and improvement. Here we examine whether epistatic and gene-by-environment interactions have shaped how domestication traits have evolved. We review available evidence from the literature, and we analyze two domestication-related traits, shattering and flowering time, in a mapping population derived from a cross between domesticated foxtail millet and its wild progenitor. We find that compared with wild progenitor alleles, those favored during domestication often have large phenotypic effects and are relatively insensitive to genetic background and environmental effects. Consistent selection should thus be able to rapidly change traits during domestication. We conclude that if phenotypic evolution was slow during crop domestication, this is more likely due to cultural or historical factors than epistatic or environmental constraints.

Current perspectives and the future of domestication studies.

It is difficult to overstate the cultural and biological impacts that the domestication of plants and animals has had on our species. Fundamental questions regarding where, when, and how many times domestication took place have been of primary interest within a wide range of academic disciplines. Within the last two decades, the advent of new archaeological and genetic techniques has revolutionized our understanding of the pattern and process of domestication and agricultural origins that led to our modern way of life. In the spring of 2011, 25 scholars with a central interest in domestication representing the fields of genetics, archaeobotany, zooarchaeology, geoarchaeology, and archaeology met at the National Evolutionary Synthesis Center to discuss recent domestication research progress and identify challenges for the future. In this introduction to the resulting Special Feature, we present the state of the art in the field by discussing what is known about the spatial and temporal patterns of domestication, and controversies surrounding the speed, intentionality, and evolutionary aspects of the domestication process. We then highlight three key challenges for future research. We conclude by arguing that although recent progress has been impressive, the next decade will yield even more substantial insights not only into how domestication took place, but also when and where it did, and where and why it did not.

Changes in light quality alter physiological responses of soybean to thiamethoxam.

MAIN CONCLUSION: The interaction between neighboring weed-induced far-red enriched light and thiamethoxam can significantly alter soybean seedling morphology, nodulation, isoflavone levels, UV-absorbing phenolics, and carbon and nitrogen content. Neonicotinoid insecticides that are widely used on major crop plants can enhance plant growth and yield. Although the underlying mechanism of this enhanced growth and yield is not clear, recent studies suggest that neonicotinoids such as thiamethoxam (TMX) may exert their effects at least in part via signals that involve salicylic acid (SA) and jasmonic acid (JA). In the current research, effects of TMX on morphological and physiological responses of soybean have been compared under far-red-depleted (FR-D) and far-red-enriched (FR-E) light reflected by neighboring weeds. TMX significantly enhanced shoot and root growth but did not prevent stem elongation under FR-E light. Also, TMX did not prevent reductions in shoot carbon content and shoot carbon to nitrogen ratio under FR-E light. Despite similarities between these TMX effects in soybean and those known for SA and JA in other plant species, TMX significantly enhanced root-nodule numbers per plant and levels of root isoflavones malonyl-daidzin and malonyl-genistin under FR-E light only. These results suggest that the combined effect of FR-E light and TMX triggers a mechanism that operates concomitantly to enhance root isoflavones and nodulation in soybean.

Brevis plant1, a putative inositol polyphosphate 5-phosphatase, is required for internode elongation in maize.

In maize (Zea mays L.), as in other grass species, stem elongation occurs during growth and most noticeably upon the transition to flowering. Genes that reduce stem elongation have been important to reduce stem breakage, or lodging. Stem elongation has been mediated by dwarf and brachytic/brevis plant mutants that affect giberellic acid and auxin pathways, respectively. Maize brevis plant1 (bv1) mutants, first identified over 80 years ago, strongly resemble brachytic2 mutants that have shortened internodes, short internode cells, and are deficient in auxin transport. Here, we characterized two novel bv1 maize mutants. We found that an inositol polyphosphate 5-phosphatase orthologue of the rice gene dwarf50 was the molecular basis for the bv1 phenotype, implicating auxin-mediated inositol polyphosphate and/or phosphoinositide signalling in stem elongation. We suggest that auxin-mediated internode elongation involves processes that also contribute to stem gravitropism. Genes misregulated in bv1 mutants included genes important for cell wall synthesis, transmembrane transport, and cytoskeletal function. Mutant and wild-type plants were indistinguishable early in development, responded similarly to changes in light quality, had unaltered flowering times, and had normal flower development. These attributes suggest that breeding could utilize bv1 alleles to increase crop grain yields.

Evidence for selection on gene expression in cultivated rice (Oryza sativa).

Artificial selection has been used throughout plant domestication and breeding to develop crops that are adapted to diverse environments. Here, we investigate whether gene regulatory changes have been widespread targets of lineage-specific selection in cultivated lines Minghui 63 and Zhenshan 97 of rice, Oryza sativa. A line experiencing positive selection for either an increase or a decrease in genes' transcript abundances is expected to have an overabundance of expression quantitative trait locus (eQTL) alleles that increase or decrease those genes' expression, respectively. Results indicate that several genes that share Gene Ontology terms or are members of the same coexpression module have eQTL alleles from one parent that consistently increase gene expression relative to the second parent. A second line of evidence for lineage-specific selection is an overabundance of cis-trans pairs of eQTL alleles that affect gene expression in the same direction (are reinforcing). Across all cis-trans pairs of eQTL, including pairs that both weakly and strongly affect gene expression, there is no evidence for selection. However, the frequency of genes with reinforcing eQTL increases with eQTL strength. Therefore, there is evidence that eQTL with strong effects were positively selected during rice cultivation. Among 41 cis-trans pairs with strong trans eQTL, 31 have reinforcing eQTL. Several of the candidate genes under positive selection accurately predict phenotypic differences between Minghui 63 and Zhenshan 97. Overall, our results suggest that positive selection for regulatory alleles may be a key factor in plant improvement.

High throughput RNA sequencing of a hybrid maize and its parents shows different mechanisms responsive to nitrogen limitation.

BACKGROUND: Development of crop varieties with high nitrogen use efficiency (NUE) is crucial for minimizing N loss, reducing environmental pollution and decreasing input cost. Maize is one of the most important crops cultivated worldwide and its productivity is closely linked to the amount of fertilizer used. A survey of the transcriptomes of shoot and root tissues of a maize hybrid line and its two parental inbred lines grown under sufficient and limiting N conditions by mRNA-Seq has been conducted to have a better understanding of how different maize genotypes respond to N limitation. RESULTS: A different set of genes were found to be N-responsive in the three genotypes. Many biological processes important for N metabolism such as the cellular nitrogen compound metabolic process and the cellular amino acid metabolic process were enriched in the N-responsive gene list from the hybrid shoots but not from the parental lines' shoots. Coupled to this, sugar, carbohydrate, monosaccharide, glucose, and sorbitol transport pathways were all up-regulated in the hybrid, but not in the parents under N limitation. Expression patterns also differed between shoots and roots, such as the up-regulation of the cytokinin degradation pathway in the shoots of the hybrid and down-regulation of that pathway in the roots. The change of gene expression under N limitation in the hybrid resembled the parent with the higher NUE trait. The transcript abundances of alleles derived from each parent were estimated using polymorphic sites in mapped reads in the hybrid. While there were allele abundance differences, there was no correlation between these and the expression differences seen between the hybrid and the two parents. CONCLUSIONS: Gene expression in two parental inbreds and the corresponding hybrid line in response to N limitation was surveyed using the mRNA-Seq technology. The data showed that the three genotypes respond very differently to N-limiting conditions, and the hybrid clearly has a unique expression pattern compared to its parents. Our results expand our current understanding of N responses and will help move us forward towards effective strategies to improve NUE and enhance crop production.

Protein-coding cis-natural antisense transcripts have high and broad expression in Arabidopsis.

Pairs of genes within eukaryotic genomes are often located on opposite DNA strands such that transcription generates cis-natural sense antisense transcripts (cis-NATs). This orientation of genes has been associated with the biogenesis of splice variants and natural antisense small RNAs. Here, in an analysis of currently available data, we report that within Arabidopsis (Arabidopsis thaliana), protein-coding cis-NATs are also characterized by high abundance, high coexpression, and broad expression. Our results suggest that a permissive chromatin environment may have led to the proximity of these genes. Compared with other genes, cis-NAT-encoding genes have enriched low-nucleosome-density regions, high levels of histone H3 lysine-9 acetylation, and low levels of H3 lysine-27 trimethylation. Promoters associated with broadly expressed genes are preferentially found in the 5' regulatory sequences of cis-NAT-encoding genes. Our results further suggest that natural antisense small RNA production from cis-NATs is limited. Small RNAs sequenced from natural antisense small RNA biogenesis mutants including dcl1, dcl2, dcl3, and rdr6 map to cis-NATs as frequently as small RNAs sequenced from wild-type plants. Future work will investigate if the positive transcriptional regulation of overlapping protein-coding genes contributes to the prevalence of these genes within other eukaryotic genomes.

A developmental transcriptional network for maize defines coexpression modules.

Here, we present a genome-wide overview of transcriptional circuits in the agriculturally significant crop species maize (Zea mays). We examined transcript abundance data at 50 developmental stages, from embryogenesis to senescence, for 34,876 gene models and classified genes into 24 robust coexpression modules. Modules were strongly associated with tissue types and related biological processes. Sixteen of the 24 modules (67%) have preferential transcript abundance within specific tissues. One-third of modules had an absence of gene expression in specific tissues. Genes within a number of modules also correlated with the developmental age of tissues. Coexpression of genes is likely due to transcriptional control. For a number of modules, key genes involved in transcriptional control have expression profiles that mimic the expression profiles of module genes, although the expression of transcriptional control genes is not unusually representative of module gene expression. Known regulatory motifs are enriched in several modules. Finally, of the 13 network modules with more than 200 genes, three contain genes that are notably clustered (P < 0.05) within the genome. This work, based on a carefully selected set of major tissues representing diverse stages of maize development, demonstrates the remarkable power of transcript-level coexpression networks to identify underlying biological processes and their molecular components.

Regulatory motifs identified from a maize developmental coexpression network.

Transcriptional control is an important determinant of plant development, and distinct modules of coordinated genes characterize the maize developmental transcriptome. Upstream regulatory sequences are often the primary factors that control gene expression pattern and abundance. Here, we identify 244 regulatory motifs that are significantly enriched within 24 gene expression modules previously constructed from transcript abundances of 34,876 Zea mays (maize) gene models from embryogenesis to senescence. Within modules, we identify motifs that have not been characterized. In addition, we identify motifs similar to experimentally verified motifs, and the functions of these motifs overlap with predicted module functions. This work demonstrates the power of transcript-level coexpression modules to identify both variants of known regulatory motifs and novel motifs that control a species' developmental transcriptome.

The origin of the naked grains of maize.

The most critical step in maize (Zea mays ssp. mays) domestication was the liberation of the kernel from the hardened, protective casing that envelops the kernel in the maize progenitor, teosinte. This evolutionary step exposed the kernel on the surface of the ear, such that it could readily be used by humans as a food source. Here we show that this key event in maize domestication is controlled by a single gene (teosinte glume architecture or tga1), belonging to the SBP-domain family of transcriptional regulators. The factor controlling the phenotypic difference between maize and teosinte maps to a 1-kilobase region, within which maize and teosinte show only seven fixed differences in their DNA sequences. One of these differences encodes a non-conservative amino acid substitution and may affect protein function, and the other six differences potentially affect gene regulation. Molecular evolution analyses show that this region was the target of selection during maize domestication. Our results demonstrate that modest genetic changes in single genes can induce dramatic changes in phenotype during domestication and evolution.