Application Of Genomic Selection Research Articles

306 Application of genomic selection for enhancing health, welfare, efficiency and production traits in turkeys.

The implementation of single-step genomic BLUP (ssGBLUP) to selection programs can considerably increase the accuracy of breeding values, and therefore the ability to estimate the genetic merit, of livestock. In particular, traits with low heritability, such as livability and other health and welfare traits, are amenable to genomic selection. Effective application of ssGBLUP requires genomic information, population-based pedigree data, and high-quality phenotypes. The aim of this research is to apply ssGBLUP to commercial turkey lines (15,000 turkeys genotyped at 65K). This population will be used to assess the increase in accuracy of selection using genomic information, and to monitor inbreeding within the different lines. In addition to the high-quality phenotypes available related to fertility, growth, production, and carcass composition, auxiliary health and behaviour phenotypes related to livability will be investigated. Furthermore, novel meat quality (e.g. water holding capacity, pH, sheer force and colour) and total carcass composition phenotypes will also be identified and collected.With improved methodology, more detailed phenotypic information, and comprehensive data collection and integration, more accurate selection of parent stock will be possible.

Karyomapping for simultaneous genomic evaluation and aneuploidy screening of preimplantation bovine embryos: The first live-born calves

In cattle breeding, the development of genomic selection strategies based on single nucleotide polymorphism (SNP) interrogation has led to improved rates of genetic gain. Additionally, the application of genomic selection to in-vitro produced (IVP) embryos is expected to bring further benefits thanks to the ability to test a greater number of individuals before establishing a pregnancy and to ensure only carriers of desirable traits are born. However, aneuploidy, a leading cause of developmental arrest, is known to be common in IVP embryos. Karyomapping is a comprehensive screening test based on SNP typing that can be used for simultaneous genomic selection and aneuploidy detection, offering the potential to maximize pregnancy rates. Moreover, Karyomapping can be used to characterize the frequency and parental origin of aneuploidy in bovine IVP embryos, which have remained underexplored to date. Here, we report the use of Karyomapping to characterize the frequency and parental origin of aneuploidy in IVP bovine embryos in order to establish an estimate of total aneuploidy rates in each parental germline. We report an estimate of genome wide recombination rate in cattle and demonstrate, for the first time, a proof of principle for the application of Karyomapping to cattle breeding, with the birth of five calves after screening. This combined genomic selection and aneuploidy screening approach was highly reliable, with calves showing 98% concordance with their respective embryo biopsies for SNP typing and 100% concordance with their respective biopsies for aneuploidy screening. This approach has the potential to simultaneously improve pregnancy rates following embryo transfer and the rate of genetic gain in cattle breeding, and is applicable to basic research to investigate meiosis and aneuploidy.

Genomic selection efficiency and a priori estimation of accuracy in a structured dent maize panel.

Population structure affects genomic selection efficiency as well as the ability to forecast accuracy using standard GBLUP. Genomic prediction models usually assume that the individuals used for calibration belong to the same population as those to be predicted. Most of the a priori indicators of precision, such as the coefficient of determination (CD), were derived from those same models. But genetic structure is a common feature in plant species, and it may impact genomic selection efficiency and the ability to forecast prediction accuracy. We investigated the impact of genetic structure in a dent maize panel ("Amaizing Dent") using different scenarios including within- or across-group predictions. For a given training set size, the best accuracies were achieved when predicting individuals using a model calibrated on the same genetic group. Nevertheless, a diverse training set representing all the groups had a certain predictive efficiency for all the validation sets, and adding extra-group individuals was almost always beneficial. It underlines the potential of sucha generic training set for dent maize genomic selection applications. Alternative prediction models, taking genetic structure explicitly into account, did not improve the prediction accuracy compared to GBLUP. We also investigated the ability of different indicators of precision to forecast accuracy in the within- or across-group scenarios. There was a global encouraging trend of the CD to differentiate scenarios, although there were specific combinations of target populations and traits where the efficiency of this indicator proved to be null. One hypothesis to explain such erratic performances is the impact of genetic structure through group-specific allele diversity at QTLs rather than group-specific allele effects.

Genomic selection prediction models comparing sequence capture and SNP array genotyping methods

The successful application of genomic selection (GS) approaches is dependent on genetic makers derived from high-throughput and low-cost genotyping methods. Recent GS studies in trees have predominantly relied on SNP arrays as the source of genotyping, though this technology has a high entry cost. The recent development of alternative genotyping platforms, tailored to specific species and with low entry cost, has become possible due to advances in next-generation sequencing and genome complexity reduction methods such as sequence capture. However, the performance of these new platforms in GS models has not yet been evaluated, or compared to models developed from SNP arrays. Here, we evaluate the impact of these genotyping technologies on the development of GS prediction models for a Eucalyptus breeding population composed of 739 trees phenotyped for 13 wood quality and growth traits. Genotyping data obtained with both methods were compared for linkage disequilibrium, minor allele frequency, and missing data. Phenotypic prediction methods RR-BLUP and BayesB were employed, while predictive ability using cross validation was used to evaluate the performance of GS models derived from the different genotyping platforms. Differences in linkage disequilibrium patterns, minor allele frequency, missing data, and marker distribution were detected between sequence capture and SNP arrays. However, RR-BLUP and BayesB GS models resulted in similar predictive abilities. These results demonstrate that both genotyping methods are equivalent for genomic prediction of the traits evaluated. Sequence capture offers an alternative for species where SNP arrays are not available, or for when the initial development cost is too high.

A comparison between genotyping-by-sequencing and array-based scoring of SNPs for genomic prediction accuracy in winter wheat

The utilization of DNA molecular markers in plant breeding to maximize selection response via marker-assisted selection (MAS) and genomic selection (GS) has revolutionized plant breeding. A key factor affecting GS applicability is the choice of molecular marker platform. Genotyping-by-sequencing scored SNPs (GBS-scored SNPs) provides a large number of markers, albeit with high rates of missing data. Array scored SNPs are of high quality, but the cost per sample is substantially higher. The objectives of this study were 1) compare GBS-scored SNPs, and array scored SNPs for genomic selection applications, and 2) compare estimates of genomic kinship and population structure calculated using the two marker platforms. SNPs were compared in a diversity panel consisting of 299 hard winter wheat (Triticum aestivum L.) accessions that were part of a multi-year, multi-environments association mapping study. The panel was phenotyped in Ithaca, Nebraska for heading date, plant height, days to physiological maturity and grain yield in 2012 and 2013. The panel was genotyped using GBS-scored SNPs, and array scored SNPs. Results indicate that GBS-scored SNPs is comparable to or better than Array-scored SNPs for genomic prediction application.Both platforms identified the same genetic patterns in the panel where 90% of the lines were classified to common genetic groups. Overall, we concluded that GBS-scored SNPs have the potential to be the marker platform of choice for genetic diversity and genomic selection in winter wheat.

Genomic prediction using DArT-Seq technology for yellowtail kingfish Seriola lalandi

BackgroundGenomic prediction using Diversity Arrays Technology (DArT) genotype by sequencing platform has not been reported in yellowtail kingfish (Seriola lalandi). The principal aim of this study was to address this knowledge gap and to assess predictive ability of genomic Best Linear Unbiased Prediction (gBLUP) for traits of commercial importance in a yellowtail kingfish population comprising 752 individuals that had DNA sequence and phenotypic records for growth traits (body weight, fork length and condition index). The gBLUP method was used due to its computational efficiency and it showed similar predictive performance to other approaches, especially for traits whose variation is of polygenic nature, such as body traits analysed in this study. The accuracy or predictive ability of the gBLUP model was estimated for three growth traits: body weight, folk length and condition index.ResultsThe prediction accuracy was moderate to high (0.44 to 0.69) for growth-related traits. The predictive ability for body weight increased by 17.0% (from 0.69 to 0.83) when missing genotype was imputed. Within population prediction using five-fold across validation approach showed that the gBLUP model performed well for growth traits (weight, length and condition factor), with the coefficient of determination (R2) from linear regression analysis ranging from 0.49 to 0.71.ConclusionsCollectively our results demonstrated, for the first time in yellowtail kingfish, the potential application of genomic selection for growth-related traits in the future breeding program for this species, S. lalandi.

Linkage disequilibrium levels and allele frequency distribution in Blanco Orejinegro and Romosinuano Creole cattle using medium density SNP chip data.

The linkage disequilibrium (LD) between molecular markers affects the accuracy of genome-wide association studies and genomic selection application. High-density genotyping platforms allow identifying the genotype of thousands of single nucleotide polymorphisms (SNPs) distributed throughout the animal genomes, which increases the resolution of LD evaluations. This study evaluated the distribution of minor allele frequencies (MAF) and the level of LD in the Colombian Creole cattle breeds Blanco Orejinegro (BON) and Romosinuano (ROMO) using a medium density SNP panel (BovineSNP50K_v2). The LD decay in these breeds was lower than those reported for other taurine breeds, achieving optimal LD values (r2 ≥ 0.3) up to a distance of 70 kb in BON and 100 kb in ROMO, which is possibly associated with the conservation status of these cattle populations and their effective population size. The average MAF for both breeds was 0.27 ± 0.14 with a higher SNP proportion having high MAF values (≥ 0.3). The LD levels and distribution of allele frequencies found in this study suggest that it is possible to have adequate coverage throughout the genome of these breeds using the BovineSNP50K_v2, capturing the effect of most QTL related with productive traits, and ensuring an adequate prediction capacity in genomic analysis.

Application of genomic selection to crop breeding

Application of genomic selection to crop breeding

Evaluation of the potential use of a meta-population for genomic selection in autochthonous beef cattle populations

This study investigated the potential application of genomic selection under a multi-breed scheme in the Spanish autochthonous beef cattle populations using a simulation study that replicates the structure of linkage disequilibrium obtained from a sample of 25 triplets of sire/dam/offspring per population and using the BovineHD Beadchip. Purebred and combined reference sets were used for the genomic evaluation and several scenarios of different genetic architecture of the trait were investigated. The single-breed evaluations yielded the highest within-breed accuracies. Across breed accuracies were found low but positive on average confirming the genetic connectedness between the populations. If the same genotyping effort is split in several populations, the accuracies were lower when compared with single-breed evaluation, but showed a small advantage over small-sized purebred reference sets over the accuracies of subsequent generations. Besides, the genetic architecture of the trait did not show any relevant effect on the accuracy with the exception of rare variants, which yielded slightly lower results and higher loss of predictive ability over the generations.

Breeding Differently—the Digital Revolution: High-Throughput Phenotyping and Genotyping

A conventional potato breeding strategy uses targeted outcrossing, followed by phenotypic recurrent selection over a series of generations to identify improved cultivars. This paper reviews recent research in Australia aimed at improving the efficiency of such breeding. To develop marker-assisted selection (MAS) for traits of interest, our initial targets were qualitative disease resistances for potato cyst nematode (Globodera rostochiensis Ro1), Potato virus Y and Potato virus X. We undertook a cost analysis comparison between MAS and conventional screening, confirming that MAS would be cost-effective within a breeding programme. Then, as the majority of target traits are quantitative in nature, we also looked at methods to address these traits, including progeny testing and a quantitative genetic analysis technique to develop estimated breeding values (EBVs). We found the markers were useful for detecting the disease resistance characters, while the EBVs improved the analysis of the complex traits. Using a combination of MAS, EBVs and conventional screening methods, we then designed a breeding scheme for rapid selection of cultivars with multiple desirable traits, reducing the breeding cycle from over 10 to 4 years. We then explored the factors that will affect the application of genomic selection in potato and investigated strategies to incorporate genomic selection in potato breeding, as we found that it would accelerate genetic gain as the breeding cycle can be reduced to 1 year. Improvements in computational power are also flowing on to research capabilities such as sequencing, high-throughput phenotyping and data analysis, which will accelerate germplasm improvement and breeding. High-throughput phenotyping facilities are being developed that include automated glasshouse systems equipped with imaging sensors and in-field high-throughput phenotyping systems with sensors mounted on ground- or aerial-based vehicles. Using these technological improvements in phenotypic and genotypic analysis will reduce the breeding cycle in a cost-effective manner and means that we can now breed differently.

Efficiency of genomic prediction of non-assessed single crosses.

An important application of genomic selection in plant breeding is predicting untested single crosses (SCs). Most investigations on the prediction efficiency were based on tested SCs using cross-validation. The main objective was to assess the prediction efficiency by correlating the predicted and true genotypic values of untested SCs (accuracy) and measuring the efficacy of identification of the best 300 untested SCs (coincidence) using simulated data. We assumed 10,000 SNPs, 400 QTLs, two groups of 70 selected DH lines, and 4900 SCs. The heritabilities for the assessed SCs were 30, 60, and 100%. The scenarios included three sampling processes of DH lines, two sampling processes of SCs for testing, two SNP densities, DH lines from distinct and the same populations, DH lines from populations with lower LD, two genetic models, three statistical models, and three statistical approaches. We derived a model for genomic prediction based on SNP average effects of substitution and dominance deviations. The prediction accuracy is not affected by the linkage phase. The prediction of untested SCs is very efficient. The accuracies and coincidences ranged from ~0.8 and 0.5 at low heritability to 0.9 and 0.7 at high heritability, respectively. We also highlight the relevance of the overall LD and demonstrate that efficient prediction of untested SCs can be achieved for crops that show no heterotic pattern, for reduced training set size (10%), for SNP density of 1 cM, and for distinct sampling processes of DH lines based on random choice of the SCs for testing.

Application of genomic selection in farm animal breeding.

Genomic selection (GS) has become a widely accepted method in animal breeding to genetically improve economic traits. With the declining costs of high-density SNP chips and next-generation sequencing, GS has been applied in dairy cattle, swine, poultry and other animals and gained varying degrees of success. Currently, major challenges in GS studies include further reducing the cost of genome-wide SNP genotyping and improving the predictive accuracy of genomic estimated breeding value (GEBV). In this review, we summarize various methods for genome-wide SNP genotyping and GEBV prediction, and give a brief introduction of GS in livestock and poultry breeding. This review will provide a reference for further implementation of GS in farm animal breeding.

Effect of Trait Heritability, Training Population Size and Marker Density on Genomic Prediction Accuracy Estimation in 22 bi-parental Tropical Maize Populations.

Genomic selection is being used increasingly in plant breeding to accelerate genetic gain per unit time. One of the most important applications of genomic selection in maize breeding is to predict and select the best un-phenotyped lines in bi-parental populations based on genomic estimated breeding values. In the present study, 22 bi-parental tropical maize populations genotyped with low density SNPs were used to evaluate the genomic prediction accuracy (rMG) of the six trait-environment combinations under various levels of training population size (TPS) and marker density (MD), and assess the effect of trait heritability (h2), TPS and MD on rMG estimation. Our results showed that: (1) moderate rMG values were obtained for different trait-environment combinations, when 50% of the total genotypes was used as training population and ~200 SNPs were used for prediction; (2) rMG increased with an increase in h2, TPS and MD, both correlation and variance analyses showed that h2 is the most important factor and MD is the least important factor on rMG estimation for most of the trait-environment combinations; (3) predictions between pairwise half-sib populations showed that the rMG values for all the six trait-environment combinations were centered around zero, 49% predictions had rMG values above zero; (4) the trend observed in rMG differed with the trend observed in rMG/h, and h is the square root of heritability of the predicted trait, it indicated that both rMG and rMG/h values should be presented in GS study to show the accuracy of genomic selection and the relative accuracy of genomic selection compared with phenotypic selection, respectively. This study provides useful information to maize breeders to design genomic selection workflow in their breeding programs.

Genome-Wide Association and Prediction of Grain and Semolina Quality Traits in Durum Wheat Breeding Populations.

Grain yield and semolina quality traits are essential selection criteria in durum wheat breeding. However, high phenotypic screening costs limit selection to relatively few breeding lines in late generations. This selection paradigm confers relatively low selection efficiency due to the advancement of undesirable lines into expensive yield trials for grain yield and quality trait testing. Marker-aided selection can enhance selection efficiency, especially for traits that are difficult or costly to phenotype. The aim of this study was to identify major quality trait quantitative trait loci (QTL) for marker-assisted selection (MAS) and to explore potential application of genomic selection (GS) in a durum wheat breeding program. In this study, genome-wide association mapping was conducted for five quality traits using 1184 lines from the North Dakota State University (NDSU) durum wheat breeding program. Several QTL associated with test weight, semolina color, and gluten strength were identified. Genomic selection models were developed and forward prediction accuracies of 0.27 to 0.66 were obtained for the five quality traits. Our results show the potential for grain and semolina quality traits to be selected more efficiently through MAS and GS with further refinement. Considerable opportunity exists to extend these techniques to other traits such as grain yield and agronomic characteristics, further improving breeding efficiency in durum cultivar development.

Practical application of genomic selection in a doubled-haploid winter wheat breeding program

Crop improvement is a long-term, expensive institutional endeavor. Genomic selection (GS), which uses single nucleotide polymorphism (SNP) information to estimate genomic breeding values, has proven efficient to increasing genetic gain by accelerating the breeding process in animal breeding programs. As for crop improvement, with few exceptions, GS applicability remains in the evaluation of algorithm performance. In this study, we examined factors related to GS applicability in line development stage for grain yield using a hard red winter wheat (Triticum aestivum L.) doubled-haploid population. The performance of GS was evaluated in two consecutive years to predict grain yield. In general, the semi-parametric reproducing kernel Hilbert space prediction algorithm outperformed parametric genomic best linear unbiased prediction. For both parametric and semi-parametric algorithms, an upward bias in predictability was apparent in within-year cross-validation, suggesting the prerequisite of cross-year validation for a more reliable prediction. Adjusting the training population’s phenotype for genotype by environment effect had a positive impact on GS model’s predictive ability. Possibly due to marker redundancy, a selected subset of SNPs at an absolute pairwise correlation coefficient threshold value of 0.4 produced comparable results and reduced the computational burden of considering the full SNP set. Finally, in the context of an ongoing breeding and selection effort, the present study has provided a measure of confidence based on the deviation of line selection from GS results, supporting the implementation of GS in wheat variety development.

Accelerating wheat breeding for end-use quality with multi-trait genomic predictions incorporating near infrared and nuclear magnetic resonance-derived phenotypes

Using NIR and NMR predictions of quality traits overcomes a major barrier for the application of genomic selection to accelerate improvement in grain end-use quality traits of wheat. Grain end-use quality traits are among the most important in wheat breeding. These traits are difficult to breed for, as their assays require flour quantities only obtainable late in the breeding cycle, and are expensive. These traits are therefore an ideal target for genomic selection. However, large reference populations are required for accurate genomic predictions, which are challenging to assemble for these traits for the same reasons they are challenging to breed for. Here, we use predictions of end-use quality derived from near infrared (NIR) or nuclear magnetic resonance (NMR), that require very small amounts of flour, as well as end-use quality measured by industry standard assay in a subset of accessions, in a multi-trait approach for genomic prediction. The NIR and NMR predictions were derived for 19 end-use quality traits in 398 accessions, and were then assayed in 2420 diverse wheat accessions. The accessions were grown out in multiple locations and multiple years, and were genotyped for 51208 SNP. Incorporating NIR and NMR phenotypes in the multi-trait approach increased the accuracy of genomic prediction for most quality traits. The accuracy ranged from 0 to 0.47 before the addition of the NIR/NMR data, while after these data were added, it ranged from 0 to 0.69. Genomic predictions were reasonably robust across locations and years for most traits. Using NIR and NMR predictions of quality traits overcomes a major barrier for the application of genomic selection for grain end-use quality traits in wheat breeding.

Multiple-trait QTL mapping and genomic prediction for wool traits in sheep

BackgroundThe application of genomic selection to sheep breeding could lead to substantial increases in profitability of wool production due to the availability of accurate breeding values from single nucleotide polymorphism (SNP) data. Several key traits determine the value of wool and influence a sheep’s susceptibility to fleece rot and fly strike. Our aim was to predict genomic estimated breeding values (GEBV) and to compare three methods of combining information across traits to map polymorphisms that affect these traits.MethodsGEBV for 5726 Merino and Merino crossbred sheep were calculated using BayesR and genomic best linear unbiased prediction (GBLUP) with real and imputed 510,174 SNPs for 22 traits (at yearling and adult ages) including wool production and quality, and breech conformation traits that are associated with susceptibility to fly strike. Accuracies of these GEBV were assessed using fivefold cross-validation. We also devised and compared three approximate multi-trait analyses to map pleiotropic quantitative trait loci (QTL): a multi-trait genome-wide association study and two multi-trait methods that use the output from BayesR analyses. One BayesR method used local GEBV for each trait, while the other used the posterior probabilities that a SNP had an effect on each trait.ResultsBayesR and GBLUP resulted in similar average GEBV accuracies across traits (~0.22). BayesR accuracies were highest for wool yield and fibre diameter (>0.40) and lowest for skin quality and dag score (<0.10). Generally, accuracy was higher for traits with larger reference populations and higher heritability. In total, the three multi-trait analyses identified 206 putative QTL, of which 20 were common to the three analyses. The two BayesR multi-trait approaches mapped QTL in a more defined manner than the multi-trait GWAS. We identified genes with known effects on hair growth (i.e. FGF5, STAT3, KRT86, and ALX4) near SNPs with pleiotropic effects on wool traits.ConclusionsThe mean accuracy of genomic prediction across wool traits was around 0.22. The three multi-trait analyses identified 206 putative QTL across the ovine genome. Detailed phenotypic information helped to identify likely candidate genes.

Genome-enabled predictions for fruit weight and quality from repeated records in European peach progenies

BackgroundHighly polygenic traits such as fruit weight, sugar content and acidity strongly influence the agroeconomic value of peach varieties. Genomic Selection (GS) can accelerate peach yield and quality gain if predictions show higher levels of accuracy compared to phenotypic selection. The available IPSC 9K SNP array V1 allows standardized and highly reliable genotyping, preparing the ground for GS in peach.ResultsA repeatability model (multiple records per individual plant) for genome-enabled predictions in eleven European peach populations is presented. The analysis included 1147 individuals derived from both commercial and non-commercial peach or peach-related accessions. Considered traits were average fruit weight (FW), sugar content (SC) and titratable acidity (TA). Plants were genotyped with the 9K IPSC array, grown in three countries (France, Italy, Spain) and phenotyped for 3–5 years. An analysis of imputation accuracy of missing genotypic data was conducted using the software Beagle, showing that two of the eleven populations were highly sensitive to increasing levels of missing data. The regression model produced, for each trait and each population, estimates of heritability (FW:0.35, SC:0.48, TA:0.53, on average) and repeatability (FW:0.56, SC:0.63, TA:0.62, on average). Predictive ability was estimated in a five-fold cross validation scheme within population as the correlation of true and predicted phenotypes. Results differed by populations and traits, but predictive abilities were in general high (FW:0.60, SC:0.72, TA:0.65, on average).ConclusionsThis study assessed the feasibility of Genomic Selection in peach for highly polygenic traits linked to yield and fruit quality. The accuracy of imputing missing genotypes was as high as 96%, and the genomic predictive ability was on average 0.65, but could be as high as 0.84 for fruit weight or 0.83 for titratable acidity. The estimated repeatability may prove very useful in the management of the typical long cycles involved in peach productions. All together, these results are very promising for the application of genomic selection to peach breeding programmes.

Applications of Genomic Selection in Breeding Wheat for Rust Resistance.

There are a lot of methods developed to predict untested phenotypes in schemes commonly used in genomic selection (GS) breeding. The use of GS for predicting disease resistance has its own particularities: (a) most populations shows additivity in quantitative adult plant resistance (APR); (b) resistance needs effective combinations of major and minor genes; and (c) phenotype is commonly expressed in ordinal categorical traits, whereas most parametric applications assume that the response variable is continuous and normally distributed. Machine learning methods (MLM) can take advantage of examples (data) that capture characteristics of interest from an unknown underlying probability distribution (i.e., data-driven). We introduce some state-of-the-art MLM capable to predict rust resistance in wheat. We also present two parametric R packages for the reader to be able to compare.

The optimal number of partial least squares components in genomic selection for pork pH

ABSTRACT: The main application of genomic selection (GS) is the early identification of genetically superior animals for traits difficult-to-measure or lately evaluated, such as meat pH (measured after slaughter). Because the number of markers in GS is generally larger than the number of genotyped animals and these markers are highly correlated owing to linkage disequilibrium, statistical methods based on dimensionality reduction have been proposed. Among them, the partial least squares (PLS) technique stands out, because of its simplicity and high predictive accuracy. However, choosing the optimal number of components remains a relevant issue for PLS applications. Thus, we applied PLS (and principal component and traditional multiple regression) techniques to GS for pork pH traits (with pH measured at 45min and 24h after slaughter) and also identified the optimal number of PLS components based on the degree-of-freedom (DoF) and cross-validation (CV) methods. The PLS method out performs the principal component and traditional multiple regression techniques, enabling satisfactory predictions for pork pH traits using only genotypic data (low-density SNP panel). Furthermore, the SNP marker estimates from PLS revealed a relevant region on chromosome 4, which may affect these traits. The DoF and CV methods showed similar results for determining the optimal number of components in PLS analysis; thus, from the statistical viewpoint, the DoF method should be preferred because of its theoretical background (based on the "statistical information theory"), whereas CV is an empirical method based on computational effort.