Investigating the genetic control of plant development in spring barley under speed breeding conditions

KEY MESSAGE

This study found that the genes, PPD-H1 and ELF3, control the acceleration of plant development under speed breeding, with important implications for optimizing the delivery of climate-resilient crops. Speed breeding is a tool to accelerate breeding and research programmes. Despite its success and growing popularity with breeders, the genetic basis of plant development under speed breeding remains unknown. This study explored the developmental advancements of barley genotypes under different photoperiod regimes. A subset of the HEB-25 Nested Association Mapping population was evaluated for days to heading and maturity under two contrasting photoperiod conditions: (1) Speed breeding (SB) consisting of 22 h of light and 2 h of darkness, and (2) normal breeding (NB) consisting of 16 h of light and 8 h of darkness. GWAS revealed that developmental responses under both conditions were largely controlled by two loci: PPDH-1 and ELF3. Allelic variants at these genes determine whether plants display early flowering and maturity under both conditions. At key QTL regions, domesticated alleles were associated with late flowering and maturity in NB and early flowering and maturity in SB, whereas wild alleles were associated with early flowering under both conditions. We hypothesize that this is related to the dark-dependent repression of PPD-H1 by ELF3 which might be more prominent in NB conditions. Furthermore, by comparing development under two photoperiod regimes, we derived an estimate of plasticity for the two traits. Interestingly, plasticity in development was largely attributed to allelic variation at ELF3. Our results have important implications for our understanding and optimization of speed breeding protocols particularly for introgression breeding and the design of breeding programmes to support the delivery of climate-resilient crops.

Phenomic and genomic prediction of yield on multiple locations in winter wheat

Genomic selection has recently become an established part of breeding strategies in cereals. However, a limitation of linear genomic prediction models for complex traits such as yield is that these are unable to accommodate Genotype by Environment effects, which are commonly observed over trials on multiple locations. In this study, we investigated how this environmental variation can be captured by the collection of a large number of phenomic markers using high-throughput field phenotyping and whether it can increase GS prediction accuracy. For this purpose, 44 winter wheat (Triticum aestivum L.) elite populations, comprising 2,994 lines, were grown on two sites over 2 years, to approximate the size of trials in a practical breeding programme. At various growth stages, remote sensing data from multi- and hyperspectral cameras, as well as traditional ground-based visual crop assessment scores, were collected with approximately 100 different data variables collected per plot. The predictive power for grain yield was tested for the various data types, with or without genome-wide marker data sets. Models using phenomic traits alone had a greater predictive value (R² = 0.39-0.47) than genomic data (approximately R² = 0.1). The average improvement in predictive power by combining trait and marker data was 6%-12% over the best phenomic-only model, and performed best when data from one full location was used to predict the yield on an entire second location. The results suggest that genetic gain in breeding programmes can be increased by utilisation of large numbers of phenotypic variables using remote sensing in field trials, although at what stage of the breeding cycle phenomic selection could be most profitably applied remains to be answered.

Genome-wide association mapping of Hagberg falling number, protein content, test weight, and grain yield in U.K. wheat

Association mapping using crop cultivars allows identification of genetic loci of direct relevance to breeding. Here, 150 U.K. wheat (Triticum aestivum L.) cultivars genotyped with 23,288 single nucleotide polymorphisms (SNPs) were used for genome-wide association studies (GWAS) using historical phenotypic data for grain protein content, Hagberg falling number (HFN), test weight, and grain yield. Power calculations indicated experimental design would enable detection of quantitative trait loci (QTL) explaining ≥20% of the variation (PVE) at a relatively high power of >80%, falling to 40% for detection of a SNP with an R² ≥ .5 with the same QTL. Genome-wide association studies identified marker-trait associations for all four traits. For HFN (h ² = .89), six QTL were identified, including a major locus on chromosome 7B explaining 49% PVE and reducing HFN by 44 s. For protein content (h ² = 0.86), 10 QTL were found on chromosomes 1A, 2A, 2B, 3A, 3B, and 6B, together explaining 48.9% PVE. For test weight, five QTL were identified (one on 1B and four on 3B; 26.3% PVE). Finally, 14 loci were identified for grain yield (h ² = 0.95) on eight chromosomes (1A, 2A, 2B, 2D, 3A, 5B, 6A, 6B; 68.1% PVE), of which five were located within 16 Mbp of genetic regions previously identified as under breeder selection in European wheat. Our study demonstrates the utility of exploiting historical crop datasets, identifying genomic targets for independent validation, and ultimately for wheat genetic improvement.

The evolving battle between yellow rust and wheat: implications for global food security

Wheat (Triticum aestivum L.) is a global commodity, and its production is a key component underpinning worldwide food security. Yellow rust, also known as stripe rust, is a wheat disease caused by the fungus Puccinia striiformis Westend f. sp. tritici (Pst), and results in yield losses in most wheat growing areas. Recently, the rapid global spread of genetically diverse sexually derived Pst races, which have now largely replaced the previous clonally propagated slowly evolving endemic populations, has resulted in further challenges for the protection of global wheat yields. However, advances in the application of genomics approaches, in both the host and pathogen, combined with classical genetic approaches, pathogen and disease monitoring, provide resources to help increase the rate of genetic gain for yellow rust resistance via wheat breeding while reducing the carbon footprint of the crop. Here we review key elements in the evolving battle between the pathogen and host, with a focus on solutions to help protect future wheat production from this globally important disease.

Trends of genetic changes uncovered by Env- and Eigen-GWAS in wheat and barley

Variety age and population structure detect novel QTL for yield and adaptation in wheat and barley without the need to phenotype. The process of crop breeding over the last century has delivered new varieties with increased genetic gains, resulting in higher crop performance and yield. However, in many cases, the alleles and genomic regions underpinning this success remain unknown. This is partly due to the difficulty of generating sufficient phenotypic data on large numbers of historical varieties to enable such analyses. Here we demonstrate the ability to circumvent such bottlenecks by identifying genomic regions selected over 100 years of crop breeding using age of a variety as a surrogate for yield. Rather than collecting phenotype data, we deployed 'environmental genome-wide association scans' (EnvGWAS) based on variety age in two of the world's most important crops, wheat and barley, and detected strong signals of selection across both genomes. EnvGWAS identified 16 genomic regions in barley and 10 in wheat with contrasting patterns between spring and winter types of the two crops. To further examine changes in genome structure, we used the genomic relationship matrix of the genotypic data to derive eigenvectors for analysis in EigenGWAS. This detected seven major chromosomal introgressions that contributed to adaptation in wheat. EigenGWAS and EnvGWAS based on variety age avoid costly phenotyping and facilitate the identification of genomic tracts that have been under selection during breeding. Our results demonstrate the potential of using historical cultivar collections coupled with genomic data to identify chromosomal regions under selection and may help guide future plant breeding strategies to maximise the rate of genetic gain and adaptation.

Wheat genetic loci conferring resistance to stripe rust in the face of genetically diverse races of the fungus Puccinia striiformis f. sp. tritici

KEY MESSAGE

Analysis of a wheat multi-founder population identified 14 yellow rust resistance QTL. For three of the four most significant QTL, haplotype analysis indicated resistance alleles were rare in European wheat. Stripe rust, or yellow rust (YR), is a major fungal disease of wheat (Triticum aestivum) caused by Puccinia striiformis Westend f. sp. tritici (Pst). Since 2011, the historically clonal European Pst races have been superseded by the rapid incursion of genetically diverse lineages, reducing the resistance of varieties previously showing durable resistance. Identification of sources of genetic resistance to such races is a high priority for wheat breeding. Here we use a wheat eight-founder multi-parent population genotyped with a 90,000 feature single nucleotide polymorphism array to genetically map YR resistance to such new Pst races. Genetic analysis of five field trials at three UK sites identified 14 quantitative trait loci (QTL) conferring resistance. Of these, four highly significant loci were consistently identified across all test environments, located on chromosomes 1A (QYr.niab-1A.1), 2A (QYr.niab-2A.1), 2B (QYr.niab-2B.1) and 2D (QYr.niab-2D.1), together explaining ~ 50% of the phenotypic variation. Analysis of these four QTL in two-way and three-way combinations showed combinations conferred greater resistance than single QTL, and genetic markers were developed that distinguished resistant and susceptible alleles. Haplotype analysis in a collection of wheat varieties found that the haplotypes associated with YR resistance at three of these four major loci were rare (≤ 7%) in European wheat, highlighting their potential utility for future targeted improvement of disease resistance. Notably, the physical interval for QTL QYr.niab-2B.1 contained five nucleotide-binding leucine-rich repeat candidate genes with integrated BED domains, of which two corresponded to the cloned resistance genes Yr7 and Yr5/YrSp.

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.

Understanding the classics: the unifying concepts of transgressive segregation, inbreeding depression and heterosis and their central relevance for crop breeding

Transgressive segregation and heterosis are the reasons that plant breeding works. Molecular explanations for both phenomena have been suggested and play a contributing role. However, it is often overlooked by molecular genetic researchers that transgressive segregation and heterosis are most simply explained by dispersion of favorable alleles. Therefore, advances in molecular biology will deliver the most impact on plant breeding when integrated with sources of heritable trait variation - and this will be best achieved within a quantitative genetics framework. An example of the power of quantitative approaches is the implementation of genomic selection, which has recently revolutionized animal breeding. Genomic selection is now being applied to both hybrid and inbred crops and is likely to be the major source of improvement in plant breeding practice over the next decade. Breeders' ability to efficiently apply genomic selection methodologies is due to recent technology advances in genotyping and sequencing. Furthermore, targeted integration of additional molecular data (such as gene expression, gene copy number and methylation status) into genomic prediction models may increase their performance. In this review, we discuss and contextualize a suite of established quantitative genetics themes relating to hybrid vigour, transgressive segregation and their central relevance to plant breeding, with the aim of informing crop researchers outside of the quantitative genetics discipline of their relevance and importance to crop improvement. Better understanding between molecular and quantitative disciplines will increase the potential for further improvements in plant breeding methodologies and so help underpin future food security.

Multi-parent populations in crops: a toolbox integrating genomics and genetic mapping with breeding

Crop populations derived from experimental crosses enable the genetic dissection of complex traits and support modern plant breeding. Among these, multi-parent populations now play a central role. By mixing and recombining the genomes of multiple founders, multi-parent populations combine many commonly sought beneficial properties of genetic mapping populations. For example, they have high power and resolution for mapping quantitative trait loci, high genetic diversity and minimal population structure. Many multi-parent populations have been constructed in crop species, and their inbred germplasm and associated phenotypic and genotypic data serve as enduring resources. Their utility has grown from being a tool for mapping quantitative trait loci to a means of providing germplasm for breeding programmes. Genomics approaches, including de novo genome assemblies and gene annotations for the population founders, have allowed the imputation of rich sequence information into the descendent population, expanding the breadth of research and breeding applications of multi-parent populations. Here, we report recent successes from crop multi-parent populations in crops. We also propose an ideal genotypic, phenotypic and germplasm 'package' that multi-parent populations should feature to optimise their use as powerful community resources for crop research, development and breeding.

Reference Genome Anchoring of High-Density Markers for Association Mapping and Genomic Prediction in European Winter Wheat

In this study, we anchored genotyping-by-sequencing data to the International Wheat Genome Sequencing Consortium Reference Sequence v1.0 assembly to generate over 40,000 high quality single nucleotide polymorphism markers on a panel of 376 elite European winter wheat varieties released between 1946 and 2007. We compared association mapping and genomic prediction accuracy for a range of productivity traits with previous results based on lower density dominant DArT markers. The results demonstrate that the availability of RefSeq v1.0 supports higher precision trait mapping and provides the density of markers required to obtain accurate predictions of traits controlled by multiple small effect loci, including grain yield.

Genomic index selection provides a pragmatic framework for setting and refining multi-objective breeding targets in Miscanthus

BACKGROUND

Miscanthus has potential as a biomass crop but the development of varieties that are consistently superior to the natural hybrid M. × giganteus has been challenging, presumably because of strong G × E interactions and poor knowledge of the complex genetic architectures of traits underlying biomass productivity and climatic adaptation. While linkage and association mapping studies are starting to generate long lists of candidate regions and even individual genes, it seems unlikely that this information can be translated into effective marker-assisted selection for the needs of breeding programmes. Genomic selection has emerged as a viable alternative, and prediction accuracies are moderate across a range of phenological and morphometric traits in Miscanthus, though relatively low for biomass yield per se.

METHODS

We have previously proposed a combination of index selection and genomic prediction as a way of overcoming the limitations imposed by the inherent complexity of biomass yield. Here we extend this approach and illustrate its potential to achieve multiple breeding targets simultaneously, in the absence of a priori knowledge about their relative economic importance, while also monitoring correlated selection responses for non-target traits. We evaluate two hypothetical scenarios of increasing biomass yield by 20 % within a single round of selection. In the first scenario, this is achieved in combination with delaying flowering by 44 d (roughly 20 %), whereas, in the second, increased yield is targeted jointly with reduced lignin (-5 %) and increased cellulose (+5 %) content, relative to current average levels in the breeding population.

KEY RESULTS

In both scenarios, the objectives were achieved efficiently (selection intensities corresponding to keeping the best 20 and 4 % of genotypes, respectively). However, the outcomes were strikingly different in terms of correlated responses, and the relative economic values (i.e. value per unit of change in each trait compared with that for biomass yield) of secondary traits included in selection indices varied considerably.

CONCLUSIONS

Although these calculations rely on multiple assumptions, they highlight the need to evaluate breeding objectives and explicitly consider correlated responses in silico, prior to committing extensive resources. The proposed approach is broadly applicable for this purpose and can readily incorporate high-throughput phenotyping data as part of integrated breeding platforms.

The effects of training population design on genomic prediction accuracy in wheat

Genomic selection offers several routes for increasing the genetic gain or efficiency of plant breeding programmes. In various species of livestock, there is empirical evidence of increased rates of genetic gain from the use of genomic selection to target different aspects of the breeder's equation. Accurate predictions of genomic breeding value are central to this, and the design of training sets is in turn central to achieving sufficient levels of accuracy. In summary, small numbers of close relatives and very large numbers of distant relatives are expected to enable predictions with higher accuracy. To quantify the effect of some of the properties of training sets on the accuracy of genomic selection in crops, we performed an extensive field-based winter wheat trial. In summary, this trial involved the construction of 44 F_2:4 bi- and tri-parental populations, from which 2992 lines were grown on four field locations and yield was measured. For each line, genotype data were generated for 25 K segregating SNP markers. The overall heritability of yield was estimated to 0.65, and estimates within individual families ranged between 0.10 and 0.85. Genomic prediction accuracies of yield BLUEs were 0.125-0.127 using two different cross-validation approaches and generally increased with training set size. Using related crosses in training and validation sets generally resulted in higher prediction accuracies than using unrelated crosses. The results of this study emphasise the importance of the training panel design in relation to the genetic material to which the resulting prediction model is to be applied.

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.

Assessing European Wheat Sensitivities to Parastagonospora nodorum Necrotrophic Effectors and Fine-Mapping the Snn3-B1 Locus Conferring Sensitivity to the Effector SnTox3

Parastagonospora nodorum is a necrotrophic fungal pathogen of wheat (Triticum aestivum L.), one of the world's most important crops. P. nodorum mediates host cell death using proteinaceous necrotrophic effectors, presumably liberating nutrients that allow the infection process to continue. The identification of pathogen effectors has allowed host genetic resistance mechanisms to be separated into their constituent parts. In P. nodorum, three proteinaceous effectors have been cloned: SnToxA, SnTox1, and SnTox3. Here, we survey sensitivity to all three effectors in a panel of 480 European wheat varieties, and fine-map the wheat SnTox3 sensitivity locus Snn3-B1 using genome-wide association scans (GWAS) and an eight-founder wheat multi-parent advanced generation inter-cross (MAGIC) population. Using a Bonferroni corrected P ≤ 0.05 significance threshold, GWAS identified 10 significant markers defining a single locus, Snn3-B1, located on the short arm of chromosome 5B explaining 32% of the phenotypic variation [peak single nucleotide polymorphisms (SNPs), Excalibur_c47452_183 and GENE-3324_338, -log₁₀P = 20.44]. Single marker analysis of SnTox3 sensitivity in the MAGIC population located Snn3-B1 via five significant SNPs, defining a 6.2-kb region that included the two peak SNPs identified in the association mapping panel. Accordingly, SNP Excalibur_c47452_183 was converted to the KASP genotyping system, and validated by screening a subset of 95 wheat varieties, providing a valuable resource for marker assisted breeding and for further genetic investigation. In addition, composite interval mapping in the MAGIC population identified six minor SnTox3 sensitivity quantitative trait loci, on chromosomes 2A (QTox3.niab-2A.1, P-value = 9.17^-7), 2B (QTox3.niab-2B.1, P = 0.018), 3B (QTox3.niab-3B.1, P = 48.51^-4), 4D (QTox3.niab-4D.1, P = 0.028), 6A (QTox3.niab-6A.1, P = 8.51^-4), and 7B (QTox3.niab-7B.1, P = 0.020), each accounting for between 3.1 and 6.0 % of the phenotypic variance. Collectively, the outcomes of this study provides breeders with knowledge and resources regarding the sensitivity of European wheat germplasm to P. nodorum effectors, as well as simple diagnostic markers for determining allelic state at Snn3-B1.

Assessing European Wheat Sensitivities to Parastagonospora nodorum Necrotrophic Effectors and Fine-Mapping the Snn3-B1 Locus Conferring Sensitivity to the Effector SnTox3

Parastagonospora nodorum is a necrotrophic fungal pathogen of wheat (Triticum aestivum L.), one of the world's most important crops. P. nodorum mediates host cell death using proteinaceous necrotrophic effectors, presumably liberating nutrients that allow the infection process to continue. The identification of pathogen effectors has allowed host genetic resistance mechanisms to be separated into their constituent parts. In P. nodorum, three proteinaceous effectors have been cloned: SnToxA, SnTox1, and SnTox3. Here, we survey sensitivity to all three effectors in a panel of 480 European wheat varieties, and fine-map the wheat SnTox3 sensitivity locus Snn3-B1 using genome-wide association scans (GWAS) and an eight-founder wheat multi-parent advanced generation inter-cross (MAGIC) population. Using a Bonferroni corrected P ≤ 0.05 significance threshold, GWAS identified 10 significant markers defining a single locus, Snn3-B1, located on the short arm of chromosome 5B explaining 32% of the phenotypic variation [peak single nucleotide polymorphisms (SNPs), Excalibur_c47452_183 and GENE-3324_338, -log₁₀P = 20.44]. Single marker analysis of SnTox3 sensitivity in the MAGIC population located Snn3-B1 via five significant SNPs, defining a 6.2-kb region that included the two peak SNPs identified in the association mapping panel. Accordingly, SNP Excalibur_c47452_183 was converted to the KASP genotyping system, and validated by screening a subset of 95 wheat varieties, providing a valuable resource for marker assisted breeding and for further genetic investigation. In addition, composite interval mapping in the MAGIC population identified six minor SnTox3 sensitivity quantitative trait loci, on chromosomes 2A (QTox3.niab-2A.1, P-value = 9.17^-7), 2B (QTox3.niab-2B.1, P = 0.018), 3B (QTox3.niab-3B.1, P = 48.51^-4), 4D (QTox3.niab-4D.1, P = 0.028), 6A (QTox3.niab-6A.1, P = 8.51^-4), and 7B (QTox3.niab-7B.1, P = 0.020), each accounting for between 3.1 and 6.0 % of the phenotypic variance. Collectively, the outcomes of this study provides breeders with knowledge and resources regarding the sensitivity of European wheat germplasm to P. nodorum effectors, as well as simple diagnostic markers for determining allelic state at Snn3-B1.

Maximizing the potential of multi-parental crop populations

Most agriculturally significant crop traits are quantitatively inherited which limits the ease and efficiency of trait dissection. Multi-parent populations overcome the limitations of traditional trait mapping and offer new potential to accurately define the genetic basis of complex crop traits. The increasing popularity and use of nested association mapping (NAM) and multi-parent advanced generation intercross (MAGIC) populations raises questions about the optimal design and allocation of resources in their creation. In this paper we review strategies for the creation of multi-parent populations and describe two complementary in silico studies addressing the design and construction of NAM and MAGIC populations. The first simulates the selection of diverse founder parents and the second the influence of multi-parent crossing schemes (and number of founders) on haplotype creation and diversity. We present and apply two open software resources to simulate alternate strategies for the development of multi-parent populations.

A highly recombined, high-density, eight-founder wheat MAGIC map reveals extensive segregation distortion and genomic locations of introgression segments

Multiparent Advanced Generation Intercross (MAGIC) mapping populations offer unique opportunities and challenges for marker and QTL mapping in crop species. We have constructed the first eight-parent MAGIC genetic map for wheat, comprising 18 601 SNP markers. We validated the accuracy of our map against the wheat genome sequence and found an improvement in accuracy compared to published genetic maps. Our map shows a notable increase in precision resulting from the three generations of intercrossing required to create the population. This is most pronounced in the pericentromeric regions of the chromosomes. Sixteen percent of mapped markers exhibited segregation distortion (SD) with many occurring in long (>20 cM) blocks. Some of the longest and most distorted blocks were collinear with noncentromeric high-marker-density regions of the genome, suggesting they were candidates for introgression fragments introduced into the bread wheat gene pool from other grass species. We investigated two of these linkage blocks in detail and found strong evidence that one on chromosome 4AL, showing SD against the founder Robigus, is an interspecific introgression fragment. The completed map is available from http://www.niab.com/pages/id/326/Resources.

Determining Phenological Patterns Associated with the Onset of Senescence in a Wheat MAGIC Mapping Population

The appropriate timing of developmental transitions is critical for adapting many crops to their local climatic conditions. Therefore, understanding the genetic basis of different aspects of phenology could be useful in highlighting mechanisms underpinning adaptation, with implications in breeding for climate change. For bread wheat (Triticum aestivum), the transition from vegetative to reproductive growth, the start and rate of leaf senescence and the relative timing of different stages of flowering and grain filling all contribute to plant performance. In this study we screened under Smart house conditions a large, multi-founder "NIAB elite MAGIC" wheat population, to evaluate the genetic elements that influence the timing of developmental stages in European elite varieties. This panel of recombinant inbred lines was derived from eight parents that are or recently have been grown commercially in the UK and Northern Europe. We undertook a detailed temporal phenotypic analysis under Smart house conditions of the population and its parents, to try to identify known or novel Quantitative Trait Loci associated with variation in the timing of key phenological stages in senescence. This analysis resulted in the detection of QTL interactions with novel traits such the time between "half of ear emergence above flag leaf ligule" and the onset of senescence at the flag leaf as well as traits associated with plant morphology such as stem height. In addition, strong correlations between several traits and the onset of senescence of the flag leaf were identified. This work establishes the value of systematically phenotyping genetically unstructured populations to reveal the genetic architecture underlying morphological variation in commercial wheat.

Applying association mapping and genomic selection to the dissection of key traits in elite European wheat

We show the application of association mapping and genomic selection for key breeding targets using a large panel of elite winter wheat varieties and a large volume of agronomic data. The heightening urgency to increase wheat production in line with the needs of a growing population, and in the face of climatic uncertainty, mean new approaches, including association mapping (AM) and genomic selection (GS) need to be validated and applied in wheat breeding. Key adaptive responses are the cornerstone of regional breeding. There is evidence that new ideotypes for long-standing traits such as flowering time may be required. In order to detect targets for future marker-assisted improvement and validate the practical application of GS for wheat breeding we genotyped 376 elite wheat varieties with 3,046 DArT, single nucleotide polymorphism and gene markers and measured seven traits in replicated yield trials over 2 years in France, Germany and the UK. The scale of the phenotyping exceeds the breadth of previous AM and GS studies in these key economic wheat production regions of Northern Europe. Mixed-linear modelling (MLM) detected significant marker-trait associations across and within regions. Genomic prediction using elastic net gave low to high prediction accuracies depending on the trait, and could be experimentally increased by modifying the constituents of the training population (TP). We also tested the use of differentially penalised regression to integrate candidate gene and genome-wide markers to predict traits, demonstrating the validity and simplicity of this approach. Overall, our results suggest that whilst AM offers potential for application in both research and breeding, GS represents an exciting opportunity to select key traits, and that optimisation of the TP is crucial to its successful implementation.

An eight-parent multiparent advanced generation inter-cross population for winter-sown wheat: creation, properties, and validation

MAGIC populations represent one of a new generation of crop genetic mapping resources combining high genetic recombination and diversity. We describe the creation and validation of an eight-parent MAGIC population consisting of 1091 F7 lines of winter-sown wheat (Triticum aestivum L.). Analyses based on genotypes from a 90,000-single nucleotide polymorphism (SNP) array find the population to be well-suited as a platform for fine-mapping quantitative trait loci (QTL) and gene isolation. Patterns of linkage disequilibrium (LD) show the population to be highly recombined; genetic marker diversity among the founders was 74% of that captured in a larger set of 64 wheat varieties, and 54% of SNPs segregating among the 64 lines also segregated among the eight founder lines. In contrast, a commonly used reference bi-parental population had only 54% of the diversity of the 64 varieties with 27% of SNPs segregating. We demonstrate the potential of this MAGIC resource by identifying a highly diagnostic marker for the morphological character "awn presence/absence" and independently validate it in an association-mapping panel. These analyses show this large, diverse, and highly recombined MAGIC population to be a powerful resource for the genetic dissection of target traits in wheat, and it is well-placed to efficiently exploit ongoing advances in phenomics and genomics. Genetic marker and trait data, together with instructions for access to seed, are available at http://www.niab.com/MAGIC/.

Rht-1 and Ppd-D1 associations with height, GA sensitivity, and days to heading in a worldwide bread wheat collection

Reduced height (Rht)-1 and Photoperiod (Ppd) have major effects on the adaptability of bread wheat (Triticum aestivum) to specific environments. Ppd-D1a is a photoperiod insensitive allele that reduces time to flowering. The gibberellin (GA) insensitive alleles Rht-B1b and Rht-D1b shorten plant stature and were important components of the 'green revolution'. Two additional Rht-B1 alleles were recently identified that contain a 160 or 197 bp insertion upstream of the coding region and may affect plant height or GA sensitivity Wilhelm et al. (Theor Appl Gen doi: 10.1007/s00122-013-2088-7 , 2013b). We determined the frequency of the five alleles in a worldwide core collection of 372 wheat accessions (372CC) and estimated their effects on height, days to heading, and GA sensitivity when the collection was grown in pots outdoors or in the glasshouse. This revealed that each allele was widespread geographically with frequencies ranging from 0.12 to 0.25. Ppd-D1a was associated with significant (p ≤ 0.05) reductions in days to heading and height relative to photoperiod sensitive Ppd-D1b. Relative to wild type, Rht-B1b and Rht-D1b each resulted in significant reductions in height (approximately 30 %) and GA sensitivity. The 160 and 197 bp alleles were associated with significant height reductions of 18 and 12 %, respectively, and with non-significant reductions in GA sensitivity relative to wild type. Two statistical methods were developed and used to estimate GA sensitivity of the 372CC accessions, but novel GA insensitive alleles were not identified. Further characterization of the Rht-B1 insertion alleles is required, but our results suggest these may enable fine adjustments in plant height.

Haplotype dictionary for the Rht-1 loci in wheat

The introduction of Reduced height (Rht)-B1b and Rht-D1b into bread wheat (Triticum aestivum) varieties was a key component of the 'green revolution' and today these alleles are the primary sources of semi-dwarfism in wheat. The Rht-1 loci encode DELLA proteins, which are transcription factors that affect plant growth and stress tolerance. In bread wheat, Rht-D1b and Rht-B1b influence resistance to the disease Fusarium Head Blight. To identify Rht-1 variants, locus specific primers were developed and used to sequence the entire open reading frame (ORF) and 1.7 kb of the 5' and 0.5 kb of the 3' flanking regions of Rht-A1 (Rht-A1+f), Rht-B1 (Rht-B1+f), and Rht-D1 (Rht-D1+f) in bread wheat (36 sequences from each genome) and tetraploid and diploid wheat (TDW) (one to three sequences from each genome). Among the bread wheat accessions, the Rht-A1+f and Rht-D1+f sequences contained relatively low genetic diversity and few haplotypes relative to the Rht-B1+f sequences. The TDW accessions were relatively rich in genetic diversity and contained the majority of the polymorphic sites. Novel polymorphisms, relative to 'Chinese Spring', discovered among the accessions include 160 and 197 bp insertions 5' of Rht-B1 and a frameshift in the Rht-B1 ORF. Quantitative real-time PCR using shoot and leaf tissue from 5-day-old seedlings of genotypes lacking or containing the 5' insertions revealed no major effect on Rht-B1 transcript accumulation. This research provides insights into the genetic diversity present at the Rht-1 loci in modern bread wheat and in relation to ancestral wheat accessions.

Short, natural, and extended photoperiod response in BC2F4 lines of bread wheat with different photoperiod-1 (Ppd-1) alleles

Flowering is a critical period in the life cycle of flowering plant species, resulting in an irreversible commitment of significant resources. Wheat is photoperiod sensitive, flowering only when daylength surpasses a critical length; however, photoperiod insensitivity (PI) has been selected by plant breeders for >40 years to enhance yield in certain environments. Control of flowering time has been greatly facilitated by the development of molecular markers for the Photoperiod-1 (Ppd-1) homeoloci, on the group 2 chromosomes. In the current study, an allelic series of BC2F4 lines in the winter wheat cultivars 'Robigus' and 'Alchemy' was developed to elucidate the influence on flowering of eight gene variants from the B- and D-genomes of bread wheat and the A-genome of durum wheat. Allele effects were tested in short, natural, and extended photoperiods in the field and controlled environments. Across genetic background and treatment, the D-genome PI allele, Ppd-D1a, had a more potent effect on reducing flowering time than Ppd-B1a. However, there was significant donor allele effect for both Ppd-D1a and Ppd-B1a, suggesting the presence of linked modifier genes and/or additional sources of latent sensitivity. Development of Ppd-A1a BC2F4 lines derived from synthetic hexaploid wheat provided an opportunity to compare directly the flowering time effect of the A-genome allele from durum with the B- and D-genome variants from bread wheat for the first time. Analyses indicated that the reducing effect of Ppd-A1a is comparable with that of Ppd-D1a, confirming it as a useful alternative source of PI.

Association mapping of partitioning loci in barley

Background

Association mapping, initially developed in human disease genetics, is now being applied to plant species. The model species Arabidopsis provided some of the first examples of association mapping in plants, identifying previously cloned flowering time genes, despite high population sub-structure. More recently, association genetics has been applied to barley, where breeding activity has resulted in a high degree of population sub-structure. A major genotypic division within barley is that between winter- and spring-sown varieties, which differ in their requirement for vernalization to promote subsequent flowering. To date, all attempts to validate association genetics in barley by identifying major flowering time loci that control vernalization requirement (VRN-H1 and VRN-H2) have failed. Here, we validate the use of association genetics in barley by identifying VRN-H1 and VRN-H2, despite their prominent role in determining population sub-structure.

Results

By taking barley as a typical inbreeding crop, and seasonal growth habit as a major partitioning phenotype, we develop an association mapping approach which successfully identifies VRN-H1 and VRN-H2, the underlying loci largely responsible for this agronomic division. We find a combination of Structured Association followed by Genomic Control to correct for population structure and inflation of the test statistic, resolved significant associations only with VRN-H1 and the VRN-H2 candidate genes, as well as two genes closely linked to VRN-H1 (HvCSFs1 and HvPHYC).

Conclusion

We show that, after employing appropriate statistical methods to correct for population sub-structure, the genome-wide partitioning effect of allelic status at VRN-H1 and VRN-H2 does not result in the high levels of spurious association expected to occur in highly structured samples. Furthermore, we demonstrate that both VRN-H1 and the candidate VRN-H2 genes can be identified using association mapping. Discrimination between intragenic VRN-H1 markers was achieved, indicating that candidate causative polymorphisms may be discerned and prioritised within a larger set of positive associations. This proof of concept study demonstrates the feasibility of association mapping in barley, even within highly structured populations. A major advantage of this method is that it does not require large numbers of genome-wide markers, and is therefore suitable for fine mapping and candidate gene evaluation, especially in species for which large numbers of genetic markers are either unavailable or too costly.

Genomewide linkage study in 1,176 affected sister pair families identifies a significant susceptibility locus for endometriosis on chromosome 10q26

Endometriosis is a common gynecological disease that affects up to 10% of women in their reproductive years. It causes pelvic pain, severe dysmenorrhea, and subfertility. The disease is defined as the presence of tissue resembling endometrium in sites outside the uterus. Its cause remains uncertain despite >50 years of hypothesis-driven research, and thus the therapeutic options are limited. Disease predisposition is inherited as a complex genetic trait, which provides an alternative route to understanding the disease. We seek to identify susceptibility loci, using a positional-cloning approach that starts with linkage analysis to identify genomic regions likely to harbor these genes. We conducted a linkage study of 1,176 families (931 from an Australian group and 245 from a U.K. group), each with at least two members--mainly affected sister pairs--with surgically diagnosed disease. We have identified a region of significant linkage on chromosome 10q26 (maximum LOD score [MLS] of 3.09; genomewide P = .047) and another region of suggestive linkage on chromosome 20p13 (MLS = 2.09). Minor peaks (with MLS > 1.0) were found on chromosomes 2, 6, 7, 8, 12, 14, 15, and 17. This is the first report of linkage to a major locus for endometriosis. The findings will facilitate discovery of novel positional genetic variants that influence the risk of developing this debilitating disease. Greater understanding of the aberrant cellular and molecular mechanisms involved in the etiology and pathophysiology of endometriosis should lead to better diagnostic methods and targeted treatments.

Genome-wide linkage scan on estimated breeding values for a quantitative trait

BACKGROUND

A genome-wide linkage scan was performed on Replicate 1 of the simulated data for fasting triglyceride levels. The aim of this study was to implement mixed-model methodology to estimate breeding values for each individual for this trait and to assess the merit of these breeding values in linkage analysis. These breeding values utilize all the pedigree information, and the genetic and phenotypic correlations with other measured traits across the two cohorts. A genome-wide linkage scan was run on both the new breeding value traits and the original traits.

RESULTS

Using breeding values, a maximum LOD of 7.78 was found on chromosome 5 at a position very close to a gene underlying the triglyceride levels. This effect was not detected using the original trait.

CONCLUSION

The results imply that estimating breeding values may be a suitable method of deriving traits for use in genome-wide scans.

Evidence for a major gene for bone mineral density in idiopathic osteoporotic families

Although there have been a number of studies indicating a heritable component for osteoporosis in middle to late adulthood, the etiology of osteoporosis in young people is uncertain. The present study aims to evaluate the extent to which genetic factors influence familial resemblance for bone mineral density (BMD) in families ascertained on the basis of young osteoporotic probands. The sample comprises eight families (74 total individuals) that were identified through a proband under the age of 35 years with a history of two or more fractures and a spinal bone density of at least 2.5 SDs below the mean for age and sex (Z score). Secondary causes of osteoporosis were excluded in the probands. In total, 27% (18/66) of the probands' relatives had osteoporosis and an additional 30% (20/66) had osteopenia. Classical segregation analysis was performed to evaluate the extent to which a genetic etiology could account for familial resemblance in these families. The results indicate a major gene of codominant inheritance for spinal BMD. Model-fitting comparisons revealed no support for environmental effects or for polygenic inheritance.

The effect of gametic-phase disequilibrium on the prediction of response to recurrent selection in plants

Selection reduces additive genetic variation by generating gametic-phase disequilibrium, a phenomenon largely ignored when predicting response in plant breeding programs. The development of gametic-phase disequilibrium is here taken into account when predicting the response to selection for various schemes of recurrent selection applicable to plant populations. A general program permits prediction of response to selection from schemes of recurrent selection in which two or more rounds of selection occur in each cycle. An example from Sugar Beet, with alternate rounds of half-sib and S1 family selection, is illustrated. It is shown that failure to take into account the effects of gametic-phase disequilibrium can result in substantial overestimation of the response to selection as well as to changes in rank of the merits of alternative breeding schemes. For a given scheme, ignoring gametic-phase disequilibrium has only small effects on defining the optimum allocation of plots and the numbers of families tested.