The genomic signature of crop-wild introgression in maize

A genomic journey across the past, present, and future of South American maize

Maize (Zea mays ssp. mays) diverged from one of its wild relatives, the teosinte Zea mays ssp. parviglumis, in the lowlands of southwest Mexico ~9000 years ago. Following this divergence, maize rapidly expanded throughout the Americas becoming a staple food. This dispersal was accompanied by significant demographic and selective changes, leading to the development of numerous local varieties with a complex evolutionary history that remains incompletely understood. In recent years, genomic advances have challenged traditional models of maize domestication and spread to South America. At least three distinct genetic lineages associated with different migratory waves have been described: ancestral Andean, ancestral Lowland, and Pan-American. Additionally, the significant role of the teosinte Zea mays ssp. mexicana in the evolution of modern maize has been recently uncovered. Genomic studies have shed light on highland adaptation processes, revealing largely independent adaptation events in Meso and South America. As new evidence emerges, the regional complexity underlying maize diversity and the need for comprehensive, multiscale approaches have become evident. In the face of climate change and evolving agricultural landscapes, the conservation of native maize in South America is of growing interest, with genomics serving as an invaluable tool for identifying and preserving the genetic variability of locally adapted germplasm.

The trichome pattern diversity of Cardamine shares genetic mechanisms with Arabidopsis but differs in environmental drivers

Natural variation in trichome pattern (amount and distribution) is prominent among populations of many angiosperms. However, the degree of parallelism in the genetic mechanisms underlying this diversity and its environmental drivers in different species remain unclear. To address these questions, we analyzed the genomic and environmental bases of leaf trichome pattern diversity in Cardamine hirsuta, a relative of Arabidopsis (Arabidopsis thaliana). We characterized 123 wild accessions for their genomic diversity, leaf trichome patterns at different temperatures, and environmental adjustments. Nucleotide diversities and biogeographical distribution models identified two major genetic lineages with distinct demographic and adaptive histories. Additionally, C. hirsuta showed substantial variation in trichome pattern and plasticity to temperature. Trichome amount in C. hirsuta correlated positively with spring precipitation but negatively with temperature, which is opposite to climatic patterns in A. thaliana. Contrastingly, genetic analysis of C. hirsuta glabrous accessions indicated that, like for A. thaliana, glabrousness is caused by null mutations in ChGLABRA1 (ChGL1). Phenotypic genome-wide association studies (GWAS) further identified a ChGL1 haplogroup associated with low trichome density and ChGL1 expression. Therefore, a ChGL1 series of null and partial loss-of-function alleles accounts for the parallel evolution of leaf trichome pattern in C. hirsuta and A. thaliana. Finally, GWAS also detected other candidate genes (e.g. ChETC3, ChCLE17) that might affect trichome pattern. Accordingly, the evolution of this trait in C. hirsuta and A. thaliana shows partially conserved genetic mechanisms but is likely involved in adaptation to different environments.

Teosinte populations exhibit weak local adaptation to their rhizosphere biota despite strong effects of biota source on teosinte fitness and traits

While biotic interactions often impose selection, species and populations vary in whether they are locally adapted to biotic interactions. Evolutionary theory predicts that environmental conditions drive this variable local adaptation by altering the fitness impacts of species interactions. To investigate the influence of an environmental gradient on adaptation between a plant and its associated rhizosphere biota, we cross-combined teosinte (Zea mays ssp. mexicana) and rhizosphere biota collected across a gradient of decreasing temperature, precipitation, and nutrients in a greenhouse common garden experiment. We measured both fitness and phenotypes expected to be influenced by biota, including concentrations of nutrients in leaves. Independent, main effects of teosinte and biota source explained most variation in teosinte fitness and traits. For example, biota from warmer sites provided population-independent fitness benefits across teosinte hosts. Effects of biota that depended on teosinte genotype were often not specific to their local hosts, and most traits had similar relationships to fitness across biota treatments. However, we found weak patterns of local adaptation between teosinte and biota from colder sites, suggesting environmental gradients may alter the importance of local adaptation in teosinte-biota interactions, as evolutionary theory predicts.

Genome-Wide Diversity in Lowland and Highland Maize Landraces From Southern South America: Population Genetics Insights to Assist Conservation

Maize (Zea mays ssp. mays L.) landraces are traditional American crops with high genetic variability that conform a source of original alleles for conventional maize breeding. Northern Argentina, one the southernmost regions of traditional maize cultivation in the Americas, harbours around 57 races traditionally grown in two regions with contrasting environmental conditions, namely, the Andean mountains in the Northwest and the tropical grasslands and Atlantic Forest in the Northeast. These races encounter diverse threats to their genetic diversity and persistence in their regions of origin, with climate change standing out as one of the major challenges. In this work, we use genome-wide SNPs derived from ddRADseq to study the genetic diversity of individuals representing the five groups previously described for this area. This allowed us to distinguish two clearly differentiated gene pools, the highland northwestern maize (HNWA) and the floury northeastern maize (FNEA). Subsequently, we employed essential biodiversity variables at the genetic level, as proposed by the Group on Earth Observations Biodiversity Observation Network (GEO BON), to evaluate the conservation status of these two groups. This assessment encompassed genetic diversity (Pi), inbreeding coefficient (F) and effective population size (Ne). FNEA showed low Ne values and high F values, while HNWA showed low Ne values and low Pi values, indicating that further genetic erosion is imminent for these landraces. Outlier detection methods allowed identification of putative adaptive genomic regions, consistent with previously reported flowering-time loci and chromosomal regions displaying introgression from the teosinte Zea mays ssp. mexicana. Finally, species distribution models were obtained for two future climate scenarios, showing a notable reduction in the potential planting area of HNWA and a shift in the cultivation areas of FNEA. These results suggest that maize landraces from Northern Argentina may be unable to cope with climate change. Therefore, active conservation policies are advisable.

Unveiling the Genomic Symphony: Identification Cultivar-Specific Genes and Enhanced Insights on Sweet Sorghum Genomes Through Comprehensive superTranscriptomic Analysis

Sorghum (Sorghum bicolor (L.) Moench) is a multipurpose crop grown for food, fodder, and bioenergy production. Its cultivated varieties, along with their wild counterparts, contribute to the core genetic pool. Despite the availability of several re-sequenced sorghum genomes, a variable portion of sorghum genomes is not reported during reference genome assembly and annotation. The present analysis used 223 publicly available RNA-seq datasets from seven sweet sorghum cultivars to construct superTranscriptome. This approach yielded 45,864 Representative Transcript Assemblies (RTAs) that showcased intriguing Presence/Absence Variation (PAV) across 15 published sorghum genomes. We found 301 superTranscripts were exclusive to sweet sorghum, including 58 de novo genes encoded core and linker histones, zinc finger domains, glucosyl transferases, cellulose synthase, etc. The superTranscriptome added 2,802 new protein-coding genes to the Sweet Sorghum Reference Genome (SSRG), of which 559 code for different transcription factors (TFs). Our analysis revealed that MULE-like transposases were abundant in the sweet sorghum genome and could play a hidden role in the evolution of sweet sorghum. We observed large deletions in the D locus and terminal deletions in four other NAC encoding loci in the SSRG compared to its wild progenitor (353) suggesting non-functional NAC genes contributed to trait development in sweet sorghum. Moreover, superTranscript-based methods for Differential Exon Usage (DEU) and Differential Gene Expression (DGE) analyses were more accurate than those based on the SSRG. This study demonstrates that the superTranscriptome can enhance our understanding of fundamental sorghum mechanisms, improve genome annotations, and potentially even replace the reference genome.

Teosinte Pollen Drive guides maize diversification and domestication by RNAi

Selfish genetic elements contribute to hybrid incompatibility and bias or 'drive' their own transmission1,2. Chromosomal drive typically functions in asymmetric female meiosis, whereas gene drive is normally post-meiotic and typically found in males. Here, using single-molecule and single-pollen genome sequencing, we describe Teosinte Pollen Drive, an instance of gene drive in hybrids between maize (Zea mays ssp. mays) and teosinte mexicana (Z. mays ssp. mexicana) that depends on RNA interference (RNAi). 22-nucleotide small RNAs from a non-coding RNA hairpin in mexicana depend on Dicer-like 2 (Dcl2) and target Teosinte Drive Responder 1 (Tdr1), which encodes a lipase required for pollen viability. Dcl2, Tdr1 and the hairpin are in tight pseudolinkage on chromosome 5, but only when transmitted through the male. Introgression of mexicana into early cultivated maize is thought to have been critical to its geographical dispersal throughout the Americas3, and a tightly linked inversion in mexicana spans a major domestication sweep in modern maize4. A survey of maize traditional varieties and sympatric populations of teosinte mexicana reveals correlated patterns of admixture among unlinked genes required for RNAi on at least four chromosomes that are also subject to gene drive in pollen from synthetic hybrids. Teosinte Pollen Drive probably had a major role in maize domestication and diversification, and offers an explanation for the widespread abundance of 'self' small RNAs in the germ lines of plants and animals.

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.

DNA methylation variations underlie lettuce domestication and divergence

BACKGROUND

Lettuce (Lactuca sativa L.) is an economically important vegetable crop worldwide. Lettuce is believed to be domesticated from a single wild ancestor Lactuca serriola and subsequently diverged into two major morphologically distinct vegetable types: leafy lettuce and stem lettuce. However, the role of epigenetic variation in lettuce domestication and divergence remains largely unknown.

RESULTS

To understand the genetic and epigenetic basis underlying lettuce domestication and divergence, we generate single-base resolution DNA methylomes from 52 Lactuca accessions, including major lettuce cultivars and wild relatives. We find a significant increase of DNA methylation during lettuce domestication and uncover abundant epigenetic variations associated with lettuce domestication and divergence. Interestingly, DNA methylation variations specifically associated with leafy and stem lettuce are related to regulation and metabolic processes, respectively, while those associated with both types are enriched in stress responses. Moreover, we reveal that domestication-induced DNA methylation changes could influence expression levels of nearby and distal genes possibly through affecting chromatin accessibility and chromatin loop.

CONCLUSION

Our study provides population epigenomic insights into crop domestication and divergence and valuable resources for further domestication for diversity and epigenetic breeding to boost crop improvement.

Crop-to-wild gene flow in wild coffee species: the case of Coffea canephora in the Democratic Republic of the Congo

BACKGROUND AND AIMS

Plant breeders are increasingly turning to crop wild relatives (CWRs) to ensure food security in a rapidly changing environment. However, CWR populations are confronted with various human-induced threats, including hybridization with their nearby cultivated crops. This might be a particular problem for wild coffee species, which often occur near coffee cultivation areas. Here, we briefly review the evidence for wild Coffea arabica (cultivated as Arabica coffee) and Coffea canephora (cultivated as Robusta coffee) and then focused on C. canephora in the Yangambi region in the Democratic Republic of the Congo. There, we examined the geographical distribution of cultivated C. canephora and the incidence of hybridization between cultivated and wild individuals within the rainforest.

METHODS

We collected 71 C. canephora individuals from home gardens and 12 C. canephora individuals from the tropical rainforest in the Yangambi region and genotyped them using genotyping-by-sequencing (GBS). We compared the fingerprints with existing GBS data from 388 C. canephora individuals from natural tropical rainforests and the INERA Coffee Collection, a Robusta coffee field gene bank and the most probable source of cultivated genotypes in the area. We then established robust diagnostic fingerprints that genetically differentiate cultivated from wild coffee, identified cultivated-wild hybrids and mapped their geographical position in the rainforest.

KEY RESULTS

We identified cultivated genotypes and cultivated-wild hybrids in zones with clear anthropogenic activity, and where cultivated C. canephora in home gardens may serve as a source for crop-to-wild gene flow. We found relatively few hybrids and backcrosses in the rainforests.

CONCLUSIONS

The cultivation of C. canephora in close proximity to its wild gene pool has led to cultivated genotypes and cultivated-wild hybrids appearing within the natural habitats of C. canephora. Yet, given the high genetic similarity between the cultivated and wild gene pool, together with the relatively low incidence of hybridization, our results indicate that the overall impact in terms of risk of introgression remains limited so far.

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.

Teosinte Pollen Drive guides maize diversification and domestication by RNAi

Meiotic drivers subvert Mendelian expectations by manipulating reproductive development to bias their own transmission. Chromosomal drive typically functions in asymmetric female meiosis, while gene drive is normally postmeiotic and typically found in males. Using single molecule and single-pollen genome sequencing, we describe Teosinte Pollen Drive, an instance of gene drive in hybrids between maize (Zea mays ssp. mays) and teosinte mexicana (Zea mays ssp. mexicana), that depends on RNA interference (RNAi). 22nt small RNAs from a non-coding RNA hairpin in mexicana depend on Dicer-Like 2 (Dcl2) and target Teosinte Drive Responder 1 (Tdr1), which encodes a lipase required for pollen viability. Dcl2, Tdr1, and the hairpin are in tight pseudolinkage on chromosome 5, but only when transmitted through the male. Introgression of mexicana into early cultivated maize is thought to have been critical to its geographical dispersal throughout the Americas, and a tightly linked inversion in mexicana spans a major domestication sweep in modern maize. A survey of maize landraces and sympatric populations of teosinte mexicana reveals correlated patterns of admixture among unlinked genes required for RNAi on at least 4 chromosomes that are also subject to gene drive in pollen from synthetic hybrids. Teosinte Pollen Drive likely played a major role in maize domestication and diversification, and offers an explanation for the widespread abundance of "self" small RNAs in the germlines of plants and animals.

Introgression and persistence of cultivar alleles in wild carrot (Daucus carota) populations in the United States

PREMISE

Cultivated species and their wild relatives often hybridize in the wild, and the hybrids can survive and reproduce in some environments. However, it is unclear whether cultivar alleles are permanently incorporated into the wild genomes or whether they are purged by natural selection. This question is key to accurately assessing the risk of escape and spread of cultivar genes into wild populations.

METHODS

We used genomic data and population genomic methods to study hybridization and introgression between cultivated and wild carrot (Daucus carota) in the United States. We used single nucleotide polymorphisms (SNPs) obtained via genotyping by sequencing for 450 wild individuals from 29 wild georeferenced populations in seven states and 144 cultivars from the United States, Europe, and Asia.

RESULTS

Cultivated and wild carrot formed two genetically differentiated groups, and evidence of crop-wild admixture was detected in several but not all wild carrot populations in the United States. Two regions were identified where cultivar alleles were present in wild carrots: California and Nantucket Island (Massachusetts). Surprisingly, there was no evidence of introgression in some populations with a long-known history of sympatry with the crop, suggesting that post-hybridization barriers might prevent introgression in some areas.

CONCLUSIONS

Our results provide support for the introgression and long-term persistence of cultivar alleles in wild carrots populations. We thus anticipate that the release of genetically engineered (GE) cultivars would lead to the introduction and spread of GE alleles in wild carrot populations.

Fluctuation of ecological niches and geographic range shifts along chile pepper's domestication gradient

Domestication is an ongoing well-described process. However, while many have studied the changes domestication causes in plant genetics, few have explored its impact on the portion of the geographic landscape in which the plants exist. Therefore, the goal of this study was to understand how the process of domestication changed the geographic space suitable for chile pepper (Capsicum annuum) in its center of origin (domestication). C. annuum is a major crop species globally whose center of domestication, Mexico, has been well-studied. It provides a unique opportunity to explore the degree to which ranges of different domestication classes diverged and how these ranges might be altered by climate change. To this end, we created ecological niche models for four domestication classes (wild, semiwild, landrace, modern cultivar) based on present climate and future climate scenarios for 2050, 2070, and 2090. Considering present environment, we found substantial overlap in the geographic niches of all the domestication classes. Yet, environmental and geographic aspects of the current ranges did vary among classes. Wild and commercial varieties could grow in desert conditions, while landraces could not. With projections into the future, habitat was lost asymmetrically, with wild, semiwild, and landraces at greater risk of territorial declines than modern cultivars. Further, we identified areas where future suitability overlap between landraces and wilds is expected to be lost. While range expansion is widely associated with domestication, we found little support of a constant niche expansion (either in environmental or geographical space) throughout the domestication gradient in chile peppers in Mexico. Instead, particular domestication transitions resulted in loss, followed by capturing or recapturing environmental or geographic space. The differences in environmental characterization among domestication gradient classes and their future potential range shifts increase the need for conservation efforts to preserve landraces and semiwild genotypes.

Genomic diversity and population structure of teosinte (Zea spp.) and its conservation implications

The wild species of the genus Zea commonly named teosintes, comprise nine different taxa, distributed from northern Mexico to Costa Rica. Although this genus of plants has been extensively studied from a morphological, ecogeographical and genetic point of view, most contributions have been limited to the study of a few populations and taxa. To understand the great variability that exists between and within teosinte species, it is necessary to include the vast majority of known populations. In this context, the objective of this work was to evaluate the diversity and genomic structure of 276 teosinte populations. Molecular analyzes were performed with 3,604 plants and with data from 33,929 SNPs. The levels of genetic diversity by taxonomic group show a marked difference between species, races and sections, where the highest values of genomic diversity was found in ssp. parviglumis and ssp. mexicana. The lower values were obtained for the Luxuriantes section as well as ssp. huehuetenagensis of the section Zea. The results of structure show that there is a great genetic differentiation in all the taxonomic groups considered. For ssp. parviglumis and mexicana, which are the taxa with the largest number of populations, a marked genomic differentiation was found that is consistent with their geographic distribution patterns. These results showed a loss of diversity in several teosinte populations, making a strong case for further collection, and ex situ and in situ conservation. Also, this study highlights the importance of integrating genomic diversity and structure for the applications of conservation and management.

Genetic Incompatibilities and Evolutionary Rescue by Wild Relatives Shaped Grain Amaranth Domestication

Crop domestication and the subsequent expansion of crops have long been thought of as a linear process from a wild ancestor to a domesticate. However, evidence of gene flow from locally adapted wild relatives that provided adaptive alleles into crops has been identified in multiple species. Yet, little is known about the evolutionary consequences of gene flow during domestication and the interaction of gene flow and genetic load in crop populations. We study the pseudo-cereal grain amaranth that has been domesticated three times in different geographic regions of the Americas. We quantify the amount and distribution of gene flow and genetic load along the genome of the three grain amaranth species and their two wild relatives. Our results show ample gene flow between crop species and between crops and their wild relatives. Gene flow from wild relatives decreased genetic load in the three crop species. This suggests that wild relatives could provide evolutionary rescue by replacing deleterious alleles in crops. We assess experimental hybrids between the three crop species and found genetic incompatibilities between one Central American grain amaranth and the other two crop species. These incompatibilities might have created recent reproductive barriers and maintained species integrity today. Together, our results show that gene flow played an important role in the domestication and expansion of grain amaranth, despite genetic species barriers. The domestication of plants was likely not linear and created a genomic mosaic by multiple contributors with varying fitness effects for today's crops.

Genomic Resequencing Unravels the Genetic Basis of Domestication, Expansion, and Trait Improvement in Morus Atropurpurea

Mulberry is an economically important plant in the sericulture industry and traditional medicine. However, the genetic and evolutionary history of mulberry remains largely unknown. Here, this work presents the chromosome-level genome assembly of Morus atropurpurea (M. atropurpurea), originating from south China. Population genomic analysis using 425 mulberry accessions reveal that cultivated mulberry is classified into two species, M. atropurpurea and M. alba, which may have originated from two different mulberry progenitors and have independent and parallel domestication in north and south China, respectively. Extensive gene flow is revealed between different mulberry populations, contributing to genetic diversity in modern hybrid cultivars. This work also identifies the genetic architecture of the flowering time and leaf size. In addition, the genomic structure and evolution of sex-determining regions are identified. This study significantly advances the understanding of the genetic basis and domestication history of mulberry in the north and south, and provides valuable molecular markers of desirable traits for mulberry breeding.

Diamonds in the not-so-rough: Wild relative diversity hidden in crop genomes

Crop production is becoming an increasing challenge as the global population grows and the climate changes. Modern cultivated crop species are selected for productivity under optimal growth environments and have often lost genetic variants that could allow them to adapt to diverse, and now rapidly changing, environments. These genetic variants are often present in their closest wild relatives, but so are less desirable traits. How to preserve and effectively utilize the rich genetic resources that crop wild relatives offer while avoiding detrimental variants and maladaptive genetic contributions is a central challenge for ongoing crop improvement. This Essay explores this challenge and potential paths that could lead to a solution.

Adaptive and maladaptive introgression in grapevine domestication

Domesticated grapevines spread to Europe around 3,000 years ago. Previous studies have revealed genomic signals of introgression from wild to cultivated grapes in Europe, but the time, mode, genomic pattern, and biological effects of these introgression events have not been investigated. Here, we studied resequencing data from 345 samples spanning the distributional range of wild (Vitis vinifera ssp. sylvestris) and cultivated (V. vinifera ssp. vinifera) grapes. Based on machine learning-based population genetic analyses, we detected evidence for a single domestication of grapevine, followed by continuous gene flow between European wild grapes (EU) and cultivated grapes over the past ~2,000 y, especially from EU to wine grapes. We also inferred that soft-selective sweeps were the dominant signals of artificial selection. Gene pathways associated with the synthesis of aromatic compounds were enriched in regions that were both selected and introgressed, suggesting EU wild grapes were an important resource for improving the flavor of cultivated grapes. Despite the potential benefits of introgression in grape improvement, the introgressed fragments introduced a higher deleterious burden, with most deleterious SNPs and structural variants hidden in a heterozygous state. Cultivated wine grapes have benefited from adaptive introgression with wild grapes, but introgression has also increased the genetic load. In general, our study of beneficial and harmful effects of introgression is critical for genomic breeding of grapevine to take advantage of wild resources.

Population genomics unravels the Holocene history of bread wheat and its relatives

Deep knowledge of crop biodiversity is essential to improving global food security. Despite bread wheat serving as a keystone crop worldwide, the population history of bread wheat and its relatives, both cultivated and wild, remains elusive. By analysing whole-genome sequences of 795 wheat accessions, we found that bread wheat originated from the southwest coast of the Caspian Sea and underwent a slow speciation process, lasting ~3,300 yr owing to persistent gene flow from its relatives. Soon after, bread wheat spread across Eurasia and reached Europe, South Asia and East Asia ~7,000 to ~5,000 yr ago, shaping a diversified but occasionally convergent adaptive landscape in novel environments. By contrast, the cultivated relatives of bread wheat experienced a population decline by ~82% over the past ~2,000 yr due to the food choice shift of humans. Further biogeographical modelling predicted a continued population shrinking of many bread wheat relatives in the coming decades because of their vulnerability to the changing climate. These findings will guide future efforts in protecting and utilizing wheat biodiversity to enhance global wheat production.

Large haplotypes highlight a complex age structure within the maize pan-genome

The genomes of maize and other eukaryotes contain stable haplotypes in regions of low recombination. These regions, including centromeres, long heterochromatic blocks, and rDNA arrays, have been difficult to analyze with respect to their diversity and origin. Greatly improved genome assemblies are now available that enable comparative genomics over these and other nongenic spaces. Using 26 complete maize genomes, we developed methods to align intergenic sequences while excluding genes and regulatory regions. The centromere haplotypes (cenhaps) extend for megabases on either side of the functional centromere regions and appear as evolutionary strata, with haplotype divergence/coalescence times dating as far back as 450 thousand years ago (kya). Application of the same methods to other low recombination regions (heterochromatic knobs and rDNA) and all intergenic spaces revealed that deep coalescence times are ubiquitous across the maize pan-genome. Divergence estimates vary over a broad timescale with peaks at ∼16 and 300 kya, reflecting a complex history of gene flow among diverging populations and changes in population size associated with domestication. Cenhaps and other long haplotypes provide vivid displays of this ancient diversity.

Genome-wide signatures of the geographic expansion and breeding of soybean

Soybean is a leguminous crop that provides oil and protein. Exploring the genomic signatures of soybean evolution is crucial for breeding varieties with improved adaptability to environmental extremes. We analyzed the genome sequences of 2,214 soybeans and proposed a soybean evolutionary route, i.e., the expansion of annual wild soybean (Glycine soja Sieb. & Zucc.) from southern China and its domestication in central China, followed by the expansion and local breeding selection of its landraces (G. max (L.) Merr.). We observed that the genetic introgression in soybean landraces was mostly derived from sympatric rather than allopatric wild populations during the geographic expansion. Soybean expansion and breeding were accompanied by the positive selection of flowering time genes, including GmSPA3c. Our study sheds light on the evolutionary history of soybean and provides valuable genetic resources for its future breeding.

Distinct traces of mixed ancestry in western commercial pig genomes following gene flow from Chinese indigenous breeds

Studying gene flow between different livestock breeds will benefit the discovery of genes related to production traits and provide insight into human historical breeding. Chinese pigs have played an indispensable role in the breeding of Western commercial pigs. However, the differences in the timing and volume of the contribution of pigs from different Chinese regions to Western pigs are not yet apparent. In this paper, we combine the whole-genome sequencing data of 592 pigs from different studies and illustrate patterns of gene flow from Chinese pigs into Western commercial pigs. We describe introgression patterns from four distinct Chinese indigenous groups into five Western commercial groups. There were considerable differences in the number and length of the putative introgressed segments from Chinese pig groups that contributed to Western commercial pig breeds. The contribution of pigs from different Chinese geographical locations to a given western commercial breed varied more than that from a specific Chinese pig group to different Western commercial breeds, implying admixture within Europe after introgression. Within different Western commercial lines from the same breed, the introgression patterns from a given Chinese pig group seemed highly conserved, suggesting that introgression of Chinese pigs into Western commercial pig breeds mainly occurred at an early stage of breed formation. Finally, based on analyses of introgression signals, allele frequencies, and selection footprints, we identified a ∼2.65 Mb Chinese-derived haplotype under selection in Duroc pigs (CHR14: 95.68-98.33 Mb). Functional and phenotypic studies demonstrate that this PRKG1 haplotype is related to backfat and loin depth in Duroc pigs. Overall, we demonstrate that the introgression history of domestic pigs is complex and that Western commercial pigs contain distinct traces of mixed ancestry, likely derived from various Chinese pig breeds.

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.

Genetic introgression between different groups reveals the differential process of Asian cultivated rice

The Asian cultivated rice consists of two major subspecies: indica and japonica. There are already many reports about the existence of genetic introgression between the two subspecies. They propose some possible introgression-related genes from the comparison of population parameters. This study uses the genome-wide variation data of Asia cultivated rice to investigate their genetic introgression on the whole genome level. We detect a total of 13 significantly high introgression loci between the tropical japonica and indica populations. Two different methods are used to identify the genetic introgression regions. For most of the detected introgression regions, they generally get consistent results. Some previously known introgression genes are detected in the identified introgression loci, such as heat resistance gene TT1 and GLW7. The biological functions for these genetic introgression regions are annotated by the published QTL mapping results. We find that genetic introgression plays a vital role in the determination of both the phenotype and the domestication process of different groups. Our study also provides useful information and resources for the study of rice gene function and the domestication process.

Conflict over fertilization underlies the transient evolution of reinforcement

When two species meet in secondary contact, the production of low fitness hybrids may be prevented by the adaptive evolution of increased prezygotic isolation, a process known as reinforcement. Theoretical challenges to the evolution of reinforcement are generally cast as a coordination problem, i.e., "how can statistical associations between traits and preferences be maintained in the face of recombination?" However, the evolution of reinforcement also poses a potential conflict between mates. For example, the opportunity costs to hybridization may differ between the sexes or species. This is particularly likely for reinforcement based on postmating prezygotic (PMPZ) incompatibilities, as the ability to fertilize both conspecific and heterospecific eggs is beneficial to male gametes, but heterospecific mating may incur a cost for female gametes. We develop a population genetic model of interspecific conflict over reinforcement inspired by "gametophytic factors", which act as PMPZ barriers among Zea mays subspecies. We demonstrate that this conflict results in the transient evolution of reinforcement-after females adaptively evolve to reject gametes lacking a signal common in conspecific gametes, this gamete signal adaptively introgresses into the other population. Ultimately, the male gamete signal fixes in both species, and isolation returns to pre-reinforcement levels. We interpret geographic patterns of isolation among Z. mays subspecies considering these findings and suggest when and how this conflict can be resolved. Our results suggest that sexual conflict over fertilization may pose an understudied obstacle to the evolution of reinforcement.

An identical-by-descent segment harbors a 12-bp insertion determining fruit softening during domestication and speciation in Pyrus

BACKGROUND

Although the wild relatives of pear originated in southwest China, this fruit crop was independently domesticated and improved in Asia and Europe, and there are major phenotypic differences (e.g., maturity and fruit firmness) between Asian and European pears.  RESULTS: In this study, we examined the genomes of 113 diverse pear accessions using an identity-by-descent (IBD) approach to investigate how historical gene flow has shaped fruit firmness traits in Asian and European pears. We found a 3-Mbp IBD-enriched region (IBD-ER) that has undergone "convergent domestication" in both the Asian and European pear lineages, and a genome-wide association study (GWAS) of fruit firmness phenotypes strongly implicated the TRANSLOCON AT THE INNER CHLOROPLAST ENVELOPE55 (TIC55) locus within this 3-Mbp IBD-ER. Furthermore, we identified a tandem duplication that includes a 12-bp insertion located in the first exon of TIC55 that is uniquely present in Asian pears, and expression analysis showed that the pear TIC55 gene is highly expressed in Asian pear, while it is weakly or not expressed in European pear; this could contribute to the differences in fruit firmness between Asian and European pear fruits.

CONCLUSIONS

Our findings provide insights into how pear fruit softening has been impacted during domestication, and we identified candidate genes associated with fruit softening that can contribute to the breeding and improvement of pear and other fruit crops.

The genomic signature of wild-to-crop introgression during the domestication of scarlet runner bean (Phaseolus coccineus L.)

The scarlet runner bean (Phaseolus coccineus) is one of the five domesticated Phaseolus species. It is cultivated in small-scale agriculture in the highlands of Mesoamerica for its dry seeds and immature pods, and unlike the other domesticated beans, P. coccineus is an open-pollinated legume. Contrasting with its close relative, the common bean, few studies focusing on its domestication history have been conducted. Demographic bottlenecks associated with domestication might reduce genetic diversity and facilitate the accumulation of deleterious mutations. Conversely, introgression from wild relatives could be a source of variation. Using Genotyping by Sequencing data (79,286 single-nucleotide variants) from 237 cultivated and wild samples, we evaluated the demographic history of traditional varieties from different regions of Mexico and looked for evidence of introgression between sympatric wild and cultivated populations. Traditional varieties have high levels of diversity, even though there is evidence of a severe initial genetic bottleneck followed by a population expansion. Introgression from wild to domesticated populations was detected, which might contribute to the recovery of the genetic variation. Introgression has occurred at different times: constantly in the center of Mexico; recently in the North West; and anciently in the South. Several factors are acting together to increase and maintain genetic diversity in P. coccineus cultivars, such as demographic expansion and introgression. Wild relatives represent a valuable genetic resource and have played a key role in scarlet runner bean evolution via introgression into traditional varieties.

Alternative modes of introgression-mediated selection shaped crop adaptation to novel climates

Recent plant genomic studies provide fine-grained details on the evolutionary consequences of adaptive introgression during crop domestication. Modern genomic approaches and analytical methods now make it possible to better separate the introgression signal from the demographic signal thus providing a more comprehensive and complex picture of the role of introgression in local adaptation. Adaptive introgression has been fundamental for crop expansion and has involved complex patterns of gene flow. In addition to providing new and more favorable alleles of large effect, introgression during the early stages of domestication also increased allelic diversity at adaptive loci. Previous studies have largely underestimated the effect of such increased diversity following introgression. Recent genomic studies in wheat, potato, maize, grapevine, and ryegrass show that introgression of multiple genes, of as yet unknown effect, increased the effectiveness of purifying selection, and promoted disruptive or fluctuating selection in early cultivars and landraces. Historical selection processes associated with introgression from crop wild relatives provide an instructive analog for adaptation to current climate change and offer new avenues for crop breeding research that are expected to be instrumental for strengthening food security in the coming years.

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.

Genotype calling and haplotype inference from low coverage sequence data in heterozygous plant genome using HetMap

This study developed a new genotyping method that can accurately infer heterozygous genotype information from the complex plant genome sequence data, which helped discover new alleles in the association studies. Many software packages and pipelines had been developed to handle the sequence data of the model species. However, Genotyping from complex heterozygous plant genome needs further improvement on the previous methods. Here we present a new pipeline available at https://github.com/Ncgrhg/HetMapv1 ) for variant calling and missing genotype imputation from low coverage sequence data for heterozygous plant genomes. To check the performance of the HetMap on the real sequence data, HetMap was applied to both the F₁ hybrid rice population, which consists of 1495 samples and the wild rice population with 446 samples. The high coverage sequence data of four hybrid rice accessions and two wild rice accessions, which were also included in low coverage sequence data, were used to validate the accuracy of genotype inference. The validation results showed that HetMap archieved significant improvement in heterozygous genotype inference accuracy (13.65% for hybrid rice, 26.05% for wild rice) and total accuracy compared with similar software packages. The application of the new genotype with the genome-wide association study also showed improvement of association power in wild rice awn length phenotype. It could archive high genotype inference accuracy in low sequence coverage in a small population with both the natural and constructed recombination population. HetMap provided a powerful tool for the heterozygous plant genome sequence data analysis, which may help to discover new phenotype regions for the plant species with the complex heterozygous genome.

Maize dispersal patterns associated with different types of endosperm and migration of indigenous groups in lowland South America

BACKGROUND AND AIMS

The lowlands of South America appear to be remarkably important in the evolutionary history of maize, due to new evidence that suggests that maize dispersed from Mexico and arrived in this region in a state of partial domestication. This study aimed to identify dispersal patterns of maize genetic diversity in this part of the continent.

METHODS

A total of 170 maize accessions were characterized with 4398 single nucleotide polymorphisms (SNPs) and analysed to determine if maize dispersal was associated with types of endosperm and indigenous language families.

KEY RESULTS

Four genetic groups were identified in the discriminant analysis of principal components and five groups in the cluster analysis (neighbour-joining method). The groups were structured according to the predominance of endosperm types (popcorn, floury, flint/semi-flint). Spatial principal component analysis of genetic variation suggests different dispersal patterns for each endosperm type and can be associated with hypotheses of expansions of different indigenous groups.

CONCLUSIONS

From a possible origin in Southwestern Amazonia, different maize dispersal routes emerged: (1) towards Northern Amazonia, which continued towards the Caatinga and south-eastern Atlantic Forest (Floury); (2) towards Southern Brazil, passing through the Cerrado and Southern Atlantic Forest reaching the Pampa region (Floury); and (3) along the Atlantic Coast, following Tupi movements originating from two separate expansions: one (Tupinamba) from north to south, and the other (Guarani) in the opposite direction, from south to north (flint, floury and popcorn).

Demonstration of local adaptation in maize landraces by reciprocal transplantation

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

Phylogeny and disparate selection signatures suggest two genetically independent domestication events in pea (Pisum L.)

Domestication is considered a model of adaptation that can be used to draw conclusions about the modus operandi of selection in natural systems. Investigating domestication may give insights into how plants react to different intensities of human manipulation, which has direct implication for the continuing efforts of crop improvement. Therefore, scientists of various disciplines study domestication-related questions to understand the biological and cultural bases of the domestication process. We employed restriction site-associated DNA sequencing (RAD-seq) of 494 Pisum sativum (pea) samples from all wild and domesticated groups to analyze the genetic structure of the collection. Patterns of ancient admixture were investigated by analysis of admixture graphs. We used two complementary approaches, one diversity based and one based on differentiation, to detect the selection signatures putatively associated with domestication. An analysis of the subpopulation structure of wild P. sativum revealed five distinct groups with a notable geographic pattern. Pisum abyssinicum clustered unequivocally within the P. sativum complex, without any indication of hybrid origin. We detected 32 genomic regions putatively subjected to selection: 29 in P. sativum ssp. sativum and three in P. abyssinicum. The two domesticated groups did not share regions under selection and did not display similar haplotype patterns within those regions. Wild P. sativum is structured into well-diverged subgroups. Although Pisum sativum ssp. humile is not supported as a taxonomic entity, the so-called 'southern humile' is a genuine wild group. Introgression did not shape the variation observed within the sampled germplasm. The two domesticated pea groups display distinct genetic bases of domestication, suggesting two genetically independent domestication events.

A B73×Palomero Toluqueño mapping population reveals local adaptation in Mexican highland maize

Generations of farmer selection in the central Mexican highlands have produced unique maize varieties adapted to the challenges of the local environment. In addition to possessing great agronomic and cultural value, Mexican highland maize represents a good system for the study of local adaptation and acquisition of adaptive phenotypes under cultivation. In this study, we characterize a recombinant inbred line population derived from the B73 reference line and the Mexican highland maize variety Palomero Toluqueño. B73 and Palomero Toluqueño showed classic rank-changing differences in performance between lowland and highland field sites, indicative of local adaptation. Quantitative trait mapping identified genomic regions linked to effects on yield components that were conditionally expressed depending on the environment. For the principal genomic regions associated with ear weight and total kernel number, the Palomero Toluqueño allele conferred an advantage specifically in the highland site, consistent with local adaptation. We identified Palomero Toluqueño alleles associated with expression of characteristic highland traits, including reduced tassel branching, increased sheath pigmentation and the presence of sheath macrohairs. The oligogenic architecture of these three morphological traits supports their role in adaptation, suggesting they have arisen from consistent directional selection acting at distinct points across the genome. We discuss these results in the context of the origin of phenotypic novelty during selection, commenting on the role of de novo mutation and the acquisition of adaptive variation by gene flow from endemic wild relatives.

Genetic Architecture of Grain Yield-Related Traits in Sorghum and Maize

Grain size, grain number per panicle, and grain weight are crucial determinants of yield-related traits in cereals. Understanding the genetic basis of grain yield-related traits has been the main research object and nodal in crop science. Sorghum and maize, as very close C4 crops with high photosynthetic rates, stress tolerance and large biomass characteristics, are extensively used to produce food, feed, and biofuels worldwide. In this review, we comprehensively summarize a large number of quantitative trait loci (QTLs) associated with grain yield in sorghum and maize. We placed great emphasis on discussing 22 fine-mapped QTLs and 30 functionally characterized genes, which greatly hinders our deep understanding at the molecular mechanism level. This review provides a general overview of the comprehensive findings on grain yield QTLs and discusses the emerging trend in molecular marker-assisted breeding with these QTLs.

Extent of wild-to-crop interspecific introgression in grapevine (Vitis vinifera) as a consequence of resistance breeding and implications for the crop species definition

Over the past two centuries, introgression through repeated backcrossing has introduced disease resistance from wild grape species into the domesticated lineage Vitis vinifera subsp. sativa. Introgression lines are being cultivated over increasing vineyard surface areas, as their wines now rival in quality those obtained from preexisting varieties. There is, however, a lot of debate about whether and how wine laws defining commercial product categories, which are based on the classification of V. vinifera and interspecific hybrid grapes, should be revised to accommodate novel varieties that do not fit either category. Here, we developed a method of multilocus genotype analysis using short-read resequencing to identify haplotypic blocks of wild ancestry in introgression lines and quantify the physical length of chromosome segments free-of-introgression or with monoallelic and biallelic introgression. We used this genomic data to characterize species, hybrids and introgression lines and show that newly released resistant varieties contain 76.5-94.8% of V. vinifera DNA. We found that varietal wine ratings are not always commensurate with the percentage of V. vinifera ancestry and linkage drag of wild alleles around known resistance genes persists over at least 7.1-11.5 Mb, slowing down the recovery of the recurrent parental genome. This method also allowed us to identify the donor species of known resistance haplotypes, define the ancestry of wild genetic background in introgression lines with complex pedigrees, validate the ancestry of the historic varieties Concord and Norton, and unravel sample curation errors in public databases.

Demographic history and natural selection shape patterns of deleterious mutation load and barriers to introgression across Populus genome

Hybridization and resulting introgression are important processes shaping the tree of life and appear to be far more common than previously thought. However, how the genome evolution was shaped by various genetic and evolutionary forces after hybridization remains unresolved. Here we used whole-genome resequencing data of 227 individuals from multiple widespread Populus species to characterize their contemporary patterns of hybridization and to quantify genomic signatures of past introgression. We observe a high frequency of contemporary hybridization and confirm that multiple previously ambiguous species are in fact F₁ hybrids. Seven species were identified, which experienced different demographic histories that resulted in strikingly varied efficacy of selection and burdens of deleterious mutations. Frequent past introgression has been found to be a pervasive feature throughout the speciation of these Populus species. The retained introgressed regions, more generally, tend to contain reduced genetic load and to be located in regions of high recombination. We also find that in pairs of species with substantial differences in effective population size, introgressed regions are inferred to have undergone selective sweeps at greater than expected frequencies in the species with lower effective population size, suggesting that introgression likely have higher potential to provide beneficial variation for species with small populations. Our results, therefore, illustrate that demography and recombination have interplayed with both positive and negative selection in determining the genomic evolution after hybridization.

Plasticity of root anatomy during domestication of a maize-teosinte derived population

Maize (Zea mays L.) has undergone profound changes in root anatomy for environmental adaptation during domestication. However, the genetic mechanism of plasticity of maize root anatomy during the domestication process remains unclear. In this study, high-resolution mapping was performed for nine root anatomical traits using a maize-teosinte population (mexicana × Mo17) across three environments. Large genetic variations were detected for different root anatomical traits. The cortex, stele, aerenchyma areas, xylem vessel number, and cortical cell number had large variations across three environments, indicating high plasticity. Sixteen quantitative trait loci (QTL) were identified, including seven QTL with QTL × environment interaction (EIQTL) for high plasticity traits and nine QTL without QTL × environment interaction (SQTL). Most of the root loci were consistent with shoot QTL depicting domestication signals. Combining transcriptome and genome-wide association studies revealed that AUXIN EFFLUX CARRIER PIN-FORMED LIKE 4 (ZmPILS4) serves as a candidate gene underlying a major QTL of xylem traits. The near-isogenic lines (NILs) with lower expression of ZmPILS4 had 18-24% more auxin concentration in the root tips and 8-15% more xylem vessels. Nucleotide diversity values analysis in the promoter region suggested that ZmPILS4 was involved in maize domestication and adaptation. These results revealed the potential genetic basis of root anatomical plasticity during domestication.

Domestication reshaped the genetic basis of inbreeding depression in a maize landrace compared to its wild relative, teosinte

Inbreeding depression is the reduction in fitness and vigor resulting from mating of close relatives observed in many plant and animal species. The extent to which the genetic load of mutations contributing to inbreeding depression is due to large-effect mutations versus variants with very small individual effects is unknown and may be affected by population history. We compared the effects of outcrossing and self-fertilization on 18 traits in a landrace population of maize, which underwent a population bottleneck during domestication, and a neighboring population of its wild relative teosinte. Inbreeding depression was greater in maize than teosinte for 15 of 18 traits, congruent with the greater segregating genetic load in the maize population that we predicted from sequence data. Parental breeding values were highly consistent between outcross and selfed offspring, indicating that additive effects determine most of the genetic value even in the presence of strong inbreeding depression. We developed a novel linkage scan to identify quantitative trait loci (QTL) representing large-effect rare variants carried by only a single parent, which were more important in teosinte than maize. Teosinte also carried more putative juvenile-acting lethal variants identified by segregation distortion. These results suggest a mixture of mostly polygenic, small-effect partially recessive effects in linkage disequilibrium underlying inbreeding depression, with an additional contribution from rare larger-effect variants that was more important in teosinte but depleted in maize following the domestication bottleneck. Purging associated with the maize domestication bottleneck may have selected against some large effect variants, but polygenic load is harder to purge and overall segregating mutational burden increased in maize compared to teosinte.

Inbreeding depression is the reduction in fitness and vigor resulting from mating of close relatives observed in many plant and animal species. Mating of close relatives increases the probability that an individual inherits two non-functioning mutations at the same gene, resulting in lower fitness of such matings. We do not know the extent to which inbreeding depression is due to mutations with large-effects versus small-effect polygenic variants. We compared the effects of outcrossing and self-fertilization on 18 traits in a landrace population of maize, which underwent a population bottleneck during domestication, and a neighboring population of its wild relative teosinte. Inbreeding depression was greater in maize than teosinte for 15 of 18 traits and we found that this was consistent with higher predicted ‘genetic load’ in maize based solely on the evolutionary conservation of the sequence variants observed in the population. We also mapped genome positions associated with inbreeding depression, identifying more and larger-effect genetic variants in teosinte than maize. These results suggest that during domestication, some of the rare large-effect variants in teosinte were bred out, but many genetic variants of small effects on inbreeding depression increased in frequency maize.

Wheat speciation and adaptation: perspectives from reticulate evolution

Reticulate evolution through the interchanging of genetic components across organisms can impact significantly on the fitness and adaptation of species. Bread wheat (Triticum aestivum subsp. aestivum) is one of the most important crops in the world. Allopolyploid speciation, frequent hybridization, extensive introgression, and occasional horizontal gene transfer (HGT) have been shaping a typical paradigm of reticulate evolution in bread wheat and its wild relatives, which is likely to have a substantial influence on phenotypic traits and environmental adaptability of bread wheat. In this review, we outlined the evolutionary history of bread wheat and its wild relatives with a highlight on the interspecific hybridization events, demonstrating the reticulate relationship between species/subspecies in the genera Triticum and Aegilops. Furthermore, we discussed the genetic mechanisms and evolutionary significance underlying the introgression of bread wheat and its wild relatives. An in-depth understanding of the evolutionary process of Triticum species should be beneficial to future genetic study and breeding of bread wheat.

Defense mechanisms involving secondary metabolism in the grass family

Plants synthesize and accumulate a wide variety of compounds called secondary metabolites. Secondary metabolites serve as chemical barriers to protect plants from pathogens and herbivores. Antimicrobial secondary metabolites are accumulated to prevent pathogen infection. These metabolites are classified into phytoalexins (induced in response to pathogen attack) and phytoanticipins (present prior to pathogen infection). The antimicrobial compounds in the grass family (Poaceae) were studied from the viewpoint of evolution. The studies were performed at three hierarchies, families, genera, and species and include the following: 1) the distribution of benzoxazinoids (Bxs) in the grass family, 2) evolutionary replacement of phytoanticipins from Bxs to hydroxycinnamic acid amide dimers in the genus Hordeum, and 3) chemodiversity of flavonoid and diterpenoid phytoalexins in rice. These studies demonstrated dynamic changes in secondary metabolism during evolution, indicating the adaptation of plants to their environment by repeating scrap-and-build cycles.

Pervasive hybridization with local wild relatives in Western European grapevine varieties

Grapevine (Vitis vinifera L.) diversity richness results from a complex domestication history over multiple historical periods. Here, we used whole-genome resequencing to elucidate different aspects of its recent evolutionary history. Our results support a model in which a central domestication event in grapevine was followed by postdomestication hybridization with local wild genotypes, leading to the presence of an introgression signature in modern wine varieties across Western Europe. The strongest signal was associated with a subset of Iberian grapevine varieties showing large introgression tracts. We targeted this study group for further analysis, demonstrating how regions under selection in wild populations from the Iberian Peninsula were preferentially passed on to the cultivated varieties by gene flow. Examination of underlying genes suggests that environmental adaptation played a fundamental role in both the evolution of wild genotypes and the outcome of hybridization with cultivated varieties, supporting a case of adaptive introgression in grapevine.

Anchoring alien chromosome segment substitutions bearing gene(s) for resistance to mustard aphid in Brassica juncea-B. fruticulosa introgression lines and their possible disruption through gamma irradiation

KEY MESSAGE

Heavy doses of gamma irradiation can reduce linkage drag by disrupting large sized alien translocations and promoting exchanges between crop and wild genomes. Resistance to mustard aphid (Lipaphis erysimi) infestation was significantly improved in Brassica juncea through B. juncea-B. fruticulosa introgression. However, linkage drag caused by introgressed chromatin fragments has so far prevented the deployment of this resistance source in commercial cultivars. We investigated the patterns of donor chromatin segment substitutions in the introgression lines (ILs) through genomic in situ hybridization (GISH) coupled with B. juncea chromosome-specific oligonucleotide probes. These allowed identification of large chromosome translocations from B. fruticulosa in the terminal regions of chromosomes A05, B02, B03 and B04 in three founder ILs (AD-64, 101 and 104). Only AD-101 carried an additional translocation at the sub-terminal to intercalary position in both homologues of chromosome A01. We validated these translocations with a reciprocal blast hit analysis using shotgun sequencing of three ILs and species-specific contigs/scaffolds (kb sized) from a de novo assembly of B. fruticulosa. Alien segment substitution on chromosome A05 could not be validated. Current studies also endeavoured to break linkage drag by exposing seeds to a heavy dose (200kR) of gamma radiation. Reduction in the size of introgressed chromatin fragments was observed in many M3 plants. There was a complete loss of the alien chromosome fragment in one instance. A few M3 plants with novel patterns of chromosome segment substitutions displayed improved agronomic performance coupled with resistance to mustard aphid. SNPs in such genomic spaces should aid the development of markers to track introgressed DNA and allow application in plant breeding.

Selective sorting of ancestral introgression in maize and teosinte along an elevational cline

While often deleterious, hybridization can also be a key source of genetic variation and pre-adapted haplotypes, enabling rapid evolution and niche expansion. Here we evaluate these opposing selection forces on introgressed ancestry between maize (Zea mays ssp. mays) and its wild teosinte relative, mexicana (Zea mays ssp. mexicana). Introgression from ecologically diverse teosinte may have facilitated maize's global range expansion, in particular to challenging high elevation regions (> 1500 m). We generated low-coverage genome sequencing data for 348 maize and mexicana individuals to evaluate patterns of introgression in 14 sympatric population pairs, spanning the elevational range of mexicana, a teosinte endemic to the mountains of Mexico. While recent hybrids are commonly observed in sympatric populations and mexicana demonstrates fine-scale local adaptation, we find that the majority of mexicana ancestry tracts introgressed into maize over 1000 generations ago. This mexicana ancestry seems to have maintained much of its diversity and likely came from a common ancestral source, rather than contemporary sympatric populations, resulting in relatively low FST between mexicana ancestry tracts sampled from geographically distant maize populations. Introgressed mexicana ancestry in maize is reduced in lower-recombination rate quintiles of the genome and around domestication genes, consistent with pervasive selection against introgression. However, we also find mexicana ancestry increases across the sampled elevational gradient and that high introgression peaks are most commonly shared among high-elevation maize populations, consistent with introgression from mexicana facilitating adaptation to the highland environment. In the other direction, we find patterns consistent with adaptive and clinal introgression of maize ancestry into sympatric mexicana at many loci across the genome, suggesting that maize also contributes to adaptation in mexicana, especially at the lower end of its elevational range. In sympatric maize, in addition to high introgression regions we find many genomic regions where selection for local adaptation maintains steep gradients in introgressed mexicana ancestry across elevation, including at least two inversions: the well-characterized 14 Mb Inv4m on chromosome 4 and a novel 3 Mb inversion Inv9f surrounding the macrohairless1 locus on chromosome 9. Most outlier loci with high mexicana introgression show no signals of sweeps or local sourcing from sympatric populations and so likely represent ancestral introgression sorted by selection, resulting in correlated but distinct outcomes of introgression in different contemporary maize populations.

The evolutionary genomics of species' responses to climate change

Climate change is a threat to biodiversity. One way that this threat manifests is through pronounced shifts in the geographical range of species over time. To predict these shifts, researchers have primarily used species distribution models. However, these models are based on assumptions of niche conservatism and do not consider evolutionary processes, potentially limiting their accuracy and value. To incorporate evolution into the prediction of species' responses to climate change, researchers have turned to landscape genomic data and examined information about local genetic adaptation using climate models. Although this is an important advancement, this approach currently does not include other evolutionary processes-such as gene flow, population dispersal and genomic load-that are critical for predicting the fate of species across the landscape. Here, we briefly review the current practices for the use of species distribution models and for incorporating local adaptation. We next discuss the rationale and theory for considering additional processes, reviewing how they can be incorporated into studies of species' responses to climate change. We summarize with a conceptual framework of how manifold layers of information can be combined to predict the potential response of specific populations to climate change. We illustrate all of the topics using an exemplar dataset and provide the source code as potential tutorials. This Perspective is intended to be a step towards a more comprehensive integration of population genomics with climate change science.

Rewilding crops for climate resilience: economic analysis and de novo domestication strategies

To match predicted population growth, annual food production should be doubled by 2050. This is not achievable by current agronomical and breeding practices, due to the impact of climate changes and associated abiotic stresses on agricultural production systems. Here, we analyze the impact of global climate trends on crop productivity and show that the overall loss in crop production from climate-driven abiotic stresses may exceed US$170 billion year-1 and represents a major threat to global food security. We also show that abiotic stress tolerance had been present in wild progenitors of modern crops but was lost during their domestication. We argue for a major shift in our paradigm of crop breeding, focusing on climate resilience, and call for a broader use of wild relatives as a major tool in this process. We argue that, while molecular tools are currently in place to harness the potential of climate-resilient genes present in wild relatives, the complex polygenic nature of tolerance traits remains a major bottleneck in this process. Future research efforts should be focused not only on finding appropriate wild relatives but also on development of efficient cell-based high-throughput phenotyping platforms allowing assessment of the in planta operation of key genes.

Introgression among North American wild grapes (Vitis) fuels biotic and abiotic adaptation

BACKGROUND

Introgressive hybridization can reassort genetic variants into beneficial combinations, permitting adaptation to new ecological niches. To evaluate evolutionary patterns and dynamics that contribute to introgression, we investigate six wild Vitis species that are native to the Southwestern United States and useful for breeding grapevine (V. vinifera) rootstocks.

RESULTS

By creating a reference genome assembly from one wild species, V. arizonica, and by resequencing 130 accessions, we focus on identifying putatively introgressed regions (pIRs) between species. We find six species pairs with signals of introgression between them, comprising up to ~ 8% of the extant genome for some pairs. The pIRs tend to be gene poor, located in regions of high recombination and enriched for genes implicated in disease resistance functions. To assess potential pIR function, we explore SNP associations to bioclimatic variables and to bacterial levels after infection with the causative agent of Pierce's disease (Xylella fastidiosa). pIRs are enriched for SNPs associated with both climate and bacterial levels, suggesting that introgression is driven by adaptation to biotic and abiotic stressors.

CONCLUSIONS

Altogether, this study yields insights into the genomic extent of introgression, potential pressures that shape adaptive introgression, and the evolutionary history of economically important wild relatives of a critical crop.

Molecular Parallelism Underlies Convergent Highland Adaptation of Maize Landraces

Convergent phenotypic evolution provides some of the strongest evidence for adaptation. However, the extent to which recurrent phenotypic adaptation has arisen via parallelism at the molecular level remains unresolved, as does the evolutionary origin of alleles underlying such adaptation. Here, we investigate genetic mechanisms of convergent highland adaptation in maize landrace populations and evaluate the genetic sources of recurrently selected alleles. Population branch excess statistics reveal substantial evidence of parallel adaptation at the level of individual single-nucleotide polymorphism (SNPs), genes, and pathways in four independent highland maize populations. The majority of convergently selected SNPs originated via migration from a single population, most likely in the Mesoamerican highlands, while standing variation introduced by ancient gene flow was also a contributor. Polygenic adaptation analyses of quantitative traits reveal that alleles affecting flowering time are significantly associated with elevation, indicating the flowering time pathway was targeted by highland adaptation. In addition, repeatedly selected genes were significantly enriched in the flowering time pathway, indicating their significance in adapting to highland conditions. Overall, our study system represents a promising model to study convergent evolution in plants with potential applications to crop adaptation across environmental gradients.

The genomic consequences of hybridization

In the past decade, advances in genome sequencing have allowed researchers to uncover the history of hybridization in diverse groups of species, including our own. Although the field has made impressive progress in documenting the extent of natural hybridization, both historical and recent, there are still many unanswered questions about its genetic and evolutionary consequences. Recent work has suggested that the outcomes of hybridization in the genome may be in part predictable, but many open questions about the nature of selection on hybrids and the biological variables that shape such selection have hampered progress in this area. We synthesize what is known about the mechanisms that drive changes in ancestry in the genome after hybridization, highlight major unresolved questions, and discuss their implications for the predictability of genome evolution after hybridization.