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.

De Novo Genome Assembly of the Japanese Wheat Cultivar Norin 61 Highlights Functional Variation in Flowering Time and Fusarium-Resistant Genes in East Asian Genotypes.

Bread wheat is a major crop that has long been the focus of basic and breeding research. Assembly of its genome has been difficult because of its large size and allohexaploid nature (AABBDD genome). Following the first reported assembly of the genome of the experimental strain Chinese Spring (CS), the 10+ Wheat Genomes Project was launched to produce multiple assemblies of worldwide modern cultivars. The only Asian cultivar in the project is Norin 61, a representative Japanese cultivar adapted to grow across a broad latitudinal range, mostly characterized by a wet climate and a short growing season. Here, we characterize the key aspects of its chromosome-scale genome assembly spanning 15 Gb with a raw scaffold N50 of 22 Mb. Analysis of the repetitive elements identified chromosomal regions unique to Norin 61 that encompass a tandem array of the pathogenesis-related 13 family. We report novel copy-number variations in the B homeolog of the florigen gene FT1/VRN3, pseudogenization of its D homeolog and the association of its A homeologous alleles with the spring/winter growth habit. Furthermore, the Norin 61 genome carries typical East Asian functional variants different from CS, ranging from a single nucleotide to multi-Mb scale. Examples of such variation are the Fhb1 locus, which confers Fusarium head-blight resistance, Ppd-D1a, which confers early flowering, Glu-D1f for Asian noodle quality and Rht-D1b, which introduced semi-dwarfism during the green revolution. The adoption of Norin 61 as a reference assembly for functional and evolutionary studies will enable comprehensive characterization of the underexploited Asian bread wheat diversity.

Resistance and susceptibility QTL identified in a rice MAGIC population by screening with a minor-effect virulence factor from Xanthomonas oryzae pv. oryzae.

Effective and durable disease resistance for bacterial blight (BB) of rice is a continuous challenge due to the evolution and adaptation of the pathogen, Xanthomonas oryzae pv. oryzae (Xoo), on cultivated rice varieties. Fundamental to this pathogens' virulence is transcription activator-like (TAL) effectors that activate transcription of host genes and contribute differently to pathogen virulence, fitness or both. Host plant resistance is predicted to be more durable if directed at strategic virulence factors that impact both pathogen virulence and fitness. We characterized Tal7b, a minor-effect virulence factor that contributes incrementally to pathogen virulence in rice, is a fitness factor to the pathogen and is widely present in geographically diverse strains of Xoo. To identify sources of resistance to this conserved effector, we used a highly virulent strain carrying a plasmid borne copy of Tal7b to screen an indica multi-parent advanced generation inter-cross (MAGIC) population. Of 18 QTL revealed by genome-wide association studies and interval mapping analysis, six were specific to Tal7b (qBB-tal7b). Overall, 150 predicted Tal7b gene targets overlapped with qBB-tal7b QTL. Of these, 21 showed polymorphisms in the predicted effector binding element (EBE) site and 23 lost the EBE sequence altogether. Inoculation and bioinformatics studies suggest that the Tal7b target in one of the Tal7b-specific QTL, qBB-tal7b-8, is a disease susceptibility gene and that the resistance mechanism for this locus may be through loss of susceptibility. Our work demonstrates that minor-effect virulence factors significantly contribute to disease and provide a potential new approach to identify effective disease resistance.

Allelic variation for broad-spectrum resistance and susceptibility to bacterial pathogens identified in a rice MAGIC population.

Quantitative trait loci (QTL) that confer broad-spectrum resistance (BSR), or resistance that is effective against multiple and diverse plant pathogens, have been elusive targets of crop breeding programmes. Multiparent advanced generation intercross (MAGIC) populations, with their diverse genetic composition and high levels of recombination, are potential resources for the identification of QTL for BSR. In this study, a rice MAGIC population was used to map QTL conferring BSR to two major rice diseases, bacterial leaf streak (BLS) and bacterial blight (BB), caused by Xanthomonas oryzae pathovars (pv.) oryzicola (Xoc) and oryzae (Xoo), respectively. Controlling these diseases is particularly important in sub-Saharan Africa, where no sources of BSR are currently available in deployed varieties. The MAGIC founders and lines were genotyped by sequencing and phenotyped in the greenhouse and field by inoculation with multiple strains of Xoc and Xoo. A combination of genomewide association studies (GWAS) and interval mapping analyses revealed 11 BSR QTL, effective against both diseases, and three pathovar-specific QTL. The most promising BSR QTL (qXO-2-1, qXO-4-1 and qXO-11-2) conferred resistance to more than nine Xoc and Xoo strains. GWAS detected 369 significant SNP markers with distinguishable phenotypic effects, allowing the identification of alleles conferring disease resistance and susceptibility. The BSR and susceptibility QTL will improve our understanding of the mechanisms of both resistance and susceptibility in the long term and will be immediately useful resources for rice breeding programmes.

Representing true plant genomes: haplotype-resolved hybrid pepper genome with trio-binning.

As sequencing costs decrease and availability of high fidelity long-read sequencing increases, generating experiment specific de novo genome assemblies becomes feasible. In many crop species, obtaining the genome of a hybrid or heterozygous individual is necessary for systems that do not tolerate inbreeding or for investigating important biological questions, such as hybrid vigor. However, most genome assembly methods that have been used in plants result in a merged single sequence representation that is not a true biologically accurate representation of either haplotype within a diploid individual. The resulting genome assembly is often fragmented and exhibits a mosaic of the two haplotypes, referred to as haplotype-switching. Important haplotype level information, such as causal mutations and structural variation is therefore lost causing difficulties in interpreting downstream analyses. To overcome this challenge, we have applied a method developed for animal genome assembly called trio-binning to an intra-specific hybrid of chili pepper (Capsicum annuum L. cv. HDA149 x Capsicum annuum L. cv. HDA330). We tested all currently available softwares for performing trio-binning, combined with multiple scaffolding technologies including Bionano to determine the optimal method of producing the best haplotype-resolved assembly. Ultimately, we produced highly contiguous biologically true haplotype-resolved genome assemblies for each parent, with scaffold N50s of 266.0 Mb and 281.3 Mb, with 99.6% and 99.8% positioned into chromosomes respectively. The assemblies captured 3.10 Gb and 3.12 Gb of the estimated 3.5 Gb chili pepper genome size. These assemblies represent the complete genome structure of the intraspecific hybrid, as well as the two parental genomes, and show measurable improvements over the currently available reference genomes. Our manuscript provides a valuable guide on how to apply trio-binning to other plant genomes.

High molecular weight glutenin gene diversity in Aegilops tauschii demonstrates unique origin of superior wheat quality.

Central to the diversity of wheat products was the origin of hexaploid bread wheat, which added the D-genome of Aegilops tauschii to tetraploid wheat giving rise to superior dough properties in leavened breads. The polyploidization, however, imposed a genetic bottleneck, with only limited diversity introduced in the wheat D-subgenome. To understand genetic variants for quality, we sequenced 273 accessions spanning the known diversity of Ae. tauschii. We discovered 45 haplotypes in Glu-D1, a major determinant of quality, relative to the two predominant haplotypes in wheat. The wheat allele 2 + 12 was found in Ae. tauschii Lineage 2, the donor of the wheat D-subgenome. Conversely, the superior quality wheat allele 5 + 10 allele originated in Lineage 3, a recently characterized lineage of Ae. tauschii, showing a unique origin of this important allele. These two wheat alleles were also quite similar relative to the total observed molecular diversity in Ae. tauschii at Glu-D1. Ae. tauschii is thus a reservoir for unique Glu-D1 alleles and provides the genomic resource to begin utilizing new alleles for end-use quality improvement in wheat breeding programs.

Genome-based prediction of multiple wheat quality traits in multiple years.

Wheat quality improvement is an important objective in all wheat breeding programs. However, due to the cost, time and quantity of seed required, wheat quality is typically analyzed only in the last stages of the breeding cycle on a limited number of samples. The use of genomic prediction could greatly help to select for wheat quality more efficiently by reducing the cost and time required for this analysis. Here were evaluated the prediction performances of 13 wheat quality traits under two multi-trait models (Bayesian multi-trait multi-environment [BMTME] and multi-trait ridge regression [MTR]) using five data sets of wheat lines evaluated in the field during two consecutive years. Lines in the second year (testing) were predicted using the quality information obtained in the first year (training). For most quality traits were found moderate to high prediction accuracies, suggesting that the use of genomic selection could be feasible. The best predictions were obtained with the BMTME model in all traits and the worst with the MTR model. The best predictions with the BMTME model under the mean arctangent absolute percentage error (MAAPE) were for test weight across the five data sets, whereas the worst predictions were for the alveograph trait ALVPL. In contrast, under Pearson's correlation, the best predictions depended on the data set. The results obtained suggest that the BMTME model should be preferred for multi-trait prediction analyses. This model allows to obtain not only the correlation among traits, but also the correlation among environments, helping to increase the prediction accuracy.