Strong linkage disequilibrium near the selected Yr17 resistance gene in a wheat experimental population

Analysis of historical selection in winter wheat

KEY MESSAGE

Modeling of the distribution of allele frequency over year of variety release identifies major loci involved in historical breeding of winter wheat. Winter wheat is a major crop with a rich selection history in the modern era of crop breeding. Genetic gains across economically important traits like yield have been well characterized and are the major force driving its production. Winter wheat is also an excellent model for analyzing historical genetic selection. As a proof of concept, we analyze two major collections of winter wheat varieties that were bred in Western Europe from 1916 to 2010, namely the Triticeae Genome (TG) and WAGTAIL panels, which include 333 and 403 varieties, respectively. We develop and apply a selection mapping approach, Regression of Alleles on Years (RALLY), in these panels, as well as in simulated populations. RALLY maps loci under sustained historical selection by using a simple logistic model to regress allele counts on years of variety release. To control for drift-induced allele frequency change, we develop a hybrid approach of genomic control and delta control. Within the TG panel, we identify 22 significant RALLY quantitative selection loci (QSLs) and estimate the local heritabilities for 12 traits across these QSLs. By correlating predicted marker effects with RALLY regression estimates, we show that alleles whose frequencies have increased over time are heavily biased toward conferring positive yield effect, but negative effects in flowering time, lodging, plant height and grain protein content. Altogether, our results (1) demonstrate the use of RALLY to identify selected genomic regions while controlling for drift, and (2) reveal key patterns in the historical selection in winter wheat and guide its future breeding.

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.

A large-scale pedigree resource of wheat reveals evidence for adaptation and selection by breeders

Information on crop pedigrees can be used to help maximise genetic gain in crop breeding and allow efficient management of genetic resources. We present a pedigree resource of 2,657 wheat (Triticum aestivum L.) genotypes originating from 38 countries, representing more than a century of breeding and variety development. Visualisation of the pedigree enables illustration of the key developments in United Kingdom wheat breeding, highlights the wide genetic background of the UK wheat gene pool, and facilitates tracing the origin of beneficial alleles. A relatively high correlation between pedigree- and marker-based kinship coefficients was found, which validated the pedigree and enabled identification of errors in the pedigree or marker data. Using simulations with a combination of pedigree and genotype data, we found evidence for significant effects of selection by breeders. Within crosses, genotypes are often more closely related than expected by simulations to one of the parents, which indicates selection for favourable alleles during the breeding process. Selection across the pedigree was demonstrated on a subset of the pedigree in which 110 genotyped varieties released before the year 2000 were used to simulate the distribution of marker alleles of 45 genotyped varieties released after the year 2000, in the absence of selection. Allelic diversity in the 45 varieties was found to deviate significantly from the simulated distributions at a number of loci, indicating regions under selection over this period. The identification of one of these regions as coinciding with a strong yield component quantitative trait locus (QTL) highlights both the potential of the remaining loci as wheat breeding targets for further investigation, as well as the utility of this pedigree-based methodology to identify important breeding targets in other crops. Further evidence for selection was found as greater linkage disequilibrium (LD) for observed versus simulated genotypes within all chromosomes. This difference was greater at shorter genetic distances, indicating that breeder selections have conserved beneficial linkage blocks. Collectively, this work highlights the benefits of generating detailed pedigree resources for crop species. The wheat pedigree database developed here represents a valuable community resource and will be updated as new varieties are released at https://www.niab.com/pages/id/501/UK_Wheat_varieties_Pedigree.

Breeding activities undertaken in the world’s most important crop species have resulted in large increases in yield potential over the last century. Bread wheat is a key crop for both human and animal nutrition worldwide. To help inform future breeding and research activities, we have developed a pedigree resource of over 2,600 bread wheat accessions, originating from 38 countries and representing more than a century of breeding and variety development. Pedigree-based relationships between lines are largely confirmed by genetic marker data. By combining the genetic and pedigree data sets, we are able to identify genetic signatures of selection across the pedigree, identifying genomic regions selected for via modern breeding activities. The resource developed here will serve as an evolving platform to inform and manage wheat genetic diversity in breeding programmes around the world and highlights the utility of developing and exploiting similar resources in other crop species.

Characterization of a major QTL for adult plant resistance to stripe rust in US soft red winter wheat

Stripe rust, caused by Puccinia striiformis f. sp. tritici, is an important disease of soft red winter wheat in the eastern region of the USA. Pioneer 26R61 has provided effective resistance to stripe rust for 10 years. To elucidate the genetic basis of the resistance, a mapping population of 178 recombinant inbred lines (RILs) was developed using single-seed descent from a cross between Pioneer 26R61 and the susceptible cultivar AGS 2000. A genetic map with 895 markers covering all 21 chromosomes was used for QTL analysis. One major QTL was detected, explaining up to 56.0% of the mean phenotypic variation, flanked by markers Xbarc124 and Xgwm359, and assigned to the distal 22% of the short arm of wheat chromosome 2A. Evidence showed that it was different from Yr17 derived from Ae. ventricosa, the only formally named Yr gene in 2AS, and the QTL was temporarily designated as YrR61. In addition, a minor QTL, QYr.uga-6AS, probably conditioned high-temperature adult plant resistance. The QTL explained 6-7% of the trait variation. Preliminary test of the flanking markers for YrR61, in two cultivars and two promising breeding lines with Pioneer 26R61 in their pedigree, indicated that YrR61 was present in these cultivars and lines, and these markers could therefore be used in marker-assisted selection.

Dynamic management of crop diversity: From an experimental approach to on-farm conservation

In agricultural systems, biodiversity includes diversity within species and among species and provides many benefits for production, resilience and conservation. This article addresses the effects of a strategy of in situ conservation called dynamic management (DM) on population evolution, adaptation and diversity. Two French DM initiatives are considered, the first one corresponding to an experimental context, the second to an on-farm management. Results from a study over 26 years of experimental DM of bread wheat (Triticum aestivum L.) are first presented, including the evolution of agronomic traits and genetic diversity at neutral and fitness related loci. While this experiment greatly increased scientific knowledge of the effects of natural selection on cultivated populations, it also showed that population conservation cannot rely only on a network of experimental stations. In collaboration with a farmers' network in France, researchers have begun studying the effects of on-farm DM (conservation and selection) on diversity and adaptation. Results from these studies show that on-farm DM is a key element for the long-term conservation and use of agricultural biodiversity. This method of in situ conservation deserves more attention in industrialised countries.

Evolution of flowering time in experimental wheat populations: a comprehensive approach to detect genetic signatures of natural selection

In annual plant species, flowering time is a major adaptive trait that synchronizes the initiation of reproduction with favorable environmental conditions. Here, we aimed at studying the evolution of flowering time in three experimental populations of bread wheat, grown in contrasting environments (Northern to Southern France) for 12 generations. By comparing the distribution of phenotypic and presumably neutral variation, we first showed that flowering time responded to selection during the 12 generations of the experiment. To get insight into the genetic architecture of that trait, we then tested whether the distribution of genetic polymorphisms at six candidate genes, presumably involved in the trait expression, departed from neutral expectation. To that end, we focused on the temporal variation during the course of the experiment, and on the spatial differentiation at the end of the experiment, using previously published methods adapted to our experimental design. Only those genes that were strongly associated with flowering time variation were detected as responding to selection. For genes that had low-to-moderate phenotypic effects, or when there was interaction across different genes, we did not find evidence of selection using methods based on the distribution of temporal or spatial variation. In such cases, it might be more informative to consider multilocus and multiallelic combinations across genes, which could be the targets of selection.

Association mapping of agronomic traits on chromosome 2A of wheat

Association mapping is a method to test the association between molecular markers and quantitative trait loci (QTL) based on linkage disequilibrium (LD). In this study, the collection of 108 wheat germplasm accessions form China were evaluated for their plant heights, spike length, spikelets per spike, grains per spike, thousand kernel weight and spikelets density in 3 years at three locations. And they were genotyped with 85 SSR markers and 40 EST-SSR markers. The population structure was inferred on the basis of unlinked 48 SSR markers and 40 EST-SSR markers. The extent of LD on chromosome 2A was 2.3 cM. Association of 37 SSR loci on chromosomes 2A with six agronomic traits was analysed with a mixed linear model. A total of 14 SSR loci were significantly associated with agronomic traits. Some of the associated markers were located in the QTL region detected in previous linkage mapping analysis. Our results demonstrated that association mapping can enhance QTL information and achieves higher resolution with short LD extent.

Insight into the genetic bases of climatic adaptation in experimentally evolving wheat populations

Experimental populations evolving under natural selection represent an interesting tool to study genetic bases of adaptation. Evolution of genes possibly involved in adaptive response can be followed together with the corresponding phenotypic traits. Using experimental populations of hexaploid wheat, we studied the evolution of flowering time, a major adaptive trait that synchronizes the initiation of reproduction and the occurrence of favourable environmental conditions. During 12 generations, three populations were grown in contrasted environments (Vervins North France, Le Moulon near Paris, Toulouse South France) under the influence of natural selection, drift, mutation and recombination. Evolution of diversity at the major gene VRN-1 involved in wheat vernalization response has been analysed jointly with earliness estimated in controlled conditions. Whatever the population, rapid phenotypic changes as well as parallel genotypic variations were observed in the first seven generations, probably as the result of selection acting on this major gene which explains 80% of the trait variation overall. Different allelic combinations at physically unlinked copies of VRN-1 located on distinct genomes (A, B and D) were selected between populations. As theoretically expected, due to population differentiation, a high level of genetic diversity was maintained overall in generation 12. Surprisingly, in two populations out of three, the emergence of new alleles by mutation or migration, coupled with temporal variable selection or frequency-dependent selection, allowed to maintain within-population diversity despite local genetic drift and natural selection. This result may plead for an evolutionary approach of wheat genetic resource conservation.

Soft selective sweep near a gene that increases plant height in wheat

Strong selection within a given population locally reduces genetic variability not only in the selected gene itself but also in neighbouring loci. This so-called hitch-hiking effect is related to the initial linkage disequilibrium between markers and the selected gene, and depends mainly on the number of copies of the beneficial allele at the start of the selection phase. Contrary to the classical case, in which selection acts on a single, newly arisen beneficial mutation, we considered selection from standing variation (soft selective sweeps) on a gene (Rht-B1) with a major effect on plant height, a selected trait in an experimental wheat population grown for 17 generations, and we documented the evolution of gene diversity and linkage disequilibrium near this gene. As expected, Rht-B1 was found to be under strong selection (s = 0.15) and its variation in frequency accounted for 15% of the total trait evolution. This led to a smaller genetic effective population size at Rht-B1 (N(eg) = 18) compared to the whole genome estimation (N(eg) = 167). When compared with expectations under genetic drift only, no significant decrease in gene diversity was found at the closest loci. We computed expected di-locus frequencies for any linked marker-Rht-B1 pair due to hitch-hiking effects. We found that hitch-hiking was expected to affect the two most closely linked loci, but expected reduction in gene diversity was not greater than that due to genetic drift, which was consistent with the observations. Such limited effect was attributed to the low level of linkage disequilibrium (0.16) estimated after parental intercrosses, together with a relatively high initial frequency of the gene. This situation is favourable to candidate gene approaches where small linkage disequilibrium around selected genes is expected.