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

An eight-founder wheat MAGIC population allows fine-mapping of flowering time loci and provides novel insights into the genetic control of flowering time

Flowering time synchronizes reproductive development with favorable environmental conditions to optimize yield. Improved understanding of the genetic control of flowering will help optimize varietal adaptation to future agricultural systems under climate change. Here, we investigate the genetic basis of flowering time in winter wheat (Triticum aestivum L.) using an eight-founder multi-parent advanced generation intercross (MAGIC) population. Flowering time data was collected from field trials across six growing seasons in the United Kingdom, followed by genetic analysis using a combination of linear modelling, simple interval mapping and composite interval mapping, using either single markers or founder haplotype probabilities. We detected 57 quantitative trait loci (QTL) across three growth stages linked to flowering time, of which 17 QTL were identified only when the major photoperiod response locus Ppd-D1 was included as a covariate. Of the 57 loci, ten were identified using all genetic mapping approaches and classified as 'major' QTL, including homoeologous loci on chromosomes 1B and 1D, and 4A and 4B. Additional Earliness per se flowering time QTL were identified, along with growth stage- and year-specific effects. Furthermore, six of the main-effect QTL were found to interact epistatically with Ppd-D1. Finally, we exploited residual heterozygosity in the MAGIC recombinant inbred lines to Mendelize the Earliness per se QTL QFt.niab-5A.03, which was confirmed to modulate flowering time by at least four days. This work provides detailed understanding of the genetic control of phenological variation within varieties relevant to the north-western European wheat genepool, aiding informed manipulation of flowering time in wheat breeding.

Diversification of quantitative morphological traits in wheat

BACKGROUND AND AIMS

The development and morphology of crop plants have been profoundly altered by evolution under cultivation, initially through unconscious selection, without deliberate foresight, and later by directed breeding. Wild wheats remain an important potential source of variation for modern breeders; however, the sequence and timing of morphological changes during domestication are not fully resolved.

METHODS

We grew and measured 142 wheat accessions representing different stages in wheat evolution, including three independent domestication events, and compared their morphological traits to define the morphospace of each group.

KEY RESULTS

The results show that wild and domesticated wheats have overlapping morphospaces, but each also occupies a distinct area of morphospace from one another. Polyploid formation in wheat increased leaf biomass and seed weight but had its largest effects on tiller loss. Domestication continued to increase the sizes of wheat leaves and seeds and made wheat grow taller, with more erect architecture. Associated changes to the biomass of domesticated wheats generated more grains and achieved higher yields. Landrace improvement subsequently decreased the numbers of tillers and spikes, to focus resource allocation to the main stem, accompanied by a thicker main stem and larger flag leaves. During the Green Revolution, wheat height was reduced to increase the harvest index and therefore yield. Modern wheats also have more erect leaves and larger flower biomass proportions than landraces.

CONCLUSIONS

Quantitative trait history in wheat differs by trait. Some trait values show progressive changes in the same direction (e.g. leaf size, grain weight), whereas others change in a punctuated way at particular stages (e.g. canopy architecture), and other trait values switch directions during wheat evolution (e.g. plant height, flower biomass proportion). Agronomically valued domestication traits arose during different stages of wheat history, such that modern wheats are the product of >10 000 years of morphological evolution.

Genotyping of SNPs in bread wheat at reduced cost from pooled experiments and imputation

KEY MESSAGE

Pooling and imputation are computational methods that can be combined for achieving cost-effective and accurate high-density genotyping of both common and rare variants, as demonstrated in a MAGIC wheat population. The plant breeding industry has shown growing interest in using the genotype data of relevant markers for performing selection of new competitive varieties. The selection usually benefits from large amounts of marker data, and it is therefore crucial to dispose of data collection methods that are both cost-effective and reliable. Computational methods such as genotype imputation have been proposed earlier in several plant science studies for addressing the cost challenge. Genotype imputation methods have though been used more frequently and investigated more extensively in human genetics research. The various algorithms that exist have shown lower accuracy at inferring the genotype of genetic variants occurring at low frequency, while these rare variants can have great significance and impact in the genetic studies that underlie selection. In contrast, pooling is a technique that can efficiently identify low-frequency items in a population, and it has been successfully used for detecting the samples that carry rare variants in a population. In this study, we propose to combine pooling and imputation and demonstrate this by simulating a hypothetical microarray for genotyping a population of recombinant inbred lines in a cost-effective and accurate manner, even for rare variants. We show that with an adequate imputation model, it is feasible to accurately predict the individual genotypes at lower cost than sample-wise genotyping and time-effectively. Moreover, we provide code resources for reproducing the results presented in this study in the form of a containerized workflow.

Fast integration and accumulation of beneficial breeding alleles through an AB-NAMIC strategy in wheat

Wheat (Triticum aestivum) is among the most important staple crops for safeguarding the food security of the growing world population. To bridge the gap between genebank diversity and breeding programs, we developed an advanced backcross-nested association mapping plus inter-crossed population (AB-NAMIC) by crossing three popular wheat cultivars as recurrent founders to 20 germplasm lines from a mini core collection. Selective backcrossing combined with selection against undesirable traits and extensive crossing within and between sub-populations created new opportunities to detect unknown genes and increase the frequency of beneficial alleles in the AB-NAMIC population. We performed phenotyping of 590 AB-NAMIC lines and a natural panel of 476 cultivars for six consecutive growing seasons and genotyped these 1066 lines with a 660K SNP array. Genome-wide association studies of both panels for plant development and yield traits demonstrated improved power to detect rare alleles and loci with medium genetic effects in AB-NAMIC. Notably, genome-wide association studies in AB-NAMIC detected the candidate gene TaSWEET6-7B (TraesCS7B03G1216700), which has high homology to the rice SWEET6b gene and exerts strong effects on adaptation and yield traits. The commercial release of two derived AB-NAMIC lines attests to its direct applicability in wheat improvement. Valuable information on genome-wide association study mapping, candidate genes, and their haplotypes for breeding traits are available through WheatGAB. Our research provides an excellent framework for fast-tracking exploration and accumulation of beneficial alleles stored in genebanks.

Identifying genomic regions determining shoot and root traits related to nitrogen uptake efficiency in a multiparent advanced generation intercross (MAGIC) winter wheat population in a high-throughput phenotyping facility

In the context of a continuously increasing human population that needs to be fed, with environmental protection in mind, nitrogen use efficiency (NUE) improvement is becoming very important. To understand the natural variation of traits linked to nitrogen uptake efficiency (UPE), one component of NUE, the multiparent advanced generation intercross (MAGIC) winter wheat population WM-800 was phenotyped under two contrasting nitrogen (N) levels in a high-throughput phenotyping facility for six weeks. Three biomass-related, three root-related, and two reflectance-related traits were measured weekly under each treatment. Subsequently, the population was genetically analysed using a total of 13,060 polymorphic haplotypes and singular SNPs for a genome-wide association study (GWAS). In total, we detected 543 quantitative trait loci (QTL) across all time points and traits, which were pooled into 42 stable QTL (sQTL; present in at least three of the six weeks). Besides Rht-B1 and Rht-D1, candidate genes playing a role in gibberellic acid-regulated growth and nitrate transporter genes from the NPF gene family, like NRT 1.1, were linked to sQTL. Two novel sQTL on chromosomes 5 A and 6D showed pleiotropic effects on several traits. The high number of N-specific sQTL indicates that selection for UPE is useful specifically under N-limited conditions.

Longer epidermal cells underlie a quantitative source of variation in wheat flag leaf size

The wheat flag leaf is the main contributor of photosynthetic assimilates to developing grains. Understanding how canopy architecture strategies affect source strength and yield will aid improved crop design. We used an eight-founder population to investigate the genetic architecture of flag leaf area, length, width and angle in European wheat. For the strongest genetic locus identified, we subsequently created a near-isogenic line (NIL) pair for more detailed investigation across seven test environments. Genetic control of traits investigated was highly polygenic, with colocalisation of replicated quantitative trait loci (QTL) for one or more traits identifying 24 loci. For QTL QFll.niab-5A.1 (FLL5A), development of a NIL pair found the FLL5A+ allele commonly conferred a c. 7% increase in flag and second leaf length and a more erect leaf angle, resulting in higher flag and/or second leaf area. Increased FLL5A-mediated flag leaf length was associated with: (1) longer pavement cells and (2) larger stomata at lower density, with a trend for decreased maximum stomatal conductance (G_smax ) per unit leaf area. For FLL5A, cell size rather than number predominantly determined leaf length. The observed trade-offs between leaf size and stomatal morphology highlight the need for future studies to consider these traits at the whole-leaf level.

Wheat EARLY FLOWERING 3 affects heading date without disrupting circadian oscillations

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

A Haplotype-Based GWAS Identified Trait-Improving QTL Alleles Controlling Agronomic Traits under Contrasting Nitrogen Fertilization Treatments in the MAGIC Wheat Population WM-800

The multi-parent-advanced-generation-intercross (MAGIC) population WM-800 was developed by intercrossing eight modern winter wheat cultivars to enhance the genetic diversity present in breeding populations. We cultivated WM-800 during two seasons in seven environments under two contrasting nitrogen fertilization treatments. WM-800 lines exhibited highly significant differences between treatments, as well as high heritabilities among the seven agronomic traits studied. The highest-yielding WM-line achieved an average yield increase of 4.40 dt/ha (5.2%) compared to the best founder cultivar Tobak. The subsequent genome-wide-association-study (GWAS), which was based on haplotypes, located QTL for seven agronomic traits including grain yield. In total, 40, 51, and 46 QTL were detected under low, high, and across nitrogen treatments, respectively. For example, the effect of QYLD_3A could be associated with the haplotype allele of cultivar Julius increasing yield by an average of 4.47 dt/ha (5.2%). A novel QTL on chromosome 2B exhibited pleiotropic effects, acting simultaneously on three-grain yield components (ears-per-square-meter, grains-per-ear, and thousand-grain-weight) and plant-height. These effects may be explained by a member of the nitrate-transporter-1 (NRT1)/peptide-family, TaNPF5.34, located 1.05 Mb apart. The WM-800 lines and favorable QTL haplotypes, associated with yield improvements, are currently implemented in wheat breeding programs to develop advanced nitrogen-use efficient wheat cultivars.

statgenMPP: an R package implementing an IBD-based mixed model approach for QTL mapping in a wide range of multi-parent populations

MOTIVATION

Multi-parent populations (MPPs) are popular for QTL mapping because they combine wide genetic diversity in parents with easy control of population structure, but a limited number of software tools for QTL mapping are specifically developed for general MPP designs.

RESULTS

We developed an R package called statgenMPP, adopting a unified identity-by-descent (IBD)-based mixed model approach for QTL analysis in MPPs. The package offers easy-to-use functionalities of IBD calculations, mixed model solutions and visualizations for QTL mapping in a wide range of MPP designs, including diallele, nested-association mapping populations, multi-parent advanced genetic inter-cross populations and other complicated MPPs with known crossing schemes.

AVAILABILITY AND IMPLEMENTATION

The R package statgenMPP is open-source and freely available on CRAN at https://CRAN.R-project.org/package=statgenMPP.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Stomata on the abaxial and adaxial leaf surfaces contribute differently to leaf gas exchange and photosynthesis in wheat

Although stomata are typically found in greater numbers on the abaxial surface, wheat flag leaves have greater densities on the adaxial surface. We determine the impact of this less common stomatal patterning on gaseous fluxes using a novel chamber that simultaneously measures both leaf surfaces. Using a combination of differential illuminations and CO₂ concentrations at each leaf surface, we found that mesophyll cells associated with the adaxial leaf surface have a higher photosynthetic capacity than those associated with the abaxial leaf surface, which is supported by an increased stomatal conductance (driven by differences in stomatal density). When vertical gas flux at the abaxial leaf surface was blocked, no compensation by adaxial stomata was observed, suggesting each surface operates independently. Similar stomatal kinetics suggested some co-ordination between the two surfaces, but factors other than light intensity played a role in these responses. Higher photosynthetic capacity on the adaxial surface facilitates greater carbon assimilation, along with higher adaxial stomatal conductance, which would also support greater evaporative leaf cooling to maintain optimal leaf temperatures for photosynthesis. Furthermore, abaxial gas exchange contributed c. 50% to leaf photosynthesis and therefore represents an important contributor to overall leaf gas exchange.

Overexpression of the WAPO-A1 gene increases the number of spikelets per spike in bread wheat

Two homoeologous QTLs for number of spikelets per spike (SPS) were mapped on chromosomes 7AL and 7BL using two wheat MAGIC populations. Sets of lines contrasting for the QTL on 7AL were developed which allowed for the validation and fine mapping of the 7AL QTL and for the identification of a previously described candidate gene, WHEAT ORTHOLOG OF APO1 (WAPO1). Using transgenic overexpression in both a low and a high SPS line, we provide a functional validation for the role of this gene in determining SPS also in hexaploid wheat. We show that the expression levels of this gene positively correlate with SPS in multiple MAGIC founder lines under field conditions as well as in transgenic lines grown in the greenhouse. This work highlights the potential use of WAPO1 in hexaploid wheat for further yield increases. The impact of WAPO1 and SPS on yield depends on other genetic and environmental factors, hence, will require a finely balanced expression level to avoid the development of detrimental pleiotropic phenotypes.

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

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

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

Image 1.

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 R2 ≥ .5 with the same QTL. Genome-wide association studies identified marker-trait associations for all four traits. For HFN (h 2 = .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 2 = 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 2 = 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.

Genetic mapping and prediction of flowering time and plant height in a maize Stiff Stalk MAGIC population

The Stiff Stalk heterotic pool is a foundation of US maize seed parent germplasm and has been heavily utilized by both public and private maize breeders since its inception in the 1930's. Flowering time and plant height are critical characteristics for both inbred parents and their test crossed hybrid progeny. To study these traits, a six parent multiparent advanced generation intercross (MAGIC) population was developed including maize inbred lines B73, B84, PHB47 (B37 type), LH145 (B14 type), PHJ40 (novel early Stiff Stalk), and NKH8431 (B73/B14 type). A set of 779 doubled haploid lines were evaluated for flowering time and plant height in two field replicates in 2016 and 2017, and a subset of 689 and 561 doubled haploid lines were crossed to two testers, respectively, and evaluated as hybrids in two locations in 2018 and 2019 using an incomplete block design. Markers were derived from a Practical Haplotype Graph built from the founder whole genome assemblies and genotype-by-sequencing and exome capture-based sequencing of the population. Genetic mapping utilizing an update to R/qtl2 revealed differing profiles of significant loci for both traits between 635 of the DH lines and two sets of 570 and 471 derived hybrids. Genomic prediction was used to test the feasibility of predicting hybrid phenotypes based on the per se data. Predictive abilities were highest on direct models trained using the data they would predict (0.55 to 0.63), and indirect models trained using per se data to predict hybrid traits had slightly lower predictive abilities (0.49 to 0.55). Overall, this finding is consistent with the overlapping and non-overlapping significant QTL found within the per se and hybrid populations and suggests that selections for phenology traits can be made effectively on doubled haploid lines before hybrid data is available.

Mapping of dwarfing QTL of Ari1327, a semi-dwarf mutant of upland cotton

BACKGROUND

Upland Cotton (Gossypium hirsutum L.) has few cotton varieties suitable for mechanical harvesting. The plant height of the cultivar is one of the key features that need to modify. Hence, this study was planned to locate the QTL for plant height in a 60Co γ treated upland cotton semi-dwarf mutant Ari1327.

RESULTS

Interestingly, bulk segregant analysis (BSA) and genotyping by sequencing (GBS) methods exhibited that candidate QTL was co-located in the region of 5.80-9.66 Mb at D01 chromosome in two F2 populations. Using three InDel markers to genotype a population of 1241 individuals confirmed that the offspring's phenotype is consistent with the genotype. Comparative analysis of RNA-seq between the mutant and wild variety exhibited that Gh_D01G0592 was identified as the source of dwarfness from 200 genes. In addition, it was also revealed that the appropriate use of partial separation markers in QTL mapping can escalate linkage information.

CONCLUSIONS

Overwhelmingly, the results will provide the basis to reveal the function of candidate genes and the utilization of excellent dwarf genetic resources in the future.

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.

Meta-QTLs and candidate genes for stripe rust resistance in wheat

In bread wheat, meta-QTL analysis was conducted using 353 QTLs that were available from earlier studies. When projected onto a dense consensus map comprising 76,753 markers, only 184 QTLs with the required information, could be utilized leading to identification of 61 MQTLs spread over 18 of the 21 chromosomes (barring 5D, 6D and 7D). The range for mean R2 (PVE %) was 1.9% to 48.1%, and that of CI was 0.02 to 11.47 cM; these CIs also carried 37 Yr genes. Using these MQTLs, 385 candidate genes (CGs) were also identified. Out of these CGs, 241 encoded known R proteins and 120 showed differential expression due to stripe rust infection at the seedling stage; the remaining 24 CGs were common in the sense that they encoded R proteins as well as showed differential expression. The proteins encoded by CGs carried the following widely known domains: NBS-LRR domain, WRKY domains, ankyrin repeat domains, sugar transport domains, etc. Thirteen breeders' MQTLs (PVE > 20%) including four pairs of closely linked MQTLs are recommended for use in wheat molecular breeding, for future studies to understand the molecular mechanism of stripe rust resistance and for gene cloning.

An IBD-based mixed model approach for QTL mapping in multiparental populations

The identity-by-descent (IBD)-based mixed model approach introduced in this study can detect quantitative trait loci (QTLs) referring to the parental origin and simultaneously account for multilevel relatedness of individuals within and across families. This unified approach is proved to be a powerful approach for all kinds of multiparental population (MPP) designs. Multiparental populations (MPPs) have become popular for quantitative trait loci (QTL) detection. Tools for QTL mapping in MPPs are mostly developed for specific MPPs and do not generalize well to other MPPs. We present an IBD-based mixed model approach for QTL mapping in all kinds of MPP designs, e.g., diallel, Nested Association Mapping (NAM), and Multiparental Advanced Generation Intercross (MAGIC) designs. The first step is to compute identity-by-descent (IBD) probabilities using a general Hidden Markov model framework, called reconstructing ancestry blocks bit by bit (RABBIT). Next, functions of IBD information are used as design matrices, or genetic predictors, in a mixed model approach to estimate variance components for multiallelic genetic effects associated with parents. Family-specific residual genetic effects are added, and a polygenic effect is structured by kinship relations between individuals. Case studies of simulated diallel, NAM, and MAGIC designs proved that the advanced IBD-based multi-QTL mixed model approach incorporating both kinship relations and family-specific residual variances (IBD.MQMkin_F) is robust across a variety of MPP designs and allele segregation patterns in comparison to a widely used benchmark association mapping method, and in most cases, outperformed or behaved at least as well as other tools developed for specific MPP designs in terms of mapping power and resolution. Successful analyses of real data cases confirmed the wide applicability of our IBD-based mixed model methodology.

Crafting for a better MAGIC: systematic design and test for Multiparental Advanced Generation Inter-Cross population

Multiparental Advanced Generation Inter-Cross (MAGIC) populations are valuable crop resources with a wide array of research uses including genetic mapping of complex traits, management of genetic resources and breeding of new varieties. Multiple founders are crossed to create a rich mosaic of highly recombined founder genomes in the MAGIC recombinant inbred lines (RILs). Many variations of MAGIC population designs exist; however, a large proportion of the currently available populations have been created empirically and based on similar designs. In our evaluations of five MAGIC populations, we found that the choice of designs has a large impact on the recombination landscape in the RILs. The most popular design used in many MAGIC populations has been shown to have a bias in recombinant haplotypes and low level of unique recombinant haplotypes, and therefore is not recommended. To address this problem and provide a remedy for the future, we have developed the “magicdesign” R package for creating and testing any MAGIC population design via simulation. A Shiny app version of the package is available as well. Our “magicdesign” package provides a unifying tool and a framework for creativity and innovation in MAGIC population designs. For example, using this package, we demonstrate that MAGIC population designs can be found which are very effective in creating haplotype diversity without the requirement for very large crossing programs. Furthermore, we show that interspersing cycles of crossing with cycles of selfing is effective in increasing haplotype diversity. These approaches are applicable in species that are hard to cross or in which resources are limited.

Construction of Consensus Genetic Map With Applications in Gene Mapping of Wheat (Triticum aestivum L.) Using 90K SNP Array

Wheat is one of the most important cereal crops worldwide. A consensus map combines genetic information from multiple populations, providing an effective alternative to improve the genome coverage and marker density. In this study, we constructed a consensus map from three populations of recombinant inbred lines (RILs) of wheat using a 90K single nucleotide polymorphism (SNP) array. Phenotypic data on plant height (PH), spike length (SL), and thousand-kernel weight (TKW) was collected in six, four, and four environments in the three populations, and then used for quantitative trait locus (QTL) mapping. The mapping results obtained using the constructed consensus map were compared with previous results obtained using individual maps and previous studies on other populations. A simulation experiment was also conducted to assess the performance of QTL mapping with the consensus map. The constructed consensus map from the three populations spanned 4558.55 cM in length, with 25,667 SNPs, having high collinearity with physical map and individual maps. Based on the consensus map, 21, 27, and 19 stable QTLs were identified for PH, SL, and TKW, much more than those detected with individual maps. Four PH QTLs and six SL QTLs were likely to be novel. A putative gene called TraesCS4D02G076400 encoding gibberellin-regulated protein was identified to be the candidate gene for one major PH QTL located on 4DS, which may enrich genetic resources in wheat semi-dwarfing breeding. The simulation results indicated that the length of the confidence interval and standard errors of the QTLs detected using the consensus map were much smaller than those detected using individual maps. The consensus map constructed in this study provides the underlying genetic information for systematic mapping, comparison, and clustering of QTL, and gene discovery in wheat genetic study. The QTLs detected in this study had stable effects across environments and can be used to improve the wide adaptation of wheat cultivars through marker-assisted breeding.

The effects of crop attributes, selection, and recombination on Canadian bread wheat molecular variation

Cultivated germplasm provides an opportunity to investigate how crop agronomic traits, selection for major genes, and differences in crossing-over rates drive patterns of allelic variation. To identify how these factors correlated with allelic variation within a collection of cultivated bread wheat (Triticum aestivum L.), we generated genotypes for 388 accessions grown in Canada over the past 170 yr using filtered single nucleotide polymorphism (SNP) calls from an Illumina Wheat iSelect 90K SNP-array. Entries' breeding program, era of release, grain texture, kernel color, and growth habit contributed to allelic differentiation. Allelic diversity and linkage disequilibrium (LD) of markers flanking some major loci known to affect traits such as gluten strength, growth habit, and grain color were consistent with selective sweeps. Nonetheless, some flanking markers of major loci had low LD and high allelic diversity. Positive selection may have acted upon homoeologous genes that had significant enrichment for the gene ontology terms 'response-to-auxin' and 'response-to-wounding.' Long regions of LD, spanning approximately one-third the length of entire chromosomes, were associated with many pericentromeric regions. These regions were also characterized by low diversity. Enhancing recombination across these regions could generate novel allele combinations to accelerate Canadian wheat improvement.

Genetic control of iron bioavailability is independent from iron concentration in a diverse winter wheat mapping population

BACKGROUND

Anemia is thought to affect up to 1.6 billion people worldwide. One of the major contributors to low iron (Fe) absorption is a higher proportion of cereals compared to meats and pulse crops in people's diets. This has now become a problem in both the developed and developing world, as a result of both modern food choice and food availability. Bread wheat accounts for 20 % of the calories consumed by humans and is an important source of protein, vitamins and minerals meaning it could be a major vehicle for bringing more bioavailable Fe into the diet.

RESULTS

To investigate whether breeding for higher concentrations of Fe in wheat grains could help increase Fe absorption, a multiparent advanced generation intercross (MAGIC) population, encompassing more than 80 % of UK wheat polymorphism, was grown over two seasons in the UK. The population was phenotyped for both Fe concentration and Fe bioavailability using an established Caco-2 cell bioassay. It was found that increasing Fe concentrations in the grains was not correlated with higher Fe bioavailability and that the underlying genetic regions controlling grain Fe concentrations do not co-localise with increased Fe absorption. Furthermore, we show that phytate concentrations do not correlate with Fe bioavailability in our wheat population and thus phytate-binding is insufficient to explain the lack of correlation between Fe bioavailability and Fe concentrations in the wheat grain. Finally, we observed no (Fe bioavailability) or low (Fe concentration) correlation between years for these traits, confirming that both are under strong environmental influence.

CONCLUSIONS

This suggests that breeders will have to select not only for Fe concentrations directly in grains, but also increased bioavailability. However the use of numerous controls and replicated trials limits the practicality of adoption of screening by Caco-2 cells by many breeders.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

The 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.

Multiparental Population in Crops: Methods of Development and Dissection of Genetic Traits

Multiparental populations are located midway between association mapping that relies on germplasm collections and classic linkage analysis, based upon biparental populations. They provide several key advantages such as the possibility to include a higher number of alleles and increased level of recombination with respect to biparental populations, and more equilibrated allelic frequencies than association mapping panels. Moreover, in these populations new allele's combinations arise from recombination that may reveal transgressive phenotypes and make them a useful pre-breeding material. Here we describe the strategies for working with multiparental populations, focusing on nested association mapping populations (NAM) and multiparent advanced generation intercross populations (MAGIC). We provide details from the selection of founders, population development, and characterization to the statistical methods for genetic mapping and quantitative trait detection.

QTL mapping of adult plant and seedling resistance to leaf rust (Puccinia triticina Eriks.) in a multiparent advanced generation intercross (MAGIC) wheat population

The Bavarian MAGIC Wheat population, comprising 394 F6:8 recombinant inbred lines was phenotyped for Puccinia triticina resistance in multi-years' field trials at three locations and in a controlled environment seedling test. Simple intervall mapping revealed 19 QTL, corresponding to 11 distinct chromosomal regions. The biotrophic rust fungus Puccinia triticina is one of the most important wheat pathogens with the potential to cause yield losses up to 70%. Growing resistant cultivars is the most cost-effective and environmentally friendly way to encounter this problem. The emergence of leaf rust races being virulent against common resistance genes increases the demand for wheat varieties with novel resistances. In the past decade, the use of complex experimental populations, like multiparent advanced generation intercross (MAGIC) populations, has risen and offers great advantages for mapping resistances. The genetic diversity of multiple parents, which has been recombined over several generations, leads to a broad phenotypic diversity, suitable for high-resolution mapping of quantitative traits. In this study, interval mapping was performed to map quantitative trait loci (QTL) for leaf rust resistance in the Bavarian MAGIC Wheat population, comprising 394 F_6:8 recombinant inbred lines (RILs). Phenotypic evaluation of the RILs for adult plant resistance was carried out in field trials at three locations and two years, as well as in a controlled-environment seedling inoculation test. In total, interval mapping revealed 19 QTL, which corresponded to 11 distinct chromosomal regions controlling leaf rust resistance. Six of these regions may represent putative new QTL. Due to the elite parental material, RILs identified to be resistant to leaf rust can be easily introduced in breeding programs.

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.

Genetic Analysis Using a Multi-Parent Wheat Population Identifies Novel Sources of Septoria Tritici Blotch Resistance

Zymoseptoria tritici is the causative fungal pathogen of septoria tritici blotch (STB) disease of wheat (Triticum aestivum L.) that continuously threatens wheat crops in Ireland and throughout Europe. Under favorable conditions, STB can cause up to 50% yield losses if left untreated. STB is commonly controlled with fungicides; however, a combination of Z. tritici populations developing fungicide resistance and increased restrictions on fungicide use in the EU has led to farmers relying on fewer active substances. Consequently, this serves to drive the emergence of Z. tritici resistance against the remaining chemistries. In response, the use of resistant wheat varieties provides a more sustainable disease management strategy. However, the number of varieties offering an adequate level of resistance against STB is limited. Therefore, new sources of resistance or improved stacking of existing resistance loci are needed to develop varieties with superior agronomic performance. Here, we identified quantitative trait loci (QTL) for STB resistance in the eight-founder "NIAB Elite MAGIC" winter wheat population. The population was screened for STB response in the field under natural infection for three seasons from 2016 to 2018. Twenty-five QTL associated with STB resistance were identified in total. QTL either co-located with previously reported QTL or represent new loci underpinning STB resistance. The genomic regions identified and the linked genetic markers serve as useful resources for STB resistance breeding, supporting rapid selection of favorable alleles for the breeding of new wheat cultivars with improved STB resistance.

Accurate calling of homeoallelic genotypes of iSelect markers using inbred structured populations

MOTIVATION

Genetic map construction is a foundational step in the analysis of structured experimental populations. For markers that hybridize to several genetically similar locations, or where several alleles are present (such as in multiparental populations), current methods often discard the marker or incorrectly call the genotypes. These errors result in information loss, or incorrect genotypes that can corrupt map construction.

RESULTS

We present a new approach for simultaneously performing genetic map construction and marker calling. Our new approach allows the calling of a larger number of markers, a larger number of unique alleles per marker and the correct use of markers which hybridize to multiple genetically similar locations. We demonstrate our new approach using simulations, a biparental wheat population and an eight-parent population of spring bread wheat. Applying our method to the eight-parent population increased the number of mapped markers by 71%. We show that the new genetic map allows the investigation of synteny in ways that were not previously possible in that dataset.

AVAILABILITY AND IMPLEMENTATION

The method described in this article has been incorporated into R package mpMap2. It is available from CRAN and also from https://github.com/rohan-shah/mpMap2.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A roadmap for gene functional characterisation in crops with large genomes: Lessons from polyploid wheat

Understanding the function of genes within staple crops will accelerate crop improvement by allowing targeted breeding approaches. Despite their importance, a lack of genomic information and resources has hindered the functional characterisation of genes in major crops. The recent release of high-quality reference sequences for these crops underpins a suite of genetic and genomic resources that support basic research and breeding. For wheat, these include gene model annotations, expression atlases and gene networks that provide information about putative function. Sequenced mutant populations, improved transformation protocols and structured natural populations provide rapid methods to study gene function directly. We highlight a case study exemplifying how to integrate these resources. This review provides a helpful guide for plant scientists, especially those expanding into crop research, to capitalise on the discoveries made in Arabidopsis and other plants. This will accelerate the improvement of crops of vital importance for food and nutrition security.

Genetic analysis of wheat sensitivity to the ToxB fungal effector from Pyrenophora tritici-repentis, the causal agent of tan spot

Genetic mapping of sensitivity to the Pyrenophora tritici-repentis effector ToxB allowed development of a diagnostic genetic marker, and investigation of wheat pedigrees allowed transmission of sensitive alleles to be tracked. Tan spot, caused by the necrotrophic fungal pathogen Pyrenophora tritici-repentis, is a major disease of wheat (Triticum aestivum). Secretion of the P. tritici-repentis effector ToxB is thought to play a part in mediating infection, causing chlorosis of plant tissue. Here, genetic analysis using an association mapping panel (n = 480) and a multiparent advanced generation intercross (MAGIC) population (n founders = 8, n progeny = 643) genotyped with a 90,000 feature single nucleotide polymorphism (SNP) array found ToxB sensitivity to be highly heritable (h² ≥ 0.9), controlled predominantly by the Tsc2 locus on chromosome 2B. Genetic mapping of Tsc2 delineated a 1921-kb interval containing 104 genes in the reference genome of ToxB-insensitive variety 'Chinese Spring'. This allowed development of a co-dominant genetic marker for Tsc2 allelic state, diagnostic for ToxB sensitivity in the association mapping panel. Phenotypic and genotypic analysis in a panel of wheat varieties post-dated the association mapping panel further supported the diagnostic nature of the marker. Combining ToxB phenotype and genotypic data with wheat pedigree datasets allowed historic sources of ToxB sensitivity to be tracked, finding the variety 'Maris Dove' to likely be the historic source of sensitive Tsc2 alleles in the wheat germplasm surveyed. Exploration of the Tsc2 region gene space in the ToxB-sensitive line 'Synthetic W7984' identified candidate genes for future investigation. Additionally, a minor ToxB sensitivity QTL was identified on chromosome 2A. The resources presented here will be of immediate use for marker-assisted selection for ToxB insensitivity and the development of germplasm with additional genetic recombination within the Tsc2 region.

Genetic mapping using a wheat multi-founder population reveals a locus on chromosome 2A controlling resistance to both leaf and glume blotch caused by the necrotrophic fungal pathogen Parastagonospora nodorum

KEY MESSAGE

A locus on wheat chromosome 2A was found to control field resistance to both leaf and glume blotch caused by the necrotrophic fungal pathogen Parastagonospora nodorum. The necrotrophic fungal pathogen Parastagonospora nodorum is the causal agent of Septoria nodorum leaf blotch and glume blotch, which are common wheat (Triticum aestivum L.) diseases in humid and temperate areas. Susceptibility to Septoria nodorum leaf blotch can partly be explained by sensitivity to corresponding P. nodorum necrotrophic effectors (NEs). Susceptibility to glume blotch is also quantitative; however, the underlying genetics have not been studied in detail. Here, we genetically map resistance/susceptibility loci to leaf and glume blotch using an eight-founder wheat multiparent advanced generation intercross population. The population was assessed in six field trials across two sites and 4 years. Seedling infiltration and inoculation assays using three P. nodorum isolates were also carried out, in order to compare quantitative trait loci (QTL) identified under controlled conditions with those identified in the field. Three significant field resistance QTL were identified on chromosomes 2A and 6A, while four significant seedling resistance QTL were detected on chromosomes 2D, 5B and 7D. Among these, QSnb.niab-2A.3 for field resistance to both leaf blotch and glume blotch was detected in Norway and the UK. Colocation with a QTL for seedling reactions against culture filtrate from a Norwegian P. nodorum isolate indicated the QTL could be caused by a novel NE sensitivity. The consistency of this QTL for leaf blotch at the seedling and adult plant stages and culture filtrate infiltration was confirmed by haplotype analysis. However, opposite effects for the leaf blotch and glume blotch reactions suggest that different genetic mechanisms may be involved.

Segregation distortion: Utilizing simulated genotyping data to evaluate statistical methods

Segregation distortion is the phenomenon in which genotypes deviate from expected Mendelian ratios in the progeny of a cross between two varieties or species. There is not currently a widely used consensus for the appropriate statistical test, or more specifically the multiple testing correction procedure, used to detect segregation distortion for high-density single-nucleotide polymorphism (SNP) data. Here we examine the efficacy of various multiple testing procedures, including chi-square test with no correction for multiple testing, false-discovery rate correction and Bonferroni correction using an in-silico simulation of a biparental mapping population. We find that the false discovery rate correction best approximates the traditional p-value threshold of 0.05 for high-density marker data. We also utilize this simulation to test the effect of segregation distortion on the genetic mapping process, specifically on the formation of linkage groups during marker clustering. Only extreme segregation distortion was found to effect genetic mapping. In addition, we utilize replicate empirical mapping populations of wheat varieties Avalon and Cadenza to assess how often segregation distortion conforms to the same pattern between closely related wheat varieties.

Effect of marker segregation distortion on high density linkage map construction and QTL mapping in Soybean (Glycine max L.)

Marker segregation distortion is a natural phenomenon. Severely distorted markers are usually excluded in the construction of linkage maps. We investigated the effect of marker segregation distortion on linkage map construction and quantitative trait locus (QTL) mapping. A total of 519 recombinant inbred lines of soybean from orthogonal and reciprocal crosses between LSZZH and NN493-1 were genotyped by specific length amplified fragment markers and seed linoleic acid content was measured in three environments. As a result, twenty linkage groups were constructed with 11,846 markers, including 1513 (12.77%) significantly distorted markers, on 20 chromosomes, and the map length was 2475.86 cM with an average marker-interval of 0.21 cM. The inclusion of distorted markers in the analysis was shown to not only improve the grouping of the markers from the same chromosomes, and the consistency of linkage maps with genome, but also increase genome coverage by markers. Combining genotypic data from both orthogonal and reciprocal crosses decreased the proportion of distorted markers and then improved the quality of linkage maps. Validation of the linkage maps was confirmed by the high collinearity between positions of markers in the soybean reference genome and in linkage maps and by the high consistency of 24 QTL regions in this study compared with the previously reported QTLs and lipid metabolism related genes. Additionally, linkage maps that include distorted markers could add more information to the outputs from QTL mapping. These results provide important information for linkage mapping, gene cloning and marker-assisted selection in soybean.

Prospects of orphan crops in climate change

MAIN CONCLUSION

Orphan crops can contribute to building resilience of marginal cropping systems as a climate chnage adaptation strategy. Orphan crops play an important role in global food and nutrition security, and may have potential to contribute to sustainable food systems under climate change. Owing to reports of their potential under water scarcity, there is an argument to promote them to sustainably address challenges such as increasing drought and water scarcity, food and nutrition insecurity, environmental degradation, and employment creation under climate change. We conducted a scoping review using online databases to identify the prospects of orphan crops to contribute to (1) sustainable and healthy food systems, (2) genetic resources for future crop improvement, and (3) improving agricultural sustainability under climate change. The review found that, as a product of generations of landrace agriculture, several orphan crops are nutritious, resilient, and adapted to niche marginal agricultural environments. Including such orphan crops in the existing monocultural cropping systems could support more sustainable, nutritious, and diverse food systems in marginalised agricultural environments. Orphan crops also represent a broad gene pool for future crop improvement. The reduction in arable land due to climate change offers opportunities to expand the area under their production. Their suitability to marginal niche and low-input environments offers opportunities for low greenhouse gas (GHG) emissions from an agro-ecosystems, production, and processing perspective. This, together with their status as a sub-set of agro-biodiversity, offers opportunities to address socio-economic and environmental challenges under climate change. With research and development, and policy to support them, orphan crops could play an important role in climate-change adaptation, especially in the global south.

Plant Genomics: Unlocking the Genome of Wheat's Progenitor

The genome sequence and analysis of wheat's progenitor provides a roadmap to enhance genomics-assisted breeding and improvement of modern wheat varieties.

Genotypic, Developmental and Environmental Effects on the Rapidity of gs in Wheat: Impacts on Carbon Gain and Water-Use Efficiency

Stomata are the primary gatekeepers for CO2 uptake for photosynthesis and water loss via transpiration and therefore play a central role in crop performance. Although stomatal conductance (gs ) and assimilation rate (A) are often highly correlated, studies have demonstrated an uncoupling between A and gs that can result in sub-optimal physiological processes in dynamic light environments. Wheat (Triticum aestivum L.) is exposed to changes in irradiance due to leaf self-shading, moving clouds and shifting sun angle to which both A and gs respond. However, stomatal responses are generally an order of magnitude slower than photosynthetic responses, leading to non-synchronized A and gs responses that impact CO2 uptake and water use efficiency ( iWUE). Here we phenotyped a panel of eight wheat cultivars (estimated to capture 80% of the single nucleotide polymorphism variation in North-West European bread wheat) for differences in the speed of stomatal responses (to changes in light intensity) and photosynthetic performance at different stages of development. The impact of water stress and elevated [CO2] on stomatal kinetics was also examined in a selected cultivar. Significant genotypic variation was reported for the time constant for stomatal opening (Ki, P = 0.038) and the time to reach 95% steady state A (P = 0.045). Slow gs opening responses limited A by ∼10% and slow closure reduced iWUE, with these impacts found to be greatest in cultivars Soissons, Alchemy and Xi19. A decrease in stomatal rapidity (and thus an increase in the limitation of photosynthesis) (P < 0.001) was found during the post-anthesis stage compared to the early booting stage. Reduced water availability triggered stomatal closure and asymmetric stomatal opening and closing responses, while elevated atmospheric [CO2] conditions reduced the time for stomatal opening during a low to high light transition, thus suggesting a major environmental effect on dynamic stomatal kinetics. We discuss these findings in terms of exploiting various traits to develop ideotypes for specific environments, and suggest that intraspecific variation in the rapidity of stomatal responses could provide a potential unexploited breeding target to optimize the physiological responses of wheat to dynamic field conditions.

CropSNPdb: a database of SNP array data for Brassica crops and hexaploid bread wheat

Advances in sequencing technology have led to a rapid rise in the genomic data available for plants, driving new insights into the evolution, domestication and improvement of crops. Single nucleotide polymorphisms (SNPs) are a major component of crop genomic diversity, and are invaluable as genetic markers in research and breeding programs. High-throughput SNP arrays, or 'SNP chips', can generate reproducible sets of informative SNP markers and have been broadly adopted. Although there are many public repositories for sequencing data, which are routinely uploaded, there are no formal repositories for crop SNP array data. To make SNP array data more easily accessible, we have developed CropSNPdb (http://snpdb.appliedbioinformatics.com.au), a database for SNP array data produced by the Illumina Infinium™ hexaploid bread wheat (Triticum aestivum) 90K and Brassica 60K arrays. We currently host SNPs from datasets covering 526 Brassica lines and 309 bread wheat lines, and provide search, download and upload utilities for users. CropSNPdb provides a useful repository for these data, which can be applied for a range of genomics and molecular crop-breeding activities.

From markers to genome-based breeding in wheat

Recent technological advances in wheat genomics provide new opportunities to uncover genetic variation in traits of breeding interest and enable genome-based breeding to deliver wheat cultivars for the projected food requirements for 2050. There has been tremendous progress in development of whole-genome sequencing resources in wheat and its progenitor species during the last 5 years. High-throughput genotyping is now possible in wheat not only for routine gene introgression but also for high-density genome-wide genotyping. This is a major transition phase to enable genome-based breeding to achieve progressive genetic gains to parallel to projected wheat production demands. These advances have intrigued wheat researchers to practice less pursued analytical approaches which were not practiced due to the short history of genome sequence availability. Such approaches have been successful in gene discovery and breeding applications in other crops and animals for which genome sequences have been available for much longer. These strategies include, (i) environmental genome-wide association studies in wheat genetic resources stored in genbanks to identify genes for local adaptation by using agroclimatic traits as phenotypes, (ii) haplotype-based analyses to improve the statistical power and resolution of genomic selection and gene mapping experiments, (iii) new breeding strategies for genome-based prediction of heterosis patterns in wheat, and (iv) ultimate use of genomics information to develop more efficient and robust genome-wide genotyping platforms to precisely predict higher yield potential and stability with greater precision. Genome-based breeding has potential to achieve the ultimate objective of ensuring sustainable wheat production through developing high yielding, climate-resilient wheat cultivars with high nutritional quality.

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.

Dense genotyping-by-sequencing linkage maps of two Synthetic W7984×Opata reference populations provide insights into wheat structural diversity

Wheat (Triticum aestivum) genetic maps are a key enabling tool for genetic studies. We used genotyping-by-sequencing-(GBS) derived markers to map recombinant inbred line (RIL) and doubled haploid (DH) populations from crosses of W7984 by Opata, and used the maps to explore features of recombination control. The RIL and DH populations, SynOpRIL and SynOpDH, were composed of 906 and 92 individuals, respectively. Two high-density genetic linkage framework maps were constructed of 2,842 and 2,961 cM, harboring 3,634 and 6,580 markers, respectively. Using imputation, we added 43,013 and 86,042 markers to the SynOpRIL and SynOpDH maps. We observed preferential recombination in telomeric regions and reduced recombination in pericentromeric regions. Recombination rates varied between subgenomes, with the D genomes of the two populations exhibiting the highest recombination rates of 0.26-0.27 cM/Mb. QTL mapping identified two additive and three epistatic loci associated with crossover number. Additionally, we used published POPSEQ data from SynOpDH to explore the structural variation in W7984 and Opata. We found that chromosome 5AS is missing from W7984. We also found 2,332 variations larger than 100 kb. Structural variants were more abundant in distal regions, and overlapped 9,196 genes. The two maps provide a resource for trait mapping and genomic-assisted breeding.

Resistance gene cloning from a wild crop relative by sequence capture and association genetics

Disease resistance (R) genes from wild relatives could be used to engineer broad-spectrum resistance in domesticated crops. We combined association genetics with R gene enrichment sequencing (AgRenSeq) to exploit pan-genome variation in wild diploid wheat and rapidly clone four stem rust resistance genes. AgRenSeq enables R gene cloning in any crop that has a diverse germplasm panel.

Usefulness of a Multiparent Advanced Generation Intercross Population With a Greatly Reduced Mating Design for Genetic Studies in Winter Wheat

Multiparent advanced generation intercross (MAGIC) populations were recently developed to allow the high-resolution mapping of quantitative traits. We present a genetic linkage map of an elite but highly diverse eight-founder MAGIC population in common wheat (Triticum aestivum L.). Our MAGIC population is composed of 394 F6:8 recombinant inbred lines lacking significant signatures of population structure. The linkage map included 5435 SNP markers distributed over 2804 loci and spanning 5230 cM. The analysis of population parameters, including genetic structure, kinship, founder probabilities, and linkage disequilibrium and congruency to other maps indicated appropriate construction of both the population and the genetic map. It was shown that eight-founder MAGIC populations exhibit a greater number of loci and higher recombination rates, especially in the pericentromeric regions, compared to four-founder MAGIC, and biparental populations. In addition, our greatly simplified eight-parental MAGIC mating design with an additional eight-way intercross step was found to be equivalent to a MAGIC design with all 210 possible four-way crosses regarding the levels of missing founder assignments and the number of recombination events. Furthermore, the MAGIC population captured 71.7% of the allelic diversity available in the German wheat breeding gene pool. As a proof of principle, we demonstrated the application of the resource for quantitative trait loci mapping analyzing seedling resistance to powdery mildew. As wheat is a crop with many breeding objectives, this resource will allow scientists and breeders to carry out genetic studies for a wide range of breeder-relevant parameters in a single genetic background and reveal possible interactions between traits of economic importance.

The genetic architecture of genome-wide recombination rate variation in allopolyploid wheat revealed by nested association mapping

Recombination affects the fate of alleles in populations by imposing constraints on the reshuffling of genetic information. Understanding the genetic basis of these constraints is critical for manipulating the recombination process to improve the resolution of genetic mapping, and reducing the negative effects of linkage drag and deleterious genetic load in breeding. Using sequence-based genotyping of a wheat nested association mapping (NAM) population of 2,100 recombinant inbred lines created by crossing 29 diverse lines, we mapped QTL affecting the distribution and frequency of 102 000 crossovers (CO). Genome-wide recombination rate variation was mostly defined by rare alleles with small effects together explaining up to 48.6% of variation. Most QTL were additive and showed predominantly trans-acting effects. The QTL affecting the proximal COs also acted additively without increasing the frequency of distal COs. We showed that the regions with decreased recombination carry more single nucleotide polymorphisms (SNPs) with possible deleterious effects than the regions with a high recombination rate. Therefore, our study offers insights into the genetic basis of recombination rate variation in wheat and its effect on the distribution of deleterious SNPs across the genome. The identified trans-acting additive QTL can be utilized to manipulate CO frequency and distribution in the large polyploid wheat genome opening the possibility to improve the efficiency of gene pyramiding and reducing the deleterious genetic load in the low-recombining pericentromeric regions of chromosomes.

Adaptive selection of founder segments and epistatic control of plant height in the MAGIC winter wheat population WM-800

BACKGROUND

Multi-parent advanced generation intercross (MAGIC) populations are a newly established tool to dissect quantitative traits. We developed the high resolution MAGIC wheat population WM-800, consisting of 910 F4:6 lines derived from intercrossing eight recently released European winter wheat cultivars.

RESULTS

Genotyping WM-800 with 7849 SNPs revealed a low mean genetic similarity of 59.7% between MAGIC lines. WM-800 harbours distinct genomic regions exposed to segregation distortion. These are mainly located on chromosomes 2 to 6 of the wheat B genome where founder specific DNA segments were positively or negatively selected. This suggests adaptive selection of individual founder alleles during population development. The application of a genome-wide association study identified 14 quantitative trait loci (QTL) controlling plant height in WM-800, including the known semi-dwarf genes Rht-B1 and Rht-D1 and a potentially novel QTL on chromosome 5A. Additionally, epistatic effects controlled plant height. For example, two loci on chromosomes 2B and 7B gave rise to an additive epistatic effect of 13.7 cm.

CONCLUSION

The present study demonstrates that plant height in the MAGIC-WHEAT population WM-800 is mainly determined by large-effect QTL and di-genic epistatic interactions. As a proof of concept, our study confirms that WM-800 is a valuable tool to dissect the genetic architecture of important agronomic traits.

Meta-analysis of chromosome-scale crossover rate variation in eukaryotes and its significance to evolutionary genomics

Understanding the distribution of crossovers along chromosomes is crucial to evolutionary genomics because the crossover rate determines how strongly a genome region is influenced by natural selection on linked sites. Nevertheless, generalities in the chromosome-scale distribution of crossovers have not been investigated formally. We fill this gap by synthesizing joint information on genetic and physical maps across 62 animal, plant and fungal species. Our quantitative analysis reveals a strong and taxonomically widespread reduction of the crossover rate in the centre of chromosomes relative to their peripheries. We demonstrate that this pattern is poorly explained by the position of the centromere, but find that the magnitude of the relative reduction in the crossover rate in chromosome centres increases with chromosome length. That is, long chromosomes often display a dramatically low crossover rate in their centre, whereas short chromosomes exhibit a relatively homogeneous crossover rate. This observation is compatible with a model in which crossover is initiated from the chromosome tips, an idea with preliminary support from mechanistic investigations of meiotic recombination. Consequently, we show that organisms achieve a higher genome-wide crossover rate by evolving smaller chromosomes. Summarizing theory and providing empirical examples, we finally highlight that taxonomically widespread and systematic heterogeneity in crossover rate along chromosomes generates predictable broad-scale trends in genetic diversity and population differentiation by modifying the impact of natural selection among regions within a genome. We conclude by emphasizing that chromosome-scale heterogeneity in crossover rate should urgently be incorporated into analytical tools in evolutionary genomics, and in the interpretation of resulting patterns.

TaAAP6-3B, a regulator of grain protein content selected during wheat improvement

BACKGROUND

The content of grain protein (GPC) in cereals is an important part of total protein in human food. Exploring and utilizing new GPC genes is one of the most effective approaches for wheat quality breeding.

RESULTS

Three homoeologues of TaAAP6(-3A, 3B, 3D)were cloned by homology cloning from OsAAP6.Temporal and spatial expression analysis showed that TaAAP6homoeologues were preferentially expressed in developing grains, and TaAAP6-3B may play a major role in regulating GPC in wheat. Association analysis indicated thatTaAAP6-3B-I is significantly correlated with higher GPC than that of TaAAP6-3B-II for 115 wheat lines in all five environments. TaAAP6-3B-I, the favored allele of TaAAP6-3B, was preferentially expressed in preliminary developing grain stage. Two functional markers were developed to discriminate 197F2populations and the result showed that TaAAP6-3B-I (high-protein content) was completely dominant. Two cis-regulatory elements appear to be associated with high GPC were found in the 5'UTR of TaAAP6-3B-I.The change of the TaAAP6-3B locus types indicated that the gene was subjected to selection pressures during long process of artificial selection.

CONCLUSIONS

TaAAP6-3B is a regulator of GPC and its favored allele TaAAP6-3B-I exhibits an obvious potential application in wheat high-GPC breeding.

TEOSINTE BRANCHED1 Regulates Inflorescence Architecture and Development in Bread Wheat (Triticum aestivum)

The flowers of major cereals are arranged on reproductive branches known as spikelets, which group together to form an inflorescence. Diversity for inflorescence architecture has been exploited during domestication to increase crop yields, and genetic variation for this trait has potential to further boost grain production. Multiple genes that regulate inflorescence architecture have been identified by studying alleles that modify gene activity or dosage; however, little is known in wheat. Here, we show TEOSINTE BRANCHED1 (TB1) regulates inflorescence architecture in bread wheat (Triticum aestivum) by investigating lines that display a form of inflorescence branching known as "paired spikelets." We show that TB1 interacts with FLOWERING LOCUS T1 and that increased dosage of TB1 alters inflorescence architecture and growth rate in a process that includes reduced expression of meristem identity genes, with allelic diversity for TB1 found to associate genetically with paired spikelet development in modern cultivars. We propose TB1 coordinates formation of axillary spikelets during the vegetative to floral transition and that alleles known to modify dosage or function of TB1 could help increase wheat yields.

SLAF-based high-density genetic map construction and QTL mapping for major economic traits in sea urchin Strongylocentrotus intermedius

Sea urchin (Strongylocentrotus intermedius) has long been a model species for developmental and evolutionary research, but only a few studies have focused on gene mapping. Here, we reported a high-density genetic map containing 4,387 polymorphism specific-length amplified fragment (SLAF) markers spanning 21 linkage groups (LG) for sea urchin. Based on this genetic map and phenotyping data for eight economic traits, 33 potentially significant QTLs were detected on ten different LGs with explanations ranging from 9.90% to 46.30%, partly including 10 QTLs for test diameter, six QTLs for body weight and eight QTLs for Aristotle's lantern weight. Moreover, we found a QTL enrichment LG, LG15, gathering QTLs for test diameter, body weight, gonad weight, light orange-yellow color difference (≥E₁) and light yellow color difference (≥E₂). Among all QTLs, we genotyped four QTLs for test diameter, Aristotle's lantern weight and body weight using High Resolution Melting (HRM) technology. Finally, we used the verified SNP marker (detected using SLAF sequencing) to explore their marker-assisted selection (MAS) breeding application potential and found that SNP-29 associated tightly with body weight and that heterozygous genotype was a dominant genotype, indicating that SNP-29 was a promising marker for MAS.