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

Association mapping and genomic prediction for processing and end-use quality traits in wheat (Triticum aestivum L.)

End-use and processing traits in wheat (Triticum aestivum L.) are crucial for varietal development but are often evaluated only in the advanced stages of the breeding program due to the amount of grain needed and the labor-intensive phenotyping assays. Advances in genomic resources have provided new tools to address the selection for these complex traits earlier in the breeding process. We used association mapping to identify key variants underlying various end-use quality traits and evaluate the usefulness of genomic prediction for these traits in hard red spring wheat from the Northern United States. A panel of 383 advanced breeding lines and cultivars representing the diversity of the University of Minnesota wheat breeding program was genotyped using the Illumina 90K single nucleotide polymorphism array and evaluated in multilocation trials using standard assessments of end-use quality. Sixty-three associations for grain or flour characteristics, mixograph, farinograph, and baking traits were identified. The majority of these associations were mapped in the vicinity of glutenin/gliadin or other known loci. In addition, a putative novel multi-trait association was identified on chromosome 6AL, and candidate gene analysis revealed eight genes of interest. Further, genomic prediction had a high predictive ability (PA) for mixograph and farinograph traits, with PA up to 0.62 and 0.50 in cross-validation and forward prediction, respectively. The deployment of 46 markers from GWAS to predict dough-rheology traits yielded low to moderate PA for various traits. The results of this study suggest that genomic prediction for end-use traits in early generations can be effective for mixograph and farinograph assays but not baking assays.

Genomic prediction of regional-scale performance in switchgrass (Panicum virgatum) by accounting for genotype-by-environment variation and yield surrogate traits

Switchgrass is a potential crop for bioenergy or carbon capture schemes, but further yield improvements through selective breeding are needed to encourage commercialization. To identify promising switchgrass germplasm for future breeding efforts, we conducted multisite and multitrait genomic prediction with a diversity panel of 630 genotypes from 4 switchgrass subpopulations (Gulf, Midwest, Coastal, and Texas), which were measured for spaced plant biomass yield across 10 sites. Our study focused on the use of genomic prediction to share information among traits and environments. Specifically, we evaluated the predictive ability of cross-validation (CV) schemes using only genetic data and the training set (cross-validation 1: CV1), a subset of the sites (cross-validation 2: CV2), and/or with 2 yield surrogates (flowering time and fall plant height). We found that genotype-by-environment interactions were largely due to the north-south distribution of sites. The genetic correlations between the yield surrogates and the biomass yield were generally positive (mean height r = 0.85; mean flowering time r = 0.45) and did not vary due to subpopulation or growing region (North, Middle, or South). Genomic prediction models had CV predictive abilities of -0.02 for individuals using only genetic data (CV1), but 0.55, 0.69, 0.76, 0.81, and 0.84 for individuals with biomass performance data from 1, 2, 3, 4, and 5 sites included in the training data (CV2), respectively. To simulate a resource-limited breeding program, we determined the predictive ability of models provided with the following: 1 site observation of flowering time (0.39); 1 site observation of flowering time and fall height (0.51); 1 site observation of fall height (0.52); 1 site observation of biomass (0.55); and 5 site observations of biomass yield (0.84). The ability to share information at a regional scale is very encouraging, but further research is required to accurately translate spaced plant biomass to commercial-scale sward biomass performance.

Multi-trait multi-environment genomic prediction of preliminary yield trial in pulse crop

Phenotypic selection of complex traits such as seed yield and protein in the preliminary yield trial (PYT) is often constrained by limited seed availability, resulting in trials with few environments and minimal to no replications. Multi-trait multi-environment enabled genomic prediction (MTME-GP) offers a valuable alternative to predict missing phenotypes of selection candidates for multiple traits and diverse environments. In this study, we assessed the efficiency of MTME-GP for improving seed protein and seed yield in field pea, the top two breeding targets but highly antagonistic traits in pulse crop. We utilized a set of 300 selection candidates in the PYT that virtually represented all possible families of the North Dakota State University field pea breeding program. Selection candidates were evaluated in three diverse, contrasting environments, as indicated by a range of heritability. Using whole- and split-environment cross validation schemes, MTME-GP had higher predictive ability than a standard additive G-BLUP model. Integrating a range of overlapping genotypes in between environments showed improvement on the predictive ability of the MTME-GP model but tends to plateau at 50%-80% training set size. Regardless of the cross-validation scheme, accuracy was among the lowest in stressed environments, presumably due to low heritability for seed protein and yield. This study provided insights into the potential of MTME-GP in a public pulse crop breeding program. The MTME-GP framework can be further improved with more testing environments and integration of additional orthogonal information in the early stages of the breeding pipeline.

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.

Multivariate analysis of phenotypic diversity elite bread wheat (Triticum aestivum L.) genotypes from ICARDA in Ethiopia

Wheat is an important crop for food security, providing a source of protein and energy for the growing population in Ethiopia. However, both biotic and abiotic factors limit national wheat productivity. The availability of genetically diverse wheat genotypes is crucial for developing new wheat varieties that are both high-yielding and resilient to stress. Therefore, this field trial aimed to assess phenotypic variation and relationship among ICARDA-derived bread wheat genotypes using multivariate analysis techniques. The trial was conducted at three locations: Enewari, Wogere, and Kulumsa using an alpha lattice design with two replications during the main cropping seasons of 2022 and 2023. Phenotypic data on eight agronomic traits and the severity of yellow rust were collected and R programming was used for data analysis. Individual and combined location data analysis of variance showed significant differences (p ≤ 0.05) among genotypes for most of the studied traits. The highest heritability and genetic advance as a percentage of the mean were observed in days to heading (90.8, 21.29), plant height (72.4, 28.6), seeds per spike (61.7, 28), thousand kernel weight (61.9, 12), and area under the disease progress curve (67, 39.8), suggesting a predominance of additive gene action. Grain yield showed a strong positive correlation with days to maturity, plant height, spike length, spikelet per spike, and thousand kernel weight for each location. Dendrogram and phylogenetic tree methods were used to group genotypes into four genetically distinct clusters. Cluster II and III had the greatest inter-cluster distance, indicating higher diversity among their genotypes. This study identified new candidate genotypes with superior agronomic performance, high grain yield traits, and robust resistance to yellow rust, making them valuable for both current and future wheat breeding programs. Additionally, the comprehensive dataset produced in this study could facilitate the identification of genetic variations influencing desirable traits through genome-wide association analysis.

Enhancing winter wheat prediction with genomics, phenomics and environmental data

In the realm of multi-environment prediction, when the goal is to predict a complete environment using the others as a training set, the efficiency of genomic selection (GS) falls short of expectations. Genotype by environment interaction poses a challenge in achieving high prediction accuracies. Consequently, current efforts are focused on enhancing efficiency by integrating various types of inputs, such as phenomics data, environmental information, and other omics data. In this study, we sought to evaluate the impact of incorporating environmental information into the modeling process, in addition to genomic and phenomics information. Our evaluation encompassed five data sets of soft white winter wheat, and the results revealed a significant improvement in prediction accuracy, as measured by the normalized root mean square error (NRMSE), through the integration of environmental information. Notably, there was an average gain in prediction accuracy of 49.19% in terms of NRMSE across the data sets. Moreover, the observed prediction accuracy ranged from 5.68% (data set 3) to 60.36% (data set 4), underscoring the substantial effect of integrating environmental information. By including genomic, phenomic, and environmental data in prediction models, plant breeding programs can improve selection efficiency across locations.

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.

Genomic selection in plant breeding: Key factors shaping two decades of progress

Genomic selection, the application of genomic prediction (GP) models to select candidate individuals, has significantly advanced in the past two decades, effectively accelerating genetic gains in plant breeding. This article provides a holistic overview of key factors that have influenced GP in plant breeding during this period. We delved into the pivotal roles of training population size and genetic diversity, and their relationship with the breeding population, in determining GP accuracy. Special emphasis was placed on optimizing training population size. We explored its benefits and the associated diminishing returns beyond an optimum size. This was done while considering the balance between resource allocation and maximizing prediction accuracy through current optimization algorithms. The density and distribution of single-nucleotide polymorphisms, level of linkage disequilibrium, genetic complexity, trait heritability, statistical machine-learning methods, and non-additive effects are the other vital factors. Using wheat, maize, and potato as examples, we summarize the effect of these factors on the accuracy of GP for various traits. The search for high accuracy in GP-theoretically reaching one when using the Pearson's correlation as a metric-is an active research area as yet far from optimal for various traits. We hypothesize that with ultra-high sizes of genotypic and phenotypic datasets, effective training population optimization methods and support from other omics approaches (transcriptomics, metabolomics and proteomics) coupled with deep-learning algorithms could overcome the boundaries of current limitations to achieve the highest possible prediction accuracy, making genomic selection an effective tool in plant breeding.

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.

Multi-trait and multi-environment genomic prediction for flowering traits in maize: a deep learning approach

Maize (Zea mays L.), the third most widely cultivated cereal crop in the world, plays a critical role in global food security. To improve the efficiency of selecting superior genotypes in breeding programs, researchers have aimed to identify key genomic regions that impact agronomic traits. In this study, the performance of multi-trait, multi-environment deep learning models was compared to that of Bayesian models (Markov Chain Monte Carlo generalized linear mixed models (MCMCglmm), Bayesian Genomic Genotype-Environment Interaction (BGGE), and Bayesian Multi-Trait and Multi-Environment (BMTME)) in terms of the prediction accuracy of flowering-related traits (Anthesis-Silking Interval: ASI, Female Flowering: FF, and Male Flowering: MF). A tropical maize panel of 258 inbred lines from Brazil was evaluated in three sites (Cambira-2018, Sabaudia-2018, and Iguatemi-2020 and 2021) using approximately 290,000 single nucleotide polymorphisms (SNPs). The results demonstrated a 14.4% increase in prediction accuracy when employing multi-trait models compared to the use of a single trait in a single environment approach. The accuracy of predictions also improved by 6.4% when using a single trait in a multi-environment scheme compared to using multi-trait analysis. Additionally, deep learning models consistently outperformed Bayesian models in both single and multiple trait and environment approaches. A complementary genome-wide association study identified associations with 26 candidate genes related to flowering time traits, and 31 marker-trait associations were identified, accounting for 37%, 37%, and 22% of the phenotypic variation of ASI, FF and MF, respectively. In conclusion, our findings suggest that deep learning models have the potential to significantly improve the accuracy of predictions, regardless of the approach used and provide support for the efficacy of this method in genomic selection for flowering-related traits in tropical maize.

Genomic prediction and selection response for grain yield in safflower

In plant breeding programs, multiple traits are recorded in each trial, and the traits are often correlated. Correlated traits can be incorporated into genomic selection models, especially for traits with low heritability, to improve prediction accuracy. In this study, we investigated the genetic correlation between important agronomic traits in safflower. We observed the moderate genetic correlations between grain yield (GY) and plant height (PH, 0.272-0.531), and low correlations between grain yield and days to flowering (DF, -0.157-0.201). A 4%-20% prediction accuracy improvement for grain yield was achieved when plant height was included in both training and validation sets with multivariate models. We further explored the selection responses for grain yield by selecting the top 20% of lines based on different selection indices. Selection responses for grain yield varied across sites. Simultaneous selection for grain yield and seed oil content (OL) showed positive gains across all sites with equal weights for both grain yield and oil content. Combining g×E interaction into genomic selection (GS) led to more balanced selection responses across sites. In conclusion, genomic selection is a valuable breeding tool for breeding high grain yield, oil content, and highly adaptable safflower varieties.

GWAS and genomic prediction for pre-harvest sprouting tolerance involving sprouting score and two other related traits in spring wheat

In wheat, a genome-wide association study (GWAS) and genomic prediction (GP) analysis were conducted for pre-harvest sprouting (PHS) tolerance and two of its related traits. For this purpose, an association panel of 190 accessions was phenotyped for PHS (using sprouting score), falling number, and grain color over two years and genotyped with 9904 DArTseq based SNP markers. GWAS for main-effect quantitative trait nucleotides (M-QTNs) using three different models (CMLM, SUPER, and FarmCPU) and epistatic QTNs (E-QTNs) using PLINK were performed. A total of 171 M-QTNs (CMLM, 47; SUPER, 70; FarmCPU, 54) for all three traits, and 15 E-QTNs involved in 20 first-order epistatic interactions were identified. Some of the above QTNs overlapped the previously reported QTLs, MTAs, and cloned genes, allowing delineating 26 PHS-responsive genomic regions that spread over 16 wheat chromosomes. As many as 20 definitive and stable QTNs were considered important for use in marker-assisted recurrent selection (MARS). The gene, TaPHS1, for PHS tolerance (PHST) associated with one of the QTNs was also validated using the KASP assay. Some of the M-QTNs were shown to have a key role in the abscisic acid pathway involved in PHST. Genomic prediction accuracies (based on the cross-validation approach) using three different models ranged from 0.41 to 0.55, which are comparable to the results of previous studies. In summary, the results of the present study improved our understanding of the genetic architecture of PHST and its related traits in wheat and provided novel genomic resources for wheat breeding based on MARS and GP.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01357-5.

Identification of leaf rust resistance loci in a geographically diverse panel of wheat using genome-wide association analysis

Leaf rust, caused by Puccinia triticina (Pt) is among the most devastating diseases posing a significant threat to global wheat production. The continuously evolving virulent Pt races in North America calls for exploring new sources of leaf rust resistance. A diversity panel of 365 bread wheat accessions selected from a worldwide population of landraces and cultivars was evaluated at the seedling stage against four Pt races (TDBJQ, TBBGS, MNPSD and, TNBJS). A wide distribution of seedling responses against the four Pt races was observed. Majority of the genotypes displayed a susceptible response with only 28 (9.8%), 59 (13.5%), 45 (12.5%), and 29 (8.1%) wheat accessions exhibiting a highly resistant response to TDBJQ, TBBGS, MNPSD and, TNBJS, respectively. Further, we conducted a high-resolution multi-locus genome-wide association study (GWAS) using a set of 302,524 high-quality single nucleotide polymorphisms (SNPs). The GWAS analysis identified 27 marker-trait associations (MTAs) for leaf rust resistance on different wheat chromosomes of which 20 MTAs were found in the vicinity of known Lr genes, MTAs, or quantitative traits loci (QTLs) identified in previous studies. The remaining seven significant MTAs identified represent genomic regions that harbor potentially novel genes for leaf rust resistance. Furthermore, the candidate gene analysis for the significant MTAs identified various genes of interest that may be involved in disease resistance. The identified resistant lines and SNPs linked to the QTLs in this study will serve as valuable resources in wheat rust resistance breeding programs.

Wheat end-use quality: State of art, genetics, genomics-assisted improvement, future challenges, and opportunities

Wheat is the most important source of food, feed, and nutrition for humans and livestock around the world. The expanding population has increasing demands for various wheat products with different quality attributes requiring the development of wheat cultivars that fulfills specific demands of end-users including millers and bakers in the international market. Therefore, wheat breeding programs continually strive to meet these quality standards by screening their improved breeding lines every year. However, the direct measurement of various end-use quality traits such as milling and baking qualities requires a large quantity of grain, traits-specific expensive instruments, time, and an expert workforce which limits the screening process. With the advancement of sequencing technologies, the study of the entire plant genome is possible, and genetic mapping techniques such as quantitative trait locus mapping and genome-wide association studies have enabled researchers to identify loci/genes associated with various end-use quality traits in wheat. Modern breeding techniques such as marker-assisted selection and genomic selection allow the utilization of these genomic resources for the prediction of quality attributes with high accuracy and efficiency which speeds up crop improvement and cultivar development endeavors. In addition, the candidate gene approach through functional as well as comparative genomics has facilitated the translation of the genomic information from several crop species including wild relatives to wheat. This review discusses the various end-use quality traits of wheat, their genetic control mechanisms, the use of genetics and genomics approaches for their improvement, and future challenges and opportunities for wheat breeding.

Multi-trait genomic prediction improves selection accuracy for enhancing seed mineral concentrations in pea

Multi-trait genomic selection (MT-GS) has the potential to improve predictive ability by maximizing the use of information across related genotypes and genetically correlated traits. In this study, we extended the use of sparse phenotyping method into the MT-GS framework by split testing of entries to maximize borrowing of information across genotypes and predict missing phenotypes for targeted traits without additional phenotyping expenditure. Using 300 advanced breeding lines from North Dakota State University (NDSU) pulse breeding program and ∼200 USDA accessions that were evaluated for 10 nutritional traits, our results show that the proposed sparse phenotyping aided MT-GS can further improve predictive ability by >12% across traits compared with univariate (UNI) genomic selection. The proposed strategy departed from the previous reports that weak genetic correlation is a limitation to the advantage of MT-GS over UNI genomic selection, which was evident in the partially balanced phenotyping-enabled MT-GS. Our results point to heritability and genetic correlation between traits as possible metrics to optimize and further improve the estimation of model parameters, and ultimately, prediction performance. Overall, our study offers a new approach to optimize the prediction performance using the MT-GS and further highlight strategy to maximize the efficiency of GS in a plant breeding program. The sparse-testing-aided MT-GS proposed in this study can be further extended to multi-environment, multi-trait GS to improve prediction performance and further reduce the cost of phenotyping and time-consuming data collection process.

Utilizing genomics and historical data to optimize gene pools for new breeding programs: A case study in winter wheat

With the rapid generation and preservation of both genomic and phenotypic information for many genotypes within crops and across locations, emerging breeding programs have a valuable opportunity to leverage these resources to 1) establish the most appropriate genetic foundation at program inception and 2) implement robust genomic prediction platforms that can effectively select future breeding lines. Integrating genomics-enabled¹ breeding into cultivar development can save costs and allow resources to be reallocated towards advanced (i.e., later) stages of field evaluation, which can facilitate an increased number of testing locations and replicates within locations. In this context, a reestablished winter wheat breeding program was used as a case study to understand best practices to leverage and tailor existing genomic and phenotypic resources to determine optimal genetics for a specific target population of environments. First, historical multi-environment phenotype data, representing 1,285 advanced breeding lines, were compiled from multi-institutional testing as part of the SunGrains cooperative and used to produce GGE biplots and PCA for yield. Locations were clustered based on highly correlated line performance among the target population of environments into 22 subsets. For each of the subsets generated, EMMs and BLUPs were calculated using linear models with the ‘lme4’ R package. Second, for each subset, TPs representative of the new SC breeding lines were determined based on genetic relatedness using the ‘STPGA’ R package. Third, for each TP, phenotypic values and SNP data were incorporated into the ‘rrBLUP’ mixed models for generation of GEBVs of YLD, TW, HD and PH. Using a five-fold cross-validation strategy, an average accuracy of r = 0.42 was obtained for yield between all TPs. The validation performed with 58 SC elite breeding lines resulted in an accuracy of r = 0.62 when the TP included complete historical data. Lastly, QTL-by-environment interaction for 18 major effect genes across three geographic regions was examined. Lines harboring major QTL in the absence of disease could potentially underperform (e.g., Fhb1 R-gene), whereas it is advantageous to express a major QTL under biotic pressure (e.g., stripe rust R-gene). This study highlights the importance of genomics-enabled breeding and multi-institutional partnerships to accelerate cultivar development.

Ascertaining yield and grain protein content stability in wheat genotypes having the Gpc-B1 gene using univariate, multivariate, and correlation analysis

The high performance and stability of wheat genotypes for yield, grain protein content (GPC), and other desirable traits are critical for varietal development and food and nutritional security. Likewise, the genotype by environment (G × E) interaction (GEI) should be thoroughly investigated and favorably utilized whenever genotype selection decisions are made. The present study was planned with the following two major objectives: 1) determination of GEI for some advanced wheat genotypes across four locations (Ludhiana, Ballowal, Patiala, and Bathinda) of Punjab, India; and 2) selection of the best genotypes with high GPC and yield in various environments. Different univariate [Eberhart and Ruessll's models; Perkins and Jinks' models; Wrike's Ecovalence; and Francis and Kannenberg's models], multivariate (AMMI and GGE biplot), and correlation analyses were used to interpret the data from the multi-environmental trial (MET). Consequently, both the univariate and multivariate analyses provided almost similar results regarding the top-performing and stable genotypes. The analysis of variance revealed that variation due to environment, genotype, and GEI was highly significant at the 0.01 and 0.001 levels of significance for all studied traits. The days to flowering, plant height, spikelets per spike, grain per spike, days to maturity, and 1000-grain weight were specifically affected by the environment, whereas yield was mainly affected by the environment and GEI. Genotypes, on the other hand, had a greater impact on the GPC than environmental conditions. As a result, a multi-environmental investigation was necessary to identify the GEI for wheat genotype selection because the GEI was very significant for all of the evaluated traits. Yield, 1000-grain weight, spikelet per spike, and days to maturity were observed to have positive correlations, implying the feasibility of their simultaneous selection for yield enhancement. However, GPC was observed to have a negative correlation with yield. Patiala was found to be the most discriminating environment for both yield and GPC and also the most effective representative environment for GPC, whereas Ludhiana was found to be the most effective representative environment for yield. Eventually, two NILs (BWL7508, and BWL7511) were selected as the top across all environments for both yield and GPC.

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.

Consensus genomic regions associated with multiple abiotic stress tolerance in wheat and implications for wheat breeding

In wheat, a meta-analysis was performed using previously identified QTLs associated with drought stress (DS), heat stress (HS), salinity stress (SS), water-logging stress (WS), pre-harvest sprouting (PHS), and aluminium stress (AS) which predicted a total of 134 meta-QTLs (MQTLs) that involved at least 28 consistent and stable MQTLs conferring tolerance to five or all six abiotic stresses under study. Seventy-six MQTLs out of the 132 physically anchored MQTLs were also verified with genome-wide association studies. Around 43% of MQTLs had genetic and physical confidence intervals of less than 1 cM and 5 Mb, respectively. Consequently, 539 genes were identified in some selected MQTLs providing tolerance to 5 or all 6 abiotic stresses. Comparative analysis of genes underlying MQTLs with four RNA-seq based transcriptomic datasets unravelled a total of 189 differentially expressed genes which also included at least 11 most promising candidate genes common among different datasets. The promoter analysis showed that the promoters of these genes include many stress responsiveness cis-regulatory elements, such as ARE, MBS, TC-rich repeats, As-1 element, STRE, LTR, WRE3, and WUN-motif among others. Further, some MQTLs also overlapped with as many as 34 known abiotic stress tolerance genes. In addition, numerous ortho-MQTLs among the wheat, maize, and rice genomes were discovered. These findings could help with fine mapping and gene cloning, as well as marker-assisted breeding for multiple abiotic stress tolerances in wheat.

Combining NDVI and Bacterial Blight Score to Predict Grain Yield in Field Pea

Field pea is the most commonly grown temperate pulse crop, with close to 15 million tons produced globally in 2020. Varieties improved through breeding are important to ensure ongoing improvements in yield and disease resistance. Genomic selection (GS) is a modern breeding approach that could substantially improve the rate of genetic gain for grain yield, and its deployment depends on the prediction accuracy (PA) that can be achieved. In our study, four yield trials representing breeding lines' advancement stages of the breeding program (S0, S1, S2, and S3) were assessed with grain yield, aerial high-throughput phenotyping (normalized difference vegetation index, NDVI), and bacterial blight disease scores (BBSC). Low-to-moderate broad-sense heritability (0.31-0.71) and narrow-sense heritability (0.13-0.71) were observed, as the estimated additive and non-additive genetic components for the three traits varied with the different models fitted. The genetic correlations among the three traits were high, particularly in the S0-S2 stages. NDVI and BBSC were combined to investigate the PA for grain yield by univariate and multivariate GS models, and multivariate models showed higher PA than univariate models in both cross-validation and forward prediction methods. A 6-50% improvement in PA was achieved when multivariate models were deployed. The highest PA was indicated in the forward prediction scenario when the training population consisted of early generation breeding stages with the multivariate models. Both NDVI and BBSC are commonly used traits that could be measured in the early growth stage; however, our study suggested that NDVI is a more useful trait to predict grain yield with high accuracy in the field pea breeding program, especially in diseased trials, through its incorporation into multivariate models.

Multi-Trait Multi-Environment Genomic Prediction for End-Use Quality Traits in Winter Wheat

Soft white wheat is a wheat class used in foreign and domestic markets to make various end products requiring specific quality attributes. Due to associated cost, time, and amount of seed needed, phenotyping for the end-use quality trait is delayed until later generations. Previously, we explored the potential of using genomic selection (GS) for selecting superior genotypes earlier in the breeding program. Breeders typically measure multiple traits across various locations, and it opens up the avenue for exploring multi-trait-based GS models. This study's main objective was to explore the potential of using multi-trait GS models for predicting seven different end-use quality traits using cross-validation, independent prediction, and across-location predictions in a wheat breeding program. The population used consisted of 666 soft white wheat genotypes planted for 5 years at two locations in Washington, United States. We optimized and compared the performances of four uni-trait- and multi-trait-based GS models, namely, Bayes B, genomic best linear unbiased prediction (GBLUP), multilayer perceptron (MLP), and random forests. The prediction accuracies for multi-trait GS models were 5.5 and 7.9% superior to uni-trait models for the within-environment and across-location predictions. Multi-trait machine and deep learning models performed superior to GBLUP and Bayes B for across-location predictions, but their advantages diminished when the genotype by environment component was included in the model. The highest improvement in prediction accuracy, that is, 35% was obtained for flour protein content with the multi-trait MLP model. This study showed the potential of using multi-trait-based GS models to enhance prediction accuracy by using information from previously phenotyped traits. It would assist in speeding up the breeding cycle time in a cost-friendly manner.

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.