Genomic Prediction from Multi-Environment Trials of Wheat Breeding.

Genomic prediction relates a set of markers to variability in observed phenotypes of cultivars and allows for the prediction of phenotypes or breeding values of genotypes on unobserved individuals. Most genomic prediction approaches predict breeding values based solely on additive effects. However, the economic value of wheat lines is not only influenced by their additive component but also encompasses a non-additive part (e.g., additive × additive epistasis interaction). In this study, genomic prediction models were implemented in three target populations of environments (TPE) in South Asia. Four models that incorporate genotype × environment interaction (G × E) and genotype × genotype (GG) were tested: Factor Analytic (FA), FA with genomic relationship matrix (FA + G), FA with epistatic relationship matrix (FA + GG), and FA with both genomic and epistatic relationship matrices (FA + G + GG). Results show that the FA + G and FA + G + GG models displayed the best and a similar performance across all tests, leading us to infer that the FA + G model effectively captures certain epistatic effects. The wheat lines tested in sites in different TPE were predicted with different precisions depending on the cross-validation employed. In general, the best prediction accuracy was obtained when some lines were observed in some sites of particular TPEs and the worse genomic prediction was observed when wheat lines were never observed in any site of one TPE.

Deep learning methods improve genomic prediction of wheat breeding.

In the field of plant breeding, various machine learning models have been developed and studied to evaluate the genomic prediction (GP) accuracy of unseen phenotypes. Deep learning has shown promise. However, most studies on deep learning in plant breeding have been limited to small datasets, and only a few have explored its application in moderate-sized datasets. In this study, we aimed to address this limitation by utilizing a moderately large dataset. We examined the performance of a deep learning (DL) model and compared it with the widely used and powerful best linear unbiased prediction (GBLUP) model. The goal was to assess the GP accuracy in the context of a five-fold cross-validation strategy and when predicting complete environments using the DL model. The results revealed the DL model outperformed the GBLUP model in terms of GP accuracy for two out of the five included traits in the five-fold cross-validation strategy, with similar results in the other traits. This indicates the superiority of the DL model in predicting these specific traits. Furthermore, when predicting complete environments using the leave-one-environment-out (LOEO) approach, the DL model demonstrated competitive performance. It is worth noting that the DL model employed in this study extends a previously proposed multi-modal DL model, which had been primarily applied to image data but with small datasets. By utilizing a moderately large dataset, we were able to evaluate the performance and potential of the DL model in a context with more information and challenging scenario in plant breeding.

A marker weighting approach for enhancing within-family accuracy in genomic prediction.

Genomic selection is revolutionizing plant breeding. However, its practical implementation is still very challenging, since predicted values do not necessarily have high correspondence to the observed phenotypic values. When the goal is to predict within-family, it is not always possible to obtain reasonable accuracies, which is of paramount importance to improve the selection process. For this reason, in this research, we propose the Adversaria-Boruta (AB) method, which combines the virtues of the adversarial validation (AV) method and the Boruta feature selection method. The AB method operates primarily by minimizing the disparity between training and testing distributions. This is accomplished by reducing the weight assigned to markers that display the most significant differences between the training and testing sets. Therefore, the AB method built a weighted genomic relationship matrix that is implemented with the genomic best linear unbiased predictor (GBLUP) model. The proposed AB method is compared using 12 real data sets with the GBLUP model that uses a nonweighted genomic relationship matrix. Our results show that the proposed AB method outperforms the GBLUP by 8.6, 19.7, and 9.8% in terms of Pearson's correlation, mean square error, and normalized root mean square error, respectively. Our results support that the proposed AB method is a useful tool to improve the prediction accuracy of a complete family, however, we encourage other investigators to evaluate the AB method to increase the empirical evidence of its potential.

Simultaneous improvement of grain yield and grain protein concentration in durum wheat by using association tests and weighted GBLUP.

KEY MESSAGE: Simultaneous improvement for GY and GPC by using GWAS and GBLUP suggested a significant application in durum wheat breeding. Despite the importance of grain protein concentration (GPC) in determining wheat quality, its negative correlation with grain yield (GY) is still one of the major challenges for breeders. Here, a durum wheat panel of 200 genotypes was evaluated for GY, GPC, and their derived indices (GPD and GYD), under eight different agronomic conditions. The plant material was genotyped with the Illumina 25 k iSelect array, and a genome-wide association study was performed. Two statistical models revealed dozens of marker-trait associations (MTAs), each explaining up to 30%. phenotypic variance. Two markers on chromosomes 2A and 6B were consistently identified by both models and were found to be significantly associated with GY and GPC. MTAs identified for phenological traits co-mapped to well-known genes (i.e., Ppd-1, Vrn-1). The significance values (p-values) that measure the strength of the association of each single nucleotide polymorphism marker with the target traits were used to perform genomic prediction by using a weighted genomic best linear unbiased prediction model. The trained models were ultimately used to predict the agronomic performances of an independent durum wheat panel, confirming the utility of genomic prediction, although environmental conditions and genetic backgrounds may still be a challenge to overcome. The results generated through our study confirmed the utility of GPD and GYD to mitigate the inverse GY and GPC relationship in wheat, provided novel markers for marker-assisted selection and opened new ways to develop cultivars through genomic prediction approaches.

Multivariate Genomic Hybrid Prediction with Kernels and Parental Information.

Genomic selection (GS) plays a pivotal role in hybrid prediction. It can enhance the selection of parental lines, accurately predict hybrid performance, and harness hybrid vigor. Likewise, it can optimize breeding strategies by reducing field trial requirements, expediting hybrid development, facilitating targeted trait improvement, and enhancing adaptability to diverse environments. Leveraging genomic information empowers breeders to make informed decisions and significantly improve the efficiency and success rate of hybrid breeding programs. In order to improve the genomic ability performance, we explored the incorporation of parental phenotypic information as covariates under a multi-trait framework. Approach 1, referred to as Pmean, directly utilized parental phenotypic information without any preprocessing. While approach 2, denoted as BV, replaced the direct use of phenotypic values of both parents with their respective breeding values. While an improvement in prediction performance was observed in both approaches, with a minimum 4.24% reduction in the normalized root mean square error (NRMSE), the direct incorporation of parental phenotypic information in the Pmean approach slightly outperformed the BV approach. We also compared these two approaches using linear and nonlinear kernels, but no relevant gain was observed. Finally, our results increase empirical evidence confirming that the integration of parental phenotypic information helps increase the prediction performance of hybrids.

A novel method for genomic-enabled prediction of cultivars in new environments.

Introduction: Genomic selection (GS) has gained global importance due to its potential to accelerate genetic progress and improve the efficiency of breeding programs. Objectives of the research: In this research we proposed a method to improve the prediction accuracy of tested lines in new (untested) environments. Method-1: The new method trained the model with a modified response variable (a difference of response variables) that decreases the lack of a non-stationary distribution between the training and testing and improved the prediction accuracy. Comparing new and conventional method: We compared the prediction accuracy of the conventional genomic best linear unbiased prediction (GBLUP) model (M1) including (or not) genotype × environment interaction (GE) (M1_GE; M1_NO_GE) versus the proposed method (M2) on several data sets. Results and discussion: The gain in prediction accuracy of M2, versus M1_GE, M1_NO_GE in terms of Pearson´s correlation was of at least 4.3%, while in terms of percentage of top-yielding lines captured when was selected the 10% (Best10) and 20% (Best20) of lines was at least of 19.5%, while in terms of Normalized Root Mean Squared Error (NRMSE) was of at least of 42.29%.

Unlocking the molecular basis of wheat straw composition and morphological traits through multi-locus GWAS.

BACKGROUND: Rapid reductions in emissions from fossil fuel burning are needed to curb global climate change. Biofuel production from crop residues can contribute to reducing the energy crisis and environmental deterioration. Wheat is a renewable source for biofuels owing to the low cost and high availability of its residues. Thus, identifying candidate genes controlling these traits is pivotal for efficient biofuel production. Here, six multi-locus genome-wide association (ML-GWAS) models were applied using 185 tetraploid wheat accessions to detect quantitative trait nucleotides (QTNs) for fifteen traits associated with biomass composition. RESULTS: Among the 470 QTNs, only 72 identified by at least two models were considered as reliable. Among these latter, 16 also showed a significant effect on the corresponding trait (p.value < 0.05). Candidate genes survey carried out within 4 Mb flanking the QTNs, revealed putative biological functions associated with lipid transfer and metabolism, cell wall modifications, cell cycle, and photosynthesis. Four genes encoded as Cellulose Synthase (CeSa), Anaphase promoting complex (APC/C), Glucoronoxylan 4-O Methyltransferase (GXM) and HYPONASTIC LEAVES1 (HYL1) might be responsible for an increase in cellulose, and natural and acid detergent fiber (NDF and ADF) content in tetraploid wheat. In addition, the SNP marker RFL_Contig3228_2154 associated with the variation in stem solidness (Q.Scsb-3B) was validated through two molecular methods (High resolution melting; HRM and RNase H2-dependent PCR; rhAMP). CONCLUSIONS: The study provides new insights into the genetic basis of biomass composition traits on tetraploid wheat. The application of six ML-GWAS models on a panel of diverse wheat genotypes represents an efficient approach to dissect complex traits with low heritability such as wheat straw composition. The discovery of genes/genomic regions associated with biomass production and straw quality parameters is expected to accelerate the development of high-yielding wheat varieties useful for biofuel production.

QTL Analysis of Five Morpho-Physiological Traits in Bread Wheat Using Two Mapping Populations Derived from Common Parents.

Traits such as plant height (PH), juvenile growth habit (GH), heading date (HD), and tiller number are important for both increasing yield potential and improving crop adaptation to climate change. In the present study, these traits were investigated by using the same bi-parental population at early (F2 and F2-derived F3 families) and late (F6 and F7, recombinant inbred lines, RILs) generations to detect quantitative trait loci (QTLs) and search for candidate genes. A total of 176 and 178 lines were genotyped by the wheat Illumina 25K Infinium SNP array. The two genetic maps spanned 2486.97 cM and 3732.84 cM in length, for the F2 and RILs, respectively. QTLs explaining the highest phenotypic variation were found on chromosomes 2B, 2D, 5A, and 7D for HD and GH, whereas those for PH were found on chromosomes 4B and 4D. Several QTL detected in the early generations (i.e., PH and tiller number) were not detected in the late generations as they were due to dominance effects. Some of the identified QTLs co-mapped to well-known adaptive genes (i.e., Ppd-1, Vrn-1, and Rht-1). Other putative candidate genes were identified for each trait, of which PINE1 and PIF4 may be considered new for GH and TTN in wheat. The use of a large F2 mapping population combined with NGS-based genotyping techniques could improve map resolution and allow closer QTL tagging.