Comparing the Potential of Marker-Assisted Selection and Genomic Prediction for Improving Rust Resistance in Hybrid Wheat

Toward combining qualitative race-specific and quantitative race-nonspecific disease resistance by genomic selection

A novel genomic selection strategy offers the unique opportunity to develop qualitative race-specific resistant varieties that possess high levels of the more durable quantitative race-nonspecific resistance in their genetic background. Race-specific qualitative resistance genes (R-genes) are conferring complete resistance in many pathosystems, but are frequently overcome by new virulent pathogen races. Once the deployed R-genes are overcome, a wide variation of quantitative disease resistance (QDR) can be observed in a set of previously race-specific, i.e., completely resistant genotypes-a phenomenon known as "vertifolia effect." This race-nonspecific QDR is considered to be more durable in the long term, but provides merely a partial protection against pathogens. This simulation study aimed to detangle race-specific R-gene-mediated resistance of pending selection candidates and the QDR in their genetic background by employing different genomic selection strategies. True breeding values that reflected performance data for rust resistance in wheat were simulated, and used in a recurrent genomic selection based on several prediction models and training population designs. Using training populations that were devoid of race-specific R-genes was thereby pivotal for an efficient improvement of QDR in the long term. Marker-assisted preselection for the presence of R-genes followed by a genomic prediction for accumulating the many small to medium effect loci underlying QDR in the genetic background of race-specific resistant genotypes appeared furthermore to be a promising approach to select simultaneously for both types of resistance. The practical application of such a knowledge-driven genomic breeding strategy offers the opportunity to develop varieties with multiple layers of resistance, which have the potential to prevent intolerable crop losses under epidemic situations by displaying a high level of QDR even when race-specific R-genes have been overcome by evolving pathogen populations.

Machine learning bridges omics sciences and plant breeding

Some of the biological knowledge obtained from fundamental research will be implemented in applied plant breeding. To bridge basic research and breeding practice, machine learning (ML) holds great promise to translate biological knowledge and omics data into precision-designed plant breeding. Here, we review ML for multi-omics analysis in plants, including data dimensionality reduction, inference of gene-regulation networks, and gene discovery and prioritization. These applications will facilitate understanding trait regulation mechanisms and identifying target genes potentially applicable to knowledge-driven molecular design breeding. We also highlight applications of deep learning in plant phenomics and ML in genomic selection-assisted breeding, such as various ML algorithms that model the correlations among genotypes (genes), phenotypes (traits), and environments, to ultimately achieve data-driven genomic design breeding.

Evaluation of the resistance to Chinese predominant races of Puccinia triticina and analysis of effective leaf rust resistance genes in wheat accessions from the U.S. National Plant Germplasm System

Puccinia triticina, which is the causative agent of wheat leaf rust, is widely spread in China and most other wheat-planting countries around the globe. Cultivating resistant wheat cultivars is the most economical, effective, and environmentally friendly method for controlling leaf rust-caused yield damage. Exploring the source of resistance is very important in wheat resistance breeding programs. In order to explore more effective resistance sources for wheat leaf rust, the resistance of 112 wheat accessions introduced from the U.S. National Plant Germplasm System were identified using a mixture of pathogenic isolates of THTT, THTS, PHTT, THJT and THJS which are the most predominant races in China. As a result, all of these accessions showed high resistance at seedling stage, of which, ninety-nine accessions exhibited resistance at adult plant stage. Eleven molecular markers of eight effective leaf rust resistance genes in China were used to screen the 112 accessions. Seven effective leaf rust resistance genes Lr9, Lr19, Lr24, Lr28, Lr29, Lr38 and Lr45 were detected, except Lr47. Twenty-three accessions had only one of those seven effective leaf rust resistance gene. Eleven accessions carried Lr24+Lr38, and 7 accessions carried Lr9+Lr24+Lr38, Lr24+Lr38+Lr45, Lr24+Lr29+Lr38 and Lr19+Lr38+Lr45 respectively. The remaining seventy-one accessions had none of those eight effective leaf rust resistance genes. This study will provide theoretical guidance for rational utilization of these introduted wheat accessions directly or for breeding the resistant wheat cultivars.

Genome-wide association study and genomic prediction of resistance to stripe rust in current Central and Northern European winter wheat germplasm

We found two loci on chromosomes 2BS and 6AL that significantly contribute to stripe rust resistance in current European winter wheat germplasm. Stripe or yellow rust, caused by the fungus Puccinia striiformis Westend f. sp. tritici, is one of the most destructive wheat diseases. Sustainable management of wheat stripe rust can be achieved through the deployment of rust resistant cultivars. To detect effective resistance loci for use in breeding programs, an association mapping panel of 230 winter wheat cultivars and breeding lines from Northern and Central Europe was employed. Genotyping with the Illumina® iSelect® 25 K Infinium® single nucleotide polymorphism (SNP) genotyping array yielded 8812 polymorphic markers. Structure analysis revealed two subpopulations with 92 Austrian breeding lines and cultivars, which were separated from the other 138 genotypes from Germany, Norway, Sweden, Denmark, Poland, and Switzerland. Genome-wide association study for adult plant stripe rust resistance identified 12 SNP markers on six wheat chromosomes which showed consistent effects over several testing environments. Among these, two marker loci on chromosomes 2BS (RAC875_c1226_652) and 6AL (Tdurum_contig29607_413) were highly predictive in three independent validation populations of 1065, 1001, and 175 breeding lines. Lines with the resistant haplotype at both loci were nearly free of stipe rust symptoms. By using mixed linear models with those markers as fixed effects, we could increase predictive ability in the three populations by 0.13-0.46 compared to a standard genomic best linear unbiased prediction approach. The obtained results facilitate an efficient selection for stripe rust resistance against the current pathogen population in the Northern and Central European winter wheat gene pool.

Comparison of single-trait and multi-trait genomic predictions on agronomic and disease resistance traits in spring wheat

This study performed comprehensive analyses on the predictive abilities of single-trait and two multi-trait models in three populations. Our results demonstrated the superiority of multi-traits over single-trait models across seven agronomic and four to seven disease resistance traits of different genetic architecture. The predictive ability of multi-trait and single-trait prediction models has not been investigated on diverse traits evaluated under organic and conventional management systems. Here, we compared the predictive abilities of 25% of a testing set that has not been evaluated for a single trait (ST), not evaluated for multi-traits (MT1), and evaluated for some traits but not others (MT2) in three spring wheat populations genotyped either with the wheat 90K single nucleotide polymorphisms array or DArTseq. Analyses were performed on seven agronomic traits evaluated under conventional and organic management systems, four to seven disease resistance traits, and all agronomic and disease resistance traits simultaneously. The average prediction accuracies of the ST, MT1, and MT2 models varied from 0.03 to 0.78 (mean 0.41), from 0.05 to 0.82 (mean 0.47), and from 0.05 to 0.92 (mean 0.67), respectively. The predictive ability of the MT2 model was significantly greater than the ST model in all traits and populations except common bunt with the MT1 model being intermediate between them. The MT2 model increased prediction accuracies over the ST and MT1 models in all traits by 9.0-82.4% (mean 37.3%) and 2.9-82.5% (mean 25.7%), respectively, except common bunt that showed up to 7.7% smaller accuracies in two populations. A joint analysis of all agronomic and disease resistance traits further improved accuracies within the MT1 and MT2 models on average by 21.4% and 17.4%, respectively, as compared to either the agronomic or disease resistance traits, demonstrating the high potential of the multi-traits models in improving prediction accuracies.

Genomic Prediction Accuracy of Stripe Rust in Six Spring Wheat Populations by Modeling Genotype by Environment Interaction

Some previous studies have assessed the predictive ability of genome-wide selection on stripe (yellow) rust resistance in wheat, but the effect of genotype by environment interaction (GEI) in prediction accuracies has not been well studied in diverse genetic backgrounds. Here, we compared the predictive ability of a model based on phenotypic data only (M1), the main effect of phenotype and molecular markers (M2), and a model that incorporated GEI (M3) using three cross-validations (CV1, CV2, and CV0) scenarios of interest to breeders in six spring wheat populations. Each population was evaluated at three to eight field nurseries and genotyped with either the DArTseq technology or the wheat 90K single nucleotide polymorphism arrays, of which a subset of 1,058- 23,795 polymorphic markers were used for the analyses. In the CV1 scenario, the mean prediction accuracies of the M1, M2, and M3 models across the six populations varied from −0.11 to −0.07, from 0.22 to 0.49, and from 0.19 to 0.48, respectively. Mean accuracies obtained using the M3 model in the CV1 scenario were significantly greater than the M2 model in two populations, the same in three populations, and smaller in one population. In both the CV2 and CV0 scenarios, the mean prediction accuracies of the three models varied from 0.53 to 0.84 and were not significantly different in all populations, except the Attila/CDC Go in the CV2, where the M3 model gave greater accuracy than both the M1 and M2 models. Overall, the M3 model increased prediction accuracies in some populations by up to 12.4% and decreased accuracy in others by up to 17.4%, demonstrating inconsistent results among genetic backgrounds that require considering each population separately. This is the first comprehensive genome-wide prediction study that investigated details of the effect of GEI on stripe rust resistance across diverse spring wheat populations.

Discovery of Resistance Genes in Rye by Targeted Long-Read Sequencing and Association Genetics

The majority of released rye cultivars are susceptible to leaf rust because of a low level of resistance in the predominant hybrid rye-breeding gene pools Petkus and Carsten. To discover new sources of leaf rust resistance, we phenotyped a diverse panel of inbred lines from the less prevalent Gülzow germplasm using six distinct isolates of Puccinia recondita f. sp. secalis and found that 55 out of 92 lines were resistant to all isolates. By performing a genome-wide association study using 261,406 informative SNP markers, we identified five resistance-associated QTLs on chromosome arms 1RS, 1RL, 2RL, 5RL and 7RS. To identify candidate Puccinia recondita (Pr) resistance genes in these QTLs, we sequenced the rye nucleotide-binding leucine-rich repeat (NLR) intracellular immune receptor complement using a Triticeae NLR bait-library and PacBio^® long-read single-molecule high-fidelity (HiFi) sequencing. Trait-genotype correlations across 10 resistant and 10 susceptible lines identified four candidate NLR-encoding Pr genes. One of these physically co-localized with molecular markers delimiting Pr3 on chromosome arm 1RS and the top-most resistance-associated QTL in the panel.

Hybrid wheat: quantitative genetic parameters and heterosis for quality and rheological traits as well as baking volume

Heterosis effects for dough quality and baking volume were close to zero. However, hybrids have a higher grain yield at a given level of bread making quality compared to their parental lines. Bread wheat cultivars have been selected according to numerous quality traits to fulfill the requirements of the bread making industry. These include beside protein content and quality also rheological traits and baking volume. We evaluated 35 male and 73 female lines and 119 of their single-cross hybrids at three different locations for grain yield, protein content, sedimentation value, extensograph traits and baking volume. No significant differences (p < 0.05) were found in the mean comparisons of males, females and hybrids, except for higher grain yield and lower protein content in the hybrids. Mid-parent and better-parent heterosis values were close to zero and slightly negative, respectively, for baking volume and extensograph traits. However, the majority of heterosis values resulted in the finding that hybrids had higher grain yield than lines for a given level of baking volume, sedimentation value or energy value of extensograph. Due to the high correlation with the mid-parent values (r > 0.70), an initial prediction of hybrid performance based on line per se performance for protein content, sedimentation value, most traits of the extensograph and baking volume is possible. The low variance due to specific combining ability effects for most quality traits points toward an additive gene action requires quality selection within both heterotic groups. Consequently, hybrid wheat can combine high grain yield with high bread making quality. However, the future use of wheat hybrids strongly depends on the establishment of a cost-efficient and reliable seed production system.

Genomic Predictions for Common Bunt, FHB, Stripe Rust, Leaf Rust, and Leaf Spotting Resistance in Spring Wheat

Some studies have investigated the potential of genomic selection (GS) on stripe rust, leaf rust, Fusarium head blight (FHB), and leaf spot in wheat, but none of them have assessed the effect of the reaction norm model that incorporated GE interactions. In addition, the prediction accuracy on common bunt has not previously been studied. Here, we investigated within-population prediction accuracies using the baseline M1 model and two reaction norm models (M2 and M3) with three random cross-validation (CV1, CV2, and CV0) schemes. Three Canadian spring wheat populations were evaluated in up to eight field environments and genotyped with 3158, 5732, and 23,795 polymorphic markers. The M3 model that incorporated GE interactions reduced residual variance by an average of 10.2% as compared with the main effect M2 model and increased prediction accuracies on average by 2–6%. In some traits, the M3 model increased prediction accuracies up to 54% as compared with the M2 model. The average prediction accuracies of the M3 model with CV1, CV2, and CV0 schemes varied from 0.02 to 0.48, from 0.25 to 0.84, and from 0.14 to 0.87, respectively. In both CV2 and CV0 schemes, stripe rust in all three populations, common bunt and leaf rust in two populations, as well as FHB severity, FHB index, and leaf spot in one population had high to very high (0.54–0.87) prediction accuracies. This is the first comprehensive genomic selection study on five major diseases in spring wheat.

Genome-based prediction of agronomic traits in spring wheat under conventional and organic management systems

Using phenotype data of three spring wheat populations evaluated at 6-15 environments under two management systems, we found moderate to very high prediction accuracies across seven traits. The phenotype data collected under an organic management system effectively predicted the performance of lines in the conventional management and vice versa. There is growing interest in developing wheat cultivars specifically for organic agriculture, but we are not aware of the effect of organic management on the predictive ability of genomic selection (GS). Here, we evaluated within populations prediction accuracies of four GS models, four combinations of training and testing sets, three reaction norm models, and three random cross-validations (CV) schemes in three populations phenotyped under organic and conventional management systems. Our study was based on a total of 578 recombinant inbred lines and varieties from three spring wheat populations, which were evaluated for seven traits at 3-9 conventionally and 3-6 organically managed field environments and genotyped either with the wheat 90 K SNP array or DArTseq. We predicted the management systems (CV0_M) or environments (CV0), a subset of lines that have been evaluated in either management (CV2_M) or some environments (CV2), and the performance of newly developed lines in either management (CV1_M) or environments (CV1). The average prediction accuracies of the model that incorporated genotype × environment interactions with CV0 and CV2 schemes varied from 0.69 to 0.97. In the CV1 and CV1M schemes, prediction accuracies ranged from - 0.12 to 0.77 depending on the reaction norm models, the traits, and populations. In most cases, grain protein showed the highest prediction accuracies. The phenotype data collected under the organic management effectively predicted the performance of lines under conventional management and vice versa. This is the first comprehensive GS study that investigated the effect of the organic management system in wheat.

Prospectus of Genomic Selection and Phenomics in Cereal, Legume and Oilseed Breeding Programs

The last decade witnessed an unprecedented increase in the adoption of genomic selection (GS) and phenomics tools in plant breeding programs, especially in major cereal crops. GS has demonstrated the potential for selecting superior genotypes with high precision and accelerating the breeding cycle. Phenomics is a rapidly advancing domain to alleviate phenotyping bottlenecks and explores new large-scale phenotyping and data acquisition methods. In this review, we discuss the lesson learned from GS and phenomics in six self-pollinated crops, primarily focusing on rice, wheat, soybean, common bean, chickpea, and groundnut, and their implementation schemes are discussed after assessing their impact in the breeding programs. Here, the status of the adoption of genomics and phenomics is provided for those crops, with a complete GS overview. GS’s progress until 2020 is discussed in detail, and relevant information and links to the source codes are provided for implementing this technology into plant breeding programs, with most of the examples from wheat breeding programs. Detailed information about various phenotyping tools is provided to strengthen the field of phenomics for a plant breeder in the coming years. Finally, we highlight the benefits of merging genomic selection, phenomics, and machine and deep learning that have resulted in extraordinary results during recent years in wheat, rice, and soybean. Hence, there is a potential for adopting these technologies into crops like the common bean, chickpea, and groundnut. The adoption of phenomics and GS into different breeding programs will accelerate genetic gain that would create an impact on food security, realizing the need to feed an ever-growing population.

Discovery of a Novel Leaf Rust (Puccinia recondita) Resistance Gene in Rye (Secale cereale L.) Using Association Genomics

Leaf rust constitutes one of the most important foliar diseases in rye (Secale cereale L.). To discover new sources of resistance, we phenotyped 180 lines belonging to a less well-characterized Gülzow germplasm at three field trial locations in Denmark and Northern Germany in 2018 and 2019. We observed lines with high leaf rust resistance efficacy at all locations in both years. A genome-wide association study using 261,406 informative single-nucleotide polymorphisms revealed two genomic regions associated with resistance on chromosome arms 1RS and 7RS, respectively. The most resistance-associated marker on chromosome arm 1RS physically co-localized with molecular markers delimiting Pr3. In the reference genomes Lo7 and Weining, the genomic region associated with resistance on chromosome arm 7RS contained a large number of nucleotide-binding leucine-rich repeat (NLR) genes. Residing in close proximity to the most resistance-associated marker, we identified a cluster of NLRs exhibiting close protein sequence similarity with the wheat leaf rust Lr1 gene situated on chromosome arm 5DL in wheat, which is syntenic to chromosome arm 7RS in rye. Due to the close proximity to the most resistance-associated marker, our findings suggest that the considered leaf rust R gene, provisionally denoted Pr6, could be a Lr1 ortholog in rye.