A haplotype-led approach to increase the precision of wheat breeding

Omics-driven utilization of wild relatives for empowering pre-breeding in pearl millet

MAIN CONCLUSION

Pearl millet wild relatives harbour novel alleles which could be utilized to broaden genetic base of cultivated species. Genomics-informed pre-breeding is needed to speed up introgression from wild to cultivated gene pool in pearl millet. Rising episodes of intense biotic and abiotic stresses challenge pearl millet production globally. Wild relatives provide a wide spectrum of novel alleles which could address challenges posed by climate change. Pre-breeding holds potential to introgress novel diversity in genetically narrow cultivated Pennisetum glaucum from diverse gene pool. Practical utilization of gene pool diversity remained elusive due to genetic intricacies. Harnessing promising traits from wild pennisetum is limited by lack of information on underlying candidate genes/QTLs. Next-Generation Omics provide vast scope to speed up pre-breeding in pearl millet. Genomic resources generated out of draft genome sequence and improved genome assemblies can be employed to utilize gene bank accessions effectively. The article highlights genetic richness in pearl millet and its utilization with a focus on harnessing next-generation Omics to empower pre-breeding.

Climate Change-The Rise of Climate-Resilient Crops

Climate change disrupts food production in many regions of the world. The accompanying extreme weather events, such as droughts, floods, heat waves, and cold snaps, pose threats to crops. The concentration of carbon dioxide also increases in the atmosphere. The United Nations is implementing the climate-smart agriculture initiative to ensure food security. An element of this project involves the breeding of climate-resilient crops or plant cultivars with enhanced resistance to unfavorable environmental conditions. Modern agriculture, which is currently homogeneous, needs to diversify the species and cultivars of cultivated plants. Plant breeding programs should extensively incorporate new molecular technologies, supported by the development of field phenotyping techniques. Breeders should closely cooperate with scientists from various fields of science.

Deciphering the evolution and complexity of wheat germplasm from a genomic perspective

Bread wheat provides an essential fraction of the daily calorific intake for humanity. Due to its huge and complex genome, progress in studying on the wheat genome is substantially trailed behind those of the other two major crops, rice and maize, for at least a decade. With rapid advances in genome assembling and reduced cost of high-throughput sequencing, emerging de novo genome assemblies of wheat and whole-genome sequencing data are leading to a paradigm shift in wheat research. Here, we review recent progress in dissecting the complex genome and germplasm evolution of wheat since the release of the first high-quality wheat genome. New insights have been gained in the evolution of wheat germplasm during domestication and modern breeding progress, genomic variations at multiple scales contributing to the diversity of wheat germplasm, and complex transcriptional and epigenetic regulations of functional genes in polyploid wheat. Genomics databases and bioinformatics tools meeting the urgent needs of wheat genomics research are also summarized. The ever-increasing omics data, along with advanced tools and well-structured databases, are expected to accelerate deciphering the germplasm and gene resources in wheat for future breeding advances.

GreenHill: a de novo chromosome-level scaffolding and phasing tool using Hi-C

Chromosome-level haplotype-resolved genome assembly is an important resource in molecular biology. However, current de novo haplotype assemblers require parental data or reference genomes and often fail to provide chromosome-level results. We present GreenHill, a novel scaffolding and phasing tool that considers various assemblers' contigs as input to reconstruct chromosome-level haplotypes using Hi-C without parental or reference data. Its unique functions include new error correction based on Hi-C contacts and the simultaneous use of Hi-C and long reads. Benchmarks reveal that GreenHill outperforms other approaches in contiguity and phasing accuracy, and the majority of chromosome arms are entirely phased.

Genetic Analysis and Physical Mapping of Oat Adult Plant Resistance Loci Against Puccinia coronata f. sp. avenae

Six quantitative trait loci (QTLs) for adult plant resistance against oat crown rust (Puccinia coronata f. sp. avenae) were identified from mapping three recombinant inbred populations. Using genotyping-by-sequencing with markers called against the OT3098 v1 reference genome, the QTLs were mapped on six different chromosomes: Chr1D, Chr4D, Chr5A, Chr5D, Chr7A, and Chr7C. Composite interval mapping with marker cofactor selection showed that the phenotypic variance explained by all identified QTLs for coefficient of infection range from 12.2 to 46.9%, whereas heritability estimates ranged from 0.11 to 0.38. The significant regions were narrowed down to intervals of 3.9 to 25 cM, equivalent to physical distances of 11 to 133 Mb. At least two flanking single-nucleotide polymorphism markers were identified within 10 cM of each QTL that could be used in marker-assisted introgression, pyramiding, and selection. The additive effects of the QTLs in each population were determined using single-nucleotide polymorphism haplotype data, which showed a significantly lower coefficient of infection in lines homozygous for the resistant alleles. Analysis of pairwise linkage disequilibrium also revealed high correlation of markers and presence of linkage blocks in the significant regions. To further facilitate marker-assisted breeding, polymerase chain reaction allelic competitive extension (PACE) markers for the adult plant resistance loci were developed. Putative candidate genes were also identified in each of the significant regions, which include resistance gene analogs that encode for kinases, ligases, and predicted receptors of avirulence proteins from pathogens.

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.

Wheat EARLY FLOWERING 3 affects heading date without disrupting circadian oscillations

Plant breeders have indirectly selected for variation at circadian-associated loci in many of the world's major crops, when breeding to increase yield and improve crop performance. Using an eight-parent Multiparent Advanced Generation Inter-Cross (MAGIC) population, we investigated how variation in circadian clock-associated genes contributes to the regulation of heading date in UK and European winter wheat (Triticum aestivum) varieties. We identified homoeologues of EARLY FLOWERING 3 (ELF3) as candidates for the Earliness per se (Eps) D1 and B1 loci under field conditions. We then confirmed a single-nucleotide polymorphism within the coding region of TaELF3-B1 as a candidate polymorphism underlying the Eps-B1 locus. We found that a reported deletion at the Eps-D1 locus encompassing TaELF3-D1 is, instead, an allele that lies within an introgression region containing an inversion relative to the Chinese Spring D genome. Using Triticum turgidum cv. Kronos carrying loss-of-function alleles of TtELF3, we showed that ELF3 regulates heading, with loss of a single ELF3 homoeologue sufficient to alter heading date. These studies demonstrated that ELF3 forms part of the circadian oscillator; however, the loss of all homoeologues was required to affect circadian rhythms. Similarly, loss of functional LUX ARRHYTHMO (LUX) in T. aestivum, an orthologue of a protein partner of Arabidopsis (Arabidopsis thaliana) ELF3, severely disrupted circadian rhythms. ELF3 and LUX transcripts are not co-expressed at dusk, suggesting that the structure of the wheat circadian oscillator might differ from that of Arabidopsis. Our demonstration that alterations to ELF3 homoeologues can affect heading date separately from effects on the circadian oscillator suggests a role for ELF3 in cereal photoperiodic responses that could be selected for without pleiotropic deleterious alterations to circadian rhythms.

A 'wiring diagram' for sink strength traits impacting wheat yield potential

Identifying traits for improving sink strength is a bottleneck to increasing wheat yield. The interacting processes determining sink strength and yield potential are reviewed and visualized in a set of 'wiring diagrams', covering critical phases of development (and summarizing known underlying genetics). Using this framework, we reviewed and assembled the main traits determining sink strength and identified research gaps and potential hypotheses to be tested for achieving gains in sink strength. In pre-anthesis, grain number could be increased through: (i) enhanced spike growth associated with optimized floret development and/or a reduction in specific stem-internode lengths and (ii) improved fruiting efficiency through an accelerated rate of floret development, improved partitioning between spikes, or optimized spike cytokinin levels. In post-anthesis, grain, sink strength could be augmented through manipulation of grain size potential via ovary size and/or endosperm cell division and expansion. Prospects for improving spike vascular architecture to support all rapidly growing florets, enabling the improved flow of assimilate, are also discussed. Finally, we considered the prospects for enhancing grain weight realization in relation to genetic variation in stay-green traits as well as stem carbohydrate remobilization. The wiring diagrams provide a potential workspace for breeders and crop scientists to achieve yield gains in wheat and other field crops.

An autoactive NB-LRR gene causes Rht13 dwarfism in wheat

Semidwarfing genes have greatly increased wheat yields globally, yet the widely used gibberellin (GA)-insensitive genes Rht-B1b and Rht-D1b have disadvantages for seedling emergence. Use of the GA-sensitive semidwarfing gene Rht13 avoids this pleiotropic effect. Here, we show that Rht13 encodes a nucleotide-binding site/leucine-rich repeat (NB-LRR) gene. A point mutation in the semidwarf Rht-B13b allele autoactivates the NB-LRR gene and causes a height reduction comparable with Rht-B1b and Rht-D1b in diverse genetic backgrounds. The autoactive Rht-B13b allele leads to transcriptional up-regulation of pathogenesis-related genes including class III peroxidases associated with cell wall remodeling. Rht13 represents a new class of reduced height (Rht) gene, unlike other Rht genes, which encode components of the GA signaling or metabolic pathways. This discovery opens avenues to use autoactive NB-LRR genes as semidwarfing genes in a range of crop species, and to apply Rht13 in wheat breeding programs using a perfect genetic marker.

Deciphering diversity at er loci for diversification of powdery mildew resistance in pea

Agricultural biotechnology aims to scrutinize the field crops which feed half of the world's population by improving their agronomic traits using various biotechnological tools. Pea- an important cash crop, rich in nutrients, but frequently infected with powdery mildew (fungal disease caused by Erysiphe pisi) that destroys the whole crop and causes economic loss for growers. We, therefore, targeted this research to find the pathogen-resistant pea lines and further decipher the diversity at er locus among resistant pea lines. Screening for resistant pea lines was done with Erysiphe pisi isolates (Genebank submission: KX455922.1) under the net house and greenhouse conditions. Molecular studies revealed that the Erysiphe resistant (er1) gene was present in 40 lines out of selected 50 pea lines and the mutational character was conferred up to 36 genotypes with 11 haplotype groups. The haplotype (gene) diversity (Hd) was found to be 0.5571 ± 0.099 SD and the nucleotide diversity (Pi) was 0.0160 ± 0.0042 SD Majority of resistant lines (67%) occurred in Hap-1, other remaining haplotypes (Hap 2-10) having 33% resistant lines, each showing characteristic nucleotide substitutions with respect to reference PsMLO1 gene; genotypes from these divergent haplotypes can be used in pea resistance breeding to avoid genetic homogeneity and genetic vulnerability.

Discovery and Validation of Grain Shape Loci in U.S. Rice Germplasm Through Haplotype Characterization

Rice grain shape is a major determinant of rice market value and the end-use. We mapped quantitative trait loci (QTL) for grain shape traits in a bi-parental recombinant inbred line population (Trenasse/Jupiter) and discovered two major grain length QTLs—qGL3.1 and qGL7.1. Previously, a major grain shape gene GS3 was reported in the qGL3.1 region and grain length gene GL7 was reported to be encompassing qGL7.1 locus. The re-sequencing SNP data on the International Rice Research Institute (IRRI) 3K Rice Genome Project (RGP) panel were obtained from the IRRI SNP-Seek database for both genes and haplotype diversity was characterized for each gene in this diverse panel. United States rice germplasm was not well represented in the IRRI 3K RGP database. Therefore, a minimum SNP set was identified for each gene that could differentiate all the characterized haplotypes. These haplotypes in the 3K RGP panel were screened across 323 elite U.S. genotypes using the minimum SNP set. The screening of haplotypes and phenotype association confirmed the role of GS3 under qGL3.1. However, screening of the GL7 haplotypes in the U.S. germplasm panel showed that GL7 did not play a role in qGL7.1, and in addition, GL7.1 did not segregate in the Trenasse/Jupiter RIL population. This concluded that qGL7.1 is a novel QTL discovered on chr7 for grain shape in the Trenasse/Jupiter RIL population. A high-throughput KASP-based SNP marker for each locus (GS3 and qGL7.1) was identified and validated in elite U.S. rice germplasm to be used in an applied rice breeding program.

TaNAC020 homoeologous genes are associated with higher thousand kernel weight and kernel length in Chinese wheat

NAC proteins constitute one of the largest plant-specific transcription factor (TF) families and play significant roles in plant growth and development. In the present study, three TaNAC020 homoeologous genes located on chromosomes 7A, 7B, and 7D were isolated from wheat (Triticum aestivum L.). TaNAC020s were predominantly expressed in developing grains. The developed transgenic rice lines for TaNAC020-B showed higher starch density and lower amylose contents than those of the wild type (WT). Sequence polymorphism studies showed seven and eight SNPs in TaNAC020-A/B, making three and two haplotypes, respectively. No sequence polymorphism was identified in TaNAC020-D. Association analysis revealed that HAP-2 of TaNAC020-A and TaNAC020-B was the favored haplotype for higher thousand kernel weight and length. Geographic distribution and allelic frequency showed that our favored haplotype experienced strong selection in China, and likewise, diversity increased in TaNAC020s during wheat polyploidization. The results obtained in this study demonstrate that TaNAC020s positively influence starch synthesis and accumulation and are one of the key regulators of the kernel (seed) size and kernel number and have the potential for utilization in wheat breeding to improve grain yield. Molecular markers developed in this study stand instrumental in marker-assisted selection for genetic improvement and germplasm enhancement in wheat.

Dispersed emergence and protracted domestication of polyploid wheat uncovered by mosaic ancestral haploblock inference

Major crops are all survivors of domestication bottlenecks. Studies have focused on the genetic loci related to the domestication syndrome, while the contribution of ancient haplotypes remains largely unknown. Here, an ancestral genomic haploblock dissection method is developed and applied to a resequencing dataset of 386 tetraploid/hexaploid wheat accessions, generating a pan-ancestry haploblock map. Together with cytoplastic evidences, we reveal that domesticated polyploid wheat emerged from the admixture of six founder wild emmer lineages, which contributed the foundation of ancestral mosaics. The key domestication-related loci, originated over a wide geographical range, were gradually pyramided through a protracted process. Diverse stable-inheritance ancestral haplotype groups of the chromosome central zone are identified, revealing the expanding routes of wheat and the trends of modern wheat breeding. Finally, an evolution model of polyploid wheat is proposed, highlighting the key role of wild-to-crop and interploidy introgression, that increased genomic diversity following bottlenecks introduced by domestication and polyploidization.

The contribution of ancient haplotypes to domestication is largely unknown. Here, the authors develop an ancestral genomic haploblock dissection method to generate a mosaic pan-ancestry genomic map and reveal that the domesticated polyploidy wheat emerged from the admixture of six founder wild emmer linages.

Characterization of Expression and Epigenetic Features of Core Genes in Common Wheat

The availability of multiple wheat genome sequences enables us to identify core genes and characterize their genetic and epigenetic features, thereby advancing our understanding of their biological implications within individual plant species. It is, however, largely understudied in wheat. To this end, we reanalyzed genome sequences from 16 different wheat varieties and identified 62,299 core genes. We found that core and non-core genes have different roles in subgenome differentiation. Meanwhile, according to their expression profiles, these core genes can be classified into genes related to tissue development and stress responses, including 3376 genes highly expressed in both spikelets and at high temperatures. After associating with six histone marks and open chromatin, we found that these core genes can be divided into eight sub-clusters with distinct epigenomic features. Furthermore, we found that ca. 51% of the expressed transcription factors (TFs) were marked with both H3K27me3 and H3K4me3, indicative of the bivalency feature, which can be involved in tissue development through the TF-centered regulatory network. Thus, our study provides a valuable resource for the functional characterization of core genes in stress responses and tissue development in wheat.

Chromosome-Level Haplotype Assembly for Equus asinu

Background: Haplotype provides significant insights into understanding genomes at both individual and population levels. However, research on many non-model organisms is still based on independent genetic variations due to the lack of haplotype.

Results: We conducted haplotype assembling for Equus asinu, a non-model organism that plays a vital role in human civilization. We described the hybrid single individual assembled haplotype of the Dezhou donkey based on the high-depth sequencing data from single-molecule real-time sequencing (×30), Illumina short-read sequencing (×211), and high-throughput chromosome conformation capture (×56). We assembled a near-complete haplotype for the high-depth sequenced Dezhou donkey individual and a phased cohort for the resequencing data of the donkey population.

Conclusion: Here, we described the complete chromosome-scale haplotype of the Dezhou donkey with more than a 99.7% phase rate. We further phased a cohort of 156 donkeys to form a donkey haplotype dataset with more than 39 million genetic variations.

High-resolution detection of quantitative trait loci for seven important yield-related traits in wheat (Triticum aestivum L.) using a high-density SLAF-seq genetic map

BACKGROUND

Yield-related traits including thousand grain weight (TGW), grain number per spike (GNS), grain width (GW), grain length (GL), plant height (PH), spike length (SL), and spikelet number per spike (SNS) are greatly associated with grain yield of wheat (Triticum aestivum L.). To detect quantitative trait loci (QTL) associated with them, 193 recombinant inbred lines derived from two elite winter wheat varieties Chuanmai42 and Chuanmai39 were employed to perform QTL mapping in six/eight environments.

RESULTS

A total of 30 QTLs on chromosomes 1A, 1B, 1D, 2A, 2B, 2D, 3A, 4A, 5A, 5B, 6A, 6D, 7A, 7B and 7D were identified. Among them, six major QTLs QTgw.cib-6A.1, QTgw.cib-6A.2, QGw.cib-6A, QGl.cib-3A, QGl.cib-6A, and QSl.cib-2D explaining 5.96-23.75% of the phenotypic variance were detected in multi-environments and showed strong and stable effects on corresponding traits. Three QTL clusters on chromosomes 2D and 6A containing 10 QTLs were also detected, which showed significant pleiotropic effects on multiple traits. Additionally, three Kompetitive Allele Specific PCR (KASP) markers linked with five of these major QTLs were developed. Candidate genes of QTgw.cib-6A.1/QGl.cib-6A and QGl.cib-3A were analyzed based on the spatiotemporal expression patterns, gene annotation, and orthologous search.

CONCLUSIONS

Six major QTLs for TGW, GL, GW and SL were detected. Three KASP markers linked with five of these major QTLs were developed. These QTLs and KASP markers will be useful for elucidating the genetic architecture of grain yield and developing new wheat varieties with high and stable yield in wheat.

Long-read and chromosome-scale assembly of the hexaploid wheat genome achieves high resolution for research and breeding

BACKGROUND

The sequencing of the wheat (Triticum aestivum) genome has been a methodological challenge for many years owing to its large size (15.5 Gb), repeat content, and hexaploidy. Many initiatives aiming at obtaining a reference genome of cultivar Chinese Spring have been launched in the past years and it was achieved in 2018 as the result of a huge effort to combine short-read sequencing with many other resources. Reference-quality genome assemblies were then produced for other accessions, but the rapid evolution of sequencing technologies offers opportunities to reach high-quality standards at lower cost.

RESULTS

Here, we report on an optimized procedure based on long reads produced on the Oxford Nanopore Technology PromethION device to assemble the genome of the French bread wheat cultivar Renan.

CONCLUSIONS

We provide the most contiguous chromosome-scale assembly of a bread wheat genome to date. Coupled with an annotation based on RNA-sequencing data, this resource will be valuable for the crop community and will facilitate the rapid selection of agronomically important traits. We also provide a framework to generate high-quality assemblies of complex genomes using ONT.

Trade-offs in the genetic control of functional and nutritional quality traits in UK winter wheat

A complex network of trade-offs exists between wheat quality and nutritional traits. We investigated the correlated relationships among several milling and baking traits as well as mineral density in refined white and whole grain flour. Our aim was to determine their pleiotropic genetic control in a multi-parent population over two trial years with direct application to practical breeding. Co-location of major quantitative trait loci (QTL) and principal component based multi-trait QTL mapping increased the power to detect QTL and revealed pleiotropic effects explaining many complementary and antagonistic trait relationships. High molecular weight glutenin subunit genes explained much of the heritable variation in important dough rheology traits, although additional QTL were detected. Several QTL, including one linked to the TaGW2 gene, controlled grain size and increased flour extraction rate. The semi-dwarf Rht-D1b allele had a positive effect on Hagberg falling number, but reduced grain size, specific weight, grain protein content and flour water absorption. Mineral nutrient concentrations were lower in Rht-D1b lines for many elements, in wholemeal and white flour, but potassium concentration was higher in Rht-D1b lines. The presence of awns increased calcium content without decreasing extraction rate, despite the negative correlation between these traits. QTL were also found that affect the relative concentrations of key mineral nutrients compared to phosphorus which may help increase bioavailability without associated anti-nutritional effects of phytic acid. Taken together these results demonstrate the potential for marker-based selection to optimise trait trade-offs and enhance wheat nutritional value by considering pleiotropic genetic effects across multiple traits.

Genome sequences of three Aegilops species of the section Sitopsis reveal phylogenetic relationships and provide resources for wheat improvement

Aegilops is a close relative of wheat (Triticum spp.), and Aegilops species in the section Sitopsis represent a rich reservoir of genetic diversity for the improvement of wheat. To understand their diversity and advance their utilization, we produced whole-genome assemblies of Aegilops longissima and Aegilops speltoides. Whole-genome comparative analysis, along with the recently sequenced Aegilops sharonensis genome, showed that the Ae. longissima and Ae. sharonensis genomes are highly similar and are most closely related to the wheat D subgenome. By contrast, the Ae. speltoides genome is more closely related to the B subgenome. Haplotype block analysis supported the idea that Ae. speltoides genome is closest to the wheat B subgenome, and highlighted variable and similar genomic regions between the three Aegilops species and wheat. Genome-wide analysis of nucleotide-binding leucine-rich repeat (NLR) genes revealed species-specific and lineage-specific NLR genes and variants, demonstrating the potential of Aegilops genomes for wheat improvement.

ggComp enables dissection of germplasm resources and construction of a multiscale germplasm network in wheat

Accurate germplasm characterization is a vital step for accelerating crop genetic improvement, which remains largely infeasible for crops such as bread wheat (Triticum aestivum L.), which has a complex genome that undergoes frequent introgression and contains many structural variations. Here, we propose a genomic strategy called ggComp, which integrates resequencing data with copy number variations and stratified single-nucleotide polymorphism densities to enable unsupervised identification of pairwise germplasm resource-based Identity-By-Descent (gIBD) blocks. The reliability of ggComp was verified in wheat cultivar Nongda5181 by dissecting parental-descent patterns represented by inherited genomic blocks. With gIBD blocks identified among 212 wheat accessions, we constructed a multi-scale genomic-based germplasm network. At the whole-genome level, the network helps to clarify pedigree relationship, demonstrate genetic flow, and identify key founder lines. At the chromosome level, we were able to trace the utilization of 1RS introgression in modern wheat breeding by hitchhiked segments. At the single block scale, the dissected germplasm-based haplotypes nicely matched with previously identified alleles of "Green Revolution" genes and can guide allele mining and dissect the trajectory of beneficial alleles in wheat breeding. Our work presents a model-based framework for precisely evaluating germplasm resources with genomic data. A database, WheatCompDB (http://wheat.cau.edu.cn/WheatCompDB/), is available for researchers to exploit the identified gIBDs with a multi-scale network.

Genomic prediction of morphometric and colorimetric traits in Solanaceous fruits

Selection of high-performance lines with respect to traits of interest is a key step in plant breeding. Genomic prediction allows to determine the genomic estimated breeding values of unseen lines for trait of interest using genetic markers, e.g. single-nucleotide polymorphisms (SNPs), and machine learning approaches, which can therefore shorten breeding cycles, referring to genomic selection (GS). Here, we applied GS approaches in two populations of Solanaceous crops, i.e. tomato and pepper, to predict morphometric and colorimetric traits. The traits were measured by using scoring-based conventional descriptors (CDs) as well as by Tomato Analyzer (TA) tool using the longitudinally and latitudinally cut fruit images. The GS performance was assessed in cross-validations of classification-based and regression-based machine learning models for CD and TA traits, respectively. The results showed the usage of TA traits and tag SNPs provide a powerful combination to predict morphology and color-related traits of Solanaceous fruits. The highest predictability of 0.89 was achieved for fruit width in pepper, with an average predictability of 0.69 over all traits. The multi-trait GS models are of slightly better predictability than single-trait models for some colorimetric traits in pepper. While model validation performs poorly on wild tomato accessions, the usage as many as one accession per wild species in the training set can increase the transferability of models to unseen populations for some traits (e.g. fruit shape for which predictability in unseen scenario increased from zero to 0.6). Overall, GS approaches can assist the selection of high-performance Solanaceous fruits in crop breeding.

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.

The 3366 chickpea genomes for research and breeding

Genome sequences provide an unprecedented resource to rapidly develop modern crops. A recent paper by Varshney et al. provides genome variation maps of 3366 chickpea accessions. Here, we highlight how this breakthrough research can fundamentally change breeding practices of chickpea and potentially other crops.

Wheat Breeding, Transcription Factories, and Genetic Interactions: New Perspectives

Epistatic interactions and negative heterosis have been shown to be associated with interchromosomal interactions in wheat. Physical gene-gene interactions between co-regulated genes clustered in "transcription factories" have been documented, and a genome-wide atlas of functionally paired, interacting regulatory elements and genes of wheat recently produced. Integration of these new studies on gene and regulatory element interactions, co-regulation of gene expression in "transcription factories," and epigenetics generates new perspectives for wheat breeding and trait enhancement.

Trend, population structure, and trait mapping from 15 years of national varietal trials of UK winter wheat

There are now a rich variety of genomic and genotypic resources available to wheat researchers and breeders. However, the generation of high-quality and field-relevant phenotyping data which is required to capture the complexities of gene × environment interactions remains a major bottleneck. Historical datasets from national variety performance trials (NVPT) provide sufficient dimensions, in terms of numbers of years and locations, to examine phenotypic trends and study gene × environment interactions. Using NVPT for winter wheat varieties grown in the United Kingdom between 2002 and 2017, we examined temporal trends for eight traits related to yield, adaptation, and grain quality performance. We show a non-stationary linear trend for yield, grain protein content, Hagberg Falling Number (HFN), and days to ripening. Our data also show high environmental stability for yield, grain protein content, and specific weight in UK winter wheat varieties and high environmental sensitivity for HFN. We also show that UK varieties released within this period cluster into four main population groups. Using the historical NVPT data in a genome-wide association analysis, we uncovered a significant marker-trait association peak on wheat chromosome 6A spanning the NAM-A1 gene that have been previously associated with early senescence. Together, our results show the value of utilizing the data routinely collected during national variety evaluation process for examining breeding progress and the genetic architecture of important traits.

Different Shades of Kale-Approaches to Analyze Kale Variety Interrelations

Brassica oleracea is a vegetable crop with an amazing morphological diversity. Among the various crops derived from B. oleracea, kale has been in the spotlight globally due to its various health-benefitting compounds and many different varieties. Knowledge of the existing genetic diversity is essential for the improved breeding of kale. Here, we analyze the interrelationships, population structures, and genetic diversity of 72 kale and cabbage varieties by extending our previous diversity analysis and evaluating the use of summed potential lengths of shared haplotypes (SPLoSH) as a new method for such analyses. To this end, we made use of the high-density Brassica 60K SNP array, analyzed SNPs included in an available Brassica genetic map, and used these resources to generate and evaluate the information from SPLoSH data. With our results we could consistently differentiate four groups of kale across all analyses: the curly kale varieties, Italian, American, and Russian varieties, as well as wild and cultivated types. The best results were achieved by using SPLoSH information, thus validating the use of this information in improving analyses of interrelations in kale. In conclusion, our definition of kale includes the curly varieties as the kales in a strict sense, regardless of their origin. These results contribute to a better understanding of the huge diversity of kale and its interrelations.

Delaying or delivering: identification of novel NAM-1 alleles that delay senescence to extend wheat grain fill duration

Senescence is a complex trait under genetic and environmental control, in which resources are remobilized from vegetative tissue into grain. Delayed senescence, or 'staygreen' traits, can confer stress tolerance, with extended photosynthetic activity hypothetically sustaining grain filling. The genetics of senescence regulation are largely unknown, with senescence variation often correlated with phenological traits. Here, we confirm staygreen phenotypes of two Triticum aestivum cv. Paragon ethyl methane sulfonate mutants previously identified during a forward genetic screen and selected for their agronomic performance, similar phenology, and differential senescence phenotypes. Grain filling experiments confirmed a positive relationship between onset of senescence and grain fill duration, reporting an associated ~14% increase in final dry grain weight for one mutant (P<0.05). Recombinant inbred line (RIL) populations segregating for the timing of senescence were developed for trait mapping purposes and phenotyped over multiple years under field conditions. Quantification and comparison of senescence metrics aided RIL selection, facilitating exome capture-enabled bulk segregant analysis (BSA). Using BSA we mapped our two staygreen traits to two independent, dominant, loci of 4.8 and 16.7 Mb in size encompassing 56 and 142 genes, respectively. Combining association analysis with variant effect prediction, we identified single nucleotide polymorphisms encoding self-validating mutations located in NAM-1 homoeologues, which we propose as gene candidates.

A Genome-Wide Association Study Pinpoints Quantitative Trait Genes for Plant Height, Heading Date, Grain Quality, and Yield in Rye (Secale cereale L.)

Rye is the only cross-pollinating Triticeae crop species. Knowledge of rye genes controlling complex-inherited traits is scarce, which, currently, largely disables the genomics assisted introgression of untapped genetic variation from self-incompatible germplasm collections in elite inbred lines for hybrid breeding. We report on the first genome-wide association study (GWAS) in rye based on the phenotypic evaluation of 526 experimental hybrids for plant height, heading date, grain quality, and yield in 2 years and up to 19 environments. We established a cross-validated NIRS calibration model as a fast, effective, and robust analytical method to determine grain quality parameters. We observed phenotypic plasticity in plant height and tiller number as a resource use strategy of rye under drought and identified increased grain arabinoxylan content as a striking phenotype in osmotically stressed rye. We used DArTseq™ as a genotyping-by-sequencing technology to reduce the complexity of the rye genome. We established a novel high-density genetic linkage map that describes the position of almost 19k markers and that allowed us to estimate a low genome-wide LD based on the assessed genetic diversity in elite germplasm. We analyzed the relationship between plant height, heading date, agronomic, as well as grain quality traits, and genotype based on 20k novel single-nucleotide polymorphism markers. In addition, we integrated the DArTseq™ markers in the recently established 'Lo7' reference genome assembly. We identified cross-validated SNPs in 'Lo7' protein-coding genes associated with all traits studied. These include associations of the WUSCHEL-related homeobox transcription factor DWT1 and grain yield, the DELLA protein gene SLR1 and heading date, the Ethylene overproducer 1-like protein gene ETOL1 and thousand-grain weight, protein and starch content, as well as the Lectin receptor kinase SIT2 and plant height. A Leucine-rich repeat receptor protein kinase and a Xyloglucan alpha-1,6-xylosyltransferase count among the cross-validated genes associated with water-extractable arabinoxylan content. This study demonstrates the power of GWAS, hybrid breeding, and the reference genome sequence in rye genetics research to dissect and identify the function of genes shaping genetic diversity in agronomic and grain quality traits of rye. The described links between genetic causes and phenotypic variation will accelerate genomics-enabled rye improvement.

Genome-wide association study uncovers major genetic loci associated with seed flooding tolerance in soybean

BACKGROUND

Seed flooding stress is one of the threatening environmental stressors that adversely limits soybean at the germination stage across the globe. The knowledge on the genetic basis underlying seed-flooding tolerance is limited. Therefore, we performed a genome-wide association study (GWAS) using 34,718 single nucleotide polymorphism (SNPs) in a panel of 243 worldwide soybean collections to identify genetic loci linked to soybean seed flooding tolerance at the germination stage.

RESULTS

In the present study, GWAS was performed with two contrasting models, Mixed Linear Model (MLM) and Multi-Locus Random-SNP-Effect Mixed Linear Model (mrMLM) to identify significant SNPs associated with electrical conductivity (EC), germination rate (GR), shoot length (ShL), and root length (RL) traits at germination stage in soybean. With MLM, a total of 20, 40, 4, and 9 SNPs associated with EC, GR, ShL and RL, respectively, whereas in the same order mrMLM detected 27, 17, 13, and 18 SNPs. Among these SNPs, two major SNPs, Gm_08_11971416, and Gm_08_46239716 were found to be consistently connected with seed-flooding tolerance related traits, namely EC and GR across two environments. We also detected two SNPs, Gm_05_1000479 and Gm_01_53535790 linked to ShL and RL, respectively. Based on Gene Ontology enrichment analysis, gene functional annotations, and protein-protein interaction network analysis, we predicted eight candidate genes and three hub genes within the regions of the four SNPs with Cis-elements in promoter regions which may be involved in seed-flooding tolerance in soybeans and these warrant further screening and functional validation.

CONCLUSIONS

Our findings demonstrate that GWAS based on high-density SNP markers is an efficient approach to dissect the genetic basis of complex traits and identify candidate genes in soybean. The trait associated SNPs could be used for genetic improvement in soybean breeding programs. The candidate genes could help researchers better understand the molecular mechanisms underlying seed-flooding stress tolerance in soybean.

Rewilding crops for climate resilience: economic analysis and de novo domestication strategies

To match predicted population growth, annual food production should be doubled by 2050. This is not achievable by current agronomical and breeding practices, due to the impact of climate changes and associated abiotic stresses on agricultural production systems. Here, we analyze the impact of global climate trends on crop productivity and show that the overall loss in crop production from climate-driven abiotic stresses may exceed US$170 billion year-1 and represents a major threat to global food security. We also show that abiotic stress tolerance had been present in wild progenitors of modern crops but was lost during their domestication. We argue for a major shift in our paradigm of crop breeding, focusing on climate resilience, and call for a broader use of wild relatives as a major tool in this process. We argue that, while molecular tools are currently in place to harness the potential of climate-resilient genes present in wild relatives, the complex polygenic nature of tolerance traits remains a major bottleneck in this process. Future research efforts should be focused not only on finding appropriate wild relatives but also on development of efficient cell-based high-throughput phenotyping platforms allowing assessment of the in planta operation of key genes.

Fast-forward breeding for a food-secure world

Crop production systems need to expand their outputs sustainably to feed a burgeoning human population. Advances in genome sequencing technologies combined with efficient trait mapping procedures accelerate the availability of beneficial alleles for breeding and research. Enhanced interoperability between different omics and phenotyping platforms, leveraged by evolving machine learning tools, will help provide mechanistic explanations for complex plant traits. Targeted and rapid assembly of beneficial alleles using optimized breeding strategies and precise genome editing techniques could deliver ideal crops for the future. Realizing desired productivity gains in the field is imperative for securing an adequate future food supply for 10 billion people.

Crossover-active regions of the wheat genome are distinguished by DMC1, the chromosome axis, H3K27me3, and signatures of adaptation

The hexaploid bread wheat genome comprises over 16 gigabases of sequence across 21 chromosomes. Meiotic crossovers are highly polarized along the chromosomes, with elevation in the gene-dense distal regions and suppression in the Gypsy retrotransposon-dense centromere-proximal regions. We profiled the genomic landscapes of the meiotic recombinase DMC1 and the chromosome axis protein ASY1 in wheat and investigated their relationships with crossovers, chromatin state, and genetic diversity. DMC1 and ASY1 chromatin immunoprecipitation followed by sequencing (ChIP-seq) revealed strong co-enrichment in the distal, crossover-active regions of the wheat chromosomes. Distal ChIP-seq enrichment is consistent with spatiotemporally biased cytological immunolocalization of DMC1 and ASY1 close to the telomeres during meiotic prophase I. DMC1 and ASY1 ChIP-seq peaks show significant overlap with genes and transposable elements in the Mariner and Mutator superfamilies. However, DMC1 and ASY1 ChIP-seq peaks were detected along the length of each chromosome, including in low-crossover regions. At the fine scale, crossover elevation at DMC1 and ASY1 peaks and genes correlates with enrichment of the Polycomb histone modification H3K27me3. This indicates a role for facultative heterochromatin, coincident with high DMC1 and ASY1, in promoting crossovers in wheat and is reflected in distalized H3K27me3 enrichment observed via ChIP-seq and immunocytology. Genes with elevated crossover rates and high DMC1 and ASY1 ChIP-seq signals are overrepresented for defense-response and immunity annotations, have higher sequence polymorphism, and exhibit signatures of selection. Our findings are consistent with meiotic recombination promoting genetic diversity, shaping host–pathogen co-evolution, and accelerating adaptation by increasing the efficiency of selection.

Data-driven decentralized breeding increases prediction accuracy in a challenging crop production environment

Crop breeding must embrace the broad diversity of smallholder agricultural systems to ensure food security to the hundreds of millions of people living in challenging production environments. This need can be addressed by combining genomics, farmers' knowledge, and environmental analysis into a data-driven decentralized approach (3D-breeding). We tested this idea as a proof-of-concept by comparing a durum wheat (Triticum durum Desf.) decentralized trial distributed as incomplete blocks in 1,165 farmer-managed fields across the Ethiopian highlands with a benchmark representing genomic prediction applied to conventional breeding. We found that 3D-breeding could double the prediction accuracy of the benchmark. 3D-breeding could identify genotypes with enhanced local adaptation providing superior productive performance across seasons. We propose this decentralized approach to leverage the diversity in farmer fields and complement conventional plant breeding to enhance local adaptation in challenging crop production environments.

Positional-based cloning 'fail-safe' approach is overpowered by wheat chromosome structural variation

Positional-based cloning is a foundational method for understanding the genes and gene networks that control valuable agronomic traits such as grain yield components. In this study, we sought to positionally clone the causal genetic variant of a 1000-grain weight (TGW) quantitative trait loci (QTL) on wheat (Triticum aestivum L.) chromosome arm 5AL. We developed heterogenous inbred families (HIFs) (>5,000 plants) for enhanced genotypic resolution and fine-mapped the QTL to a 10-Mbp region. The transcriptome of developing grains from positive and negative control HIF haplotypes revealed presence-absence chromosome arm 5AS structural variation and unexpectedly no differential expression of genes within the chromosome arm 5AL candidate region. Evaluation of genomic, transcriptomic, and phenotypic data, and predicted function of genes, identified that the 5AL QTL was the result of strong linkage disequilibrium (LD) with chromosome arm 5AS presence or absence (HIF r² = 0.91). Structural variation is common in wheat, and our results highlight that the redundant polyploid genome's masking of such variation is a significant barrier to positional cloning. We propose recommendations for more efficient and robust detection of structural variation, including transitioning from a single nucleotide polymorphism (SNP) to a haplotype-based approach to identify positional cloning targets. We also present nine candidate genes for grain yield components based on chromosome arm 5AS presence or absence, which may unveil hidden variation of homoeolog dosage-dependent genes across the group five chromosome short arms. Taken together, our discovery demonstrates the phenotypic resiliency of polyploid genomic structural variation and highlights a considerable challenge to routine positional cloning in wheat.

Utilization of a Wheat50K SNP Microarray-Derived High-Density Genetic Map for QTL Mapping of Plant Height and Grain Traits in Wheat

Plant height is significantly correlated with grain traits, which is a component of wheat yield. The purpose of this study is to investigate the main quantitative trait loci (QTLs) that control plant height and grain-related traits in multiple environments. In this study, we constructed a high-density genetic linkage map using the Wheat50K SNP Array to map QTLs for these traits in 198 recombinant inbred lines (RILs). The two ends of the chromosome were identified as recombination-rich areas in all chromosomes except chromosome 1B. Both the genetic map and the physical map showed a significant correlation, with a correlation coefficient between 0.63 and 0.99. However, there was almost no recombination between 1RS and 1BS. In terms of plant height, 1RS contributed to the reduction of plant height by 3.43 cm. In terms of grain length, 1RS contributed to the elongation of grain by 0.11 mm. A total of 43 QTLs were identified, including eight QTLs for plant height (PH), 11 QTLs for thousand grain weight (TGW), 15 QTLs for grain length (GL), and nine QTLs for grain width (GW), which explained 1.36–33.08% of the phenotypic variation. Seven were environment-stable QTLs, including two loci (Qph.nwafu-4B and Qph.nwafu-4D) that determined plant height. The explanation rates of phenotypic variation were 7.39–12.26% and 20.11–27.08%, respectively. One QTL, Qtgw.nwafu-4B, which influenced TGW, showed an explanation rate of 3.43–6.85% for phenotypic variation. Two co-segregating KASP markers were developed, and the physical locations corresponding to KASP_AX-109316968 and KASP_AX-109519968 were 25.888344 MB and 25.847691 MB, respectively. Qph.nwafu-4B, controlling plant height, and Qtgw.nwafu-4B, controlling TGW, had an obvious linkage relationship, with a distance of 7–8 cM. Breeding is based on molecular markers that control plant height and thousand-grain weight by selecting strains with low plant height and large grain weight. Another QTL, Qgw.nwafu-4D, which determined grain width, had an explanation rate of 3.43–6.85%. Three loci that affected grain length were Qgl.nwafu-5A, Qgl.nwafu-5D.2, and Qgl.nwafu-6B, illustrating the explanation rates of phenotypic variation as 6.72–9.59%, 5.62–7.75%, and 6.68–10.73%, respectively. Two QTL clusters were identified on chromosomes 4B and 4D.

Designing Future Crops: Genomics-Assisted Breeding Comes of Age

Over the past decade, genomics-assisted breeding (GAB) has been instrumental in harnessing the potential of modern genome resources and characterizing and exploiting allelic variation for germplasm enhancement and cultivar development. Sustaining GAB in the future (GAB 2.0) will rely upon a suite of new approaches that fast-track targeted manipulation of allelic variation for creating novel diversity and facilitate their rapid and efficient incorporation in crop improvement programs. Genomic breeding strategies that optimize crop genomes with accumulation of beneficial alleles and purging of deleterious alleles will be indispensable for designing future crops. In coming decades, GAB 2.0 is expected to play a crucial role in breeding more climate-smart crop cultivars with higher nutritional value in a cost-effective and timely manner.

Crop breeding for a changing climate: integrating phenomics and genomics with bioinformatics

Safeguarding crop yields in a changing climate requires bioinformatics advances in harnessing data from vast phenomics and genomics datasets to translate research findings into climate smart crops in the field. Climate change and an additional 3 billion mouths to feed by 2050 raise serious concerns over global food security. Crop breeding and land management strategies will need to evolve to maximize the utilization of finite resources in coming years. High-throughput phenotyping and genomics technologies are providing researchers with the information required to guide and inform the breeding of climate smart crops adapted to the environment. Bioinformatics has a fundamental role to play in integrating and exploiting this fast accumulating wealth of data, through association studies to detect genomic targets underlying key adaptive climate-resilient traits. These data provide tools for breeders to tailor crops to their environment and can be introduced using advanced selection or genome editing methods. To effectively translate research into the field, genomic and phenomic information will need to be integrated into comprehensive clade-specific databases and platforms alongside accessible tools that can be used by breeders to inform the selection of climate adaptive traits. Here we discuss the role of bioinformatics in extracting, analysing, integrating and managing genomic and phenomic data to improve climate resilience in crops, including current, emerging and potential approaches, applications and bottlenecks in the research and breeding pipeline.

Limited haplotype diversity underlies polygenic trait architecture across 70 years of wheat breeding

Background

Selection has dramatically shaped genetic and phenotypic variation in bread wheat. We can assess the genomic basis of historical phenotypic changes, and the potential for future improvement, using experimental populations that attempt to undo selection through the randomizing effects of recombination.

Results

We bred the NIAB Diverse MAGIC multi-parent population comprising over 500 recombinant inbred lines, descended from sixteen historical UK bread wheat varieties released between 1935 and 2004. We sequence the founders’ genes and promoters by capture, and the MAGIC population by low-coverage whole-genome sequencing. We impute 1.1 M high-quality SNPs that are over 99% concordant with array genotypes. Imputation accuracy only marginally improves when including the founders’ genomes as a haplotype reference panel. Despite capturing 73% of global wheat genetic polymorphism, 83% of genes cluster into no more than three haplotypes. We phenotype 47 agronomic traits over 2 years and map 136 genome-wide significant associations, concentrated at 42 genetic loci with large and often pleiotropic effects. Around half of these overlap known quantitative trait loci. Most traits exhibit extensive polygenicity, as revealed by multi-locus shrinkage modelling.

Conclusions

Our results are consistent with a gene pool of low haplotypic diversity, containing few novel loci of large effect. Most past, and projected future, phenotypic changes arising from existing variation involve fine-scale shuffling of a few haplotypes to recombine dozens of polygenic alleles of small effect. Moreover, extensive pleiotropy means selection on one trait will have unintended consequences, exemplified by the negative trade-off between yield and protein content, unless selection and recombination can break unfavorable trait-trait associations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02354-7.

The wheat Seven in absentia gene is associated with increases in biomass and yield in hot climates

Wheat (Triticum aestivum L.) productivity is severely reduced by high temperatures. Breeding of heat-tolerant cultivars can be achieved by identifying genes controlling physiological and agronomical traits when high temperatures occur and using these to select superior genotypes, but no gene underlying genetic variation for heat tolerance has previously been described. We advanced the positional cloning of qYDH.3BL, a quantitative trait locus (QTL) on bread wheat chromosome 3B associated with increased yield in hot and dry climates. The delimited genomic region contained 12 putative genes and a sequence variant in the promoter region of one gene, Seven in absentia, TaSINA. This was associated with the QTL's effects on early vigour, root growth, plant biomass, and yield components in two distinct wheat populations grown under various growth conditions. Near isogenic lines carrying the positive allele at qYDH.3BL underexpressed TaSINA and had increased vigour and water use efficiency early in development, as well as increased biomass, grain number, and grain weight following heat stress. A survey of worldwide distribution indicated that the positive allele became widespread from the 1950s through the CIMMYT wheat breeding programme but, to date, has been selected only in breeding programmes in Mexico and Australia.

Computational methods for chromosome-scale haplotype reconstruction

High-quality chromosome-scale haplotype sequences of diploid genomes, polyploid genomes, and metagenomes provide important insights into genetic variation associated with disease and biodiversity. However, whole-genome short read sequencing does not yield haplotype information spanning whole chromosomes directly. Computational assembly of shorter haplotype fragments is required for haplotype reconstruction, which can be challenging owing to limited fragment lengths and high haplotype and repeat variability across genomes. Recent advancements in long-read and chromosome-scale sequencing technologies, alongside computational innovations, are improving the reconstruction of haplotypes at the level of whole chromosomes. Here, we review recent and discuss methodological progress and perspectives in these areas.

A Duplicated Copy of the Meiotic Gene ZIP4 Preserves up to 50% Pollen Viability and Grain Number in Polyploid Wheat

Simple Summary

On wheat polyploidisation, the major meiotic gene ZIP4, duplicated and diverged, such that tetraploid and hexaploid wheat each carry three and four copies of ZIP4, respectively. Surprisingly, this study demonstrates that, in hexaploid wheat, despite the presence of the other three ZIP4 copies, the duplicated ZIP4 copy is required to prevent major abnormalities during meiosis. Although there is greater disruption of subsequent male rather than female fertility, the duplicated ZIP4 copy preserves up to 50% of the grain number. High wheat fertility is important since it is consumed by over 4.5 billion people on the planet, of whom 2.5 billion are dependent on it. This study highlights the potentially extraordinary value of the wheat ZIP4 duplication, mandating further studies to unravel the complexity of the ZIP4 phenotype in this global crop.

Abstract

Although most flowering plants are polyploid, little is known of how the meiotic process evolves after polyploidisation to stabilise and preserve fertility. On wheat polyploidisation, the major meiotic gene ZIP4 on chromosome 3B duplicated onto 5B and diverged (TaZIP4-B2). TaZIP4-B2 was recently shown to promote homologous pairing, synapsis and crossover, and suppress homoeologous crossover. We therefore suspected that these meiotic stabilising effects could be important for preserving wheat fertility. A CRISPR Tazip4-B2 mutant was exploited to assess the contribution of the 5B duplicated ZIP4 copy in maintaining pollen viability and grain setting. Analysis demonstrated abnormalities in 56% of meiocytes in the Tazip4-B2 mutant, with micronuclei in 50% of tetrads, reduced size in 48% of pollen grains and a near 50% reduction in grain number. Further studies showed that most of the reduced grain number occurred when Tazip4-B2 mutant plants were pollinated with the less viable Tazip4-B2 mutant pollen rather than with wild type pollen, suggesting that the stabilising effect of TaZIP4-B2 on meiosis has a greater consequence in subsequent male, rather than female gametogenesis. These studies reveal the extraordinary value of the wheat chromosome 5B TaZIP4-B2 duplication to agriculture and human nutrition. Future studies should further investigate the role of TaZIP4-B2 on female fertility and assess whether different TaZIP4-B2 alleles exhibit variable effects on meiotic stabilisation and/or resistance to temperature change.

Small "Nested" Introgressions from Wild Thinopyrum Species, Conferring Effective Resistance to Fusarium Diseases, Positively Impact Durum Wheat Yield Potential

Today wheat cultivation is facing rapidly changing climate scenarios and yield instability, aggravated by the spreading of severe diseases such as Fusarium head blight (FHB) and Fusarium crown rot (FCR). To obtain productive genotypes resilient to stress pressure, smart breeding approaches must be envisaged, including the exploitation of wild relatives. Here we report on the assessment of the breeding potential of six durum wheat-Thinopyrum spp. recombinant lines (RLs) obtained through chromosome engineering. They are characterized by having 23% or 28% of their 7AL chromosome arm replaced by a "nested" alien segment, composed of homoeologous group 7 chromosome fractions from Th. ponticum and Th. elongatum (=7el1L + 7EL) or from different Th. ponticum accessions (=7el1L + 7el2L). In addition to the 7el1L genes Lr19 + Yp (leaf rust resistance, and yellow pigment content, respectively), these recombinant lines (RLs) possess a highly effective QTL for resistance to FHB and FCR within their 7el2L or 7EL portion. The RLs, their null segregants and well-adapted and productive durum wheat cultivars were evaluated for 16 yield-related traits over two seasons under rainfed and irrigated conditions. The absence of yield penalties and excellent genetic stability of RLs was revealed in the presence of all the alien segment combinations. Both 7el2L and 7EL stacked introgressions had positive impacts on source and sink yield traits, as well as on the overall performance of RLs in conditions of reduced water availability. The four "nested" RLs tested in 2020 were among the top five yielders, overall representing good candidates to be employed in breeding programs to enhance crop security and safety.

QTL Mapping and Validation for Kernel Area and Circumference in Common Wheat via High-Density SNP-Based Genotyping

As an important component, 1,000 kernel weight (TKW) plays a significant role in the formation of yield traits of wheat. Kernel size is significantly positively correlated to TKW. Although numerous loci for kernel size in wheat have been reported, our knowledge on loci for kernel area (KA) and kernel circumference (KC) remains limited. In the present study, a recombinant inbred lines (RIL) population containing 371 lines genotyped using the Wheat55K SNP array was used to map quantitative trait loci (QTLs) controlling the KA and KC in multiple environments. A total of 54 and 44 QTLs were mapped by using the biparental population or multienvironment trial module of the inclusive composite interval mapping method, respectively. Twenty-two QTLs were considered major QTLs. BLAST analysis showed that major and stable QTLs QKc.sau-6A.1 (23.12-31.64 cM on 6A) for KC and QKa.sau-6A.2 (66.00-66.57 cM on 6A) for KA were likely novel QTLs, which explained 22.25 and 20.34% of the phenotypic variation on average in the 3 year experiments, respectively. Two Kompetitive allele-specific PCR (KASP) markers, KASP-AX-109894590 and KASP-AX-109380327, were developed and tightly linked to QKc.sau-6A.1 and QKa.sau-6A.2, respectively, and the genetic effects of the different genotypes in the RIL population were successfully confirmed. Furthermore, in the interval where QKa.sau-6A.2 was located on Chinese Spring and T. Turgidum ssp. dicoccoides reference genomes, only 11 genes were found. In addition, digenic epistatic QTLs also showed a significant influence on KC and KA. Altogether, the results revealed the genetic basis of KA and KC and will be useful for the marker-assisted selection of lines with different kernel sizes, laying the foundation for the fine mapping and cloning of the gene(s) underlying the stable QTLs detected in this study.

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 genome¹, and the lack of genome-assembly data for multiple wheat lines^2,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 stresses^4,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 Sm1⁶, 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.

Comparison of multiple genome assemblies from wheat reveals extensive diversity that results from the complex breeding history of wheat and provides a basis for further potential improvements to this important food crop.