Post-Domestication Selection in the Maize Starch Pathway

Modern crops have usually experienced domestication selection and subsequent genetic improvement (post-domestication selection). Chinese waxy maize, which originated from non-glutinous domesticated maize (Zea mays ssp. mays), provides a unique model for investigating the post-domestication selection of maize. In this study, the genetic diversity of six key genes in the starch pathway was investigated in a glutinous population that included 55 Chinese waxy accessions, and a selective bottleneck that resulted in apparent reductions in diversity in Chinese waxy maize was observed. Significant positive selection in waxy (wx) but not amylose extender1 (ae1) was detected in the glutinous population, in complete contrast to the findings in non-glutinous maize, which indicated a shift in the selection target from ae1 to wx during the improvement of Chinese waxy maize. Our results suggest that an agronomic trait can be quickly improved into a target trait with changes in the selection target among genes in a crop pathway.     

Genetic Diversity and Molecular Evolution of Chinese Waxy Maize Germplasm

Waxy maize (Zea mays L. var. certaina Kulesh), with many excellent characters in terms of starch composition and economic value, has grown in China for a long history and its production has increased dramatically in recent decades. However, the evolution and origin of waxy maize still remains unclear. We studied the genetic diversity of Chinese waxy maize including typical landraces and inbred lines by SSR analysis and the results showed a wide genetic diversity in the Chinese waxy maize germplasm. We analyzed the origin and evolution of waxy maize by sequencing 108 samples, and downloading 52 sequences from GenBank for the waxy locus in a number of accessions from genus Zea. A sharp reduction of nucleotide diversity and significant neutrality tests (Tajima’s D and Fu and Li’s F*) were observed at the waxy locus in Chinese waxy maize but not in nonglutinous maize. Phylogenetic analysis indicated that Chinese waxy maize originated from the cultivated flint maize and most of the modern waxy maize inbred lines showed a distinct independent origin and evolution process compared with the germplasm from Southwest China. The results indicated that an agronomic trait can be quickly improved to meet production demand by selection.     

The molecular basis of mutations at the Waxy locus from Amaranthus caudatus L.: evolution of the waxy phenotype in three species of grain amaranth

Polymorphisms at the Waxy locus of Amaranthus caudatus L. collected from a wide range of regions were used to investigate genetic diversity and mutation sites. A comparison of the Waxy locus revealed a very high level of sequence conservation. This result clearly showed low environmental and evolutionary variability in the Waxy gene. We also performed screening to confirm the mutation sites in the coding sequences of all accessions. The results indicate that one insertion in the coding region of Waxy genes was responsible for the change in perisperm starch leading to the waxy phenotype in all accessions of this species, and thus that a single mutation event altered the regulation of the Waxy gene during the domestication of this crop. In addition, phylogenetic analysis showed that waxy phenotypes within each of three species, A. caudatus, A. cruentus and A. hypochondriacus, originated separately or differentiated from nonwaxy phenotypes of each species through a single mutational event (i.e., a frame shift or base substitution). We also compared obvious structural features of the coding sequence of waxy and nonwaxy phenotypes with those of low-amylose phenotypes in A. caudatus. The Waxy coding sequences of low-amylose phenotypes do not show polymorphisms and are identical with those of waxy phenotypes. This could mean that there is another gene that encodes a key enzyme responsible for amylose synthesis as the elementary quantity in tissues other than perisperm in A. caudatus.     

Target enrichment sequencing of 307 germplasm accessions identified ancestry of ancient and modern hybrids and signatures of adaptation and selection in sugarcane (Saccharum spp.), a ‘sweet’ crop with ‘bitter’ genomes

Sugarcane (Saccharum spp.) is a highly energy‐efficient crop primarily for sugar and bio‐ethanol production. Sugarcane genetics and cultivar improvement have been extremely challenging largely due to its complex genomes with high polyploidy levels. In this study, we deeply sequenced the coding regions of 307 sugarcane germplasm accessions. Nearly five million sequence variations were catalogued. The average of 98× sequence depth enabled different allele dosages of sequence variation to be differentiated in this polyploid collection. With selected high‐quality genome‐wide SNPs, we performed population genomic studies and environmental association analysis. Results illustrated that the ancient sugarcane hybrids, S. barberi and S. sinense, and modern sugarcane hybrids are significantly different in terms of genomic compositions, hybridization processes and their potential ancestry contributors. Linkage disequilibrium (LD) analysis showed a large extent of LD in sugarcane, with 962.4 Kbp, 2739.2 Kbp and 3573.6 Kbp for S. spontaneum, S. officinarum and modern S. hybrids respectively. Candidate selective sweep regions and genes were identified during domestication and historical selection processes of sugarcane in addition to genes associated with environmental variables at the original locations of the collection. This research provided an extensive amount of genomic resources for sugarcane community and the in‐depth population genomic analyses shed light on the breeding and evolution history of sugarcane, a highly polyploid species.     

Domestication and geographic origin of Oryza sativa in China: insights from multilocus analysis of nucleotide variation of O. sativa and O. rufipogon

Previous studies have indicated that China is one of the domestication centres of Asian cultivated rice (Oryza sativa), and common wild rice (O. rufipogon) is the progenitor of O. sativa. However, the number of domestication times and the geographic origin of Asian cultivated rice in China are still under debate. In this study, 100 accessions of Asian cultivated rice and 111 accessions of common wild rice in China were selected to examine the relationship between O. sativa and O. rufipogon and thereby infer the domestication and evolution of O. sativa in China through sequence analyses of six gene regions, trnC‐ycf6 in chloroplast genomes, cox3 in mitochondrial genomes and ITS, Ehd1, Waxy, Hd1 in nuclear genomes. The results indicated that the two subspecies of O. sativa (indica and japonica) were domesticated independently from different populations of O. rufipogon with gene flow occurring later from japonica to indica; Southern China was the genetic diversity centre of O. rufipogon, and the Pearl River basin near the Tropic of Cancer was the domestication centre of O. sativa in China.     

Selection under domestication: evidence for a sweep in the rice waxy genomic region

Rice (Oryza sativa) was cultivated by Asian Neolithic farmers >11,000 years ago, and different cultures have selected for divergent starch qualities in the rice grain during and after the domestication process. An intron 1 splice donor site mutation of the Waxy gene is responsible for the absence of amylose in glutinous rice varieties. This mutation appears to have also played an important role in the origin of low amylose, nonglutinous temperate japonica rice varieties, which form a primary component of Northeast Asian cuisines. Waxy DNA sequence analyses indicate that the splice donor mutation is prevalent in temperate japonica rice varieties, but rare or absent in tropical japonica, indica, aus, and aromatic varieties. Sequence analysis across a 500-kb genomic region centered on Waxy reveals patterns consistent with a selective sweep in the temperate japonicas associated with the mutation. The size of the selective sweep (>250 kb) indicates very strong selection in this region, with an inferred selection coefficient that is higher than similar estimates from maize domestication genes or wild species. These findings demonstrate that selection pressures associated with crop domestication regimes can exceed by one to two orders of magnitude those observed for genes under even strong selection in natural systems.     

Origin and evolution of the waxy phenotype in Amaranthus hypochondriacus: evidence from the genetic diversity in the Waxy locus

The existence of polymorphism in the Waxy locus in a large gene pool of 53 strains with various waxy phenotypes from samples of Amaranthus hypochondriacus collected from different regions was investigated in an origin-and-evolution study. First, we screened all strains for a mutation point (G–A polymorphism in exon 6) by using PCR–RFLP and/or direct sequence analysis. The results showed that the nonsense mutation in the coding region (exon 6) of the Waxy gene was responsible for the change in perisperm starch, leading to a waxy phenotype in all strains. Second, phylogenetic analysis, which was based on the Waxy variation, indicated diverse waxy types occurring separately and independently in certain domesticated regions in Mexico. Finally, we designated nine molecular types by comparing obvious structural variations in the coding region of the Waxy gene. Among the molecular types, A. hypochondriacus contained Type III in three subtypes with the waxy phenotype, with evolutionary routes that could originate from Type II in accordance with G–A polymorphism. In addition, these types had the same mutation points by which the Waxy gene was converted into the waxy phenotype. Therefore, the present results showed that the nonsense mutation is a unique event in the evolution of waxy phenotypes in this crop. This study will provide useful information for understanding the evolutionary process of the waxy phenotype.     

Genetic diversity and domestication history of African rice (<Emphasis Type="Italic">Oryza glaberrima</Emphasis>) as inferred from multiple gene sequences

Nucleotide variation in 14 unlinked nuclear genes was investigated in species-wide samples of African rice (Oryza glaberrima) and its wild progenitor (O. barthii). Average estimates of nucleotide diversity were extremely low in both species (θ _sil = 0.0007 for O. glaberrima; θ _sil = 0.0024 for O. barthii). About 70% less diversity was found in O. glaberrima than in its progenitor O. barthii. Coalescent simulation indicated that such dramatic reduction of nucleotide diversity in African rice could be explained mainly by a severe bottleneck during its domestication. The progenitor of African rice maintained also low genetic diversity, which may be attributed to small effective population size in O. barthii. Self-pollinating would be another factor leading to the unusually low diversity in both species. Genealogical analyses showed that all O. glaberrima accessions formed a strongly supported cluster with seven O. barthii individuals that were sampled exclusively from the proposed domestication centers of African rice. Population structure and principal component analyses found that the O. glaberrima group was homogeneous with no obvious genetic subdivision, in contrast to the heterogeneous O. barthii cluster. These findings support a single domestication origin of African rice in areas of the Upper Niger and Sahelian Rivers.     

Cloning and sequence diversity analysis of <Emphasis Type="Italic">GmHs1</Emphasis><Superscript><Emphasis Type="Italic">pro-1</Emphasis></Superscript> in Chinese domesticated and wild soybeans

Soybean cyst nematode (SCN) (Heterodera glycines Ichinohe) is the most important pathogen in soybean production worldwide and causes substantial yield losses. An apparent narrow genetic base of SCN resistance was observed in current elite soybean cultivars, and searching for novel SCN resistance genes as well as novel resistance sources rather than focusing on the two important genes rhg1 and Rhg4 has become another major objective in soybean research. In the present paper we report a 1,477 bp Hs1 ^pro-1 homolog, named GmHs1 ^pro-1 . This gene was cloned from soybean variety Wenfeng 7 based on information for Hs1 ^pro-1 , a beet cyst nematode resistance gene in sugar beet. It has two domains, Hs1pro-1_N and Hs1pro-1_C, both of which are believed to confer resistance to nematodes. Of the 1,477 bp sequence in GmHs1 ^pro-1 , an open reading frame of 1,314 bp, encoding a protein with 437 amino acids, was flanked by a 5′-untranslated region of 27 bp and a 3′-untranslated region of 135 bp. Fourteen single-nucleotide polymorphisms (SNPs) were observed in 44 soybean accessions including 23 wild soybeans, 8 landraces, and 13 soybean varieties (or lines), among which 5 in wild soybeans and 3 in landrace accessions were unique. Sequence diversity analysis on the 44 soybean accessions showed π = 0.00168 and θ = 0.00218 for GmHs1 ^pro-1 ; landraces had the highest diversity, followed by wild soybeans, with varieties (or lines) having the lowest. Although we did not detect a significant effect of selection on GmHs1 ^pro-1 in the three populations, sequence diversity, unique SNPs, and phylogenetic analysis indicated a slight domestication bottleneck and an intensive selection bottleneck. High sequence diversity, more unique SNPs, and broader representation across the phylogenetic tree in wild soybeans and landraces indicated that wild collections and landrace accessions are invaluable germplasm for broadening the genetic base of elite soybean varieties resistant to SCN. C. Yuan and G. Zhou contributed to this paper equally.     

Contrasting Patterns of Sequence Evolution at the Functionally Redundant <Emphasis Type="Italic">bric à brac</Emphasis> Paralogs in <Emphasis Type="Italic">Drosophila melanogaster</Emphasis>

Genes with overlapping expression and function may gradually diverge despite retaining some common functions. To test whether such genes show distinct patterns of molecular evolution within species, we examined sequence variation at the bric à brac (bab) locus of Drosophila melanogaster. This locus is composed of two anciently duplicated paralogs, bab1 and bab2, which are involved in patterning the adult abdomen, legs, and ovaries. We have sequenced the 148 kb genomic region spanning the bab1 and bab2 genes from 94 inbred lines of D. melanogaster sampled from a single location. Two non-coding regions, one in each paralog, appear to be under selection. The strongest evidence of directional selection is found in a region of bab2 that has no known functional role. The other region is located in the bab1 paralog and is known to contain a cis-regulatory element that controls sex-specific abdominal pigmentation. The coding region of bab1 appears to be under stronger functional constraint than the bab2 coding sequences. Thus, the two paralogs are evolving under different selective regimes in the same natural population, illuminating the different evolutionary trajectories of partially redundant duplicate genes.     

Genomic and phenotypic evidence for an incomplete domestication of South American grain amaranth (Amaranthus caudatus)

The domestication syndrome comprises phenotypic changes that differentiate crops from their wild ancestors. We compared the genomic variation and phenotypic differentiation of the two putative domestication traits seed size and seed colour of the grain amaranth Amaranthus caudatus, which is an ancient crop of South America, and its two close wild relatives and putative ancestors A. hybridus and A. quitensis. Genotyping 119 accessions of the three species from the Andean region using genotyping by sequencing (GBS) resulted in 9485 SNPs that revealed a strong genetic differentiation of cultivated A. caudatus from its two relatives. A. quitensis and A. hybridus accessions did not cluster by their species assignment but formed mixed groups according to their geographic origin in Ecuador and Peru, respectively. A. caudatus had a higher genetic diversity than its close relatives and shared a high proportion of polymorphisms with their wild relatives consistent with the absence of a strong bottleneck or a high level of recent gene flow. Genome sizes and seed sizes were not significantly different between A. caudatus and its relatives, although a genetically distinct group of A. caudatus from Bolivia had significantly larger seeds. We conclude that despite a long history of human cultivation and selection for white grain colour, A. caudatus shows a weak genomic and phenotypic domestication syndrome and proposes that it is an incompletely domesticated crop species either because of weak selection or high levels of gene flow from its sympatric close undomesticated relatives that counteracted the fixation of key domestication traits.     

Physical leaf defenses – altered by Zea life‐history evolution,domestication, and breeding – mediate oviposition preference of a specialist leafhopper

Plant anti‐herbivore defenses are known to be affected by life‐history evolution, as well as by domestication and breeding in the case of crop species. A suite of plants from the maize genus Zea (Poaceae) and the specialist herbivore Dalbulus maidis (DeLong & Wolcott) (Hemiptera: Cicadellidae) were used to test the hypothesis that anti‐herbivore defenses are affected by plant life‐history evolution and human intervention through domestication and breeding for high yield. The suite of plants included a maize (Zea mays ssp. mays L.) commercial hybrid, a maize landrace, two populations of the annual Balsas teosinte (Z. mays ssp. parviglumis Iltis & Doebley), and perennial teosinte (Z. diploperennis Iltis, Doebley & Guzman). Leaf toughness, pubescence, and oviposition preference were compared among the suite of host plants looking for effects of transitions in life history (i.e., from perennial to annual life cycle), domestication (i.e., from wild annual to domesticated annual), and breeding (i.e., from landrace to hybrid maize) on defense against D. maidis. Results on leaf toughness suggested that the life‐history and domestication transitions weakened the plant's resistance to penetration by the mouthparts and ovipositor of D. maidis, whereas results on pubescence suggested that this putative defense was strengthened with the breeding transition, contrary to expectations. Results on oviposition preference of D. maidis coincided with the expectation that life‐history and domestication transitions would lead to preference for Balsas teosinte over perennial teosinte, and of landrace maize over Balsas teosinte. Also, a negative correlation suggested that oviposition preference is significantly influenced by leaf toughness. Overall, the results suggested that Zea defenses against the specialist herbivore D. maidis were variably affected by plant life‐history evolution, domestication, and breeding, and that chemical defense may play a role in Zea defense against D. maidis because leaf toughness and pubescence only partially explained its host preferences.     

Identification of Genetic Differentiation between Waxy and Common Maize by SNP Genotyping

Waxy maize (Zea mays L. var. ceratina) is an important vegetable and economic crop that is thought to have originated from cultivated flint maize and most recently underwent divergence from common maize. In this study, a total of 110 waxy and 110 common maize inbred lines were genotyped with 3072 SNPs to evaluate the genetic diversity, population structure, and linkage disequilibrium decay as well as identify putative loci that are under positive selection. The results revealed abundant genetic diversity in the studied panel and that genetic diversity was much higher in common than in waxy maize germplasms. Principal coordinate analysis and neighbor-joining cluster analysis consistently classified the 220 accessions into two major groups and a mixed group with mixed ancestry. Subpopulation structure in both waxy and common maize sets were associated with the germplasm origin and corresponding heterotic groups. The LD decay distance (1500–2000 kb) in waxy maize was lower than that in common maize. Fourteen candidate loci were identified as under positive selection between waxy and common maize at the 99% confidence level. The information from this study can assist waxy maize breeders by enhancing parental line selection and breeding program design.     

Genome‐wide analysis reveals artificial selection on coat colour and reproductive traits in Chinese domestic pigs

Pigs from Asia and Europe were independently domesticated from c. 9000 years ago. During this period, strong artificial selection has led to dramatic phenotypic changes in domestic pigs. However, the genetic basis underlying these morphological and behavioural adaptations is relatively unknown, particularly for indigenous Chinese pigs. Here, we performed a genome‐wide analysis to screen 196 regions with selective sweep signals in Tongcheng pigs, which are a typical indigenous Chinese breed. Genes located in these regions have been found to be involved in lipid metabolism, melanocyte differentiation, neural development and other biological processes, which coincide with the evolutionary phenotypic changes in this breed. A synonymous substitution, c.669T>C, in ESR1, which colocalizes with a major quantitative trait locus for litter size, shows extreme differences in allele frequency between Tongcheng pigs and wild boars. Notably, the variant C allele in this locus exhibits high allele frequency in most Chinese populations, suggesting a consequence of positive selection. Five genes (PRM1, PRM2, TNP2, GPR149 and JMJD1C) related to reproductive traits were found to have high haplotype similarity in Chinese breeds. Two selected genes, MITF and EDNRB, are implied to shape the two‐end black colour trait in Tongcheng pig. Subsequent SNP microarray studies of five Chinese white‐spotted breeds displayed a concordant signature at both loci, suggesting that these two genes are responsible for colour variations in Chinese breeds. Utilizing massively parallel sequencing, we characterized the candidate sites that adapt to artificial and environmental selections during the Chinese pig domestication. This study provides fundamental proof for further research on the evolutionary adaptation of Chinese pigs.     

Maize Domestication and Anti-Herbivore Defences: Leaf-Specific Dynamics during Early Ontogeny of Maize and Its Wild Ancestors

As a consequence of artificial selection for specific traits, crop plants underwent considerable genotypic and phenotypic changes during the process of domestication. These changes may have led to reduced resistance in the cultivated plant due to shifts in resource allocation from defensive traits to increased growth rates and yield. Modern maize (Zea mays ssp. mays) was domesticated from its ancestor Balsas teosinte (Z. mays ssp. parviglumis) approximately 9000 years ago. Although maize displays a high genetic overlap with its direct ancestor and other annual teosintes, several studies show that maize and its ancestors differ in their resistance phenotypes with teosintes being less susceptible to herbivore damage. However, the underlying mechanisms are poorly understood. Here we addressed the question to what extent maize domestication has affected two crucial chemical and one physical defence traits and whether differences in their expression may explain the differences in herbivore resistance levels. The ontogenetic trajectories of 1,4-benzoxazin-3-ones, maysin and leaf toughness were monitored for different leaf types across several maize cultivars and teosinte accessions during early vegetative growth stages. We found significant quantitative and qualitative differences in 1,4-benzoxazin-3-one accumulation in an initial pairwise comparison, but we did not find consistent differences between wild and cultivated genotypes during a more thorough examination employing several cultivars/accessions. Yet, 1,4-benzoxazin-3-one levels tended to decline more rapidly with plant age in the modern maize cultivars. Foliar maysin levels and leaf toughness increased with plant age in a leaf-specific manner, but were also unaffected by domestication. Based on our findings we suggest that defence traits other than the ones that were investigated are responsible for the observed differences in herbivore resistance between teosinte and maize. Furthermore, our results indicate that single pairwise comparisons may lead to false conclusions regarding the effects of domestication on defensive and possibly other traits.     

Construction and comparison of three reference‐quality genome assemblies for soybean

We report reference‐quality genome assemblies and annotations for two accessions of soybean (Glycine max) and for one accession of Glycine soja, the closest wild relative of G. max. The G. max assemblies provided are for widely used US cultivars: the northern line Williams 82 (Wm82) and the southern line Lee. The Wm82 assembly improves the prior published assembly, and the Lee and G. soja assemblies are new for these accessions. Comparisons among the three accessions show generally high structural conservation, but nucleotide difference of 1.7 single‐nucleotide polymorphisms (snps) per kb between Wm82 and Lee, and 4.7 snps per kb between these lines and G. soja. snp distributions and comparisons with genotypes of the Lee and Wm82 parents highlight patterns of introgression and haplotype structure. Comparisons against the US germplasm collection show placement of the sequenced accessions relative to global soybean diversity. Analysis of a pan‐gene collection shows generally high conservation, with variation occurring primarily in genomically clustered gene families. We found approximately 40–42 inversions per chromosome between either Lee or Wm82v4 and G. soja, and approximately 32 inversions per chromosome between Wm82 and Lee. We also investigated five domestication loci. For each locus, we found two different alleles with functional differences between G. soja and the two domesticated accessions. The genome assemblies for multiple cultivated accessions and for the closest wild ancestor of soybean provides a valuable set of resources for identifying causal variants that underlie traits for the domestication and improvement of soybean, serving as a basis for future research and crop improvement efforts for this important crop species.     

Largely unlinked gene sets targeted by selection for domestication syndrome phenotypes in maize and sorghum

The domestication of diverse grain crops from wild grasses was a result of artificial selection for a suite of overlapping traits producing changes referred to in aggregate as ‘domestication syndrome’. Parallel phenotypic change can be accomplished by either selection on orthologous genes or selection on non‐orthologous genes with parallel phenotypic effects. To determine how often artificial selection for domestication traits in the grasses targeted orthologous genes, we employed resequencing data from wild and domesticated accessions of Zea (maize) and Sorghum (sorghum). Many ‘classic’ domestication genes identified through quantitative trait locus mapping in populations resulting from wild/domesticated crosses indeed show signatures of parallel selection in both maize and sorghum. However, the overall number of genes showing signatures of parallel selection in both species is not significantly different from that expected by chance. This suggests that while a small number of genes will extremely large phenotypic effects have been targeted repeatedly by artificial selection during domestication, the optimization part of domestication targeted small and largely non‐overlapping subsets of all possible genes which could produce equivalent phenotypic alterations.     

Higher expression of induced defenses in teosintes (Zea spp.) is correlated with greater resistance to fall armyworm,Spodoptera frugiperda

Selection for plant traits important for agriculture can come at a high cost to plant defenses. While selecting for increased growth rate and yield, domestication and subsequent breeding may lead to weakened defenses and greater susceptibility of plants to herbivores. We tested whether expression of defense genes differed among maize, Zea mays ssp. mays L. (Poaceae), and its wild relatives Zea mays ssp. parviglumis Iltis & Doebley and Zea diploperennis Iltis et al. We used two populations of Z. mays ssp. parviglumis: one expected to express high levels of an herbivore resistance gene, wound‐inducible protein (wip1), and another expected to have low expression of wip1. To test whether maize and wild Zea differed in induction of defenses against Spodoptera frugiperda (Smith) (Lepidoptera: Noctuidae), we quantified expression of several genes involved in plant defense: wip1, maize protease inhibitor (mpi), pathogenesis‐related protein (PR‐1), and chitinase. Moreover, we compared growth, development, and survival of caterpillars on maize and wild Zea plants. We found that maize expressed low levels of all but one of the genes when attacked by caterpillars, whereas the wild relatives of maize expressed induced defense genes at high levels. Expression of wip1, in particular, was much greater in the Z. mays ssp. parviglumis population that we expected to naturally express high levels of wip1, with expression levels 29‐fold higher than in herbivore‐free plants. Elevated expression of defenses in wild plants was correlated with higher resistance to caterpillars. Larvae were 15–20% smaller on wild Zea compared with maize, developed 20% slower, and only 22% of them survived to pupation on Z. mays ssp. parviglumis with high levels of wip1. Our results suggest that domestication has inadvertently reduced the resistance of maize, and it is likely that expression of wip1 and other genes associated with defenses play an important role in this reduction in resistance.     

Signatures of local adaptation in lowland and highland teosintes from whole‐genome sequencing of pooled samples

Spatially varying selection triggers differential adaptation of local populations. Here, we mined the determinants of local adaptation at the genomewide scale in the two closest maize wild relatives, the teosintes Zea mays ssp parviglumis and ssp. mexicana. We sequenced 120 individuals from six populations: two lowland, two intermediate and two highland populations sampled along two altitudinal gradients. We detected 8 479 581 single nucleotide polymorphisms (SNPs) covered in the six populations with an average sequencing depth per site per population ranging from 17.0× to 32.2×. Population diversity varied from 0.10 to 0.15, and linkage disequilibrium decayed very rapidly. We combined two differentiation‐based methods, and correlation of allele frequencies with environmental variables to detect outlier SNPs. Outlier SNPs displayed significant clustering. From clusters, we identified 47 candidate regions. We further modified a haplotype‐based method to incorporate genotype uncertainties in haplotype calling, and applied it to candidate regions. We retrieved evidence for selection at the haplotype level in 53% of our candidate regions, and in 70% of the cases the same haplotype was selected in the two lowland or the two highland populations. We recovered a candidate region located within a previously characterized inversion on chromosome 1. We found evidence of a soft sweep at a locus involved in leaf macrohair variation. Finally, our results revealed frequent colocalization between our candidate regions and loci involved in the variation of traits associated with plant–soil interactions such as root morphology, aluminium and low phosphorus tolerance. Soil therefore appears to be a major driver of local adaptation in teosintes.