Chloroplast DNA non-coding sequences variation in Aegilops tauschii Coss.: evolutionary history of the species

Sequences of four chloroplast DNA non-coding regions, about 3,000?bp in total, were analysed in 112 Aegilops tauschii accessions, 56 of ssp. tauschii and 56 of ssp. strangulata, representing all of the species range. One inversion, 8 insertions/deletions, 18 base pair substitutions and 5 microsatellite loci were found. The data revealed that Ae. tauschii originated in Caucasia. Neither of the two Ae. tauschii subspecies was an ancestor to one another. Aegilops tauschii divided into ssp. tauschii and ssp. strangulata at the very beginning of its existence as a species. Subspecies tauschii was the first to start geographic expansion and relatively rapidly occupied a vast area from Caucasia—eastward up to central Tien Shan and western Himalayas. In contrast to ssp. tauschii, geographic spread of ssp. strangulata was a complicated, multi-stage and slow process. At the beginning of ssp. strangulata evolutionary history its major phylogenetic lineage for a lengthy time span had existed as a small isolated population. Several forms of ssp. strangulata, better adapted to relatively moister and cooler habitats, had originated. Each of these forms has gradually forced out ssp. tauschii from some part of its area in the west, up to central Kopet-Dag.     

Distribution and diversity of Aegilops tauschii in Iran

The wide morphological variation of Aegilops tauschii has led to the distinction of different subspecies; a typical ssp. tauschii and a second ssp. strangulata. However some researchers pointed out the existance of the intermediate form among morphologically distinguished subspecies. Distribution, diversity and the relationship between different subspecies and the intermediate form were evaluated in the Iranian Ae. tauschii collection. This collection was classified to 15 different populations according to morphological similarities and the collecting origin of accessions. The highest variation was found in tauschii population of Golestan followed by tauschii populations of Gilan and Ardebill, whereas the lowest variation was observed in tauschii populations of central Iran. Two discriminant functions suggested that the length of rachis node and spikelet glume, particularly, the length/width ratios of these traits had the highest impact on identification of different forms. Mahalanobis distances (D ² ) between the two subspecies along with intermediate form on the multidimensional scaling plot showed that the intermediate form is more similar to ssp. tauschii than ssp. strangulata. Although, the diversity within the ssp. strangulata was not very high, it widely affected the diversity of Iranian accessions of Ae. tauschii through continues crossing with the more diversed subspecies, tauschii, during thousands of years. This fact had lead to expansion of its distribution from its origin to Northern Khorasan, Northern Semnan and Eastern Ardebill by producing the intermediate form.     

Stripe rust resistance in <Emphasis Type="Italic">Aegilops tauschii</Emphasis> and its genetic analysis

Aegilops tauschii Coss., the D-genome progenitor of common wheat (Triticum aestivum L.) includes two subspecies, tauschii and strangulata (Eig) Tzvel. Subspecies tauschii has a wide geographic distribution spreading westwards to Turkey and eastwards to Afghanistan and China, while ssp. strangulata has a narrower distribution occurring only in two disjoined regions, southeastern Caspian Iran and Transcaucasia. A collection of 56 Ae. tauschii accessions was screened at adult stage against a mixture of pathotypes of stripe rust prevalent in the current wheat production in China. The results for three crop seasons indicated that among the 38 ssp. tauschii accessions, 37 were susceptible and only one was resistant, while all the 18 ssp. strangulata accessions were resistant. These results indicated that stripe rust resistance was related to taxonomic origin. Further genetic analysis revealed the resistance of stripe rust in ssp. strangulata accession AS2388 was conferred by a single dominant gene.     

Allelic diversity of high molecular weight glutenin subunits (HMW-GS) in Iranian Aegilops tauschii Coss. accessions by sodium dodecyl sulphate polyacrylamide gel electrophoresis (SDS-PAGE)

Variation of high molecular weight glutenin subunits (HMW-GS) in 28 Iranian Aegilops tauschii (2n = 2x = 14, DD) accessions studied by sodium dodecyl sulphate electrophoresis method (SDS-PAGE). The results showed high variation of HMW-GS in the accessions. The range of frequency in 14 HMW-GS combinations was 3.57–25 % in the accessions. AMOVA showed the molecular variance between the geographic areas was lower than within the geographic areas. According to Nei’s genetic diversity, the highest diversity levels were in Semnan, Golestan and Azarbayjan, on the other hand the lowest levels of diversity were found in Khorasan, Gilan and Mazandaran accessions. Hence, the Caspian Sea South East accessions also Azerbayjan in Iran have more diversity. AMOVA did not show variance between strangulata and tauschii but there was more genetic diversity in ssp. tauschii subspecies in comparison of ssp. strangulata according to Nei’s gene diversity and Shannon information index. It showed Iranian Ae. tauschii have a good potential for bread making quality improvement in bread wheat.     

Geographic patterns of Got1, Got2, Got3 and Est2 enzyme-encoding genes allelic variation in Aegilops tauschii Coss.

Geographic patterns of Got1, Got2, Got3 and Est2 enzyme-encoding genes allelic variation were investigated among 322 accessions of Aegilops tauschii Coss., 161 accessions of ssp. tauschii and 161 accessions of ssp. strangulata, representing all the species area. It was found that: (1) in the two ecologically different subspecies, ssp. tauschii and ssp. strangulata, the patterns of allelic variation of the four genes differ greatly; (2) the same allozymes have originated several times independently in different Ae. tauschii local populations; (3) allelic variation of Got1, Got3 and Est2 in ssp. strangulata corresponds to climatic conditions. The data obtained reflected that Ae. tauschii has been inhabiting its area from ancient times; and allelic variation patterns of Got1, Got2, Got3, Est2 loci were mostly formed by natural selection. Further investigations of these loci with molecular genetic methods are prospective for understanding the peculiarities of Ae. tauschii evolution.     

Microsatellite analysis of genetic diversity and population genetic structure of <Emphasis Type="Italic">Aegilops tauschii</Emphasis> Coss. in Northern Iran

Genetic diversity and population genetic structure of Aegilops tauschii in Northern Iran were studied based on nine microsatellite loci. A high level of genetic diversity was observed from the accessions collected from six regions (provinces). These accessions include 79 samples of the two subspecies (tauschii and strangulata), the intermediate form (among morphologically distinguished subspecies) and ten accessions of Triticum aestivum. The nine microsatellites revealed a total of 141 alleles, with an average of 15.7 alleles per locus. A comparison of the parameters showing genetic diversity, including the observed heterozygosity (Ho), gene diversity (He) and Shannon’s information index (I) of Ae. tauschii accessions from different provinces in Northern Iran, indicated that subsp. tauschii possesses the highest genetic diversity, followed by intermediate form. Genetic distance between subsp. strangulata and subsp. tauschii was low, confirming high gene flow between these two subspecies. However, intermediate form was more distinct from both of them. It was also found that the genetic diversity of T. aestivum is obviously lower than that of Ae. tauschii accessions. Moreover, the level of genetic diversity for Gilan, Golestan and Mazanderan provinces was higher than for Ardebil, Ghazvin and Semnan provinces, suggesting that these regions may provide a readily available source of potentially useful variation for wheat improvement.     

Novel and ancient HMW glutenin genes from Aegilops tauschii and their phylogenetic positions

A pair of novel high-molecular-weight glutenin subunits (HMW-GS) 1Dx3.1^t and 1Dy11*^t were revealed and characterized from Aegilops tauschii Coss. subspecies tauschii accession AS60. SDS-PAGE band of 1Dx3.1^t was between those of 1Dx2 and 1Dx3, while 1Dy11*^t was between 1Dy11 and 1Dy12. The lengths of 1Dx3.1 ^t and 1Dy11* ^t were 2,514?bp and 1,968?bp, encoding 836 and 654 amino acid residues, respectively. Their authenticity was confirmed by successful expression of the coding regions in Escherichia coli. Network analysis indicated that 1Dx3.1 ^t together with other five rare alleles only detected in Asia common wheat populations represented the ancestral sequences in Glu-D1 locus. Neighbor-joining tree analysis of previously cloned x-type and y-type alleles in the Glu-D1 locus supported the hypothesis that more than one Ae. tauschii genotypes were involved in the origin of hexaploid wheat and that different Ae. tauschii accessions contributed the D genome to common wheat and Ae. cylindrical Host, respectively. An Ae. tauschii accession with 1Dx3.1 ^t or a closely related allele probably have involved in the origin of common wheat. Since accession AS60 used in this study belonged to typical ssp. tauschii, present results suggested the possibility that ssp. tauschii was involved in the evolution of common wheat.     

Spatial patterns of adenylate kinase,catalase, endopeptidase and fructose-1,6-diphosphatase encoding genes allelic variation in <Emphasis Type="Italic">Aegilops tauschii</Emphasis> Coss.

Investigation of spatial patterns of adenylate kinase, catalase, endopeptidase and fructose-1,6-diphosphatase encoding genes (Ak, Cat1, Cat2, Ep, Fdp) allelic variation in Aegilops tauschii was carried out. About 300 accessions, representing all the species range were taken for the study. Cat2 and Fdp loci are completely monomorphic in ssp. strangulata and in the western part of ssp. tauschii range, as well. Both Cat2 and Fdp are highly polymorphic in the eastern part of ssp. tauschii range, with the patterns of this polymorphism being discordant in these two loci. Ak ¹⁰⁸, a rare allele with sporadical spatial occurrence, was found in ssp. tauschii only. On the contrary, Ak ⁹² is absent in ssp. tauschii: it is the most common Ak allele in ssp. strangulata in Precaspian Iran, the most moist part of the area, and is very rare in other parts of ssp. strangulata area. Ep is a highly polymorphic locus with the highest level of variation in the west of Ae. tauschii area, where this species had originated. Ep allele variation patterns are rather similar in ssp. tauschii and ssp. strangulata. The data reveal the adaptive nature of Ak, Cat2, and Fdp allele variation, while Ep polymorphism seems to be mostly neutral.     

Aegilops tauschii: Genetic Variation in Iran

All the 79 Aegilops tauschii Coss. accessions of Iranian origin from Prof. Kihara’s collection were analyzed electrophoretically. Of 23 enzyme-encoding loci studied, 11 were polymorphic. In Iran Ae. tauschii is presented by ssp. tauschii and ssp. strangulata which distinctly differ genetically, morphologically and ecologically. Variation patterns of low polymorphic locus Aco2 and highly polymorphic Ep are similar in both subspecies. In contrast, variation of Acph1, Ak, Est2, Est5, Got1, Got2, Got3 and Lap is a set of diverse patterns which markedly differ between subspecies and natural regions also, implying that natural selection is involved.     

Multivariate analysis of genetic variation in Aegilops tauschii from the world germplasm collection

A study of Aegilops tauschii subspecies constitution was undertaken. The data on allozyme and morphologic variation among 308 plants from 154 accessions were used for multivariate analysis. ACPH1 and (glume width)/(rachis segment width) ratio were found to be reliable criteria to distinguish between sspp. tauschii and strangulata.     

Polymorphism of <Emphasis Type="Italic">Got2</Emphasis> DNA sequences sheds light on <Emphasis Type="Italic">Aegilops tauschii</Emphasis> Coss. intraspecies divergence and origin of <Emphasis Type="Italic">Triticum aestivum</Emphasis> L.

DNA sequences of nuclear gene Got2 was studied in 60 accessions of Aegilops tauschii, 29 of subsp. tauschii and 31 of subsp. strangulata. It was found that Got2 allozyme polymorphism in Ae. tauschii is due to a single, unique, mutation which led to replacement of glutamic acid by isoleucine in residue 256 of the enzyme molecule, encoded by Got2. As revealed by Got2 DNA sequences variation, initially in its history Ae. tauschii was presented by subsp. strangulata, and among phylogenetic lineages of subsp. strangulata, the lineage “t-9¹s” (TauL3) is the most ancient, a relict one. Subspecies tauschii is relatively “young”. Initially it was presented by the lineage marked by combination of allozyme alleles Got2 ¹⁰⁵ and Acph1 ¹⁰⁰. In the past it inhabited the Continental area from Caucasia to Pakistan, but later on it was forced out by newly originated, now—a major lineage of subsp. tauschii, marked by Got2 ¹⁰⁰. This lineage extended the Continental area of the species up to Kirgizstan, but actually failed to penetrate into pre-Caspian area, occupied by subsp. strangulata. These results essentially differ from those obtained previously, using chloroplast DNA (cpDNA) sequences polymorphism. As revealed by cpDNA, the major, “usual”, subsp. strangulata (TauL2) is “younger” than subsp. tauschii, which resided on phylogenetic tree between relict lineage “t-9¹s”of subsp. strangulata—and major subsp. strangulata. But both cpDNA and Got2 DNA sequences indicate that the level of genetic variation in subsp. tauschii is much lower than in subsp. strangulata. According to Got2 DNA sequences variation, it was Ae. tauschii subsp. strangulata lineage “k-109″ which donated genome D to Triticum aestivum L. This lineage includes accessions: k-109 from South-Eastern Precaspian Azerbaijan; KU-2105, KU-2159 from Western Precaspian Iran; KU-2080 from Eastern Precaspian Iran.     

Polymorphism of high molecular weight glutenin subunits in three species of Aegilops

A collection of 136 accessions of Aegilops umbellulata (39), Ae. comosa (75) and Ae. markgrafii (22) was analysed for high-molecular-weight (HMW) glutenin subunits composition. The homogeneity of the accessions was studied and 55.1% of the collection was homogeneous for HMW glutenin subunits (29 Ae. umbellulata, 33 Ae. comosa and 14 Ae. markgrafii). The HMW glutenin subunits of Ae. umbellulata are encoded by the Glu-U1 locus; in Ae. comosa results showed that this proteins are encoded at the 1M chromosome, and the locus was named Glu-M1. In Ae. markgrafii it was assumed that HMW glutenin subunits were encoded by an homoeologous locus and it was named Glu-C1. All the accessions of Ae. umbellulata and Ae. markgrafii expressed both, x-type and y-type subunits. Among the Ae. comosa accessions, only one expressed an x-type subunit alone. All the accessions of Ae. umbellulata and some of Ae. comosa had x-type glutenins of higher molecular weights than those commonly present in bread wheat. A total of 8 alleles were detected at the Glu-U1 locus, 11 at the Glu-M1 and 4 at the Glu-C1. The new HMW glutenin variation found in this work suggests their possible utilisation in breeding for wheat quality.     

Aegilops tauschii Coss.: allelic variation of enzyme-encoding genes and ecological differentiation of the species

Three hundred and seven accessions of Aegilops tauschii Coss., including 160 of subsp. tauschii and 147 of subsp. strangulata, representing all the species range—from Turkey to Kirgizstan, were analyzed electrophoretically. Twenty polymorphic enzyme-encoding loci were studied, 10 of which were essentially polymorphic in Ae. tauschii. Climatic data for each of the 307 Ae. tauschii habitats were taken from WORLDCLIM database of computer system ArcGIS. Forty-nine climatic parameters were considered: precipitation, minimum, mean and maximum temperatures for each month, and also the total annual level of precipitation. The data were analyzed with multivariate statistical methods, such as Principal Components Analysis (PCA), Multiple Correspondence Analysis (MCA) and Two-Block Partial Least Squares. Variability of climatic conditions among Ae. tauschii habitats is reflected by the two approximately orthogonal “vectors”. The “first vector” is mostly determined by negative impact of precipitation and minimum temperatures during winter. The “second vector” is mostly determined by negative impact of maximum temperatures during summer, and positive impact of precipitation during late spring and summer. Aegilops tauschii is essentially variable along the “second vector”, and especially high level of variation is characteristic for subsp. tauschii. This variation reflects that Ae. tauschii is very tolerable to the climatic variation during summer season. Aegilops tauschii subsp. strangulata is also characterized by the high level of variation along the “first vector”. Moreover, all the habitats of subsp. strangulata fall into the two distinct separate clusters: the habitats in Precaspian Iran, which have the highest minimum temperatures in winter,—and all the other habitats. In the plot of the first two factors of PCA, the “cluster of Precaspian Iran” can be further divided into “Western Precaspian Iran (WPI)”, having relatively higher level of annual rainfall, and relatively dryer “Eastern Precaspian Iran (EPI)”. This three groups of subsp. strangulata accessions, from WPI, EPI and other areas, are also distinctly differed in enzyme-encoding genes allelic variation, as revealed on the plot of the first two axes of MCA. In contrast to subsp. strangulata, the level of variation of subsp. tauschii along the “first vector” is rather low. It was pointed out that variation along “the first vector” reflects adaptive intraspecies divergence of Ae. tauschii: its subspecies strangulata “prefers” the habitats of seaside climate, with warm and moist winter; while subsp. tauschii mostly occupies the habitats with rather continental climate, with relatively cold and dry winter. Allelic variation of enzyme-encoding genes Acph1, Ak, Est2, Est5, Got1, Got2, and Got3 correlate with climate along “the first vector”. Apparently, polymorphism of these loci were involved into the process of Ae. tauschii intraspecies adaptive divergence. Allelic variation of Cat2 and Fdp loci correspond to climatic variation along “the second vector” in subsp. tauschii. Therefore Cat2 and Fdp are likely to be among the genes which polymorphism “helped” subsp. tauschii to succeed in vast geographical expansion far to the east from Caspian Sea.     

Polymorphism of gliadins in <Emphasis Type="Italic">Aegilops tauschii</Emphasis> Coss. local populations in two primary habitats in Dagestan

Polymorphism of gliadins was investigated in Aegilops tauschii from primary habitats “4”, near Hily, and “6”, near Rukel, in Dagestan, Russia 205 individual plants were analysed (53/50 and 54/48 plants of subsp. tauschii/subsp. strangulata from the habitats “4” and “6”, respectively) and 1/7 and 18/14 different haplotypes were found among the plants of subsp. tauschii/subsp. strangulata from the habitats “4” and “6”, respectively. No direct evidences of cross-pollination were pointed out, although gliadins electrophoretic phenotypes obtained allowed to suggest that it occur in Ae. tauschii with very low frequency. The data obtained revealed that during Ae. tauschii evolutionary history a local habitat could be populated many times by different phylogenetic lineages of the species. It was found that in Dagestan, at the very edge of the species area, several different lineages belonging to different subspecies could for a long time co-exist together in a local habit, and in such case a very high level of genetic variation in Ae. tauschii could be accumulated on a square of less than one hectare. The further studies of genetic variation in Ae. tauschii local populations, based on molecular genetic methods seems to be very prospective for understanding of peculiarities of the species evolution.     

Purification and Characterization of the Glutenin Subunits of Triticum tauschii,Progenitor of the D Genome in Hexaploid Bread Wheat

N-terminal amino acid sequences and sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) molecular weights have been determined for high-performance liquid chromatography (HPLC)-purified high molecular weight (HMW) and low molecular weight (LMW) glutenin subunits (GS) of Triticum tauschii ssp. strangulata, contributor of the D genome to hexaploid bread wheat. The use of three different extraction procedures resulted in similar glutenin preparations. On the basis of N-terminal sequences, the same types of glutenin subunits that have been reported in bread and durum wheats (HMW-GS of both the x and y types and LMW-GS of the LMW-s, LMW-m, α-, and γ-types) were found in T. tauschii. However, the HMW-GS in T. tauschii were in greater proportion relative to LMW-GS when compared to reported values for a bread and durum wheat. Our results support the likelihood that differences in the proportions of the various subunits contributed by the A, B, and D genomes, rather than qualitative differences in the types of subunits, are responsible for the major differences in quality characteristics between bread wheat and durum wheat.     

Evolution of tetraploid wheat based on variations in 5′ UTR regions of Ppd-A1: evidence of gene flow between emmer and timopheevi wheat

Previous study showed that tetraploid wheat was divided into two groups (Type AI and Type AII) based on sequences around Ppd-A1 gene (Takenaka and Kawahara in Theor Appl Genet 125(5):999–1014, 2012). That study focused on domesticated emmer wheat and used only 19 wild emmer wheats, so could not be clear the evolutional relationship between Type AI and Type AII. Here, a total of 669 accessions comprising 65 einkorn wheats, 185 wild emmer wheats, 107 hulled emmer wheats, 204 free-threshing (FT) emmer wheats, and 108 timopheevii wheats were studied by PCR assay and DNA sequencing for Type AI/AII. Type AII was an older type than Type AI because all einkorn accessions had Type AII. In wild emmer, Type AI was distributed in the northeast regions of its distribution and Type AII was found to be centered on Israel. A total of 37.4 % of hulled emmer accessions were Type AI, while 92.2 % of FT emmer accessions were Type AI. Differences in the proportion of Type AI/AII in domesticated emmer suggested a strong bottle-neck effect. We also found two MITE-like sequence deletion patterns from a part of Type AII accessions (dic-del and ara-del). Dic-del was found from only Israeli wild emmer accessions and ara-del was found from almost all timopheevii wheat accessions. Only three timopheevii accessions did not have ara-del, and one wild emmer accession and ten hulled emmer accessions had ara-del. These accessions suggested gene flow between emmer and timopheevii wheat.     

Re-sequencing of vrs1 and int-c loci shows that labile barleys (Hordeum vulgare convar. labile) have a six-rowed genetic background

Labile-barleys (Hordeum vulgare convar. labile (Schiem.) Mansf.) are found in the highlands of Ethiopia, Eretria and North India-Pakistan districts. They represent a distinct spike form showing row-type alterations even within individual spikes of the same genotypes. Variation at the six-rowed spike 1 (vrs1) locus is sufficient to control barley lateral spikelet fertility, which is also modified by alleles at the intermedium-c (int-c) locus. This study aimed at re-sequencing these two loci to investigate whether labile-barleys have a two-rowed genetic background, resulting in increased lateral spikelet fertility, or show reduced lateral spikelet fertility if they possess a six-rowed genetic background. The Vrs1 re-sequencing results of 221 supposedly labile-barley accessions from Ethiopia revealed 13 accessions with two novel vrs1.a1 haplotypes. Following the current nomenclature of vrs1 haplotypes, the new haplotypes were named as haplotypes 66 and 67. Re-sequencing at the int-c locus showed that 118 of the labile-barleys possessed the previously described Int-c.a allele but only one accession was found having a novel Int-c.a haplotype in the homozygous state (termed Int-c.a haplotype1; Hap_1). Interestingly, 101 labile-barleys carried the Int-c.a allele and Int-c.a haplotype1 simultaneously, suggesting maintained heterozygosity or recent gene duplication at this locus. Only one accession had a two-rowed haplotype (Vrs1.b3, int-c.b1) and one accession possessed the Vrs1.t (deficiens) and Int-c.a alleles (six-rowed). These two accessions were considered as misclassified labile genotypes and not included in further analysis. Thus, these results confirmed that all of the 219 labile accessions studied in this work showed six-rowed alleles at vrs1 but reduced lateral spikelet fertility. This reduction is most likely caused by the recessive labile (lab) locus which we are in the process to characterize further.     

Single-copy nuclear PolA1 gene sheds light on the origin of S genome with relationships to B and G genomes of polyploid wheat species

PolA1, a single-copy nuclear gene encoding the largest subunit of RNA polymerase I, comprises highly polymorphic intron 19 and nucleotide tag (Ntag) sequences. We analyzed these sequences in 42 accessions, which differed in ploidy, of Triticum–Aegilops and Hordeum species. The lengths of the intron 19 sequences were ca. 110?bp long in Triticum–Aegilops species, except in four Sitopsis species, Ae. longissima, Ae. searsii, Ae. sharonensis, Ae. speltoides, which had introns similar in length to those of Hordeum species, i.e., ca. 240?bp long. Phylogenetic analyses of the Ntag sequences showed that the four Sitopsis and remaining Triticum–Aegilops species were classified into two discrete Hordeum and Triticum clades, respectively. The A and D genome-specific Ntag sequences of polyploid wheats were highly homologous with those of T. urartu and Ae. tauschii, respectively. In Ae. bicornis, another Sitopsis species, two accessions had the short intron 19 and Triticum–type Ntag sequence, which were highly homologous with those of the B genome in polyploid wheats, whereas one accession contained the long intron 19 and Hordeum–type Ntag sequences. In contrast, partial sequence analyses revealed that the three accessions of Ae. bicornis shared highly homology to single-copy DMC1 and EF-G genes. The discrepancy between these results indicates that the Sitopsis species were probably established by hybrid speciation including ancient introgressive hybridization between progenitors of Triticum–Aegilops and Hordeum. Although many researchers have proposed Ae. speltoides as a candidate for the B genome donor, our data suggest the existence of diploid B genome species in the past that were responsible for the origin of both polyploid wheats and Sitopsis species, including Ae. speltoides.     

Biodiversity of Diploid D-Genome Aegilops Tauschii Coss. in Iran Measured Using Microsatellites

Microsatellite markers were used to analyse the biodiversity of 57 accessions of different subspecies and varieties of wild Aegilops tauschii (2n = 2x = 14; D genome) collected across the major areas where it grows in Iran. Levels of diversity were high, with numbers of alleles averaging 7.3 (ranging up to 12) and polymorphism information contents averaging 0.6591. One accession was notably more similar to two of the D genome in hexaploid wheats (Triticum aestivum) used as outgroups. Within the Ae. tauschii accessions, no markers were characteristic for taxa or geographical origin, suggesting high gene flow between the subspecies and varieties, although some groupings, which could be related to geographical origin, were evident. This survey demonstrates the high diversity present in wild goatgrass in Iran, and indicates that there is value in sampling for useful genes for wheat breeding.