Adaptation and the Geographic Spread of Crop Species

Crops are plant species that were domesticated starting about 11,000 years ago from several centers of origin, most prominently the Fertile Crescent, East Asia, and Mesoamerica. From their domestication centers, these crops spread across the globe and had to adapt to differing environments as a result of this dispersal. We discuss broad patterns of crop spread, including the early diffusion of crops associated with the rise and spread of agriculture, the later movement via ancient trading networks, and the exchange between the Old and New Worlds over the last ∼550 years after the European colonization of the Americas. We also examine the various genetic mechanisms associated with the evolutionary adaptation of crops to their new environments after dispersal, most prominently seasonal adaptation associated with movement across latitudes, as well as altitudinal, temperature, and other environmental factors.

Seedling root system adaptation to water availability during maize domestication and global expansion

The maize root system has been reshaped by indirect selection during global adaptation to new agricultural environments. In this study, we characterized the root systems of more than 9,000 global maize accessions and its wild relatives, defining the geographical signature and genomic basis of variation in seminal root number. We demonstrate that seminal root number has increased during maize domestication followed by a decrease in response to limited water availability in locally adapted varieties. By combining environmental and phenotypic association analyses with linkage mapping, we identified genes linking environmental variation and seminal root number. Functional characterization of the transcription factor ZmHb77 and in silico root modeling provides evidence that reshaping root system architecture by reducing the number of seminal roots and promoting lateral root density is beneficial for the resilience of maize seedlings to drought.

Are cereal grasses a single genetic system?

In 1993, a passionate and provocative call to arms urged cereal researchers to consider the taxon they study as a single genetic system and collaborate with each other. Since then, that group of scientists has seen their discipline blossom. In an attempt to understand what unity of genetic systems means and how the notion was borne out by later research, we survey the progress and prospects of cereal genomics: sequence assemblies, population-scale sequencing, resistance gene cloning and domestication genetics. Gene order may not be as extraordinarily well conserved in the grasses as once thought. Still, several recurring themes have emerged. The same ancestral molecular pathways defining plant architecture have been co-opted in the evolution of different cereal crops. Such genetic convergence as much as cross-fertilization of ideas between cereal geneticists has led to a rich harvest of genes that, it is hoped, will lead to improved varieties.

Unraveling the diversity and cultural heritage of fruit crops through paleogenomics

Abundant and plentiful fruit crops are threatened by the loss of diverse legacy cultivars which are being replaced by a limited set of high-yielding ones. This article delves into the potential of paleogenomics that utilizes ancient DNA analysis to revive lost diversity. By focusing on grapevines, date palms, and tomatoes, recent studies showcase the effectiveness of paleogenomic techniques in identifying and understanding genetic traits crucial for crop resilience, disease resistance, and nutritional value. The approach not only tracks landrace dispersal and introgression but also sheds light on domestication events. In the face of major future environmental challenges, integrating paleogenomics with modern breeding strategies emerges as a promising avenue to significantly bolster fruit crop sustainability.

ZmELF6-ZmPRR37 module regulates maize flowering and salt response

The control of flowering time in maize is crucial for reproductive success and yield, and it can be influenced by environmental stresses. Using the approaches of Ac/Ds transposon and transposable element amplicon sequencing techniques, we identified a Ds insertion mutant in the ZmPRR37 gene. The Ds insertion showed a significant correlation with days to anthesis. Further research indicated that ZmPRR37-CR knockout mutants exhibited early flowering, whereas ZmPRR37-overexpression lines displayed delayed flowering compared to WT under long-day (LD) conditions. We demonstrated that ZmPRR37 repressed the expression of ZmNF-YC2 and ZmNF-YA3 to delay flowering. Association analysis revealed a significant correlation between flowering time and a SNP2071-C/T located upstream of ZmPRR37. The SNP2071-C/T impacted the binding capacity of ZmELF6 to the promoter of ZmPRR37. ZmELF6 also acted as a flowering suppressor in maize under LD conditions. Notably, our study unveiled that ZmPRR37 can enhance salt stress tolerance in maize by directly regulating the expression of ABA-responsive gene ZmDhn1. ZmDhn1 negatively regulated maize salt stress resistance. In summary, our findings proposed a novel pathway for regulating photoperiodic flowering and responding to salt stress based on ZmPRR37 in maize, providing novel insights into the integration of abiotic stress signals into floral pathways.

Cytoplasmic genome contributions to domestication and improvement of modern maize

BACKGROUND

Studies on maize evolution and domestication are largely limited to the nuclear genomes, and the contribution of cytoplasmic genomes to selection and domestication of modern maize remains elusive. Maize cytoplasmic genomes have been classified into fertile (NA and NB) and cytoplasmic-nuclear male-sterility (CMS-S, CMS-C, and CMS-T) groups, but their contributions to modern maize breeding have not been systematically investigated.

RESULTS

Here we report co-selection and convergent evolution between nuclear and cytoplasmic genomes by analyzing whole genome sequencing data of 630 maize accessions modern maize and its relatives, including 24 fully assembled mitochondrial and chloroplast genomes. We show that the NB cytotype is associated with the expansion of modern maize to North America, gradually replaces the fertile NA cytotype probably through unequal division, and predominates in over 90% of modern elite inbred lines. The mode of cytoplasmic evolution is increased nucleotypic diversity among the genes involved in photosynthesis and energy metabolism, which are driven by selection and domestication. Furthermore, genome-wide association study reveals correlation of cytoplasmic nucleotypic variation with key agronomic and reproductive traits accompanied with the diversification of the nuclear genomes.

CONCLUSIONS

Our results indicate convergent evolution between cytoplasmic and nuclear genomes during maize domestication and breeding. These new insights into the important roles of mitochondrial and chloroplast genomes in maize domestication and improvement should help select elite inbred lines to improve yield stability and crop resilience of maize hybrids.

Genomic insight into the origin, domestication, dispersal, diversification and human selection of Tartary buckwheat

BACKGROUND

Tartary buckwheat, Fagopyrum tataricum, is a pseudocereal crop with worldwide distribution and high nutritional value. However, the origin and domestication history of this crop remain to be elucidated.

RESULTS

Here, by analyzing the population genomics of 567 accessions collected worldwide and reviewing historical documents, we find that Tartary buckwheat originated in the Himalayan region and then spread southwest possibly along with the migration of the Yi people, a minority in Southwestern China that has a long history of planting Tartary buckwheat. Along with the expansion of the Mongol Empire, Tartary buckwheat dispersed to Europe and ultimately to the rest of the world. The different natural growth environments resulted in adaptation, especially significant differences in salt tolerance between northern and southern Chinese Tartary buckwheat populations. By scanning for selective sweeps and using a genome-wide association study, we identify genes responsible for Tartary buckwheat domestication and differentiation, which we then experimentally validate. Comparative genomics and QTL analysis further shed light on the genetic foundation of the easily dehulled trait in a particular variety that was artificially selected by the Wa people, a minority group in Southwestern China known for cultivating Tartary buckwheat specifically for steaming as a staple food to prevent lysine deficiency.

CONCLUSIONS

This study provides both comprehensive insights into the origin and domestication of, and a foundation for molecular breeding for, Tartary buckwheat.

Tracking population structure and phenology through time using ancient genomes from waterlogged white oak wood

Whole genome characterizations of crop plants based on ancient DNA have provided unique keys for a better understanding of the evolutionary origins of modern cultivars, the pace and mode of selection underlying their adaptation to new environments and the production of phenotypes of interest. Although forests are among the most biologically rich ecosystems on earth and represent a fundamental resource for human societies, no ancient genome sequences have been generated for trees. This contrasts with the generation of multiple ancient reference genomes for important crops. Here, we sequenced the first ancient tree genomes using two white oak wood remains from Germany dating to the Last Little Ice Age (15th century CE, 7.3× and 4.0×) and one from France dating to the Bronze Age (1700 BCE, 3.4×). We assessed the underlying species and identified one medieval remains as a hybrid between two common oak species (Quercus robur and Q. petraea) and the other two remains as Q. robur. We found that diversity at the global genome level had not changed over time. However, exploratory analyses suggested that a reduction of diversity took place at different time periods. Finally, we determined the timing of leaf unfolding for ancient trees for the first time. The study extends the application of ancient wood beyond the classical proxies of dendroclimatology, dendrochronology, dendroarchaeology and dendroecology, thereby enhancing resolution of inferences on the responses of forest ecosystems to past environmental changes, epidemics and silvicultural practices.

Conservation and diversity of the pollen microbiome of Pan-American maize using PacBio and MiSeq

Pollen is a vector for diversification, fitness-selection, and transmission of plant genetic material. The extent to which the pollen microbiome may contribute to host diversification is largely unknown, because pollen microbiome diversity within a plant species has not been reported, and studies have been limited to conventional short-read 16S rRNA gene sequencing (e.g., V4-MiSeq) which suffers from poor taxonomic resolution. Here we report the pollen microbiomes of 16 primitive and traditional accessions of maize (corn) selected by indigenous peoples across the Americas, along with the modern U.S. inbred B73. The maize pollen microbiome has not previously been reported. The pollen microbiomes were identified using full-length (FL) 16S rRNA gene PacBio SMRT sequencing compared to V4-MiSeq. The Pan-American maize pollen microbiome encompasses 765 taxa spanning 39 genera and 46 species, including known plant growth promoters, insect-obligates, plant pathogens, nitrogen-fixers and biocontrol agents. Eleven genera and 13 species composed the core microbiome. Of 765 taxa, 63% belonged to only four genera: 28% were Pantoea, 15% were Lactococcus, 11% were Pseudomonas, and 10% were Erwinia. Interestingly, of the 215 Pantoea taxa, 180 belonged to a single species, P. ananatis. Surprisingly, the diversity within P. ananatis ranged nearly 10-fold amongst the maize accessions analyzed (those with ≥3 replicates), despite being grown in a common field. The highest diversity within P. ananatis occurred in accessions that originated near the center of diversity of domesticated maize, with reduced diversity associated with the north-south migration of maize. This sub-species diversity was revealed by FL-PacBio but missed by V4-MiSeq. V4-MiSeq also mis-identified some dominant genera captured by FL-PacBio. The study, though limited to a single season and common field, provides initial evidence that pollen microbiomes reflect evolutionary and migratory relationships of their host plants.

Temporal collections to study invasion biology

Biological invasions represent an extraordinary opportunity to study evolution. This is because accidental or deliberate species introductions have taken place for centuries across large geographical scales, frequently prompting rapid evolutionary transitions in invasive populations. Until recently, however, the utility of invasions as evolutionary experiments has been hampered by limited information on the makeup of populations that were part of earlier invasion stages. Now, developments in ancient and historical DNA technologies, as well as the quickening pace of digitization for millions of specimens that are housed in herbaria and museums globally, promise to help overcome this obstacle. In this review, we first introduce the types of temporal data that can be used to study invasions, highlighting the timescale captured by each approach and their respective limitations. We then discuss how ancient and historical specimens as well as data available from prior invasion studies can be used to answer questions on mechanisms of (mal)adaptation, rates of evolution, or community-level changes during invasions. By bridging the gap between contemporary and historical invasive populations, temporal data can help us connect pattern to process in invasion science. These data will become increasingly important if invasions are to achieve their full potential as experiments of evolution in nature.

Two teosintes made modern maize

The origins of maize were the topic of vigorous debate for nearly a century, but neither the current genetic model nor earlier archaeological models account for the totality of available data, and recent work has highlighted the potential contribution of a wild relative, Zea mays ssp. mexicana. Our population genetic analysis reveals that the origin of modern maize can be traced to an admixture between ancient maize and Zea mays ssp. mexicana in the highlands of Mexico some 4000 years after domestication began. We show that variation in admixture is a key component of maize diversity, both at individual loci and for additive genetic variation underlying agronomic traits. Our results clarify the origin of modern maize and raise new questions about the anthropogenic mechanisms underlying dispersal throughout the Americas.

Genetic analysis of phenotypic plasticity identifies BBX6 as the candidate gene for maize adaptation to temperate regions

Introduction

Climate changes pose a significant threat to crop adaptation and production. Dissecting the genetic basis of phenotypic plasticity and uncovering the responsiveness of regulatory genes to environmental factors can significantly contribute to the improvement of climate- resilience in crops.

Methods

We established a BC1F3:4 population using the elite inbred lines Zheng58 and PH4CV and evaluated plant height (PH) across four environments characterized by substantial variations in environmental factors. Then, we quantified the correlation between the environmental mean of PH (the mean performance in each environment) and the environmental parameters within a specific growth window. Furthermore, we performed GWAS analysis of phenotypic plasticity, and identified QTLs and candidate gene that respond to key environment index. After that, we constructed the coexpression network involving the candidate gene, and performed selective sweep analysis of the candidate gene.

Results

We found that the environmental parameters demonstrated substantial variation across the environments, and genotype by environment interaction contributed to the variations of PH. Then, we identified PTT(35-48) (PTT is the abbreviation for photothermal units), the mean PTT from 35 to 48 days after planting, as the pivotal environmental index that closely correlated with environmental mean of PH. Leveraging the slopes of the response of PH to both the environmental mean and PTT(35-48), we successfully pinpointed QTLs for phenotypic plasticity on chromosomes 1 and 2. Notably, the PH4CV genotypes at these two QTLs exhibited positive contributions to phenotypic plasticity. Furthermore, our analysis demonstrated a direct correlation between the additive effects of each QTL and PTT(35-48). By analyzing transcriptome data of the parental lines in two environments, we found that the 1009 genes responding to PTT(35-48) were enriched in the biological processes related to environmental sensitivity. BBX6 was the prime candidate gene among the 13 genes in the two QTL regions. The coexpression network of BBX6 contained other genes related to flowering time and photoperiod sensitivity. Our investigation, including selective sweep analysis and genetic differentiation analysis, suggested that BBX6 underwent selection during maize domestication.

Discussion

Th is research substantially advances our understanding of critical environmental factors influencing maize adaptation while simultaneously provides an invaluable gene resource for the development of climate-resilient maize hybrid varieties.

Wealth inequality in the prehispanic northern US Southwest: from Malthus to Tyche

Persistent differences in wealth and power among prehispanic Pueblo societies are visible from the late AD 800s through the late 1200s, after which large portions of the northern US Southwest were depopulated. In this paper we measure these differences in wealth using Gini coefficients based on house size, and show that high Ginis (large wealth differences) are positively related to persistence in settlements and inversely related to an annual measure of the size of the unoccupied dry-farming niche. We argue that wealth inequality in this record is due first to processes inherent in village life which have internally different distributions of the most productive maize fields, exacerbated by the dynamics of systems of balanced reciprocity; and second to decreasing ability to escape village life owing to shrinking availability of unoccupied places within the maize dry-farming niche as villages get enmeshed in regional systems of tribute or taxation. We embed this analytical reconstruction in the model of an 'Abrupt imposition of Malthusian equilibrium in a natural-fertility, agrarian society' proposed by Puleston et al. (Puleston C, Tuljapurkar S, Winterhalder B. 2014 PLoS ONE 9, e87541 (doi:10.1371/journal.pone.0087541)), but show that the transition to Malthusian dynamics in this area is not abrupt but extends over centuries This article is part of the theme issue 'Evolutionary ecology of inequality'.

Trans-Holocene Bayesian chronology for tree and field crop use from El Gigante rockshelter, Honduras

El Gigante rockshelter in western Honduras provides a deeply stratified archaeological record of human-environment interaction spanning the entirety of the Holocene. Botanical materials are remarkably well preserved and include important tree (e.g., ciruela (Spondias), avocado (Persea americana)) and field (maize (Zea mays), beans (Phaseolus), and squash (Cucurbita)) crops. Here we provide a major update to the chronology of tree and field crop use evident in the sequence. We report 375 radiocarbon dates, a majority of which are for short-lived botanical macrofossils (e.g., maize cobs, avocado seeds, or rinds). Radiocarbon dates were used in combination with stratigraphic details to establish a Bayesian chronology for ~9,800 identified botanical samples spanning the last 11,000 years. We estimate that at least 16 discrete intervals of use occurred during this time, separated by gaps of ~100-2,000 years. The longest hiatus in rockshelter occupation was between ~6,400 and 4,400 years ago and the deposition of botanical remains peaked at ~2,000 calendar years before present (cal BP). Tree fruits and squash appeared early in the occupational sequence (~11,000 cal BP) with most other field crops appearing later in time (e.g., maize at ~4,400 cal BP; beans at ~2,200 cal BP). The early focus on tree fruits and squash is consistent with early coevolutionary partnering with humans as seed dispersers in the wake of megafaunal extinction in Mesoamerica. Tree crops predominated through much of the Holocene, and there was an overall shift to field crops after 4,000 cal BP that was largely driven by increased reliance on maize farming.

Genotyping of DNA pools identifies untapped landraces and genomic regions to develop next-generation varieties

Landraces, that is, traditional varieties, have a large diversity that is underexploited in modern breeding. A novel DNA pooling strategy was implemented to identify promising landraces and genomic regions to enlarge the genetic diversity of modern varieties. As proof of concept, DNA pools from 156 American and European maize landraces representing 2340 individuals were genotyped with an SNP array to assess their genome-wide diversity. They were compared to elite cultivars produced across the 20th century, represented by 327 inbred lines. Detection of selective footprints between landraces of different geographic origin identified genes involved in environmental adaptation (flowering times, growth) and tolerance to abiotic and biotic stress (drought, cold, salinity). Promising landraces were identified by developing two novel indicators that estimate their contribution to the genome of inbred lines: (i) a modified Roger's distance standardized by gene diversity and (ii) the assignation of lines to landraces using supervised analysis. It showed that most landraces do not have closely related lines and that only 10 landraces, including famous landraces as Reid's Yellow Dent, Lancaster Surecrop and Lacaune, cumulated half of the total contribution to inbred lines. Comparison of ancestral lines directly derived from landraces with lines from more advanced breeding cycles showed a decrease in the number of landraces with a large contribution. New inbred lines derived from landraces with limited contributions enriched more the haplotype diversity of reference inbred lines than those with a high contribution. Our approach opens an avenue for the identification of promising landraces for pre-breeding.

Heritage genetics for adaptation to marginal soils in barley

Future crops need to be sustainable in the face of climate change. Modern barley varieties have been bred for high productivity and quality; however, they have suffered considerable genetic erosion, losing crucial genetic diversity. This renders modern cultivars vulnerable to climate change and stressful environments. We highlight the potential to tailor crops to a specific environment by utilising diversity inherent in an adapted landrace population. Tapping into natural biodiversity, while incorporating information about local environmental and climatic conditions, allows targeting of key traits and genotypes, enabling crop production in marginal soils. We outline future directions for the utilisation of genetic resources maintained in landrace collections to support sustainable agriculture through germplasm development via the use of genomics technologies and big data.

A role for heritable transcriptomic variation in maize adaptation to temperate environments

Background

Transcription bridges genetic information and phenotypes. Here, we evaluated how changes in transcriptional regulation enable maize (Zea mays), a crop originally domesticated in the tropics, to adapt to temperate environments.

Result

We generated 572 unique RNA-seq datasets from the roots of 340 maize genotypes. Genes involved in core processes such as cell division, chromosome organization and cytoskeleton organization showed lower heritability of gene expression, while genes involved in anti-oxidation activity exhibited higher expression heritability. An expression genome-wide association study (eGWAS) identified 19,602 expression quantitative trait loci (eQTLs) associated with the expression of 11,444 genes. A GWAS for alternative splicing identified 49,897 splicing QTLs (sQTLs) for 7614 genes. Genes harboring both cis-eQTLs and cis-sQTLs in linkage disequilibrium were disproportionately likely to encode transcription factors or were annotated as responding to one or more stresses. Independent component analysis of gene expression data identified loci regulating co-expression modules involved in oxidation reduction, response to water deprivation, plastid biogenesis, protein biogenesis, and plant-pathogen interaction. Several genes involved in cell proliferation, flower development, DNA replication, and gene silencing showed lower gene expression variation explained by genetic factors between temperate and tropical maize lines. A GWAS of 27 previously published phenotypes identified several candidate genes overlapping with genomic intervals showing signatures of selection during adaptation to temperate environments.

Conclusion

Our results illustrate how maize transcriptional regulatory networks enable changes in transcriptional regulation to adapt to temperate regions.

Supplementary information

The online version contains supplementary material available at 10.1186/s13059-023-02891-3.

Exploring the sorghum race level diversity utilizing 272 sorghum accessions genomic resources

Due to evolutionary divergence, sorghum race populations exhibit significant genetic and morphological variation. A k-mer-based sorghum race sequence comparison identified the conserved k-mers of all 272 accessions from sorghum and the race-specific genetic signatures identified the gene variability in 10,321 genes (PAVs). To understand sorghum race structure, diversity and domestication, a deep learning-based variant calling approach was employed in a set of genotypic data derived from a diverse panel of 272 sorghum accessions. The data resulted in 1.7 million high-quality genome-wide SNPs and identified selective signature (both positive and negative) regions through a genome-wide scan with different (iHS and XP-EHH) statistical methods. We discovered 2,370 genes associated with selection signatures including 179 selective sweep regions distributed over 10 chromosomes. Co-localization of these regions undergoing selective pressure with previously reported QTLs and genes revealed that the signatures of selection could be related to the domestication of important agronomic traits such as biomass and plant height. The developed k-mer signatures will be useful in the future to identify the sorghum race and for trait and SNP markers for assisting in plant breeding programs.

A common resequencing-based genetic marker data set for global maize diversity

Maize (Zea mays ssp. mays) populations exhibit vast ranges of genetic and phenotypic diversity. As sequencing costs have declined, an increasing number of projects have sought to measure genetic differences between and within maize populations using whole-genome resequencing strategies, identifying millions of segregating single-nucleotide polymorphisms (SNPs) and insertions/deletions (InDels). Unlike older genotyping strategies like microarrays and genotyping by sequencing, resequencing should, in principle, frequently identify and score common genetic variants. However, in practice, different projects frequently employ different analytical pipelines, often employ different reference genome assemblies and consistently filter for minor allele frequency within the study population. This constrains the potential to reuse and remix data on genetic diversity generated from different projects to address new biological questions in new ways. Here, we employ resequencing data from 1276 previously published maize samples and 239 newly resequenced maize samples to generate a single unified marker set of approximately 366 million segregating variants and approximately 46 million high-confidence variants scored across crop wild relatives, landraces as well as tropical and temperate lines from different breeding eras. We demonstrate that the new variant set provides increased power to identify known causal flowering-time genes using previously published trait data sets, as well as the potential to track changes in the frequency of functionally distinct alleles across the global distribution of modern maize.

The evening complex promotes maize flowering and adaptation to temperate regions

Maize (Zea mays) originated in southern Mexico and has spread over a wide latitudinal range. Maize expansion from tropical to temperate regions has necessitated a reduction of its photoperiod sensitivity. In this study, we cloned a quantitative trait locus (QTL) regulating flowering time in maize and show that the maize ortholog of Arabidopsis thaliana EARLY FLOWERING3, ZmELF3.1, is the causal locus. We demonstrate that ZmELF3.1 and ZmELF3.2 proteins can physically interact with ZmELF4.1/4.2 and ZmLUX1/2, to form evening complex(es; ECs) in the maize circadian clock. Loss-of-function mutants for ZmELF3.1/3.2 and ZmLUX1/2 exhibited delayed flowering under long-day and short-day conditions. We show that EC directly represses the expression of several flowering suppressor genes, such as the CONSTANS, CONSTANS-LIKE, TOC1 (CCT) genes ZmCCT9 and ZmCCT10, ZmCONSTANS-LIKE 3, and the PSEUDORESPONSE REGULATOR (PRR) genes ZmPRR37a and ZmPRR73, thus alleviating their inhibition, allowing florigen gene expression and promoting flowering. Further, we identify two closely linked retrotransposons located in the ZmELF3.1 promoter that regulate the expression levels of ZmELF3.1 and may have been positively selected during postdomestication spread of maize from tropical to temperate regions during the pre-Columbian era. These findings provide insights into circadian clock-mediated regulation of photoperiodic flowering in maize and new targets of genetic improvement for breeding.

The herbarium of the future

The ~400 million specimens deposited across ~3000 herbaria are essential for: (i) understanding where plants have lived in the past, (ii) forecasting where they may live in the future, and (iii) delineating their conservation status. An open access 'global metaherbarium' is emerging as these specimens are digitized, mobilized, and interlinked online. This virtual biodiversity resource is attracting new users who are accelerating traditional applications of herbaria and generating basic and applied scientific innovations, including e-monographs and floras produced by diverse, interdisciplinary, and inclusive teams; robust machine-learning algorithms for species identification and phenotyping; collection and synthesis of ecological trait data at large spatiotemporal and phylogenetic scales; and exhibitions and installations that convey the beauty of plants and the value of herbaria in addressing broader societal issues.

Genome sequencing reveals evidence of adaptive variation in the genus Zea

Maize is a globally valuable commodity and one of the most extensively studied genetic model organisms. However, we know surprisingly little about the extent and potential utility of the genetic variation found in wild relatives of maize. Here, we characterize a high-density genomic variation map from 744 genomes encompassing maize and all wild taxa of the genus Zea, identifying over 70 million single-nucleotide polymorphisms. The variation map reveals evidence of selection within taxa displaying novel adaptations. We focus on adaptive alleles in highland teosinte and temperate maize, highlighting the key role of flowering-time-related pathways in their adaptation. To show the utility of variants in these data, we generate mutant alleles for two flowering-time candidate genes. This work provides an extensive sampling of the genetic diversity of Zea, resolving questions on evolution and identifying adaptive variants for direct use in modern breeding.

Cultivation has selected for a wider niche and large range shifts in maize

Background

Maize (Zea mays L.) is a staple crop cultivated on a global scale. However, its ability to feed the rapidly growing human population may be impaired by climate change, especially if it has low climatic niche and range lability. One important question requiring clarification is therefore whether maize shows high niche and range lability.

Methods

We used the COUE scheme (a unified terminology representing niche centroid shift, overlap, unfilling and expansion) and species distribution models to study the niche and range changes between maize and its wild progenitors using occurrence records of maize, lowland teosinte (Zea mays ssp. parviglumis) and highland teosinte (Zea mays ssp. mexicana), respectively, as well as explore the mechanisms underlying the niche and range changes.

Results

In contrast to maize in Mexico, maize did not conserve its niche inherited from lowland and highland teosinte at the global scale. The niche breadth of maize at the global scale was wider than that of its wild progenitors (ca. 5.21 and 3.53 times wider compared with lowland and highland teosinte, respectively). Compared with its wild progenitors, maize at global scale can survive in regions with colder, wetter climatic conditions, as well as with wider ranges of climatic variables (ca. 4.51 and 2.40 times wider compared with lowland and highland teosinte, respectively). The niche changes of maize were largely driven by human introduction and cultivation, which have exposed maize to climatic conditions different from those experienced by its wild progenitors. Small changes in niche breadth had large effects on the magnitude of range shifts; changes in niche breadth thus merit increased attention.

Discussion

Our results demonstrate that maize shows wide climatic niche and range lability, and this substantially expanded its realized niche and potential range. Our findings also suggest that niche and range shifts probably triggered by natural and artificial selection in cultivation may enable maize to become a global staple crop to feed the growing population and adapting to changing climatic conditions. Future analyses are needed to determine the limits of the novel conditions that maize can tolerate, especially relative to projected climate change.

A pan-Zea genome map for enhancing maize improvement

Background

Maize (Zea mays L.) is at the vanguard facing the upcoming breeding challenges. However, both a super pan-genome for the Zea genus and a comprehensive genetic variation map for maize breeding are still lacking.

Results

Here, we construct an approximately 6.71-Gb pan-Zea genome that contains around 4.57-Gb non-B73 reference sequences from fragmented de novo assemblies of 721 pan-Zea individuals. We annotate a total of 58,944 pan-Zea genes and find around 44.34% of them are dispensable in the pan-Zea population. Moreover, 255,821 common structural variations are identified and genotyped in a maize association mapping panel. Further analyses reveal gene presence/absence variants and their potential roles during domestication of maize. Combining genetic analyses with multi-omics data, we demonstrate how structural variants are associated with complex agronomic traits.

Conclusions

Our results highlight the underexplored role of the pan-Zea genome and structural variations to further understand domestication of maize and explore their potential utilization in crop improvement.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02742-7.

An adaptive teosinte mexicana introgression modulates phosphatidylcholine levels and is associated with maize flowering time

Native Americans domesticated maize (Zea mays ssp. mays) from lowland teosinte parviglumis (Zea mays ssp. parviglumis) in the warm Mexican southwest and brought it to the highlands of Mexico and South America where it was exposed to lower temperatures that imposed strong selection on flowering time. Phospholipids are important metabolites in plant responses to low-temperature and phosphorus availability and have been suggested to influence flowering time. Here, we combined linkage mapping with genome scans to identify High PhosphatidylCholine 1 (HPC1), a gene that encodes a phospholipase A1 enzyme, as a major driver of phospholipid variation in highland maize. Common garden experiments demonstrated strong genotype-by-environment interactions associated with variation at HPC1, with the highland HPC1 allele leading to higher fitness in highlands, possibly by hastening flowering. The highland maize HPC1 variant resulted in impaired function of the encoded protein due to a polymorphism in a highly conserved sequence. A meta-analysis across HPC1 orthologs indicated a strong association between the identity of the amino acid at this position and optimal growth in prokaryotes. Mutagenesis of HPC1 via genome editing validated its role in regulating phospholipid metabolism. Finally, we showed that the highland HPC1 allele entered cultivated maize by introgression from the wild highland teosinte Zea mays ssp. mexicana and has been maintained in maize breeding lines from the Northern United States, Canada, and Europe. Thus, HPC1 introgressed from teosinte mexicana underlies a large metabolic QTL that modulates phosphatidylcholine levels and has an adaptive effect at least in part via induction of early flowering time.

Polygenic score accuracy in ancient samples: Quantifying the effects of allelic turnover

Polygenic scores link the genotypes of ancient individuals to their phenotypes, which are often unobservable, offering a tantalizing opportunity to reconstruct complex trait evolution. In practice, however, interpretation of ancient polygenic scores is subject to numerous assumptions. For one, the genome-wide association (GWA) studies from which polygenic scores are derived, can only estimate effect sizes for loci segregating in contemporary populations. Therefore, a GWA study may not correctly identify all loci relevant to trait variation in the ancient population. In addition, the frequencies of trait-associated loci may have changed in the intervening years. Here, we devise a theoretical framework to quantify the effect of this allelic turnover on the statistical properties of polygenic scores as functions of population genetic dynamics, trait architecture, power to detect significant loci, and the age of the ancient sample. We model the allele frequencies of loci underlying trait variation using the Wright-Fisher diffusion, and employ the spectral representation of its transition density to find analytical expressions for several error metrics, including the expected sample correlation between the polygenic scores of ancient individuals and their true phenotypes, referred to as polygenic score accuracy. Our theory also applies to a two-population scenario and demonstrates that allelic turnover alone may explain a substantial percentage of the reduced accuracy observed in cross-population predictions, akin to those performed in human genetics. Finally, we use simulations to explore the effects of recent directional selection, a bias-inducing process, on the statistics of interest. We find that even in the presence of bias, weak selection induces minimal deviations from our neutral expectations for the decay of polygenic score accuracy. By quantifying the limitations of polygenic scores in an explicit evolutionary context, our work lays the foundation for the development of more sophisticated statistical procedures to analyze both temporally and geographically resolved polygenic scores.

More opportunities more species: Pleistocene differentiation and northward expansion of an evergreen broad-leaved tree species Machilus thunbergii (Lauraceae) in Southeast China

BACKGROUND

The broad continuum between tropical and temperate floras in Eastern Asia (EAS) are thought to be one of the main factors responsible for a prominent species diversity anomaly of temperate plants between EAS and eastern North America (ENS). However, how the broad continuum and niche evolution between tropical and temperate floras in EAS contributes to lineage divergence and species diversity remains largely unknown.

RESULTS

Population genetic structure, demography, and determinants of genetic structure [i.e., isolation-by-distance (IBD), isolation-by-resistance (IBR), and isolation-by-environment (IBE)] of Machilus thunbergii Sieb. et Zucc. (Lauraceae) were evaluated by examining sequence variation of ten low-copy nuclear genes across 43 populations in southeast China. Climatic niche difference and potential distributions across four periods (Current, mid-Holocene, the last glacial maximum, the last interglacial) of two genetic clusters were determined by niche modelling. North and south clusters of populations in M. thunbergii were revealed and their demarcation line corresponds well with the northern boundary of tropical zone in China of Zhu & Wan. The divergence time between the clusters and demographic expansion of M. thunbergii occurred after the mid-Pleistocene climate transition (MPT, 0.8-1.2 Ma). Migration rates between clusters were asymmetrical, being much greater from north to south than the reverse. Significant effects of IBE, but non-significant effects of IBD and IBR on population genetic divergence were detected. The two clusters have different ecological niches and require different temperature regimes.

CONCLUSIONS

The north-south genetic differentiation may be common across the temperate-tropical boundary in southeast China. Divergent selection under different temperature regimes (possibly above and below freezing temperature in winter) could account for this divergence pattern. The broad continuum between tropical and temperate floras in EAS may have provided ample opportunities for tropical plant lineages to acquire freezing tolerance and to colonize the temperate regions during the late-Cenozoic global cooling. Our findings shed deeper insights into the high temperate plant species diversity in EAS.

Discovering the Secrets of Ancient Plants: Recovery of DNA from Museum and Archaeological Plant Specimens

Plant DNA preserved in ancient specimens has recently gained importance as a tool in comparative genomics, allowing the investigation of evolutionary processes in plant genomes through time. However, recovering the genomic information contained in such specimens is challenging owing to the presence of secondary substances that limit DNA retrieval. In this chapter, we provide a DNA extraction protocol optimized for the recovery of DNA from degraded plant materials. The protocol is based on a commercially available DNA extraction kit that does not require handling of hazardous reagents.

Emerging evidence of plant domestication as a landscape-level process

The evidence from ancient crops over the past decade challenges some of our most basic assumptions about the process of domestication. The emergence of crops has been viewed as a technologically progressive process in which single or multiple localized populations adapt to human environments in response to cultivation. By contrast, new genetic and archaeological evidence reveals a slow process that involved large populations over wide areas with unexpectedly sustained cultural connections in deep time. We review evidence that calls for a new landscape framework of crop origins. Evolutionary processes operate across vast distances of landscape and time, and the origins of domesticates are complex. The domestication bottleneck is a redundant concept and the progressive nature of domestication is in doubt.

The arches and spandrels of maize domestication, adaptation, and improvement

People living in the Balsas River basin in southwest México domesticated maize from the bushy grass teosinte. Nine thousand years later, in 2021, Ms. Deb Haaland - a member of the Pueblo of Laguna tribe of New Mexico - wore a dress adorned with a cornstalk when she was sworn in as the Secretary of Interior of the United States of America. This choice of garment highlights the importance of the coevolution of maize and the farmers who, through careful selection over thousands of years, domesticated maize and adapted the physiology and shoot architecture of maize to fit local environments and growth habits. Some traits such as tillering were directly selected on (arches), and others such as tassel size are the by-products (spandrels) of maize evolution. Here, we review current knowledge of the underlying cellular, developmental, physiological, and metabolic processes that were selected by farmers and breeders, which have positioned maize as a top global staple crop.

The genetic architecture and population genomic signatures of glyphosate resistance in Amaranthus tuberculatus

Much of what we know about the genetic basis of herbicide resistance has come from detailed investigations of monogenic adaptation at known target-sites, despite the increasingly recognized importance of polygenic resistance. Little work has been done to characterize the broader genomic basis of herbicide resistance, including the number and distribution of genes involved, their effect sizes, allele frequencies and signatures of selection. In this work, we implemented genome-wide association (GWA) and population genomic approaches to examine the genetic architecture of glyphosate (Round-up) resistance in the problematic agricultural weed Amaranthus tuberculatus. A GWA was able to correctly identify the known target-gene but statistically controlling for two causal target-site mechanisms revealed an additional 250 genes across all 16 chromosomes associated with non-target-site resistance (NTSR). The encoded proteins had functions that have been linked to NTSR, the most significant of which is response to chemicals, but also showed pleiotropic roles in reproduction and growth. Compared to an empirical null that accounts for complex population structure, the architecture of NTSR was enriched for large effect sizes and low allele frequencies, suggesting the role of pleiotropic constraints on its evolution. The enrichment of rare alleles also suggested that the genetic architecture of NTSR may be population-specific and heterogeneous across the range. Despite their rarity, we found signals of recent positive selection on NTSR-alleles by both window- and haplotype-based statistics, and an enrichment of amino acid changing variants. In our samples, genome-wide single nucleotide polymorphisms explain a comparable amount of the total variation in glyphosate resistance to monogenic mechanisms, even in a collection of individuals where 80% of resistant individuals have large-effect TSR mutations, indicating an underappreciated polygenic contribution to the evolution of herbicide resistance in weed populations.

Reference genomes of the two cultivated jute species

Cultivated jute, which comprises the two species Corchorus capsularis and C. olitorius, is the second most important natural fibre source after cotton. Here we describe chromosome-level assemblies of the genomes of both cultivated species. The C. capsularis and C. olitorius assemblies are each comprised of seven pseudo-chromosomes, with the C. capsularis assembly consisting of 336 Mb with 25,874 genes and the C. olitorius assembly containing 361 Mb with 28 479 genes. Although the two Corchorus genomes exhibit collinearity, the genome of C. olitorius contains 25 Mb of additional sequences than that of C. capsularis with 13 putative inversions, which might give a hint to the difference of phenotypic variants between the two cultivated jute species. Analysis of gene expression in isolated fibre tissues reveals candidate genes involved in fibre development. Our analysis of the population structures of 242 cultivars from C. capsularis and 57 cultivars from C. olitorius by whole-genome resequencing resulted in post-domestication bottlenecks occurred ~2000 years ago in these species. We identified hundreds of putative significant marker-trait associations (MTAs) controlling fibre fineness, cellulose content and lignin content of fibre by integrating data from genome-wide association studies (GWAS) with data from analyses of selective sweeps due to natural and artificial selection in these two jute species. Among them, we further validated that CcCOBRA1 and CcC4H1 regulate fibre quality in transgenic plants via improving the biosynthesis of the secondary cell wall. Our results yielded important new resources for functional genomics research and genetic improvement in jute and allied fibre crops.

Phenotypic Plasticity Contributes to Maize Adaptation and Heterosis

Plant phenotypic plasticity describes altered phenotypic performance of an individual when grown in different environments. Exploring genetic architecture underlying plant plasticity variation may help mitigate the detrimental effects of a rapidly changing climate on agriculture, but little research has been done in this area to date. In the present study, we established a population of 976 maize F₁ hybrids by crossing 488 diverse inbred lines with two elite testers. Genome-wide association study identified hundreds of quantitative trait loci associated with phenotypic plasticity variation across diverse F₁ hybrids, the majority of which contributed very little variance, in accordance with the polygenic nature of these traits. We identified several quantitative trait locus regions that may have been selected during the tropical-temperate adaptation process. We also observed heterosis in terms of phenotypic plasticity, in addition to the traditional genetic value differences measured between hybrid and inbred lines, and the pattern of which was affected by genetic background. Our results demonstrate a landscape of phenotypic plasticity in maize, which will aid in the understanding of its genetic architecture, its contribution to adaptation and heterosis, and how it may be exploited for future maize breeding in a rapidly changing environment.

Local adaptation contributes to gene expression divergence in maize

Gene expression links genotypes to phenotypes, so identifying genes whose expression is shaped by selection will be important for understanding the traits and processes underlying local adaptation. However, detecting local adaptation for gene expression will require distinguishing between divergence due to selection and divergence due to genetic drift. Here, we adapt a QST−FST framework to detect local adaptation for transcriptome-wide gene expression levels in a population of diverse maize genotypes. We compare the number and types of selected genes across a wide range of maize populations and tissues, as well as selection on cold-response genes, drought-response genes, and coexpression clusters. We identify a number of genes whose expression levels are consistent with local adaptation and show that genes involved in stress response show enrichment for selection. Due to its history of intense selective breeding and domestication, maize evolution has long been of interest to researchers, and our study provides insight into the genes and processes important for in local adaptation of maize.

Molecular clocks and archaeogenomics of a Late Period Egyptian date palm leaf reveal introgression from wild relatives and add timestamps on the domestication

The date palm, Phoenix dactylifera, has been a cornerstone of Middle Eastern and North African agriculture for millennia. It was first domesticated in the Persian Gulf, and its evolution appears to have been influenced by gene flow from two wild relatives, P. theophrasti, currently restricted to Crete and Turkey, and P. sylvestris, widespread from Bangladesh to the West Himalayas. Genomes of ancient date palm seeds show that gene flow from P. theophrasti to P. dactylifera may have occurred by ∼2,200 years ago, but traces of P. sylvestris could not be detected. We here integrate archeogenomics of a ∼2,100-year-old P. dactylifera leaf from Saqqara (Egypt), molecular-clock dating, and coalescence approaches with population genomic tests, to probe the hybridization between the date palm and its two closest relatives and provide minimum and maximum timestamps for its reticulated evolution. The Saqqara date palm shares a close genetic affinity with North African date palm populations, and we find clear genomic admixture from both P. theophrasti, and P. sylvestris, indicating that both had contributed to the date palm genome by 2,100 years ago. Molecular-clocks placed the divergence of P. theophrasti from P. dactylifera/P. sylvestris and that of P. dactylifera from P. sylvestris in the Upper Miocene, but strongly supported, conflicting topologies point to older gene flow between P. theophrasti and P. dactylifera, and P. sylvestris and P. dactylifera. Our work highlights the ancient hybrid origin of the date palms, and prompts the investigation of the functional significance of genetic material introgressed from both close relatives, which in turn could prove useful for modern date palm breeding.

Archaeology and agriculture: plants, people, and past land-use

As a specialised branch of archaeology requiring specific field and laboratory methodologies, the contributions of archaeobotany have often been overlooked by the ecological research community. Developments in the fields of botany, chemistry, and ancient DNA analyses have greatly increased the potential for archaeobotany to contribute to topical questions relating to the Anthropocene and landscape transformations. We review the role of archaeobotany in identifying and describing past arable land use. Analytical techniques are illustrated with examples at both local and regional scales, demonstrating how archaeobotany can provide unique details of the wide array of past subsistence and land-use strategies. These data and their potential should be better recognised as important information that could underpin models seeking to evaluate or predict the effects of socioenvironmental interactions.

Meta-Analysis Identifies Pleiotropic Loci Controlling Phenotypic Trade-offs in Sorghum

Community association populations are composed of phenotypically and genetically diverse accessions. Once these populations are genotyped, the resulting marker data can be reused by different groups investigating the genetic basis of different traits. Because the same genotypes are observed and scored for a wide range of traits in different environments, these populations represent a unique resource to investigate pleiotropy. Here we assembled a set of 234 separate trait datasets for the Sorghum Association Panel, a group of 406 sorghum genotypes widely employed by the sorghum genetics community. Comparison of genome wide association studies conducted with two independently generated marker sets for this population demonstrate that existing genetic marker sets do not saturate the genome and likely capture only 35-43% of potentially detectable loci controlling variation for traits scored in this population. While limited evidence for pleiotropy was apparent in cross-GWAS comparisons, a multivariate adaptive shrinkage approach recovered both known pleiotropic effects of existing loci and new pleiotropic effects, particularly significant impacts of known dwarfing genes on root architecture. In addition, we identified new loci with pleiotropic effects consistent with known trade-offs in sorghum development. These results demonstrate the potential for mining existing trait datasets from widely used community association populations to enable new discoveries from existing trait datasets as new, denser genetic marker datasets are generated for existing community association populations.

Joint analysis of days to flowering reveals independent temperate adaptations in maize

Domesticates are an excellent model for understanding biological consequences of rapid climate change. Maize (Zea mays ssp. mays) was domesticated from a tropical grass yet is widespread across temperate regions today. We investigate the biological basis of temperate adaptation in diverse structured nested association mapping (NAM) populations from China, Europe (Dent and Flint) and the United States as well as in the Ames inbred diversity panel, using days to flowering as a proxy. Using cross-population prediction, where high prediction accuracy derives from overall genomic relatedness, shared genetic architecture, and sufficient diversity in the training population, we identify patterns in predictive ability across the five populations. To identify the source of temperate adapted alleles in these populations, we predict top associated genome-wide association study (GWAS) identified loci in a Random Forest Classifier using independent temperate-tropical North American populations based on lines selected from Hapmap3 as predictors. We find that North American populations are well predicted (AUC equals 0.89 and 0.85 for Ames and USNAM, respectively), European populations somewhat well predicted (AUC equals 0.59 and 0.67 for the Dent and Flint panels, respectively) and that the Chinese population is not predicted well at all (AUC is 0.47), suggesting an independent adaptation process for early flowering in China. Multiple adaptations for the complex trait days to flowering in maize provide hope for similar natural systems under climate change.

Eleven biosynthetic genes explain the majority of natural variation in carotenoid levels in maize grain

Vitamin A deficiency remains prevalent in parts of Asia, Latin America, and sub-Saharan Africa where maize (Zea mays) is a food staple. Extensive natural variation exists for carotenoids in maize grain. Here, to understand its genetic basis, we conducted a joint linkage and genome-wide association study of the US maize nested association mapping panel. Eleven of the 44 detected quantitative trait loci (QTL) were resolved to individual genes. Six of these were correlated expression and effect QTL (ceeQTL), showing strong correlations between RNA-seq expression abundances and QTL allelic effect estimates across six stages of grain development. These six ceeQTL also had the largest percentage of phenotypic variance explained, and in major part comprised the three to five loci capturing the bulk of genetic variation for each trait. Most of these ceeQTL had strongly correlated QTL allelic effect estimates across multiple traits. These findings provide an in-depth genome-level understanding of the genetic and molecular control of carotenoids in plants. In addition, these findings provide a roadmap to accelerate breeding for provitamin A and other priority carotenoid traits in maize grain that should be readily extendable to other cereals.

Maize metabolome and proteome responses to controlled cold stress partly mimic early-sowing effects in the field and differ from those of Arabidopsis

In Northern Europe, sowing maize one-month earlier than current agricultural practices may lead to moderate chilling damage. However, studies of the metabolic responses to low, non-freezing, temperatures remain scarce. Here, genetically-diverse maize hybrids (Zea mays, dent inbred lines crossed with a flint inbred line) were cultivated in a growth chamber at optimal temperature and then three decreasing temperatures for 2 days each, as well as in the field. Leaf metabolomic and proteomic profiles were determined. In the growth chamber, 50% of metabolites and 18% of proteins changed between 20 and 16°C. These maize responses, partly differing from those of Arabidopsis to short-term chilling, were mapped on genome-wide metabolic maps. Several metabolites and proteins showed similar variation for all temperature decreases: seven MS-based metabolite signatures and two proteins involved in photosynthesis decreased continuously. Several increasing metabolites or proteins in the growth-chamber chilling conditions showed similar trends in the early-sowing field experiment, including trans-aconitate, three hydroxycinnamate derivatives, a benzoxazinoid, a sucrose synthase, lethal leaf-spot 1 protein, an allene oxide synthase, several glutathione transferases and peroxidases. Hybrid groups based on field biomass were used to search for the metabolite or protein responses differentiating them in growth-chamber conditions, which could be of interest for breeding.

Evolution: Replicated mutation of COL2 contributed long-day flowering in common bean

The ability to flower without strict daylength constraints has helped spread cultivation of crop plants to new locations. The generation of daylength-insensitive common bean accessions in central and South America involved the repeated selection of mutant alleles for a key transcription factor that suppresses long-day flowering.

Genotyping by Sequencing Reveals Genetic Relatedness of Southwestern U.S. Blue Maize Landraces

Maize has played a key role in the sustenance and cultural traditions of the inhabitants of the southwestern USA for many centuries. Blue maize is an important component of the diverse landraces still cultivated in the region but the degree to which they are related is unknown. This research was designed to ascertain the genotypic, morphological, and phenotypic diversity of six representative southwestern blue maize landraces. Their genotypic diversity was examined using tunable genotyping-by-sequencing (tGBS™). A total of 81,038 high quality SNPs were identified and obtained through tGBS. A total of 45 morphological and biochemical traits were evaluated at two locations in New Mexico. The varieties Los Lunas High and Flor del Rio were genetically less related with other southwestern landraces whereas diffusion between Navajo Blue, Hopi Blue, Yoeme Blue, and Taos Blue demonstrated that these landraces were genetically related. Phenotypic variability was highest for kernel traits and least for plant traits. Plant, ear, and kernel traits were fairly consistent within and across locations. Principal component analysis and tGBS showed that Corn Belt variety ‘Ohio Blue’ was distinctly different from southwestern landraces. Genotypic analysis displayed that southwestern landraces are genetically closely related, but selection has resulted in differing phenotypes. This study has provided additional insight into the genetic relatedness of southwestern blue maize landraces.

Genome assembly and population genomic analysis provide insights into the evolution of modern sweet corn

Sweet corn is one of the most important vegetables in the United States and Canada. Here, we present a de novo assembly of a sweet corn inbred line Ia453 with the mutated shrunken2-reference allele (Ia453-sh2). This mutation accumulates more sugar and is present in most commercial hybrids developed for the processing and fresh markets. The ten pseudochromosomes cover 92% of the total assembly and 99% of the estimated genome size, with a scaffold N50 of 222.2 Mb. This reference genome completely assembles the large structural variation that created the mutant sh2-R allele. Furthermore, comparative genomics analysis with six field corn genomes highlights differences in single-nucleotide polymorphisms, structural variations, and transposon composition. Phylogenetic analysis of 5,381 diverse maize and teosinte accessions reveals genetic relationships between sweet corn and other types of maize. Our results show evidence for a common origin in northern Mexico for modern sweet corn in the U.S. Finally, population genomic analysis identifies regions of the genome under selection and candidate genes associated with sweet corn traits, such as early flowering, endosperm composition, plant and tassel architecture, and kernel row number. Our study provides a high-quality reference-genome sequence to facilitate comparative genomics, functional studies, and genomic-assisted breeding for sweet corn.

3D morphology of nematode encapsulation in snail shells, revealed by micro-CT imaging

Many parasites and hosts are embroiled in an on-going arms race that affects the evolution of each participant. One such battle is between parasitic nematodes and terrestrial gastropods which have co-evolved for 90–130 MY. Recently, snails have been shown to encase and kill invading nematodes using their shell as a defence mechanism. However, there is remarkably little known about this process in terms of understanding where, when and how nematodes are fixed within the shell. Also there has never been any attempt to observe this process using methods other than light microscopy. Therefore, we used micro CT scanning of a Cepaea nemoralis shell (a common host for nematodes) to 3D visualise encased nematode parasites and quantify morphological parameters. By taking this approach future studies could use micro CT scanning of fossil shells in conchology collections to understand nematode/snail co-evolution.

Archaeological Central American maize genomes suggest ancient gene flow from South America

Maize is a global food staple with great economic and cultural importance. Archaeogenomic studies have revealed a process of protracted maize domestication and multiple waves of human-mediated dispersal in the Americas. Maize first arrived in South America as a partial domesticate, where the domestication syndrome became independently fixed and improved varieties developed away from the influence of wild gene flow. We demonstrate that hybrids of some of these improved varieties were likely reintroduced back to Central America. We hypothesize that this backflow of South American genetic material may have contributed to the development of a more productive staple, which was related to the growth and aggregation of human populations, and the formation of more complex social and political structures regionally.

Maize (Zea mays ssp. mays) domestication began in southwestern Mexico ∼9,000 calendar years before present (cal. BP) and humans dispersed this important grain to South America by at least 7,000 cal. BP as a partial domesticate. South America served as a secondary improvement center where the domestication syndrome became fixed and new lineages emerged in parallel with similar processes in Mesoamerica. Later, Indigenous cultivators carried a second major wave of maize southward from Mesoamerica, but it has been unclear until now whether the deeply divergent maize lineages underwent any subsequent gene flow between these regions. Here we report ancient maize genomes (2,300–1,900 cal. BP) from El Gigante rock shelter, Honduras, that are closely related to ancient and modern maize from South America. Our findings suggest that the second wave of maize brought into South America hybridized with long-established landraces from the first wave, and that some of the resulting newly admixed lineages were then reintroduced to Central America. Direct radiocarbon dates and cob morphological data from the rock shelter suggest that more productive maize varieties developed between 4,300 and 2,500 cal. BP. We hypothesize that the influx of maize from South America into Central America may have been an important source of genetic diversity as maize was becoming a staple grain in Central and Mesoamerica.

Isolation, Library Preparation, and Bioinformatic Analysis of Historical and Ancient Plant DNA

The ability to sequence DNA retrieved from ancient and historical material plays a crucial role in reinforcing evolutionary and anthropological inference. While the focus of the field is largely on analyzing DNA from ancient hominids and other animals, we have also learned from plant ancient DNA (aDNA), in particular, about human farming practices, crop domestication, environment management, species invasion, and adaptation to various environmental conditions. In the following protocols, we outline best practices for plant aDNA isolation, preparation for sequencing, bioinformatic processing, and authentication. We describe the process all the way from processing of archaeological or historical plant material to characterizing and authenticating sequencing reads. In alternative protocols, we include modifications to this process that are tailored to strongly degraded DNA. Throughout, we stress the importance of precautionary measures to successfully analyze aDNA. Finally, we discuss the evolution of the archaeogenomics field and the development of new methods, which both shaped this protocol. © 2020 Wiley Periodicals LLC. Basic Protocol 1: Isolation of aDNA Alternate Protocol 1: Isolation of ultra-short DNA (Dabney modification) Support Protocol 1: Preparation of PTB-based mix Support Protocol 2: Preparation of binding buffer Basic Protocol 2: Preparation of genomic libraries Alternate Protocol 2: Preparation of genomic libraries with uracil removal Basic Protocol 3: Bioinformatic processing and authentication of aDNA.

Evolutionary and functional genomics of DNA methylation in maize domestication and improvement

DNA methylation is a ubiquitous chromatin feature, present in 25% of cytosines in the maize genome, but variation and evolution of the methylation landscape during maize domestication remain largely unknown. Here, we leverage whole-genome sequencing (WGS) and whole-genome bisulfite sequencing (WGBS) data on populations of modern maize, landrace, and teosinte (Zea mays ssp. parviglumis) to estimate epimutation rates and selection coefficients. We find weak evidence for direct selection on DNA methylation in any context, but thousands of differentially methylated regions (DMRs) are identified population-wide that are correlated with recent selection. For two trait-associated DMRs, vgt1-DMR and tb1-DMR, HiChIP data indicate that the interactive loops between DMRs and respective downstream genes are present in B73, a modern maize line, but absent in teosinte. Our results enable a better understanding of the evolutionary forces acting on patterns of DNA methylation and suggest a role of methylation variation in adaptive evolution.

Variation and evolution of DNA methylation during maize domestication remain largely unknown. Here, the authors generate genome and methylome sequencing data as well as HiChIP-based interactome data to investigate the adaptive and phenotypic consequences of methylation variations in maize.

Predicting Thermal Adaptation by Looking Into Populations' Genomic Past

Molecular evolution offers an insightful theory to interpret the genomic consequences of thermal adaptation to previous events of climate change beyond range shifts. However, disentangling often mixed footprints of selective and demographic processes from those due to lineage sorting, recombination rate variation, and genomic constrains is not trivial. Therefore, here we condense current and historical population genomic tools to study thermal adaptation and outline key developments (genomic prediction, machine learning) that might assist their utilization for improving forecasts of populations' responses to thermal variation. We start by summarizing how recent thermal-driven selective and demographic responses can be inferred by coalescent methods and in turn how quantitative genetic theory offers suitable multi-trait predictions over a few generations via the breeder's equation. We later assume that enough generations have passed as to display genomic signatures of divergent selection to thermal variation and describe how these footprints can be reconstructed using genome-wide association and selection scans or, alternatively, may be used for forward prediction over multiple generations under an infinitesimal genomic prediction model. Finally, we move deeper in time to comprehend the genomic consequences of thermal shifts at an evolutionary time scale by relying on phylogeographic approaches that allow for reticulate evolution and ecological parapatric speciation, and end by envisioning the potential of modern machine learning techniques to better inform long-term predictions. We conclude that foreseeing future thermal adaptive responses requires bridging the multiple spatial scales of historical and predictive environmental change research under modern cohesive approaches such as genomic prediction and machine learning frameworks.