Next Article in Journal
Age Related Differences in Monocyte Subsets and Cytokine Pattern during Acute COVID-19—A Prospective Observational Longitudinal Study
Next Article in Special Issue
Analyses of Lysin-motif Receptor-like Kinase (LysM-RLK) Gene Family in Allotetraploid Brassica napus L. and Its Progenitor Species: An In Silico Study
Previous Article in Journal
Growth Hormone and IGF1 Actions in Kidney Development and Function
Previous Article in Special Issue
Metabolite Diversity and Metabolic Genome-Wide Marker Association Studies (Mgwas) for Health Benefiting Nutritional Traits in Pearl Millet Grains
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Strategies to Increase Prediction Accuracy in Genomic Selection of Complex Traits in Alfalfa (Medicago sativa L.)

1
United States Department of Agriculture-Agricultural Research Service, Plant Germplasm Introduction and Testing Research, Prosser, WA 99350, USA
2
Department of Plant and Environmental Sciences, New Mexico State University, Las Cruces, NM 88003, USA
*
Author to whom correspondence should be addressed.
Cells 2021, 10(12), 3372; https://doi.org/10.3390/cells10123372
Submission received: 19 October 2021 / Revised: 19 November 2021 / Accepted: 24 November 2021 / Published: 30 November 2021
(This article belongs to the Special Issue Omics in Plant Genetics and Breeding)

Abstract

:
Agronomic traits such as biomass yield and abiotic stress tolerance are genetically complex and challenging to improve through conventional breeding approaches. Genomic selection (GS) is an alternative approach in which genome-wide markers are used to determine the genomic estimated breeding value (GEBV) of individuals in a population. In alfalfa (Medicago sativa L.), previous results indicated that low to moderate prediction accuracy values (<70%) were obtained in complex traits, such as yield and abiotic stress resistance. There is a need to increase the prediction value in order to employ GS in breeding programs. In this paper we reviewed different statistic models and their applications in polyploid crops, such as alfalfa and potato. Specifically, we used empirical data affiliated with alfalfa yield under salt stress to investigate approaches that use DNA marker importance values derived from machine learning models, and genome-wide association studies (GWAS) of marker-trait association scores based on different GWASpoly models, in weighted GBLUP analyses. This approach increased prediction accuracies from 50% to more than 80% for alfalfa yield under salt stress. Finally, we expended the weighted GBLUP approach to potato and analyzed 13 phenotypic traits and obtained similar results. This is the first report on alfalfa to use variable importance and GWAS-assisted approaches to increase the prediction accuracy of GS, thus helping to select superior alfalfa lines based on their GEBVs.

1. Introduction

Alfalfa (Medicago sativa L.) is an autotetraploid (2n = 4x = 32) perennial forage crop with a genome size of 800–1000 Mb [1]. However, alfalfa breeding is complicated by its high heterozygosity, polysomic inheritance, and out-crossing nature, which hinder the creation of inbred lines. Alfalfa breeding goals target improvement of forage yield, quality, and tolerance to biotic and abiotic stresses. This process requires the selection of perennial plants that can maintain biomass productivity and quality over several years. Therefore, traits must be evaluated over multiple harvests each year for several years. Consequently, genetic gain is slower compared to annual crops. In addition, alfalfa breeding programs have largely focused on recurrent phenotypic selection (PS) in field environments to improve quantitative traits of interest. However, this approach is constrained by breeding population size, genotype × environment interactions, or low heritability of the trait, thus hindering the development of superior varieties.
One promising alternative to recurrent PS is indirect selection based on the use of molecular markers generated, for example, via genotyping by sequencing (GBS) [2]. Markers closely linked to quantitative trait loci (QTL) can then be used for marker-assisted selection (MAS) in breeding programs. Initially, QTLs are detected through genetic mapping or genome-wide association studies (GWAS), where marker-trait associations that exceed specific thresholds are declared statistically significant (Figure 1a). However, MAS is primarily effective for traits controlled by relatively few genes with large effects. For complex traits (e.g., stress tolerance or yield) in elite populations, it can be difficult to clearly identify QTL with major effect because the trait is often controlled by multiple loci possessing small effects.
Soil salinity is one of the major environmental factors limiting the productivity of alfalfa. Genomic tools have been applied to identify important genetic factors that influence salt tolerance in alfalfa. Genetic loci associated with salt stress tolerance have been identified from different studies. Yu et al. (2016) identified 23 markers and 14 functional genes associated with germination under salt stress [3]. Liu et al. (2017) identified 42 markers associated with forage yield, plant height, leaf chlorophyll content, and stomatal conductance under salt stress [4]. Liu et al. (2019) identified 49 markers associated with forage yield, plant height, leaf relative water content, and stomatal conductance under salt stress [5]. Most recently, Medina et al. (2020) identified 27 markers associated with biomass yield and plant vigor under salt stress [6]. Those results highlighted the genetic complexity of salt stress response in alfalfa.
Conventional breeding strategy for improving salt tolerance of alfalfa is time consuming and less effective. Genomic selection (GS) offers the potential to shorten alfalfa breeding and selection cycles. GS is a promising alternative to determine the genetic potential or breeding value of an individual based on whole-genome markers (Figure 1a). This method follows the infinitesimal model, which assumes that a quantitative trait is determined by an infinite number of unlinked and non-epistatic loci, each one with a very small effect that satisfies normality and linearity [7]. This technique uses both parametric and non-parametric statistical models to determine associations of phenotypic trait values with genome-wide molecular markers. This information is subsequently used to predict future breeding values (i.e., genomic-estimated breeding values, GEBVs) for each individual in a population based on their genome-wide marker profile/genotype [8]. Hence, rapid marker-based selection cycles can replace some time-intensive phenotypic selection cycles to accelerate genetic gain. In this paper, we review different GS models and their application to polyploid crops. We also demonstrate the implementation of GS models on a real dataset of alfalfa and potato to identify improved approaches to implement GS in different breeding programs.

2. Statistical Methods in GS

There is a wide repertoire of parametric and non-parametric models to obtain GEBVs that differ in complexity, accuracy, and computational requirements (Figure 1b). Some phenotypic traits are highly complex and more difficult to predict using their genetic information, therefore, accuracy in GS modeling is a cornerstone. Model accuracy metric is calculated as the Pearson’s correlation coefficient ( r G E B V : y ) between GEBVs in a training population and observed phenotypes in a testing population. Determining a GEBV can be solved as a regression:
y = X β + e  
where y is a vector ( n × 1 ) of phenotypic outcomes in n observations, X is a matrix ( n × p ) with p number of markers or predictors in n observations, β is the vector ( p × 1 ) of marker effects and e is a vector of residual effects. In GS, however, molecular markers or predictors (p) are greater than observations (n), generating a large p small n problem ( p n ). Therefore, estimation of marker effects via multiple regression by ordinary least squares is not possible. To resolve this issue, multiple methods have been developed to handle the high dimensionality of the genomic data. Shrinkage models, such as best linear unbiased prediction using ridge-regression (RRBLUP) [9] or genomic best linear unbiased prediction (GBLUP) [10], are the most popular models used in GS. Both of these models assume that the effect of all single nucleotide polymorphism (SNP) markers is normally distributed with equal variance [9]. Bayesian and least absolute shrinkage and selection operator (LASSO) models assume that some SNPs have large or moderate effects, and others have small or null effects [11]. Finally, machine learning (ML) models like random forest (RF), support vector machine (SVM) and deep learning (DL) algorithms do not assume linearity in the model. ML models can use nonlinear kernels to capture complex SNP–SNP interactions and nonlinear relationships.

2.1. Ridge-Regression Best Linear Unbiased Prediction (RRBLUP)

The RRBLUP is a shrinkage method to obtain GEBVs by incorporating genomic information into BLUP using ridge regression (RR). This model has been implemented in the R package RRBLUP [9]. Prediction equations used by RRBLUP assume a priori that all loci explain equal amounts of the genetic variation. The core of the RRBLUP package is the function mixed.solve, which solves any mixed model of the form:
y = X β + Z u + e  
where y is a vector ( n × 1 ) of phenotypic outcomes in n observations, X is a matrix ( n × p ) with p number of markers or predictors in n observations, β is a vector ( p × 1 ) of fixed effects, u ~ N ( 0 , K σ u 2 ) is a vector ( n × 1 ) of random effects distributed normally with mean zero and variance σ u 2 and K is a positive semidefinite matrix, Z is a design matrix ( n × p ) for the random effects, and e ~ N ( 0 , I σ e 2 ) is a vector ( n × 1 ) of residual effects distributed normally with mean zero and variance σ e 2 and I is the identity matrix.

2.2. Genomic Best Linear Unbiased Prediction (GBLUP)

GBLUP measures the relationship between individuals with the aid of marker data. The difference with RRBLUP is the use of a marker-based relationship matrix named the genomic relationship matrix (GRM) or G matrix [12]. The G matrix defines the covariance between known relatives in a population, based on DNA marker information. The mixed model for GBLUP analysis uses the following formula:
y = 1 μ + X β + Z g + e  
where y , X , β and e were defined in Equation (2), μ is the overall mean, g ~ N ( 0 , G σ g 2 ) is a vector ( n × 1 ) of random effects distributed normally with mean zero and variance σ g 2 and G is the G matrix, which can be obtained according to the approach of VanRaden [13]:
G = Z Z 2 Σ p i ( 1 p i )  
where Z is an identity matrix for the markers and p i is the observed minimum allele frequency (MAF) of all individuals genotyped. However, Yang et al. (2010) combined the information on all SNPs (i) coded as 0 = AA, 1 = BB, 2 = AB according to alternative allele dosage to calculate the relationship between individuals j and k into a GRM ( G i j k ) using a weighting scheme based on allele frequencies:
G j k = 1 N i G i j k = { 1 N i ( w i j 2 p i ) ( w i k 2 p i ) 2 p i ( 1 p i ) , j k   1 + 1 N i w i j 2 ( 1 + 2 p i ) w i j + 2 p i 2 2 p i ( 1 p i ) , j = k  
where G j k is the G matrix averaged over all SNP positions in the genome, N is the number of markers, w i j is the element of W pertaining to marker i and individual j, and w i k is the element of W pertaining to marker i and individual k. The G i j k or GRM produces the off-diagonal ( j k ) and diagonal ( j = k ) elements [14]. Based on this approach, Slater et al. (2016) proposed a full autotetraploid model to obtain the G matrix:
G j k = { 1 + 1 M i ( w i j p i ) ( w i k p i ) p i ( 1 p i ) , j k   1 + 1 M i w i j 2 2 p i w i j + p i 2 p i ( 1 p i ) , j = k  
where M is the number of markers × 5 and p i is the frequency of each genotype. The genomic relationship matrices described are based on identity-by-state and simply measure the similarity of alleles between individuals [15].
Analysis results from RRBLUP and GBLUP can be similar; however, GBLUP is more computationally efficient than RRBLUP. GBLUP requires a G matrix of dimensions n × n (where n is the number of individuals in the population), whereas RRBLUP requires a genotypic matrix n × m (where m is the number of markers) with high dimensionality. In summary, GBLUP does not provide marker effects but is more time/memory efficient than RRBLUP.

Weighted Genomic Best Linear Unbiased Prediction (WGBLUP)

The GBLUP method usually assumes that all SNPs explain the same fraction of genetic variance. However, traits are affected by different genetic architectures associated with SNPs that possess varying effects (e.g., major SNPs). To account for varying effects of different SNP alleles, the weighted GBLUP (WGBLUP) method was developed to incorporate unequal weights for all SNPs [16]. The G * matrix is constructed as follows:
G * = Z D Z 2 Σ p i ( 1 p i )  
where the asterisk symbol ( G * ) is used to differentiate the weighted G matrix from the regular G matrix, Z is an identity matrix for the markers, D is a diagonal matrix, where each element of the diagonal corresponds to SNP weights, and p i is the observed minimum allele frequency (MAF) of all genotyped individuals. To obtain the D matrix, each element of this matrix is defined as:
D = ( w 1 0 0 w n )
where w is weight based on SNP effects from different methods. SNP weights can be obtained from Bayesian Regressions with good results [17]. Other methods to determine weights for SNPs include prioritization based on Wright’s fixation scores ( F s t ) [18]. The F s t score measures the level of genetic differentiation between populations based on a change in allele frequencies. When F s t scores were used to compute relative weights, prediction accuracy increased up to 5%. Ren et al. (2021) developed several methods to optimize WGBLUP by the generation of different weighted G matrices. They noted that the choice of an optimal GBLUP matrix will depend on the number of loci controlling the trait. Results indicated that estimated marker-variance-weighted (EVW)-GBLUP was superior for traits controlled by loci of a large effect, and absolute value of the estimated marker-effect-weighted (AEW)-GBLUP was better for traits controlled by loci with moderate effect [19].

2.3. Bayesian Models

Bayesian models applied to GS do not assume a normal distribution of marker effects. Instead, they assume that few markers will have large effects on the trait, allowing markers to have different effects and variances. Bayesian models impose stronger shrinkage towards zero on small SNP effects and less shrinkage on relatively large SNP effects. The BGLR R package implements a large collection of Bayesian models [20]. The Bayesian models for continuous variables are represented by the equation:
y i = 1 μ + j = 1 m x i j β j + e i
where y i is the vector of adjusted phenotypic observations { y 1 , , y n } , μ is the overall mean for the trait, m is the number of markers or SNPs, x i j is the i t h genotype for j t h SNP, β j is a vector for the effect of the j t h SNP, and e i is a vector of residual effects with assumed normal distribution e ~ N ( 0 , I σ e 2 ) , where σ e 2 is the residual variance and I is the identity matrix. Bayes A, B, Cπ, and Bayesian LASSO (BL) are the most common models used. All models assume different prior distributions for SNP effects (Table 1).

2.4. Machine Learning Models

Machine learning is a field that involves the application of computer algorithms and statistical models to interpret and predict large datasets. Algorithmic modeling is a rapidly developing discipline with strong potential to provide accurate and informative analyses or predictions using large and complex data sets [24]. These models are widely used to solve problems across different disciplines, such as medicine, genomics, natural language processing, and stock market forecasting. Compared with classical statistical models, ML models have fewer assumptions about normality and distribution of data. One important remark is that ML models are being developed much faster than their interpretability, developing a new field to be explored. The most common problem with ML algorithms is data overfitting, which results in models that poorly predict the behavior of future data. To avoid this problem, it is necessary to use a robust validation method, such as cross-validation, which provides an indication of performance on new data. SVM and RF are the most common ML models for classification and regression in GS [25,26].

2.4.1. Support Vector Machine (SVM)

SVM is a machine learning algorithm used in classification or regression problems [27]. The objective of SVM is to find the best hyperplane with the maximal margin in an n-dimensional space (genotypic matrix) with respect to a given collection of data (phenotypic values) and predict the correct classification/regression of unseen examples. Support vector regression (SVR) is an application of the SVM. In SVR, each n -dimensional input vector ( x i ) of p SNP markers is associated with a y i as response variable (e.g., yield), where x i p and y i . Then, linear regression f ( x ) is performed using the following equation [28]:
f ( x ) = w x + b
where, w is a vector of unknown weights (i.e., regression coefficients) and b is the bias. The training data is used to learn w . The coefficients w and b are estimated by minimizing the following regularized loss function R ( C ) [10]:
R ( C ) = 1 2 w 2 + C i = 1 n L ε ( y i f ( x i ) )
where w 2 = w w, represents model complexity, and C is a positive cost parameter specified by the user. C determines the trade-off between model complexity and training error, y i f ( x i ) is the error associated with i t h training data point and L ε is the empirical error measured by ε -intensive loss-function:
L ε ( y i f ( x i ) ) = { 0   i f   | y i f ( x i ) | < ε | y i f ( x i ) | ε   otherwise
where the loss function is zero (“insensitive”) for any absolute error smaller than a predefined value ε . For an error value larger than ε , the loss function is the difference between the absolute error and ε . ε -SVR solutions is sparseness, with a fraction of errors equal to zero and thereby vanishing in the final model f ( x ) and only absolute errors > 0 are relevant and used as “support vectors”. f ( x ) can be assumed as linear or non-linear. For non-linear functions, the data can be mapped into a higher dimensionality space using a kernel space. In nonlinear SVR modeling, different kernels can be used to increase the predictive power of the model. The kernel function provides a solution to the classification/regression dataset by adding an additional dimension to the data. Different kernel functions can be selected to transform input data to feature space. Commonly used kernels in SVM include linear, polynomial, radial basis function (RBF), and sigmoidal kernels (Table 2). In high-dimensional data (i.e., microarrays or GS), the RBF kernel is preferred [28].

2.4.2. Random Forest (RF)

The RF method is a machine learning model for classification and regression problems based on the identification of an objective function and its optimization [29]. The objective function measures the distance between the RF output and desired scores to modify internal parameters to reduce this error. Random forests consist of numerous independent decision trees that are independently trained using a random subset of data. The final prediction is calculated as the average values over all the trees. The RF model attempts to reduce the computational cost to train the model, capture complex interactions and reduce the over-fitting risk in the data [26]. In each decision tree, multiple binary filters are applied to create a bifurcation generating branches and a treelike structure. Every point where the samples are filtered is called a decision node. Optimization in RF consists of determining the best way to split samples at decision nodes based on the predictors. RF can optimize four different hyperparameters to increase the predictive power of the model: total number of observations (N), total number of predictors (M), subset of predictors chosen for determining a decision tree (mtry), and total number of decision trees to generate the RF (ntree). Subsequently, RF creates a series of filters based on the predictor variables. Gini impurity score (as described below) and mean squared error are used to select the best variables in decision nodes.
The RF approach can provide accurate predictions with complex genomic datasets. A very useful feature in the RF model is a function-designated, variable importance metrics, which ranks each SNP’s impact according to the trait. Two approaches to computing variable importance include mean prediction accuracy decrease when a variable/marker is removed, and the mean decrease in impurity (or Gini importance). Gini impurity measures how well a potential split is able to separate the samples of two classes at a particular node. Important limitations when using RF algorithms for GS involve slower processing time when a large number of trees are chosen for the model, or the number of SNP markers is too high. For RF analysis, we observed that a processing limit was reached when the genotypic matrix was composed of more than 10,000 markers.

2.4.3. Deep Learning (DL)

Deep learning is a subfield of machine learning with great success in natural language processing, image recognition, or virtual assistance [30]. The DL architecture uses several layers of nonlinear processing units called hidden layers. Hidden layers allow the network to capture higher-order interactions from the data. This method uses artificial neural network architecture where the perceptron is the fundamental unit for comparison as a neuron in a biological neural network. The implementation of DL in genomic selection is recent, and some studies have reported a modest increase in prediction accuracy in comparison with parametric and non-parametric models [31,32,33,34]. In theory, deep learning could perform better for traits with large epistatic effects and low narrow-sense heritability, a concept which is reinforced by the high predictive ability of mixed models as prediction machinery [35]. Numerically encoded SNPs are the inputs to the first layer to produce a centered vector of phenotypes. The most common DL models used in GS are multilayer perceptron neural network (MLP), convolutional neural network (CNN), and recurrent neural networks (RNN).
The formal description of MLP is a feed-forward DL model composed of multiple perceptrons ordered in hidden layers in a directed graph. In MLP, each layer is fully connected with the next one by nonlinear activation functions, such as rectified linear unit activation function (ReLU), to minimize the mean square error. MLP is flexible because no assumption is made about the joint distribution of inputs and outputs. CNN is a special case of a neural network that uses convolution instead of full matrix multiplication in the hidden layers. The convolution is a function that can be defined as an “integral transform” to reduce the number of hyperparameters to be estimated. CNN was proposed to accommodate situations where input variables are distributed along a space pattern resembling an SNP matrix. CNN seems to perform best in GS because it can detect patterns in the genotypic matrix, discovering correlations between adjacent SNPs. In addition, CNN appears to perform better when epistatic components are important and the narrow-sense heritability is low [34]. It is important to note that DL depends on an adequate hyperparameter choice and high-performance computing with graphics processing units (GPUs) architecture, which can be challenging to implement in small breeding programs.

2.5. Other Models

Klápště et al. (2020) presented a strategy to generate a marker-based relationship matrix that prioritized markers using Partial Least Squares (PLS). This approach downweighs noisy predictors, but does not remove them from the model. The advantage of PLS is that it deals with multicollinearity and can handle several response variables at a time. The authors used PLS-Canonical Analysis (PLS-CA) for constructing marker-based relationship matrices with different numbers of markers. This strategy attempts to improve the accuracy of traits with low heritability by taking advantage of the genetic covariance common across all investigated traits. In order to perform the marker selection by PLS-CA, all individuals in the training population must be phenotyped for all traits that will be included in the analysis [36].
The incorporation of staking ensemble ML (SEML) in GS is a promising alternative to increase the predictive ability of ML models. The SEML method uses a meta-learning algorithm to determine how to best combine the predictions from two or more base ML models. Hence, SEML has the potential to generate predictions with better performance than any single model [37]. Liang et al. (2021) tested the prediction accuracy of the SEML approach using three ML models: SVM, kernel ridge regression, and elastic net in Loblolly pine, beef, and dairy cattle. On average, there was an increase of 7.70% in prediction accuracy in nine traits tested. However, Bayes B demonstrated higher prediction accuracies for some traits, including milk fat percentage or tree stiffness [38].

3. Genomic Selection in Polyploids

GS requires high-quality genome-wide markers to determine GEBVs. Two types of high-throughput genotyping methods can be employed: SNP arrays and GBS. There are SNPs arrays with different marker densities in potato [39] and wheat [40]. Alfalfa also has an array with 9277 SNPs [41]. However, its use has not been widely adopted, and GBS is currently the best option to obtain genome-wide markers. During the genotyping process by GBS, different types of markers, such as single nucleotide polymorphisms (SNPs), insertions/deletions (indels), or short tandem repeats (STRs) can be obtained. Genome-wide markers can then be arrayed in a genotypic matrix of m samples and n markers. The genotypic matrix can be filtered to retain only biallelic SNPs, which are the most abundant and stable markers for identifying QTLs associated with traits of interest [42].
Allele dosage counts alternative allele frequency for each biallelic SNP. In diploid species the genotypic matrix is coded as { 0 ,   1 ,   2 } , reflecting if a given marker is present in the homozygous reference (AA), heterozygous (AB), or homozygous alternate (BB) allelic state. For biallelic SNPs in polyploid species with ploidy N , the biallelic dosage is N + 1 and the genotypic matrix is coded as { 0 , , N } . Genotype calling in autotetraploids requires bioinformatics tools to distinguish among five possible genotypes (AAAA, AAAB, AABB, ABBB, BBBB) with biallelic SNPs coded as {0, 1, 2, 3, 4}. There are several R packages, such as polyRAD [43], superMASSA [44], FitTetra 2.0 [45] or Updog [46], with which to obtain allele dosage in numeric format from variant call format file (vcf) format. Some of these R packages, such as Updog, require users to specify genotype priors [46] to accurately calculate the allele dosage and distinguish between all possible genotypes. However, the most common option is to use high depth sequence reads (e.g., ~60×) which leads to 98.4% accuracy in genotypic calls [47] The effects of marker allele dosage on phenotype for genomic selection have been reported previously. Slater et al. (2016) described three different models for GS in autopolyploids: additive autotetraploid, pseudodiploid, and full autotetraploid. In the additive autotetraploid model, the allele dosage has an additive effect, and 0, 1, 2, 3, 4 corresponds to AAAA, AAAB, AABB, ABBB, BBBB, respectively. In the pseudodiploid model, all heterozygous genotypes (AAAB, AABB, ABBB) have the same effect of 1 on the genotypic variation, while the two homozygotes AAAA and BBBB have an effect of 0 and 2, respectively. Finally, the full autotetraploid model assumes that each genotype has its own effect with five possible effects per marker, assuming that the markers are fitted as random effects [15]. In addition, Rosyara et al. (2016) developed GWASpoly, a software for genome-wide association studies in autopolyploids. GWASpoly has different assumptions over allele dosages and conducts the hypothesis tests for each marker using six models (general, diploidized general, diploidized additive, additive, simplex dominant, and duplex dominant models) (Table 3).
Amadeu et al. (2019) evaluated the inclusion of dominance effects for genomic prediction in autotetraploid crops. They reported that a full autotetraploid model, including additive and dominance effects jointly modeled into a unique general effect, increased the total genetic variance explained [48]. In potato, different covariance genomic marker- and pedigree-based matrices, designated G and A, respectively, were tested to identify additive and nonadditive genetic effects and to improve the accuracy in GS. A matrix (also known as numerator relationship matrix) was calculated from a 13-generation pedigree. The A matrix was defined as a matrix containing kinship coefficients among all individuals in the population, multiplied by four. They reported that the G matrix was superior to the A matrix and adding allele dosage information increased the prediction accuracy. Finally, the use of a pseudodiploid matrix reduced the prediction accuracy by 0.13, on average [49].
In the autotetraploid forage grass Panicum maximum, de C. Lara et al. (2019) compared the predictive ability of six GS models in six traits using tetraploid and pseudodiploid allele dosages and a minimum depth of 25 reads. Additionally, multiple harvests were modeled with a variance-covariance matrix for genotypes nested across harvests, treating the genotypes as a random factor. The incorporation of correlations among harvests provided a better fit for the traits analyzed. Including the tetraploid dosage also produced higher predictive accuracy compared with pseudodiploid dosages. In autotetraploids with highly mixed ploidy, such as sugarcane and sweet potato, the incorporation of allele dosage information increased model predictive abilities up to 140% in comparison to using diploidized markers [50]. The accuracy of different models showed few changes to ploidy or allele dosage information in sugarcane on sweet potato. In sugarcane, the Brix trait possessed the highest mean predictive accuracy (0.24) using the GBLUP model that included allele dosage. In sweet potato, prediction accuracies were moderate to high. For example, color saturation had the highest mean predictive ability (0.75) using a G model with allele dosage information.
In blueberry (Vaccinium spp.), several approaches have been evaluated to improve the genomic selection process, because the conventional breeding pipeline takes up to 12 years [35,51,52]. Implementation of GS in the early stages of the breeding program could shorten the cycle time to nine years and increase the expected genetic gain by 86%. De Bem Oliveira et al. (2019) compared diploid, tetraploid, and continuous allele dosages at the individual plant level for the application of genomic selection in potato and blueberry. In general, there was no difference among the models tested, but continuous genotypes resulted in a better predictive ability for some traits, such as fruit firmness, fruit scar, and fruit diameter. Furthermore, the use of a marker-based relationship matrix generated better predictions than a pedigree-based relationship matrix (A matrix). Ferrão et al. (2021) reported similar prediction accuracies of GBLUP for four traits using two genotype calling approaches (dosage and ratio) and two read-depth scenarios (6× and 60×). They also observed that combining allele dosage for low to mid sequencing depths (6×–12×) produced similar accuracies to that obtained by high read-depth (60×). The use of mid sequencing depths will allow modifying economic resource allocation to increase the number of individuals genotyped.
Enormous progress has been made during the last few years in the application of GS approaches to polyploids. In alfalfa, GS has been tested in different traits using parametric and non-parametric models. Table 4 summarizes the progress that has been made towards applying GS in multiple polyploids, including alfalfa. Although yield is the main trait in alfalfa breeding programs, other agronomic traits, such as forage quality and plant regrowth, have also been tested [53]. Furthermore, use of an allele dosage genotype matrix has been reported to improve prediction accuracies of forage quality [54] and yield under salt stress [6]. The current challenge is to implement GS in breeding programs and to evaluate increases in GS-derived genetic gain in comparison with PS-derived materials.

4. Case Study: Logan 2020 Population

In this review, we tested some models of GS using the dataset of alfalfa previously published [6]. Datasets were collected from a multi-parental population generated to select lines tolerant to salt stress. Forage yield under salt stress was measured over seven harvests in 265 individuals for two years. Each harvest was spatially corrected by a two-dimensional P-spline mixed model with the Mr.Bean web application [65] using the SpATS package [66]. Multiple best linear unbiased estimator (BLUE) values were adjusted in multi-environmental trials using Factor Analytic II covariance structure [67] with ASReml R software [68]. A genotypic matrix of 6796 high-quality GBS-derived SNP markers was obtained using NGSEP v4.0.0 software [69] with parameters previously reported in [6]. SNPs were coded from 0 to 4 according to allele types using the Sommer R package [70].
A GS approach using regression analysis between phenotypes ( y ) and a genotypic matrix was transformed by eight different models ( f ( X ) ). These included: RRBLUP, Bayes A, B, C, Bayesian Lasso (BL), GBLUP using two G matrices (VanRaden [VR] [Equation (4)] and full-autotetraploid [FA] [Equation (6)]), RF and SVM (Figure 2a). All models were compared for execution time and Pearson’s correlation using ten-fold cross-validation with the GROAN R package [71]. Execution time is an important factor to consider for GS modeling when computing power is limited. Consequently, system time (seconds) was measured for each model with cross-validation. The fastest models were GBLUP and RRBLUP with an average of 0.06 s, whereas SVM and RF required 10.57 and 12.99 s, respectively. More time was required for ML models when a grid search was used to estimate the best values for hyperparameters such as cost and sigma in SVM, or mtry in RF. However, time was reduced to 1.90 s with cross-validation in the SVM training model when hyperparameters were previously defined. Prediction accuracy was approximately 0.3 among the GS models RRBLUP, Bayes A, B, C, BL, RF, GBLUP-VR, and GBLUP-FA. The SVM model possessed the highest accuracy (0.46), in agreement with previous reports [6,59].
Variable importance values or SNP weights were obtained using SVM and RF models with the Caret R package [72], or by retrieving −log10 p-values resulting from six models of the GWASpoly R package (i.e., general, diploidized general, diploidized additive, additive, simplex dominant and duplex dominant models [Table 3]) [73]. SNP weights were used as input in a D diagonal matrix [Equation (8)] for the construction of a G * matrix [Equation (7)] in the WGBLUP model (Figure 2b,c). Pearson’s correlation among variable importance values of different models was measured to identify models with similar SNPs weights. Diploidized additive and diploidized general models had the highest Pearson’s correlation (0.87), followed by additive and diploidized additive models (0.74). Variable importance values derived from RF had low correlations across all models tested (Figure 2c). Prediction accuracies for GBLUP with two G matrix and 10 WGBLUP models were compared by measuring Pearson’s correlation 10 times with ten-fold cross-validation. The incorporation of variable importance values in WGBLUP increased prediction accuracies. Pearson correlations ranged from a low of 0.32 in GBLUP-VR (no variable importance values) to 0.63 in WGBLUP-SVM, with the highest prediction accuracy (0.83) achieved when −log10 p-values from the additive GWASpoly model were used as a weight vector (Figure 2d). Thus, incorporation of a diagonal matrix D with variable importance values to the G matrix increased GS predictive ability almost three times without increasing computational time. This is the first report using SNPs weights to increase the prediction accuracy of GS in alfalfa. Our results suggest that including SNP marker −log10 p-values derived from the additive GWASpoly model in a WGBLUP model may benefit prediction accuracy and selection for improvement of complex traits in alfalfa breeding programs.

5. Case Study: Potato Diversity Panel

To ensure that our approach is useful in other polyploid crops, we extended GS analysis in potato. RRBLUP, GBLUP and, WGBLUP models were evaluated in potato using supplemental files from [73]. Phenotypic data consisted of 13 agronomic traits evaluated in 187 lines of the SolCAP potato diversity panel were used. The genotypic matrix was generated by an Infinium SNP array containing 3521 markers with allele dosage. The genotypic matrix was transformed to digital format using the sommer R package [70]. The GBLUP and WGBLUP models were tested using the VanRaden G matrix [Equation (4)] and the WGBLUP equation [Equation (7)], respectively. Six different D matrices were generated according to the SNP −log10 p-values from six models of the GWASpoly (i.e., general, diploidized general, diploidized additive, additive, simplex dominant, and duplex dominant models [Table 3]). G matrices were constructed using the function Gmatrix from the AGHmatrix R package [74]. Prediction accuracies for RRBLUP, GBLUP, and WGBLUP models were compared by measuring Pearson’s correlation 10 times with 10-fold cross-validation using the GROAN R package (Table 5) [71].
Higher prediction accuracies in all 13 agronomic traits of potato were obtained using the WGBLUP model with SNP −log10 p-values derived from the additive GWASpoly model. Traits of glucose, tuber length, or tuber shape showed accuracies higher than 0.9. It is important to point out that traits of tuber length or tuber shape had high accuracies (0.82 and 0.78 respectively using RRBLUP and GBLUP models) and the use of the WGBLUP model increased the prediction accuracy up to 0.93. Total yield had low prediction accuracies with RRBLUP or GBLUP models (0.132 and 0.117, respectively), and the use of the WGBLUP model increased prediction accuracy by almost five times (Table 5). These results agree with our previous results in alfalfa (Figure 2).

6. Conclusions

Genomic selection is a breeding strategy that predicts the genomic estimated breeding value (GEBV) of individuals in a population using genomic-wide genetic markers. A significant advantage of GS is the ability to select superior individuals in a population at very early stages in a breeding cycle based on their genotype, versus conducting lengthy/expensive phenotyping trials prior to each selection cycle. For instance, in a simulated wheat breeding program, selection based on GEBVs for grain yield tripled the genetic gain compared with PS [75]. In alfalfa breeding programs, GS can be implemented in elite large germplasm panels with genotyped individuals to decrease PS efforts, thus reducing overall selection cycle time and accelerating variety development. However, determining correct GEBVs relies on prediction accuracy, which can vary according to a combination of phenotypic trait, genotypic information, and the statistical model. Consequently, GS research efforts have focused heavily on evaluating prediction accuracies of multiple parametric and non-parametric models to develop robust strategies that can be used in testing populations. Once robust modeling strategies are developed, GS has the potential to accumulate thousands of favorable alleles to develop climate-resilient crops with high yield potential. Additionally, as genotyping is required only once for a given population, multiple traits can be associated with the same genotypic matrix to determine GEBVs for each trait, thus making GS a valuable approach in multi-trait selection [76].
There is a need to explore new methodologies to improve molecular and bioinformatic tools for the application of GS in polyploid crops. The development of new approaches to obtain high-quality genome-wide markers will help to resolve the genetic architecture of complex traits. For example, the PRINCESS platform uses long-read sequencing to detect SNPs, indels, or methylation sites with high accuracy [77]. In GWAS and GS, the number of individuals in a population is crucial to maximize statistical power. Therefore, researchers search for genotyping methods that optimize the balance between cost, sample size, and the number of SNPs. In this regard, GBS is a relatively affordable genotyping methodology. However, for polyploid crops, there is a need for high-coverage sequencing (i.e., read-depth) to accurately estimate allele dosage, which increases genotyping cost. Biallelic SNPs are commonly used in polyploid GS because they are the most abundant type and are easier to transform into a numerical format for developing a genotypic matrix. However, during the construction of the vcf file, ~20% of high-quality SNP markers are discarded because they are not biallelic. Additionally, indels or simple sequence repeats (SSR) could add important information to GS model to increase prediction accuracy.
Parametric models such as RRBLUP assume that all SNPs have an effect on a specific trait, but the actual effect of each SNP is very small (heavy shrinkage) [9]. Although multiple loci have an effect on a complex trait, they often have different weights. Thus, identifying trait-specific weights for SNP marker alleles should increase prediction accuracy in GS models such as WGBLUP. Such outcomes have been demonstrated in animal breeding research [78,79], but not in polyploid crops. In this regard, utilizing a variable importance approach based on −log10 p-values for the additive GWASpoly model for alfalfa yield under salt stress was the best strategy to generate a diagonal matrix D . For this complex trait, we present empirical evidence demonstrating that the WGBLUP model increased prediction accuracy by almost 3 times compared to RRBLUP, GBLUP, Bayesian or ML models. Finally, we expanded the use of the WGBLUP model to 13 agronomic traits of potato and demonstrated the increase in prediction accuracy up to five times in complex traits, such as yield. The WGBLUP approach is simple, does not require high-performance computing, and can be applied to different crops to predict breeding values and accelerate selection cycles.

Author Contributions

C.A.M.: Manuscript preparation and data analysis; H.K.: Manuscript preparation and data analysis; I.R.: Conception and outline of this research and manuscript preparation; L.-X.Y.: Conception and outline of this research and manuscript preparation. All authors have read and agreed to the published version of the manuscript.

Funding

The authors would like to thank USDA-NIFA for funding support (2015-70005-24071).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Alfalfa raw data of GBS were submitted to the NCBI Sequence Read Archive with bioproject ID: PRJNA611554 and biosample # SAMN14336867 and phenotypic data can be retrieved from [6] (Table S2). Genotypic and phenotypic potato data were obtained from SolCAP potato diversity panel and have been published in [73] as Table S1.csv and Table S2.csv.

Acknowledgments

The authors would like to acknowledge Martha Rivera for technical help. They would further like to thank USDA-NIFA for funding support (2015-70005-24071).

Conflicts of Interest

The authors do not have conflict of interest.

References

  1. Blondon, F.; Marie, D.; Brown, S.; Kondorosi, A. Genome size and base composition in Medicago sativa and M. truncatula species. Genome 1994, 37, 264–270. [Google Scholar] [CrossRef]
  2. Elshire, R.J.; Glaubitz, J.C.; Sun, Q.; Poland, J.A.; Kawamoto, K.; Buckler, E.S.; Mitchell, S.E. A Robust, Simple Genotyping-by-Sequencing (GBS) Approach for High Diversity Species. PLoS ONE 2011, 6, e19379. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Yu, L.-X.; Liu, X.; Boge, W.; Liu, X.-P. Genome-Wide Association Study Identifies Loci for Salt Tolerance during Germination in Autotetraploid Alfalfa (Medicago sativa L.) Using Genotyping-by-Sequencing. Front. Plant Sci. 2016, 7, 956. [Google Scholar] [CrossRef] [Green Version]
  4. Liu, X.-P.; Yu, L.-X. Genome-Wide Association Mapping of Loci Associated with Plant Growth and Forage Production under Salt Stress in Alfalfa (Medicago sativa L.). Front. Plant Sci. 2017, 8, 853. [Google Scholar] [CrossRef]
  5. Liu, X.; Hawkins, C.; Peel, M.D.; Yu, L. Genetic Loci Associated with Salt Tolerance in Advanced Breeding Populations of Tetraploid Alfalfa Using Genome-Wide Association Studies. Plant Genome 2019, 12, 180026. [Google Scholar] [CrossRef]
  6. Medina, C.A.; Hawkins, C.; Liu, X.-P.; Peel, M.; Yu, L.-X. Genome-Wide Association and Prediction of Traits Related to Salt Tolerance in Autotetraploid Alfalfa (Medicago sativa L.). Int. J. Mol. Sci. 2020, 21, 3361. [Google Scholar] [CrossRef]
  7. Bulmer, M.G. The Effect of Selection on Genetic Variability. Am. Nat. 1971, 105, 201–211. [Google Scholar] [CrossRef]
  8. Meuwissen, T.H.E.; Hayes, B.J.; Goddard, M.E. Prediction of total genetic value using genome-wide dense marker maps. Genetics 2001, 157, 1819–1829. [Google Scholar] [CrossRef]
  9. Endelman, J.B. Ridge Regression and Other Kernels for Genomic Selection with R Package rrBLUP. Plant Genome 2011, 4, 250–255. [Google Scholar] [CrossRef] [Green Version]
  10. Hayes, B.J.; Bowman, P.J.; Chamberlain, A.C.; Verbyla, K.; Goddard, M.E. Accuracy of genomic breeding values in multi-breed dairy cattle populations. Genet. Sel. Evol. 2009, 41, 51. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Pérez, P.; Campos, G.D.L.; Crossa, J.; Gianola, D. Genomic-Enabled Prediction Based on Molecular Markers and Pedigree Using the Bayesian Linear Regression Package in R. Plant Genome 2010, 3, 106–116. [Google Scholar] [CrossRef]
  12. Vanraden, P.M. Genomic measures of relationship and inbreeding. Interbull Bull. 2007, 25, 33. [Google Scholar]
  13. VanRaden, P.M. Efficient Methods to Compute Genomic Predictions. J. Dairy Sci. 2008, 91, 4414–4423. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Yang, J.; Benyamin, B.; McEvoy, B.P.; Gordon, S.; Henders, A.; Nyholt, D.; Madden, P.A.; Heath, A.C.; Martin, N.; Montgomery, G.; et al. Common SNPs explain a large proportion of the heritability for human height. Nat. Genet. 2010, 42, 565–569. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Slater, A.T.; Cogan, N.O.; Forster, J.W.; Hayes, B.; Daetwyler, H.D. Improving Genetic Gain with Genomic Selection in Autotetraploid Potato. Plant Genome 2016, 9, 1–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Zhang, Z.; Liu, J.; Ding, X.; Bijma, P.; de Koning, D.-J.; Zhang, Q. Best Linear Unbiased Prediction of Genomic Breeding Values Using a Trait-Specific Marker-Derived Relationship Matrix. PLoS ONE 2010, 5, e12648. [Google Scholar] [CrossRef] [Green Version]
  17. Legarra, A.; Robert-Granié, C.; Croiseau, P.; Guillaume, F.; Fritz, S. Improved Lasso for genomic selection. Genet. Res. 2010, 93, 77–87. [Google Scholar] [CrossRef] [Green Version]
  18. Chang, L.-Y.; Toghiani, S.; Hay, E.H.; Aggrey, S.E.; Rekaya, R. A Weighted Genomic Relationship Matrix Based on Fixation Index (FST) Prioritized SNPs for Genomic Selection. Genes 2019, 10, 922. [Google Scholar] [CrossRef] [Green Version]
  19. Ren, D.; An, L.; Li, B.; Qiao, L.; Liu, W. Efficient weighting methods for genomic best linear-unbiased prediction (BLUP) adapted to the genetic architectures of quantitative traits. Heredity 2020, 126, 320–334. [Google Scholar] [CrossRef]
  20. Pérez, P.; de los Campos, G. Genome-Wide Regression and Prediction with the BGLR Statistical Package. Genetics 2014, 198, 483–495. [Google Scholar] [CrossRef]
  21. Meuwissen, T.H.; Solberg, T.R.; Shepherd, R.; Woolliams, J.A. A fast algorithm for BayesB type of prediction of genome-wide estimates of genetic value. Genet. Sel. Evol. 2009, 41, 2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Habier, D.; Fernando, R.L.; Kizilkaya, K.; Garrick, D.J. Extension of the bayesian alphabet for genomic selection. BMC Bioinform. 2011, 12, 186. [Google Scholar] [CrossRef] [Green Version]
  23. de los Campos, G.; Naya, H.; Gianola, D.; Crossa, J.; Legarra, A.; Manfredi, E.; Weigel, K.; Cotes, J.M. Predicting Quantitative Traits with Regression Models for Dense Molecular Markers and Pedigree. Genetics 2009, 182, 375–385. [Google Scholar] [CrossRef] [Green Version]
  24. Breiman, L. Statistical Modeling: The Two Cultures (with comments and a rejoinder by the author). Stat. Sci. 2001, 16, 199–231. [Google Scholar] [CrossRef]
  25. Drucker, H.; Surges, C.J.C.; Kaufman, L.; Smola, A.; Vapnik, V. Support vector regression machines. In Advances in Neural Information Processing Systems; Mozer, M.C., Jordan, M.I., Petsche, T., Eds.; MIT Press: Cambridge, MA, USA, 1997; pp. 155–161. [Google Scholar]
  26. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  27. Cortes, C.; Vapnik, V. Support-vector networks. Mach. Learn. 1995, 20, 273–297. [Google Scholar] [CrossRef]
  28. Liu, W.; Meng, X.; Xu, Q.; Flower, D.R.; Li, T. Quantitative prediction of mouse class I MHC peptide binding affinity using support vector machine regression (SVR) models. BMC Bioinform. 2006, 7, 182. [Google Scholar] [CrossRef] [Green Version]
  29. Sun, S.; Cao, Z.; Zhu, H.; Zhao, J. A Survey of Optimization Methods from a Machine Learning Perspective. IEEE Trans. Cybern. 2019, 50, 3668–3681. [Google Scholar] [CrossRef] [Green Version]
  30. LeCun, Y.; Bengio, Y.; Hinton, G. Deep learning. Nature 2015, 521, 436–444. [Google Scholar] [CrossRef]
  31. Pook, T.; Freudenthal, J.; Korte, A.; Simianer, H. Using Local Convolutional Neural Networks for Genomic Prediction. Front. Genet. 2020, 11, 1366. [Google Scholar] [CrossRef]
  32. Sandhu, K.S.; Lozada, D.N.; Zhang, Z.; Pumphrey, M.O.; Carter, A.H. Deep Learning for Predicting Complex Traits in Spring Wheat Breeding Program. Front. Plant Sci. 2021, 11, 2084. [Google Scholar] [CrossRef] [PubMed]
  33. Sandhu, K.; Patil, S.S.; Pumphrey, M.; Carter, A. Multitrait machine- and deep-learning models for genomic selection using spectral information in a wheat breeding program. Plant Genome 2021, 11, 1–17. [Google Scholar] [CrossRef] [PubMed]
  34. Zingaretti, L.M.; Gezan, S.A.; Ferrão, L.F.V.; Osorio, L.F.; Monfort, A.; Muñoz, P.R.; Whitaker, V.M.; Pérez-Enciso, M. Exploring Deep Learning for Complex Trait Genomic Prediction in Polyploid Outcrossing Species. Front. Plant Sci. 2020, 11, 25. [Google Scholar] [CrossRef] [Green Version]
  35. Ferrão, L.F.V.; Amadeu, R.R.; Benevenuto, J.; Oliveira, I.D.B.; Munoz, P.R. Genomic Selection in an Outcrossing Autotetraploid Fruit Crop: Lessons from Blueberry Breeding. Front. Plant Sci. 2021, 12, 1–13. [Google Scholar] [CrossRef] [PubMed]
  36. Klápště, J.; Dungey, H.S.; Telfer, E.J.; Suontama, M.; Graham, N.J.; Li, Y.; McKinley, R. Marker Selection in Multivariate Genomic Prediction Improves Accuracy of Low Heritability Traits. Front. Genet. 2020, 11, 499094. [Google Scholar] [CrossRef]
  37. Kyriakides, G.; Margaritis, K.G. Hands-on Ensemble Learning with Python: Build Highly Optimized Ensemble Machine Learning Models Using Scikit-Learn and Keras LK, 1st ed.; Packt Publishing Ltd.: Birmingham, UK, 2019; Available online: https://www.packtpub.com/product/hands-on-ensemble-learning-with-python/9781789612851 (accessed on 1 November 2021).
  38. Liang, M.; Chang, T.; An, B.; Duan, X.; Du, L.; Wang, X.; Miao, J.; Xu, L.; Gao, X.; Zhang, L.; et al. A Stacking Ensemble Learning Framework for Genomic Prediction. Front. Genet. 2021, 12, 600040. [Google Scholar] [CrossRef]
  39. Vos, P.G.; Uitdewilligen, J.G.A.M.L.; Voorrips, R.E.; Visser, R.G.F.; Van Eck, H.J. Development and analysis of a 20K SNP array for potato (Solanum tuberosum): An insight into the breeding history. Theor. Appl. Genet. 2015, 128, 2387–2401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Winfield, M.O.; Allen, A.M.; Burridge, A.; Barker, G.L.A.; Benbow, H.R.; Wilkinson, P.A.; Coghill, J.; Waterfall, C.; Davassi, A.; Scopes, G.; et al. High-density SNP genotyping array for hexaploid wheat and its secondary and tertiary gene pool. Plant Biotechnol. J. 2016, 14, 1195–1206. [Google Scholar] [CrossRef]
  41. Li, X.; Han, Y.; Wei, Y.; Acharya, A.; Farmer, A.D.; Ho, J.; Monteros, M.; Brummer, E.C. Development of an Alfalfa SNP Array and Its Use to Evaluate Patterns of Population Structure and Linkage Disequilibrium. PLoS ONE 2014, 9, e84329. [Google Scholar] [CrossRef] [Green Version]
  42. Perkel, J. SNP genotyping: Six technologies that keyed a revolution. Nat. Chem. Biol. 2008, 5, 447–453. [Google Scholar] [CrossRef]
  43. Clark, L.V.; Lipka, A.E.; Sacks, E.J. polyRAD: Genotype Calling with Uncertainty from Sequencing Data in Polyploids and Diploids. G3 Genes Genomes Genet. 2019, 9, 663–673. [Google Scholar] [CrossRef] [Green Version]
  44. Pereira, G.S.; Garcia, A.A.F.; Margarido, G.R.A. A fully automated pipeline for quantitative genotype calling from next generation sequencing data in autopolyploids. BMC Bioinform. 2018, 19, 398. [Google Scholar] [CrossRef] [Green Version]
  45. Zych, K.; Gort, G.; Maliepaard, C.A.; Jansen, R.C.; Voorrips, R.E. FitTetra 2.0—Improved genotype calling for tetraploids with multiple population and parental data support. BMC Bioinform. 2019, 20, 148. [Google Scholar] [CrossRef] [Green Version]
  46. Gerard, D.; Ferrão, L.F.V.; Garcia, A.A.F.; Stephens, M. Genotyping Polyploids from Messy Sequencing Data. Genetics 2018, 210, 789–807. [Google Scholar] [CrossRef] [Green Version]
  47. Uitdewilligen, J.G.A.M.L.; Wolters, A.-M.; D’Hoop, B.B.; Borm, T.J.A.; Visser, R.G.F.; Van Eck, H.J. A Next-Generation Sequencing Method for Genotyping-by-Sequencing of Highly Heterozygous Autotetraploid Potato. PLoS ONE 2013, 8, e62355. [Google Scholar] [CrossRef] [Green Version]
  48. Amadeu, R.R.; Ferrão, L.F.V.; Oliveira, I.D.B.; Benevenuto, J.; Endelman, J.B.; Munoz, P.R. Impact of Dominance Effects on Autotetraploid Genomic Prediction. Crop Sci. 2019, 60, 656–665. [Google Scholar] [CrossRef]
  49. Endelman, J.B.; Carley, C.A.S.; Bethke, P.C.; Coombs, J.J.; Clough, M.E.; Silva, W.L.; De Jong, W.S.; Douches, D.S.; Frederick, C.M.; Haynes, K.G.; et al. Genetic Variance Partitioning and Genome-Wide Autotetraploid Potato. Genetics 2018, 209, 77–87. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Batista, L.G.; Mello, V.H.; Souza, A.P.; Margarido, G.R.A. Genomic prediction with allele dosage information in highly polyploid species. BioRxiv 2021. [Google Scholar] [CrossRef]
  51. Oliveira, I.D.B.; Resende, M.F.R.; Ferrão, L.F.V.; Amadeu, R.R.; Endelman, J.B.; Kirst, M.; Coelho, A.S.G.; Muñoz, P.R. Genomic Prediction of Autotetraploids; Influence of Relationship Matrices, Allele Dosage, and Continuous Genotyping Calls in Phenotype Prediction. G3 Genes Genomes Genet. 2019, 9, 1189–1198. [Google Scholar] [CrossRef] [Green Version]
  52. Oliveira, I.D.B.; Amadeu, R.R.; Ferrão, L.F.V.; Muñoz, P.R. Optimizing whole-genomic prediction for autotetraploid blueberry breeding. Heredity 2020, 125, 437–448. [Google Scholar] [CrossRef]
  53. Jia, C.; Zhao, F.; Wang, X.; Han, J.; Zhao, H.; Liu, G.; Wang, Z. Genomic Prediction for 25 Agronomic and Quality Traits in Alfalfa (Medicago sativa). Front. Plant Sci. 2018, 9, 1220. [Google Scholar] [CrossRef]
  54. Biazzi, E.; Nazzicari, N.; Pecetti, L.; Brummer, E.C.; Palmonari, A.; Tava, A.; Annicchiarico, P. Genome-Wide Association Mapping and Genomic Selection for Alfalfa (Medicago sativa) Forage Quality Traits. PLoS ONE 2017, 12, e0169234. [Google Scholar] [CrossRef] [PubMed]
  55. Campbell, M.T.; Hu, H.; Yeats, T.H.; Brzozowski, L.J.; Caffe-Treml, M.; Gutiérrez, L.; Smith, K.P.; Sorrells, M.E.; Gore, M.A.; Jannink, J.-L. Improving Genomic Prediction for Seed Quality Traits in Oat (Avena sativa L.) Using Trait-Specific Relationship Matrices. Front. Genet. 2021, 12, 437. [Google Scholar] [CrossRef]
  56. Fikere, M.; Barbulescu, D.M.; Malmberg, M.M.; Maharjan, P.; Salisbury, P.A.; Kant, S.; Panozzo, J.; Norton, S.; Spangenberg, G.C.; Cogan, N.O.I.; et al. Genomic Prediction and Genetic Correlation of Agronomic, Blackleg Disease, and Seed Quality Traits in Canola (Brassica napus L.). Plants 2020, 9, 719. [Google Scholar] [CrossRef]
  57. Sousa, T.; Caixeta, E.T.; Alkimim, E.; Oliveira, A.C.B.; Pereira, A.A.; Sakiyama, N.S.; Zambolim, L.; Resende, M.D.V. Early Selection Enabled by the Implementation of Genomic Selection in Coffea arabica Breeding. Front. Plant Sci. 2019, 9, 1934. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Li, X.; Wei, Y.; Acharya, A.; Hansen, J.L.; Crawford, J.L.; Viands, D.R.; Michaud, R.; Claessens, A.; Brummer, E.C. Genomic Prediction of Biomass Yield in Two Selection Cycles of a Tetraploid Alfalfa Breeding Population. Plant Genome 2015, 8, 1–10. [Google Scholar] [CrossRef]
  59. Annicchiarico, P.; Nazzicari, N.; Li, X.; Wei, Y.; Pecetti, L.; Brummer, E.C. Accuracy of genomic selection for alfalfa biomass yield in different reference populations. BMC Genom. 2015, 16, 1–13. [Google Scholar] [CrossRef] [Green Version]
  60. Lara, L.A.C.; Santos, M.F.; Jank, L.; Chiari, L.; Vilela, M.M.; Amadeu, R.R.; dos Santos, J.P.R.; Pereira, G.A.; Zeng, Z.-B.; Garcia, A.A.F. Genomic Selection with Allele Dosage in Panicum maximum Jacq. G3 Genes Genomes Genet. 2019, 9, 2463–2475. [Google Scholar] [CrossRef] [Green Version]
  61. Wilson, S.; Zheng, C.; Maliepaard, C.; Mulder, H.A.; Visser, R.G.F.; van der Burgt, A.; van Eeuwijk, F. Understanding the Effectiveness of Genomic Prediction in Tetraploid Potato. Front. Plant Sci. 2021, 12, 1634. [Google Scholar] [CrossRef]
  62. Yadav, S.; Wei, X.; Joyce, P.; Atkin, F.; Deomano, E.; Sun, Y.; Nguyen, L.T.; Ross, E.M.; Cavallaro, T.; Aitken, K.S.; et al. Improved genomic prediction of clonal performance in sugarcane by exploiting non-additive genetic effects. Theor. Appl. Genet. 2021, 134, 2235–2252. [Google Scholar] [CrossRef] [PubMed]
  63. Michel, S.; Löschenberger, F.; Ametz, C.; Pachler, B.; Sparry, E.; Bürstmayr, H. Simultaneous selection for grain yield and protein content in genomics-assisted wheat breeding. Theor. Appl. Genet. 2019, 132, 1745–1760. [Google Scholar] [CrossRef] [PubMed]
  64. Sehgal, D.; Rosyara, U.; Mondal, S.; Singh, R.; Poland, J.; Dreisigacker, S. Incorporating Genome-Wide Association Mapping Results Into Genomic Prediction Models for Grain Yield and Yield Stability in CIMMYT Spring Bread Wheat. Front. Plant Sci. 2020, 11, 197. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Aparicio Arce, J.S. Mr. Bean. 2018. Available online: https://apariciojohan.shinyapps.io/Mrbean/ (accessed on 18 March 2020).
  66. Rodríguez-Álvarez, M.X.; Boer, M.P.; van Eeuwijk, F.A.; Eilers, P.H.C. Spatial Models for Field Trials. arXiv 2016, arXiv:1607.08255. [Google Scholar]
  67. Isik, F.; Holland, J.; Maltecca, C. Multi Environmental Trials. In Genetic Data Analysis for Plant and Animal Breeding; Isik, F., Holland, J., Maltecca, C., Eds.; Springer International Publishing: Cham, Switzerland, 2017; pp. 227–262. [Google Scholar] [CrossRef]
  68. Butler, D.G.; Cullis, B.R.; Gilmour, A.R.; Gogel, B.J.; Thompson, R. ASReml-R Reference Manual Version 4, ASReml-R Ref. Man. 2018. Available online: http://www.homepages.ed.ac.uk/iwhite/asreml/uop (accessed on 1 November 2021).
  69. Duitama, J.; Quintero, J.C.; Cruz, D.F.; Quintero, C.; Hubmann, G.; Foulquié-Moreno, M.R.; Verstrepen, K.; Thevelein, J.; Tohme, J. An integrated framework for discovery and genotyping of genomic variants from high-throughput sequencing experiments. Nucleic Acids Res. 2014, 42, e44. [Google Scholar] [CrossRef] [PubMed]
  70. Covarrubias-Pazaran, G. Genome-Assisted Prediction of Quantitative Traits Using the R Package sommer. PLoS ONE 2016, 11, e0156744. [Google Scholar] [CrossRef] [Green Version]
  71. Bernardo, R.; Yu, J. Prospects for Genomewide Selection for Quantitative Traits in Maize. Crop Sci. 2007, 47, 1082–1090. [Google Scholar] [CrossRef] [Green Version]
  72. Kuhn, M.; Wing, J.; Weston, S.; Williams, A.; Keefer, C.; Engelhardt, A.; Cooper, T.; Mayer, Z.; Kenkel, B.; The R.C. Team; et al. Caret: Classification and Regression Training. 2019. Available online: https://cran.r-project.org/package=caret (accessed on 20 March 2020).
  73. Rosyara, U.R.; De Jong, W.S.; Douches, D.S.; Endelman, J.B. Software for Genome-Wide Association Studies in Autopolyploids and Its Application to Potato. Plant Genome 2016, 9, 1–10. [Google Scholar] [CrossRef] [Green Version]
  74. Amadeu, R.R.; Cellon, C.; Olmstead, J.W.; Garcia, A.A.F.; Resende, M.F.R.; Muñoz, P.R. AGHmatrix: R Package to Construct Relationship Matrices for Autotetraploid and Diploid Species: A Blueberry Example. Plant Genome 2016, 9, 9 . [Google Scholar] [CrossRef]
  75. Tessema, B.B.; Liu, H.; Sørensen, A.C.; Andersen, J.R.; Jensen, J. Strategies Using Genomic Selection to Increase Genetic Gain in Breeding Programs for Wheat. Front. Genet. 2020, 11, 578123. [Google Scholar] [CrossRef]
  76. Moeinizade, S.; Kusmec, A.; Hu, G.; Wang, L.; Schnable, P. Multi-trait Genomic Selection Methods for Crop Improvement. Genetics 2020, 215, 931–945. [Google Scholar] [CrossRef] [PubMed]
  77. Mahmoud, M.; Doddapaneni, H.; Timp, W.; Sedlazeck, F.J. PRINCESS: Comprehensive detection of haplotype resolved SNVs, SVs, and methylation. Genome Biol. 2021, 22, 268. [Google Scholar] [CrossRef] [PubMed]
  78. Lopez, B.; Lee, S.-H.; Park, J.-E.; Shin, D.-H.; Oh, J.-D.; Heras-Saldana, S.D.L.; Van Der Werf, J.; Chai, H.-H.; Park, W.; Lim, D. Correction: Weighted Genomic Best Linear Unbiased Prediction for Carcass Traits in Hanwoo Cattle. Genes 2020, 11, 1013. [Google Scholar] [CrossRef] [PubMed]
  79. Zhang, X.; Lourenco, D.; Aguilar, I.; Legarra, A.; Misztal, I. Weighting Strategies for Single-Step Genomic BLUP: An Iterative Approach for Accurate Calculation of GEBV and GWAS. Front. Genet. 2016, 7, 151. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Indirect selection based on molecular markers. (a) Generalized Manhattan plots illustrating a comparison of GWAS effectiveness in simple (left) vs. complex traits (right). Note: Bold dashed line indicates minimum threshold to select significant markers. A significant signal (i.e., QTL) was identified in the simple trait (left panel), while no defined QTL was identified for the complex trait. Therefore, genomic selection (GS) is more appropriate and practical for complex traits. (b) Common parametric and non-parametric models used in GS and their computational requirements. GBLUP, genomic best linear unbiased prediction; RRBLUP, ridge-regression BLUP; RF, random forest; SVM, support vector machine; MLP, multilayer perceptron; CNN, convolutional neural network; RNN, recurrent neural network.
Figure 1. Indirect selection based on molecular markers. (a) Generalized Manhattan plots illustrating a comparison of GWAS effectiveness in simple (left) vs. complex traits (right). Note: Bold dashed line indicates minimum threshold to select significant markers. A significant signal (i.e., QTL) was identified in the simple trait (left panel), while no defined QTL was identified for the complex trait. Therefore, genomic selection (GS) is more appropriate and practical for complex traits. (b) Common parametric and non-parametric models used in GS and their computational requirements. GBLUP, genomic best linear unbiased prediction; RRBLUP, ridge-regression BLUP; RF, random forest; SVM, support vector machine; MLP, multilayer perceptron; CNN, convolutional neural network; RNN, recurrent neural network.
Cells 10 03372 g001
Figure 2. Optimization of GS models. (a) GS model accuracy measured as Pearson’s correlation after 10-fold cross-validation for biomass yield under salt stress. Computing time was measured as system time in seconds to run one cross-validation. (b) Example of variable importance values derived from SVM for 10 randomly chosen SNPs. (c) Pearson’s correlation for 6796 SNPs weights obtained by variable importance (SVM, RF) or by −log10 p-values of different GWASpoly models. (d) Accuracy of GBLUP (GBLUP VR and GBLUP FA) and WGBLUP models. Accuracy was measured 10 times using Pearson’s correlation with 10-fold cross-validation. SNP weights for WGBLUP were obtained from variable importance values (SVM, RF) or −log10 p-values of different GWASpoly models. RRBLUP, best linear unbiased prediction using ridge-regression; BL Bayes LASSO; GBLUP, genomic best linear unbiased prediction; VR, VanRaden G matrix; FA, full autotetraploid G matrix; RF, random forest; SVM, support vector machine; WGBLUP, weighted GBLUP; 1-dom-alt and 1-dom-ref, simplex dominant models; 2-dom-alt and 2-dom-ref, duplex dominant models; diplo-general, diploidized general; diplo-additive, diploidized additive.
Figure 2. Optimization of GS models. (a) GS model accuracy measured as Pearson’s correlation after 10-fold cross-validation for biomass yield under salt stress. Computing time was measured as system time in seconds to run one cross-validation. (b) Example of variable importance values derived from SVM for 10 randomly chosen SNPs. (c) Pearson’s correlation for 6796 SNPs weights obtained by variable importance (SVM, RF) or by −log10 p-values of different GWASpoly models. (d) Accuracy of GBLUP (GBLUP VR and GBLUP FA) and WGBLUP models. Accuracy was measured 10 times using Pearson’s correlation with 10-fold cross-validation. SNP weights for WGBLUP were obtained from variable importance values (SVM, RF) or −log10 p-values of different GWASpoly models. RRBLUP, best linear unbiased prediction using ridge-regression; BL Bayes LASSO; GBLUP, genomic best linear unbiased prediction; VR, VanRaden G matrix; FA, full autotetraploid G matrix; RF, random forest; SVM, support vector machine; WGBLUP, weighted GBLUP; 1-dom-alt and 1-dom-ref, simplex dominant models; 2-dom-alt and 2-dom-ref, duplex dominant models; diplo-general, diploidized general; diplo-additive, diploidized additive.
Cells 10 03372 g002
Table 1. Different prior distributions for Bayesian models.
Table 1. Different prior distributions for Bayesian models.
ModelPrior Distribution Ref.
Bayes A β j ~ t ( d f β , S β ) [8]
Bayes B β j = { 1 / 2 γ λ exp ( λ | β j | ) ( 1 γ ) f o r   β j 0 f o r   β j = 0 [21]
Bayes Cπ β j | π ,   σ β j 2 { β j ~ 0 β j ~ N ( 0 , σ β j 2 ) w i t h   p r o b   π                         w i t h   p r o b   ( 1 π ) [22]
Bayesian LASSO β j ~ D E ( λ 2 , σ e 2 ) [23]
; β j , is the additive effect of the j t h ; t , scaled-t distribution; d f β , degree of freedom; S β , scale parameters; γ , fraction of the SNPs that are in linkage disequilibrium with a quantitative trait locus; SNP; π , probability of the marker effect equal to zero; D E , double exponential; λ , parameter of exponential distribution.
Table 2. Kernels used in support vector machine (SVM) model. Meta-parameters used for tuning include gamma ( γ ), degree of polynomial ( d ) and intercept ( α ).
Table 2. Kernels used in support vector machine (SVM) model. Meta-parameters used for tuning include gamma ( γ ), degree of polynomial ( d ) and intercept ( α ).
KernelFormula
Linear K ( x i , y j ) = x i T y j
Polynomial K ( x i , y j ) = γ ( x i T y j + α ) d
Radial basis function K ( x i , y j ) = e γ x i y j 2
Sigmoidal K ( x i , y j ) = tan h ( γ x i T y j + α )
; x i , y j are two vectors in the n-dimensional space.
Table 3. Coding effect assumptions of GWASpoly models according to allele dosage in biallelic SNPs.
Table 3. Coding effect assumptions of GWASpoly models according to allele dosage in biallelic SNPs.
Allele Dosage AAAAAAABAABBABBBBBBB
Numerical Code01234
GWASpoly ModelsPhenotypic Effect §
Diplo-additive0.000.501.00
Diplo-general 0.000.00 < x <1.001.00
Additive0.000.250.500.751.00
1-dom-ref (A > B simplex)1.001.001.001.000.00
2-dom-ref (A > B duplex)1.001.001.000.000.00
1-dom-alt (B > A simplex)0.001.001.001.001.00
2-dom-alt (B > A duplex)0.000.001.001.001.00
General No restrictions
, allele dosage A is coded as the reference allele and B is coded as the alternative allele; §, phenotypic effects are scaled from 0.00 to 1.00; , for the diplo-general model all heterozygotes have the same effect (x), but x is not constrained to be halfway between the homozygous effects; , the general model has no restrictions on the effects of the different dosage levels.
Table 4. Recent achievements in genomic selection (GS) in polyploid crops.
Table 4. Recent achievements in genomic selection (GS) in polyploid crops.
CropPloidyTrait §GS MethodAcc NotesAuthor
Avena sativaAllohexaploidSeed lipid contentMK-BLUP0.48Use of additive marker effects of Bayesian models during the construction of G matrix[55]
Brassica napusAlloteteraploidSeed yieldGBLUP0.69Several agronomic and seed quality traits were tested[56]
Coffea arabicaAllotetraploidCanopy diameterGBLUP0.4018 agronomic traits were tested. Diploid dosage assumed[57]
Eucalyptus nitensPaleotetraploidWood densityMVGLUP 0.77Marker selection in multivariate analysis. Requires uses multiple traits highly correlated[36]
Medicago sativaAutotetraploidYieldRRBLUP0.66Multi-environment trials over two generations. First report of GS in alfalfa.[58]
Medicago sativaAutotetraploidYieldSVM0.35Six GS models were tested. First report of machine learning models in alfalfa[59]
Medicago sativaAutotetraploidLeaf crude proteinRRBLUP0.40Nine alfalfa forage quality traits were tested by five GS models[54]
Medicago sativaAutotetraploidFall plant heightBayes B0.6515 quality traits and 10 agronomic traits were tested using three GS models[53]
Medicago sativaAutotetraploidYield under salt stressSVM0.50Multi-environment trials with seven yield measurements. Eight GS models were tested[6]
Panicum maximumAutotetraploidOrganic matterBayes B-TD0.39Genomic selection using tetraploid dosage (GS-TD) vs. diploid dosage (GS-DD)[60]
Solanum tuberosumAutopolyploidYieldGBLUP0.55Incorporation of additive and digenic dominant G covariance matrix[49]
Solanum tuberosumAutopolyploidTuber weightRKHS0.59Four agronomic tuber traits were tested by eight GS models[61]
SugarcaneOctaploid and decaploidFiberGBLUP0.44Inclusion of additive and non-additive genetic components for GS[62]
Triticum aestivumAllohexaploidGrain yieldGBLUP0.47Multi-trait selection for grain yield and protein content[63]
Triticum aestivumAllohexaploidGrain yieldGBLUP0.53GWAS markers as fixed effects in GS models.[64]
Vaccinium corymbosumAutotetraploidWeightGBLUP0.49Comparison of allele dosage with depth sequencing: 6×–60×)[35]
§ For multiple traits, the trait with the highest predictive accuracy was selected; , predictive accuracy measured as Pearson’s correlation; MK-BLUP, multi-kernel trait-specific BLUP; MVGLUP, Multi-trait model GBLUP; SVM, support vector machine; Bayes B-TD, Bayes B with tetraploid allele dosage; RKHS, Reproducing Kernel Hilbert Space; In multi-trait genomic selection (MT-GS) a secondary trait that is genetically correlated with the primary trait is incorporated in the prediction model, to predict the primary trait with higher accuracy.
Table 5. Comparison of genomic selection (GS) models in 13 phenotypic traits collected in the SolCAP potato diversity panel. Mean and standard deviation of Pearson’s correlation obtained by 10-fold cross validation in 10 replicates. SNP weights for WGBLUP were obtained from −log10 p-values of different GWASpoly models.
Table 5. Comparison of genomic selection (GS) models in 13 phenotypic traits collected in the SolCAP potato diversity panel. Mean and standard deviation of Pearson’s correlation obtained by 10-fold cross validation in 10 replicates. SNP weights for WGBLUP were obtained from −log10 p-values of different GWASpoly models.
TraitRRBLUPGBLUPWGBLUP
1-d-a1-d-r2-d-a2-d-rGenerald-Gend-AddAdditive
Chip color0.7230.7210.8260.7980.8590.8500.8670.8490.8550.896
(±0.014)(±0.015)(±0.009)(±0.011)(±0.007)(±0.013)(±0.008)(±0.009)(±0.007)(±0.007)
log10 fructose0.6820.6760.8190.7850.8450.8330.8680.8390.8550.895
(±0.024)(±0.025)(±0.014)(±0.017)(±0.007)(±0.011)(±0.011)(±0.015)(±0.003)(±0.008)
log10 glucose0.6780.6680.7960.8090.8550.8490.8750.8440.8480.91
(±0.017)(±0.030)(±0.009)(±0.016)(±0.009)(±0.009)(±0.009)(±0.011)(±0.013)(±0.007)
Malic acid0.6020.5980.7510.7450.8020.8010.8380.8080.8260.876
(±0.016)(±0.027)(±0.021)(±0.022)(±0.021)(±0.016)(±0.011)(±0.016)(±0.009)(±0.007)
Sucrose0.5390.5190.6760.6750.7020.7160.7250.7220.7390.806
(±0.024)(±0.034)(±0.011)(±0.022)(±0.019)(±0.015)(±0.023)(±0.011)(±0.019)(±0.011)
Total yield0.1320.1170.4010.4130.4180.4280.4700.4920.5040.584
(±0.023)(±0.041)(±0.026)(±0.030)(±0.031)(±0.017)(±0.029)(±0.030)(±0.030)(±0.028)
Tuber eye depth0.4950.4780.6050.6550.6930.7170.7400.6930.7360.812
(±0.026)(±0.019)(±0.029)(±0.016)(±0.025)(±0.014)(±0.020)(±0.020)(±0.018)(±0.007)
Tuber length0.8260.8210.8910.8840.8990.8890.9040.9080.9120.928
(±0.012)(±0.014)(±0.006)(±0.009)(±0.006)(±0.012)(±0.008)(±0.008)(±0.005)(±0.009)
Tuber shape0.7750.7800.8650.8530.8860.8630.8960.890.8910.922
(±0.018)(±0.017)(±0.010)(±0.013)(±0.008)(±0.005)(±0.010)(±0.008)(±0.009)(±0.006)
Tuber size0.5010.4990.6410.6500.6790.6630.6660.6610.6790.742
(±0.024)(±0.027)(±0.019)(±0.020)(±0.020)(±0.022)(±0.024)(±0.022)(±0.019)(±0.021)
Tuber width0.6350.6380.7520.7490.7820.7720.8050.7890.8030.847
(±0.023)(±0.021)(±0.020)(±0.021)(±0.016)(±0.018)(±0.012)(±0.015)(±0.013)(±0.017)
Vine maturity 95 days0.2880.2860.5500.5380.6030.5890.6680.6320.650.746
(±0.035)(±0.042)(±0.028)(±0.020)(±0.022)(±0.028)(±0.022)(±0.019)(±0.025)(±0.017)
Vine maturity 120 days0.3210.3230.4950.5690.6360.6330.6690.6160.6660.755
(±0.047)(±0.024)(±0.026)(±0.021)(±0.021)(±0.013)(±0.025)(±0.023)(±0.026)(±0.019)
RRBLUP, best linear unbiased prediction using ridge-regression; GBLUP, genomic best linear unbiased prediction using VanRaden G matrix; WGBLUP, weighted GBLUP; 1-d-a and 1-d-r, simplex dominant models; 2-d-a and 2-d-r, duplex dominant models; d-gen, diploidized general; d-add, diploidized additive.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Medina, C.A.; Kaur, H.; Ray, I.; Yu, L.-X. Strategies to Increase Prediction Accuracy in Genomic Selection of Complex Traits in Alfalfa (Medicago sativa L.). Cells 2021, 10, 3372. https://doi.org/10.3390/cells10123372

AMA Style

Medina CA, Kaur H, Ray I, Yu L-X. Strategies to Increase Prediction Accuracy in Genomic Selection of Complex Traits in Alfalfa (Medicago sativa L.). Cells. 2021; 10(12):3372. https://doi.org/10.3390/cells10123372

Chicago/Turabian Style

Medina, Cesar A., Harpreet Kaur, Ian Ray, and Long-Xi Yu. 2021. "Strategies to Increase Prediction Accuracy in Genomic Selection of Complex Traits in Alfalfa (Medicago sativa L.)" Cells 10, no. 12: 3372. https://doi.org/10.3390/cells10123372

APA Style

Medina, C. A., Kaur, H., Ray, I., & Yu, L. -X. (2021). Strategies to Increase Prediction Accuracy in Genomic Selection of Complex Traits in Alfalfa (Medicago sativa L.). Cells, 10(12), 3372. https://doi.org/10.3390/cells10123372

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop