Novel PHOTOPERIOD-1 gene variants associate with yield-related and root-angle traits in European bread wheat

KEY MESSAGE

PHOTOPERIOD-1 homoeologous gene copies play a pivotal role in regulation of flowering time in wheat. Here, we show that their influence also extends to spike and shoot architecture and even impacts root development. The sequence diversity of three homoeologous copies of the PHOTOPERIOD-1 gene in European winter wheat was analyzed by Oxford Nanopore amplicon-based multiplex sequencing and molecular markers in a panel of 194 cultivars representing breeding progress over the past 5 decades. A strong, consistent association with an average 8% increase in grain yield was observed for the PpdA1-Hap1 haplotype across multiple environments. This haplotype was found to be linked in 51% of cultivars to the 2NS/2AS translocation, originally introduced from Aegilops ventricosa, which leads to an overestimation of its effect. However, even in cultivars without the 2NS/2AS translocation, PpdA1-Hap1 was significantly associated with increased grain yield, kernel per spike and kernel per m2 under optimal growth conditions, conferring a 4% yield advantage compared to haplotype PpdA1-Hap4. In contrast to Ppd-B1 and Ppd-D1, the Ppd-A1 gene exhibits novel structural variations and a high number of SNPs, highlighting the evolutionary changes that have occurred in this region over the course of wheat breeding history. Additionally, cultivars carrying the photoperiod-insensitive Ppd-D1a allele not only exhibit earlier heading, but also deeper roots compared to those with photoperiod-sensitive alleles under German conditions. PCR and KASP assays have been developed that can be effectively employed in marker-assisted breeding programs to introduce these favorable haplotypes.

Use of continuous genotypes for genomic prediction in sugarcane

Genomic selection in sugarcane faces challenges due to limited genomic tools and high genomic complexity, particularly because of its high and variable ploidy. The classification of genotypes for single nucleotide polymorphisms (SNPs) becomes difficult due to the wide range of possible allele dosages. Previous genomic studies in sugarcane used pseudo-diploid genotyping, grouping all heterozygotes into a single class. In this study, we investigate the use of continuous genotypes as a proxy for allele-dosage in genomic prediction models. The hypothesis is that continuous genotypes could better reflect allele dosage at SNPs linked to mutations affecting target traits, resulting in phenotypic variation. The dataset included genotypes of 1318 clones at 58K SNP markers, with about 26K markers filtered using standard quality controls. Predictions for tonnes of cane per hectare (TCH), commercial cane sugar (CCS), and fiber content (Fiber) were made using parametric, non-parametric, and Bayesian methods. Continuous genotypes increased accuracy by 5%-7% for CCS and Fiber. The pseudo-diploid parametrization performed better for TCH. Reproducing kernel Hilbert spaces model with Gaussian kernel and AK4 (arc-cosine kernel with hidden layer 4) kernel outperformed other methods for TCH and CCS, suggesting that non-additive effects might influence these traits. The prevalence of low-dosage markers in the study may have limited the benefits of approximating allele-dosage information with continuous genotypes in genomic prediction models. Continuous genotypes simplify genomic prediction in polyploid crops, allowing additional markers to be used without adhering to pseudo-diploid inheritance. The approach can particularly benefit high ploidy species or emerging crops with unknown ploidy.

Optimising clonal performance in sugarcane: leveraging non-additive effects via mate-allocation strategies

Mate-allocation strategies in breeding programs can improve progeny performance by harnessing non-additive genetic effects. These approaches prioritise predicted progeny merit over parental breeding value, making them particularly appealing for clonally propagated crops such as sugarcane. We conducted a comparative analysis of mate-allocation strategies, exploring utilising non-additive and heterozygosity effects to maximise clonal performance with schemes that solely consider additive effects to optimise breeding value. Using phenotypic and genotypic data from a population of 2,909 clones evaluated in final assessment trials of Australian sugarcane breeding programs, we focused on three important traits: tonnes of cane per hectare (TCH), commercial cane sugar (CCS), and Fibre. By simulating families from all possible crosses (1,225) with 50 progenies each, we predicted the breeding and clonal values of progeny using two models: GBLUP (considering additive effects only) and extended-GBLUP (incorporating additive, non-additive, and heterozygosity effects). Integer linear programming was used to identify the optimal mate-allocation among selected parents. Compared to breeding value-based approaches, mate-allocation strategies based on clonal performance yielded substantial improvements, with predicted progeny values increasing by 57% for TCH, 12% for CCS, and 16% for fibre. Our simulation study highlights the effectiveness of mate-allocation approaches that exploit non-additive and heterozygosity effects, resulting in superior clonal performance. However, there was a notable decline in additive gain, particularly for TCH, likely due to significant epistatic effects. When selecting crosses based on clonal performance for TCH, the inbreeding coefficient of progeny was significantly lower compared to random mating, underscoring the advantages of leveraging non-additive and heterozygosity effects in mitigating inbreeding depression. Thus, mate-allocation strategies are recommended in clonally propagated crops to enhance clonal performance and reduce the negative impacts of inbreeding.

Dissecting the genetic architecture of root-related traits in a grafted wild Vitis berlandieri population for grapevine rootstock breeding

In woody perennial plants, quantitative genetics and association studies remain scarce for root-related traits, due to the time required to obtain mature plants and the complexity of phenotyping. In grapevine, a grafted cultivated plant, most of the rootstocks used are hybrids between American Vitis species (V. rupestris, V. riparia, and V. berlandieri). In this study, we used a wild population of an American Vitis species (V. berlandieri) to analyze the genetic architecture of the root-related traits of rootstocks in a grafted context. We studied a population consisting of 211 genotypes, with one to five replicates each (n = 846 individuals), plus four commercial rootstocks as control genotypes (110R, 5BB, Börner, and SO4). After two independent years of experimentation, the best linear unbiased estimates method revealed root-related traits with a moderate-to-high heritability (0.36-0.82) and coefficient of genetic variation (0.15-0.45). A genome-wide association study was performed with the BLINK model, leading to the detection of 11 QTL associated with four root-related traits (one QTL was associated with the total number of roots, four were associated with the number of small roots (< 1 mm in diameter), two were associated with the number of medium-sized roots (1 mm < diameter < 2 mm), and four were associated with mean diameter) accounting for up to 25.1% of the variance. Three genotypes were found to have better root-related trait performances than the commercial rootstocks and therefore constitute possible new candidates for use in grapevine rootstock breeding programs.

Haplotype blocks for genomic prediction: a comparative evaluation in multiple crop datasets

In modern plant breeding, genomic selection is becoming the gold standard for selection of superior genotypes. The basis for genomic prediction models is a set of phenotyped lines along with their genotypic profile. With high marker density and linkage disequilibrium (LD) between markers, genotype data in breeding populations tends to exhibit considerable redundancy. Therefore, interest is growing in the use of haplotype blocks to overcome redundancy by summarizing co-inherited features. Moreover, haplotype blocks can help to capture local epistasis caused by interacting loci. Here, we compared genomic prediction methods that either used single SNPs or haplotype blocks with regards to their prediction accuracy for important traits in crop datasets. We used four published datasets from canola, maize, wheat and soybean. Different approaches to construct haplotype blocks were compared, including blocks based on LD, physical distance, number of adjacent markers and the algorithms implemented in the software "Haploview" and "HaploBlocker". The tested prediction methods included Genomic Best Linear Unbiased Prediction (GBLUP), Extended GBLUP to account for additive by additive epistasis (EGBLUP), Bayesian LASSO and Reproducing Kernel Hilbert Space (RKHS) regression. We found improved prediction accuracy in some traits when using haplotype blocks compared to SNP-based predictions, however the magnitude of improvement was very trait- and model-specific. Especially in settings with low marker density, haplotype blocks can improve genomic prediction accuracy. In most cases, physically large haplotype blocks yielded a strong decrease in prediction accuracy. Especially when prediction accuracy varies greatly across different prediction models, prediction based on haplotype blocks can improve prediction accuracy of underperforming models. However, there is no "best" method to build haplotype blocks, since prediction accuracy varied considerably across methods and traits. Hence, criteria used to define haplotype blocks should not be viewed as fixed biological parameters, but rather as hyperparameters that need to be adjusted for every dataset.

Towards grapevine root architectural models to adapt viticulture to drought

To sustainably adapt viticultural production to drought, the planting of rootstock genotypes adapted to a changing climate is a promising means. Rootstocks contribute to the regulation of scion vigor and water consumption, modulate scion phenological development and determine resource availability by root system architecture development. There is, however, a lack of knowledge on spatio-temporal root system development of rootstock genotypes and its interactions with environment and management that prevents efficient knowledge transfer into practice. Hence, winegrowers take only limited advantage of the large variability of existing rootstock genotypes. Models of vineyard water balance combined with root architectural models, using both static and dynamic representations of the root system, seem promising tools to match rootstock genotypes to frequently occurring future drought stress scenarios and address scientific knowledge gaps. In this perspective, we discuss how current developments in vineyard water balance modeling may provide the background for a better understanding of the interplay of rootstock genotypes, environment and management. We argue that root architecture traits are key drivers of this interplay, but our knowledge on rootstock architectures in the field remains limited both qualitatively and quantitatively. We propose phenotyping methods to help close current knowledge gaps and discuss approaches to integrate phenotyping data into different models to advance our understanding of rootstock x environment x management interactions and predict rootstock genotype performance in a changing climate. This could also provide a valuable basis for optimizing breeding efforts to develop new grapevine rootstock cultivars with optimal trait configurations for future growing conditions.

genomicSimulation: fast R functions for stochastic simulation of breeding programs

Simulation tools are key to designing and optimizing breeding programs that are multiyear, high-effort endeavors. Tools that operate on real genotypes and integrate easily with other analysis software can guide users toward crossing decisions that best balance genetic gains and genetic diversity required to maintain gains in the future. Here, we present genomicSimulation, a fast and flexible tool for the stochastic simulation of crossing and selection based on real genotypes. It is fully written in C for high execution speeds, has minimal dependencies, and is available as an R package for the integration with R’s broad range of analysis and visualization tools. Comparisons of a simulated recreation of a breeding program to a real data set demonstrate the simulated offspring from the tool correctly show key population features, such as genomic relationships and approximate linkage disequilibrium patterns. Both versions of genomicSimulation are freely available on GitHub: The R package version at https://github.com/vllrs/genomicSimulation/ and the C library version at https://github.com/vllrs/genomicSimulationC/.

Mining the Vavilov wheat diversity panel for new sources of adult plant resistance to stripe rust

Multi-year evaluation of the Vavilov wheat diversity panel identified new sources of adult plant resistance to stripe rust. Genome-wide association studies revealed the key genomic regions influencing resistance, including seven novel loci. Wheat stripe rust (YR) caused by Puccinia striiformis f. sp. tritici (Pst) poses a significant threat to global food security. Resistance genes commonly found in many wheat varieties have been rendered ineffective due to the rapid evolution of the pathogen. To identify novel sources of adult plant resistance (APR), 292 accessions from the N.I. Vavilov Institute of Plant Genetic Resources, Saint Petersburg, Russia, were screened for known APR genes (i.e. Yr18, Yr29, Yr46, Yr33, Yr39 and Yr59) using linked polymerase chain reaction (PCR) molecular markers. Accessions were evaluated against Pst (pathotype 134 E16 A + Yr17 + Yr27) at seedling and adult plant stages across multiple years (2014, 2015 and 2016) in Australia. Phenotypic analyses identified 132 lines that potentially carry novel sources of APR to YR. Genome-wide association studies (GWAS) identified 68 significant marker-trait associations (P < 0.001) for YR resistance, representing 47 independent quantitative trait loci (QTL) regions. Fourteen genomic regions overlapped with previously reported Yr genes, including Yr29, Yr56, Yr5, Yr43, Yr57, Yr30, Yr46, Yr47, Yr35, Yr36, Yrxy1, Yr59, Yr52 and YrYL. In total, seven QTL (positioned on chromosomes 1D, 2A, 3A, 3D, 5D, 7B and 7D) did not collocate with previously reported genes or QTL, indicating the presence of promising novel resistance factors. Overall, the Vavilov diversity panel provides a rich source of new alleles which could be used to broaden the genetic bases of YR resistance in modern wheat varieties.

Environmental variation effects fertility in tropical beef cattle

The northern Australia beef cattle industry operates in harsh environmental conditions which consistently suppress female fertility. To better understand the environmental effect on cattle raised extensively in northern Australia, new environmental descriptors were defined for 54 commercial herds located across the region. Three fertility traits, based on the presence of a corpus luteum at 600 d of age, indicating puberty, (CL Presence, n = 25,176), heifer pregnancy (n = 20,989) and first lactation pregnancy (n = 10,072) were recorded. Temperature, humidity, and rainfall were obtained from publicly available data based on herd location. Being pubertal at 600 d (i.e. CL Presence) increased the likelihood of success at heifer pregnancy and first lactation pregnancy (P < 0.05), underscoring the importance of early puberty in reproductive success. A temperature humidity index (THI) of 65–70 had a significant (P < 0.05) negative effect on first lactation pregnancy rate, heifer pregnancy and puberty at 600 d of age. Area under the curve of daily THI was significant (P < 0.05) and reduced the likelihood of pregnancy at first lactation and puberty at 600 days. Deviation from long-term average rainfall was not significant (P < 0.05) for any trait. Average daily weight gain had a significant and positive relationship (P < 0.05) for heifer and first lactation pregnancy. The results indicate that chronic or cumulative heat load is more determinantal to reproductive performance than acute heat stress. The reason for the lack of a clear relationship between acute heat stress and reproductive performance is unclear but may be partially explained by peak THI and peak nutrition coinciding at the same time. Sufficient evidence was found to justify the use of average daily weight gain and chronic heat load as descriptors to define an environmental gradient.

A toolkit to rapidly modify root systems through single plant selection

BACKGROUND

The incorporation of root traits into elite germplasm is typically a slow process. Thus, innovative approaches are required to accelerate research and pre-breeding programs targeting root traits to improve yield stability in different environments and soil types. Marker-assisted selection (MAS) can help to speed up the process by selecting key genes or quantitative trait loci (QTL) associated with root traits. However, this approach is limited due to the complex genetic control of root traits and the limited number of well-characterised large effect QTL. Coupling MAS with phenotyping could increase the reliability of selection. Here we present a useful framework to rapidly modify root traits in elite germplasm. In this wheat exemplar, a single plant selection (SPS) approach combined three main elements: phenotypic selection (in this case for seminal root angle); MAS using KASP markers (targeting a root biomass QTL); and speed breeding to accelerate each cycle.

RESULTS

To develop a SPS approach that integrates non-destructive screening for seminal root angle and root biomass, two initial experiments were conducted. Firstly, we demonstrated that transplanting wheat seedlings from clear pots (for seminal root angle assessment) into sand pots (for root biomass assessment) did not impact the ability to differentiate genotypes with high and low root biomass. Secondly, we demonstrated that visual scores for root biomass were correlated with root dry weight (r = 0.72), indicating that single plants could be evaluated for root biomass in a non-destructive manner. To highlight the potential of the approach, we applied SPS in a backcrossing program which integrated MAS and speed breeding for the purpose of rapidly modifying the root system of elite bread wheat line Borlaug100. Bi-directional selection for root angle in segregating generations successfully shifted the mean root angle by 30° in the subsequent generation (P ≤ 0.05). Within 18 months, BC₂F₄:F₅ introgression lines were developed that displayed a full range of root configurations, while retaining similar above-ground traits to the recurrent parent. Notably, the seminal root angle displayed by introgression lines varied more than 30° compared to the recurrent parent, resulting in lines with both narrow and wide root angles, and high and low root biomass phenotypes.

CONCLUSION

The SPS approach enables researchers and plant breeders to rapidly manipulate root traits of future crop varieties, which could help improve productivity in the face of increasing environmental fluctuations. The newly developed elite wheat lines with modified root traits provide valuable materials to study the value of different root systems to support yield in different environments and soil types.

A linkage disequilibrium-based approach to position unmapped SNPs in crop species

BACKGROUND

High-density SNP arrays are now available for a wide range of crop species. Despite the development of many tools for generating genetic maps, the genome position of many SNPs from these arrays is unknown. Here we propose a linkage disequilibrium (LD)-based algorithm to allocate unassigned SNPs to chromosome regions from sparse genetic maps. This algorithm was tested on sugarcane, wheat, and barley data sets. We calculated the algorithm's efficiency by masking SNPs with known locations, then assigning their position to the map with the algorithm, and finally comparing the assigned and true positions.

RESULTS

In the 20-fold cross-validation, the mean proportion of masked mapped SNPs that were placed by the algorithm to a chromosome was 89.53, 94.25, and 97.23% for sugarcane, wheat, and barley, respectively. Of the markers that were placed in the genome, 98.73, 96.45 and 98.53% of the SNPs were positioned on the correct chromosome. The mean correlations between known and new estimated SNP positions were 0.97, 0.98, and 0.97 for sugarcane, wheat, and barley. The LD-based algorithm was used to assign 5920 out of 21,251 unpositioned markers to the current Q208 sugarcane genetic map, representing the highest density genetic map for this species to date.

CONCLUSIONS

Our LD-based approach can be used to accurately assign unpositioned SNPs to existing genetic maps, improving genome-wide association studies and genomic prediction in crop species with fragmented and incomplete genome assemblies. This approach will facilitate genomic-assisted breeding for many orphan crops that lack genetic and genomic resources.

Genetic characterization of adult-plant resistance to tan spot (syn, yellow spot) in wheat

QTL mapping identified key genomic regions associated with adult-plant resistance to tan spot, which are effective even in the presence of the sensitivity gene Tsn1, thus serving as a new genetic solution to develop disease-resistant wheat cultivars. Improving resistance to tan spot (Pyrenophora tritici-repentis; Ptr) in wheat by eliminating race-specific susceptibility genes is a common breeding approach worldwide. The potential to exploit variation in quantitative forms of resistance, such as adult-plant resistance (APR), offers an alternative approach that could lead to broad-spectrum protection. We previously identified wheat landraces in the Vavilov diversity panel that exhibited high levels of APR despite carrying the sensitivity gene Tsn1. In this study, we characterised the genetic control of APR by developing a recombinant inbred line population fixed for Tsn1, but segregating for the APR trait. Linkage mapping using DArTseq markers and disease response phenotypes identified a QTL associated with APR to Ptr race 1 (producing Ptr ToxA- and Ptr ToxC) on chromosome 2B (Qts.313-2B), which was consistently detected in multiple adult-plant experiments. Additional loci were also detected on chromosomes 2A, 3D, 5A, 5D, 6A, 6B and 7A at the seedling stage, and on chromosomes 1A and 5B at the adult stage. We demonstrate that Qts.313-2B can be combined with other adult-plant QTL (i.e. Qts.313-1A and Qts.313-5B) to strengthen resistance levels. The APR QTL reported in this study provide a new genetic solution to tan spot in Australia and could be deployed in wheat cultivars, even in the presence of Tsn1, to decrease production losses and reduce the application of fungicides.

Improved genomic prediction of clonal performance in sugarcane by exploiting non-additive genetic effects

Non-additive genetic effects seem to play a substantial role in the expression of complex traits in sugarcane. Including non-additive effects in genomic prediction models significantly improves the prediction accuracy of clonal performance. In the recent decade, genetic progress has been slow in sugarcane. One reason might be that non-additive genetic effects contribute substantially to complex traits. Dense marker information provides the opportunity to exploit non-additive effects in genomic prediction. In this study, a series of genomic best linear unbiased prediction (GBLUP) models that account for additive and non-additive effects were assessed to improve the accuracy of clonal prediction. The reproducible kernel Hilbert space model, which captures non-additive genetic effects, was also tested. The models were compared using 3,006 genotyped elite clones measured for cane per hectare (TCH), commercial cane sugar (CCS), and Fibre content. Three forward prediction scenarios were considered to investigate the robustness of genomic prediction. By using a pseudo-diploid parameterization, we found significant non-additive effects that accounted for almost two-thirds of the total genetic variance for TCH. Average heterozygosity also had a major impact on TCH, indicating that directional dominance may be an important source of phenotypic variation for this trait. The extended-GBLUP model improved the prediction accuracies by at least 17% for TCH, but no improvement was observed for CCS and Fibre. Our results imply that non-additive genetic variance is important for complex traits in sugarcane, although further work is required to better understand the variance component partitioning in a highly polyploid context. Genomics-based breeding will likely benefit from exploiting non-additive genetic effects, especially in designing crossing schemes. These findings can help to improve clonal prediction, enabling a more accurate identification of variety candidates for the sugarcane industry.

Perspectives on Applications of Hierarchical Gene-To-Phenotype (G2P) Maps to Capture Non-stationary Effects of Alleles in Genomic Prediction

Genomic prediction of complex traits across environments, breeding cycles, and populations remains a challenge for plant breeding. A potential explanation for this is that underlying non-additive genetic (GxG) and genotype-by-environment (GxE) interactions generate allele substitution effects that are non-stationary across different contexts. Such non-stationary effects of alleles are either ignored or assumed to be implicitly captured by most gene-to-phenotype (G2P) maps used in genomic prediction. The implicit capture of non-stationary effects of alleles requires the G2P map to be re-estimated across different contexts. We discuss the development and application of hierarchical G2P maps that explicitly capture non-stationary effects of alleles and have successfully increased short-term prediction accuracy in plant breeding. These hierarchical G2P maps achieve increases in prediction accuracy by allowing intermediate processes such as other traits and environmental factors and their interactions to contribute to complex trait variation. However, long-term prediction remains a challenge. The plant breeding community should undertake complementary simulation and empirical experiments to interrogate various hierarchical G2P maps that connect GxG and GxE interactions simultaneously. The existing genetic correlation framework can be used to assess the magnitude of non-stationary effects of alleles and the predictive ability of these hierarchical G2P maps in long-term, multi-context genomic predictions of complex traits in plant breeding.

Wheat root systems as a breeding target for climate resilience

In the coming decades, larger genetic gains in yield will be necessary to meet projected demand, and this must be achieved despite the destabilizing impacts of climate change on crop production. The root systems of crops capture the water and nutrients needed to support crop growth, and improved root systems tailored to the challenges of specific agricultural environments could improve climate resiliency. Each component of root initiation, growth and development is controlled genetically and responds to the environment, which translates to a complex quantitative system to navigate for the breeder, but also a world of opportunity given the right tools. In this review, we argue that it is important to know more about the 'hidden half' of crop plants and hypothesize that crop improvement could be further enhanced using approaches that directly target selection for root system architecture. To explore these issues, we focus predominantly on bread wheat (Triticum aestivum L.), a staple crop that plays a major role in underpinning global food security. We review the tools available for root phenotyping under controlled and field conditions and the use of these platforms alongside modern genetics and genomics resources to dissect the genetic architecture controlling the wheat root system. To contextualize these advances for applied wheat breeding, we explore questions surrounding which root system architectures should be selected for, which agricultural environments and genetic trait configurations of breeding populations are these best suited to, and how might direct selection for these root ideotypes be implemented in practice.

Tackling G × E × M interactions to close on-farm yield-gaps: creating novel pathways for crop improvement by predicting contributions of genetics and management to crop productivity

KEY MESSAGE

Climate change and Genotype-by-Environment-by-Management interactions together challenge our strategies for crop improvement. Research to advance prediction methods for breeding and agronomy is opening new opportunities to tackle these challenges and overcome on-farm crop productivity yield-gaps through design of responsive crop improvement strategies. Genotype-by-Environment-by-Management (G × E × M) interactions underpin many aspects of crop productivity. An important question for crop improvement is "How can breeders and agronomists effectively explore the diverse opportunities within the high dimensionality of the complex G × E × M factorial to achieve sustainable improvements in crop productivity?" Whenever G × E × M interactions make important contributions to attainment of crop productivity, we should consider how to design crop improvement strategies that can explore the potential space of G × E × M possibilities, reveal the interesting Genotype-Management (G-M) technology opportunities for the Target Population of Environments (TPE), and enable the practical exploitation of the associated improved levels of crop productivity under on-farm conditions. Climate change adds additional layers of complexity and uncertainty to this challenge, by introducing directional changes in the environmental dimension of the G × E × M factorial. These directional changes have the potential to create further conditional changes in the contributions of the genetic and management dimensions to future crop productivity. Therefore, in the presence of G × E × M interactions and climate change, the challenge for both breeders and agronomists is to co-design new G-M technologies for a non-stationary TPE. Understanding these conditional changes in crop productivity through the relevant sciences for each dimension, Genotype, Environment, and Management, creates opportunities to predict novel G-M technology combinations suitable to achieve sustainable crop productivity and global food security targets for the likely climate change scenarios. Here we consider critical foundations required for any prediction framework that aims to move us from the current unprepared state of describing G × E × M outcomes to a future responsive state equipped to predict the crop productivity consequences of G-M technology combinations for the range of environmental conditions expected for a complex, non-stationary TPE under the influences of climate change.

Strategies and considerations for implementing genomic selection to improve traits with additive and non-additive genetic architectures in sugarcane breeding

Simulations highlight the potential of genomic selection to substantially increase genetic gain for complex traits in sugarcane. The success rate depends on the trait genetic architecture and the implementation strategy. Genomic selection (GS) has the potential to increase the rate of genetic gain in sugarcane beyond the levels achieved by conventional phenotypic selection (PS). To assess different implementation strategies, we simulated two different GS-based breeding strategies and compared genetic gain and genetic variance over five breeding cycles to standard PS. GS scheme 1 followed similar routines like conventional PS but included three rapid recurrent genomic selection (RRGS) steps. GS scheme 2 also included three RRGS steps but did not include a progeny assessment stage and therefore differed more fundamentally from PS. Under an additive trait model, both simulated GS schemes achieved annual genetic gains of 2.6-2.7% which were 1.9 times higher compared to standard phenotypic selection (1.4%). For a complex non-additive trait model, the expected annual rates of genetic gain were lower for all breeding schemes; however, the rates for the GS schemes (1.5-1.6%) were still greater than PS (1.1%). Investigating cost-benefit ratios with regard to numbers of genotyped clones showed that substantial benefits could be achieved when only 1500 clones were genotyped per 10-year breeding cycle for the additive genetic model. Our results show that under a complex non-additive genetic model, the success rate of GS depends on the implementation strategy, the number of genotyped clones and the stage of the breeding program, likely reflecting how changes in QTL allele frequencies change additive genetic variance and therefore the efficiency of selection. These results are encouraging and motivate further work to facilitate the adoption of GS in sugarcane breeding.

Regional association analysis coupled with transcriptome analyses reveal candidate genes affecting seed oil accumulation in Brassica napus

Regional association analysis of 50 re-sequenced Chinese semi-winter rapeseed accessions in combination with co-expression analysis reveal candidate genes affecting oil accumulation in Brassica napus. One of the breeding goals in rapeseed production is to enhance the seed oil content to cater to the increased demand for vegetable oils due to a growing global population. To investigate the genetic basis of variation in seed oil content, we used 60 K Brassica Infinium SNP array along with phenotype data of 203 Chinese semi-winter rapeseed accessions to perform a genome-wide analysis of haplotype blocks associated with the oil content. Nine haplotype regions harbouring lipid synthesis/transport-, carbohydrate metabolism- and photosynthesis-related genes were identified as significantly associated with the oil content and were mapped to chromosomes A02, A04, A05, A07, C03, C04, C05, C08 and C09, respectively. Regional association analysis of 50 re-sequenced Chinese semi-winter rapeseed accessions combined with transcriptome datasets from 13 accessions was further performed on these nine haplotype regions. This revealed natural variation in the BnTGD3-A02 and BnSSE1-A05 gene regions correlated with the phenotypic variation of the oil content within the A02 and A04 chromosome haplotype regions, respectively. Moreover, co-expression network analysis revealed that BnTGD3-A02 and BnSSE1-A05 were directly linked with fatty acid beta-oxidation-related gene BnKAT2-C04, thus forming a molecular network involved in the potential regulation of seed oil accumulation. The results of this study could be used to combine favourable haplotype alleles for further improvement of the seed oil content in rapeseed.

Accuracy of genomic prediction of complex traits in sugarcane

Complex traits in sugarcane can be accurately predicted using genome-wide DNA markers. Genomic single-step prediction is an attractive method for genomic selection in commercial breeding programs. Sugarcane breeding programs have achieved up to 1% genetic gain in key traits such as tonnes of cane per hectare (TCH), commercial cane sugar (CCS) and Fibre content over the past decades. Here, we assess the potential of genomic selection to increase the rate of genetic gain for these traits by deriving genomic estimated breeding values (GEBVs) from a reference population of 3984 clones genotyped for 26 K SNP. We evaluated the three different genomic prediction approaches GBLUP, genomic single step (GenomicSS), and BayesR. GenomicSS combining pedigree and SNP information from historic and recent breeding programs achieved the most accurate predictions for most traits (0.3-0.44). This method is attractive for routine genetic evaluation because it requires relatively little modification to the existing evaluation and results in breeding value estimates for all individuals, not only those genotyped. Adding information from early-stage trials added up to 5% accuracy for CCS and Fibre, but 0% for TCH, reflecting the importance of competition effects for TCH. These GEBV accuracies are sufficiently high that, combined with the right breeding strategy, a doubling of the rate of genetic gain could be achieved. We also assessed the flowering traits days to flowering, gender and pollen viability and found high heritabilities of 0.57, 0.78 and 0.72, respectively. The GEBV accuracies indicated that genomic selection could be used to improve these traits. This could open new avenues for breeders to manage their breeding programs, for example, by synchronising flowering time and selecting males with high pollen viability.

Need for speed: manipulating plant growth to accelerate breeding cycles

To develop more productive and resilient crops that are capable of feeding 10 billion people by 2050, we must accelerate the rate of genetic improvement in plant breeding programs. Speed breeding manipulates the growing environment by regulating light and temperature for the purpose of rapid generation advance. Protocols are now available for a range of short-day and long-day species and the approach is highly compatible with other cutting-edge breeding tools such as genomic selection. Here, we highlight how speed breeding hijacks biological processes for applied plant breeding outcomes and provide a case study examining wheat growth and development under speed breeding conditions. The establishment of speed breeding facilities worldwide is expected to provide benefits for capacity building, discovery research, pre-breeding, and plant breeding to accelerate the development of productive and robust crops.

Overcoming polyploidy pitfalls: a user guide for effective SNP conversion into KASP markers in wheat

Conversion of SNP chip assays into locus-specific KASP markers requires adapted strategies in polyploid species with high genome homeology. Procedures are exemplified by QTL-associated SNPs in hexaploid wheat. Kompetitive allele-specific PCR (KASP) markers are commonly used in marker-assisted commercial plant breeding due to their cost-effectiveness and throughput for high sample volumes. However, conversion of trait-linked SNP markers from array-based SNP detection technologies into KASP markers is particularly challenging in polyploid crop species, due to the presence of highly similar homeologous and paralogous genome sequences. We evaluated strategies and identified key requirements for successful conversion of Illumina Infinium assays from the wheat 90 K SNP array into robust locus-specific KASP markers. Numerous examples showed that commonly used software for semiautomated KASP primer design frequently fails to achieve locus-specificity of KASP assays in wheat. Instead, alignment of SNP probes with multiple reference genomes and Sanger sequencing of relevant genotypes, followed by visual KASP primer placement, was critical for locus-specificity. To identify KASP assays resulting in false calling of heterozygous individuals, validation of KASP assays using extended reference genotype sets including heterozygous genotypes is strongly advised for polyploid crop species. Applying this strategy, we developed highly reproducible, stable KASP assays that are predictive for root biomass QTL haplotypes from highly homoeologous wheat chromosome regions. Due to their locus-specificity, these assays predicted root biomass considerably better than the original trait-associated markers from the Illumina array.

Adaptive Traits to Improve Durum Wheat Yield in Drought and Crown Rot Environments

Durum wheat (Triticum turgidum L. ssp. durum) production can experience significant yield losses due to crown rot (CR) disease. Losses are usually exacerbated when disease infection coincides with terminal drought. Durum wheat is very susceptible to CR, and resistant germplasm is not currently available in elite breeding pools. We hypothesize that deploying physiological traits for drought adaptation, such as optimal root system architecture to reduce water stress, might minimize losses due to CR infection. This study evaluated a subset of lines from a nested association mapping population for stay-green traits, CR incidence and yield in field experiments as well as root traits under controlled conditions. Weekly measurements of normalized difference vegetative index (NDVI) in the field were used to model canopy senescence and to determine stay-green traits for each genotype. Genome-wide association studies using DArTseq molecular markers identified quantitative trait loci (QTLs) on chromosome 6B (qCR-6B) associated with CR tolerance and stay-green. We explored the value of qCR-6B and a major QTL for root angle QTL qSRA-6A using yield datasets from six rainfed environments, including two environments with high CR disease pressure. In the absence of CR, the favorable allele for qSRA-6A provided an average yield advantage of 0.57 t·ha⁻¹, whereas in the presence of CR, the combination of favorable alleles for both qSRA-6A and qCR-6B resulted in a yield advantage of 0.90 t·ha⁻¹. Results of this study highlight the value of combining above- and belowground physiological traits to enhance yield potential. We anticipate that these insights will assist breeders to design improved durum varieties that mitigate production losses due to water deficit and CR.

GWAS and co-expression network combination uncovers multigenes with close linkage effects on the oleic acid content accumulation in Brassica napus

Background

Strong artificial and natural selection causes the formation of highly conserved haplotypes that harbor agronomically important genes. GWAS combination with haplotype analysis has evolved as an effective method to dissect the genetic architecture of complex traits in crop species.

Results

We used the 60 K Brassica Infinium SNP array to perform a genome-wide analysis of haplotype blocks associated with oleic acid (C18:1) in rapeseed. Six haplotype regions were identified as significantly associated with oleic acid (C18:1) that mapped to chromosomes A02, A07, A08, C01, C02, and C03. Additionally, whole-genome sequencing of 50 rapeseed accessions revealed three genes (BnmtACP2-A02, BnABCI13-A02 and BnECI1-A02) in the A02 chromosome haplotype region and two genes (BnFAD8-C02 and BnSDP1-C02) in the C02 chromosome haplotype region that were closely linked to oleic acid content phenotypic variation. Moreover, the co-expression network analysis uncovered candidate genes from these two different haplotype regions with potential regulatory interrelationships with oleic acid content accumulation.

Conclusions

Our results suggest that several candidate genes are closely linked, which provides us with an opportunity to develop functional haplotype markers for the improvement of the oleic acid content in rapeseed.

Revolutionizing agriculture with synthetic biology

Synthetic biology is here to stay and will transform agriculture if given the chance. The huge challenges facing food, fuel and chemical production make it vital to give synthetic biology that chance-notwithstanding the shifts in mindset, training and infrastructure investment this demands. Here, we assess opportunities for agricultural synthetic biology and ways to remove barriers to their realization.

High-resolution mapping of rachis nodes per rachis, a critical determinant of grain yield components in wheat

Exploring large genomic data sets based on the latest reference genome assembly identifies the rice ortholog APO1 as a key candidate gene for number of rachis nodes per spike in wheat. Increasing grain yield in wheat is a key breeding objective worldwide. Several component traits contribute to grain yield with spike attributes being among the most important. In this study, we performed a genome-wide association analysis for 12 grain yield and component traits measured in field trials with contrasting agrochemical input levels in a panel of 220 hexaploid winter wheats. A highly significant, environmentally consistent QTL was detected for number of rachis nodes per rachis (NRN) on chromosome 7AL. The five most significant SNPs formed a strong linkage disequilibrium (LD) block and tagged a 2.23 Mb region. Using pairwise LD for exome SNPs located across this interval in a large worldwide hexaploid wheat collection, we reduced the genomic region for NRN to a 258 Kb interval containing four of the original SNP and six high-confidence genes. The ortholog of one (TraesCS7A01G481600) of these genes in rice was ABBERANT PANICLE ORGANIZATION1 (APO1), which is known to have significant effects on panicle attributes. The APO1 ortholog was the best candidate for NRN and was associated with a 115 bp promoter deletion and two amino acid (C47F and D384 N) changes. Using a large worldwide collection of tetraploid and hexaploid wheat, we found 12 haplotypes for the NRN QTL and evidence for positive enrichment of two haplotypes in modern germplasm. Comparison of five QTL haplotypes in Australian yield trials revealed their relative, context-dependent contribution to grain yield. Our study provides diagnostic SNPs and value propositions to support deployment of the NRN trait in wheat breeding.

Breeding improves wheat productivity under contrasting agrochemical input levels

The world cropping area for wheat exceeds that of any other crop, and high grain yields in intensive wheat cropping systems are essential for global food security. Breeding has raised yields dramatically in high-input production systems; however, selection under optimal growth conditions is widely believed to diminish the adaptive capacity of cultivars to less optimal cropping environments. Here, we demonstrate, in a large-scale study spanning five decades of wheat breeding progress in western Europe, where grain yields are among the highest worldwide, that breeding for high performance in fact enhances cultivar performance not only under optimal production conditions but also in production systems with reduced agrochemical inputs. New cultivars incrementally accumulated genetic variants conferring favourable effects on key yield parameters, disease resistance, nutrient use efficiency, photosynthetic efficiency and grain quality. Combining beneficial, genome-wide haplotypes could help breeders to more efficiently exploit available genetic variation, optimizing future yield potential in more sustainable production systems.

Breeding improves wheat productivity under contrasting agrochemical input levels

The world cropping area for wheat exceeds that of any other crop, and high grain yields in intensive wheat cropping systems are essential for global food security. Breeding has raised yields dramatically in high-input production systems; however, selection under optimal growth conditions is widely believed to diminish the adaptive capacity of cultivars to less optimal cropping environments. Here, we demonstrate, in a large-scale study spanning five decades of wheat breeding progress in western Europe, where grain yields are among the highest worldwide, that breeding for high performance in fact enhances cultivar performance not only under optimal production conditions but also in production systems with reduced agrochemical inputs. New cultivars incrementally accumulated genetic variants conferring favourable effects on key yield parameters, disease resistance, nutrient use efficiency, photosynthetic efficiency and grain quality. Combining beneficial, genome-wide haplotypes could help breeders to more efficiently exploit available genetic variation, optimizing future yield potential in more sustainable production systems.

A Major Root Architecture QTL Responding to Water Limitation in Durum Wheat

The optimal root system architecture (RSA) of a crop is context dependent and critical for efficient resource capture in the soil. Narrow root growth angle promoting deeper root growth is often associated with improved access to water and nutrients in deep soils during terminal drought. RSA, therefore is a drought-adaptive trait that could minimize yield losses in regions with limited rainfall. Here, GWAS for seminal root angle (SRA) identified seven marker-trait associations clustered on chromosome 6A, representing a major quantitative trait locus (qSRA-6A) which also displayed high levels of pairwise LD (r 2 = 0.67). Subsequent haplotype analysis revealed significant differences between major groups. Candidate gene analysis revealed loci related to gravitropism, polar growth and hormonal signaling. No differences were observed for root biomass between lines carrying hap1 and hap2 for qSRA-6A, highlighting the opportunity to perform marker-assisted selection for the qSRA-6A locus and directly select for wide or narrow RSA, without influencing root biomass. Our study revealed that the genetic predisposition for deep rooting was best expressed under water-limitation, yet the root system displayed plasticity producing root growth in response to water availability in upper soil layers. We discuss the potential to deploy root architectural traits in cultivars to enhance yield stability in environments that experience limited rainfall.

Q&A: modern crop breeding for future food security

Farmers around the world have recently experienced significant crop losses due to severe heat and drought. Such extreme weather events and the need to feed a rapidly growing population have raised concerns for global food security. While plant breeding has been very successful and has delivered today’s highly productive crop varieties, the rate of genetic improvement must double to meet the projected future demands. Here we discuss basic principles and features of crop breeding and how modern technologies could efficiently be explored to boost crop improvement in the face of increasingly challenging production conditions.

Discovering new alleles for yellow spot resistance in the Vavilov wheat collection

GWAS detected 11 yellow spot resistance QTL in the Vavilov wheat collection. Promising adult-plant resistance loci could provide a sustainable genetic solution to yellow spot in modern wheat varieties. Yellow spot, caused by the fungal pathogen Pyrenophora tritici-repentis (Ptr), is the most economically damaging foliar disease of wheat in Australia. Genetic resistance is considered to be the most sustainable means for disease management, yet the genomic regions underpinning resistance to Ptr, particularly adult-plant resistance (APR), remain vastly unknown. In this study, we report results of a genome-wide association study using 295 accessions from the Vavilov wheat collection which were extensively tested for response to Ptr infections in glasshouse and field trials at both seedling an adult growth stages. Combining phenotypic datasets from multiple experiments in Australia and Russia with 25,286 genome-wide, high-quality DArTseq markers, we detected a total of 11 QTL, of which 5 were associated with seedling resistance, 3 with all-stage resistance, and 3 with APR. Interestingly, the novel APR QTL were effective even in the presence of host sensitivity gene Tsn1. These genomic regions could offer broad-spectrum yellow spot protection, not just to ToxA but also other pathogenicity or virulence factors. Vavilov wheat accessions carrying APR QTL combinations displayed enhanced levels of resistance highlighting the potential for QTL stacking through breeding. We propose that the APR genetic factors discovered in our study could be used to improve resistance levels in modern wheat varieties and contribute to the sustainable control of yellow spot.

Designer Roots for Future Crops

Despite the importance of roots, they have largely been ignored by modern crop research and breeding. We discuss important progress in crop root research and highlight how the context-dependent optimisation of below- and above-ground plant components provides opportunities to improve future crops in the face of increasing environmental fluctuations.

Genetic insights into underground responses to Fusarium graminearum infection in wheat

The ongoing global intensification of wheat production will likely be accompanied by a rising pressure of Fusarium diseases. While utmost attention was given to Fusarium head blight (FHB) belowground plant infections of the pathogen have largely been ignored. The current knowledge about the impact of soil borne Fusarium infection on plant performance and the underlying genetic mechanisms for resistance remain very limited. Here, we present the first large-scale investigation of Fusarium root rot (FRR) resistance using a diverse panel of 215 international wheat lines. We obtained data for a total of 21 resistance-related traits, including large-scale Real-time PCR experiments to quantify fungal spread. Association mapping and subsequent haplotype analyses discovered a number of highly conserved genomic regions associated with resistance, and revealed a significant effect of allele stacking on the stembase discoloration. Resistance alleles were accumulated in European winter wheat germplasm, implying indirect prior selection for improved FRR resistance in elite breeding programs. Our results give first insights into the genetic basis of FRR resistance in wheat and demonstrate how molecular parameters can successfully be explored in genomic prediction. Ongoing work will help to further improve our understanding of the complex interactions of genetic factors influencing FRR resistance.

Effective Genomic Selection in a Narrow-Genepool Crop with Low-Density Markers: Asian Rapeseed as an Example

Genomic selection (GS) has revolutionized breeding for quantitative traits in plants, offering potential to optimize resource allocation in breeding programs and increase genetic gain per unit of time. Modern high-density single nucleotide polymorphism (SNP) arrays comprising up to several hundred thousand markers provide a user-friendly technology to characterize the genetic constitution of whole populations and for implementing GS in breeding programs. However, GS does not build upon detailed genotype profiling facilitated by maximum marker density. With extensive genome-wide linkage disequilibrium (LD) being a common characteristic of breeding pools, fewer representative markers from available high-density genotyping platforms could be sufficient to capture the association between a genomic region and a phenotypic trait. To examine the effects of reduced marker density on genomic prediction accuracy, we collected data on three traits across 2 yr in a panel of 203 homozygous Chinese semiwinter rapeseed ( L.) inbred lines, broadly encompassing allelic variability in the Asian genepool. We investigated two approaches to selecting subsets of markers: a trait-dependent strategy based on genome-wide association study (GWAS) significance thresholds and a trait-independent method to detect representative tag SNPs. Prediction accuracies were evaluated using cross-validation with ridge-regression best linear unbiased predictions (rrBLUP). With semiwinter rapeseed as a model species, we demonstrate that low-density marker sets comprising a few hundred to a few thousand markers enable high prediction accuracies in breeding populations with strong LD comparable to those achieved with high-density arrays. Our results are valuable for facilitating routine application of cost-efficient GS in breeding programs.

Genome-wide regression models considering general and specific combining ability predict hybrid performance in oilseed rape with similar accuracy regardless of trait architecture

Genomic prediction using the Brassica 60 k genotyping array is efficient in oilseed rape hybrids. Prediction accuracy is more dependent on trait complexity than on the prediction model. In oilseed rape breeding programs, performance prediction of parental combinations is of fundamental importance. Due to the phenomenon of heterosis, per se performance is not a reliable indicator for F₁-hybrid performance, and selection of well-paired parents requires the testing of large quantities of hybrid combinations in extensive field trials. However, the number of potential hybrids, in general, dramatically exceeds breeding capacity and budget. Integration of genomic selection (GS) could substantially increase the number of potential combinations that can be evaluated. GS models can be used to predict the performance of untested individuals based only on their genotypic profiles, using marker effects previously predicted in a training population. This allows for a preselection of promising genotypes, enabling a more efficient allocation of resources. In this study, we evaluated the usefulness of the Illumina Brassica 60 k SNP array for genomic prediction and compared three alternative approaches based on a homoscedastic ridge regression BLUP and three Bayesian prediction models that considered general and specific combining ability (GCA and SCA, respectively). A total of 448 hybrids were produced in a commercial breeding program from unbalanced crosses between 220 paternal doubled haploid lines and five male-sterile testers. Predictive ability was evaluated for seven agronomic traits. We demonstrate that the Brassica 60 k genotyping array is an adequate and highly valuable platform to implement genomic prediction of hybrid performance in oilseed rape. Furthermore, we present first insights into the application of established statistical models for prediction of important agronomical traits with contrasting patterns of polygenic control.

Speed breeding for multiple quantitative traits in durum wheat

BACKGROUND

Plant breeding requires numerous generations to be cycled and evaluated before an improved cultivar is released. This lengthy process is required to introduce and test multiple traits of interest. However, a technology for rapid generation advance named 'speed breeding' was successfully deployed in bread wheat (Triticum aestivum L.) to achieve six generations per year while imposing phenotypic selection for foliar disease resistance and grain dormancy. Here, for the first time the deployment of this methodology is presented in durum wheat (Triticum durum Desf.) by integrating selection for key traits, including above and below ground traits on the same set of plants. This involved phenotyping for seminal root angle (RA), seminal root number (RN), tolerance to crown rot (CR), resistance to leaf rust (LR) and plant height (PH). In durum wheat, these traits are desirable in environments where yield is limited by in-season rainfall with the occurrence of CR and epidemics of LR. To evaluate this multi-trait screening approach, we applied selection to a large segregating F₂ population (n = 1000) derived from a bi-parental cross (Outrob4/Caparoi). A weighted selection index (SI) was developed and applied. The gain for each trait was determined by evaluating F₃ progeny derived from 100 'selected' and 100 'unselected' F₂ individuals.

RESULTS

Transgressive segregation was observed for all assayed traits in the Outrob4/Caparoi F₂ population. Application of the SI successfully shifted the population mean for four traits, as determined by a significant mean difference between 'selected' and 'unselected' F₃ families for CR tolerance, LR resistance, RA and RN. No significant shift for PH was observed.

CONCLUSIONS

The novel multi-trait phenotyping method presents a useful tool for rapid selection of early filial generations or for the characterization of fixed lines out-of-season. Further, it offers efficient use of resources by assaying multiple traits on the same set of plants. Results suggest that when performed in parallel with speed breeding in early generations, selection will enrich recombinant inbred lines with desirable alleles and will reduce the length and number of years required to combine these traits in elite breeding populations and therefore cultivars.

Unlocking new alleles for leaf rust resistance in the Vavilov wheat collection

Thirteen potentially new leaf rust resistance loci were identified in a Vavilov wheat diversity panel. We demonstrated the potential of allele stacking to strengthen resistance against this important pathogen. Leaf rust (LR) caused by Puccinia triticina is an important disease of wheat (Triticum aestivum L.), and the deployment of genetically resistant cultivars is the most viable strategy to minimise yield losses. In this study, we evaluated a diversity panel of 295 bread wheat accessions from the N. I. Vavilov Institute of Plant Genetic Resources (St Petersburg, Russia) for LR resistance and performed genome-wide association studies (GWAS) using 10,748 polymorphic DArT-seq markers. The diversity panel was evaluated at seedling and adult plant growth stages using three P. triticina pathotypes prevalent in Australia. GWAS was applied to 11 phenotypic data sets which identified a total of 52 significant marker-trait associations representing 31 quantitative trait loci (QTL). Among them, 29 QTL were associated with adult plant resistance (APR). Of the 31 QTL, 13 were considered potentially new loci, whereas 4 co-located with previously catalogued Lr genes and 14 aligned to regions reported in other GWAS and genomic prediction studies. One seedling LR resistance QTL located on chromosome 3A showed pronounced levels of linkage disequilibrium among markers (r ² = 0.7), suggested a high allelic fixation. Subsequent haplotype analysis for this region found seven haplotype variants, of which two were strongly associated with LR resistance at seedling stage. Similarly, analysis of an APR QTL on chromosome 7B revealed 22 variants, of which 4 were associated with resistance at the adult plant stage. Furthermore, most of the tested lines in the diversity panel carried 10 or more combined resistance-associated marker alleles, highlighting the potential of allele stacking for long-lasting resistance.

Genome-wide regression models considering general and specific combining ability predict hybrid performance in oilseed rape with similar accuracy regardless of trait architecture

KEY MESSAGE

Genomic prediction using the Brassica 60 k genotyping array is efficient in oilseed rape hybrids. Prediction accuracy is more dependent on trait complexity than on the prediction model. In oilseed rape breeding programs, performance prediction of parental combinations is of fundamental importance. Due to the phenomenon of heterosis, per se performance is not a reliable indicator for F₁-hybrid performance, and selection of well-paired parents requires the testing of large quantities of hybrid combinations in extensive field trials. However, the number of potential hybrids, in general, dramatically exceeds breeding capacity and budget. Integration of genomic selection (GS) could substantially increase the number of potential combinations that can be evaluated. GS models can be used to predict the performance of untested individuals based only on their genotypic profiles, using marker effects previously predicted in a training population. This allows for a preselection of promising genotypes, enabling a more efficient allocation of resources. In this study, we evaluated the usefulness of the Illumina Brassica 60 k SNP array for genomic prediction and compared three alternative approaches based on a homoscedastic ridge regression BLUP and three Bayesian prediction models that considered general and specific combining ability (GCA and SCA, respectively). A total of 448 hybrids were produced in a commercial breeding program from unbalanced crosses between 220 paternal doubled haploid lines and five male-sterile testers. Predictive ability was evaluated for seven agronomic traits. We demonstrate that the Brassica 60 k genotyping array is an adequate and highly valuable platform to implement genomic prediction of hybrid performance in oilseed rape. Furthermore, we present first insights into the application of established statistical models for prediction of important agronomical traits with contrasting patterns of polygenic control.

Exploring and Harnessing Haplotype Diversity to Improve Yield Stability in Crops

In order to meet future food, feed, fiber, and bioenergy demands, global yields of all major crops need to be increased significantly. At the same time, the increasing frequency of extreme weather events such as heat and drought necessitates improvements in the environmental resilience of modern crop cultivars. Achieving sustainably increase yields implies rapid improvement of quantitative traits with a very complex genetic architecture and strong environmental interaction. Latest advances in genome analysis technologies today provide molecular information at an ultrahigh resolution, revolutionizing crop genomic research, and paving the way for advanced quantitative genetic approaches. These include highly detailed assessment of population structure and genotypic diversity, facilitating the identification of selective sweeps and signatures of directional selection, dissection of genetic variants that underlie important agronomic traits, and genomic selection (GS) strategies that not only consider major-effect genes. Single-nucleotide polymorphism (SNP) markers today represent the genotyping system of choice for crop genetic studies because they occur abundantly in plant genomes and are easy to detect. SNPs are typically biallelic, however, hence their information content compared to multiallelic markers is low, limiting the resolution at which SNP-trait relationships can be delineated. An efficient way to overcome this limitation is to construct haplotypes based on linkage disequilibrium, one of the most important features influencing genetic analyses of crop genomes. Here, we give an overview of the latest advances in genomics-based haplotype analyses in crops, highlighting their importance in the context of polyploidy and genome evolution, linkage drag, and co-selection. We provide examples of how haplotype analyses can complement well-established quantitative genetics frameworks, such as quantitative trait analysis and GS, ultimately providing an effective tool to equip modern crops with environment-tailored characteristics.

Linkage drag constrains the roots of modern wheat

Roots, the hidden half of crop plants, are essential for resource acquisition. However, knowledge about the genetic control of below-ground plant development in wheat, one of the most important small-grain crops in the world, is very limited. The molecular interactions connecting root and shoot development and growth, and thus modulating the plant's demand for water and nutrients along with its ability to access them, are largely unexplored. Here, we demonstrate that linkage drag in European bread wheat, driven by strong selection for a haplotype variant controlling heading date, has eliminated a specific combination of two flanking, highly conserved haplotype variants whose interaction confers increased root biomass. Reversing this inadvertent consequence of selection could recover root diversity that may prove essential for future food production in fluctuating environments. Highly conserved synteny to rice across this chromosome segment suggests that adaptive selection has shaped the diversity landscape of this locus across different, globally important cereal crops. By mining wheat gene expression data, we identified root-expressed genes within the region of interest that could help breeders to select positive variants adapted to specific target soil environments.