Genomic Selection Models Research Articles

Statistical Methods for the Quantitative Genetic Analysis of High-Throughput Phenotyping Data.

The advent of plant phenomics, coupled with the wealth of genotypic data generated by next-generation sequencing technologies, provides exciting new resources for investigations into and improvement of complex traits. However, these new technologies also bring new challenges in quantitative genetics, namely, a need for the development of robust frameworks that can accommodate these high-dimensional data. In this chapter, we describe methods for the statistical analysis of high-throughput phenotyping (HTP) data with the goal of enhancing the prediction accuracy of genomic selection (GS). Following the Introduction in Sec. 1, Sec.2 discusses field-based HTP, including the use of unoccupied aerial vehicles and light detection and ranging, as well as how we can achieve increased genetic gain by utilizing image data derived from HTP. Section 3 considers extending commonly used GS models to integrate HTP data as covariates associated with the principal trait response, such as yield. Particular focus is placed on single-trait, multi-trait, and genotype by environment interaction models. One unique aspect of HTP data is that phenomics platforms often produce large-scale data with high spatial and temporal resolution for capturing dynamic growth, development, and stress responses. Section 4 discusses the utility of a random regression model for performing longitudinal modeling. The chapter concludes with a discussion of some standing issues.

Overview of Genomic Prediction Methods and the Associated Assumptions on the Variance of Marker Effect, and on the Architecture of the Target Trait.

Genomic selection (GS) is a methodology that revolutionized the process of breeding improved genetic materials in plant and animal breeding programs. It uses predicted genomic values of the potential of untested/unobserved genotypes as surrogates of phenotypes during the selection process. Such that the predicted genomic values are obtained using exclusively the marker profiles of the untested genotypes, and these potentially can be used by breeders for screening the genotypes to be advanced in the breeding pipeline, to identify potential parents for next improvement cycles, or to find optimal crosses for targeting genotypes among others. Conceptually, GS initially requires a set of genotypes with both molecular marker information and phenotypic data for model calibration and then the performance of untested genotypes is predicted using their marker profiles only. Hence, it is expected that breeders would look at these values in order to conduct selections. Even though the concept of GS seems trivial, due to the high dimensional nature of the data delivered from modern sequencing technologies where the number of molecular markers (p) excess by far the number of data points available for model fitting (n; p≫n) a complete renovated set of prediction models was needed to cope with this challenge. In this chapter, we provide a conceptual framework for comparing statistical models to overcome the "large p, small n problem." Given the very large diversity of GS models only the most popular are presented here; mainly we focused on linear regression-based models and nonparametric models that predict the genetic estimated breeding values (GEBV) in a single environment considering a single trait only, mainly in the context of plant breeding.

Genome and EnvironmentBased Prediction Models and Methods of Complex Traits Incorporating Genotype × Environment Interaction.

Genomic-enabled prediction models are of paramount importance for the successful implementation of genomic selection (GS) based on breeding values. As opposed to animal breeding, plant breeding includes extensive multienvironment and multiyear field trial data. Hence, genomic-enabled prediction models should include genotype × environment (G×E) interaction, which most of the time increases the prediction performance when the response of lines are different from environment to environment. In this chapter, we describe a historical timeline since 2012 related to advances of the GS models that take into account G×E interaction. We describe theoretical and practical aspects of those GS models, including the gains in prediction performance when including G×E structures for both complex continuous and categorical scale traits. Then, we detailed and explained the main G×E genomic prediction models for complex traits measured in continuous and noncontinuous (categorical) scale. Related to G×E interaction models this review also examine the analyses of the information generated with high-throughput phenotype data (phenomic) and the joint analyses of multitrait and multienvironment field trial data that is also employed in the general assessment of multitrait G×E interaction. The inclusion of nongenomic data in increasing the accuracy and biological reliability of the G×E approach is also outlined. We show the recent advances in large-scale envirotyping (enviromics), and how the use of mechanistic computational modeling can derive the crop growth and development aspects useful for predicting phenotypes and explaining G×E.

Comprehensive evaluation of mapping complex traits in wheat using genome-wide association studies.

Genome-wide association studies (GWAS) are effectively applied to detect the marker trait associations (MTAs) using whole genome-wide variants for complex quantitative traits in different crop species. GWAS has been applied in wheat for different quality, biotic and abiotic stresses, and agronomic and yield-related traits. Predictions for marker-trait associations are controlled with the development of better statistical models taking population structure and familial relatedness into account. In this review, we have provided a detailed overview of the importance of association mapping, population design, high-throughput genotyping and phenotyping platforms, advancements in statistical models and multiple threshold comparisons, and recent GWA studies conducted in wheat. The information about MTAs utilized for gene characterization and adopted in breeding programs is also provided. In the literature that we surveyed, as many as 86,122 wheat lines have been studied under various GWA studies reporting 46,940 loci. However, further utilization of these is largely limited. The future breakthroughs in area of genomic selection, multi-omics-based approaches, machine, and deep learning models in wheat breeding after exploring the complex genetic structure with the GWAS are also discussed. This is a most comprehensive study of a large number of reports on wheat GWAS and gives a comparison and timeline of technological developments in this area. This will be useful to new researchers or groups who wish to invest in GWAS.

Solution of Mixed-Integer Optimization Problems in Bioinformatics with Differential Evolution Method

The solution of the so-called mixed-integer optimization problem is an important challenge for modern life sciences. A wide range of methods has been developed for its solution, including metaheuristics approaches. Here, a modification is proposed of the differential evolution entirely parallel (DEEP) method introduced recently that was successfully applied to mixed-integer optimization problems. The triangulation recombination rule was implemented and the recombination coefficients were included in the evolution process in order to increase the robustness of the optimization. The deduplication step included in the procedure ensures the uniqueness of individual integer-valued parameters in the solution vectors. The developed algorithms were implemented in the DEEP software package and applied to three bioinformatic problems. The application of the method to the optimization of predictors set in the genomic selection model in wheat resulted in dimensionality reduction such that the phenotype can be predicted with acceptable accuracy using a selected subset of SNP markers. The method was also successfully used to optimize the training set of samples for such a genomic selection model. According to the obtained results, the developed algorithm was capable of constructing a non-linear phenomenological regression model of gene expression in developing a Drosophila eye with almost the same average accuracy but significantly less standard deviation than the linear models obtained earlier.

Genomic Selection and Genome-Wide Association Studies for Grain Protein Content Stability in a Nested Association Mapping Population of Wheat

Grain protein content (GPC) is controlled by complex genetic systems and their interactions and is an important quality determinant for hard spring wheat as it has a positive effect on bread and pasta quality. GPC is variable among genotypes and strongly influenced by the environment. Thus, understanding the genetic control of wheat GPC and identifying genotypes with improved stability is an important breeding goal. The objectives of this research were to identify genetic backgrounds with less variation for GPC across environments and identify quantitative trait loci (QTLs) controlling the stability of GPC. A spring wheat nested association mapping (NAM) population of 650 recombinant inbred lines (RIL) derived from 26 diverse founder parents crossed to one common parent, ‘Berkut’, was phenotyped over three years of field trials (2014–2016). Genomic selection models were developed and compared based on predictions of GPC and GPC stability. After observing variable genetic control of GPC within the NAM population, seven RIL families displaying reduced marker-by-environment interaction were selected based on a stability index derived from a Finlay–Wilkinson regression. A genome-wide association study identified eighteen significant QTLs for GPC stability with a Bonferroni-adjusted p-value < 0.05 using four different models and out of these eighteen QTLs eight were identified by two or more GWAS models simultaneously. This study also demonstrated that genome-wide prediction of GPC with ridge regression best linear unbiased estimates reached up to r = 0.69. Genomic selection can be used to apply selection pressure for GPC and improve genetic gain for GPC.

Genotyping-by-sequencing and genomic selection applications in hexaploid triticale.

Triticale, a hybrid species between wheat and rye, is one of the newest additions to the plant kingdom with a very short history of improvement. It has very limited genomic resources because of its large and complex genome. Objectives of this study were to generate dense marker data, understand genetic diversity, population structure, linkage disequilibrium (LD), and estimate accuracies of commonly used genomic selection (GS) models on forage yield of triticale. Genotyping-by-sequencing (GBS), using PstI and MspI restriction enzymes for reducing genome complexity, was performed on a triticale diversity panel (n = 289). After filtering for biallelic loci with more than 70% genome coverage, and minor allele frequency (MAF) > 0.05, de novo variant calling identified 16,378 single nucleotide polymorphism (SNP) markers. Sequences of these variants were mapped to wheat and rye reference genomes to infer their homologous groups and chromosome positions. About 45% (7430), and 58% (9500) of the de novo identified SNPs were mapped to the wheat and rye reference genomes, respectively. Interestingly, 28.9% (2151) of the 7430 SNPs were mapped to the D genome of hexaploid wheat, indicating substantial substitution of the R genome with D genome in cultivated triticale. About 27% of marker pairs were in significant LD with an average r2 > 0.18 (P < 0.05). Genome-wide LD declined rapidly to r2 < 0.1 beyond 10 kb physical distance. The three sub-genomes (A, B, and R) showed comparable LD decay patterns. Genetic diversity and population structure analyses identified five distinct clusters. Genotype grouping did not follow prior winter vs spring-type classification. However, one of the clusters was largely dominated by winter triticale. GS accuracies were estimated for forage yield using three commonly used models with different training population sizes and marker densities. GS accuracy increased with increasing training population size while gain in accuracy tended to plateau with marker densities of 2000 SNPs or more. Average GS accuracy was about 0.52, indicating the potential of using GS in triticale forage yield improvement.

Genomic prediction with allele dosage information in highly polyploid species.

Including allele, dosage can improve genomic selection in highly polyploid species under higher frequency of different heterozygous genotypic classes and high dominance degree levels. Several studies have shown how to leverage allele dosage information to improve the accuracy of genomic selection models in autotetraploid. In this study, we expanded the methodology used for genomic selection in autotetraploid to higher (and mixed) ploidy levels. We adapted the models to build covariance matrices of both additive and digenic dominance effects that are subsequently used in genomic selection models. We applied these models using estimates of ploidy and allele dosage to sugarcane and sweet potato datasets and validated our results by also applying the models in simulated data. For the simulated datasets, including allele dosage information led up to 140% higher mean predictive abilities in comparison to using diploidized markers. Including dominance effects were highly advantageous when using diploidized markers, leading to mean predictive abilities which were up to 115% higher in comparison to only including additive effects. When the frequency of heterozygous genotypes in the population was low, such as in the sugarcane and sweet potato datasets, there was little advantage in including allele dosage information in the models. Overall, we show that including allele dosage can improve genomic selection in highly polyploid species under higher frequency of different heterozygous genotypic classes and high dominance degree levels.

FMixFN: A Fast Big Data-Oriented Genomic Selection Model Based on an Iterative Conditional Expectation algorithm.

Genomic selection is an approach to select elite breeding stock based on the use of dense genetic markers and that has led to the development of various models to derive a predictive equation. However, the current genomic selection software faces several issues such as low prediction accuracy, low computational efficiency, or an inability to handle large-scale sample data. We report the development of a genomic prediction model named FMixFN with four zero-mean normal distributions as the prior distributions to optimize the predictive ability and computing efficiency. The variance of the prior distributions in our model is precisely determined based on an F2 population, and genomic estimated breeding values (GEBV) can be obtained accurately and quickly in combination with an iterative conditional expectation algorithm. We demonstrated that FMixFN improves computational efficiency and predictive ability compared to other methods, such as GBLUP, SSgblup, MIX, BayesR, BayesA, and BayesB. Most importantly, FMixFN may handle large-scale sample data, and thus should be able to meet the needs of large breeding companies or combined breeding schedules. Our study developed a Bayes genomic selection model called FMixFN, which combines stable predictive ability and high computational efficiency, and is a big data-oriented genomic selection model that has potential in the future. The FMixFN method can be freely accessed at https://zenodo.org/record/5560913 (DOI: 10.5281/zenodo.5560913).

Accelerating wheat breeding for end-use quality through association mapping and multivariate genomic prediction.

In hard-winter wheat (Triticum aestivum L.) breeding, the evaluation of end-use quality is expensive and time-consuming, being relegated to the final stages of the breeding program after selection for many traits including disease resistance, agronomic performance, and grain yield. In this study, our objectives were to identify genetic variants underlying baking quality traits through genome-wide association study (GWAS) and develop improved genomic selection (GS) models for the quality traits in hard-winter wheat. Advanced breeding lines (n = 462) from 2015-2017 were genotyped using genotyping-by-sequencing (GBS) and evaluated for baking quality. Significant associations were detected for mixograph mixing time and bake mixing time, most of which were within or in tight linkage to glutenin and gliadin loci and could be suitable for marker-assisted breeding. Candidate genes for newly associated loci are phosphate-dependent decarboxylase and lipid transfer protein genes, which are believed to affect nitrogen metabolism and dough development, respectively. The use of GS can both shorten the breeding cycle time and significantly increase the number of lines that could be selected for quality traits, thus we evaluated various GS models for end-use quality traits. As a baseline, univariate GS models had 0.25-0.55 prediction accuracy in cross-validation and from 0 to 0.41 in forward prediction. By including secondary traits as additional predictor variables (univariate GS with covariates) or correlated response variables (multivariate GS), the prediction accuracies were increased relative to the univariate model using only genomic information. The improved genomic prediction models have great potential to further accelerate wheat breeding for end-use quality.

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 90K SNP array or DArTseq. We predicted the management systems (CV0M) or environments (CV0), a subset of lines that have been evaluated in either management (CV2M) or some environments (CV2), and the performance of newly developed lines in either management (CV1M) 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.

Comparison of genomic selection models for exploring predictive ability of complex traits in breeding programs.

Traits with a complex unknown genetic architecture are common in breeding programs. However, they pose a challenge for selection due to a combination of complex environmental and pleiotropic effects that impede the ability to create mapping populations to characterize the trait's genetic basis. One such trait, seedling emergence of wheat (Triticum aestivum L.) from deep planting, presents a unique opportunity to explore the best method to use and implement genetic selection (GS) models to predict a complex trait. Seventeen GS models were compared using two training populations, consisting of 473 genotypes from a diverse association mapping panel phenotyped from 2015 to 2019 and the other training population consisting of 643 breeding lines phenotyped in 2015 and 2020 in Lind, WA, with 40,368 markers. There were only a few significant differences between GS models, with support vector machines reaching the highest accuracy of 0.56 in a single breeding line trial using cross-validations. However, the consistent moderate accuracy of the parametric models indicates little advantage of using nonparametric models within individual years, but the nonparametric models show a slight increase in accuracy when combing years for complex traits. There was an increase in accuracy using cross-validations from 0.40 to 0.41 using diversity panels lines to breeding lines. Overall, our study showed that breeders can accurately predict and implement GS for a complex trait by using nonparametric machine learning models within their own breeding programs with increased accuracy as they combine training populations over the years.

Quantitative trait locus mapping and genomic selection of tobacco (Nicotiana tabacum L.) based on high-density genetic map

Tobacco (Nicotiana tabacum L.) is an economic crop and a model organism for studies of plant biology and genetics. As an allotetraploid plant generated from interspecific hybridization, tobacco has a massive genome (4.5 Gb). Recently, a genetic map with 45,081 single nucleotide polymorphism (SNP) markers was constructed using whole-genome sequencing data for a tobacco population including 274 individuals. This provides a basis for quantitative trait locus (QTL) mapping and genomic selection, which have been widely applied to other crops but have not been feasible in tobacco. Based on this high-density genetic map, we identified QTLs associated with important agronomic traits, chemical compounds in dry leaves, and hazardous substances in processed cigarettes. The LOD values for major QTLs were highest for agronomic traits, followed by chemical compounds and hazardous substances. In addition to the identification of molecular markers, we evaluated genomic selection models and found that BayesB had the highest prediction accuracy for the recombinant inbred line population. Our results offer new insights into the genetic mechanism underlying important traits, such as agronomic traits and quality-related chemical compounds in tobacco, and will be able to support the application of molecular breeding to tobacco.

Predicting Fusarium Head Blight Resistance for Advanced Trials in a Soft Red Winter Wheat Breeding Program With Genomic Selection.

Many studies have evaluated the effectiveness of genomic selection (GS) using cross-validation within training populations; however, few have looked at its performance for forward prediction within a breeding program. The objectives for this study were to compare the performance of naïve GS (NGS) models without covariates and multi-trait GS (MTGS) models by predicting two years of F4:7 advanced breeding lines for three Fusarium head blight (FHB) resistance traits, deoxynivalenol (DON) accumulation, Fusarium damaged kernels (FDK), and severity (SEV) in soft red winter wheat and comparing predictions with phenotypic performance over two years of selection based on selection accuracy and response to selection. On average, for DON, the NGS model correctly selected 69.2% of elite genotypes, while the MTGS model correctly selected 70.1% of elite genotypes compared with 33.0% based on phenotypic selection from the advanced generation. During the 2018 breeding cycle, GS models had the greatest response to selection for DON, FDK, and SEV compared with phenotypic selection. The MTGS model performed better than NGS during the 2019 breeding cycle for all three traits, whereas NGS outperformed MTGS during the 2018 breeding cycle for all traits except for SEV. Overall, GS models were comparable, if not better than phenotypic selection for FHB resistance traits. This is particularly helpful when adverse environmental conditions prohibit accurate phenotyping. This study also shows that MTGS models can be effective for forward prediction when there are strong correlations between traits of interest and covariates in both training and validation populations.

Genome optimization via virtual simulation to accelerate maize hybrid breeding.

The employment of doubled-haploid (DH) technology in maize has vastly accelerated the efficiency of developing inbred lines. The selection of superior lines has to rely on genotypes with genomic selection (GS) model, rather than phenotypes due to the high expense of field phenotyping. In this work, we implemented 'genome optimization via virtual simulation (GOVS)' using the genotype and phenotype data of 1404 maize lines and their F1 progeny. GOVS simulates a virtual genome encompassing the most abundant 'optimal genotypes' or 'advantageous alleles' in a genetic pool. Such a virtually optimized genome, although can never be developed in reality, may help plot the optimal route to direct breeding decisions. GOVS assists in the selection of superior lines based on the genomic fragments that a line contributes to the simulated genome. The assumption is that the more fragments of optimal genotypes a line contributes to the assembly, the higher the likelihood of the line favored in the F1 phenotype, e.g. grain yield. Compared to traditional GS method, GOVS-assisted selection may avoid using an arbitrary threshold for the predicted F1 yield to assist selection. Additionally, the selected lines contributed complementary sets of advantageous alleles to the virtual genome. This feature facilitates plotting the optimal route for DH production, whereby the fewest lines and F1 combinations are needed to pyramid a maximum number of advantageous alleles in the new DH lines. In summary, incorporation of DH production, GS and genome optimization will ultimately improve genomically designed breeding in maize. Short abstract: Doubled-haploid (DH) technology has been widely applied in maize breeding industry, as it greatly shortens the period of developing homozygous inbred lines via bypassing several rounds of self-crossing. The current challenge is how to efficiently screen the large volume of inbred lines based on genotypes. We present the toolbox of genome optimization via virtual simulation (GOVS), which complements the traditional genomic selection model. GOVS simulates a virtual genome encompassing the most abundant 'optimal genotypes' in a breeding population, and then assists in selection of superior lines based on the genomic fragments that a line contributes to the simulated genome. Availability of GOVS (https://govs-pack.github.io/) to the public may ultimately facilitate genomically designed breeding in maize.

Increased Predictive Accuracy of Multi-Environment Genomic Prediction Model for Yield and Related Traits in Spring Wheat (Triticum aestivum L.).

Genomic selection (GS) has the potential to improve the selection gain for complex traits in crop breeding programs from resource-poor countries. The GS model performance in multi-environment (ME) trials was assessed for 141 advanced breeding lines under four field environments via cross-predictions. We compared prediction accuracy (PA) of two GS models with or without accounting for the environmental variation on four quantitative traits of significant importance, i.e., grain yield (GRYLD), thousand-grain weight, days to heading, and days to maturity, under North and Central Indian conditions. For each trait, we generated PA using the following two different ME cross-validation (CV) schemes representing actual breeding scenarios: (1) predicting untested lines in tested environments through the ME model (ME_CV1) and (2) predicting tested lines in untested environments through the ME model (ME_CV2). The ME predictions were compared with the baseline single-environment (SE) GS model (SE_CV1) representing a breeding scenario, where relationships and interactions are not leveraged across environments. Our results suggested that the ME models provide a clear advantage over SE models in terms of robust trait predictions. Both ME models provided 2–3 times higher prediction accuracies for all four traits across the four tested environments, highlighting the importance of accounting environmental variance in GS models. While the improvement in PA from SE to ME models was significant, the CV1 and CV2 schemes did not show any clear differences within ME, indicating the ME model was able to predict the untested environments and lines equally well. Overall, our results provide an important insight into the impact of environmental variation on GS in smaller breeding programs where these programs can potentially increase the rate of genetic gain by leveraging the ME wheat breeding trials.

Tomato genomic prediction for good performance under high-temperature and identification of loci involved in thermotolerance response

Many studies showed that few degrees above tomato optimum growth temperature threshold can lead to serious loss in production. Therefore, the development of innovative strategies to obtain tomato cultivars with improved yield under high temperature conditions is a main goal both for basic genetic studies and breeding activities. In this paper, a F4 segregating population was phenotypically evaluated for quantitative and qualitative traits under heat stress conditions. Moreover, a genotyping by sequencing (GBS) approach has been employed for building up genomic selection (GS) models both for yield and soluble solid content (SCC). Several parameters, including training population size, composition and marker quality were tested to predict genotype performance under heat stress conditions. A good prediction accuracy for the two analyzed traits (0.729 for yield production and 0.715 for SCC) was obtained. The predicted models improved the genetic gain of selection in the next breeding cycles, suggesting that GS approach is a promising strategy to accelerate breeding for heat tolerance in tomato. Finally, the annotation of SNPs located in gene body regions combined with QTL analysis allowed the identification of five candidates putatively involved in high temperatures response, and the building up of a GS model based on calibrated panel of SNP markers.

Accuracy of genomic prediction for seed oil concentration in high‐oleic soybean populations using a low‐density marker panel

AbstractInsoybean [Glycine max (L.) Merr.], seed oil concentration is a complex quantitative trait, and genomic selection (GS) has been shown to be a valuable tool for performing selection on such traits. The objectives of this study were to evaluate multiple GS models for seed oil concentration using a low‐density marker panel in four biparental soybean populations and to assess predictive ability of the models using six unique training populations (TPs). Individuals were grown as BC1F4:F5 progeny rows in 2014. Genomic estimated breeding values (GEBVs) were calculated for each genotype within a population using genomic best linear unbiased predictor (GBLUP), BayesA, and BayesB models in a biparental specific context. In 2015, 60 individuals from each population were randomly selected and grown at six locations with two replications each to generate a “true” phenotypic value for each genotype. Prediction accuracies for each estimation set were generated by correlating the GEBVs with the “true” phenotypic value. Across all populations, the GS prediction accuracy was greatest using GBLUP; however, no GS prediction model showed a significant advantage in accuracy over the phenotypic values. Generally, TPs consisting of more individuals had higher prediction accuracies; however, variations were observed across populations and models. The results show that GS in a biparental context and with low marker densities can be a valuable tool for breeders focused on making gains for oil concentration; however, consideration must be given as to how to apply these methods to each situation.

Experimental evaluation of genomic selection prediction for rust resistance in sugarcane

The total sugarcane (SaccharumL.) production has increased worldwide; however, the rate of growth is lower compared with other major crops, mainly due to a plateauing of genetic gain. Genomic selection (GS) has proven to substantially increase the rate of genetic gain in many crops. To investigate the utility of GS in future sugarcane breeding, a field trial was conducted using 432 sugarcane clones using an augmented design with two replications. Two major diseases in sugarcane, brown and orange rust (BR and OR), were screened artificially using whorl inoculation method in the field over two crop cycles. The genotypic data were generated through target enrichment sequencing technologies. After filtering, a set of 8,825 single nucleotide polymorphic markers were used to assess the prediction accuracy of multiple GS models. Using fivefold cross-validation, we observed GS prediction accuracies for BR and OR that ranged from 0.28 to 0.43 and 0.13 to 0.29, respectively, across two crop cycles and combined cycles. The prediction ability further improved by including a known major gene for resistance to BR as a fixed effect in the GS model. It also substantially reduced the minimum number of markers and training population size required for GS. The nonparametric GS models outperformed the parametric GS suggesting that nonadditive genetic effects could contribute genomic sources underlying BR and OR. This study demonstrated that GS could potentially predict the genomic estimated breeding value for selecting the desired germplasm for sugarcane breeding for disease resistance.

Genetic and phenotypic characterization of rice grain quality traits to define research strategies for improving rice milling, appearance, and cooking qualities in Latin America and the Caribbean.

Rice (Oryza sativa L.)grain quality is a set of complex interrelated traits that include grain milling, appearance, cooking, and edible properties. As consumer preferences in Latin America and the Caribbean evolve, determining what traits best capture regional grain quality preferences is fundamental for breeding and cultivar release. In this study, a genome-wide association study (GWAS), marker-assisted selection (MAS), and genomic selection (GS) were evaluated to help guide the development of new breeding strategies for rice grain quality improvement. For this purpose, 284 rice lines representing over 20 yr of breeding in Latin America and the Caribbean were genotyped and phenotyped for 10 different traits including grain milling, appearance, cooking, and edible quality traits. Genetic correlations among the 10 traits ranged from -0.83 to 0.85. A GWAS identified 19 significant marker/trait combinations associated with eight grain quality traits. Four functional markers, three located in the Waxy and one in the starch synthase IIa genes, were significantly associated with six grain-quality traits. These markers individually explained 51-75% of the phenotypic variance depending on the trait, clearly indicating their potential utility for MAS. Cross-validation studies to evaluate predictive abilities of four different GS models for each of the 10 quality traits were conducted and predictive abilities ranged from 0.3 to 0.72. Overall, the machine learning model random forest had the highest predictive abilities and was especially effective for traits where large effect quantitative trait loci were identified. This study provides the foundation for deploying effective molecular breeding strategies for grain quality in Latin American rice breeding programs.