Inference of domestication history and differentiation between early‐ and late‐flowering varieties in pearl millet

Pearl millet (Pennisetum glaucum) is a staple crop in Sahelian Africa. Farmers usually grow varieties with different cycle lengths and complementary functions in Sahelian agrosystems. Both the level of genetic differentiation of these varieties and the domestication history of pearl millet have been poorly studied. We investigated the neutral genetic diversity and population genetic structure of early‐ and late‐flowering domesticated and wild pearl millet populations using 18 microsatellite loci and 8 nucleotide sequences. Strikingly, early‐ and late‐flowering domesticated varieties were not differentiated over their whole distribution area, despite a clear difference in their isolation‐by‐distance pattern. Conversely, our data brought evidence for two well‐differentiated genetic pools in wild pearl millet, allowing us to test scenarios with different numbers and origins of domestication using approximate Bayesian computation (ABC). The ABC analysis showed the likely existence of asymmetric migration between wild and domesticated populations. The model choice procedure indicated that a single domestication from the eastern wild populations was the more likely scenario to explain the polymorphism patterns observed in cultivated pearl millet.     

Tuareg, a novel miniature-inverted repeat family of pearl millet (Pennisetum glaucum) related to the PIF superfamily of maize

Miniature-inverted repeat transposable elements (MITEs) are abundantly repeated in plant genomes and are especially found in genic regions where they could contribute regulatory elements for gene expression. We describe with molecular and cytological tools the first MITE family reported in pearl millet: Tuareg. It was initially detected in the pearl millet ortholog of Teosinte-branched1, an important developmental gene involved in the domestication of maize. The Tuareg family was amplified recently in the pearl millet genome and elements were found more abundant in wild than in domesticated plants. We found that they shared similarity in their terminal repeats with the previously described mPIF MITEs and that they are also present in other Pennisetum species, in maize and more distantly related grasses. The Tuareg family may be part of MITEs activated by PIF-like transposases and it could have been mobile since pearl millet domestication. Electronic supplementary material Electronic supplementary material is available for this article at and accessible for authorised users. O. Robin contributed the FISH and fiber-FISH hybridizations.     

Signatures of de‐domestication in autochthonous pig breeds and of domestication in wild boar populations from MC1R and NR6A1 allele distribution

Autochthonous pig breeds are usually reared in extensive or semi‐extensive production systems that might facilitate contact with wild boars and, thus, reciprocal genetic exchanges. In this study, we analysed variants in the melanocortin 1 receptor (MC1R) gene (which cause different coat colour phenotypes) and in the nuclear receptor subfamily 6 group A member 1 (NR6A1) gene (associated with increased vertebral number) in 712 pigs of 12 local pig breeds raised in Italy (Apulo‐Calabrese, Casertana, Cinta Senese, Mora Romagnola, Nero Siciliano and Sarda) and south‐eastern European countries (Kr?kopolje from Slovenia, Black Slavonian and Turopolje from Croatia, Mangalitsa and Moravka from Serbia and East Balkan Swine from Bulgaria) and compared the data with the genetic variability at these loci investigated in 229 wild boars from populations spread in the same macro‐geographic areas. None of the autochthonous pig breeds or wild boar populations were fixed for one allele at both loci. Domestic and wild‐type alleles at these two genes were present in both domestic and wild populations. Findings of the distribution of MC1R alleles might be useful for tracing back the complex genetic history of autochthonous breeds. Altogether, these results indirectly demonstrate that bidirectional introgression of wild and domestic alleles is derived and affected by the human and naturally driven evolutionary forces that are shaping the Sus scrofa genome: autochthonous breeds are experiencing a sort of ‘de‐domestication’ process, and wild resources are challenged by a ‘domestication’ drift. Both need to be further investigated and managed.     

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.     

A draft genome of sweet cherry (Prunus avium L.) reveals genome‐wide and local effects of domestication

Sweet cherry (Prunus avium L.) trees are both economically important fruit crops but also important components of natural forest ecosystems in Europe, Asia and Africa. Wild and domesticated trees currently coexist in the same geographic areas with important questions arising on their historical relationships. Little is known about the effects of the domestication process on the evolution of the sweet cherry genome. We assembled and annotated the genome of the cultivated variety “Big Star*” and assessed the genetic diversity among 97 sweet cherry accessions representing three different stages in the domestication and breeding process (wild trees, landraces and modern varieties). The genetic diversity analysis revealed significant genome‐wide losses of variation among the three stages and supports a clear distinction between wild and domesticated trees, with only limited gene flow being detected between wild trees and domesticated landraces. We identified 11 domestication sweeps and five breeding sweeps covering, respectively, 11.0 and 2.4 Mb of the P. avium genome. A considerable fraction of the domestication sweeps overlaps with those detected in the related species, Prunus persica (peach), indicating that artificial selection during domestication may have acted independently on the same regions and genes in the two species. We detected 104 candidate genes in sweep regions involved in different processes, such as the determination of fruit texture, the regulation of flowering and fruit ripening and the resistance to pathogens. The signatures of selection identified will enable future evolutionary studies and provide a valuable resource for genetic improvement and conservation programs in sweet cherry.     

High degree of genetic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (F_ST = 0.033 or 0.011, respectively), whereas average F_ST for directionally selected markers was significantly higher (F_ST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.     

Y‐chromosome and mtDNA variation confirms independent domestications and directional hybridization in South American camelids

Investigations of genetic diversity and domestication in South American camelids (SAC) have relied on autosomal microsatellite and maternally‐inherited mitochondrial data. We present the first integrated analysis of domestic and wild SAC combining male and female sex‐specific markers (male specific Y‐chromosome and female‐specific mtDNA sequence variation) to assess: (i) hypotheses about the origin of domestic camelids, (ii) directionality of introgression among domestic and/or wild taxa as evidence of hybridization and (iii) currently recognized subspecies patterns. Three male‐specific Y‐chromosome markers and control region sequences of mitochondrial DNA are studied here. Although no sequence variation was found in SRY and ZFY, there were seven variable sites in DBY generating five haplotypes on the Y‐chromosome. The haplotype network showed clear separation between haplogroups of guanaco–llama and vicuña–alpaca, indicating two genetically distinct patrilineages with near absence of shared haplotypes between guanacos and vicuñas. Although we document some examples of directional hybridization, the patterns strongly support the hypothesis that llama (Lama glama) is derived from guanaco (Lama guanicoe) and the alpaca (Vicugna pacos) from vicuña (Vicugna vicugna). Within male guanacos we identified a haplogroup formed by three haplotypes with different geographical distributions, the northernmost of which (Peru and northern Chile) was also observed in llamas, supporting the commonly held hypothesis that llamas were domesticated from the northernmost populations of guanacos (L. g. cacilensis). Southern guanacos shared the other two haplotypes. A second haplogroup, consisting of two haplotypes, was mostly present in vicuñas and alpacas. However, Y‐chromosome variation did not distinguish the two subspecies of vicuñas.     

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.     

Genomic variation in cline shape across a hybrid zone

Hybrid zones are unique biological interfaces that reveal both population level and species level evolutionary processes. A genome‐scale approach to assess gene flow across hybrid zones is vital, and now possible. In Mexican towhees (genus Pipilo), several morphological hybrid gradients exist. We completed a genome survey across one such gradient (9 populations, 140 birds) using mitochondrial DNA, 28 isozyme, and 377 AFLP markers. To assess variation in introgression among loci, cline parameters (i.e., width, center) for the 61 clinally varying loci were estimated and compiled into genomic distributions for tests against three empirical models spanning the range of observed cline shape. No single model accounts for observed variation in cline shape among loci. Numerous backcross individuals near the gradient center confirm a hybrid origin for these populations, contrary to a previous hypothesis based on social mimicry and character displacement. In addition, the observed variation does not bin into well‐defined categories of locus types (e.g., neutral vs. highly selected). Our multi‐locus analysis reveals cross‐genomic variation in selective constraints on gene flow and locus‐specific flexibility in the permeability of the interspecies membrane.     

The structure of wild and domesticated emmer wheat populations,gene flow between them,and the site of emmer domestication

The domestication of emmer wheat (Triticum turgidum spp. dicoccoides, genomes BBAA) was one of the key events during the emergence of agriculture in southwestern Asia, and was a prerequisite for the evolution of durum and common wheat. Single- and multilocus genotypes based on restriction fragment length polymorphism at 131 loci were analyzed to describe the structure of populations of wild and domesticated emmer and to generate a picture of emmer domestication and its subsequent diffusion across Asia, Europe and Africa. Wild emmer consists of two populations, southern and northern, each further subdivided. Domesticated emmer mirrors the geographic subdivision of wild emmer into the northern and southern populations and also shows an additional structure in both regions. Gene flow between wild and domesticated emmer occurred across the entire area of wild emmer distribution. Emmer was likely domesticated in the Diyarbakir region in southeastern Turkey, which was followed by subsequent hybridization and introgression from wild to domesticated emmer in southern Levant. A less likely scenario is that emmer was domesticated independently in the Diyarbakir region and southern Levant, and the Levantine genepool was absorbed into the genepool of domesticated emmer diffusing from southeastern Turkey. Durum wheat is closely related to domesticated emmer in the eastern Mediterranean and likely originated there. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The role of Bh4 in parallel evolution of hull colour in domesticated and weedy rice

The two independent domestication events in the genus Oryza that led to African and Asian rice offer an extremely useful system for studying the genetic basis of parallel evolution. This system is also characterized by parallel de‐domestication events, with two genetically distinct weedy rice biotypes in the US derived from the Asian domesticate. One important trait that has been altered by rice domestication and de‐domestication is hull colour. The wild progenitors of the two cultivated rice species have predominantly black‐coloured hulls, as does one of the two U.S. weed biotypes; both cultivated species and one of the US weedy biotypes are characterized by straw‐coloured hulls. Using Black hull 4 (Bh4) as a hull colour candidate gene, we examined DNA sequence variation at this locus to study the parallel evolution of hull colour variation in the domesticated and weedy rice system. We find that independent Bh4‐coding mutations have arisen in African and Asian rice that are correlated with the straw hull phenotype, suggesting that the same gene is responsible for parallel trait evolution. For the U.S. weeds, Bh4 haplotype sequences support current hypotheses on the phylogenetic relationship between the two biotypes and domesticated Asian rice; straw hull weeds are most similar to indica crops, and black hull weeds are most similar to aus crops. Tests for selection indicate that Asian crops and straw hull weeds deviate from neutrality at this gene, suggesting possible selection on Bh4 during both rice domestication and de‐domestication.     

Gene editing of the wheat homologs of TONNEAU1‐recruiting motif encoding gene affects grain shape and weight in wheat

Grain size and weight are important components of a suite of yield‐related traits in crops. Here, we showed that the CRISPR‐Cas9 gene editing of TaGW7, a homolog of rice OsGW7 encoding a TONNEAU1‐recruiting motif (TRM) protein, affects grain shape and weight in allohexaploid wheat. By editing the TaGW7 homoeologs in the B and D genomes, we showed that mutations in either of the two or both genomes increased the grain width and weight but reduced the grain length. The effect sizes of mutations in the TaGW7 gene homoeologs coincided with the relative levels of their expression in the B and D genomes. The effects of gene editing on grain morphology and weight traits were dosage dependent with the double‐copy mutant showing larger effect than the respective single copy mutants. The TaGW7‐centered gene co‐expression network indicated that this gene is involved in the pathways regulating cell division and organ growth, also confirmed by the cellular co‐localization of TaGW7 with α‐ and β‐tubulin proteins, the building blocks of microtubule arrays. The analyses of exome capture data in tetraploid domesticated and wild emmer, and hexaploid wheat revealed the loss of diversity around TaGW7‐associated with domestication selection, suggesting that TaGW7 is likely to play an important role in the evolution of yield component traits in wheat. Our study showed how integrating CRISPR‐Cas9 system with cross‐species comparison can help to uncover the function of a gene fixed in wheat for allelic variants targeted by domestication selection and select targets for engineering new gene variants for crop improvement.     

Footprints of selection in wild populations of Bicyclus anynana along a latitudinal cline

One of the major questions in ecology and evolutionary biology is how variation in the genome enables species to adapt to divergent environments. Here, we study footprints of thermal selection in candidate genes in six wild populations of the afrotropical butterfly Bicyclus anynana sampled along a c. 3000 km latitudinal cline. We sequenced coding regions of 31 selected genes with known functions in metabolism, pigment production, development and heat shock responses. These include genes for which we expect a priori a role in thermal adaptation and, thus, varying selection pressures along a latitudinal cline, and genes we do not expect to vary clinally and can be used as controls. We identified amino acid substitution polymorphisms in 13 genes and tested these for clinal variation by correlation analysis of allele frequencies with latitude. In addition, we used two F_ST‐based outlier methods to identify loci with higher population differentiation than expected under neutral evolution, while accounting for potentially confounding effects of population structure and demographic history. Two metabolic enzymes of the glycolytic pathway, UGP and Treh, showed clinal variation. The same loci showed elevated population differentiation and were identified as significant outliers. We found no evidence of clines in the pigmentation genes, heat shock proteins and developmental genes. However, we identified outlier loci in more localized parts of the range in the pigmentation genes yellow and black. We discuss that the observed clinal variation and elevated population divergence in UGP and Treh may reflect adaptation to a geographic thermal gradient.     

Domestication of tomato has reduced the attraction of herbivore natural enemies to pest‐damaged plants

相似文献   

Quantitative trait loci for growth and body size in the nine‐spined stickleback Pungitius pungitius L.

Body size is an ecologically important trait shown to be genetically variable both within and among different animal populations as revealed by quantitative genetic studies. However, few studies have looked into underlying genetic architecture of body size variability in the wild using genetic mapping methods. With the aid of quantitative trait loci (QTL) analyses based on 226 microsatellite markers, we mapped body size and growth rate traits in the nine‐spined stickleback (Pungitius pungitius) using an F₂‐intercross (n = 283 offspring) between size‐divergent populations. In total, 17 QTL locations were detected. The proportion of phenotypic variation explained by individual body size‐related QTL ranged from 3% to 12% and those related to growth parameters and increments from 3% to 10%. Several of the detected QTL affected either early or late growth. These results provide a solid starting point for more in depth investigations of structure and function of genomic regions involved in determination of body size in this popular model of ecological and evolutionary research.     

Genetic diversity of Chinese cattle revealed by Y‐SNP and Y‐STR markers

With its vast territory and complex natural environment, China boasts rich cattle genetic resources. To gain the further insight into the genetic diversity and paternal origins of Chinese cattle, we analyzed the polymorphism of Y‐SNPs (UTY19 and ZFY10) and Y‐STRs (INRA189 and BM861) in 34 Chinese cattle breeds/populations, including 606 males representative of 24 cattle breeds/populations collected in this study as well as previously published data for 302 bulls. Combined genotypic data identified 14 Y‐chromosome haplotypes that represented three haplogroups. Y2‐104‐158 and Y2‐102‐158 were the most common taurine haplotypes detected mainly in northern and central China, whereas the indicine haplotype Y3‐88‐156 predominates in southern China. Haplotypes Y2‐108‐158, Y2‐110‐158, Y2‐112‐158 and Y3‐92‐156 were private to Chinese cattle. The population structure revealed by multidimensional scaling analysis differentiated Tibetan cattle from the other three groups of cattle. Analysis of molecular variance showed that the majority of the genetic variation was explained by the genetic differences among groups. Overall, our study indicates that Chinese cattle retain high paternal diversity (H = 0.607 ± 0.016) and probably much of the original lineages that derived from the domestication center in the Near East without strong admixture from commercial cattle carrying Y1 haplotypes.     

SNP discovery in wild and domesticated populations of blue catfish,Ictalurus furcatus,using genotyping‐by‐sequencing and subsequent SNP validation

Blue catfish, Ictalurus furcatus, are valued in the United States as a trophy fishery for their capacity to reach large sizes, sometimes exceeding 45 kg. Additionally, blue catfish × channel catfish (I. punctatus) hybrid food fish production has recently increased the demand for blue catfish broodstock. However, there has been little study of the genetic impacts and interaction of farmed, introduced and stocked populations of blue catfish. We utilized genotyping‐by‐sequencing (GBS) to capture and genotype SNP markers on 190 individuals from five wild and domesticated populations (Mississippi River, Missouri, D&B, Rio Grande and Texas). Stringent filtering of SNP‐calling parameters resulted in 4275 SNP loci represented across all five populations. Population genetics and structure analyses revealed potential shared ancestry and admixture between populations. We utilized the Sequenom MassARRAY to validate two multiplex panels of SNPs selected from the GBS data. Selection criteria included SNPs shared between populations, SNPs specific to populations, number of reads per individual and number of individuals genotyped by GBS. Putative SNPs were validated in the discovery population and in two additional populations not used in the GBS analysis. A total of 64 SNPs were genotyped successfully in 191 individuals from nine populations. Our results should guide the development of highly informative, flexible genotyping multiplexes for blue catfish from the larger GBS SNP set as well as provide an example of a rapid, low‐cost approach to generate and genotype informative marker loci in aquatic species with minimal previous genetic information.     

Global investigation of the co‐evolution of MIRNA genes and microRNA targets during soybean domestication

Although the selection of coding genes during plant domestication has been well studied, the evolution of MIRNA genes (MIRs) and the interaction between microRNAs (miRNAs) and their targets in this process are poorly understood. Here, we present a genome‐wide survey of the selection of MIRs and miRNA targets during soybean domestication and improvement. Our results suggest that, overall, MIRs have higher evolutionary rates than miRNA targets. Nonetheless, they do demonstrate certain similar evolutionary patterns during soybean domestication: MIRs and miRNA targets with high expression and duplication status, and with greater numbers of partners, exhibit lower nucleotide divergence than their counterparts without these characteristics, suggesting that expression level, duplication status, and miRNA–target interaction are essential for evolution of MIRs and miRNA targets. Further investigation revealed that miRNA–target pairs that are subjected to strong purifying selection have greater similarities than those that exhibited genetic diversity. Moreover, mediated by domestication and improvement, the similarities of a large number of miRNA–target pairs in cultivated soybean populations were increased compared to those in wild soybeans, whereas a small number of miRNA–target pairs exhibited decreased similarity, which may be associated with the adoption of particular domestication traits. Taken together, our results shed light on the co‐evolution of MIRs and miRNA targets during soybean domestication.     

Elucidation of the origin of ‘agriocrithon’ based on domestication genes questions the hypothesis that Tibet is one of the centers of barley domestication

Wild barley forms a two‐rowed spike with a brittle rachis whereas domesticated barley has two‐ or six‐rowed spikes with a tough rachis. Like domesticated barley, ‘agriocrithon’ forms a six‐rowed spike; however, the spike is brittle as in wild barley, which makes the origin of agriocrithon obscure. Haplotype analysis of the Six‐rowed spike 1 (vrs1) and Non‐brittle rachis 1 (btr1) and 2 (btr2) genes was conducted to infer the origin of agriocrithon barley. Some agriocrithon barley accessions (eu‐agriocrithon) carried Btr1 and Btr2 haplotypes that are not found in any cultivars, implying that they are directly derived from wild barley through a mutation at the vrs1 locus. Other agriocrithon barley accessions (pseudo‐agriocrithon) carried Btr1 or Btr2 from cultivated barley, thus implying that they originated from hybridization between six‐rowed landraces carrying btr1Btr2 and Btr1btr2 genotypes followed by recombination to produce Btr1Btr2. All materials we collected from Tibet belong to pseudo‐agriocrithon and thus do not support the Tibetan Plateau as being a center of barley domestication. Tracing the evolutionary history of these allelic variants revealed that eu‐agriocrithon represents six‐rowed barley lineages that were selected by early farmers, once in south‐eastern Turkmenistan (vrs1.a1) and again in the eastern part of Uzbekistan (vrs1.a4).     

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.