Edward S. Buckler
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Edward S. Buckler has a total of 1000 co-authors. Here are the top co-authors by number of shared publications:
- Maria Cinta Romay
Institute for Genomic Diversity, Cornell University
57 shared publications - Michael Allen Gore
Plant Breeding and Genetics Section, School of Integrative Plant Science, Cornell University, Ithaca, NY 14853, USA
38 shared publications - James Brendan Holland
Department of Crop and Soil Sciences North Carolina State University Raleigh North Carolina USA
32 shared publications - Zhiwu Zhang
Washington State University, USA
28 shared publications - Sherry Flint-Garcia
USDA-ARS, University of Missouri
25 shared publications - Peter J. Bradbury
USDA
22 shared publications
Latest Publications
- Open Access
- 14 October 2026
- Published by Public Library of Science (PLoS) in PLOS ONE
- Vol. 21 (10) , e0358458
Evolutionary constraint on a core winter proteome across independently cold-adapted PACMAD grasses
Grasses in the Panicoideae, Aristidoideae, Chloridoideae, Micrairoideae, Arundinoideae, and Danthonioideae (PACMAD) clade independently colonized cold environments from warm-climate ancestors, but whether their molecular responses to freezing reflect shared evolutionary solutions or lineage-specific innovations remains poorly understood at the protein level. We used comparative proteomics to test whether protein-level cold responses show stronger cross-species conservation than previously observed at the transcript level. We quantified seasonal rhizome proteomes (winter vs summer) from five PACMAD species grown in a common-garden experiment exposed to sustained sub-zero temperatures, identified differentially abundant proteins, and compared fold-change magnitudes across species using orthogroup-based correlation analyses. We further examined LEA3 protein structure through hydropathy profiling and motif analysis. Shared cold-responsive proteins showed higher cross-species fold-change correlation ( ρ = 0.80) than background proteins ( ρ = 0.45), despite greater divergence in baseline abundance. This correlation advantage persisted (Δ ρ = 0.27) after controlling for effect-size selection by comparing shared cold-responsive proteins against magnitude-matched single-species responses. LEA3 was the only ortholog elevated across all five species, with cold-tolerant species containing more tandem 11-mer repeats than the cold-sensitive maize. Two species accumulated multiple LEA3 paralogs, increasing total LEA3 abundance. Independent evolution of freezing tolerance in PACMAD grasses is consistent with evolutionary constraints on protein-level response magnitude, reflecting the retention of an ancestral protective capacity. Structural divergence of LEA3 in maize suggests that transcriptional induction alone does not ensure freezing tolerance; functional protection likely requires intact motif architecture.
- 14 October 2026
- Published by Public Library of Science (PLoS) in PLOS ONE
- Vol. 21 (10) , e0358458
Open AccessEvolutionary constraint on a core winter proteome across independently cold-adapted PACMAD grasses
Grasses in the Panicoideae, Aristidoideae, Chloridoideae, Micrairoideae, Arundinoideae, and Danthonioideae (PACMAD) clade independently colonized cold environments from warm-climate ancestors, but whether their molecular responses to freezing reflect shared evolutionary solutions or lineage-specific innovations remains poorly understood at the protein level. We used comparative proteomics to test whether protein-level cold responses show stronger cross-species conservation than previously observed at the transcript level. We quantified seasonal rhizome proteomes (winter vs summer) from five PACMAD species grown in a common-garden experiment exposed to sustained sub-zero temperatures, identified differentially abundant proteins, and compared fold-change magnitudes across species using orthogroup-based correlation analyses. We further examined LEA3 protein structure through hydropathy profiling and motif analysis. Shared cold-responsive proteins showed higher cross-species fold-change correlation ( ρ = 0.80) than background proteins ( ρ = 0.45), despite greater divergence in baseline abundance. This correlation advantage persisted (Δ ρ = 0.27) after controlling for effect-size selection by comparing shared cold-responsive proteins against magnitude-matched single-species responses. LEA3 was the only ortholog elevated across all five species, with cold-tolerant species containing more tandem 11-mer repeats than the cold-sensitive maize. Two species accumulated multiple LEA3 paralogs, increasing total LEA3 abundance. Independent evolution of freezing tolerance in PACMAD grasses is consistent with evolutionary constraints on protein-level response magnitude, reflecting the retention of an ancestral protective capacity. Structural divergence of LEA3 in maize suggests that transcriptional induction alone does not ensure freezing tolerance; functional protection likely requires intact motif architecture.
- Open Access
An open field phenomics resource for multimodal maize yield prediction across divergent environments
Temporal drone phenotyping captures crop development, but irregular flight schedules complicate comparisons across environments. We release curated imagery from 356 flights across 19 Genomes to Fields environments containing 1,180 maize ( Zea mays L.) hybrids. To evaluate its utility, we integrated functional principal components of vegetation index and weather trajectories with genomic information. Combinined genomic and phenomic kernels improved yield prediction, reaching correlations up to r = 0.501 for held-out hybrids in environments represented in training and 0.408 when environments were also withheld. Accumulated growing degree days offered no consistent predictive advantage over days after planting, and weather contributed modest, task-dependent gains. A transformer neural process learned directly from irregular observations, serving as a novel application of neural process models in agriculture. Mapping vegetation index functional principal components identified recurrent quantitative trait loci on chromosomes 3 and 7. This resource and its reproducible analyses guide the use of temporal spectral data for crop prediction and genetic discovery.
Open AccessAn open field phenomics resource for multimodal maize yield prediction across divergent environments
Temporal drone phenotyping captures crop development, but irregular flight schedules complicate comparisons across environments. We release curated imagery from 356 flights across 19 Genomes to Fields environments containing 1,180 maize ( Zea mays L.) hybrids. To evaluate its utility, we integrated functional principal components of vegetation index and weather trajectories with genomic information. Combinined genomic and phenomic kernels improved yield prediction, reaching correlations up to r = 0.501 for held-out hybrids in environments represented in training and 0.408 when environments were also withheld. Accumulated growing degree days offered no consistent predictive advantage over days after planting, and weather contributed modest, task-dependent gains. A transformer neural process learned directly from irregular observations, serving as a novel application of neural process models in agriculture. Mapping vegetation index functional principal components identified recurrent quantitative trait loci on chromosomes 3 and 7. This resource and its reproducible analyses guide the use of temporal spectral data for crop prediction and genetic discovery.
- Open Access
- 8 September 2026
- Published by National Academy of Sciences in Proceedings of the National Academy of Sciences
- Vol. 123 (37)
Evolution of maize recombination landscape during domestication
Despite the plethora of knowledge about the benefits of meiotic recombination and numerous theoretical studies examining how recombination rates evolve, there is a general lack of empirical support and consensus across species. To fill this knowledge gap, we characterized the evolution of recombination landscape in maize during its domestication from teosinte and related the observed changes to established theoretical frameworks. Through examining recombination in experimental populations of maize and teosinte and the population genomics approach of identifying historical recombination events using ancestral recombination graph inference to generate saturated maize and teosinte recombination maps, we found that during domestication, maize experienced a 12% increase in its genome-wide recombination rate. Furthermore, maize evolved higher recombination rates on the long arms of chromosomes in regions closer to centromeres, where recombination is generally very low. The repatterning of crossover events came from changes in global crossover positioning rather than alterations in cis-acting chromatin factors. Consequently, we found evidence of selection acting on trans-acting recombination modifiers affecting crossover interference and controlling the interference-dependent class I crossover pathway. We show that CO repatterning was likely beneficial for maize fitness, as significant recombination rate increases were predominantly in gene-rich regions, which harbor domestication-related variation. This work suggests genomic and mechanistic processes leading to the evolution of meiotic recombination landscape in response to directional selection pressure and provides evidence for the evolutionary advantage of recombination.
- 8 September 2026
- Published by National Academy of Sciences in Proceedings of the National Academy of Sciences
- Vol. 123 (37)
Open AccessEvolution of maize recombination landscape during domestication
Despite the plethora of knowledge about the benefits of meiotic recombination and numerous theoretical studies examining how recombination rates evolve, there is a general lack of empirical support and consensus across species. To fill this knowledge gap, we characterized the evolution of recombination landscape in maize during its domestication from teosinte and related the observed changes to established theoretical frameworks. Through examining recombination in experimental populations of maize and teosinte and the population genomics approach of identifying historical recombination events using ancestral recombination graph inference to generate saturated maize and teosinte recombination maps, we found that during domestication, maize experienced a 12% increase in its genome-wide recombination rate. Furthermore, maize evolved higher recombination rates on the long arms of chromosomes in regions closer to centromeres, where recombination is generally very low. The repatterning of crossover events came from changes in global crossover positioning rather than alterations in cis-acting chromatin factors. Consequently, we found evidence of selection acting on trans-acting recombination modifiers affecting crossover interference and controlling the interference-dependent class I crossover pathway. We show that CO repatterning was likely beneficial for maize fitness, as significant recombination rate increases were predominantly in gene-rich regions, which harbor domestication-related variation. This work suggests genomic and mechanistic processes leading to the evolution of meiotic recombination landscape in response to directional selection pressure and provides evidence for the evolutionary advantage of recombination.
- Open Access
- 7 July 2026
- Published by Springer Nature in BMC Biology
Phylogenetic analysis and machine learning identify signatures of selection and predict deleterious mutations in common bean (Phaseolus vulgaris L.)
Mutations are a continuous source of new alleles and genetic diversity in populations. Domestication and selection influence the accumulation of alleles oc
- 7 July 2026
- Published by Springer Nature in BMC Biology
Open AccessPhylogenetic analysis and machine learning identify signatures of selection and predict deleterious mutations in common bean (Phaseolus vulgaris L.)
Mutations are a continuous source of new alleles and genetic diversity in populations. Domestication and selection influence the accumulation of alleles oc
- Open Access
Convergent genome- and gene-level constraints shape repeated environmental adaptation in grasses
Summary: Grasses (Poaceae) dominate terrestrial ecosystems and sustain global food security, yet the genomic principles enabling their repeated adaptation to extreme environments remain unresolved. Combining dense phylogenomic sampling, global environmental data, and genomic large language models (gLLMs), we characterize the mutational targets underlying environmental shift across 687 genomes from 561 species spanning 17 climate zones. We identify 17-30 phylogenetically independent transitions into extreme temperature, water, and soil environments, accompanied by convergent shifts in genome-scale molecular properties, including the Nitrogen-to-Carbon balance and the biosynthetic cost of the proteome. Our gLLMs-informed phylogenetic mixed modeling framework identifies 330 genes that repeatedly associated with distinct axes of environmental variation, highlighting the importance of protein modification and localization within extracellular and organellar compartments. Overlaying independent convergent selection tests identifies 27 high-confidence candidates for further characterization. Together, our results show that grass adaptation is canalized by layered constraints at genome-wide and gene-specific scales, producing predictable evolutionary trajectories.
Open AccessConvergent genome- and gene-level constraints shape repeated environmental adaptation in grasses
Summary: Grasses (Poaceae) dominate terrestrial ecosystems and sustain global food security, yet the genomic principles enabling their repeated adaptation to extreme environments remain unresolved. Combining dense phylogenomic sampling, global environmental data, and genomic large language models (gLLMs), we characterize the mutational targets underlying environmental shift across 687 genomes from 561 species spanning 17 climate zones. We identify 17-30 phylogenetically independent transitions into extreme temperature, water, and soil environments, accompanied by convergent shifts in genome-scale molecular properties, including the Nitrogen-to-Carbon balance and the biosynthetic cost of the proteome. Our gLLMs-informed phylogenetic mixed modeling framework identifies 330 genes that repeatedly associated with distinct axes of environmental variation, highlighting the importance of protein modification and localization within extracellular and organellar compartments. Overlaying independent convergent selection tests identifies 27 high-confidence candidates for further characterization. Together, our results show that grass adaptation is canalized by layered constraints at genome-wide and gene-specific scales, producing predictable evolutionary trajectories.
- 20 May 2026
- Published by Springer Nature in Nature Reviews Genetics
- Vol. 27 (10) , 797-812
Translating functional molecular knowledge into crop-breeding success
Historical plant breeding, which optimizes phenotypes through selective crossing guided by phenotypic evaluation and molecular markers, is limited by evolutionary constraints that hinder rapid crop improvement. A new paradigm, precision breeding, circumvents these limitations by targeting genetic variants through functional molecular knowledge. To generate this knowledge at scale, sequence-based deep learning leverages high-quality genome sequence data to predict variant effects at base-pair resolution. When linked to agronomically important traits, these predictions enable breeders to prioritize variants for precision selection or editing. Although it is still in the early stages of development, we foresee three key applications for this approach: introgressing genes from distant breeding pools, purging deleterious mutations and designing new plant ideotypes. Looking ahead, refined computational models will facilitate targeted editing and the systematic redesign of complex physiological processes to address emerging breeding goals under shifting environmental conditions.
- 20 May 2026
- Published by Springer Nature in Nature Reviews Genetics
- Vol. 27 (10) , 797-812
Translating functional molecular knowledge into crop-breeding success
Historical plant breeding, which optimizes phenotypes through selective crossing guided by phenotypic evaluation and molecular markers, is limited by evolutionary constraints that hinder rapid crop improvement. A new paradigm, precision breeding, circumvents these limitations by targeting genetic variants through functional molecular knowledge. To generate this knowledge at scale, sequence-based deep learning leverages high-quality genome sequence data to predict variant effects at base-pair resolution. When linked to agronomically important traits, these predictions enable breeders to prioritize variants for precision selection or editing. Although it is still in the early stages of development, we foresee three key applications for this approach: introgressing genes from distant breeding pools, purging deleterious mutations and designing new plant ideotypes. Looking ahead, refined computational models will facilitate targeted editing and the systematic redesign of complex physiological processes to address emerging breeding goals under shifting environmental conditions.
- Open Access
A transcriptomic atlas of grass senescence reveals divergent underground sink networks limit nitrogen recycling in annuals
Senescence enables plants to remobilize and recycle nutrients from aging organs to support growth, reproduction, and survival. In annual crops like maize, nitrogen remobilization from leaves to grain is incomplete, with 30-50% of nitrogen stranded in aboveground tissues and subject to environmental loss. Mitigating nitrogen loss in annual crops could be achieved by leveraging the physiological strategies of perennial grasses, which remobilize nitrogen and other nutrients into underground organs at the end of the growing season, thereby preventing environmental leakage. To uncover the molecular basis of perennial nitrogen recycling to underground organs, we built a transcriptomic atlas from field-grown plants, comprising 2,685 RNA-seq libraries from 14 grass species within the Panicoideae (Poaceae), utilizing maize and sorghum as annual references for comparative analyses. The atlas spans leaves, roots, stalks, and rhizomes across two seasons, from mid-growing season to senescence. Using a photosynthetic index to align the leaf’s transition from nitrogen sink to source across species, co-expression network analysis revealed that the subnetworks driving leaf nitrogen recycling are preserved across annuals and perennials. However, we discovered that the subnetworks associated with underground sink establishment, specifically those associated with seed-like dormancy and desiccation tolerance pathways, have diverged among annual crop accessions. Our work identifies conserved gene candidates and networks that could be used to reintroduce perennial-like nutrient recycling into annual crops to enhance long-term nutrient retention in the field.
Open AccessA transcriptomic atlas of grass senescence reveals divergent underground sink networks limit nitrogen recycling in annuals
Senescence enables plants to remobilize and recycle nutrients from aging organs to support growth, reproduction, and survival. In annual crops like maize, nitrogen remobilization from leaves to grain is incomplete, with 30-50% of nitrogen stranded in aboveground tissues and subject to environmental loss. Mitigating nitrogen loss in annual crops could be achieved by leveraging the physiological strategies of perennial grasses, which remobilize nitrogen and other nutrients into underground organs at the end of the growing season, thereby preventing environmental leakage. To uncover the molecular basis of perennial nitrogen recycling to underground organs, we built a transcriptomic atlas from field-grown plants, comprising 2,685 RNA-seq libraries from 14 grass species within the Panicoideae (Poaceae), utilizing maize and sorghum as annual references for comparative analyses. The atlas spans leaves, roots, stalks, and rhizomes across two seasons, from mid-growing season to senescence. Using a photosynthetic index to align the leaf’s transition from nitrogen sink to source across species, co-expression network analysis revealed that the subnetworks driving leaf nitrogen recycling are preserved across annuals and perennials. However, we discovered that the subnetworks associated with underground sink establishment, specifically those associated with seed-like dormancy and desiccation tolerance pathways, have diverged among annual crop accessions. Our work identifies conserved gene candidates and networks that could be used to reintroduce perennial-like nutrient recycling into annual crops to enhance long-term nutrient retention in the field.
- Open Access
Inter-variety competition dynamics in US inbred and hybrid maize
Variety mixtures provide a potential avenue in US cropping systems to improve yield stability and disease resistance. However, implementation of variety mixtures requires an understanding of the competitive dynamics of the crop. In this study, we examine the effects of plant competition both between and within plots through five unique experiments: 1) 5,000 diverse inbred lines in single-row plots, 2) hybrids in two-row plots developed from the above inbred lines, 3) over 4,000 hybrids measured in 141 locations in two-row plots as part of Genomes to Fields, 4) mixtures of two hybrids within a two-row plot planted across two years and five locations, and 5) mixtures of up to twenty hybrids in four-row plots in three locations. Across all experiments, we find that competitive interactions are extremely limited. Within inbred lines, height of the neighboring plot accounts for 1.2% of the variance in focal plot height. Similarly, neighbor height explains 1.7% of the variance in focal plot yield in hybrids developed from the inbred lines. The genetics of neighboring plots explains 1.55% of the variation in yield across 141 location-year environments, reinforcing the generally modest impacts of neighbor competition. In evaluating mixtures of hybrids in both two and four-row plots, we observe no yield penalty compared to conventional single hybrid plots, even with large height differentials of the hybrids included in the mixture or in mixtures of up to 20 hybrids within a plot. Finally, we observe that mixtures have more yield stability compared to conventional plots, highlighting a new avenue for increased stability in higher risk environments. The lack of yield penalty and stability benefits are promising for future investigations of mixtures that may complement each other in disease resistance or abiotic stress tolerance and increase overall yield stability in the field.
Open AccessInter-variety competition dynamics in US inbred and hybrid maize
Variety mixtures provide a potential avenue in US cropping systems to improve yield stability and disease resistance. However, implementation of variety mixtures requires an understanding of the competitive dynamics of the crop. In this study, we examine the effects of plant competition both between and within plots through five unique experiments: 1) 5,000 diverse inbred lines in single-row plots, 2) hybrids in two-row plots developed from the above inbred lines, 3) over 4,000 hybrids measured in 141 locations in two-row plots as part of Genomes to Fields, 4) mixtures of two hybrids within a two-row plot planted across two years and five locations, and 5) mixtures of up to twenty hybrids in four-row plots in three locations. Across all experiments, we find that competitive interactions are extremely limited. Within inbred lines, height of the neighboring plot accounts for 1.2% of the variance in focal plot height. Similarly, neighbor height explains 1.7% of the variance in focal plot yield in hybrids developed from the inbred lines. The genetics of neighboring plots explains 1.55% of the variation in yield across 141 location-year environments, reinforcing the generally modest impacts of neighbor competition. In evaluating mixtures of hybrids in both two and four-row plots, we observe no yield penalty compared to conventional single hybrid plots, even with large height differentials of the hybrids included in the mixture or in mixtures of up to 20 hybrids within a plot. Finally, we observe that mixtures have more yield stability compared to conventional plots, highlighting a new avenue for increased stability in higher risk environments. The lack of yield penalty and stability benefits are promising for future investigations of mixtures that may complement each other in disease resistance or abiotic stress tolerance and increase overall yield stability in the field.
- Open Access
Identification of Freezing Tolerance QTLs in Tripsacum dactyloides Using Open-Pollinated Bulk Segregant Analysis
This study investigates the genetic basis of freezing tolerance in Tripsacum dactyloides and related subspecies as a potential source of valuable traits for improving maize agriculture. Recognizing the significant economic losses in corn yields due to frost damage, we hypothesized that northern populations of T. dactyloides are enriched for freezing tolerance alleles. 40 diverse Tripsacum accessions were collected from natural populations and long-established field collections and used to generate F1 hybrids and open-pollinated F2 families. F2 seedlings were germinated then screened within a growth chamber for freezing tolerance by exposure to freezing temperatures. Seedlings were then phenotyped by tissue survival, and extremes were pooled to create tolerant and susceptible bulks. DNA sequencing was performed on founders, F1s, and tolerant/susceptible F2 bulks. To overcome challenges in traditional SNP calling in bulked samples, we developed a regression-based approach to estimate gamete frequencies and impute allele frequencies in pooled populations. The results showed genetic diversity among Tripsacum accessions, with divergence between northern and southern populations. We tracked segregation of alleles across genomic loci, and performed a joint bulk segregant analysis, identifying 9 QTLs significantly associated with freezing tolerance. These findings highlight potential loci for freezing tolerance that could inform genetic engineering of maize. Central Hypothesis: Northern populations of Tripsacum dactyloides , a wild relative of maize, are enriched for freezing tolerance alleles which can be identified by mapping.
Open AccessIdentification of Freezing Tolerance QTLs in Tripsacum dactyloides Using Open-Pollinated Bulk Segregant Analysis
This study investigates the genetic basis of freezing tolerance in Tripsacum dactyloides and related subspecies as a potential source of valuable traits for improving maize agriculture. Recognizing the significant economic losses in corn yields due to frost damage, we hypothesized that northern populations of T. dactyloides are enriched for freezing tolerance alleles. 40 diverse Tripsacum accessions were collected from natural populations and long-established field collections and used to generate F1 hybrids and open-pollinated F2 families. F2 seedlings were germinated then screened within a growth chamber for freezing tolerance by exposure to freezing temperatures. Seedlings were then phenotyped by tissue survival, and extremes were pooled to create tolerant and susceptible bulks. DNA sequencing was performed on founders, F1s, and tolerant/susceptible F2 bulks. To overcome challenges in traditional SNP calling in bulked samples, we developed a regression-based approach to estimate gamete frequencies and impute allele frequencies in pooled populations. The results showed genetic diversity among Tripsacum accessions, with divergence between northern and southern populations. We tracked segregation of alleles across genomic loci, and performed a joint bulk segregant analysis, identifying 9 QTLs significantly associated with freezing tolerance. These findings highlight potential loci for freezing tolerance that could inform genetic engineering of maize. Central Hypothesis: Northern populations of Tripsacum dactyloides , a wild relative of maize, are enriched for freezing tolerance alleles which can be identified by mapping.
- Open Access
- 10 December 2025
- Published by Oxford University Press (OUP) in Molecular Biology and Evolution
- Vol. 43 (1)
Widespread turnover of a conserved cis - regulatory code across 589 grass species
The growing availability of genomes from non-model organisms offers new opportunities to identify functional loci underlying trait variation through comparative genomics. While cis-regulatory regions drive much of phenotypic evolution, linking them to specific functions remains challenging. We identified 514 cis-regulatory motifs enriched in regulatory regions of five diverse grass species, with 73% consistently enriched across all, suggesting a deeply conserved regulatory code. Leveraging 57 new contig-level genome assemblies, we then quantified shared occupancy of specific motif instances within gene-proximal regions across 589 grass species, revealing widespread gain and loss over evolutionary time. Shared occupancy declined rapidly over the first few million years of divergence, yet ∼50% of motif instances were shared back to the origin of grasses ∼100 million years ago. We used phylogenetic mixed models to identify motif gains and losses associated with ecological niche transitions. Our models revealed significant environmental associations across 1282 motif–orthogroup combinations, including convergent gains of HSF/GARP motifs at an alpha-N-acetylglucosaminidase gene associated with occurrence in temperate environments. Our findings support a “stable motifs, variable binding sites” model in which cis-regulatory evolution involves turnover of thousands of individual binding site instances while largely preserving transcription factors’ binding preferences. Our results highlight the potential of comparative genomics and phylogenetic mixed models to reveal the genetic basis of complex traits.
- 10 December 2025
- Published by Oxford University Press (OUP) in Molecular Biology and Evolution
- Vol. 43 (1)
Open AccessWidespread turnover of a conserved cis - regulatory code across 589 grass species
The growing availability of genomes from non-model organisms offers new opportunities to identify functional loci underlying trait variation through comparative genomics. While cis-regulatory regions drive much of phenotypic evolution, linking them to specific functions remains challenging. We identified 514 cis-regulatory motifs enriched in regulatory regions of five diverse grass species, with 73% consistently enriched across all, suggesting a deeply conserved regulatory code. Leveraging 57 new contig-level genome assemblies, we then quantified shared occupancy of specific motif instances within gene-proximal regions across 589 grass species, revealing widespread gain and loss over evolutionary time. Shared occupancy declined rapidly over the first few million years of divergence, yet ∼50% of motif instances were shared back to the origin of grasses ∼100 million years ago. We used phylogenetic mixed models to identify motif gains and losses associated with ecological niche transitions. Our models revealed significant environmental associations across 1282 motif–orthogroup combinations, including convergent gains of HSF/GARP motifs at an alpha-N-acetylglucosaminidase gene associated with occurrence in temperate environments. Our findings support a “stable motifs, variable binding sites” model in which cis-regulatory evolution involves turnover of thousands of individual binding site instances while largely preserving transcription factors’ binding preferences. Our results highlight the potential of comparative genomics and phylogenetic mixed models to reveal the genetic basis of complex traits.