Population structure of wild wheat D-genome progenitor Aegilops tauschii Coss.: implications for intraspecific lineage diversification and evolution of common wheat

On the evolution and genetic diversity of the bread wheat D genome

Bread wheat (Triticum aestivum) became a globally dominant crop after incorporating the D genome from the donor species Aegilops tauschii, but the evolutionary history that shaped the D genome during this process remains to be clarified. Here, we propose a renewed evolutionary model linking Ae. tauschii and the hexaploid wheat D genome by constructing an ancestral haplotype map covering 762 Ae. tauschii and hexaploid wheat accessions. We dissected the evolutionary trajectories of Ae. tauschii lineages and reported a few independent intermediate accessions, demonstrating that low-frequency inter-sublineage gene flow had enriched the diversity of Ae. tauschii. We discovered that the D genome of hexaploid wheat was inherited from a unified ancestral template, but with a mosaic composition that was highly mixed and derived mainly from three Ae. tauschii L2 sublineages located in the Caspian coastal region. This result suggests that early agricultural activities facilitated innovations in D-genome composition and finalized the success of hexaploidization. We found that the majority (51.4%) of genetic diversity was attributed to novel mutations absent in Ae. tauschii, and we identified large Ae. tauschii introgressions from various lineages, which expanded the diversity of the wheat D genome and introduced beneficial alleles. This work sheds light on the process of wheat hexaploidization and highlights the evolutionary significance of the multi-layered genetic diversity of the bread wheat D genome.

Origin and evolution of the bread wheat D genome

Bread wheat (Triticum aestivum) is a globally dominant crop and major source of calories and proteins for the human diet. Compared with its wild ancestors, modern bread wheat shows lower genetic diversity, caused by polyploidisation, domestication and breeding bottlenecks1,2. Wild wheat relatives represent genetic reservoirs, and harbour diversity and beneficial alleles that have not been incorporated into bread wheat. Here we establish and analyse extensive genome resources for Tausch's goatgrass (Aegilops tauschii), the donor of the bread wheat D genome. Our analysis of 46 Ae. tauschii genomes enabled us to clone a disease resistance gene and perform haplotype analysis across a complex disease resistance locus, allowing us to discern alleles from paralogous gene copies. We also reveal the complex genetic composition and history of the bread wheat D genome, which involves contributions from genetically and geographically discrete Ae. tauschii subpopulations. Together, our results reveal the complex history of the bread wheat D genome and demonstrate the potential of wild relatives in crop improvement.

Rmg10, a Novel Wheat Blast Resistance Gene Derived from Aegilops tauschii

Wheat blast, caused by Pyricularia oryzae (syn. Magnaporthe oryzae) pathotype Triticum (MoT), is a devastating disease that can result in up to 100% yield loss in affected fields. To find new resistance genes against wheat blast, we screened 199 accessions of Aegilops tauschii, the D genome progenitor of common wheat (Triticum aestivum), by seedling inoculation assays with Brazilian MoT isolate Br48 and found 14 resistant accessions. A synthetic hexaploid wheat line (Ldn/KU-2097) derived from a cross between the T. turgidum 'Langdon' (Ldn) and resistant A. tauschii accession KU-2097 exhibited resistance in seedlings and spikes against Br48. In an F2 population derived from 'Chinese Spring' × Ldn/KU-2097, resistant and susceptible individuals segregated in a 3:1 ratio, suggesting that the resistance from KU-2097 is controlled by a single dominant gene. We designated this gene Rmg10. Genetic mapping using an F2:3 population from the same cross mapped the RMG10 locus to the short arm of chromosome 2D. Rmg10 was ineffective against Bangladesh isolates but effective against Brazilian isolates. Field tests in Bolivia showed increased spike resistance in a synthetic octaploid wheat line produced from a cross between common wheat cultivar 'Gladius' and KU-2097. These results suggest that Rmg10 would be beneficial in farmers' fields in South America.

Metabolome and Mycobiome of Aegilops tauschii Subspecies Differing in Susceptibility to Brown Rust and Powdery Mildew Are Diverse

The present study demonstrated the differences in the seed metabolome and mycobiome of two Aegilops tauschii Coss accessions with different resistance to brown rust and powdery mildew. We hypothesized that the seeds of resistant accession k-1958 Ae. tauschii ssp. strangulata can contain a larger number of metabolites with antifungal activity compared with the seeds of susceptible Ae. tauschii ssp meyeri k-340, which will determine differences in the seed fungal community. Our study emphasizes the differences in the seed metabolome of the studied Ae. tauschii accessions. The resistant accession k-1958 had a higher content of glucose and organic acids, including pyruvic, salicylic and azelaic acid, as well as pipecolic acids, galactinol, glycerol and sitosterol. The seeds of Ae. tauschii-resistant accession k-1958 were found to contain more active substances with antifungal activity. The genera Cladosporium and Alternaria were dominant in the seed mycobiome of the resistant accession. The genera Alternaria, Blumeria and Cladosporium dominated in seed mycobiome of susceptible accession k-340. In the seed mycobiome of the resistant k-1958, a higher occurrence of saprotrophic micromycetes was found, and many of the micromycetes were biocontrol agents. It was concluded that differences in the seed metabolome of Ae. tauschii contributed to the determination of the differences in mycobiomes.

Discrepancy of flowering time between genetically close sublineages of Aegilops umbellulata Zhuk

Aegilops umbellulata Zhuk., a wild diploid wheat-related species, has been used as a genetic resource for several important agronomic traits. However, its genetic variations have not been comprehensively studied. We sequenced RNA from 114 accessions of Ae. umbellulata to evaluate DNA polymorphisms and phenotypic variations. Bayesian clustering and phylogenetic analysis based on SNPs detected by RNA sequencing revealed two divergent lineages, UmbL1 and UmbL2. The main differences between them were in the sizes of spikes and spikelets, and culm diameter. UmbL1 is divided into two sublineages, UmbL1e and UmbL1w. These genetic differences corresponded to geographic distributions. UmbL1e, UmbL1w, and UmbL2 are found in Turkey, Iran/Iraq, and Greece, respectively. Although UmbL1e and UmbL1w were genetically similar, flowering time and other morphological traits were more distinct between these sublineages than those between the lineages. This discrepancy can be explained by the latitudinal and longitudinal differences in habitats. Specifically, latitudinal clines of flowering time were clearly observed in Ae. umbellulata, strongly correlated with solar radiation in the winter season. This observation implies that latitudinal differences are a factor in differences in the flowering times of Ae. umbellulata. Differences in flowering time could influence other morphological differences and promote genetic divergence between sublineages.

Genetic diversity and population structure of modern wheat (Triticum aestivum L.) cultivars in Henan Province of China based on SNP markers

BACKGROUND

Henan is the province with the greatest wheat production in China. Although more than 100 cultivars are used for production, many cultivars are still insufficient in quality, disease resistance, adaptability and yield potential. To overcome these limitations, it is necessary to constantly breed new cultivars to maintain the continuous and stable growth of wheat yield and quality. To improve breeding efficiency, it is important to evaluate the genetic diversity and population genetic structure of its cultivars. However, there are no such reports from Henan Province. Therefore, in this study, single nucleotide polymorphism (SNP) markers were used to study the population genetic structure and genetic diversity of 243 wheat cultivars included in a comparative test of wheat varieties in Henan Province, aiming to provide a reference for the utilization of backbone parents and the selection of hybrid combinations in the genetic improvement of wheat cultivars.

RESULTS

In this study, 243 wheat cultivars from Henan Province of China were genotyped by the Affymetrix Axiom Wheat660K SNP chip, and 21 characteristics were investigated. The cultivars were divided into ten subgroups; each subgroup had distinct characteristics and unique utilization value. Furthermore, based on principal component analysis, Zhoumai cultivars were the main hybrid parents, followed by Aikang 58, high-quality cultivars, and Shandong cultivars. Genetic diversity analysis showed that 61.3% of SNPs had a high degree of genetic differentiation, whereas 33.4% showed a moderate degree. The nucleotide diversity of subgenome B was relatively high, with an average π value of 3.91E-5; the nucleotide diversity of subgenome D was the lowest, with an average π value of 2.44E-5.

CONCLUSION

The parents used in wheat cross-breeding in Henan Province are similar, with a relatively homogeneous genetic background and low genetic diversity. These results will not only contribute to the objective evaluation and utilization of the tested cultivars but also provide insights into the current conditions and existing challenges of wheat cultivar breeding in Henan Province, thereby facilitating the scientific formulation of breeding objectives and strategies to improve breeding efficiency.

Genome-wide analysis of a collective grave from Mentesh Tepe provides insight into the population structure of early neolithic population in the South Caucasus

Despite the localisation of the southern Caucasus at the outskirt of the Fertile Crescent, the Neolithisation process started there only at the beginning of the sixth millennium with the Shomutepe-Shulaveri culture of yet unclear origins. We present here genomic data for three new individuals from Mentesh Tepe in Azerbaijan, dating back to the beginnings of the Shomutepe-Shulaveri culture. We evidence that two juveniles, buried embracing each other, were brothers. We show that the Mentesh Tepe Neolithic population is the product of a recent gene flow between the Anatolian farmer-related population and the Caucasus/Iranian population, demonstrating that population admixture was at the core of the development of agriculture in the South Caucasus. By comparing Bronze Age individuals from the South Caucasus with Neolithic individuals from the same region, including Mentesh Tepe, we evidence that gene flows between Pontic Steppe populations and Mentesh Tepe-related groups contributed to the makeup of the Late Bronze Age and modern Caucasian populations. Our results show that the high cultural diversity during the Neolithic period of the South Caucasus deserves close genetic analysis.

Seed position in spikelet as a contributing factor to the success of the winter annual invasive grass Aegilops tauschii

Seed position - dependent effects on seed dormancy/germination are well documented at the inflorescence/infructescence level, but less is known about seeds at different positions within a dispersal unit. For the invasive winter annual grass Aegilops tauschii, we quantified morphology, mass and dormancy/germination of seeds from basal (1), middle (2), and distal (3) positions in two spikelet types (Left and Right). We also investigated seedling emergence, survival, plant size and seed production of plants from seeds in different spikelet positions of two spikelet types under different soil nutrient and water conditions. We found that these seed, seedling and plant traits performed as mirror images between the Left and Right spikelet types. The middle seed was significantly the longest and had the maximum mass, while the basal seed was the shortest and had medium mass. Middle seeds had the highest increase in mass during imbibition and the highest germination percentages and rates, while basal seeds had the lowest. Seedling emergence and survival, plant size and seed production for each position of seeds were highest in the added fertilizer combined with regular watering treatment and lowest in the no fertilizer combined with natural moisture, while height of plants derived from the middle and the distal seeds was significantly higher than that of plants derived from the basal seeds under all soil nutrient and water conditions. Seedling survival, number of tillers per plant and seed production per plant from the middle and distal seeds were significantly lower than those from basal seeds under all soil nutrient and water treatments. The considerable variation in seedling emergence and survival, plant size and seed production between seeds in different positions in the spikelet results in much flexibility in all stages of the life cycle, thereby likely contributing to the invasiveness of A. tauschii.

High genetic diversity in Aegilops tauschii Coss. accessions from North Iran as revealed by IRAP and REMAP markers

BACKGROUND

Aegilops tauschii Coss. as a donor of wheat D genome has an important role in wheat breeding programs. Genetic and phylogeographic diversity of 79 Ae. tauschii accessions collected from north and northwest of Iran were analyzed based on retroelement insertional polymorphisms using inter-retrotransposon amplified polymorphism (IRAP) and retrotransposon-microsatellite amplified polymorphism (REMAP) markers.

RESULTS

In total, 306 and 151 polymorphic bands were amplified in IRAP and REMAP analyses, respectively. As a result, a high level of polymorphism was observed among the studied accessions as revealed by an average of 25.5 bands per primer/primer combination and mean PIC value of 0.47 in IRAP and an average of 25.16 bands per primer combination and mean PIC value of 0.47 in REMAP. Genetic relationships of the accessions were analyzed using distance- and model-based cluster analyses.

CONCLUSION

The result showed that genetic distance did not seem to be related to geographic distribution, and the accessions could be divided into three groups, which was further supported by principal coordinate analysis. These results on genetic diversity and population structure of Ae. tauschii in Iran should provide important knowledge on genetic resources and their applications in wheat breeding programs.

Population genomic analysis of Aegilops tauschii identifies targets for bread wheat improvement

Aegilops tauschii, the diploid wild progenitor of the D subgenome of bread wheat, is a reservoir of genetic diversity for improving bread wheat performance and environmental resilience. Here we sequenced 242 Ae. tauschii accessions and compared them to the wheat D subgenome to characterize genomic diversity. We found that a rare lineage of Ae. tauschii geographically restricted to present-day Georgia contributed to the wheat D subgenome in the independent hybridizations that gave rise to modern bread wheat. Through k-mer-based association mapping, we identified discrete genomic regions with candidate genes for disease and pest resistance and demonstrated their functional transfer into wheat by transgenesis and wide crossing, including the generation of a library of hexaploids incorporating diverse Ae. tauschii genomes. Exploiting the genomic diversity of the Ae. tauschii ancestral diploid genome permits rapid trait discovery and functional genetic validation in a hexaploid background amenable to breeding.

Assessing the Potential of Using the Langdon 5D(5B) Substitution Line for the Introgression of Aegilops tauschii Into Durum Wheat

Aegilops tauschii, the D-genome donor of hexaploid wheat, provides a source of genetic variation that could be used for tetraploid (durum) wheat improvement. In addition to the genes for wheat quality on the D-genome, which differentiate between bread and durum wheats in terms of end-use properties, genes coding for resistances to biotic and abiotic stresses are also present on the D-genome which would be useful in durum wheat. The introgression of Ae. tauschii into durum wheat, however, requires cytogenetic manipulation to induce homoeologous chromosome pairing to promote recombination. For this purpose, the introgression of Ae. tauschii into durum wheat was performed through a bridge cross of the wild species to the Langdon 5D(5B) disomic substitution line that lacks the Ph1 locus present on chromosome 5B, followed by a cross of the F₁ to the durum wheat cultivar Om Rabi 5. Subsequent generations were self-fertilized, and these were screened for D-genome introgressions using (i) D-genome-specific Kompetitive Allele-Specific PCR (KASP) markers and (ii) KASP markers polymorphic between the 5D chromosomes of wheat, present in the Langdon 5D(5B) substitution line, and of Ae. tauschii. Homozygous introgression lines were confirmed using genomic and fluorescence in situ hybridization. The results showed that the use of the Langdon 5D(5B) disomic substitution line did not promote D-genome introgression across all linkage groups with only a limited success in the introgression of Ae. tauschii 5D segments into durum wheat.

New insights into the dispersion history and adaptive evolution of taxon Aegilops tauschii in China

Aegilops tauschii, the wild progenitor of wheat D-genome and a valuable germplasm for wheat improvement, has a wide natural distribution from eastern Turkey to China. However, the phylogenetic relationship and dispersion history of Ae. tauschii in China has not been scientifically clarified. In this study, we genotyped 208 accessions (with 104 in China) using ddRAD sequencing and 55K SNP array, and classified the population into six sublineages. Three possible spreading routes or events were identified, resulting in specific distribution patterns, with four sublineages found in Xinjiang, one in Qinghai, two in Shaanxi and one in Henan. We also established the correlation of SNP-based, karyotype-based and spike-morphology-based techniques to demonstrate the internal classification of Ae. tauschii, and developed consensus dataset with 1245 putative accessions by merging data previously published. Our analysis suggested that eight inter-lineage accessions could be assigned to the putative Lineage 3 and these accessions would help to conserve the genetic diversity of the species. By developing the consensus phylogenetic relationships of Ae. tauschii, our work validated the hypothesis on the dispersal history of Ae. tauschii in China, and contributed to the efficient and comprehensive germplasm-mining of the species.

Differentially Amplified Repetitive Sequences Among Aegilops tauschii Subspecies and Genotypes

Genomic repetitive sequences commonly show species-specific sequence type, abundance, and distribution patterns, however, their intraspecific characteristics have been poorly described. We quantified the genomic repetitive sequences and performed single nucleotide polymorphism (SNP) analysis between 29 Ae. tauschii genotypes and subspecies using publicly available raw genomic Illumina sequence reads and used fluorescence in situ hybridization (FISH) to experimentally analyze some repeats. The majority of the identified repetitive sequences had similar contents and proportions between anathera, meyeri, and strangulata subspecies. However, two Ty3/gypsy retrotransposons (CL62 and CL87) showed significantly higher abundances, and CL1, CL119, CL213, CL217 tandem repeats, and CL142 retrotransposon (Ty1/copia type) showed significantly lower abundances in subspecies strangulata compared with the subspecies anathera and meyeri. One tandem repeat and 45S ribosomal DNA (45S rDNA) abundances showed a high variation between genotypes but their abundances were not subspecies specific. Phylogenetic analysis using the repeat abundances of the aforementioned clusters placed the strangulata subsp. in a distinct clade but could not discriminate anathera and meyeri. A near complete differentiation of anathera and strangulata subspecies was observed using SNP analysis; however, var. meyeri showed higher genetic diversity. FISH using major tandem repeats couldn't detect differences between subspecies, although (GAA)₁₀ signal patterns generated two different karyotype groups. Taken together, the different classes of repetitive DNA sequences have differentially accumulated between strangulata and the other two subspecies of Ae. tauschii that is generally in agreement with spike morphology, implying that factors affecting repeatome evolution are variable even among highly closely related lineages.

RNA-Seq-based DNA marker analysis of the genetics and molecular evolution of Triticeae species

The release of high-quality chromosome-level genome sequences of members of the Triticeae tribe has greatly facilitated genetic and genomic analyses of important crops such as wheat (Triticum aestivum) and barley (Hordeum vulgare). Due to the large diploid genome size of Triticeae plants (ca. 5 Gbp), transcript analysis is an important method for identifying genetic and genomic differences among Triticeae species. In this review, we summarize our results of RNA-Seq analyses of diploid wheat accessions belonging to the genera Aegilops and Triticum. We also describe studies of the molecular relationships among these accessions and provide insight into the evolution of common hexaploid wheat. DNA markers based on polymorphisms within species can be used to map loci of interest. Even though the genome sequence of diploid Aegilops tauschii, the D-genome donor of common wheat, has been released, the diploid barley genome continues to provide key information about the physical structures of diploid wheat genomes. We describe how a series of RNA-Seq analyses of wheat relatives has helped uncover the structural and evolutionary features of genomic and genetic systems in wild and cultivated Triticeae species.

Homeolog expression quantification methods for allopolyploids

Genome duplication with hybridization, or allopolyploidization, occurs in animals, fungi and plants, and is especially common in crop plants. There is an increasing interest in the study of allopolyploids because of advances in polyploid genome assembly; however, the high level of sequence similarity in duplicated gene copies (homeologs) poses many challenges. Here we compared standard RNA-seq expression quantification approaches used currently for diploid species against subgenome-classification approaches which maps reads to each subgenome separately. We examined mapping error using our previous and new RNA-seq data in which a subgenome is experimentally added (synthetic allotetraploid Arabidopsis kamchatica) or reduced (allohexaploid wheat Triticum aestivum versus extracted allotetraploid) as ground truth. The error rates in the two species were very similar. The standard approaches showed higher error rates (>10% using pseudo-alignment with Kallisto) while subgenome-classification approaches showed much lower error rates (<1% using EAGLE-RC, <2% using HomeoRoq). Although downstream analysis may partly mitigate mapping errors, the difference in methods was substantial in hexaploid wheat, where Kallisto appeared to have systematic differences relative to other methods. Only approximately half of the differentially expressed homeologs detected using Kallisto overlapped with those by any other method in wheat. In general, disagreement in low-expression genes was responsible for most of the discordance between methods, which is consistent with known biases in Kallisto. We also observed that there exist uncertainties in genome sequences and annotation which can affect each method differently. Overall, subgenome-classification approaches tend to perform better than standard approaches with EAGLE-RC having the highest precision.

Exploring the alpha-gliadin locus: the 33-mer peptide with six overlapping coeliac disease epitopes in Triticum aestivum is derived from a subgroup of Aegilops tauschii

Most alpha-gliadin genes of the Gli-D2 locus on the D genome of hexaploid bread wheat (Triticum aestivum) encode for proteins with epitopes that can trigger coeliac disease (CD), and several contain a 33-mer peptide with six partly overlapping copies of three epitopes, which is regarded as a remarkably potent T-cell stimulator. To increase genetic diversity in the D genome, synthetic hexaploid wheat lines are being made by hybridising accessions of Triticum turgidum (AB genome) and Aegilops tauschii (the progenitor of the D genome). The diversity of alpha-gliadins in A. tauschii has not been studied extensively. We analysed the alpha-gliadin transcriptome of 51 A. tauschii accessions representative of the diversity in A. tauschii. We extracted RNA from developing seeds and performed 454 amplicon sequencing of the first part of the alpha-gliadin genes. The expression profile of allelic variants of the alpha-gliadins was different between accessions, and also between accessions of the Western and Eastern clades of A. tauschii. Generally, both clades expressed many allelic variants not found in bread wheat. In contrast to earlier studies, we detected the 33-mer peptide in some A. tauschii accessions, indicating that it was introduced along with the D genome into bread wheat. In these accessions, transcripts with the 33-mer peptide were present at lower frequencies than in bread wheat varieties. In most A. tauschii accessions, however, the alpha-gliadins do not contain the epitope, and this may be exploited, through synthetic hexaploid wheats, to breed bread wheat varieties with fewer or no coeliac disease epitopes.

Genome-Wide Association Study of Morpho-Physiological Traits in Aegilops tauschii to Broaden Wheat Genetic Diversity

Aegilops tauschii, the D-genome donor of bread wheat, is a storehouse of genetic diversity that can be used for wheat improvement. This species consists of two main lineages (TauL1 and TauL2) and one minor lineage (TauL3). Its morpho-physiological diversity is large, with adaptations to a wide ecological range. Identification of allelic diversity in Ae. tauschii is of utmost importance for efficient breeding and widening of the genetic base of wheat. This study aimed at identifying markers or genes associated with morpho-physiological traits in Ae. tauschii, and at understanding the difference in genetic diversity between the two main lineages. We performed genome-wide association studies of 11 morpho-physiological traits of 343 Ae. tauschii accessions representing the entire range of habitats using 34,829 DArTseq markers. We observed a wide range of morpho-physiological variation among all accessions. We identified 23 marker-trait associations (MTAs) in all accessions, 15 specific to TauL1 and eight specific to TauL2, suggesting independent evolution in each lineage. Some of the MTAs could be novel and have not been reported in bread wheat. The markers or genes identified in this study will help reveal the genes controlling the morpho-physiological traits in Ae. tauschii, and thus in bread wheat even if the plant morphology is different.

GRAS-Di system facilitates high-density genetic map construction and QTL identification in recombinant inbred lines of the wheat progenitor Aegilops tauschii

Due to large and complex genomes of Triticeae species, skim sequencing approaches have cost and analytical advantages for detecting genetic markers and building linkage maps. Here, we develop a high-density linkage map and identify quantitative trait loci (QTLs) for recombinant inbred lines of Aegilops tauschii, a D-genome donor of bread wheat, using the recently developed genotyping by Random Amplicon Sequencing-Direct (GRAS-Di) system, which facilitates skimming of the large and complicated genome and generates a large number of genetic markers. The deduced linkage groups based on the GRAS-Di genetic markers corresponded to the chromosome number of Ae. tauschii. We successfully identified stable QTLs for flowering time and spikelet shape-related traits. Genotype differences of RILs at the QTL-linked markers were significantly associated with the trait variations. In particular, one of the QTL-linked markers for flowering time was mapped close to VRN3 (also known as FLOWERING LOCUS T), which controls flowering. The GRAS-Di system is, therefore, an efficient and useful application for genotyping and linkage mapping in species with large and complex genomes, such as Triticeae species.

The Allelic Diversity of the Gibberellin Signaling Pathway Genes in Aegilops tauschii Coss

Gibberellin-insensitive reduced height genes are widely spread in modern wheat varieties, making them resistant to lodging under conditions of intensive farming. However, the limited diversity of these genes present in wheat germplasm can limit the adaptability of newly created cultivars to the changing climate. The diversity of the gibberellin signaling pathway genes involved in plant height control- Reduced height 1 (Rht-D1), Gibberellin-insensitive dwarf 1 (Gid1‑D) and Gibberellin-insensitive dwarf 2 (Gid2-D)-was studied in the diploid wild goatgrass Aegilops tauschii Coss., one of the ancestral species of the bread wheat (Triticum aestivum L.) and the donor of its D subgenome, using high-throughput sequencing. The examination of 24 Ae. tauschii accessions of different geographical origins revealed a large number of new alleles (haplotypes) not found in bread wheat varieties. Some of the detected polymorphisms lead to changes in the amino acid sequence of proteins. Four isoforms (amino acid sequence variants) were found for the RHT-D1 protein, and two isoforms-for the GID1 and GID2 proteins, each. An analysis of the co-occurrence frequencies of various isoforms of the three proteins showed that their combinations were not random in Ae. tauschii, which may indicate the functional significance of their differences. New alleles of the Rht-D1, Gid1-D, and Gid2-D genes are promising for introgression into bread wheat and studying their effect on plant height and adaptability.

Reproductive and genetic roles of the maternal progenitor in the origin of common wheat (Triticum aestivum L.)

Common wheat (Triticum aestivum L., AABBDD genome) is thought to have emerged through natural hybridization between Triticum turgidum L. (AABB genome) and Aegilops tauschii Coss. (DD genome). Hybridization barriers and doubling of the trihaploid F1 hybrids' genome (ABD) via unreduced gamete fusion had key roles in the process. However, how T. turgidum, the maternal progenitor, was involved in these mechanisms remains unknown. An artificial cross-experiment using 46 cultivated and 31 wild T. turgidum accessions and a single Ae. tauschii tester with a very short genetic distance to the common wheat D genome was conducted. Cytological and quantitative trait locus analyses of F1 hybrid genome doubling were performed. The crossability and ability to cause hybrid inviability did not greatly differ between the cultivars and wild accessions. The ability to cause hybrid genome doubling was higher in the cultivars. Three novel T. turgidum loci for hybrid genome doubling, which influenced unreduced gamete production in F1 hybrids, were identified. Cultivated T. turgidum might have increased the probability of the emergence of common wheat through its enhanced ability to cause genome doubling in F1 hybrids with Ae. tauschii. The ability enhancement might have involved alterations at a relatively small number of loci.

Mapping QTLs for enhancing early biomass derived from Aegilops tauschii in synthetic hexaploid wheat

Strong early vigour plays a crucial role in wheat yield improvement by enhancing resource utilization efficiency. Synthetic hexaploid wheat (SHW) combines the elite genes of tetraploid wheat with Aegilops tauschii and has been widely used in wheat genetic improvement for its abundant genetic variation. The two SHWs Syn79 and Syn80 were derived from the crossing of the same tetraploid wheat DOY1 with two different Ae. tauschii accessions, AT333 and AT428, respectively. The Syn80 possessed better early vigour traits than Syn79, theretically caused by their D genome from Ae. tauschii. To dissect their genetic basis in a hexaploid background, 203 recombinant inbred lines (RILs) derived from the cross of Syn79 x Syn80 were developed to detect quantitative trait loci (QTL) for four early biomass related traits: plant height (PH), tiller number (TN), shoot fresh weight (SFW) and shoot dry weight (SDW) per plant, under five different environmental conditions. Determined from the data of SNP markers, two genome regions on 1DS and 7D were stably associated with the four early biomass related traits showing pleiotropic effects. Four stable QTLs QPh.saas-1DS, QTn.saas-1DS, QSfw.saas-1DS and QSdw.saas-1DS explaining 7.92, 15.34, 9.64 and 10.15% of the phenotypic variation, respectively, were clustered in the region of 1DS from AX-94812958 to AX-110910133. Meanwhile, QPh.saas-7D, QTn.saas-7D, QSfw.saas-7D and QSdw.saas-7D were flanked by AX-109917900 and AX-110605376 on 7D, explaining 16.12, 24.35, 15.25 and 13.37% of the phenotypic variation on average, respectively. Moreover, these genomic QTLs on 1DS and 7D enhancing biomass in the parent Syn80 were from Ae. tauschii AT428. These findings suggest that these two QTLs from Ae. tauschii can be expressed stably in a hexaploid background at the jointing stage and be used for wheat improvement.

Comparative Transcriptome Analysis of Two Aegilops tauschii with Contrasting Drought Tolerance by RNA-Seq

As the diploid progenitor of common wheat, Aegilops tauschii is considered to be a valuable resistance source to various biotic and abiotic stresses. However, little has been reported concerning the molecular mechanism of drought tolerance in Ae. tauschii. In this work, the drought tolerance of 155 Ae. tauschii accessions was firstly screened on the basis of their coleoptile lengths under simulated drought stress. Subsequently, two accessions (XJ002 and XJ098) with contrasting coleoptile lengths were selected and intensively analyzed on rate of water loss (RWL) as well as physiological characters, confirming the difference in drought tolerance at the seedling stage. Further, RNA-seq was utilized for global transcriptome profiling of the two accessions seedling leaves under drought stress conditions. A total of 6969 differentially expressed genes (DEGs) associated with drought tolerance were identified, and their functional annotations demonstrated that the stress response was mediated by pathways involving alpha-linolenic acid metabolism, starch and sucrose metabolism, peroxisome, mitogen-activated protein kinase (MAPK) signaling, carbon fixation in photosynthetic organisms, and glycerophospholipid metabolism. In addition, DEGs with obvious differences between the two accessions were intensively analyzed, indicating that the expression level of DEGs was basically in alignment with the physiological changes of Ae. tauschii under drought stress. The results not only shed fundamental light on the regulatory process of drought tolerance in Ae. tauschii, but also provide a new gene resource for improving the drought tolerance of common wheat.

Allele mining of TaGRF-2D gene 5'-UTR in Triticum aestivum and Aegilops tauschii genotypes

The low diversity of the D-subgenome of bread wheat requires the involvement of new alleles for breeding. In grasses, the allelic state of Growth Regulating Factor (GRF) gene is correlated with nitrogen uptake. In this study, we characterized the sequence of TaGRF-2D and assessed its diversity in bread wheat and goatgrass Aegilops tauschii (genome DD). In silico analysis was performed for reference sequence searching, primer pairs design and sequence assembly. The gene sequence was obtained using Illumina and Sanger sequencing. The complete sequences of TaGRF-2D were obtained for 18 varieties of wheat. The polymorphism in the presence/absence of two GCAGCC repeats in 5' UTR was revealed and the GRF-2D-SSR marker was developed. Our results showed that the alleles 5' UTR-250 and 5' UTR-238 were present in wheat varieties, 5' UTR-250 was presented in the majority of wheat varieties. In Ae. tauschii ssp. strangulata (likely donor of the D-subgenome of polyploid wheat), most accessions carried the 5' UTR-250 allele, whilst most Ae. tauschii ssp. tauschii have 5' UTR-244. The developed GRF-2D-SSR marker can be used to study the genetic diversity of wheat and Ae. tauschii.

Phenotypic effects of the U-genome variation in nascent synthetic hexaploids derived from interspecific crosses between durum wheat and its diploid relative Aegilops umbellulata

Aegilops umbellulata is a wild diploid wheat species with the UU genome that is an important genetic resource for wheat breeding. To exploit new synthetic allohexaploid lines available as bridges for wheat breeding, a total of 26 synthetic hexaploid lines were generated through crossing between the durum wheat cultivar Langdon and 26 accessions of Ae. umbellulata. In nascent synthetic hexaploids with the AABBUU genome, the presence of the set of seven U-genome chromosomes was confirmed with U-genome chromosome-specific markers developed based on RNA-seq-derived data from Ae. umbellulata. The AABBUU synthetic hexaploids showed large variations in flowering- and morphology-related traits, and these large variations transmitted well from the parental Ae. umbellulata accessions. However, the variation ranges in most traits examined were reduced under the AABBUU hexaploid background compared with under the diploid parents. The AABBUU and AABBDD synthetic hexaploids were clearly discriminated by several morphological traits, and an increase of plant height and in the number of spikes and a decrease of spike length were commonly observed in the AABBUU synthetics. Thus, interspecific differences in several morphological traits between Ae. umbellulata and A. tauschii largely affected the basic plant architecture of the synthetic hexaploids. In conclusion, the AABBUU synthetic hexaploid lines produced in the present study are useful resources for the introgression of desirable genes from Ae. umbellulata to common wheat.

Diploid genome differentiation conferred by RNA sequencing-based survey of genome-wide polymorphisms throughout homoeologous loci in Triticum and Aegilops

Background

Triticum and Aegilops diploid species have morphological and genetic diversity and are crucial genetic resources for wheat breeding. According to the chromosomal pairing-affinity of these species, their genome nomenclatures have been defined. However, evaluations of genome differentiation based on genome-wide nucleotide variations are still limited, especially in the three genomes of the genus Aegilops: Ae. caudata L. (CC genome), Ae. comosa Sibth. et Sm. (MM genome), and Ae. uniaristata Vis. (NN genome). To reveal the genome differentiation of these diploid species, we first performed RNA-seq-based polymorphic analyses for C, M, and N genomes, and then expanded the analysis to include the 12 diploid species of Triticum and Aegilops.

Results

Genetic divergence of the exon regions throughout the entire chromosomes in the M and N genomes was larger than that between A- and A^m-genomes. Ae. caudata had the second highest genetic diversity following Ae. speltoides, the putative B genome donor of common wheat. In the phylogenetic trees derived from the nuclear and chloroplast genome-wide polymorphism data, the C, D, M, N, U, and S genome species were connected with short internal branches, suggesting that these diploid species emerged during a relatively short period in the evolutionary process. The highly consistent nuclear and chloroplast phylogenetic topologies indicated that nuclear and chloroplast genomes of the diploid Triticum and Aegilops species coevolved after their diversification into each genome, accounting for most of the genome differentiation among the diploid species.

Conclusions

RNA-sequencing-based analyses successfully evaluated genome differentiation among the diploid Triticum and Aegilops species and supported the chromosome-pairing-based genome nomenclature system, except for the position of Ae. speltoides. Phylogenomic and epigenetic analyses of intergenic and centromeric regions could be essential for clarifying the mechanisms behind this inconsistency.

The complete chloroplast genomes of seventeen Aegilops tauschii: genome comparative analysis and phylogenetic inference

The D genome progenitor of bread wheat, Aegilops tauschii Cosson (DD, 2n = 2x = 14), which is naturally distributed in Central Eurasia, ranging from northern Syria and Turkey to western China, is considered a potential genetic resource for improving bread wheat. In this study, the chloroplast (cp) genomes of 17 Ae. tauschii accessions were reconstructed. The cp genome sizes ranged from 135,551 bp to 136,009 bp and contained a typical quadripartite structure of angiosperms. Within these genomes, we identified a total of 124 functional genes, including 82 protein-coding genes, 34 transfer RNA genes and eight ribosomal RNA genes, with 17 duplicated genes in the IRs. Although the comparative analysis revealed that the genomic structure (gene order, gene number and IR/SC boundary regions) is conserved, a few variant loci were detected, predominantly in the non-coding regions (intergenic spacer regions). The phylogenetic relationships determined based on the complete genome sequences were consistent with the hypothesis that Ae. tauschii populations in the Yellow River region of China originated in South Asia not Xinjiang province or Iran, which could contribute to more effective utilization of wild germplasm resources. Furthermore, we confirmed that Ae. tauschii was derived from monophyletic speciation rather than hybrid speciation at the cp genome level. We also identified four variable genomic regions, rpl32-trnL-UAG, ccsA-ndhD, rbcL-psaI and rps18-rpl20, showing high levels of nucleotide polymorphisms, which may accordingly prove useful as cpDNA markers in studying the intraspecific genetic structure and diversity of Ae. tauschii.

Unlocking the Genetic Diversity and Population Structure of a Wild Gene Source of Wheat, Aegilops biuncialis Vis., and Its Relationship With the Heading Time

Understanding the genetic diversity of Aegilops biuncialis, a valuable source of agronomical useful genes, may significantly facilitate the introgression breeding of wheat. The genetic diversity and population structure of 86 Ae. biuncialis genotypes were investigated by 32700 DArT markers with the simultaneous application of three statistical methods- neighbor-joining clustering, Principal Coordinate Analysis, and the Bayesian approach to classification. The collection of Ae. biuncialis accessions was divided into five groups that correlated well with their eco-geographic habitat: A (North Africa), B (mainly from Balkans), C (Kosovo and Near East), D (Turkey, Crimea, and Peloponnese), and E (Azerbaijan and the Levant region). The diversity between the Ae. biuncialis accessions for a phenological trait (heading time), which is of decisive importance in the adaptation of plants to different eco-geographical environments, was studied over 3 years. A comparison of the intraspecific variation in the heading time trait by means of analysis of variance and principal component analysis revealed four phenotypic categories showing association with the genetic structure and geographic distribution, except for minor differences. The detailed exploration of genetic and phenologic divergence provides an insight into the adaptation capacity of Ae. biuncialis, identifying promising genotypes that could be utilized for wheat improvement.

Genomic Analysis Confirms Population Structure and Identifies Inter-Lineage Hybrids in Aegilops tauschii

Aegilops tauschii, the D-genome donor of bread wheat, Triticum aestivum, is a storehouse of genetic diversity, and an important resource for future wheat improvement. Genomic and population analysis of 549 Ae. tauschii and 103 wheat accessions was performed by using 13,135 high quality SNPs. Population structure, principal component, and cluster analysis confirmed the differentiation of Ae. tauschii into two lineages; lineage 1 (L1) and lineage 2 (L2), the latter being the wheat D-genome donor. Lineage L1 contributes only 2.7% of the total introgression from Ae. tauschii for a set of United States winter wheat lines, confirming the great amount of untapped genetic diversity in L1. Lineage L2 accessions had overall greater allelic diversity and wheat accessions had the least allelic diversity. Both lineages also showed intra-lineage differentiation with L1 being driven by longitudinal gradient and L2 differentiated by altitude. There has previously been little reported on natural hybridization between L1 and L2. We found nine putative inter-lineage hybrids in the population structure analysis, each containing numerous lineage-specific private alleles from both lineages. One hybrid was confirmed as a recombinant inbred between the two lineages, likely artificially post collection. Of the remaining eight putative hybrids, a group of seven from Georgia carry 713 SNPs with private alleles, which points to the possibility of a novel L1-L2 hybrid lineage. To facilitate the use of Ae. tauschii in wheat improvement, a MiniCore consisting of 29 L1 and 11 L2 accessions, has been developed based on genotypic, phenotypic and geographical data. MiniCore reduces the collection size by over 10-fold and captures 84% of the total allelic diversity in the whole collection.

Developing a High-Throughput SNP-Based Marker System to Facilitate the Introgression of Traits From Aegilops Species Into Bread Wheat (Triticum aestivum)

The genus Aegilops contains a diverse collection of wild species exhibiting variation in geographical distribution, ecological adaptation, ploidy and genome organization. Aegilops is the most closely related genus to Triticum which includes cultivated wheat, a globally important crop that has a limited gene pool for modern breeding. Aegilops species are a potential future resource for wheat breeding for traits, such as adaptation to different ecological conditions and pest and disease resistance. This study describes the development and application of the first high-throughput genotyping platform specifically designed for screening wheat relative species. The platform was used to screen multiple accessions representing all species in the genus Aegilops. Firstly, the data was demonstrated to be useful for screening diversity and examining relationships within and between Aegilops species. Secondly, markers able to characterize and track introgressions from Aegilops species in hexaploid wheat were identified and validated using two different approaches.

Population structure and association analyses of the core collection of hexaploid accessions conserved ex situ in the Japanese gene bank NBRP-Wheat

In this study, we investigated the genetic diversity and population structure of the core collection of hexaploid wheat accessions in the Japanese wheat gene bank NBRP-Wheat. The core collection, consisting of 188 accessions of Triticum aestivum, T. spelta, T. compactum, T. sphaerococcum, T. macha and T. vavilovii, was intensively genotyped by DArTseq markers and consisted of 20,186 SNPs and 60,077 present and absent variations (PAVs). Polymorphic markers were distributed in all chromosomes, with a tendency for smaller numbers on the D-genome chromosomes. We examined the population structure by Bayesian clustering and principal component analysis with a general linear model. Overall, the core collection was divided into seven clusters. Non-admixture accessions in each cluster indicated that the clusters reflect the geographic distribution of the accessions. Both structure analyses strongly suggested that the cluster consisting of T. spelta and T. macha is out-grouped from other hexaploid wheat accessions. We performed genome-wide association analysis pilot studies for nine quantitative and seven qualitative traits and found marker-trait associations for all traits but one, indicating that the current core collection will be useful for detecting uncharacterized QTLs associated with phenotypes of interest.

RNA Sequencing-Based Bulked Segregant Analysis Facilitates Efficient D-genome Marker Development for a Specific Chromosomal Region of Synthetic Hexaploid Wheat

Common wheat originated from interspecific hybridization between cultivated tetraploid wheat and its wild diploid relative Aegilops tauschii followed by amphidiploidization. This evolutionary process can be reproduced artificially, resulting in synthetic hexaploid wheat lines. Here we performed RNA sequencing (RNA-seq)-based bulked segregant analysis (BSA) using a bi-parental mapping population of two synthetic hexaploid wheat lines that shared identical A and B genomes but included with D-genomes of distinct origins. This analysis permitted identification of D-genome-specific polymorphisms around the Net2 gene, a causative locus to hybrid necrosis. The resulting single nucleotide polymorphisms (SNPs) were classified into homoeologous polymorphisms and D-genome allelic variations, based on the RNA-seq results of a parental tetraploid and two Ae. tauschii accessions. The difference in allele frequency at the D-genome-specific SNP sites between the contrasting bulks (ΔSNP-index) was higher on the target chromosome than on the other chromosomes. Several SNPs with the highest ΔSNP-indices were converted into molecular markers and assigned to the Net2 chromosomal region. These results indicated that RNA-seq-based BSA can be applied efficiently to a synthetic hexaploid wheat population to permit molecular marker development in a specific chromosomal region of the D genome.

DArTseq-based analysis of genomic relationships among species of tribe Triticeae

Precise utilization of wild genetic resources to improve the resistance of their cultivated relatives to environmental growth limiting factors, such as salinity stress and diseases, requires a clear understanding of their genomic relationships. Although seriously criticized, analyzing these relationships in tribe Triticeae has largely been based on meiotic chromosome pairing in hybrids of wide crosses, a specialized and labourious strategy. In this study, DArTseq, an efficient genotyping-by-sequencing platform, was applied to analyze the genomes of 34 Triticeae species. We reconstructed the phylogenetic relationships among diploid and polyploid Aegilops and Triticum species, including hexaploid wheat. Tentatively, we have identified the diploid genomes that are likely to have been involved in the evolution of five polyploid species of Aegilops, which have remained unresolved for decades. Explanations which cast light on the progenitor of the A genomes and the complex genomic status of the B/G genomes of polyploid Triticum species in the Emmer and Timopheevi lineages of wheat have also been provided. This study has, therefore, demonstrated that DArTseq genotyping can be effectively applied to analyze the genomes of plants, especially where their genome sequence information are not available.

Fluorescence in situ hybridization karyotyping reveals the presence of two distinct genomes in the taxon Aegilops tauschii

Background

Aegilops tauschii is the donor of the bread wheat D genome. Based on spike morphology, the taxon has conventionally been subdivided into ssp. tauschii and ssp. strangulata. The present study was intended to address the poor match between this whole plant morphology-based subdivision and genetic relationships inferred from genotyping by fluorescence in situ hybridization karyotyping a set of 31 Ae. tauschii accessions.

Results

The distribution of sites hybridizing to the two probes oligo-pTa-535 and (CTT)₁₀ split the Ae. tauschii accessions into two clades, designated D^t and D^s, which corresponded perfectly with a previously assembled phylogeny based on marker genotype. The D^t cluster was populated exclusively by ssp. tauschii accessions, while the D^s cluster harbored both ssp. strangulata and morphologically intermediate accessions. As a result, it is proposed that Ae. tauschii ssp. tauschii is restricted to carriers of the D^t karyotype: their spikelets are regularly spaced along the rachis, at least in the central portion of their spike. Accessions classified as Ae. tauschii ssp. strangulata carry the D^s karyotype; their spikelets are irregularly spaced. Based on this criterion, forms formerly classified as ssp. tauschii var. meyeri have been re-designated ssp. strangulata var. meyeri.

Conclusions

According to the reworking of the taxon, the bread wheat D genome was most probably donated by ssp. strangulata var. meyeri. Chromosomal differentiation reveals intra-species taxon of Ae. tauschii. Ae. tauschii ssp. tauschii has more distant relationship with breed wheat than ssp. strangulata and can be used for breeding improving effectively.

Development and Utilization of Introgression Lines Using Synthetic Octaploid Wheat (Aegilops tauschii × Hexaploid Wheat) as Donor

As the diploid progenitor of common wheat, Aegilops tauschii Cosson (DD, 2n = 2x = 14) is considered to be a promising genetic resource for the improvement of common wheat. In this work, we demonstrated that the efficiency of transferring A. tauschii segments to common wheat was clearly improved through the use of synthetic octaploid wheat (AABBDDDD, 2n = 8x = 56) as a "bridge." The synthetic octaploid was obtained by chromosome doubling of hybrid F₁ (A. tauschii T015 × common wheat Zhoumai 18). A set of introgression lines (BC₁F₈) containing 6016 A. tauschii segments was developed and displayed significant phenotype variance among lines. Twelve agronomic traits, including growth duration, panicle traits, grain traits, and plant height (PH), were evaluated. And transgressive segregation was identified in partial lines. Additionally, better agronomic traits could be observed in some lines, compared to the recurrent parent Zhoumai 18. To verify that the significant variance of those agronomic traits was supposedly controlled by A. tauschii segments, 14 quantitative trait loci (QTLs) for three important agronomic traits (thousand kernel weight, spike length, and PH) were further located in the two environments (Huixian and Zhongmou), indicating the introgression of favorable alleles from A. tauschii into common wheat. This study provides an ameliorated strategy to improve common wheat utilizing a single A. tauschii genome.

The role of reproductive isolation in allopolyploid speciation patterns: empirical insights from the progenitors of common wheat

The ability to cause reproductive isolation often varies among individuals within a plant species. We addressed whether such polymorphism influenced speciation of the allopolyploid common wheat (Triticum aestivum L., AABBDD genome) by evaluating the expression of pre-pollination (outcrossing potential) and post-pollination (crossability) barriers in Aegilops tauschii Coss. (the D genome progenitor). In total, 201 Ae. tauschii accessions representing the entire natural habitat range of the species were used for anther length measurement and artificial crosses with a Triticum turgidum L. (the AB genome progenitor) tester. Intraspecific comparisons showed that both barriers were more strongly expressed in the TauL1 lineage than in the TauL2 lineage. The ability of Ae. tauschii to cause reproductive isolation in the hybridisation with T. turgidum might have markedly influenced common wheat’s speciation by inducing lineage-associated patterns of gene flow. The TauL2 accessions with high potential for natural hybridisation with T. turgidum clustered in the southern coastal Caspian region. This provided phenotypic support for the derivation of the D genome of common wheat from southern Caspian populations. The present study underscored the importance of approaches that incorporate the genealogical and geographic structure of the parental species’ reproductive isolation in understanding the mechanism of plant allopolyploid speciation.

Genetic mapping reveals a dominant awn-inhibiting gene related to differentiation of the variety anathera in the wild diploid wheat Aegilops tauschii

Aegilops tauschii, a wild wheat relative, is the D-genome donor of common wheat. Subspecies and varieties of Ae. tauschii are traditionally classified based on differences in their inflorescence architecture. However, the genetic information for their diversification has been quite limited in the wild wheat relatives. The variety anathera has no awn on the lemma, but the genetic basis for this diagnostic character is unknown. Wide variations in awn length traits at the top and middle spikes were found in the Ae. tauschii core collection, and the awn length at the middle spike was significantly smaller in the eastward-dispersed sublineage than in those in other sublineages. To clarify loci controlling the awnless phenotype of var. anathera, we measured awn length of an intervariety F₂ mapping population, and found that the F₂ individuals could be divided into two groups mainly based on the awn length at the middle of spike, namely short and long awn groups, significantly fitting a 3:1 segregation ratio, which indicated that a single locus controls the awnless phenotype. The awnless locus, Anathera (Antr), was assigned to the distal region of the short arm of chromosome 5D. Quantitative trait locus analysis using the awn length data of each F₂ individual showed that only one major locus was at the same chromosomal position as Antr. These results suggest that a single dominant allele determines the awnless diagnostic character in the variety anathera. The Antr dominant allele is a novel gene inhibiting awn elongation in wheat and its relatives.

Quantitative trait locus analysis for spikelet shape-related traits in wild wheat progenitor Aegilops tauschii: Implications for intraspecific diversification and subspecies differentiation

Wild diploid wheat Aegilops tauschii, the D-genome progenitor of common wheat, carries large genetic variation in spikelet and grain morphology. Two differentiated subspecies of Ae. tauschii, subspecies tauschii and strangulata, have been traditionally defined based on differences in spikelet morphology. Here, we first assessed six spikelet shape-related traits among 199 Ae. tauschii accessions, and found that the accessions belonging to TauL1major lineage produced significantly longer spikes, higher spikelet density, and shorter, narrower spikelets than another major lineage, TauL2, in which the strangulata accessions are included. Next, we performed quantitative trait locus (QTL) analysis of the spikelet and grain shape using three mapping populations derived from interlineage crosses between TauL1 and TauL2 to identify the genetic loci for the morphological variations of the spikelet and grain shape in Ae. tauschii. Three major QTL regions for the examined traits were detected on chromosomes 3D, 4D and 7D. The 3D and 4D QTL regions for several spikelet shape-related traits were conserved in the three mapping populations, which indicated that the 3D and 4D QTLs contribute to divergence of the two major lineages. The 7D QTLs were found only in a mapping population from a cross of the two subspecies, suggesting that these 7D QTLs may be closely related to subspecies differentiation in Ae. tauschii. Thus, QTL analysis for spikelet and grain morphology may provide useful information to elucidate the evolutionary processes of intraspecific differentiation.

Genome-Wide Association Study of Grain Architecture in Wild Wheat Aegilops tauschii

Aegilops tauschii, the D-genome progenitor of Triticum aestivum, encompasses huge diversity for various traits of potential economic importance such as yield, biotic and abiotic stress tolerance, quality and nutrition. In the present study, variation for grain size in Ae. tauschii germplasm was studied and its genetic basis dissected using genome-wide association study (GWAS). Grain length, width, and weight evaluated in 177 Ae. tauschii accessions over 3 years showed near normal distribution with 1.74-, 1.75-, and 2.82-fold variation, respectively. These lines were genetically characterized using genotyping-by-sequencing (GBS) protocol that produced 11,489 single nucleotide polymorphic (SNP) markers. Genetic diversity analysis revealed the presence of two distinct subgroups (designated as lineage 1 and 2) in Ae. tauschii. Based on GBS markers, the genetic similarity was calculated between the accessions and GWAS was conducted using 114 non-redundant accessions and 5,249 SNP markers. A total of 17 SNPs associated with grain size traits distributed over all the seven chromosomes were revealed with 6D, 5D, and 2D harboring most significant marker-trait associations. Some of the chromosomal regions such as 6D_66.4-71.1 cM, 1D_143.5-156.7 cM, and 2D_89.9-92.5 cM had associations with multiple traits. Candidate genes associated with cell division and differentiation were identified for some of the associated SNP markers. Further efforts to validate these loci will help to understand their role in determining grain size and allelic diversity in current germplasm and its effect on grain size upon transfer to bread wheat background.

A large-scale chromosome-specific SNP discovery guideline

Single-nucleotide polymorphisms (SNPs) are the most prevalent type of variation in genomes that are increasingly being used as molecular markers in diversity analyses, mapping and cloning of genes, and germplasm characterization. However, only a few studies reported large-scale SNP discovery in Aegilops tauschii, restricting their potential use as markers for the low-polymorphic D genome. Here, we report 68,592 SNPs found on the gene-related sequences of the 5D chromosome of Ae. tauschii genotype MvGB589 using genomic and transcriptomic sequences from seven Ae. tauschii accessions, including AL8/78, the only genotype for which a draft genome sequence is available at present. We also suggest a workflow to compare SNP positions in homologous regions on the 5D chromosome of Triticum aestivum, bread wheat, to mark single nucleotide variations between these closely related species. Overall, the identified SNPs define a density of 4.49 SNPs per kilobyte, among the highest reported for the genic regions of Ae. tauschii so far. To our knowledge, this study also presents the first chromosome-specific SNP catalog in Ae. tauschii that should facilitate the association of these SNPs with morphological traits on chromosome 5D to be ultimately targeted for wheat improvement.

Complete chloroplast genomes of Aegilops tauschii Coss. and Ae. cylindrica Host sheds light on plasmon D evolution

Hexaploid wheat (Triticum aestivum L., genomes AABBDD) originated in South Caucasus by allopolyploidization of the cultivated Emmer wheat T. dicoccum (genomes AABB) with the Caucasian Ae. tauschii ssp strangulata (genomes DD). Genetic variation of Ae. tauschii is an important natural resource, that is why it is of particular importance to investigate how this variation was formed during Ae. tauschii evolutionary history and how it is presented through the species area. The D genome is also found in tetraploid Ae. cylindrica Host (2n = 28, CCDD). The plasmon diversity that exists in Triticum and Aegilops species is of great significance for understanding the evolution of these genera. In the present investigation the complete nucleotide sequence of plasmon D (chloroplast DNA) of nine accessions of Ae. tauschii and two accessions of Ae. cylindrica are presented. Twenty-eight SNPs are characteristic for both TauL1 and TauL2 accessions of Ae. tauschii using TauL3 as a reference. Four SNPs are additionally observed for TauL2 lineage. The longest (27 bp) indel is located in the intergenic spacer Rps15-ndhF of SSC. This indel can be used for simple determination of TauL3 lineage among Ae. tauschii accessions. In the case of Ae. cylindrica additionally 7 SNPs were observed. The phylogeny tree shows that chloroplast DNA of TauL1 and TauL2 diverged from the TauL3 lineage. TauL1 lineage is relatively older then TauL2. The position of Ae. cylindrica accessions on Ae. tauschii phylogeny tree constructed on chloroplast DNA variation data is intermediate between TauL1 and TauL2. The complete nucleotide sequence of chloroplast DNA of Ae. tauschii and Ae. cylindrica allows to refine the origin and evolution of D plasmon of genus Aegilops.

Intraspecific lineage divergence and its association with reproductive trait change during species range expansion in central Eurasian wild wheat Aegilops tauschii Coss. (Poaceae)

Background

How species ranges form in landscapes is a matter of long-standing evolutionary interest. However, little is known about how natural phenotypic variations of ecologically important traits contribute to species range expansion. In this study, we examined the phylogeographic patterns of phenotypic changes in life history (seed production) and phenological (flowering time) traits during the range expansion of Aegilops tauschii Coss. from the Transcaucasus and Middle East to central Asia.

Results

Our comparative analyses of the patterns of natural variations for those traits and their association with the intraspecific lineage structure showed that (1) the eastward expansion to Asia was driven by an intraspecific sublineage (named TauL1b), (2) high seed production ability likely had an important role at the initial dispersal stage of TauL1b’s expansion to Asia, and (3) the phenological change to early flowering phenotypes was one of the key adaptation events for TauL1b to further expand its range in Asia.

Conclusions

This study provides for the first time a broad picture of the process of Ae. tauschii’s eastward range expansion in which life history and phenological traits may have had respective roles in its dispersal and adaptation in Asia. The clear association of seed production and flowering time patterns with the intraspecific lineage divergence found in this study invites further genetic research to bring the mechanistic understanding of the changes in these key functional traits during range expansion within reach.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0496-9) contains supplementary material, which is available to authorized users.

Fine mapping of shattering locus Br2 reveals a putative chromosomal inversion polymorphism between the two lineages of Aegilops tauschii

This work laid the foundation for cloning of shattering gene Br2 and provided first line of evidence that two major Aegilops tauschii lineages are differentiated by an inversion polymorphism. Chromosome inversions often accompany population differentiation and capture local adaptation during speciation. Aegilops tauschii, the D-genome donor species of hexaploid wheat, consists of two genetically isolated lineages, L1 and L2, but little is known about the genetic mechanisms underlying the population differentiation in this diploid species. During fine mapping of the shattering gene Br2 using a large F2 population derived from a cross between TA1604 (an L1 accession) and AL8/78 (an L2 accession), we found contrasting patterns of crossover distribution in the Br2 interval and neighboring regions despite the high local gene synteny with Brachypodium distachyon and rice. Br2 was localized in a 0.08-cM interval, and 13 marker loci formed a block, where single-crossovers were completely suppressed, but double-crossovers were enriched with a recombination rate of ~11 cM/Mb. In contrast, in a neighboring region no double-crossover was recovered, but single-crossover rate reached 24 cM/Mb, which is much higher than the genome-wide average. This result suggests a putative inversion polymorphism between the parental lines in the Br2 region. Genotyping using the markers from the Br2 region divided a collection of 55 randomly sampled A. tauschii accessions into two major groups, and they are largely genetically isolated. The two groups correspond to the L1 and L2 lineages based on their geographic distribution patterns. This provides first evidence that inversions may underlie the evolution of A. tauschii lineages. The presence of inter-lineage inversions may complicate map-based cloning in A. tauschii and transfer of useful traits to wheat.

The cuticular wax inhibitor locus Iw2 in wild diploid wheat Aegilops tauschii: phenotypic survey, genetic analysis, and implications for the evolution of common wheat

BACKGROUND

Cuticular wax production on plant surfaces confers a glaucous appearance and plays important roles in plant stress tolerance. Most common wheat cultivars, which are hexaploid, and most tetraploid wheat cultivars are glaucous; in contrast, a wild wheat progenitor, Aegilops tauschii, can be glaucous or non-glaucous. A dominant non-glaucous allele, Iw2, resides on the short arm of chromosome 2D, which was inherited from Ae. tauschii through polyploidization. Iw2 is one of the major causal genes related to variation in glaucousness among hexaploid wheat. Detailed genetic and phylogeographic knowledge of the Iw2 locus in Ae. tauschii may provide important information and lead to a better understanding of the evolution of common wheat.

RESULTS

Glaucous Ae. tauschii accessions were collected from a broad area ranging from Armenia to the southwestern coastal part of the Caspian Sea. Linkage analyses with five mapping populations showed that the glaucous versus non-glaucous difference was mainly controlled by the Iw2 locus in Ae. tauschii. Comparative genomic analysis of barley and Ae. tauschii was then used to develop molecular markers tightly linked with Ae. tauschii Iw2. Chromosomal synteny around the orthologous Iw2 regions indicated that some chromosomal rearrangement had occurred during the genetic divergence leading to Ae. tauschii, barley, and Brachypodium. Genetic associations between specific Iw2-linked markers and respective glaucous phenotypes in Ae. tauschii indicated that at least two non-glaucous accessions might carry other glaucousness-determining loci outside of the Iw2 locus.

CONCLUSION

Allelic differences at the Iw2 locus were the main contributors to the phenotypic difference between the glaucous and non-glaucous accessions of Ae. tauschii. Our results supported the previous assumption that the D-genome donor of common wheat could have been any Ae. tauschii variant that carried the recessive iw2 allele.

Identification of quantitative trait loci for abscisic acid responsiveness in the D-genome of hexaploid wheat

In crop species such as wheat, abiotic stresses and preharvest sprouting reduce grain yield and quality. The plant hormone abscisic acid (ABA) plays important roles in abiotic stress tolerance and seed dormancy. In previous studies, we evaluated ABA responsiveness of 67 Aegilops tauschii accessions and their synthetic hexaploid wheat lines, finding wide variation that was due to the D-genome. In this study, quantitative trait locus (QTL) analysis was performed using an F2 population derived from crosses of highly ABA-responsive and less-responsive synthetic wheat lines. A significant QTL was detected on chromosome 6D, in a similar location to that reported for ABA responsiveness using recombinant inbred lines derived from common wheat cultivars Mironovskaya 808 and Chinese Spring. A comparative map and physiological and expression analyses of the 6D QTL suggested that this locus involved in line differences among wheat synthetics is different from that involved in cultivar differences in common wheat. The common wheat 6D QTL was found to affect seed dormancy and the regulation of cold-responsive/late embryogenesis abundant genes during dehydration. However, in synthetic wheat, we failed to detect any association of ABA responsiveness with abiotic stress tolerance or seed dormancy, at least under our experimental conditions. Development of near-isogenic lines will be important for functional analyses of the synthetic wheat 6D QTL.

Genome-wide marker development for the wheat D genome based on single nucleotide polymorphisms identified from transcripts in the wild wheat progenitor Aegilops tauschii

13,347 high-confidence SNPs were discovered through transcriptome sequencing of Aegilops tauschii, which are useful for genomic analysis and molecular breeding of hexaploid wheat. In organisms with large and complex genomes, such as wheat, RNA-seq analysis is cost-effective for discovery of genome-wide single nucleotide polymorphisms (SNPs). In this study, deep sequencing of the spike transcriptome from two Aegilops tauschii accessions representing two major lineages led to the discovery of 13,347 high-confidence (HC) SNPs in 4,872 contigs. After removing redundant SNPs detected in the leaf transcriptome from the same accessions in an earlier study, 10,589 new SNPs were discovered. In total, 5,642 out of 5,808 contigs with HC SNPs were assigned to the Ae. tauschii draft genome sequence. On average, 732 HC polymorphic contigs were mapped in silico to each Ae. tauschii chromosome. Based on the polymorphic data, we developed markers to target the short arm of chromosome 2D and validated the polymorphisms using 20 Ae. tauschii accessions. Of the 29 polymorphic markers, 28 were successfully mapped to 2DS in the diploid F2 population of Ae. tauschii. Among ten hexaploid wheat lines, which included wheat synthetics and common wheat cultivars, 25 of the 43 markers were polymorphic. In the hexaploid F2 population between a common wheat cultivar and a synthetic wheat line, 23 of the 25 polymorphic markers between the parents were available for genotyping of the F2 plants and 22 markers mapped to chromosome 2DS. These results indicate that molecular markers that developed from polymorphisms between two distinct lineages of Ae. tauschii might be useful for analysis not only of the diploid, but also of the hexaploid wheat genome.