The evolution of polyploid wheats: identification of the A genome donor species

Chromosome-scale genome assembly of bread wheat's wild relative Triticum timopheevii

Wheat (Triticum aestivum) is one of the most important food crops with an urgent need for increase in its production to feed the growing world. Triticum timopheevii (2n = 4x = 28) is an allotetraploid wheat wild relative species containing the At and G genomes that has been exploited in many pre-breeding programmes for wheat improvement. In this study, we report the generation of a chromosome-scale reference genome assembly of T. timopheevii accession PI 94760 based on PacBio HiFi reads and chromosome conformation capture (Hi-C). The assembly comprised a total size of 9.35 Gb, featuring a contig N50 of 42.4 Mb and included the mitochondrial and plastid genome sequences. Genome annotation predicted 166,325 gene models including 70,365 genes with high confidence. DNA methylation analysis showed that the G genome had on average more methylated bases than the At genome. In summary, the T. timopheevii genome assembly provides a valuable resource for genome-informed discovery of agronomically important genes for food security.

Allopolyploid subgenome identification and implications for evolutionary analysis

Whole-genome duplications (WGDs) are widespread genomic events in eukaryotes that are hypothesized to contribute to the evolutionary success of many lineages, including flowering plants, Saccharomyces yeast, and vertebrates. WGDs generally can be classified into autopolyploids (ploidy increase descended from one species) or allopolyploids (ploidy increase descended from multiple species). Assignment of allopolyploid progenitor species (called subgenomes in the polyploid) is important to understanding the biology and evolution of polyploids, including the asymmetric subgenome evolution following hybridization (biased fractionation). Here, I review the different methodologies used to identify the ancestors of allopolyploid subgenomes, discuss the advantages and disadvantages of these methods, and outline the implications of how these methods affect the subsequent evolutionary analysis of these genomes.

Enrichment and Diversification of the Wheat Genome via Alien Introgression

Wheat, including durum and common wheat, respectively, is an allopolyploid with two or three homoeologous subgenomes originating from diploid wild ancestral species. The wheat genome's polyploid origin consisting of just three diploid ancestors has constrained its genetic variation, which has bottlenecked improvement. However, wheat has a large number of relatives, including cultivated crop species (e.g., barley and rye), wild grass species, and ancestral species. Moreover, each ancestor and relative has many other related subspecies that have evolved to inhabit specific geographic areas. Cumulatively, they represent an invaluable source of genetic diversity and variation available to enrich and diversify the wheat genome. The ancestral species share one or more homologous genomes with wheat, which can be utilized in breeding efforts through typical meiotic homologous recombination. Additionally, genome introgressions of distant relatives can be moved into wheat using chromosome engineering-based approaches that feature induced meiotic homoeologous recombination. Recent advances in genomics have dramatically improved the efficacy and throughput of chromosome engineering for alien introgressions, which has served to boost the genetic potential of the wheat genome in breeding efforts. Here, we report research strategies and progress made using alien introgressions toward the enrichment and diversification of the wheat genome in the genomics era.

Genome-Wide Identification and Expression Analysis of Catalase Gene Families in Triticeae

Aerobic metabolism in plants results in the production of hydrogen peroxide (H2O2), a significant and comparatively stable non-radical reactive oxygen species (ROS). H2O2 is a signaling molecule that regulates particular physiological and biological processes (the cell cycle, photosynthesis, plant growth and development, and plant responses to environmental challenges) at low concentrations. Plants may experience oxidative stress and ultimately die from cell death if excess H2O2 builds up. Triticum dicoccoides, Triticum urartu, and Triticum spelta are different ancient wheat species that present different interesting characteristics, and their importance is becoming more and more clear. In fact, due to their interesting nutritive health, flavor, and nutritional values, as well as their resistance to different parasites, the cultivation of these species is increasingly important. Thus, it is important to understand the mechanisms of plant tolerance to different biotic and abiotic stresses by studying different stress-induced gene families such as catalases (CAT), which are important H2O2-metabolizing enzymes found in plants. Here, we identified seven CAT-encoding genes (TdCATs) in Triticum dicoccoides, four genes in Triticum urartu (TuCATs), and eight genes in Triticum spelta (TsCATs). The accuracy of the newly identified wheat CAT gene members in different wheat genomes is confirmed by the gene structures, phylogenetic relationships, protein domains, and subcellular location analyses discussed in this article. In fact, our analysis showed that the identified genes harbor the following two conserved domains: a catalase domain (pfam00199) and a catalase-related domain (pfam06628). Phylogenetic analyses showed that the identified wheat CAT proteins were present in an analogous form in durum wheat and bread wheat. Moreover, the identified CAT proteins were located essentially in the peroxisome, as revealed by in silico analyses. Interestingly, analyses of CAT promoters in those species revealed the presence of different cis elements related to plant development, maturation, and plant responses to different environmental stresses. According to RT-qPCR, Triticum CAT genes showed distinctive expression designs in the studied organs and in response to different treatments (salt, heat, cold, mannitol, and ABA). This study completed a thorough analysis of the CAT genes in Triticeae, which advances our knowledge of CAT genes and establishes a framework for further functional analyses of the wheat gene family.

Rht12b, a widely used ancient allele of TaGA2oxA13, reduces plant height and enhances yield potential in wheat

KEY MESSAGE

We identified a new wheat dwarfing allele Rht12b conferring reduced height and higher grain yield, pinpointed its causal variations, developed a breeding-applicable marker, and traced its origin and worldwide distribution. Plant height control is essential to optimize lodging resistance and yield gain in crops. RHT12 is a reduced height (Rht) locus that is identified in a mutationally induced dwarfing mutant and encodes a gibberellin 2-oxidase TaGA2oxA13. However, the artificial dwarfing allele is not used in wheat breeding due to excessive height reduction. Here, we confirmed a stable Rht locus, overlapping with RHT12, in a panel of wheat cultivars and its dwarfing allele reduced plant height by 5.4-8.2 cm, equivalent to Rht12b, a new allele of RHT12. We validated the effect of Rht12b on plant height in a bi-parent mapping population. Importantly, wheat cultivars carrying Rht12b had higher grain yield than those with the contrasting Rht12a allele. Rht12b conferred higher expression level of TaGA2oxA13. Transient activation assays defined SNP-390(C/A) in the promoter of TaGA2oxA13 as the causal variation. An efficient kompetitive allele-specific PCR marker was developed to diagnose Rht12b. Conjoint analysis showed that Rht12b plus the widely used Rht-D1b, Rht8 and Rht24b was the predominant Rht combination and conferred a moderate plant height in tested wheat cultivars. Evolutionary tracking uncovered that RHT12 locus arose from a tandem duplication event with Rht12b firstly appearing in wild emmer. The frequency of Rht12b was approximately 70% (700/1005) in a worldwide wheat panel and comparable to or higher than those of other widely used Rht genes, suggesting it had been subjected to positive selection. These findings not only identify a valuable Rht gene for wheat improvement but also develop a functionally diagnostic tool for marker-assisted breeding.

Molecular cytological analysis of alien introgressions in common wheat lines created by crossing of Triticum aestivum with T. dicoccoides and T. dicoccum

Wild and domesticated emmer (ВВАА, 2n = 28) are of significant interest for expanding the genetic diversity of common wheat as sources of a high protein and microelement grain content, resistance to many biotic and abiotic factors. Particular interest in these species is also determined by their close relationship with Triticum aestivum L., which facilitates interspecific hybridization. The objective of this work was to analyze the nature of alien introgressions in hybrid lines from crossing common wheat varieties with T. dicoccoides and T. dicoccum, and to assess the effect of their genome fragments on the cytological stability of introgression lines. A C-banding technique and genotyping with SNP and SSR markers were used to determine localization and length of introgression fragments. Assessment of cytological stability was carried out on the basis of chromosome behavior in microsporogenesis. A molecular cytogenetic analysis of introgression wheat lines indicated that the inclusion of the genetic material of wild and domesticated emmer was carried out mainly in the form of whole arms or large fragments in the chromosomes of the B genome and less extended inserts in the A genome. At the same time, the highest frequency of introgressions of the emmer genome was observed in chromosomes 1A, 1B, 2B, and 3B. The analysis of the final stage of meiosis showed a high level of cytological stability in the vast majority of introgression wheat lines (meiotic index was 83.0-99.0 %), which ensures the formation of functional gametes in an amount sufficient for successful reproduction. These lines are of interest for the selection of promising material with agronomically valuable traits and their subsequent inclusion in the breeding process.

New insights on the evolution of nucleolar dominance in newly resynthesized hexaploid wheat Triticum zhukovskyi

Nucleolar dominance (ND) is a widespread epigenetic phenomenon in hybridizations where nucleolus transcription fails at the nucleolus organizer region (NOR). However, the dynamics of NORs during the formation of Triticum zhukovskyi (GGAu Au Am Am ), another evolutionary branch of allohexaploid wheat, remains poorly understood. Here, we elucidated genetic and epigenetic changes occurring at the NOR loci within the Am , G, and D subgenomes during allopolyploidization by synthesizing hexaploid wheat GGAu Au Am Am and GGAu Au DD. In T. zhukovskyi, Au genome NORs from T. timopheevii (GGAu Au ) were lost, while the second incoming NORs from T. monococcum (Am Am ) were retained. Analysis of the synthesized T. zhukovskyi revealed that rRNA genes from the Am genome were silenced in F1 hybrids (GAu Am ) and remained inactive after genome doubling and subsequent self-pollinations. We observed increased DNA methylation accompanying the inactivation of NORs in the Am genome and found that silencing of NORs in the S1 generation could be reversed by a cytidine methylase inhibitor. Our findings provide insights into the ND process during the evolutionary period of T. zhukovskyi and highlight that inactive rDNA units may serve as a 'first reserve' in the form of R-loops, contributing to the successful evolution of T. zhukovskyi.

Exploring genetic diversity of wild and related tetraploid wheat species Triticum turgidum and Triticum timopheevii

INTRODUCTION

The domestication bottleneck has reduced genetic diversity inwheat, necessitating the use of wild relatives in breeding programs. Wild tetraploid wheat are widely used in the breeding programs but with morphological characters, it is difficult to distinguish these, resulting in misclassification/mislabeling or duplication of accessions in the Gene bank.

OBJECTIVES

The study aims to exploreGenotyping by sequencing (GBS) to characterize wild and domesticated tetraploid wheat accessions to generate a core set of accessions to be used in the breeding program.

METHODS

TASSEL-GBS pipeline was used for SNP discovery, fastStructure was used to determine the population structure and PowerCore was used to generate a core sets. Nucleotide diversity matrices of Nie's and F-statistics (FST) index were used to determine the center of genetic diversity.

RESULTS

We found 65 % and 47 % duplicated accessions in Triticum timopheevii and T. turgidum respectively. Genome-wide nucleotide diversity and FST scan uncovered a lower intra and higher inter-species differentiation. Distinct FST regions were identified in genomic regions belonging to domestication genes: non-brittle rachis (Btr1) and vernalization (VRN-1).Our results suggest that Israel, Jordan, Syria, and Lebanonas the hub of genetic diversity of wild emmer;Turkey, and Georgia for T. durum; and Iraq, Azerbaijan, and Armenia for theT. timopheevii. Identified core set accessions preserved more than 93 % of the available genetic diversity. Genome wide association study (GWAS) indicated the potential chromosomal segment for resistance to leaf rust in T. timopheevii.

CONCLUSION

The present study explored the potential of GBS technology in data reduction while maintaining the significant genetic diversity of the species. Wild germplasm showed more differentiation than domesticated accessions, indicating the availability of sufficient diversity for crop improvement. With reduced complexity, the core set preserves the genetic diversity of the gene bank collections and will aid in a more robust characterization of wild germplasm.

Differential alternative polyadenylation of homoeologous genes of allohexaploid wheat ABD subgenomes during drought stress response

Because allohexaploid wheat genome contains ABD subgenomes, how the expression of homoeologous genes is coordinated remains largely unknown, particularly at the co-transcriptional level. Alternative polyadenylation (APA) is an important part of co-transcriptional regulation, which is crucial in developmental processes and stress responses. Drought stress is a major threat to the stable yield of wheat. Focusing on APA, we used poly(A) tag sequencing to track poly(A) site dynamics in wheat under drought stress. The results showed that drought stress led to extensive APA involving 37-47% of differentially expressed genes in wheat. Significant poly(A) site switching was found in stress-responsive genes. Interestingly, homoeologous genes exhibit unequal numbers of poly(A) sites, divergent APA patterns with tissue specificity and time-course dynamics, and distinct 3'-UTR length changes. Moreover, differentially expressed transcripts in leaves and roots used different poly(A) signals, the up- and downregulated isoforms had distinct preferences for non-canonical poly(A) sites. Genes that encode key polyadenylation factors showed differential expression patterns under drought stress. In summary, poly(A) signals and the changes in core poly(A) factors may widely affect the selection of poly(A) sites and gene expression levels during the response to drought stress, and divergent APA patterns among homoeologous genes add extensive plasticity to this responsive network. These results not only reveal the significant role of APA in drought stress response, but also provide a fresh perspective on how homoeologous genes contribute to adaptability through transcriptome diversity. In addition, this work provides information about the ends of transcripts for a better annotation of the wheat genome.

Genome sequencing of Sitopsis species provides insights into their contribution to the B subgenome of bread wheat

Wheat (Triticum aestivum, BBAADD) is an allohexaploid species that originated from two polyploidization events. The progenitors of the A and D subgenomes have been identified as Triticum urartu and Aegilops tauschii, respectively. Current research suggests that Aegilops speltoides is the closest but not the direct ancestor of the B subgenome. However, whether Ae. speltoides has contributed genomically to the wheat B subgenome and which chromosome regions are conserved between Ae. speltoides and the B subgenome remain unclear. Here, we assembled a high-quality reference genome for Ae. speltoides, resequenced 53 accessions from seven species (Aegilops bicornis, Aegilops longissima, Aegilops searsii, Aegilops sharonensis, Ae. speltoides, Aegilops mutica [syn. Amblyopyrum muticum], and Triticum dicoccoides) and revealed their genomic contributions to the wheat B subgenome. Our results showed that centromeric regions were particularly conserved between Aegilops and Triticum and revealed 0.17 Gb of conserved blocks between Ae. speltoides and the B subgenome. We classified five groups of conserved and non-conserved genes between Aegilops and Triticum, revealing their biological characteristics, differentiation in gene expression patterns, and collinear relationships between Ae. speltoides and the wheat B subgenome. We also identified gene families that expanded in Ae. speltoides during its evolution and 789 genes specific to Ae. speltoides. These genes can serve as genetic resources for improvement of adaptability to biotic and abiotic stress. The newly constructed reference genome and large-scale resequencing data for Sitopsis species will provide a valuable genomic resource for wheat genetic improvement and genomic studies.

Genetic dissection of the roles of β-hydroxylases in carotenoid metabolism, photosynthesis, and plant growth in tetraploid wheat (Triticum turgidum L.)

Functional redundancy and subfunctionalization of β-hydroxylases in tetraploid wheat tissues open up opportunities for manipulation of carotenoid metabolism for trait improvement. The genetic diversity provided by subgenome homoeologs in allopolyploid wheat can be leveraged for developing improved wheat varieties with modified chemical traits, including profiles of carotenoids, which play critical roles in photosynthesis, photoprotection, and growth regulation. Carotenoid profiles are greatly influenced by hydroxylation catalyzed by β-hydroxylases (HYDs). To genetically dissect the contribution of HYDs to carotenoid metabolism and wheat growth and yield, we isolated loss-of-function mutants of the two homoeologs of HYD1 (HYD-A1 and HYD-B1) and HYD2 (HYD-A2 and HYD-B2) from the sequenced ethyl methanesulfonate mutant population of the tetraploid wheat cultivar Kronos, and generated various mutant combinations. Although functional redundancy between HYD1 and HYD2 paralogs was observed in leaves, HYD1 homoeologs contributed more than HYD2 homoeologs to carotenoid β-ring hydroxylation in this tissue. By contrast, the HYD2 homoeologs functioned toward production of lutein, the major carotenoid in mature grains, whereas HYD1 homoeologs had a limited role. These results collectively suggested subfunctionalization of HYD genes and homoeologs in different tissues of tetraploid wheat. Despite reduced photoprotective responses observed in the triple hyd-A1 hyd-B1 hyd-A2 and the quadruple hyd-A1 hyd-B1 hyd-A2 hyd-B2 combinatorial mutants, comprehensive plant phenotyping analysis revealed that all mutants analyzed were comparable to the control for growth, yield, and fertility, except for a slight delay in anthesis and senescence as well as accelerated germination in the quadruple mutant. Overall, this research takes steps toward untangling the functions of HYDs in wheat and has implications for improving performance and consumer traits of this economically important global crop.

Identification of Fusarium Head Blight Resistance in Triticum timopheevii Accessions and Characterization of Wheat-T. timopheevii Introgression Lines for Enhanced Resistance

A diverse panel of wheat wild relative species was screened for resistance to Fusarium head blight (FHB) by spray inoculation. The great majority of species and accessions were susceptible or highly susceptible to FHB. Accessions of Triticum timopheevii (P95-99.1-1), Agropyron desertorum (9439957), and Elymus vaillantianus (531552) were highly resistant to FHB while additional accessions of T. timopheevii were found to be susceptible to FHB. A combination of spray and point inoculation assessments over two consecutive seasons indicated that the resistance in accession P95-99.1-1 was due to enhanced resistance to initial infection of the fungus (type 1 resistance), and not to reduction in spread (type 2 resistance). A panel of wheat-T. timopheevii (accession P95-99.1-1) introgression lines was screened for FHB resistance over two consecutive seasons using spray inoculation. Most introgression lines were similar in susceptibility to FHB as the wheat recipient (Paragon) but substitution of the terminal portion of chromosome 3BS of wheat with a similar-sized portion of 3G of T. timopheevii significantly enhanced FHB resistance in the wheat background.

Assessing the Heat Tolerance of Meiosis in Spanish Landraces of Tetraploid Wheat Triticum turgidum

Heat stress alters the number and distribution of meiotic crossovers in wild and cultivated plant species. Hence, global warming may have a negative impact on meiosis, fertility, and crop productions. Assessment of germplasm collections to identify heat-tolerant genotypes is a priority for future crop improvement. Durum wheat, Triticum turgidum, is an important cultivated cereal worldwide and given the genetic diversity of the durum wheat Spanish landraces core collection, we decided to analyse the heat stress effect on chiasma formation in a sample of 16 landraces of T. turgidum ssp. turgidum and T. turgidum ssp. durum, from localities with variable climate conditions. Plants of each landrace were grown at 18–22 °C and at 30 °C during the premeiotic temperature-sensitive stage. The number of chiasmata was not affected by heat stress in three genotypes, but decreased by 0.3–2 chiasmata in ten genotypes and more than two chiasmata in the remaining three ones. Both thermotolerant and temperature-sensitive genotypes were found in the two subspecies, and in some of the agroecological zones studied, which supports that genotypes conferring a heat tolerant meiotic phenotype are not dependent on subspecies or geographical origin. Implications of heat adaptive genotypes in future research and breeding are discussed.

Transfer of the ph1b Deletion Chromosome 5B From Chinese Spring Wheat Into a Winter Wheat Line and Induction of Chromosome Rearrangements in Wheat-Aegilops biuncialis Hybrids

Effective utilization of genetic diversity in wild relatives to improve wheat requires recombination between wheat and alien chromosomes. However, this is suppressed by the Pairing homoeologous gene, Ph1, on the long arm of wheat chromosome 5B. A deletion mutant of the Ph1 locus (ph1b) has been used widely to induce homoeologous recombination in wheat × alien hybrids. However, the original ph1b mutation, developed in Chinese Spring (CS) background has poor agronomic performance. Hence, alien introgression lines are first backcrossed with adapted wheat genotypes and after this step, alien chromosome segments are introduced into breeding lines. In this work, the ph1b mutation was transferred from two CSph1b mutants into winter wheat line Mv9kr1. Homozygous genotypes Mv9kr1 ph1b/ph1b exhibited improved plant and spike morphology compared to Chinese Spring. Flow cytometric chromosome analysis confirmed reduced DNA content of the mutant 5B chromosome in both wheat genotype relative to the wild type chromosome. The ph1b mutation in the Mv9kr1 genotype allowed wheat-alien chromosome pairing in meiosis of Mv9kr1ph1b_K × Aegilops biuncialis F₁ hybrids, predominantly with the M^b-genome chromosomes of Aegilops relative to those of the U^b genome. High frequency of wheat-Aegilops chromosome interactions resulted in rearranged chromosomes identified in the new Mv9kr1ph1b × Ae. Biuncialis amphiploids, making these lines valuable sources for alien introgressions. The new Mv9kr1ph1b mutant genotype is a unique resource to support alien introgression breeding of hexaploid wheat.

Whole-genome resequencing of the wheat A subgenome progenitor Triticum urartu provides insights into its demographic history and geographic adaptation

Triticum urartu is the progenitor of the A subgenome in tetraploid and hexaploid wheat. Uncovering the landscape of genetic variations in T. urartu will help us understand the evolutionary and polyploid characteristics of wheat. Here, we investigated the population genomics of T. urartu by genome-wide sequencing of 59 representative accessions collected around the world. A total of 42.2 million high-quality single-nucleotide polymorphisms and 3 million insertions and deletions were obtained by mapping reads to the reference genome. The ancient T. urartu population experienced a significant reduction in effective population size (Ne) from ∼3 000 000 to ∼140 000 and subsequently split into eastern Mediterranean coastal and Mesopotamian-Transcaucasian populations during the Younger Dryas period. A map of allelic drift paths displayed splits and mixtures between different geographic groups, and a strong genetic drift towards hexaploid wheat was also observed, indicating that the direct donor of the A subgenome originated from northwestern Syria. Genetic changes were revealed between the eastern Mediterranean coastal and Mesopotamian-Transcaucasian populations in genes orthologous to those regulating plant development and stress responses. A genome-wide association study identified two single-nucleotide polymorphisms in the exonic regions of the SEMI-DWARF 37 ortholog that corresponded to the different T. urartu ecotype groups. Our study provides novel insights into the origin and genetic legacy of the A subgenome in polyploid wheat and contributes a gene repertoire for genomics-enabled improvements in wheat breeding.

Ancient Caucasian Wheats: A Contribution for Sustainable Diets and Food Diversity

Through the centuries, the domestication and modern breeding of wheat led to a significant loss of genetic variation in the cultivated gene pool with a consequent decrease in food diversity. Current trends towards low-input and sustainable agriculture call for the revitalization and exploitation of ancient wheats, which represent a reservoir of biodiversity useful to ensure sustainable wheat production in the context of climate change and low-input farming systems. Ancient Caucasian wheat species, such as the hulled wheats Triticum timopheevii (tetraploid A^uA^uGG) and Triticum zhukovskyi (hexaploid A^uA^uA^mA^mGG), are still grown to a limited extent in the Caucasus for the production of traditional foods. These Caucasian wheats were grown in Italy and were analyzed for physical, nutritional and technological characteristics and compared to durum wheat. Both Caucasian species revealed a high protein content (on average 18.5%) associated with a low gluten index, mainly in T. zhukovskyi, and test weight values comparable to commercial wheats. The total antioxidant capacity was revealed to be the double of that in durum wheat, suggesting the use of ancient Caucasian wheats for the production of healthy foods. Finally, the technological and rheological results indicated that Caucasian wheats could be potential raw material for the formulation of flat breads, biscuits and pasta.

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.

Genome sequences of five Sitopsis species of Aegilops and the origin of polyploid wheat B subgenome

Common wheat (Triticum aestivum, BBAADD) is a major staple food crop worldwide. The diploid progenitors of the A and D subgenomes have been unequivocally identified; that of B, however, remains ambiguous and controversial but is suspected to be related to species of Aegilops, section Sitopsis. Here, we report the assembly of chromosome-level genome sequences of all five Sitopsis species, namely Aegilops bicornis, Ae. longissima, Ae. searsii, Ae. sharonensis, and Ae. speltoides, as well as the partial assembly of the Amblyopyrum muticum (synonym Aegilops mutica) genome for phylogenetic analysis. Our results reveal that the donor of the common wheat B subgenome is a distinct, and most probably extinct, diploid species that diverged from an ancestral progenitor of the B lineage to which the still extant Ae. speltoides and Am. muticum belong. In addition, we identified interspecific genetic introgressions throughout the evolution of the Triticum/Aegilops species complex. The five Sitopsis species have various assembled genome sizes (4.11-5.89 Gb) with high proportions of repetitive sequences (85.99%-89.81%); nonetheless, they retain high collinearity with other genomes or subgenomes of species in the Triticum/Aegilops complex. Differences in genome size were primarily due to independent post-speciation amplification of transposons. We also identified a set of Sitopsis genes pertinent to important agronomic traits that can be harnessed for wheat breeding. These newly assembled genome resources provide a new roadmap for evolutionary and genetic studies of the Triticum/Aegilops complex, as well as for wheat improvement.

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.

Improvement and Re-Evolution of Tetraploid Wheat for Global Environmental Challenge and Diversity Consumption Demand

Allotetraploid durum wheat is the second most widely cultivated wheat, following hexaploid bread wheat, and is one of the major protein and calorie sources of the human diet. However, durum wheat is encountered with a severe grain yield bottleneck due to the erosion of genetic diversity stemming from long-term domestication and especially modern breeding programs. The improvement of yield and grain quality of durum wheat is crucial when confronted with the increasing global population, changing climate environments, and the non-ignorable increasing incidence of wheat-related disorders. This review summarized the domestication and evolution process and discussed the durum wheat re-evolution attempts performed by global researchers using diploid einkorn, tetraploid emmer wheat, hexaploid wheat (particularly the D-subgenome), etc. In addition, the re-evolution of durum wheat would be promoted by the genetic enrichment process, which could diversify allelic combinations through enhancing chromosome recombination (pentaploid hybridization or pairing of homologous chromosomes gene Ph mutant line induced homoeologous recombination) and environmental adaptability via alien introgressive genes (wide cross or distant hybridization followed by embryo rescue), and modifying target genes or traits by molecular approaches, such as CRISPR/Cas9 or RNA interference (RNAi). A brief discussion of the future perspectives for exploring germplasm for the modern improvement and re-evolution of durum wheat is included.

Detecting major introgressions in wheat and their putative origins using coverage analysis

Introgressions from crop wild relatives (CWRs) have been used to introduce beneficial traits into cultivated plants. Introgressions have traditionally been detected using cytological methods. Recently, single nucleotide polymorphism (SNP)-based methods have been proposed to detect introgressions in crosses for which both parents are known. However, for unknown material, no method was available to detect introgressions and predict the putative donor species. Here, we present a method to detect introgressions and the putative donor species. We demonstrate the utility of this method using 10 publicly available wheat genome sequences and identify nine major introgressions. We show that the method can distinguish different introgressions at the same locus. We trace introgressions to early wheat cultivars and show that natural introgressions were utilised in early breeding history and still influence elite lines today. Finally, we provide evidence that these introgressions harbour resistance genes.

Detecting major introgressions in wheat and their putative origins using coverage analysis

Introgressions from crop wild relatives (CWRs) have been used to introduce beneficial traits into cultivated plants. Introgressions have traditionally been detected using cytological methods. Recently, single nucleotide polymorphism (SNP)-based methods have been proposed to detect introgressions in crosses for which both parents are known. However, for unknown material, no method was available to detect introgressions and predict the putative donor species. Here, we present a method to detect introgressions and the putative donor species. We demonstrate the utility of this method using 10 publicly available wheat genome sequences and identify nine major introgressions. We show that the method can distinguish different introgressions at the same locus. We trace introgressions to early wheat cultivars and show that natural introgressions were utilised in early breeding history and still influence elite lines today. Finally, we provide evidence that these introgressions harbour resistance genes.

Introgression of the Triticum timopheevii Genome Into Wheat Detected by Chromosome-Specific Kompetitive Allele Specific PCR Markers

Triticum timopheevii (2n = 28, A ^t A ^t GG) is a tetraploid wild relative species with great potential to increase the genetic diversity of hexaploid wheat Triticum aestivum (2n = 42, AABBDD) for various important agronomic traits. A breeding scheme that propagated advanced backcrossed populations of wheat-T. timopheevii introgression lines through further backcrossing and self-fertilisation resulted in the generation of 99 introgression lines (ILs) that carried 309 homozygous segments from the A ^t and G subgenomes of T. timopheevii. These introgressions contained 89 and 74 unique segments from the A ^t and G subgenomes, respectively. These overlapping segments covered 98.9% of the T. timopheevii genome that has now been introgressed into bread wheat cv. Paragon including the entirety of all T. timopheevii chromosomes via varying sized segments except for chromosomes 3A ^t , 4G, and 6G. Homozygous ILs contained between one and eight of these introgressions with an average of three per introgression line. These homozygous introgressions were detected through the development of a set of 480 chromosome-specific Kompetitive allele specific PCR (KASP) markers that are well-distributed across the wheat genome. Of these, 149 were developed in this study based on single nucleotide polymorphisms (SNPs) discovered through whole genome sequencing of T. timopheevii. A majority of these KASP markers were also found to be T. timopheevii subgenome specific with 182 detecting A ^t subgenome and 275 detecting G subgenome segments. These markers showed that 98% of the A ^t segments had recombined with the A genome of wheat and 74% of the G genome segments had recombined with the B genome of wheat with the rest recombining with the D genome of wheat. These results were validated through multi-colour in situ hybridisation analysis. Together these homozygous wheat-T. timopheevii ILs and chromosome-specific KASP markers provide an invaluable resource to wheat breeders for trait discovery to combat biotic and abiotic stress factors affecting wheat production due to climate change.

Wheat speciation and adaptation: perspectives from reticulate evolution

Reticulate evolution through the interchanging of genetic components across organisms can impact significantly on the fitness and adaptation of species. Bread wheat (Triticum aestivum subsp. aestivum) is one of the most important crops in the world. Allopolyploid speciation, frequent hybridization, extensive introgression, and occasional horizontal gene transfer (HGT) have been shaping a typical paradigm of reticulate evolution in bread wheat and its wild relatives, which is likely to have a substantial influence on phenotypic traits and environmental adaptability of bread wheat. In this review, we outlined the evolutionary history of bread wheat and its wild relatives with a highlight on the interspecific hybridization events, demonstrating the reticulate relationship between species/subspecies in the genera Triticum and Aegilops. Furthermore, we discussed the genetic mechanisms and evolutionary significance underlying the introgression of bread wheat and its wild relatives. An in-depth understanding of the evolutionary process of Triticum species should be beneficial to future genetic study and breeding of bread wheat.

Evolutionary rewiring of the wheat transcriptional regulatory network by lineage-specific transposable elements

More than 80% of the wheat genome consists of transposable elements (TEs), which act as major drivers of wheat genome evolution. However, their contributions to the regulatory evolution of wheat adaptations remain largely unclear. Here, we created genome-binding maps for 53 transcription factors (TFs) underlying environmental responses by leveraging DAP-seq in Triticum urartu, together with epigenomic profiles. Most TF binding sites (TFBSs) located distally from genes are embedded in TEs, whose functional relevance is supported by purifying selection and active epigenomic features. About 24% of the non-TE TFBSs share significantly high sequence similarity with TE-embedded TFBSs. These non-TE TFBSs have almost no homologous sequences in non-Triticeae species and are potentially derived from Triticeae-specific TEs. The expansion of TE-derived TFBS linked to wheat-specific gene responses, suggesting TEs are an important driving force for regulatory innovations. Altogether, TEs have been significantly and continuously shaping regulatory networks related to wheat genome evolution and adaptation.

Take-All Disease: New Insights into an Important Wheat Root Pathogen

Take-all disease, caused by the fungal root pathogen Gaeumannomyces tritici, is considered to be the most important root disease of wheat worldwide. Here we review the advances in take-all research over the last 15 years, focusing on the identification of new sources of genetic resistance in wheat relatives and the role of the microbiome in disease development. We also highlight recent breakthroughs in the molecular interactions between G. tritici and wheat, including genome and transcriptome analyses. These new findings will aid the development of novel control strategies against take-all disease. In light of this growing understanding, the G. tritici-wheat interaction could provide a model study system for root-infecting fungal pathogens of cereals.

Polyphenol oxidase genes as integral part of the evolutionary history of domesticated tetraploid wheat

Studying and understanding the genetic basis of polyphenol oxidases (PPO)-related traits plays a crucial role in genetic improvement of crops. A tetraploid wheat collection (T. turgidum ssp., TWC) was analyzed using the 90K wheat SNP iSelect assay and phenotyped for PPO activity. A total of 21,347 polymorphic SNPs were used to perform genome-wide association analysis (GWA) in TWC and durum wheat sub-groups, detecting 23 and 85 marker-trait associations (MTA). In addition, candidate genes responsible for PPO activity were predicted. Based on the 23 MTAs detected in TWC, two haplotypes associated with low and high PPO activity were identified. Four SNPs were developed and validated providing one reliable marker (IWB75732) for marker assisted selection. The 23 MTAs were used to evaluate the genetic divergence (F_ST > 0.25) between the T. turgidum subspecies, providing new information important for understanding the domestication process of Triticum turgidum ssp. and in particular of ssp. carthlicum.

Characterization of the gene family encoding metal tolerance proteins in Triticum urartu: Phylogenetic, transcriptional, and functional analyses

Homeostasis of microelements in organisms is vital for normal metabolism. In plants, the cation diffusion facilitator (CDF) protein family, also known as metal tolerance proteins (MTPs), play critical roles in maintaining trace metal homeostasis. However, little is known about these proteins in wheat. In this study, we characterized the MTP family of Triticum urartu, the donor of 'A' genome of Triticum aestivum, and analysed their phylogenetic relationships, sequence signatures, spatial expression patterns in the diploid wheat, and their transport activity when heterologously expressed. Nine MTPs were identified in the T. urartu genome database, and were classified and designated based on their sequence similarity to Arabidopsis thaliana (Arabidopsis) and Oryza sativa MTPs. Phylogenetic and sequence analyses indicated that the triticum urartu metal tolerance protein (TuMTP)s comprise three Zn-CDFs, two Fe/Zn-CDFs, and four Mn-CDFs; and can be further classified into six subgroups. Among the TuMTPs, there are no MTP2-5 and MTP9-10 counterparts but two MTP1/8/11 orthologs in relation to AtMTPs. It was also shown that members of the same cluster share similar sequence characteristic, i.e. number of introns, predicted transmembrane domains, and motifs. When expressed in yeast, TuMTP1 and TuMTP1.1 conferred tolerance to Zn and Co but not to other metal ions; while TuMTP8, TuMTP8.1, TuMTP11, and TuMTP11.1 conferred tolerance to Mn. When expressed in Arabidopsis, TuMTP1 localized to the tonoplast and significantly enhanced Zn and Co tolerance. TuMTPs showed diverse tissue-specific expression patterns. Taken together, the closely clustered TuMTPs share structural features and metal specificity but play diverse roles in the homeostasis of microelements in plant cells.

The Independent Domestication of Timopheev's Wheat: Insights from Haplotype Analysis of the Brittle rachis 1 (BTR1-A) Gene

Triticum turgidum and T. timopheevii are two tetraploid wheat species sharing T. urartu as a common ancestor, and domesticated accessions from both of these allopolyploids exhibit nonbrittle rachis (i.e., nonshattering spikes). We previously described the loss-of-function mutations in the Brittle Rachis 1 genes BTR1-A and BTR1-B in the A and B subgenomes, respectively, that are responsible for this most visible domestication trait in T. turgidum. Resequencing of a large panel of wild and domesticated T. turgidum accessions subsequently led to the identification of the two progenitor haplotypes of the btr1-A and btr1-B domesticated alleles. Here, we extended the haplotype analysis to other T. turgidum subspecies and to the BTR1 homologues in the related T. timopheevii species. Our results showed that all the domesticated wheat subspecies within T. turgidum share common BTR1-A and BTR1-B haplotypes, confirming their common origin. In T. timopheevii, however, we identified a novel loss-of-function btr1-A allele underlying a partially brittle spike phenotype. This novel recessive allele appeared fixed within the pool of domesticated Timopheev’s wheat but was also carried by one wild timopheevii accession exhibiting partial brittleness. The promoter region for BTR1-B could not be amplified in any T. timopheevii accessions with any T. turgidum primer combination, exemplifying the gene-level distance between the two species. Altogether, our results support the concept of independent domestication processes for the two polyploid, wheat-related species.

Functional gene assessment of bread wheat: breeding implications in Ningxia Province

BACKGROUND

The overall genetic distribution and divergence of cloned genes among bread wheat varieties that have occurred during the breeding process over the past few decades in Ningxia Province, China, are poorly understood. Here, we report the genetic diversities of 44 important genes related to grain yield, quality, adaptation and resistance in 121 Ningxia and 86 introduced wheat cultivars and advanced lines.

RESULTS

The population structure indicated characteristics of genetic components of Ningxia wheat, including landraces of particular genetic resources, introduced varieties with rich genetic diversities and modern cultivars in different periods. Analysis of allele frequencies showed that the dwarfing alleles Rht-B1b at Rht-B1 and Rht-D1b at Rht-D1, 1BL/1RS translocation, Hap-1 at GW2-6B and Hap-H at Sus2-2B are very frequently present in modern Ningxia cultivars and in introduced varieties from other regions but absent in landraces. This indicates that the introduced wheat germplasm with numerous beneficial genes is vital for broadening the genetic diversity of Ningxia wheat varieties. Large population differentiation between modern cultivars and landraces has occurred in adaptation genes. Founder parents carry excellent allele combinations of important genes, with a higher number of favorable alleles than modern cultivars. Gene flow analysis showed that six founder parents have greatly contributed to breeding improvement in Ningxia Province, particularly Zhou 8425B, for yield-related genes.

CONCLUSIONS

Varieties introduced from other regions with rich genetic diversity and landraces with well-adapted genetic resources have been applied to improve modern cultivars. Founder parents, particularly Zhou 8425B, for yield-related genes have contributed greatly to wheat breeding improvement in Ningxia Province. These findings will greatly benefit bread wheat breeding in Ningxia Province as well as other areas with similar ecological environments.

Generation of Doubled Haploid Wheat-Triticum urartu Introgression Lines and Their Characterisation Using Chromosome-Specific KASP Markers

Wheat is one of the most important food and protein sources in the world and although, in recent years wheat breeders have achieved yield gains, they are not sufficient to meet the demands of an ever-growing population. Development of high yielding wheat varieties, resilient to abiotic and biotic stress resulting from climate change, has been limited by wheat's narrow genetic base. In contrast to wheat, the wild relatives of wheat provide a vast reservoir of genetic variation for most, if not all, agronomic traits. Previous studies by the authors have shown the transfer of genetic variation from T. urartu into bread wheat. However, before the introgression lines can be exploited for trait analysis, they are required to have stable transmission of the introgressions to the next generation. In this work, we describe the generation of 86 doubled haploid (DH) wheat-T. urartu introgression lines that carry homozygous introgressions which are stably inherited. The DH lines were characterised using the Axiom^® Wheat Relative Genotyping Array and 151 KASP markers to identify 65 unique T. urartu introgressions in a bread wheat background. DH production has helped accelerate the breeding process and facilitated the early release of homozygous wheat-T. urartu introgression lines. Together with the KASP markers, this valuable resource could greatly advance identification of beneficial alleles that can be used in wheat improvement.

Cereal asparagine synthetase genes

Asparagine synthetase catalyses the transfer of an amino group from glutamine to aspartate to form glutamate and asparagine. The accumulation of free (nonprotein) asparagine in crops has implications for food safety because free asparagine is the precursor for acrylamide, a carcinogenic contaminant that forms during high-temperature cooking and processing. Here we review publicly available genome data for asparagine synthetase genes from species of the Pooideae subfamily, including bread wheat and related wheat species (Triticum and Aegilops spp.), barley (Hordeum vulgare) and rye (Secale cereale) of the Triticeae tribe. Also from the Pooideae subfamily: brachypodium (Brachypodium dIstachyon) of the Brachypodiae tribe. More diverse species are also included, comprising sorghum (Sorghum bicolor) and maize (Zea mays) of the Panicoideae subfamily and rice (Oryza sativa) of the Ehrhartoideae subfamily. The asparagine synthetase gene families of the Triticeae species each comprise five genes per genome, with the genes assigned to four groups: 1, 2, 3 (subdivided into 3.1 and 3.2) and 4. Each species has a single gene per genome in each group, except that some bread wheat varieties (genomes AABBDD) and emmer wheat (Triticum dicoccoides; genomes AABB) lack a group 2 gene in the B genome. This raises questions about the ancestry of cultivated pasta wheat and the B genome donor of bread wheat, suggesting that the hybridisation event that gave rise to hexaploid bread wheat occurred more than once. In phylogenetic analyses, genes from the other species cluster with the Triticeae genes, but brachypodium, sorghum and maize lack a group 2 gene, while rice has only two genes, one group 3 and one group 4. This means that TaASN2, the most highly expressed asparagine synthetase gene in wheat grain, has no equivalent in maize, rice, sorghum or brachypodium. An evolutionary pathway is proposed in which a series of gene duplications gave rise to the five genes found in modern Triticeae species.

Domestication affects the composition, diversity, and co-occurrence of the cereal seed microbiota

Introduction

The seed-associated microbiome has a strong influence on plant ecology, fitness, and productivity. Plant microbiota could be exploited for a more responsible crop management in sustainable agriculture. However, the relationships between seed microbiota and hosts related to the changes from ancestor species to breeded crops still remain poor understood.

Objectives

Our aims were i) to understand the effect of cereal domestication on seed endophytes in terms of diversity, structure and co-occurrence, by comparing four cereal crops and the respective ancestor species; ii) to test the phylogenetic coherence between cereals and their seed microbiota (clue of co-evolution).

Methods

We investigated the seed microbiota of four cereal crops (Triticum aestivum, Triticum monococcum, Triticum durum, and Hordeum vulgare), along with their respective ancestors (Aegilops tauschii, Triticum baeoticum, Triticum dicoccoides, and Hordeum spontaneum, respectively) using 16S rRNA gene metabarcoding, Randomly Amplified Polymorphic DNA (RAPD) profiling of host plants and co-evolution analysis.

Results

The diversity of seed microbiota was generally higher in cultivated cereals than in wild ancestors, suggesting that domestication lead to a bacterial diversification. On the other hand, more microbe-microbe interactions were detected in wild species, indicating a better-structured, mature community. Typical human-associated taxa, such as Cutibacterium, dominated in cultivated cereals, suggesting an interkingdom transfers of microbes from human to plants during domestication. Co-evolution analysis revealed a significant phylogenetic congruence between seed endophytes and host plants, indicating clues of co-evolution between hosts and seed-associated microbes during domestication.

Conclusion

This study demonstrates a diversification of the seed microbiome as a consequence of domestication, and provides clues of co-evolution between cereals and their seed microbiota. This knowledge is useful to develop effective strategies of microbiome exploitation for sustainable agriculture.

Production of synthetic wheat lines to exploit the genetic diversity of emmer wheat and D genome containing Aegilops species in wheat breeding

Due to the accumulation of various useful traits over evolutionary time, emmer wheat (Triticum turgidum subsp. dicoccum and dicoccoides, 2n = 4x = 28; AABB), durum wheat (T. turgidum subsp. durum, 2n = 4x = 28; AABB), T. timopheevii (2n = 4x = 28; AAGG) and D genome containing Aegilops species offer excellent sources of novel variation for the improvement of bread wheat (T. aestivum L., AABBDD). Here, we made 192 different cross combinations between diverse genotypes of wheat and Aegilops species including emmer wheat × Ae. tauschii (2n = DD or DDDD), durum wheat × Ae. tauschii, T. timopheevii × Ae. tauschii, Ae. crassa × durum wheat, Ae. cylindrica × durum wheat and Ae. ventricosa × durum wheat in the field over three successive years. We successfully recovered 56 different synthetic hexaploid and octaploid F₂ lines with AABBDD, AABBDDDD, AAGGDD, D¹D¹X^crX^crAABB, D^cD^cC^cC^cAABB and D^vD^vN^vN^vAABB genomes via in vitro rescue of F₁ embryos and spontaneous production of F₂ seeds on the F_l plants. Cytogenetic analysis of F₂ lines showed that the produced synthetic wheat lines were generally promising stable amphiploids. Contribution of D genome bearing Aegilops and the less-investigated emmer wheat genotypes as parents in the crosses resulted in synthetic amphiploids which are a valuable resource for bread wheat breeding.

Development and deployment of KASP markers for multiple alleles of Lr34 in wheat

Heterogeneous Lr34 genes for leaf rust in winter wheat cultivar 'Duster' and KASP markers for allelic variation in exon 11 and exon 22 of Lr34. Wheat, Triticum aestivum (2n = 6x = 42, AABBDD), is a hexaploid species, and each of three homoeologous genomes A, B, and D should have one copy for a gene in its ancestral form if the gene has no duplication. Previously reported leaf rust resistance gene Lr34 has one copy on the short arm of chromosome 7D in hexaploid wheat, and allelic variation in Lr34 is in intron 4, exon 11, exon 12, or exon 22. In this study, we discovered that Oklahoma hard red winter wheat cultivar 'Duster' (PI 644,016) has two copies of the Lr34 gene, the resistance allele Lr34a and the susceptibility allele Lr34b. Both Lr34a and Lr34b were mapped in the same linkage group on chromosome 7D in a doubled-haploid population generated from a cross between Duster and a winter wheat cultivar 'Billings' which carries the susceptibility allele Lr34c. A chromosomal fragment including Lr34 and at least two neighboring genes on its proximal side but excluding genes on its distal side was duplicated in Duster. The Duster Lr34ab allele was associated with tip necrosis and increased resistance against leaf rust at adult plants in the Duster × Billings DH population tested in the field, demonstrating the function of the Duster Lr34ab allele in wheat. We have developed KASP markers for allelic variation in exon 11 and exon 22 of Lr34 in wheat. These markers can be utilized to accelerate the selection of Lr34 in wheat.

Genome-Wide Discovery of G-Quadruplexes in Wheat: Distribution and Putative Functional Roles

G-quadruplexes are nucleic acid secondary structures formed by a stack of square planar G-quartets. G-quadruplexes were implicated in many biological functions including telomere maintenance, replication, transcription, and translation, in many species including humans and plants. For wheat, however, though it is one of the world's most important staple food, no G-quadruplex studies have been reported to date. Here, we computationally identify putative G4 structures (G4s) in wheat genome for the first time and compare its distribution across the genome against five other genomes (human, maize, Arabidopsis, rice, and sorghum). We identified close to 1 million G4 motifs with a density of 76 G4s/Mb across the whole genome and 93 G4s/Mb over genic regions. Remarkably, G4s were enriched around three regions, two located on the antisense and one on the sense strand at the following positions: 1) the transcription start site (TSS) (antisense), 2) the first coding domain sequence (CDS) (antisense), and 3) the start codon (sense). Functional enrichment analysis revealed that the gene models containing G4 motifs within these peaks were associated with specific gene ontology (GO) terms, such as developmental process, localization, and cellular component organization or biogenesis. We investigated genes encoding MADS-box transcription factors and showed examples of G4 motifs within critical regulatory regions in the VRN-1 genes in wheat. Furthermore, comparison with other plants showed that monocots share a similar distribution of G4s, but Arabidopsis shows a unique G4 distribution. Our study shows for the first time the prevalence and possible functional roles of G4s in wheat.

Identification and characterization of large-scale genomic rearrangements during wheat evolution

Following allopolyploidization, nascent polyploid wheat species react with massive genomic rearrangements, including deletion of transposable element-containing sequences. While such massive rearrangements are considered to be a prominent process in wheat genome evolution and speciation, their structure, extent, and underlying mechanisms remain poorly understood. In this study, we retrieved ~3500 insertions of a specific variant of Fatima, one of the most dynamic gypsy long-terminal repeat retrotransposons in wheat from the recently available high-quality genome drafts of Triticum aestivum (bread wheat) and Triticum turgidum ssp. dicoccoides or wild emmer, the allotetraploid mother of all modern wheats. The dynamic nature of Fatima facilitated the identification of large (i.e., up to ~ 1 million bases) Fatima-containing insertions/deletions (indels) upon comparison of bread wheat and wild emmer genomes. We characterized 11 such indels using computer-assisted analysis followed by PCR validation, and found that they might have occurred via unequal intra-strand recombination or double-strand break (DSB) events. Additionally, we observed one case of introgression of novel DNA fragments from an unknown source into the wheat genome. Our data thus indicate that massive large-scale DNA rearrangements might play a prominent role in wheat speciation.

Partitioning and physical mapping of wheat chromosome 3B and its homoeologue 3E in Thinopyrum elongatum by inducing homoeologous recombination

We performed homoeologous recombination-based partitioning and physical mapping of wheat chromosome 3B and Th. elongatum chromosome 3E, providing a unique physical framework of this homoeologous pair for genome studies. The wheat (Triticum aestivum, 2n = 6x = 42, AABBDD) and Thinopyrum elongatum (2n = 2x = 14, EE) genomes can be differentiated from each other by fluorescent genomic in situ hybridization (FGISH) as well as molecular markers. This has facilitated homoeologous recombination-based partitioning and engineering of their genomes for physical mapping and alien introgression. Here, we constructed a special wheat genotype, which was double monosomic for wheat chromosome 3B and Th. elongatum chromosome 3E and homozygous for the ph1b mutant, to induce 3B-3E homoeologous recombination. Totally, 81 3B-3E recombinants were recovered and detected in the primary, secondary, and tertiary homoeologous recombination cycles by FGISH. Comparing to the primary recombination, the secondary and tertiary recombination shifted toward the proximal regions due to the increase in homology between the pairing partners. The 3B-3E recombinants were genotyped by high-throughput wheat 90-K single nucleotide polymorphism (SNP) arrays and their recombination breakpoints physically mapped based on the FGISH patterns and SNP results. The 3B-3E recombination physically partitioned chromosome 3B into 38 bins, and 429 SNPs were assigned to the distinct bins. Integrative analysis of FGISH and SNP results led to the construction of a composite bin map for chromosome 3B. Additionally, we developed 22 SNP-derived semi-thermal asymmetric reverse PCR markers specific for chromosome 3E and constructed a comparative map of homoeologous chromosomes 3E, 3B, 3A, and 3D. In summary, this work provides a unique physical framework for further studies of the 3B-3E homoeologous pair and diversifies the wheat genome for wheat improvement.

Genome-wide introgression from a bread wheat × Lophopyrum elongatum amphiploid into wheat

We introgressed wheatgrass germplasm from the octoploid amphiploid Triticum aestivum× Lophopyrum elongatum into wheat by manipulating the wheat Ph1 gene and discovered and characterized 130 introgression lines harboring single or, in various combinations, complete and recombined L. elongatum chromosomes. Diploid wheatgrass Lophopyrum elongatum (genomes EE) possesses valuable traits for wheat genetics and breeding. We evaluated several strategies for introgression of this germplasm into wheat. To detect it, we developed and validated multiplexed sets of Sequenom MassARRAY single nucleotide polymorphism (SNP) markers, which differentiated disomic and monosomic L. elongatum chromosomes from wheat chromosomes. We identified 130 introgression lines (ILs), which harbored 108 complete and 89 recombined L. elongatum chromosomes. Of the latter, 59 chromosomes were recombined by one or more crossovers and 30 were involved in centromeric (Robertsonian) translocations or were telocentric. To identify wheat chromosomes substituted for or recombined with L. elongatum chromosomes, we genotyped the ILs with the wheat 90-K Infinium SNP array. We found that most of the wheat 90-K probes correctly detected their targets in the L. elongatum genome and showed that some wheat SNPs are ancient and had originated prior to the divergence of the wheat and L. elongatum lineages. Of the 130 ILs, 52% were homozygous for Ph1 deletion and thus are staged to be recombined further. We failed to detect in the L. elongatum genome the 4/5 reciprocal translocation that has been reported in Thinopyrum bessarabicum and several other Triticeae genomes.

An ankyrin-repeat and WRKY-domain-containing immune receptor confers stripe rust resistance in wheat

Perception of pathogenic effectors in plants often relies on nucleotide-binding domain (NBS) and leucine-rich-repeat-containing (NLR) proteins. Some NLRs contain additional domains that function as integrated decoys for pathogen effector targets and activation of immune signalling. Wheat stripe rust is one of the most devastating diseases of crop plants. Here, we report the cloning of YrU1, a stripe rust resistance gene from the diploid wheat Triticum urartu, the progenitor of the A genome of hexaploid wheat. YrU1 encodes a coiled-coil-NBS-leucine-rich repeat protein with N-terminal ankyrin-repeat and C-terminal WRKY domains, representing a unique NLR structure in plants. Database searches identify similar architecture only in wheat relatives. Transient expression of YrU1 in Nicotiana benthamiana does not induce cell death in the absence of pathogens. The ankyrin-repeat and coiled-coil domains of YrU1 self-associate, suggesting that homodimerisation is critical for YrU1 function. The identification and cloning of this disease resistance gene sheds light on NLR protein function and may facilitate breeding to control the devastating wheat stripe rust disease.

Genome-wide analysis of the HSP101/CLPB gene family for heat tolerance in hexaploid wheat

Heat Shock Protein 101 (HSP101), the homolog of Caseinolytic Protease B (CLPB) proteins, has functional conservation across species to play roles in heat acclimation and plant development. In wheat, several TaHSP101/CLPB genes were identified, but have not been comprehensively characterized. Given the complexity of a polyploid genome with its phenomena of homoeologous expression bias, detailed analysis on the whole TaCLPB family members is important to understand the genetic basis of heat tolerance in hexaploid wheat. In this study, a genome-wide analysis revealed thirteen members of TaCLPB gene family and their expression patterns in various tissues, developmental stages, and stress conditions. Detailed characterization of TaCLPB gene and protein structures suggested potential variations of the sub-cellular localization and their functional regulations. We revealed homoeologous specific variations among TaCLPB gene copies that have not been reported earlier. A study of the Chromosome 1 TaCLPB in four wheat genotypes demonstrated unique patterns of the homoeologous gene expression under moderate and extreme heat treatments. The results give insight into the strategies to improve heat tolerance by targeting one or some of the TaCLPB genes in wheat.

Identification of New Leaf Rust Resistance Loci in Wheat and Wild Relatives by Array-Based SNP Genotyping and Association Genetics

Leaf rust caused by Puccinia triticina is the most widespread rust disease of wheat. As pathogen populations are constantly evolving, identification of novel sources of resistance is necessary to maintain disease resistance and stay ahead of this plant-pathogen evolutionary arms race. The wild genepool of wheat is a rich source of genetic diversity, accounting for 44% of the Lr genes identified. Here we performed a genome-wide association study (GWAS) on a diverse germplasm of 385 accessions, including 27 different Triticum and Aegilops species. Genetic characterization using the wheat 90 K array and subsequent filtering identified a set of 20,501 single nucleotide polymorphic (SNP) markers. Of those, 9,570 were validated using exome capture and mapped onto the Chinese Spring reference sequence v1.0. Phylogenetic analyses illustrated four major clades, clearly separating the wild species from the T. aestivum and T. turgidum species. GWAS was conducted using eight statistical models for infection types against six leaf rust isolates and leaf rust severity rated in field trials for 3-4 years at 2-3 locations in Canada. Functional annotation of genes containing significant quantitative trait nucleotides (QTNs) identified 96 disease-related loci associated with leaf rust resistance. A total of 21 QTNs were in haplotype blocks or within flanking markers of at least 16 known Lr genes. The remaining significant QTNs were considered loci that putatively harbor new Lr resistance genes. Isolation of these candidate genes will contribute to the elucidation of their role in leaf rust resistance and promote their usefulness in marker-assisted selection and introgression.

Recombination between homoeologous chromosomes induced in durum wheat by the Aegilops speltoides Su1-Ph1 suppressor

Su1-Ph1, which we previously introgressed into wheat from Aegilops speltoides, is a potent suppressor of Ph1 and a valuable tool for gene introgression in tetraploid wheat. We previously introgressed Su1-Ph1, a suppressor of the wheat Ph1 gene, from Aegilops speltoides into durum wheat cv Langdon (LDN). Here, we evaluated the utility of the introgressed suppressor for inducing introgression of alien germplasm into durum wheat. We built LDN plants heterozygous for Su1-Ph1 that simultaneously contained a single LDN chromosome 5B and a single Ae. searsii chromosome 5S^se, which targeted them for recombination. We genotyped 28 BC₁F₁ and 84 F₂ progeny with the wheat 90-K Illumina single-nucleotide polymorphism assay and detected extensive recombination between the two chromosomes, which we confirmed by non-denaturing fluorescence in situ hybridization (ND-FISH). We constructed BC₁F₁ and F₂ genetic maps that were 65.31 and 63.71 cM long, respectively. Recombination rates between the 5B and 5S^se chromosomes were double the expected rate computed from their meiotic pairing, which we attributed to selection against aneuploid gametes. Recombination rate between 5B and 5S^se was depressed compared to that between 5B chromosomes in the proximal region of the long arm. We integrated ND-FISH signals into the genetic map and constructed a physical map, which we used to map a 172,188,453-bp Ph1 region. Despite the location of the region in a low-recombination region of the 5B chromosome, we detected three crossovers in it. Our data show that Su1-Ph1 is a valuable tool for gene introgression and gene mapping based on recombination between homoeologous chromosomes in wheat.

Genome-Wide Characterization and Expression Profiling of Squamosa Promoter Binding Protein-like (SBP) Transcription Factors in Wheat (Triticum aestivum L.)

Transcription factors (TFs) play fundamental roles in the developmental processes of all living organisms. Squamosa Promoter Binding Protein-like (SBP/SBP-Box) is a major family of plant-specific TFs, which plays important roles in multiple processes involving plant growth and development. While some work has been done, there is a lot more that is yet to be discovered in the hexaploid wheat SBP (TaSBP) family. With the completion of whole genome sequencing, genome-wide analysis of SBPs in common hexaploid wheat is now possible. In this study, we used protein–protein Basic Local Alignment Search Tool (BLASTp) to hunt the newly released reference genome sequence of hexaploid wheat (Chinese spring). Seventy-four TaSBP proteins (belonging to 56 genes) were identified and clustered into five groups. Gene structure and motif analysis indicated that most TaSBPs have relatively conserved exon–intron arrangements and motif composition. Analysis of transcriptional data showed that many TaSBP genes responded to some biological and abiotic stresses with different expression patterns. Moreover, three TaSBP genes were generally expressed in the majority of tissues throughout the wheat growth and also responded to many environmental biotic and abiotic stresses. Collectively, the detailed analyses presented here will help in understanding the roles of the TaSBP and also provide a reference for the further study of its biological function in wheat.

The novel function of the Ph1 gene to differentiate homologs from homoeologs evolved in Triticum turgidum ssp. dicoccoides via a dramatic meiosis-specific increase in the expression of the 5B copy of the C-Ph1 gene

The Ph1 gene is the principal regulator of homoeologous chromosome pairing control (HECP) that ensures the diploid-like meiotic chromosome pairing behavior of polyploid wheat. The HECP control was speculated to have evolved after the first event of polyploidization. With the objective to accurately understand the evolution of the HECP control, wild emmer wheat accessions previously known to differ for HECP control were characterized for the structure and expression of the candidate Ph1 gene, C-Ph1. The C-TdPh1-5A and 5B gene copies of emmer wheat showed 98 and 99% DNA sequence similarity respectively with the corresponding hexaploid wheat copies. Further, the C-TdPh1-5B carried the C-Ph1-5B specific structural changes and transcribed three splice variants as observed in the hexaploid wheat. Further, single nucleotide changes differentiating accessions varying for HECP control were identified. Analyzed by quantitative expression analysis, the wild emmer accessions with HECP control showed ~ 10,000-fold higher transcript abundance of the C-TdPh1-5B copy during prophase-I compared to accessions lacking the control. Differential transcriptional regulation of C-TdPh1-5B splice variants further revealed that C-Ph1-5B^alt1 variant is mainly responsible for differential accumulation of C-Ph1-5B copy in accessions with HECP control. Taken together, these results showed that the HECP control evolved via transcriptional regulation of splice variants during meiosis.

Genome-wide variation patterns between landraces and cultivars uncover divergent selection during modern wheat breeding

Genetic diversity, population structure, LD decay, and selective sweeps in 687 wheat accessions were analyzed, providing relevant guidelines to facilitate the use of the germplasm in wheat breeding. Common wheat (Triticum aestivum L.) is one of the most widely grown crops in the world. Landraces were subjected to strong human-mediated selection in developing high-yielding, good quality, and widely adapted cultivars. To investigate the genome-wide patterns of allelic variation, population structure and patterns of selective sweeps during modern wheat breeding, we tested 687 wheat accessions, including landraces (148) and cultivars (539) mainly from China and Pakistan in a wheat 90 K single nucleotide polymorphism array. Population structure analysis revealed that cultivars and landraces from China and Pakistan comprised three relatively independent genetic clusters. Cultivars displayed lower nucleotide diversity and a wider average LD decay across whole genome, indicating allelic erosion and a diversity bottleneck due to the modern breeding. Analysis of genetic differentiation between landraces and cultivars from China and Pakistan identified allelic variants subjected to selection during modern breeding. In total, 477 unique genome regions showed signatures of selection, where 109 were identified in both China and Pakistan germplasm. The majority of genomic regions were located in the B genome (225), followed by the A genome (175), and only 77 regions were located in the D genome. EigenGWAS was further used to identify key selection loci in modern wheat cultivars from China and Pakistan by comparing with global winter wheat and spring wheat diversity panels, respectively. A few known functional genes or loci found within these genome regions corresponded to known phenotypes for disease resistance, vernalization, quality, adaptability and yield-related traits. This study uncovered molecular footprints of modern wheat breeding and explained the genetic basis of polygenic adaptation in wheat. The results will be useful for understanding targets of modern wheat breeding, and in devising future breeding strategies to target beneficial alleles currently not pursued.