Multiple wheat genomes reveal global variation in modern breeding.

Advances in genomics have expedited the improvement of several agriculturally important crops but similar efforts in wheat (Triticum spp.) have been more challenging. This is largely owing to the size and complexity of the wheat genome1, and the lack of genome-assembly data for multiple wheat lines2,3. Here we generated ten chromosome pseudomolecule and five scaffold assemblies of hexaploid wheat to explore the genomic diversity among wheat lines from global breeding programs. Comparative analysis revealed extensive structural rearrangements, introgressions from wild relatives and differences in gene content resulting from complex breeding histories aimed at improving adaptation to diverse environments, grain yield and quality, and resistance to stresses4,5. We provide examples outlining the utility of these genomes, including a detailed multi-genome-derived nucleotide-binding leucine-rich repeat protein repertoire involved in disease resistance and the characterization of Sm16, a gene associated with insect resistance. These genome assemblies will provide a basis for functional gene discovery and breeding to deliver the next generation of modern wheat cultivars.

Copy number variation of TdDof controls solid-stemmed architecture in wheat.

Stem solidness is an important agronomic trait of durum (Triticum turgidum L. var. durum) and bread (Triticum aestivum L.) wheat that provides resistance to the wheat stem sawfly. This dominant trait is conferred by the SSt1 locus on chromosome 3B. However, the molecular identity and mechanisms underpinning stem solidness have not been identified. Here, we demonstrate that copy number variation of TdDof, a gene encoding a putative DNA binding with one finger protein, controls the stem solidness trait in wheat. Using map-based cloning, we localized TdDof to within a physical interval of 2.1 Mb inside the SSt1 locus. Molecular analysis revealed that hollow-stemmed wheat cultivars such as Kronos carry a single copy of TdDof, whereas solid-stemmed cultivars such as CDC Fortitude carry multiple identical copies of the gene. Deletion of all TdDof copies from CDC Fortitude resulted in the loss of stem solidness, whereas the transgenic overexpression of TdDof restored stem solidness in the TdDof deletion mutant pithless1 and conferred stem solidness in Kronos. In solid-stemmed cultivars, increased TdDof expression was correlated with the down-regulation of genes whose orthologs have been implicated in programmed cell death (PCD) in other species. Anatomical and histochemical analyses revealed that hollow-stemmed lines had stronger PCD-associated signals in the pith cells compared to solid-stemmed lines, which suggests copy number-dependent expression of TdDof could be directly or indirectly involved in the negative regulation of PCD. These findings provide opportunities to manipulate stem development in wheat and other monocots for agricultural or industrial purposes.

Dissection of a rapidly evolving wheat resistance gene cluster by long-read genome sequencing accelerated the cloning of Pm69.

Gene cloning in repeat-rich polyploid genomes remains challenging. Here, we describe a strategy for overcoming major bottlenecks in cloning of the powdery mildew resistance gene (R-gene) Pm69 derived from tetraploid wild emmer wheat. A conventional positional cloning approach was not effective owing to suppressed recombination. Chromosome sorting was compromised by insufficient purity. A Pm69 physical map, constructed by assembling Oxford Nanopore Technology (ONT) long-read genome sequences, revealed a rapidly evolving nucleotide-binding leucine-rich repeat (NLR) R-gene cluster with structural variations. A single candidate NLR was identified by anchoring RNA sequencing reads from susceptible mutants to ONT contigs and was validated by virus-induced gene silencing. Pm69 is likely a newly evolved NLR and was discovered in only one location across the wild emmer wheat distribution range in Israel. Pm69 was successfully introgressed into cultivated wheat, and a diagnostic molecular marker was used to accelerate its deployment and pyramiding with other R-genes.

Genetic architecture of rust resistance in a wheat (Triticum turgidum) diversity panel.

Introduction: Wheat rust diseases are widespread and affect all wheat growing areas around the globe. Breeding strategies focus on incorporating genetic disease resistance. However, pathogens can quickly evolve and overcome the resistance genes deployed in commercial cultivars, creating a constant need for identifying new sources of resistance. Methods: We have assembled a diverse tetraploid wheat panel comprised of 447 accessions of three Triticum turgidum subspecies and performed a genome-wide association study (GWAS) for resistance to wheat stem, stripe, and leaf rusts. The panel was genotyped with the 90K Wheat iSelect single nucleotide polymorphism (SNP) array and subsequent filtering resulted in a set of 6,410 non-redundant SNP markers with known physical positions. Results: Population structure and phylogenetic analyses revealed that the diversity panel could be divided into three subpopulations based on phylogenetic/geographic relatedness. Marker-trait associations (MTAs) were detected for two stem rust, two stripe rust and one leaf rust resistance loci. Of them, three MTAs coincide with the known rust resistance genes Sr13, Yr15 and Yr67, while the other two may harbor undescribed resistance genes. Discussion: The tetraploid wheat diversity panel, developed and characterized herein, captures wide geographic origins, genetic diversity, and evolutionary history since domestication making it a useful community resource for mapping of other agronomically important traits and for conducting evolutionary studies.

Multi-locus genome-wide association studies reveal the genetic architecture of Fusarium head blight resistance in durum wheat.

Durum wheat is more susceptible to Fusarium head blight (FHB) than other types or classes of wheat. The disease is one of the most devastating in wheat; it reduces yield and end-use quality and contaminates the grain with fungal mycotoxins such as deoxynivalenol (DON). A panel of 265 Canadian and European durum wheat cultivars, as well as breeding and experimental lines, were tested in artificially inoculated field environments (2019-2022, inclusive) and two greenhouse trials (2019 and 2020). The trials were assessed for FHB severity and incidence, visual rating index, Fusarium-damaged kernels, DON accumulation, anthesis or heading date, maturity date, and plant height. In addition, yellow pigment and protein content were analyzed for the 2020 field season. To capture loci underlying FHB resistance and related traits, GWAS was performed using single-locus and several multi-locus models, employing 13,504 SNPs. Thirty-one QTL significantly associated with one or more FHB-related traits were identified, of which nine were consistent across environments and associated with multiple FHB-related traits. Although many of the QTL were identified in regions previously reported to affect FHB, the QTL QFhb-3B.2, associated with FHB severity, incidence, and DON accumulation, appears to be novel. We developed KASP markers for six FHB-associated QTL that were consistently detected across multiple environments and validated them on the Global Durum Panel (GDP). Analysis of allelic diversity and the frequencies of these revealed that the lines in the GDP harbor between zero and six resistance alleles. This study provides a comprehensive assessment of the genetic basis of FHB resistance and DON accumulation in durum wheat. Accessions with multiple favorable alleles were identified and will be useful genetic resources to improve FHB resistance in durum breeding programs through marker-assisted recurrent selection and gene stacking.

Discovery of stripe rust resistance with incomplete dominance in wild emmer wheat using bulked segregant analysis sequencing.

Durable crop disease resistance is an essential component of global food security. Continuous pathogen evolution leads to a breakdown of resistance and there is a pressing need to characterize new resistance genes for use in plant breeding. Here we identified an accession of wild emmer wheat (Triticum turgidum ssp. dicoccoides), PI 487260, that is highly resistant to multiple stripe rust isolates. Genetic analysis revealed resistance was conferred by a single, incompletely dominant gene designated as Yr84. Through bulked segregant analysis sequencing (BSA-Seq) we identified a 52.7 Mb resistance-associated interval on chromosome 1BS. Detected variants were used to design genetic markers for recombinant screening, further refining the interval of Yr84 to a 2.3-3.3 Mb in tetraploid wheat genomes. This interval contains 34 candidate genes encoding for protein domains involved in disease resistance responses. Furthermore, KASP markers closely-linked to Yr84 were developed to facilitate marker-assisted selection for rust resistance breeding.

The effects of crop attributes, selection, and recombination on Canadian bread wheat molecular variation.

Cultivated germplasm provides an opportunity to investigate how crop agronomic traits, selection for major genes, and differences in crossing-over rates drive patterns of allelic variation. To identify how these factors correlated with allelic variation within a collection of cultivated bread wheat (Triticum aestivum L.), we generated genotypes for 388 accessions grown in Canada over the past 170 yr using filtered single nucleotide polymorphism (SNP) calls from an Illumina Wheat iSelect 90K SNP-array. Entries' breeding program, era of release, grain texture, kernel color, and growth habit contributed to allelic differentiation. Allelic diversity and linkage disequilibrium (LD) of markers flanking some major loci known to affect traits such as gluten strength, growth habit, and grain color were consistent with selective sweeps. Nonetheless, some flanking markers of major loci had low LD and high allelic diversity. Positive selection may have acted upon homoeologous genes that had significant enrichment for the gene ontology terms 'response-to-auxin' and 'response-to-wounding.' Long regions of LD, spanning approximately one-third the length of entire chromosomes, were associated with many pericentromeric regions. These regions were also characterized by low diversity. Enhancing recombination across these regions could generate novel allele combinations to accelerate Canadian wheat improvement.

Exome sequencing highlights the role of wild-relative introgression in shaping the adaptive landscape of the wheat genome.

Introgression is a potential source of beneficial genetic diversity. The contribution of introgression to adaptive evolution and improvement of wheat as it was disseminated worldwide remains unknown. We used targeted re-sequencing of 890 diverse accessions of hexaploid and tetraploid wheat to identify wild-relative introgression. Introgression, and selection for improvement and environmental adaptation, each reduced deleterious allele burden. Introgression increased diversity genome wide and in regions harboring major agronomic genes, and contributed alleles explaining a substantial proportion of phenotypic variation. These results suggest that historic gene flow from wild relatives made a substantial contribution to the adaptive diversity of modern bread wheat.

Publisher Correction: Exome sequencing highlights the role of wild-relative introgression in shaping the adaptive landscape of the wheat genome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Durum wheat genome highlights past domestication signatures and future improvement targets.

The domestication of wild emmer wheat led to the selection of modern durum wheat, grown mainly for pasta production. We describe the 10.45 gigabase (Gb) assembly of the genome of durum wheat cultivar Svevo. The assembly enabled genome-wide genetic diversity analyses revealing the changes imposed by thousands of years of empirical selection and breeding. Regions exhibiting strong signatures of genetic divergence associated with domestication and breeding were widespread in the genome with several major diversity losses in the pericentromeric regions. A locus on chromosome 5B carries a gene encoding a metal transporter (TdHMA3-B1) with a non-functional variant causing high accumulation of cadmium in grain. The high-cadmium allele, widespread among durum cultivars but undetected in wild emmer accessions, increased in frequency from domesticated emmer to modern durum wheat. The rapid cloning of TdHMA3-B1 rescues a wild beneficial allele and demonstrates the practical use of the Svevo genome for wheat improvement.

High density genetic mapping of Fusarium head blight resistance QTL in tetraploid wheat.

Breeding for Fusarium head blight (FHB) resistance in durum wheat is complicated by the quantitative trait expression and narrow genetic diversity of available resources. High-density mapping of the FHB resistance quantitative trait loci (QTL), evaluation of their co-localization with plant height and maturity QTL and the interaction among the identified QTL are the objectives of this study. Two doubled haploid (DH) populations, one developed from crosses between Triticum turgidum ssp. durum lines DT707 and DT696 and the other between T. turgidum ssp. durum cv. Strongfield and T. turgidum ssp. carthlicum cv. Blackbird were genotyped using the 90K Infinium iSelect chip and evaluated phenotypically at multiple field FHB nurseries over years. A moderate broad-sense heritability indicated a genotype-by-environment interaction for the expression of FHB resistance in both populations. Resistance QTL were identified for the DT707 × DT696 population on chromosomes 1B, 2B, 5A (two loci) and 7A and for the Strongfield × Blackbird population on chromosomes 1A, 2A, 2B, 3A, 6A, 6B and 7B with the QTL on chromosome 1A and those on chromosome 5A being more consistently expressed over environments. FHB resistance co-located with plant height and maturity QTL on chromosome 5A and with a maturity QTL on chromosome 7A for the DT707 × DT696 population. Resistance also co-located with plant height QTL on chromosomes 2A and 3A and with maturity QTL on chromosomes 1A and 7B for the Strongfield × Blackbird population. Additive × additive interactions were identified, for example between the two FHB resistance QTL on chromosome 5A for the DT707 × DT696 population and the FHB resistance QTL on chromosomes 1A and 7B for the Strongfield × Blackbird population. Application of the Single Nucleotide Polymorphic (SNP) markers associated with FHB resistance QTL identified in this study will accelerate combining genes from the two populations.

High density mapping and haplotype analysis of the major stem-solidness locus SSt1 in durum and common wheat.

Breeding for solid-stemmed durum (Triticum turgidum L. var durum) and common wheat (Triticum aestivum L.) cultivars is one strategy to minimize yield losses caused by the wheat stem sawfly (Cephus cinctus Norton). Major stem-solidness QTL have been localized to the long arm of chromosome 3B in both wheat species, but it is unclear if these QTL span a common genetic interval. In this study, we have improved the resolution of the QTL on chromosome 3B in a durum (Kofa/W9262-260D3) and common wheat (Lillian/Vesper) mapping population. Coincident QTL (LOD = 94-127, R2 = 78-92%) were localized near the telomere of chromosome 3BL in both mapping populations, which we designate SSt1. We further examined the SSt1 interval by using available consensus maps for durum and common wheat and compared genetic to physical intervals by anchoring markers to the current version of the wild emmer wheat (WEW) reference sequence. These results suggest that the SSt1 interval spans a physical distance of 1.6 Mb in WEW (positions 833.4-835.0 Mb). In addition, minor QTL were identified on chromosomes 2A, 2D, 4A, and 5A that were found to synergistically enhance expression of SSt1 to increase stem-solidness. These results suggest that developing new wheat cultivars with improved stem-solidness is possible by combining SSt1 with favorable alleles at minor loci within both wheat species.

Wild emmer genome architecture and diversity elucidate wheat evolution and domestication.

Wheat (Triticum spp.) is one of the founder crops that likely drove the Neolithic transition to sedentary agrarian societies in the Fertile Crescent more than 10,000 years ago. Identifying genetic modifications underlying wheat's domestication requires knowledge about the genome of its allo-tetraploid progenitor, wild emmer (T. turgidum ssp. dicoccoides). We report a 10.1-gigabase assembly of the 14 chromosomes of wild tetraploid wheat, as well as analyses of gene content, genome architecture, and genetic diversity. With this fully assembled polyploid wheat genome, we identified the causal mutations in Brittle Rachis 1 (TtBtr1) genes controlling shattering, a key domestication trait. A study of genomic diversity among wild and domesticated accessions revealed genomic regions bearing the signature of selection under domestication. This reference assembly will serve as a resource for accelerating the genome-assisted improvement of modern wheat varieties.