Enhancing the potential of phenomic and genomic prediction in winter wheat breeding using high-throughput phenotyping and deep learning.

Integrating high-throughput phenotyping (HTP) based traits into phenomic and genomic selection (GS) can accelerate the breeding of high-yielding and climate-resilient wheat cultivars. In this study, we explored the applicability of Unmanned Aerial Vehicles (UAV)-assisted HTP combined with deep learning (DL) for the phenomic or multi-trait (MT) genomic prediction of grain yield (GY), test weight (TW), and grain protein content (GPC) in winter wheat. Significant correlations were observed between agronomic traits and HTP-based traits across different growth stages of winter wheat. Using a deep neural network (DNN) model, HTP-based phenomic predictions showed robust prediction accuracies for GY, TW, and GPC for a single location with R2 of 0.71, 0.62, and 0.49, respectively. Further prediction accuracies increased (R2 of 0.76, 0.64, and 0.75) for GY, TW, and GPC, respectively when advanced breeding lines from multi-locations were used in the DNN model. Prediction accuracies for GY varied across growth stages, with the highest accuracy at the Feekes 11 (Milky ripe) stage. Furthermore, forward prediction of GY in preliminary breeding lines using DNN trained on multi-location data from advanced breeding lines improved the prediction accuracy by 32% compared to single-location data. Next, we evaluated the potential of incorporating HTP-based traits in multi-trait genomic selection (MT-GS) models in the prediction of GY, TW, and GPC. MT-GS, models including UAV data-based anthocyanin reflectance index (ARI), green chlorophyll index (GCI), and ratio vegetation index 2 (RVI_2) as covariates demonstrated higher predictive ability (0.40, 0.40, and 0.37, respectively) as compared to single-trait model (0.23) for GY. Overall, this study demonstrates the potential of integrating HTP traits into DL-based phenomic or MT-GS models for enhancing breeding efficiency.

Integrating genomics, phenomics, and deep learning improves the predictive ability for Fusarium head blight-related traits in winter wheat.

Fusarium head blight (FHB) remains one of the most destructive diseases of wheat (Triticum aestivum L.), causing considerable losses in yield and end-use quality. Phenotyping of FHB resistance traits, Fusarium-damaged kernels (FDK), and deoxynivalenol (DON), is either prone to human biases or resource expensive, hindering the progress in breeding for FHB-resistant cultivars. Though genomic selection (GS) can be an effective way to select these traits, inaccurate phenotyping remains a hurdle in exploiting this approach. Here, we used an artificial intelligence (AI)-based precise FDK estimation that exhibits high heritability and correlation with DON. Further, GS using AI-based FDK (FDK_QVIS/FDK_QNIR) showed a two-fold increase in predictive ability (PA) compared to GS for traditionally estimated FDK (FDK_V). Next, the AI-based FDK was evaluated along with other traits in multi-trait (MT) GS models to predict DON. The inclusion of FDK_QNIR and FDK_QVIS with days to heading as covariates improved the PA for DON by 58% over the baseline single-trait GS model. We next used hyperspectral imaging of FHB-infected wheat kernels as a novel avenue to improve the MT GS for DON. The PA for DON using selected wavebands derived from hyperspectral imaging in MT GS models surpassed the single-trait GS model by around 40%. Finally, we evaluated phenomic prediction for DON by integrating hyperspectral imaging with deep learning to directly predict DON in FHB-infected wheat kernels and observed an accuracy (R2 = 0.45) comparable to best-performing MT GS models. This study demonstrates the potential application of AI and vision-based platforms to improve PA for FHB-related traits using genomic and phenomic selection.

Multi-trait genomic selection improves the prediction accuracy of end-use quality traits in hard winter wheat.

Improvement of end-use quality remains one of the most important goals in hard winter wheat (HWW) breeding. Nevertheless, the evaluation of end-use quality traits is confined to later development generations owing to resource-intensive phenotyping. Genomic selection (GS) has shown promise in facilitating selection for end-use quality; however, lower prediction accuracy (PA) for complex traits remains a challenge in GS implementation. Multi-trait genomic prediction (MTGP) models can improve PA for complex traits by incorporating information on correlated secondary traits, but these models remain to be optimized in HWW. A set of advanced breeding lines from 2015 to 2021 were genotyped with 8725 single-nucleotide polymorphisms and was used to evaluate MTGP to predict various end-use quality traits that are otherwise difficult to phenotype in earlier generations. The MTGP model outperformed the ST model with up to a twofold increase in PA. For instance, PA was improved from 0.38 to 0.75 for bake absorption and from 0.32 to 0.52 for loaf volume. Further, we compared MTGP models by including different combinations of easy-to-score traits as covariates to predict end-use quality traits. Incorporation of simple traits, such as flour protein (FLRPRO) and sedimentation weight value (FLRSDS), substantially improved the PA of MT models. Thus, the rapid low-cost measurement of traits like FLRPRO and FLRSDS can facilitate the use of GP to predict mixograph and baking traits in earlier generations and provide breeders an opportunity for selection on end-use quality traits by culling inferior lines to increase selection accuracy and genetic gains.

Genome-wide association analysis of spike and kernel traits in the U.S. hard winter wheat.

A better understanding of the genetic control of spike and kernel traits that have higher heritability can help in the development of high-yielding wheat varieties. Here, we identified the marker-trait associations (MTAs) for various spike- and kernel-related traits in winter wheat (Triticum aestivum L.) through genome-wide association studies (GWAS). An association mapping panel comprising 297 hard winter wheat accessions from the U.S. Great Plains was evaluated for eight spike- and kernel-related traits in three different environments. A GWAS using 15,590 single-nucleotide polymorphisms (SNPs) identified a total of 53 MTAs for seven spike- and kernel-related traits, where the highest number of MTAs were identified for spike length (16) followed by the number of spikelets per spike (15) and spikelet density (11). Out of 53 MTAs, 14 were considered to represent stable quantitative trait loci (QTL) as they were identified in multiple environments. Five multi-trait MTAs were identified for various traits including the number of spikelets per spike (NSPS), spikelet density (SD), kernel width (KW), and kernel area (KA) that could facilitate the pyramiding of yield-contributing traits. Further, a significant additive effect of accumulated favorable alleles on the phenotype of four spike-related traits suggested that breeding lines and cultivars with a higher number of favorable alleles could be a valuable resource for breeders to improve yield-related traits. This study improves the understanding of the genetic basis of yield-related traits in hard winter wheat and provides reliable molecular markers that will facilitate marker-assisted selection (MAS) in wheat breeding programs.

Multi-Locus Genome-Wide Association Studies to Characterize Fusarium Head Blight (FHB) Resistance in Hard Winter Wheat.

Fusarium head blight (FHB), caused by the fungus Fusarium graminearum Schwabe is an important disease of wheat that causes severe yield losses along with serious quality concerns. Incorporating the host resistance from either wild relatives, landraces, or exotic materials remains challenging and has shown limited success. Therefore, a better understanding of the genetic basis of native FHB resistance in hard winter wheat (HWW) and combining it with major quantitative trait loci (QTLs) can facilitate the development of FHB-resistant cultivars. In this study, we evaluated a set of 257 breeding lines from the South Dakota State University (SDSU) breeding program to uncover the genetic basis of native FHB resistance in the US hard winter wheat. We conducted a multi-locus genome-wide association study (ML-GWAS) with 9,321 high-quality single-nucleotide polymorphisms (SNPs). A total of six distinct marker-trait associations (MTAs) were identified for the FHB disease index (DIS) on five different chromosomes including 2A, 2B, 3B, 4B, and 7A. Further, eight MTAs were identified for Fusarium-damaged kernels (FDK) on six chromosomes including 3B, 5A, 6B, 6D, 7A, and 7B. Out of the 14 significant MTAs, 10 were found in the proximity of previously reported regions for FHB resistance in different wheat classes and were validated in HWW, while four MTAs represent likely novel loci for FHB resistance. Accumulation of favorable alleles of reported MTAs resulted in significantly lower mean DIS and FDK score, demonstrating the additive effect of FHB resistance alleles. Candidate gene analysis for two important MTAs identified several genes with putative proteins of interest; however, further investigation of these regions is needed to identify genes conferring FHB resistance. The current study sheds light on the genetic basis of native FHB resistance in the US HWW germplasm and the resistant lines and MTAs identified in this study will be useful resources for FHB resistance breeding via marker-assisted selection.

Whole-genome analysis of hard winter wheat germplasm identifies genomic regions associated with spike and kernel traits.

Genetic dissection of yield component traits including spike and kernel characteristics is essential for the continuous improvement in wheat yield. Genome-wide association studies (GWAS) have been frequently used to identify genetic determinants for spike and kernel-related traits in wheat, though none have been employed in hard winter wheat (HWW) which represents a major class in US wheat acreage. Further, most of these studies relied on assembled diversity panels instead of adapted breeding lines, limiting the transferability of results to practical wheat breeding. Here we assembled a population of advanced/elite breeding lines and well-adapted cultivars and evaluated over four environments for phenotypic analysis of spike and kernel traits. GWAS identified 17 significant multi-environment marker-trait associations (MTAs) for various traits, representing 12 putative quantitative trait loci (QTLs), with five QTLs affecting multiple traits. Four of these QTLs mapped on three chromosomes 1A, 5B, and 7A for spike length, number of spikelets per spike (NSPS), and kernel length are likely novel. Further, a highly significant QTL was detected on chromosome 7AS that has not been previously associated with NSPS and putative candidate genes were identified in this region. The allelic frequencies of important quantitative trait nucleotides (QTNs) were deduced in a larger set of 1,124 accessions which revealed the importance of identified MTAs in the US HWW breeding programs. The results from this study could be directly used by the breeders to select the lines with favorable alleles for making crosses, and reported markers will facilitate marker-assisted selection of stable QTLs for yield components in wheat breeding.

Genomic variants affecting homoeologous gene expression dosage contribute to agronomic trait variation in allopolyploid wheat.

Allopolyploidy greatly expands the range of possible regulatory interactions among functionally redundant homoeologous genes. However, connection between the emerging regulatory complexity and expression and phenotypic diversity in polyploid crops remains elusive. Here, we use diverse wheat accessions to map expression quantitative trait loci (eQTL) and evaluate their effects on the population-scale variation in homoeolog expression dosage. The relative contribution of cis- and trans-eQTL to homoeolog expression variation is strongly affected by both selection and demographic events. Though trans-acting effects play major role in expression regulation, the expression dosage of homoeologs is largely influenced by cis-acting variants, which appear to be subjected to selection. The frequency and expression of homoeologous gene alleles showing strong expression dosage bias are predictive of variation in yield-related traits, and have likely been impacted by breeding for increased productivity. Our study highlights the importance of genomic variants affecting homoeolog expression dosage in shaping agronomic phenotypes and points at their potential utility for improving yield in polyploid crops.

Multi-Trait Multi-Environment Genomic Prediction of Agronomic Traits in Advanced Breeding Lines of Winter Wheat.

Genomic prediction is a promising approach for accelerating the genetic gain of complex traits in wheat breeding. However, increasing the prediction accuracy (PA) of genomic prediction (GP) models remains a challenge in the successful implementation of this approach. Multivariate models have shown promise when evaluated using diverse panels of unrelated accessions; however, limited information is available on their performance in advanced breeding trials. Here, we used multivariate GP models to predict multiple agronomic traits using 314 advanced and elite breeding lines of winter wheat evaluated in 10 site-year environments. We evaluated a multi-trait (MT) model with two cross-validation schemes representing different breeding scenarios (CV1, prediction of completely unphenotyped lines; and CV2, prediction of partially phenotyped lines for correlated traits). Moreover, extensive data from multi-environment trials (METs) were used to cross-validate a Bayesian multi-trait multi-environment (MTME) model that integrates the analysis of multiple-traits, such as G × E interaction. The MT-CV2 model outperformed all the other models for predicting grain yield with significant improvement in PA over the single-trait (ST-CV1) model. The MTME model performed better for all traits, with average improvement over the ST-CV1 reaching up to 19, 71, 17, 48, and 51% for grain yield, grain protein content, test weight, plant height, and days to heading, respectively. Overall, the empirical analyses elucidate the potential of both the MT-CV2 and MTME models when advanced breeding lines are used as a training population to predict related preliminary breeding lines. Further, we evaluated the practical application of the MTME model in the breeding program to reduce phenotyping cost using a sparse testing design. This showed that complementing METs with GP can substantially enhance resource efficiency. Our results demonstrate that multivariate GS models have a great potential in implementing GS in breeding programs.

Genome-wide association analysis permits characterization of Stagonospora nodorum blotch (SNB) resistance in hard winter wheat.

Stagonospora nodorum blotch (SNB) is an economically important wheat disease caused by the necrotrophic fungus Parastagonospora nodorum. SNB resistance in wheat is controlled by several quantitative trait loci (QTLs). Thus, identifying novel resistance/susceptibility QTLs is crucial for continuous improvement of the SNB resistance. Here, the hard winter wheat association mapping panel (HWWAMP) comprising accessions from breeding programs in the Great Plains region of the US, was evaluated for SNB resistance and necrotrophic effectors (NEs) sensitivity at the seedling stage. A genome-wide association study (GWAS) was performed to identify single-nucleotide polymorphism (SNP) markers associated with SNB resistance and effectors sensitivity. We found seven significant associations for SNB resistance/susceptibility distributed over chromosomes 1B, 2AL, 2DS, 4AL, 5BL, 6BS, and 7AL. Two new QTLs for SNB resistance/susceptibility at the seedling stage were identified on chromosomes 6BS and 7AL, whereas five QTLs previously reported in diverse germplasms were validated. Allele stacking analysis at seven QTLs explained the additive and complex nature of SNB resistance. We identified accessions ('Pioneer-2180' and 'Shocker') with favorable alleles at five of the seven identified loci, exhibiting a high level of resistance against SNB. Further, GWAS for sensitivity to NEs uncovered significant associations for SnToxA and SnTox3, co-locating with previously identified host sensitivity genes (Tsn1 and Snn3). Candidate region analysis for SNB resistance revealed 35 genes of putative interest with plant defense response-related functions. The QTLs identified and validated in this study could be easily employed in breeding programs using the associated markers to enhance the SNB resistance in hard winter wheat.

The acyltransferase PMAT1 malonylates brassinolide glucoside.

Brassinosteroids (BRs) are steroid hormones of plants that coordinate fundamental growth and development processes. Their homeostasis is controlled by diverse means, including glucosylation of the bioactive BR brassinolide (BL), which is catalyzed by the UDP-glycosyltransferases (UGTs) UGT73C5 and UGT73C6 and occurs mainly at the C-23 position. Additional evidence had suggested that the resultant BL-23-O-glucoside (BL-23-O-Glc) can be malonylated, but the physiological significance of and enzyme required for this reaction had remained unknown. Here, we show that in Arabidopsis thaliana malonylation of BL-23-O-Glc is catalyzed by the acyltransferase phenolic glucoside malonyl-transferase 1 (PMAT1), which is also known to malonylate phenolic glucosides and lipid amides. Loss of PMAT1 abolished BL-23-O-malonylglucoside formation and enriched BL-23-O-Glc, showing that the enzyme acts on the glucoside. An overexpression of PMAT1 in plants where UGT73C6 was also overexpressed, and thus, BL-23-O-Glc formation was promoted, enhanced the symptoms of BR-deficiency of UGT73C6oe plants, providing evidence that PMAT1 contributes to BL inactivation. Based on these results, a model is proposed in which PMAT1 acts in the conversion of both endogenous and xenobiotic glucosides to adjust metabolic homeostasis in spatial and temporal modes.

Mining and genomic characterization of resistance to tan spot, Stagonospora nodorum blotch (SNB), and Fusarium head blight in Watkins core collection of wheat landraces.

BACKGROUND: In the late 1920s, A. E. Watkins collected about 7000 landrace cultivars (LCs) of bread wheat (Triticum aestivum L.) from 32 different countries around the world. Among which 826 LCs remain viable and could be a valuable source of superior/favorable alleles to enhance disease resistance in wheat. In the present study, a core set of 121 LCs, which captures the majority of the genetic diversity of Watkins collection, was evaluated for identifying novel sources of resistance against tan spot, Stagonospora nodorum blotch (SNB), and Fusarium Head Blight (FHB). RESULTS: A diverse response was observed in 121 LCs for all three diseases. The majority of LCs were moderately susceptible to susceptible to tan spot Ptr race 1 (84%) and FHB (96%) whereas a large number of LCs were resistant or moderately resistant against tan spot Ptr race 5 (95%) and SNB (54%). Thirteen LCs were identified in this study could be a valuable source for multiple resistance to tan spot Ptr races 1 and 5, and SNB, and another five LCs could be a potential source for FHB resistance. GWAS analysis was carried out using disease phenotyping score and 8807 SNPs data of 118 LCs, which identified 30 significant marker-trait associations (MTAs) with -log10 (p-value) > 3.0. Ten, five, and five genomic regions were found to be associated with resistance to tan spot Ptr race 1, race 5, and SNB, respectively in this study. In addition to Tsn1, several novel genomic regions Q.Ts1.sdsu-4BS and Q.Ts1.sdsu-5BS (tan spot Ptr race 1) and Q.Ts5.sdsu-1BL, Q.Ts5.sdsu-2DL, Q.Ts5.sdsu-3AL, and Q.Ts5.sdsu-6BL (tan spot Ptr race 5) were also identified. Our results indicate that these putative genomic regions contain several genes that play an important role in plant defense mechanisms. CONCLUSION: Our results suggest the existence of valuable resistant alleles against leaf spot diseases in Watkins LCs. The single-nucleotide polymorphism (SNP) markers linked to the quantitative trait loci (QTLs) for tan spot and SNB resistance along with LCs harboring multiple disease resistance could be useful for future wheat breeding.

Genome-Wide Association Study Uncovers Novel Genomic Regions Associated With Coleoptile Length in Hard Winter Wheat.

Successful seedling establishment depends on the optimum depth of seed placement especially in drought-prone conditions, providing an opportunity to exploit subsoil water and increase winter survival in winter wheat. Coleoptile length is a key determinant for the appropriate depth at which seed can be sown. Thus, understanding the genetic basis of coleoptile length is necessary and important for wheat breeding. We conducted a genome-wide association study (GWAS) using a diverse panel of 298 winter wheat genotypes to dissect the genetic architecture of coleoptile length. We identified nine genomic regions associated with the coleoptile length on seven different chromosomes. Of the nine genomic regions, five have been previously reported in various studies, including one mapped to previously known Rht-B1 region. Three novel quantitative trait loci (QTLs), QCL.sdsu-2AS, QCL.sdsu-4BL, and QCL.sdsu-5BL were identified in our study. QCL.sdsu-5BL has a large substitution effect which is comparable to Rht-B1's effect and could be used to compensate for the negative effect of Rht-B1 on coleoptile length. In total, the nine QTLs explained 59% of the total phenotypic variation. Cultivars 'Agate' and 'MT06103' have the longest coleoptile length and interestingly, have favorable alleles at nine and eight coleoptile loci, respectively. These lines could be a valuable germplasm for longer coleoptile breeding. Gene annotations in the candidate regions revealed several putative proteins of specific interest including cytochrome P450-like, expansins, and phytochrome A. The QTLs for coleoptile length linked to single-nucleotide polymorphism (SNP) markers reported in this study could be employed in marker-assisted breeding for longer coleoptile in wheat. Thus, our study provides valuable insights into the genetic and molecular regulation of the coleoptile length in winter wheat.