Crossability of Triticum urartu and Triticum monococcum wheats, homoeologous recombination, and description of a panel of interspecific introgression lines

Integration of genetic and genomics resources in einkorn wheat enables precision mapping of important traits

Einkorn wheat (Triticum monococcum) is an ancient grain crop and a close relative of the diploid progenitor (T. urartu) of polyploid wheat. It is the only diploid wheat species having both domesticated and wild forms and therefore provides an excellent system to identify domestication genes and genes for traits of interest to utilize in wheat improvement. Here, we leverage genomic advancements for einkorn wheat using an einkorn reference genome assembly combined with skim-sequencing of a large genetic population of 812 recombinant inbred lines (RILs) developed from a cross between a wild and a domesticated T. monococcum accession. We identify 15,919 crossover breakpoints delimited to a median and average interval of 114 Kbp and 219 Kbp, respectively. This high-resolution mapping resource enables us to perform fine-scale mapping of one qualitative (red coleoptile) and one quantitative (spikelet number per spike) trait, resulting in the identification of small physical intervals (400 Kb to 700 Kb) with a limited number of candidate genes. Furthermore, an important domestication locus for brittle rachis is also identified on chromosome 7A. This resource presents an exciting route to perform trait discovery in diploid wheat for agronomically important traits and their further deployment in einkorn as well as tetraploid pasta wheat and hexaploid bread wheat cultivars.

A Combinatorial Q-locus and Tubulin-Based Polymorphism (TBP) Approach Helps in Discriminating Triticum Species

The simple and straightforward recognition of Triticum species is not an easy task due to their complex genetic origins. To provide a recommendation, we have compared the performance of different PCR-based methods relying on the discrimination ability of the Q- and γ-gliadin (GAG56D) genes, as well as TBP (Tubulin-Based Polymorphism), a method based on the multiple amplification of genes of the β-tubulin family. Among these approaches, the PCR-RFLP (Restriction Fragment Length Polymorphism) assay based on a single-nucleotide polymorphism (SNP) present in the Q gene is the only one capable of fully discerning hexaploid spelt and common wheat species, while both γ-gliadin and TBP fail with similar error frequencies. The Q-locus assay results in the attainment of either a single fragment or a doublet, depending on the presence of a suitable restriction site, which is affected by the mutation. This dual pattern of resolution limits both the diagnostic effectiveness, when additional Triticum species are assayed and compared to each other, and its usefulness, when commercially available flours are analyzed. These limitations are overtaken by flanking the Q-locus assay with the TBP analysis. In this way, almost all of the Triticum species can be accurately identified.

Genetic characterization and curation of diploid A-genome wheat species

A-genome diploid wheats represent the earliest domesticated and cultivated wheat species in the Fertile Crescent and include the donor of the wheat A sub-genome. The A-genome species encompass the cultivated einkorn (Triticum monococcum L. subsp. monococcum), wild einkorn (T. monococcum L. subsp. aegilopoides (Link) Thell.), and Triticum urartu. We evaluated the collection of 930 accessions in the Wheat Genetics Resource Center (WGRC) using genotyping by sequencing and identified 13,860 curated single-nucleotide polymorphisms. Genomic analysis detected misclassified and genetically identical (>99%) accessions, with most of the identical accessions originating from the same or nearby locations. About 56% (n = 520) of the WGRC A-genome species collections were genetically identical, supporting the need for genomic characterization for effective curation and maintenance of these collections. Population structure analysis confirmed the morphology-based classifications of the accessions and reflected the species geographic distributions. We also showed that T. urartu is the closest A-genome diploid to the A-subgenome in common wheat (Triticum aestivum L.) through phylogenetic analysis. Population analysis within the wild einkorn group showed three genetically distinct clusters, which corresponded with wild einkorn races α, β, and γ described previously. The T. monococcum genome-wide FST scan identified candidate genomic regions harboring a domestication selection signature at the Non-brittle rachis 1 (Btr1) locus on the short arm of chromosome 3Am at ∼70 Mb. We established an A-genome core set (79 accessions) based on allelic diversity, geographical distribution, and available phenotypic data. The individual species core set maintained at least 79% of allelic variants in the A-genome collection and constituted a valuable genetic resource to improve wheat and domesticated einkorn in breeding programs.

Chromosome-specific KASP markers for detecting Amblyopyrum muticum segments in wheat introgression lines

Many wild-relative species are being used in prebreeding programs to increase the genetic diversity of wheat (Triticum aestivum L.). Genotyping tools such as single nucleotide polymorphism (SNP)-based arrays and molecular markers have been widely used to characterize wheat-wild relative introgression lines. However, due to the polyploid nature of the recipient wheat genome, it is difficult to develop SNP-based Kompetitive allele-specific polymerase chain reaction (KASP) markers that are codominant to track the introgressions from the wild species. Previous attempts to develop KASP markers have involved both exome- and polymerase chain reaction (PCR)-amplicon-based sequencing of the wild species. But chromosome-specific KASP assays have been hindered by homoeologous SNPs within the wheat genome. This study involved whole genome sequencing of the diploid wheat wild relative Amblyopyrum muticum (Boiss.) Eig and development of a de novo SNP discovery pipeline that generated ∼38,000 SNPs in unique wheat genome sequences. New assays were designed to increase the density of Am. muticum polymorphic KASP markers. With a goal of one marker per 60 Mbp, 335 new KASP assays were validated as diagnostic for Am. muticum in a wheat background. Together with assays validated in previous studies, 498 well distributed chromosome-specific markers were used to recharacterize previously genotyped wheat-Am. muticum doubled haploid (DH) introgression lines. The chromosome-specific nature of the KASP markers allowed clarification of which wheat chromosomes were involved with recombination events or substituted with Am. muticum chromosomes and the higher density of markers allowed detection of new small introgressions in these DH lines.

Significance of the Diversification of Wheat Species for the Assembly and Functioning of the Root-Associated Microbiome

Wheat, one of the major crops in the world, has had a complex history that includes genomic hybridizations between Triticum and Aegilops species and several domestication events, which resulted in various wild and domesticated species (especially Triticum aestivum and Triticum durum), many of them still existing today. The large body of information available on wheat-microbe interactions, however, was mostly obtained without considering the importance of wheat evolutionary history and its consequences for wheat microbial ecology. This review addresses our current understanding of the microbiome of wheat root and rhizosphere in light of the information available on pre- and post-domestication wheat history, including differences between wild and domesticated wheats, ancient and modern types of cultivars as well as individual cultivars within a given wheat species. This analysis highlighted two major trends. First, most data deal with the taxonomic diversity rather than the microbial functioning of root-associated wheat microbiota, with so far a bias toward bacteria and mycorrhizal fungi that will progressively attenuate thanks to the inclusion of markers encompassing other micro-eukaryotes and archaea. Second, the comparison of wheat genotypes has mostly focused on the comparison of T. aestivum cultivars, sometimes with little consideration for their particular genetic and physiological traits. It is expected that the development of current sequencing technologies will enable to revisit the diversity of the wheat microbiome. This will provide a renewed opportunity to better understand the significance of wheat evolutionary history, and also to obtain the baseline information needed to develop microbiome-based breeding strategies for sustainable wheat farming.

Introgression of the Powdery Mildew Resistance Genes Pm60 and Pm60b from Triticum urartu to Common Wheat Using Durum as a 'Bridge'

Powdery mildew, caused by the fungus Blumeria graminis f. sp. tritici (Bgt), has limited wheat yields in many major wheat-production areas across the world. Introducing resistance genes from wild relatives into cultivated wheat can enrich the genetic resources for disease resistance breeding. The powdery mildew resistance gene Pm60 was first identified in diploid wild wheat Triticum urartu (T. urartu). In this study, we used durum as a ‘bridge’ approach to transfer Pm60 and Pm60b into hexaploid common wheat. Synthetic hexaploid wheat (SHW, AABBA^uA^u), developed by crossing T. urartu (A^uA^u) with durum (AABB), was used for crossing and backcrossing with common wheat. The Pm60 alleles were tracked by molecular markers and the resistance to powdery mildew. From BC₁F₁ backcross populations, eight recombinant types were identified based on five Pm60-flanking markers, which indicated different sizes of the introgressed chromosome segments from T. urartu. Moreover, we have selected two resistance-harboring introgression lines with high self-fertility, which could be easily used in wheat breeding system. Our results showed that the durum was an excellent ‘bridge’ for introducing the target gene from diploid T. urartu into the hexaploid cultivated wheat. Moreover, these introgression lines could be deployed in wheat resistance breeding programs, together with the assistance of the molecular markers for Pm60 alleles.

Genome-wide association study of agronomic and quality traits in a world collection of the wild wheat relative Triticum urartu

Triticum urartu (2n = 2x = 14, subgenome A^u A^u ), a wild diploid wheat progenitor, features broad allelic diversity for a number of traits of agronomic relevance. A thorough characterization of the diversity of T. urartu natural accessions may provide wheat breeders with new alleles potentially contributing to wheat improvement. In this study, we performed an extensive genotypic and phenotypic characterization of a world collection of 299 T. urartu ex situ accessions, developing 441 327 single nucleotide polymorphisms and recording trait values for agronomic and quality traits. The collection was highly diverse, with broad variation in phenology and plant architecture traits. Seed features were also varied, and analyses of flour quality reported 18 distinct patterns of glutenins, and carotenoid concentrations and sedimentation volumes in some cases surpassing those of cultivated materials. The genome-wide molecular markers developed on the collection were used to conduct a genome-wide association study reporting 25 highly significant quantitative trait nucleotides for the traits under examination, only partially overlapping loci already reported in wheat. Our data show that T. urartu may be considered a valuable allele pool to support the improvement of wheat agronomy and quality.

Development of a New Am -Genome-Specific Single Nucleotide Polymorphism Marker Set for the Molecular Characterization of Wheat-Triticum monococcum Introgression Lines

We identified 1247 polymorphic single nucleotide polymorphisms between Triticum monococcum and wheat. We identified 191 markers validated across all seven chromosomes of T. monococcum. Detected a T. monococcum introgression in leaf-rust-resistant lines. Cultivated einkorn wheat (Triticum monococcum L. subsp. monococcum, 2n = 2x = 14, A^m A^m ) and its wild relative T. monococcum subsp. aegilopoides are important sources of economically useful genes that can be exploited for wheat (Triticum aestivum L.) breeding. Einkorn has excellent resistance to fungal diseases and gene transfer is relatively simple via standard breeding methods. To fulfill the growing demand by modern prebreeding programs for a cost-effective high-throughput procedure for accurately detecting introgressed chromosomes or chromosome segments from T. monococcum into wheat, we used the Axiom Wheat-Relative Genotyping Array and developed a set of A^m genome-specific exome-based single nucleotide polymorphism (SNP) markers suitable for rapid identification of T. monococcum chromatin in a wheat background. We identified 1247 polymorphic SNPs between T. monococcum and wheat. We identified 191 markers across all seven chromosomes of T. monococcum that are also present on an existing Triticum urartu Thum. ex Gandil. genetic map and potentially ordered them on the basis of the high macrocollinearity and conservation of marker order between T. monococcum and T. urartu. The marker set has been tested on leaf-rust-resistant BC₃ F₄ progenies of wheat-T. monococcum hybrids. Two markers (AX-94492165, AX-95073542) placed on the distal end of the chromosome arm 7AL detected a T. monococcum introgression into wheat. The SNP marker set thus proved highly effective in the identification of T. monococcum chromatin in a wheat background, offering a reliable method for screening and selecting wheat-T. monococcum introgression lines, a procedure that could significantly speed up prebreeding programs.

Genome-wide polymorphisms from RNA sequencing assembly of leaf transcripts facilitate phylogenetic analysis and molecular marker development in wild einkorn wheat

A survey of genome-wide polymorphisms between closely related species is required to understand the molecular basis of the evolutionary differentiation of their genomes. Two wild diploid wheat species, namely Triticum monococcum ssp. aegilopoides and T. urartu, are closely related and harbour the A^m and A genomes, respectively. The A-genome donor of tetraploid and common wheat is T. urartu, and T. monococcum ssp. monococcum is the cultivated form derived from the wild einkorn wheat subspecies aegilopoides. Although subspecies aegilopoides has been a useful genetic resource in wheat breeding, genome-wide molecular markers for this subspecies have not been sufficiently developed. Here, we describe the detection of genome-wide polymorphisms such as single-nucleotide polymorphisms (SNPs) and insertions/deletions (indels) from RNA sequencing (RNA-seq) data of leaf transcripts in 15 accessions of the two diploid wheat species. The SNPs and indels, detected using the A genome of common wheat as the reference genome, covered the entire chromosomes of these species. The polymorphism information facilitated a comparison of the genetic diversity of einkorn wheat with that of two related diploid Aegilops species, namely, Ae. tauschii and Ae. umbellulata. Cleaved amplified polymorphic sequence (CAPS) markers converted from the SNP data were efficiently developed to confirm the addition of aegilopoides subspecies chromosomes to tetraploid wheat in nascent allohexaploid lines with AABBA^mA^m genomes. In addition, the CAPS markers permitted linkage map construction in mapping populations of aegilopoides subspecies accessions. Therefore, these RNA-seq data provide information for further breeding of closely related species with no reference genome sequence data.

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.

Detection of T. urartu Introgressions in Wheat and Development of a Panel of Interspecific Introgression Lines

Tritcum urartu (2n = 2x = 14, AuAu), the A genome donor of wheat, is an important source for new genetic variation for wheat improvement due to its high photosynthetic rate and disease resistance. By facilitating the generation of genome-wide introgressions leading to a variety of different wheat-T. urartu translocation lines, T. urartu can be practically utilized in wheat improvement. Previous studies that have generated such introgression lines have been unable to successfully use cytological methods to detect the presence of T. urartu in these lines. Many have, thus, used a variety of molecular markers with limited success due to the low-density coverage of these markers and time-consuming nature of the techniques rendering them unsuitable for large-scale breeding programs. In this study, we report the generation of a resource of single nucleotide polymorphic (SNP) markers, present on a high-throughput SNP genotyping array, that can detect the presence of T. urartu in a hexaploid wheat background making it a potentially valuable tool in wheat pre-breeding programs. A whole genome introgression approach has resulted in the transfer of different chromosome segments from T. urartu into wheat which have then been detected and characterized using these SNP markers. The molecular analysis of these wheat-T. urartu recombinant lines has resulted in the generation of a genetic map of T. urartu containing 368 SNP markers, spread across all seven chromosomes of T. urartu. Comparative analysis of the genetic map of T. urartu and the physical map of the hexaploid wheat genome showed that synteny between the two species is highly conserved at the macro-level and confirmed the presence of the 4/5 translocation in T. urartu also present in the A genome of wheat. A panel of 17 wheat-T. urartu recombinant lines, which consisted of introgressed segments that covered the whole genome of T. urartu, were also selected for self-fertilization to provide a germplasm resource for future trait analysis. This valuable resource of high-density molecular markers specifically designed for detecting wild relative chromosomes and a panel of stable interspecific introgression lines will greatly enhance the efficiency of wheat improvement through wild relative introgressions.

Molecular diversity and landscape genomics of the crop wild relative Triticum urartu across the Fertile Crescent

Modern plant breeding can benefit from the allelic variation that exists in natural populations of crop wild relatives that evolved under natural selection in varying pedoclimatic conditions. In this study, next-generation sequencing was used to generate 1.3 million genome-wide single nucleotide polymorphisms (SNPs) on ex situ collections of Triticum urartu L., the wild donor of the A^u subgenome of modern wheat. A set of 75 511 high-quality SNPs were retained to describe 298 T. urartu accessions collected throughout the Fertile Crescent. Triticum urartu showed a complex pattern of genetic diversity, with two main genetic groups distributed sequentially from west to east. The incorporation of geographical information on sampling points showed that genetic diversity was correlated to the geographical distance (R²  = 0.19) separating samples from Jordan and Lebanon, from Syria and southern Turkey, and from eastern Turkey, Iran and Iraq. The wild emmer genome was used to derive the physical positions of SNPs on the seven chromosomes of the A^u subgenome, allowing us to describe a relatively slow decay of linkage disequilibrium in the collection. Outlier loci were described on the basis of the geographic distribution of the T. urartu accessions, identifying a hotspot of directional selection on chromosome 4A. Bioclimatic variation was derived from grid data and related to allelic variation using a genome-wide association approach, identifying several marker-environment associations (MEAs). Fifty-seven MEAs were associated with altitude and temperature measures while 358 were associated with rainfall measures. The most significant MEAs and outlier loci were used to identify genomic loci with adaptive potential (some already reported in wheat), including dormancy and frost resistance loci. We advocate the application of genomics and landscape genomics on ex situ collections of crop wild relatives to efficiently identify promising alleles and genetic materials for incorporation into modern crop breeding.

Wheat miRNA ancestors: evident by transcriptome analysis of A, B, and D genome donors

MicroRNAs are critical players of post-transcriptional gene regulation with profound effects on the fundamental processes of cellular life. Their identification and characterization, together with their targets, hold great significance in exploring and exploiting their roles on a functional context, providing valuable clues into the regulation of important biological processes, such as stress tolerance or environmental adaptation. Wheat is a hardy crop, extensively harvested in temperate regions, and is a major component of the human diet. With the advent of the next generation sequencing technologies considerably decreasing sequencing costs per base-pair, genomic, and transcriptomic data from several wheat species, including the progenitors and wild relatives have become available. In this study, we performed in silico identification and comparative analysis of microRNA repertoires of bread wheat (Triticum aestivum L.) and its diploid progenitors and relatives, Aegilops sharonensis, Aegilops speltoides, Aegilops tauschii, Triticum monococcum, and Triticum urartu through the utilization of publicly available transcriptomic data. Over 200 miRNA families were identified, majority of which have not previously been reported. Ancestral relationships expanded our understanding of wheat miRNA evolution, while T. monococcum miRNAs delivered important clues on the effects of domestication on miRNA expression. Comparative analyses on wild Ae. sharonensis accessions highlighted candidate miRNAs that can be linked to stress tolerance. The miRNA repertoires of bread wheat and its diploid progenitors and relatives provide important insight into the diversification and distribution of miRNA genes, which should contribute to the elucidation of miRNA evolution of Poaceae family. A thorough understanding of the convergent and divergent expression profiles of miRNAs in different genetic backgrounds can provide unique opportunities to modulation of gene regulation for better crop performance.

Independent mis-splicing mutations in TaPHS1 causing loss of preharvest sprouting (PHS) resistance during wheat domestication

Preharvest sprouting (PHS) is one of the major constraints of wheat production in areas where prolonged rainfall occurs during harvest. TaPHS1 is a gene that regulates PHS resistance on chromosome 3A of wheat, and two causal mutations in the positions +646 and +666 of the TaPHS1 coding region result in wheat PHS susceptibility. Three competitive allele-specific PCR (KASP) markers were developed based on the two mutations in the coding region and one in the promoter region and validated in 82 wheat cultivars with known genotypes. These markers can be used to transfer TaPHS1 in breeding through marker-assisted selection. Screening of 327 accessions of wheat A genome progenitors using the three KASP markers identified different haplotypes in both diploid and tetraploid wheats. Only one Triticum monococcum accession, however, carries both causal mutations in the TaPHS1 coding region and shows PHS susceptibility. Five of 249 common wheat landraces collected from the Fertile Crescent and surrounding areas carried the mutation (C) in the promoter (-222), and one landrace carries both the causal mutations in the TaPHS1 coding region, indicating that the mis-splicing (+646) mutation occurred during common wheat domestication. PHS assay of wheat progenitor accessions demonstrated that the wild-types were highly PHS-resistant, whereas the domesticated type showed increased PHS susceptibility. The mis-splicing TaPHS1 mutation for PHS susceptibility was involved in wheat domestication and might arise independently between T. monococcum and Triticum aestivum.

Comparison of gene expression profiles and responses to zinc chloride among inter- and intraspecific hybrids with growth abnormalities in wheat and its relatives

Hybrid necrosis is a well-known reproductive isolation mechanism in plant species, and an autoimmune response is generally considered to trigger hybrid necrosis through epistatic interaction between disease resistance-related genes in hybrids. In common wheat, the complementary Ne1 and Ne2 genes control hybrid necrosis, defined as type I necrosis. Two other types of hybrid necrosis (type II and type III) have been observed in interspecific hybrids between tetraploid wheat and Aegilops tauschii. Another type of hybrid necrosis, defined here as type IV necrosis, has been reported in F1 hybrids between Triticum urartu and some accessions of Triticum monococcum ssp. aegilopoides. In types I, III and IV, cell death occurs gradually starting in older tissues, whereas type II necrosis symptoms occur only under low temperature. To compare comprehensive gene expression patterns of hybrids showing growth abnormalities, transcriptome analysis of type I and type IV necrosis was performed using a wheat 38k oligo-DNA microarray. Defense-related genes including many WRKY transcription factor genes were dramatically up-regulated in plants showing type I and type IV necrosis, similarly to other known hybrid abnormalities, suggesting an association with an autoimmune response. Reactive oxygen species generation and necrotic cell death were effectively inhibited by ZnCl2 treatment in types I, III and IV necrosis, suggesting a significant association of Ca(2+) influx in upstream signaling of necrotic cell death in wheat hybrid necrosis.

(GAA)n microsatellite as an indicator of the A genome reorganization during wheat evolution and domestication

Although the wheat A genomes have been intensively studied over past decades, many questions concerning the mechanisms of their divergence and evolution still remain unsolved. In the present study we performed comparative analysis of the A genome chromosomes in diploid (Triticum urartu Tumanian ex Gandilyan, 1972, Triticum boeoticum Boissier, 1874 and Triticum monococcum Linnaeus, 1753) and polyploid wheat species representing two evolutionary lineages, Timopheevi (Triticum timopheevii (Zhukovsky) Zhukovsky, 1934 and Triticum zhukovskyi Menabde & Ericzjan, 1960) and Emmer (Triticum dicoccoides (Körnicke ex Ascherson & Graebner) Schweinfurth, 1908, Triticum durum Desfontaines, 1798, and Triticum aestivum Linnaeus, 1753) using a new cytogenetic marker - the pTm30 probe cloned from Triticum monococcum genome and containing (GAA)56 microsatellite sequence. Up to four pTm30 sites located on 1AS, 5AS, 2AS, and 4AL chromosomes have been revealed in the wild diploid species, although most accessions contained one-two (GAA)n sites. The domesticated diploid species Triticum monococcum differs from the wild diploid species by almost complete lack of polymorphism in the distribution of (GAA)n site. Only one (GAA)n site in the 4AL chromosome has been found in Triticum monococcum. Among three wild emmer (Triticum dicoccoides) accessions we detected 4 conserved and 9 polymorphic (GAA)n sites in the A genome. The (GAA)n loci on chromosomes 2AS, 4AL, and 5AL found in of Triticum dicoccoides were retained in Triticum durum and Triticum aestivum. In species of the Timopheevi lineage, the only one, large (GAA)n site has been detected in the short arm of 6A(t) chromosome. (GAA)n site observed in Triticum monococcum are undetectable in the A(b) genome of Triticum zhukovskyi, this site could be eliminated over the course of amphiploidization, while the species was established. We also demonstrated that changes in the distribution of (GAA)n sequence on the A-genome chromosomes of diploid and polyploid wheats are associated with chromosomal rearrangements/ modifications, involving mainly the NOR (nucleolus organizer region)-bearing chromosomes, that took place during the evolution of wild and domesticated species.