Harnessing landrace diversity empowers wheat breeding.

Harnessing genetic diversity in major staple crops through the development of new breeding capabilities is essential to ensure food security1. Here we examined the genetic and phenotypic diversity of the A.E. Watkins landrace collection2 of bread wheat (Triticum aestivum), a major global cereal, through whole-genome re-sequencing (827 Watkins landraces and 208 modern cultivars) and in-depth field evaluation spanning a decade. We discovered that modern cultivars are derived from just two of the seven ancestral groups of wheat and maintain very long-range haplotype integrity. The remaining five groups represent untapped genetic sources, providing access to landrace-specific alleles and haplotypes for breeding. Linkage disequilibrium (LD) based haplotypes and association genetics analyses link Watkins genomes to the thousands of high-resolution quantitative trait loci (QTL), and significant marker-trait associations identified. Using these structured germplasm, genotyping and informatics resources, we revealed many Watkins-unique beneficial haplotypes that can confer superior traits in modern wheat. Furthermore, we assessed the phenotypic effects of 44,338 Watkins-unique haplotypes, introgressed from 143 prioritised QTL in the context of modern cultivars, bridging the gap between landrace diversity and current breeding. This study establishes a framework for systematically utilising genetic diversity in crop improvement to achieve sustainable food security.

Resolving a QTL complex for height, heading, and grain yield on chromosome 3A in bread wheat.

Crop height (Ht), heading date (Hd), and grain yield (GY) are inter-related in wheat. Independent manipulation of each is important for adaptation and performance. Validated quantitative trait loci (QTLs) for all three co-locate on chromosome 3A in the Avalon×Cadenza population, with increased Ht, Hd, and GY contributed by Cadenza. We asked if these are linked or pleiotropic effects using recombinant lines, and showed that Ht and Hd effects are independent. The Chinese Spring equivalent to the newly defined Ht interval contained a gene cluster involved in cell wall growth and displaying high levels of differential transcript expression. The Hd locus is larger and rearranged compared with the reference genome, but FT2 (Flowering Locus T2) is of particular interest. The Hd effect acted independently of photoperiod and vernalization, but did exhibit seasonal genotype×environment interaction. Recombinants were phenotyped for GY in replicated field experiments. GY was most associated with Cadenza alleles for later Hd, supporting physiological studies using the same lines proposing that 'late' alleles at this locus increase spike fertility and grain number (GN). The work has uncoupled height from heading and yield, and shown that one of very few validated GY QTLs in wheat is probably mediated by phenological variation.

Functional QTL mapping and genomic prediction of canopy height in wheat measured using a robotic field phenotyping platform.

Genetic studies increasingly rely on high-throughput phenotyping, but the resulting longitudinal data pose analytical challenges. We used canopy height data from an automated field phenotyping platform to compare several approaches to scanning for quantitative trait loci (QTLs) and performing genomic prediction in a wheat recombinant inbred line mapping population based on up to 26 sampled time points (TPs). We detected four persistent QTLs (i.e. expressed for most of the growing season), with both empirical and simulation analyses demonstrating superior statistical power of detecting such QTLs through functional mapping approaches compared with conventional individual TP analyses. In contrast, even very simple individual TP approaches (e.g. interval mapping) had superior detection power for transient QTLs (i.e. expressed during very short periods). Using spline-smoothed phenotypic data resulted in improved genomic predictive abilities (5-8% higher than individual TP prediction), while the effect of including significant QTLs in prediction models was relatively minor (<1-4% improvement). Finally, although QTL detection power and predictive ability generally increased with the number of TPs analysed, gains beyond five or 10 TPs chosen based on phenological information had little practical significance. These results will inform the development of an integrated, semi-automated analytical pipeline, which will be more broadly applicable to similar data sets in wheat and other crops.

High-density genotyping of the A.E. Watkins Collection of hexaploid landraces identifies a large molecular diversity compared to elite bread wheat.

The importance of wheat as a food crop makes it a major target for agricultural improvements. As one of the most widely grown cereal grains, together with maize and rice, wheat is the leading provider of calories in the global diet, constituting 29% of global cereal production in 2015. In the last few decades, however, yields have plateaued, suggesting that the green revolution, at least for wheat, might have run its course and that new sources of genetic variation are urgently required. The overall aim of our work was to identify novel variation that may then be used to enable the breeding process. As landraces are a potential source of such diversity, here we have characterized the A.E. Watkins Collection alongside a collection of elite accessions using two complementary high-density and high-throughput genotyping platforms. While our results show the importance of using the appropriate SNP collection to compare diverse accessions, they also show that the Watkins Collection contains a substantial amount of novel genetic diversity which has either not been captured in current breeding programmes or which has been lost through previous selection pressures. As a consequence of our analysis, we have identified a number of accessions which carry an array of novel alleles along with a number of interesting chromosome rearrangements which confirm the variable nature of the wheat genome.

The identification of new candidate genes Triticum aestivum FLOWERING LOCUS T3-B1 (TaFT3-B1) and TARGET OF EAT1 (TaTOE1-B1) controlling the short-day photoperiod response in bread wheat.

Perception of photoperiod changes enables plants to flower under optimum conditions for survival. We used doubled haploid populations of crosses among Avalon × Cadenza, Charger × Badger and Spark × Rialto and identified short-day flowering time response quantitative trait loci (QTL) on wheat chromosomes 1BS and 1BL. We used synteny between Brachypodium distachyon and wheat to identify potential candidates for both QTL. The 1BL QTL peak coincided with TaFT3-B1, a homologue of the barley gene HvFT3, the most likely candidate gene. The 1BS QTL peak coincided with homologues of Arabidopsis thaliana SENSITIVITY TO RED LIGHT REDUCED 1, WUSCHEL-like and RAP2.7, which is also known as Zea mays TARGET OF EAT1, named TaSRR1-B1, TaWUSCHELL-B1 and TaTOE1-B1, respectively. Gene expression assays suggest that TaTOE1-B1 and TaFT3-B1 are expressed more during short days. We identified four alleles of TaFT3-B1 and three alleles of TaTOE1-B1. We studied the effect of these alleles in the Watkins and GEDIFLUX diversity panels by using 936 and 431 accessions, respectively. Loss of TaFT3-B1 by deletion was associated with late flowering. Increased TaFT3-B1 copy number was associated with early flowering, suggesting that TaFT3-B1 promotes flowering. Significant association was observed in the GEDIFLUX collection for TaTOE1-B1, a putative flowering repressor.

Application of a library of near isogenic lines to understand context dependent expression of QTL for grain yield and adaptive traits in bread wheat.

BACKGROUND: Previous quantitative trait loci (QTLs) studies using the Avalon × Cadenza doubled haploid (DH) population identified eleven QTLs determining plant height, heading date and grain yield. The objectives of this study were: (i) to provide insight into the effects of these QTLs using reciprocal multiple near isogenic lines (NILs) with each pair of alleles compared in both parental backgrounds (Avalon or Cadenza), (ii) quantifying epistasis by looking at the background effects and (iii) predict favourable allelic combinations to develop superior genotypes adapted to a target environment. RESULTS: To this aim, a library of 553 BC2 NILs and their recurrent parents were tested over two growing seasons (2012/2013 and 2013/2014). The results obtained in the present study validated the plant height, heading date and grain yield QTLs previously identified. Epistatic interactions were detected for the 6B QTL for plant height and heading date, 3A QTL for heading date and grain yield and 2A QTL for grain yield. CONCLUSION: The marker assisted backcrossing strategy used provided an efficient method of resolving QTL for key agronomic traits in wheat as Mendelian factors determining possible epistatic interactions. The study shows that these QTLs are amenable to marker assisted selection, fine mapping, future positional cloning, and physiological trait dissection.

Delimitation of the Earliness per se D1 (Eps-D1) flowering gene to a subtelomeric chromosomal deletion in bread wheat (Triticum aestivum).

Earliness per se (Eps) genes account for the variation in flowering time when vernalization and photoperiod requirements are satisfied. Genomics and bioinformatics approaches were used to describe allelic variation for 40 Triticum aestivum genes predicted, by synteny with Brachypodium distachyon, to be in the 1DL Eps region. Re-sequencing 1DL genes revealed that varieties carrying early heading alleles at this locus, Spark and Cadenza, carry a subtelomeric deletion including several genes. The equivalent region in Rialto and Avalon is intact. A bimodal distribution in the segregating Spark X Rialto single seed descent (SSD) populations enabled the 1DL QTL to be defined as a discrete Mendelian factor, which we named Eps-D1. Near isogenic lines (NILs) and NIL derived key recombinants between markers flanking Eps-D1 suggest that the 1DL deletion contains the gene(s) underlying Eps-D1. The deletion spans the equivalent of the Triticum monoccocum Eps-A (m) 1 locus, and hence includes MODIFIER OF TRANSCRIPTION 1 (MOT1) and FTSH PROTEASE 4 (FTSH4), the candidates for Eps-A (m) 1. The deletion also contains T. aestivum EARLY FLOWERING 3-D1 (TaELF3-D1) a homologue of the Arabidopsis thaliana circadian clock gene EARLY FLOWERING 3. Eps-D1 is possibly a homologue of Eps-B1 on chromosome 1BL. NILs carrying the Eps-D1 deletion have significantly reduced total TaELF3 expression and altered TaGIGANTEA (TaGI) expression compared with wild type. Altered TaGI expression is consistent with an ELF3 mutant, hence we propose TaELF3-D1 as the more likely candidate for Eps-D1. This is the first direct fine mapping of Eps effect in bread wheat.

Phenotyping pipeline reveals major seedling root growth QTL in hexaploid wheat.

Seedling root traits of wheat (Triticum aestivum L.) have been shown to be important for efficient establishment and linked to mature plant traits such as height and yield. A root phenotyping pipeline, consisting of a germination paper-based screen combined with image segmentation and analysis software, was developed and used to characterize seedling traits in 94 doubled haploid progeny derived from a cross between the winter wheat cultivars Rialto and Savannah. Field experiments were conducted to measure mature plant height, grain yield, and nitrogen (N) uptake in three sites over 2 years. In total, 29 quantitative trait loci (QTLs) for seedling root traits were identified. Two QTLs for grain yield and N uptake co-localize with root QTLs on chromosomes 2B and 7D, respectively. Of the 29 root QTLs identified, 11 were found to co-localize on 6D, with four of these achieving highly significant logarithm of odds scores (>20). These results suggest the presence of a major-effect gene regulating seedling root vigour/growth on chromosome 6D.

Molecular genetic basis of pod corn (Tunicate maize).

Pod corn is a classic morphological mutant of maize in which the mature kernels of the cob are covered by glumes, in contrast to generally grown maize varieties in which kernels are naked. Pod corn, known since pre-Columbian times, is the result of a dominant gain-of-function mutation at the Tunicate (Tu) locus. Some classic articles of 20th century maize genetics reported that the mutant Tu locus is complex, but molecular details remained elusive. Here, we show that pod corn is caused by a cis-regulatory mutation and duplication of the ZMM19 MADS-box gene. Although the WT locus contains a single-copy gene that is expressed in vegetative organs only, mutation and duplication of ZMM19 in Tu lead to ectopic expression of the gene in the inflorescences, thus conferring vegetative traits to reproductive organs.

Establishing the A. E. Watkins landrace cultivar collection as a resource for systematic gene discovery in bread wheat.

KEY MESSAGE: A high level of genetic diversity was found in the A. E. Watkins bread wheat landrace collection. Genotypic information was used to determine the population structure and to develop germplasm resources. In the 1930s A. E. Watkins acquired landrace cultivars of bread wheat (Triticum aestivum L.) from official channels of the board of Trade in London, many of which originated from local markets in 32 countries. The geographic distribution of the 826 landrace cultivars of the current collection, here called the Watkins collection, covers many Asian and European countries and some from Africa. The cultivars were genotyped with 41 microsatellite markers in order to investigate the genetic diversity and population structure of the collection. A high level of genetic diversity was found, higher than in a collection of modern European winter bread wheat varieties from 1945 to 2000. Furthermore, although weak, the population structure of the Watkins collection reveals nine ancestral geographical groupings. An exchange of genetic material between ancestral groups before commercial wheat-breeding started would be a possible explanation for this. The increased knowledge regarding the diversity of the Watkins collection was used to develop resources for wheat research and breeding, one of them a core set, which captures the majority of the genetic diversity detected. The understanding of genetic diversity and population structure together with the availability of breeding resources should help to accelerate the detection of new alleles in the Watkins collection.

The wheat powdery mildew resistance gene Pm4 also confers resistance to wheat blast.

Wheat blast, caused by the fungus Magnaporthe oryzae, threatens global cereal production since its emergence in Brazil in 1985 and recently spread to Bangladesh and Zambia. Here we demonstrate that the AVR-Rmg8 effector, common in wheat-infecting isolates, is recognized by the gene Pm4, previously shown to confer resistance to specific races of Blumeria graminis f. sp. tritici, the cause of powdery mildew of wheat. We show that Pm4 alleles differ in their recognition of different AVR-Rmg8 alleles, and some confer resistance only in seedling leaves but not spikes, making it important to select for those alleles that function in both tissues. This study has identified a gene recognizing an important virulence factor present in wheat blast isolates in Bangladesh and Zambia and represents an important first step towards developing durably resistant wheat cultivars for these regions.

Nitrogen uptake and remobilization from pre- and post-anthesis stages contribute towards grain yield and grain protein concentration in wheat grown in limited nitrogen conditions.

Background: In wheat, nitrogen (N) remobilization from vegetative tissues to developing grains largely depends on genetic and environmental factors. The evaluation of genetic potential of crops under limited resource inputs such as limited N supply would provide an opportunity to identify N-efficient lines with improved N utilisation efficiency and yield potential. We assessed the genetic variation in wheat recombinant inbred lines (RILs) for uptake, partitioning, and remobilization of N towards grain, its association with grain protein concentration (GPC) and grain yield. Methods: We used the nested association mapping (NAM) population (195 lines) derived by crossing Paragon (P) with CIMMYT core germplasm (P × Cim), Baj (P × Baj), Watkins (P × Wat), and Wyalkatchem (P × Wya). These lines were evaluated in the field for two seasons under limited N supply. The plant sampling was done at anthesis and physiological maturity stages. Various physiological traits were recorded and total N uptake and other N related indices were calculated. The grain protein deviation (GPD) was calculated from the regression of grain yield on GPC. These lines were grouped into different clusters by hierarchical cluster analysis based on grain yield and N-remobilization efficiency (NRE). Results: The genetic variation in accumulation of biomass at both pre- and post-anthesis stages were correlated with grain-yield. The NRE significantly correlated with aboveground N uptake at anthesis (AGNa) and grain yield but negatively associated with AGN at post-anthesis (AGNpa) suggesting higher N uptake till anthesis favours high N remobilization during grain filling. Hierarchical cluster analysis of these RILs based on NRE and yield resulted in four clusters, efficient (31), moderately efficient (59), moderately inefficient (58), and inefficient (47). In the N-efficient lines, AGNa contributed to 77% of total N accumulated in grains, while it was 63% in N-inefficient lines. Several N-efficient lines also exhibited positive grain protein deviation (GPD), combining high grain yield and GPC. Among crosses, the P × Cim were superior and N-efficient, while P × Wya responded poorly to low N input. Conclusions: We propose that traits favouring pre- or post-anthesis biomass accumulation and pre-anthesis N uptake may be targeted for breeding to improve grain-yield under limited N. The lines with positive GPD, a first report of genotype-dependent GPD associated with both AGNpa and AGNa in wheat, may be used as varieties or genetic resources to improve grain yield with high GPC for sustainable development under limited N conditions. Supplementary Information: The online version contains supplementary material available at 10.1186/s43170-023-00153-7.

Genetic control of grain amino acid composition in a UK soft wheat mapping population.

Wheat (Triticum aestivum L.) is a major source of nutrients for populations across the globe, but the amino acid composition of wheat grain does not provide optimal nutrition. The nutritional value of wheat grain is limited by low concentrations of lysine (the most limiting essential amino acid) and high concentrations of free asparagine (precursor to the processing contaminant acrylamide). There are currently few available solutions for asparagine reduction and lysine biofortification through breeding. In this study, we investigated the genetic architecture controlling grain free amino acid composition and its relationship to other traits in a Robigus × Claire doubled haploid population. Multivariate analysis of amino acids and other traits showed that the two groups are largely independent of one another, with the largest effect on amino acids being from the environment. Linkage analysis of the population allowed identification of quantitative trait loci (QTL) controlling free amino acids and other traits, and this was compared against genomic prediction methods. Following identification of a QTL controlling free lysine content, wheat pangenome resources facilitated analysis of candidate genes in this region of the genome. These findings can be used to select appropriate strategies for lysine biofortification and free asparagine reduction in wheat breeding programs.

A wheat kinase and immune receptor form host-specificity barriers against the blast fungus.

Since emerging in Brazil in 1985, wheat blast has spread throughout South America and recently appeared in Bangladesh and Zambia. Here we show that two wheat resistance genes, Rwt3 and Rwt4, acting as host-specificity barriers against non-Triticum blast pathotypes encode a nucleotide-binding leucine-rich repeat immune receptor and a tandem kinase, respectively. Molecular isolation of these genes will enable study of the molecular interaction between pathogen effector and host resistance genes.

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.

Opinion Exploiting genomics to improve the benefits of wheat: Prospects and limitations.

Image 1.

An Integrated Linkage Map of Three Recombinant Inbred Populations of Pea (Pisum sativum L.).

Biparental recombinant inbred line (RIL) populations are sets of genetically stable lines and have a simple population structure that facilitates the dissection of the genetics of interesting traits. On the other hand, populations derived from multiparent intercrosses combine both greater diversity and higher numbers of recombination events than RILs. Here, we describe a simple population structure: a three-way recombinant inbred population combination. This structure was easy to produce and was a compromise between biparental and multiparent populations. We show that this structure had advantages when analyzing cultivar crosses, and could achieve a mapping resolution of a few genes.

Interactions between two QTLs for time to anthesis on spike development and fertility in wheat.

Earliness per se (Eps) genes are reported to be important in fine-tuning flowering time in wheat independently of photoperiod (Ppd) and vernalisation (Vrn). Unlike Ppd and Vrn genes, Eps have relatively small effects and their physiological effect along with chromosomal position are not well defined. We evaluated eight lines derived from crossing two vernalisation insensitive lines, Paragon and Baj (late and early flowering respectively), to study the detailed effects of two newly identified QTLs, Eps-7D and Eps-2B and their interactions under field conditions. The effect of both QTLs was minor and was affected by the allelic status of the other. While the magnitude of effect of these QTLs on anthesis was similar, they are associated with very different profiles of pre-anthesis development which also depends on their interaction. Eps-7D affected both duration before and after terminal spikelet while not affecting final leaf number (FLN) so Eps-7D-early had a faster rate of leaf appearance. Eps-2B acted more specifically in the early reproductive phase and slightly altered FLN without affecting the leaf appearance rate. Both QTLs affected the spike fertility by altering the rate of floret development and mortality. The effect of Eps-2B was very small but consistent in that -late allele tended to produce more fertile florets.

The role of gene flow and chromosomal instability in shaping the bread wheat genome.

Bread wheat (Triticum aestivum) is one of the world's most important crops; however, a low level of genetic diversity within commercial breeding accessions can significantly limit breeding potential. In contrast, wheat relatives exhibit considerable genetic variation and so potentially provide a valuable source of novel alleles for use in breeding new cultivars. Historically, gene flow between wheat and its relatives may have contributed novel alleles to the bread wheat pangenome. To assess the contribution made by wheat relatives to genetic diversity in bread wheat, we used markers based on single nucleotide polymorphisms to compare bread wheat accessions, created in the past 150 years, with 45 related species. We show that many bread wheat accessions share near-identical haplotype blocks with close relatives of wheat's diploid and tetraploid progenitors, while some show evidence of introgressions from more distant species and structural variation between accessions. Hence, introgressions and chromosomal rearrangements appear to have made a major contribution to genetic diversity in cultivar collections. As gene flow from relatives to bread wheat is an ongoing process, we assess the impact that introgressions might have on future breeding strategies.

Genetic variation in wheat grain quality is associated with differences in the galactolipid content of flour and the gas bubble properties of dough liquor.

Lipids affect the quality of wheat flour for breadmaking. One possible mechanism is stabilization of the gas cells which are formed during dough mixing and expanded during fermentation, leading to a greater loaf volume and evenness of texture. We therefore compared the lipidomic profiles of flour and dough liquor fractions (which contain surface-active components present at the gas bubble interface) from two sets of wheat lines differing in allelic variation at a QTL for loaf volume. Analyses of fractions from three field trials showed consistent increases in the contents of galactolipids (monogalactosyl diglyceride and digalactosyl diglyceride) in flour and dough liquor of the lines with the increasing (good quality) allele. Biophysical analysis showed that this was associated with greater elasticity of the dough liquor fraction. This is consistent with published studies reporting a relationship between galactolipids and breadmaking quality and suggests a mechanism of action for the QTL.