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.

Profiling of Fusarium head blight resistance QTL haplotypes through molecular markers, genotyping-by-sequencing, and machine learning.

KEY MESSAGE: Marker-assisted selection is important for cultivar development. We propose a system where a training population genotyped for QTL and genome-wide markers may predict QTL haplotypes in early development germplasm. Breeders screen germplasm with molecular markers to identify and select individuals that have desirable haplotypes. The objective of this research was to investigate whether QTL haplotypes can be accurately predicted using SNPs derived by genotyping-by-sequencing (GBS). In the SunGrains program during 2020 (SG20) and 2021 (SG21), 1,536 and 2,352 lines submitted for GBS were genotyped with markers linked to the Fusarium head blight QTL: Qfhb.nc-1A, Qfhb.vt-1B, Fhb1, and Qfhb.nc-4A. In parallel, data were compiled from the 2011-2020 Southern Uniform Winter Wheat Scab Nursery (SUWWSN), which had been screened for the same QTL, sequenced via GBS, and phenotyped for: visual Fusarium severity rating (SEV), percent Fusarium damaged kernels (FDK), deoxynivalenol content (DON), plant height, and heading date. Three machine learning models were evaluated: random forest, k-nearest neighbors, and gradient boosting machine. Data were randomly partitioned into training-testing splits. The QTL haplotype and 100 most correlated GBS SNPs were used for training and tuning of each model. Trained machine learning models were used to predict QTL haplotypes in the testing partition of SG20, SG21, and the total SUWWSN. Mean disease ratings for the observed and predicted QTL haplotypes were compared in the SUWWSN. For all models trained using the SG20 and SG21, the observed Fhb1 haplotype estimated group means for SEV, FDK, DON, plant height, and heading date in the SUWWSN were not significantly different from any of the predicted Fhb1 calls. This indicated that machine learning may be utilized in breeding programs to accurately predict QTL haplotypes in earlier generations.

Training population selection and use of fixed effects to optimize genomic predictions in a historical USA winter wheat panel.

KEY MESSAGE: The optimization of training populations and the use of diagnostic markers as fixed effects increase the predictive ability of genomic prediction models in a cooperative wheat breeding panel. Plant breeding programs often have access to a large amount of historical data that is highly unbalanced, particularly across years. This study examined approaches to utilize these data sets as training populations to integrate genomic selection into existing pipelines. We used cross-validation to evaluate predictive ability in an unbalanced data set of 467 winter wheat (Triticum aestivum L.) genotypes evaluated in the Gulf Atlantic Wheat Nursery from 2008 to 2016. We evaluated the impact of different training population sizes and training population selection methods (Random, Clustering, PEVmean and PEVmean1) on predictive ability. We also evaluated inclusion of markers associated with major genes as fixed effects in prediction models for heading date, plant height, and resistance to powdery mildew (caused by Blumeria graminis f. sp. tritici). Increases in predictive ability as the size of the training population increased were more evident for Random and Clustering training population selection methods than for PEVmean and PEVmean1. The selection methods based on minimization of the prediction error variance (PEV) outperformed the Random and Clustering methods across all the population sizes. Major genes added as fixed effects always improved model predictive ability, with the greatest gains coming from combinations of multiple genes. Maximum predictabilities among all prediction methods were 0.64 for grain yield, 0.56 for test weight, 0.71 for heading date, 0.73 for plant height, and 0.60 for powdery mildew resistance. Our results demonstrate the utility of combining unbalanced phenotypic records with genome-wide SNP marker data for predicting the performance of untested genotypes.

Functionality of gliadin proteins in wheat flour tortillas.

Gliadins are monomeric proteins that are encoded by the genes at the loci Gli 1 and Gli 2 present on the short arm of homologous wheat chromosomes 1 and 6, respectively. Studies have suggested that gliadins may play an important role in determining the functional properties of wheat flour. The main objective of this study was to understand the functionality of gliadins with respect to tortilla quality. The important tortilla quality attributes are diameter, opacity, and shelf stability, designated here as rollability or the ability to roll or fold the tortilla without cracking. In this study gliadin functionality in tortilla quality was studied using near-isogenic wheat lines that have deletions in either Gli A1, Gli D1, Gli A2, or Gli D2 gliadin loci. The deletion lines are designated by the same abbreviations. Dough and tortillas were prepared from the parent line used to derive these deletion lines, each individual deletion line, and a control commercial tortilla flour. Quantitative and qualitative evaluations were performed on the dough and tortillas derived from the flour from each of these lines. None of the deletions in the gliadin loci altered the shelf stability versus that found for the parent to the deletion lines or control tortilla flour. However, deletions in the Gli 2 loci, in particular Gli A2 reduced the relative proportion of alpha- and beta-gliadins with a greater cysteine amino acid content and gluten cross-link function versus the chain-terminating omega-gliadins in Gli 1, which were still present. As such, the dough and gluten matrix appeared to have greater extensibility, which improved the diameter and overall quality of the tortillas while not altering the rollability. Deletions in the Gli 1 loci had the opposite result with increased cross-linking of alpha- and beta-gliadins, polymeric protein content, and a stronger dough that decreased the diameter and overall quality of the tortillas. The data suggest that altering certain Gli 2 loci through null alleles could be a viable strategy to develop cultivars improved for the specific functionality requirements needed for the rapidly growing tortilla market.