Improving selection decisions with mating information by accounting for Mendelian sampling variances looking two generations ahead.

BACKGROUND: Breeding programs are judged by the genetic level of animals that are used to disseminate genetic progress. These animals are typically the best ones of the population. To maximise the genetic level of very good animals in the next generation, parents that are more likely to produce top performing offspring need to be selected. The ability of individuals to produce high-performing progeny differs because of differences in their breeding values and gametic variances. Differences in gametic variances among individuals are caused by differences in heterozygosity and linkage. The use of the gametic Mendelian sampling variance has been proposed before, for use in the usefulness criterion or Index5, and in this work, we extend existing approaches by not only considering the gametic Mendelian sampling variance of individuals, but also of their potential offspring. Thus, the criteria developed in this study plan one additional generation ahead. For simplicity, we assumed that the true quantitative trait loci (QTL) effects, genetic map and the haplotypes of all animals are known. RESULTS: In this study, we propose a new selection criterion, ExpBVSelGrOff, which describes the genetic level of selected grand-offspring that are produced by selected offspring of a particular mating. We compare our criterion with other published criteria in a stochastic simulation of an ongoing breeding program for 21 generations for proof of concept. ExpBVSelGrOff performed better than all other tested criteria, like the usefulness criterion or Index5 which have been proposed in the literature, without compromising short-term gains. After only five generations, when selection is strong (1%), selection based on ExpBVSelGrOff achieved 5.8% more commercial genetic gain and retained 25% more genetic variance without compromising inbreeding rate compared to selection based only on breeding values. CONCLUSIONS: Our proposed selection criterion offers a new tool to accelerate genetic progress for contemporary genomic breeding programs. It retains more genetic variance than previously published criteria that plan less far ahead. Considering future gametic Mendelian sampling variances in the selection process also seems promising for maintaining more genetic variance.

Genetic parameters, reciprocal cross differences, and age-related heterosis of egg-laying performance in chickens.

BACKGROUND: Egg-laying performance is economically important in poultry breeding programs. Crossbreeding between indigenous and elite commercial lines to exploit heterosis has been an upward trend in traditional layer breeding for niche markets. The objective of this study was to analyse the genetic background and to estimate the heterosis of longitudinal egg-laying traits in reciprocal crosses between an indigenous Beijing-You and an elite commercial White Leghorn layer line. Egg weights were measured for the first three eggs, monthly from 28 to 76 weeks of age, and at 86 and 100 weeks of age. Egg quality traits were measured at 32, 54, 72, 86, and 100 weeks of age. Egg production traits were measured from the start of lay until 43, 72, and 100 weeks of age. Heritabilities and phenotypic and genetic correlations were estimated. Heterosis was estimated as the percentage difference of performance of a crossbred from that of the parental average. Reciprocal cross differences were estimated as the difference between the reciprocal crossbreds as a percentage of the parental average. RESULTS: Estimates of heritability of egg weights ranged from 0.29 to 0.75. Estimates of genetic correlations between egg weights at different ages ranged from 0.72 to 1.00. Estimates of heritability for cumulative egg numbers until 43, 72, and 100 weeks of age were around 0.15. Estimates of heterosis for egg weight and cumulative egg number increased with age, ranging from 1.0 to 9.0% and from 1.4 to 11.6%, respectively. From 72 to 100 weeks of age, crossbreds produced more eggs per week than the superior parent White Leghorn (3.5 eggs for White Leghorn, 3.8 and 3.9 eggs for crossbreds). Heterosis for eggshell thickness ranged from 2.7 to 6.6% when using Beijing-You as the sire breed. No significant difference between reciprocal crosses was observed for the investigated traits, except for eggshell strength at 54 weeks of age. CONCLUSIONS: The heterosis was substantial for egg weight and cumulative egg number, and increased with age, suggesting that non-additive genetic effects are important in crossbreds between the indigenous and elite breeds. Generally, the crossbreds performed similar to or even outperformed the commercial White Leghorns for egg production persistency.

The long-term effects of genomic selection: 2. Changes in allele frequencies of causal loci and new mutations.

Genetic selection has been applied for many generations in animal, plant, and experimental populations. Selection changes the allelic architecture of traits to create genetic gain. It remains unknown whether the changes in allelic architecture are different for the recently introduced technique of genomic selection compared to traditional selection methods and whether they depend on the genetic architectures of traits. Here, we investigate the allele frequency changes of old and new causal loci under 50 generations of phenotypic, pedigree, and genomic selection, for a trait controlled by either additive, additive and dominance, or additive, dominance, and epistatic effects. Genomic selection resulted in slightly larger and faster changes in allele frequencies of causal loci than pedigree selection. For each locus, allele frequency change per generation was not only influenced by its statistical additive effect but also to a large extent by the linkage phase with other loci and its allele frequency. Selection fixed a large number of loci, and 5 times more unfavorable alleles became fixed with genomic and pedigree selection than with phenotypic selection. For pedigree selection, this was mainly a result of increased genetic drift, while genetic hitchhiking had a larger effect on genomic selection. When epistasis was present, the average allele frequency change was smaller (∼15% lower), and a lower number of loci became fixed for all selection methods. We conclude that for long-term genetic improvement using genomic selection, it is important to consider hitchhiking and to limit the loss of favorable alleles.

Integration of beef cattle international pedigree and genomic estimated breeding values into national evaluations, with an application to the Italian Limousin population.

BACKGROUND: International evaluations combine data from different countries allowing breeders to have access to larger panels of elite bulls and to increase the accuracy of estimated breeding values (EBV). However, international and national evaluations can use different sources of information to compute EBV (EBVINT and EBVNAT, respectively), leading to differences between them. Choosing one of these EBV results in losing the information that is contained only in the discarded EBV. Our objectives were to define and validate a procedure to integrate publishable sires' EBVINT and their associated reliabilities computed from pedigree-based or single-step international beef cattle evaluations into national evaluations to obtain "blended" EBV. The Italian (ITA) pedigree-based national evaluation was used as a case study to validate the integration procedure. METHODS: Publishable sires' international information, i.e. EBVINT and their associated reliabilities, was included in the national evaluation as pseudo-records. Data were available for 444,199 individual age-adjusted weaning weights of Limousin cattle from eight countries and 17,607 genotypes from four countries (ITA excluded). To mimic differences between international and national evaluations, international evaluations included phenotypes (and genotypes) of animals born prior to January 2019, while national evaluations included ITA phenotypes of animals born until April 2019. International evaluations using all available information were considered as reference scenarios. Publishable sires were divided into three groups: sires with ≥ 15, < 15 and no recorded offspring in ITA. RESULTS: Overall, for these three groups, integrating either pedigree-based or single-step international information into national pedigree-based evaluations improved the similarity of the blended EBV with the reference EBV compared to national evaluations without integration. For instance, the correlation with the reference EBV for direct (maternal) EBV went from 0.61 (0.79) for a national evaluation without integration to 0.97 (0.88) when integrating single-step international information, on average across all groups of publishable sires. CONCLUSIONS: Our proposed one-animal-at-a-time integration procedure yields blended EBV that are in close agreement with full international EBV for all groups of animals analysed. The procedure can be directly applied by countries since it does not rely on specific software and is computationally inexpensive, allowing straightforward integration of publishable sires' EBVINT from pedigree-based or single-step based international beef cattle evaluations into national evaluations.

Efficient large-scale single-step evaluations and indirect genomic prediction of genotyped selection candidates.

BACKGROUND: Single-step genomic best linear unbiased prediction (ssGBLUP) models allow the combination of genomic, pedigree, and phenotypic data into a single model, which is computationally challenging for large genotyped populations. In practice, genotypes of animals without their own phenotype and progeny, so-called genotyped selection candidates, can become available after genomic breeding values have been estimated by ssGBLUP. In some breeding programmes, genomic estimated breeding values (GEBV) for these animals should be known shortly after obtaining genotype information but recomputing GEBV using the full ssGBLUP takes too much time. In this study, first we compare two equivalent formulations of ssGBLUP models, i.e. one that is based on the Woodbury matrix identity applied to the inverse of the genomic relationship matrix, and one that is based on marker equations. Second, we present computationally-fast approaches to indirectly compute GEBV for genotyped selection candidates, without the need to do the full ssGBLUP evaluation. RESULTS: The indirect approaches use information from the latest ssGBLUP evaluation and rely on the decomposition of GEBV into its components. The two equivalent ssGBLUP models and indirect approaches were tested on a six-trait calving difficulty model using Irish dairy and beef cattle data that include 2.6 million genotyped animals of which about 500,000 were considered as genotyped selection candidates. When using the same computational approaches, the solving phase of the two equivalent ssGBLUP models showed similar requirements for memory and time per iteration. The computational differences between them were due to the preprocessing phase of the genomic information. Regarding the indirect approaches, compared to GEBV obtained from single-step evaluations including all genotypes, indirect GEBV had correlations higher than 0.99 for all traits while showing little dispersion and level bias. CONCLUSIONS: In conclusion, ssGBLUP predictions for the genotyped selection candidates were accurately approximated using the presented indirect approaches, which are more memory efficient and computationally fast, compared to solving a full ssGBLUP evaluation. Thus, indirect approaches can be used even on a weekly basis to estimate GEBV for newly genotyped animals, while the full single-step evaluation is done only a few times within a year.

The use of a genomic relationship matrix for breed assignment of cattle breeds: comparison and combination with a machine learning method.

To develop a breed assignment model, three main steps are generally followed: 1) The selection of breed informative single nucleotide polymorphism (SNP); 2) The training of a model, based on a reference population, that allows to classify animals to their breed of origin; and 3) The validation of the developed model on external animals i.e., that were not used in previous steps. However, there is no consensus in the literature about which methodology to follow for the first step, nor about the number of SNP to be selected. This can raise many questions when developing the model and lead to the use of sophisticated methodologies for selecting SNP (e.g., with iterative algorithms, partitions of SNP, or combination of several methods). Therefore, it may be of interest to avoid the first step by the use of all the available SNP. For this purpose, we propose the use of a genomic relationship matrix (GRM), combined or not with a machine learning method, for breed assignment. We compared it with a previously developed model based on selected informative SNP. Four methodologies were investigated: 1) The PLS_NSC methodology: selection of SNP based on a partial least square-discriminant analysis (PLS-DA) and breed assignment by classification based on the nearest shrunken centroids (NSC) method; 2) Breed assignment based on the highest mean relatedness of an animal to the reference populations of each breed (referred to mean_GRM); 3) Breed assignment based on the highest SD of the relatedness of an animal to the reference populations of each breed (referred to SD_GRM) and 4) The GRM_SVM methodology: the use of means and SD of the relatedness defined in mean_GRM and SD_GRM methodologies combined with the linear support vector machine (SVM), a machine learning method used for classification. Regarding mean global accuracies, results showed that the use of mean_GRM or GRM_SVM was not significantly different (Bonferroni corrected P > 0.0083) than the model based on a reduced SNP panel (PLS_NSC). Moreover, the mean_GRM and GRM_SVM methodology were more efficient than PLS_NSC as it was faster to compute. Therefore, it is possible to bypass the selection of SNP and, by the use of a GRM, to develop an efficient breed assignment model. In routine, we recommend the use of GRM_SVM over mean_GRM as it gave a slightly increased global accuracy, which can help endangered breeds to be maintained. The script to execute the different methodologies can be accessed on: https://github.com/hwilmot675/Breed_assignment.

Breed assignment models generally rely on three main steps: 1) Selection of markers that allow to distinguish the breeds under study; 2) Development of a classification model that assigns each animal to its breed of origin; and 3) Validation of the developed model with new animals, to verify that the developed model is not overfitted. The first step often raises several questions about the methodology to select the best markers or about the number of markers to select. That is why it can be interesting to avoid this first step and to use an appropriate methodology that performs similarly without the need for single nucleotide polymorphism (SNP) selection. In this study, we developed different methodologies based on the genomic relationship matrix (GRM), combined or not with a machine learning method, to assign animals to their breed of origin. The results showed that the model based on a GRM combined with a machine learning method showed equivalent percentage of correct assignment to a previously developed model relying on SNP selection while being substantially faster to compute. It is therefore possible to assign animals to their breed by the use of a GRM and to bypass the first step of selection of SNP.

Predicting the impact of genotype-by-genotype interaction on the purebred-crossbred genetic correlation from phenotype and genotype marker data of parental lines.

BACKGROUND: The genetic correlation between purebred (PB) and crossbred (CB) performances ([Formula: see text]) partially determines the response in CB when selection is on PB performance in the parental lines. An earlier study has derived expressions for an upper and lower bound of [Formula: see text], using the variance components of the parental purebred lines, including e.g. the additive genetic variance in the sire line for the trait expressed in one of the dam lines. How to estimate these variance components is not obvious, because animals from one parental line do not have phenotypes for the trait expressed in the other line. Thus, the aim of this study was to propose and compare three methods for approximating the required variance components. The first two methods are based on (co)variances of genomic estimated breeding values (GEBV) in the line of interest, either accounting for shrinkage (VCGEBV-S) or not (VCGEBV). The third method uses restricted maximum likelihood (REML) estimates directly from univariate and bivariate analyses (VCREML) by ignoring that the variance components should refer to the line of interest, rather than to the line in which the trait is expressed. We validated these methods by comparing the resulting predicted bounds of [Formula: see text] with the [Formula: see text] estimated from PB and CB data for five traits in a three-way cross in pigs. RESULTS: With both VCGEBV and VCREML, the estimated [Formula: see text] (plus or minus one standard error) was between the upper and lower bounds in 14 out of 15 cases. However, the range between the bounds was much smaller with VCREML (0.15-0.22) than with VCGEBV (0.44-0.57). With VCGEBV-S, the estimated [Formula: see text] was between the upper and lower bounds in only six out of 15 cases, with the bounds ranging from 0.21 to 0.44. CONCLUSIONS: We conclude that using REML estimates of variance components within and between parental lines to predict the bounds of [Formula: see text] resulted in better predictions than methods based on GEBV. Thus, we recommend that the studies that estimate [Formula: see text] with genotype data also report estimated genetic variance components within and between the parental lines.

Impact of genomic preselection on subsequent ssGBLUP evaluation of preselected animals for scarcely recorded feed intake in pigs.

We have previously shown that single-step genomic best linear unbiased prediction (ssGBLUP) estimates breeding values of genomically preselected animals without preselection bias for widely recorded traits, that is traits recorded for the majority of animals in the breeding population. This study investigated the impact of genomic preselection (GPS) on accuracy and bias in ssGBLUP evaluation of genomically preselected animals for a scarcely recorded trait, that is a trait recorded for only a small proportion of the animals, which generally has a lower prediction accuracy than widely recorded traits, mainly due to having a much smaller number of phenotypes available. We used data from a commercial pig breeding program, considering feed intake as a scarcely recorded target trait, being available for ~30% of the animals with phenotypes for any trait, and average daily gain, backfat thickness and loin depth as widely recorded predictor traits, being available for >95% of the animals with phenotypes for any trait. The data contained the routine GPS implemented by commercial animal breeding programs, and we retrospectively implemented two scenarios with additional layers of GPS by discarding pedigree, genotypes and phenotypes of animals without progeny. The ssGBLUP evaluation following GPS used records only from the target trait, only from the predictor traits, or both. Accuracy for feed intake did not differ statistically across GPS scenarios, although it tended to decrease with more intense GPS. The accuracy had average values of 0.37, 0.44, and 0.45 across all GPS scenarios when, respectively, records from only the target trait, only the predictor traits, or both were used in the ssGBLUP evaluation. Considerable deflation of the genomic breeding values for feed intake was observed in the most stringent GPS scenario, due to the variance components being underestimated as a result of the limited amount of strongly preselected data. As long as (co)variance components were unbiased, no or only marginal bias was observed. These results for accuracy and bias were observed whether records of the scarcely recorded target trait, of the predictor traits, or both were used in the ssGBLUP evaluation. Our results show that for the scarcely recorded feed intake in pigs, ssGBLUP is able to estimate breeding values of preselected animals without preselection bias, similarly as previously observed for widely recorded traits.

International single-step SNPBLUP beef cattle evaluations for Limousin weaning weight.

BACKGROUND: Compared to national evaluations, international collaboration projects further improve accuracies of estimated breeding values (EBV) by building larger reference populations or performing a joint evaluation using data (or proxy of them) from different countries. Genomic selection is increasingly adopted in beef cattle, but, to date, the benefits of including genomic information in international evaluations have not been explored. Our objective was to develop an international beef cattle single-step genomic evaluation and investigate its impact on the accuracy and bias of genomic evaluations compared to current pedigree-based evaluations. METHODS: Weaning weight records were available for 331,593 animals from seven European countries. The pedigree included 519,740 animals. After imputation and quality control, 17,607 genotypes at a density of 57,899 single nucleotide polymorphisms (SNPs) from four countries were available. We implemented two international scenarios where countries were modelled as different correlated traits: an international genomic single-step SNP best linear unbiased prediction (SNPBLUP) evaluation (ssSNPBLUPINT) and an international pedigree-based BLUP evaluation (PBLUPINT). Two national scenarios were implemented for pedigree and genomic evaluations using only nationally submitted phenotypes and genotypes. Accuracies, level and dispersion bias of EBV of animals born from 2014 onwards, and increases in population accuracies were estimated using the linear regression method. RESULTS: On average across countries, 39 and 17% of sires and maternal-grand-sires with recorded (grand-)offspring across two countries were genotyped. ssSNPBLUPINT showed the highest accuracies of EBV and, compared to PBLUPINT, led to increases in population accuracy of 13.7% for direct EBV, and 25.8% for maternal EBV, on average across countries. Increases in population accuracies when moving from national scenarios to ssSNPBLUPINT were observed for all countries. Overall, ssSNPBLUPINT level and dispersion bias remained similar or slightly reduced compared to PBLUPINT and national scenarios. CONCLUSIONS: International single-step SNPBLUP evaluations are feasible and lead to higher population accuracies for both large and small countries compared to current international pedigree-based evaluations and national evaluations. These results are likely related to the larger multi-country reference population and the inclusion of phenotypes from relatives recorded in other countries via single-step international evaluations. The proposed international single-step approach can be applied to other traits and breeds.

Adding gene transcripts into genomic prediction improves accuracy and reveals sampling time dependence.

Recent developments allowed generating multiple high-quality 'omics' data that could increase the predictive performance of genomic prediction for phenotypes and genetic merit in animals and plants. Here, we have assessed the performance of parametric and nonparametric models that leverage transcriptomics in genomic prediction for 13 complex traits recorded in 478 animals from an outbred mouse population. Parametric models were implemented using the best linear unbiased prediction, while nonparametric models were implemented using the gradient boosting machine algorithm. We also propose a new model named GTCBLUP that aims to remove between-omics-layer covariance from predictors, whereas its counterpart GTBLUP does not do that. While gradient boosting machine models captured more phenotypic variation, their predictive performance did not exceed the best linear unbiased prediction models for most traits. Models leveraging gene transcripts captured higher proportions of the phenotypic variance for almost all traits when these were measured closer to the moment of measuring gene transcripts in the liver. In most cases, the combination of layers was not able to outperform the best single-omics models to predict phenotypes. Using only gene transcripts, the gradient boosting machine model was able to outperform best linear unbiased prediction for most traits except body weight, but the same pattern was not observed when using both single nucleotide polymorphism genotypes and gene transcripts. Although the GTCBLUP model was not able to produce the most accurate phenotypic predictions, it showed the highest accuracies for breeding values for 9 out of 13 traits. We recommend using the GTBLUP model for prediction of phenotypes and using the GTCBLUP for prediction of breeding values.

Estimation of dam line composition of 3-way crossbred animals using genomic information.

BACKGROUND: In genomic prediction including data of 3- or 4-way crossbred animals, line composition is usually fitted as a regression on expected line proportions, which are 0.5, 0.25 and 0.25, respectively, for 3-way crossbred animals. However, actual line proportions for the dam lines can vary between ~ 0.1 and 0.4, and ignoring this variation may affect the genomic estimated breeding values of purebred selection candidates. Our aim was to validate a proposed gold standard to evaluate different approaches for estimating line proportions using simulated data, and to subsequently use this in actual 3-way crossbred broiler data to evaluate several other methods. RESULTS: Analysis of simulated data confirmed that line proportions computed from assigned breed-origin-of-alleles (BOA) provide a very accurate gold standard, even if the parental lines are closely related. Alternative investigated methods were linear regression of genotypes on line-specific allele frequencies, maximum likelihood estimation using the program ADMIXTURE, and the genomic relationship of crossbred animals with their maternal grandparents. The results from the simulated data showed that the genomic relationship with the maternal grandparent was most accurate, and least affected by closer relationships between the dam lines. Linear regression and ADMIXTURE performed similarly for unrelated lines, but their accuracy dropped considerably when the dam lines were more closely related. In almost all cases, estimates improved after adjusting them to ensure that the sum of dam line contributions within animals was equal to 0.5, and within dam line and across animals the average was equal to 0.25. Results from the broiler data were much more similar between methods. In both cases, stringent linkage disequilibrium pruning of genotype data led to a relatively low accuracy of predicted line proportions, due to the loss of too many single nucleotide polymorphisms. CONCLUSIONS: With relatively unrelated parental lines as typical in crosses in pigs and poultry, linear regression of crossbred genotypes on line-specific allele frequencies and ADMIXTURE are very competitive methods. Thus, linear regression may be the method of choice, as it does not require genotypes of grandparents, is computationally very efficient, and easily implemented and adapted for considering the specific nature of the crossbred animals analysed.

Impact of genomic preselection on subsequent genetic evaluations with ssGBLUP using real data from pigs.

BACKGROUND: Empirically assessing the impact of preselection on genetic evaluation of preselected animals requires comparing scenarios that take different approaches into account, including scenarios without preselection. However, preselection is almost always performed in animal breeding programs, so it is difficult to have a dataset without preselection. Hence, most studies on preselection have used simulated datasets, and have concluded that genomic estimated breeding values (GEBV) from subsequent single-step genomic best linear unbiased prediction (ssGBLUP) evaluations are unbiased. The aim of this study was to investigate the impact of genomic preselection (GPS) on accuracy and bias in subsequent ssGBLUP evaluations, using data from a commercial pig breeding program. METHODS: We used data on average daily gain during performance testing, average daily gain throughout life, backfat thickness, and loin depth from one sire line and one dam line of pigs. As these traits have different weights in the breeding goals of the two lines, we analyzed the lines separately. For each line, we implemented a reference GPS scenario that kept all available data, against which the next two scenarios were compared. We then implemented two other scenarios with additional layers of GPS by removing all animals without progeny either (i) only in the validation generation, or (ii) in all generations. We conducted subsequent ssGBLUP evaluations for each GPS scenario, using all the data remaining after implementing the GPS scenario. Accuracy and bias were computed by comparing GEBV against progeny yield deviations of validation animals. RESULTS: Results for all traits and in both lines showed a marginal loss in accuracy due to the additional layers of GPS. Average accuracies across all GPS scenarios in the two lines were 0.39, 0.47, 0.56, and 0.60, for average daily gain during performance testing and throughout life, backfat thickness, and loin depth, respectively. Biases were largely absent, and when present, did not differ greatly between the GPS scenarios. CONCLUSIONS: We conclude that the impact of preselection on accuracy and bias in subsequent ssGBLUP evaluations of selection candidates in pigs is generally minimal. We expect this conclusion to apply for other animal breeding programs as well, since preselection of any type or intensity generally has the same effect in animal breeding programs.

The long-term effects of genomic selection: 1. Response to selection, additive genetic variance, and genetic architecture.

BACKGROUND: Genomic selection has revolutionized genetic improvement in animals and plants, but little is known about its long-term effects. Here, we investigated the long-term effects of genomic selection on response to selection, genetic variance, and the genetic architecture of traits using stochastic simulations. We defined the genetic architecture as the set of causal loci underlying each trait, their allele frequencies, and their statistical additive effects. We simulated a livestock population under 50 generations of phenotypic, pedigree, or genomic selection for a single trait, controlled by either only additive, additive and dominance, or additive, dominance, and epistatic effects. The simulated epistasis was based on yeast data. RESULTS: Short-term response was always greatest with genomic selection, while response after 50 generations was greater with phenotypic selection than with genomic selection when epistasis was present, and was always greater than with pedigree selection. This was mainly because loss of genetic variance and of segregating loci was much greater with genomic and pedigree selection than with phenotypic selection. Compared to pedigree selection, selection response was always greater with genomic selection. Pedigree and genomic selection lost a similar amount of genetic variance after 50 generations of selection, but genomic selection maintained more segregating loci, which on average had lower minor allele frequencies than with pedigree selection. Based on this result, genomic selection is expected to better maintain genetic gain after 50 generations than pedigree selection. The amount of change in the genetic architecture of traits was considerable across generations and was similar for genomic and pedigree selection, but slightly less for phenotypic selection. Presence of epistasis resulted in smaller changes in allele frequencies and less fixation of causal loci, but resulted in substantial changes in statistical additive effects across generations. CONCLUSIONS: Our results show that genomic selection outperforms pedigree selection in terms of long-term genetic gain, but results in a similar reduction of genetic variance. The genetic architecture of traits changed considerably across generations, especially under selection and when non-additive effects were present. In conclusion, non-additive effects had a substantial impact on the accuracy of selection and long-term response to selection, especially when selection was accurate.

Comparison of linkage disequilibrium estimated from genotypes versus haplotypes for crossbred populations.

BACKGROUND: Linkage disequilibrium (LD) is commonly measured based on the squared coefficient of correlation [Formula: see text] between the alleles at two loci that are carried by haplotypes. LD can also be estimated as the [Formula: see text] between unphased genotype dosage at two loci when the allele frequencies and inbreeding coefficients at both loci are identical for the parental lines. Here, we investigated whether [Formula: see text] for a crossbred population (F1) can be estimated using genotype data. The parental lines of the crossbred (F1) can be purebred or crossbred. METHODS: We approached this by first showing that inbreeding coefficients for an F1 crossbred population are negative, and typically differ in size between loci. Then, we proved that the expected [Formula: see text] computed from unphased genotype data is expected to be identical to the [Formula: see text] computed from haplotype data for an F1 crossbred population, regardless of the inbreeding coefficients at the two loci. Finally, we investigated the bias and precision of the [Formula: see text] estimated using unphased genotype versus haplotype data in stochastic simulation. RESULTS: Our findings show that estimates of [Formula: see text] based on haplotype and unphased genotype data are both unbiased for different combinations of allele frequencies, sample sizes (900, 1800, and 2700), and levels of LD. In general, for any allele frequency combination and [Formula: see text] value scenarios considered, and for both methods to estimate [Formula: see text], the precision of the estimates increased, and the bias of the estimates decreased as sample size increased, indicating that both estimators are consistent. For a given scenario, the [Formula: see text] estimates using haplotype data were more precise and less biased using haplotype data than using unphased genotype data. As sample size increased, the difference in precision and biasedness between the [Formula: see text] estimates using haplotype data and unphased genotype data decreased. CONCLUSIONS: Our theoretical derivations showed that estimates of LD between loci based on unphased genotypes and haplotypes in F1 crossbreds have identical expectations. Based on our simulation results, we conclude that the LD for an F1 crossbred population can be accurately estimated from unphased genotype data. The results also apply for other crosses (F2, F3, Fn, BC1, BC2, and BCn), as long as (selected) individuals from the two parental lines mate randomly.

Prediction performance of linear models and gradient boosting machine on complex phenotypes in outbred mice.

We compared the performance of linear (GBLUP, BayesB, and elastic net) methods to a nonparametric tree-based ensemble (gradient boosting machine) method for genomic prediction of complex traits in mice. The dataset used contained genotypes for 50,112 SNP markers and phenotypes for 835 animals from 6 generations. Traits analyzed were bone mineral density, body weight at 10, 15, and 20 weeks, fat percentage, circulating cholesterol, glucose, insulin, triglycerides, and urine creatinine. The youngest generation was used as a validation subset, and predictions were based on all older generations. Model performance was evaluated by comparing predictions for animals in the validation subset against their adjusted phenotypes. Linear models outperformed gradient boosting machine for 7 out of 10 traits. For bone mineral density, cholesterol, and glucose, the gradient boosting machine model showed better prediction accuracy and lower relative root mean squared error than the linear models. Interestingly, for these 3 traits, there is evidence of a relevant portion of phenotypic variance being explained by epistatic effects. Using a subset of top markers selected from a gradient boosting machine model helped for some of the traits to improve the accuracy of prediction when these were fitted into linear and gradient boosting machine models. Our results indicate that gradient boosting machine is more strongly affected by data size and decreased connectedness between reference and validation sets than the linear models. Although the linear models outperformed gradient boosting machine for the polygenic traits, our results suggest that gradient boosting machine is a competitive method to predict complex traits with assumed epistatic effects.

The impact of direct-maternal genetic correlations on international beef cattle evaluations for Limousin weaning weight.

In beef cattle maternally influenced traits, estimates of direct-maternal genetic correlations (rdm) are usually reported to be negative. In international evaluations, rdm can differ both within countries (rdm_WC) and between countries (rdm_BC). The rdm_BC are difficult to estimate and are assumed to be zero in the current model for international beef cattle evaluations (Interbeef). Our objective was to investigate re-ranking of international estimated breeding values (IEBVs) in international beef cattle evaluations between models that either used estimated values for rdm or assumed them to be 0. Age-adjusted weaning weights and pedigree data were available for Limousin beef cattle from ten European countries. International EBVs were obtained using a multi-trait animal model with countries modeled as different traits. We compared IEBVs from a model that uses estimated rdm_BC (ranging between -0.14 and +0.14) and rdm_WC (between -0.33 and +0.40) with IEBVs obtained either from the current model that assumes rdm_BC to be 0, or from an alternative model that assumes both rdm_BC and rdm_WC to be 0. Direct and maternal IEBVs were compared across those three scenarios for different groups of animals. The ratio of population accuracies from the linear regression method was used to further investigate the impact of rdm on international evaluations, for both the whole set of animals in the evaluation and the domestic ones. Ignoring rdm_BC, i.e., replacing estimated values with 0, resulted in no (rank correlations > 0.99) or limited (between 0.98 and 0.99) re-ranking for direct and maternal IEBVs, respectively. Both rdm_BC and rdm_WC had less impact on direct IEBVs than on maternal IEBVs. Re-ranking of maternal IEBVs decreased with increasing reliability. Ignoring rdm_BC resulted in no re-ranking for sires with IEBVs that might be exchanged across countries and limited re-ranking for the top 100 sires. Using estimated rdm_BC values instead of considering them to be 0 resulted in null to limited increases in population accuracy. Ignoring both rdm_BC and rdm_WC resulted in considerable re-ranking of animals' IEBVs in all groups of animals evaluated. This study showed the limited impact of the current practice of ignoring rdm_BC in international evaluations for Limousin weaning weight, most likely because the estimated rdm_BC was close to 0. We expect that these conclusions can be extended to other traits that have reported rdm values in the range of rdm_WC values for weaning weight in Limousin.

Review: optimizing genomic selection for crossbred performance by model improvement and data collection.

Breeding programs aiming to improve the performance of crossbreds may benefit from genomic prediction of crossbred (CB) performance for purebred (PB) selection candidates. In this review, we compared genomic prediction strategies that differed in 1) the genomic prediction model used or 2) the data used in the reference population. We found 27 unique studies, two of which used deterministic simulation, 11 used stochastic simulation, and 14 real data. Differences in accuracy and response to selection between strategies depended on i) the value of the purebred crossbred genetic correlation (rpc), ii) the genetic distance between the parental lines, iii) the size of PB and CB reference populations, and iv) the relatedness of these reference populations to the selection candidates. In studies where a PB reference population was used, the use of a dominance model yielded accuracies that were equal to or higher than those of additive models. When rpc was lower than ~0.8, and was caused mainly by G × E, it was beneficial to create a reference population of PB animals that are tested in a CB environment. In general, the benefit of collecting CB information increased with decreasing rpc. For a given rpc, the benefit of collecting CB information increased with increasing size of the reference populations. Collecting CB information was not beneficial when rpc was higher than ~0.9, especially when the reference populations were small. Collecting only phenotypes of CB animals may slightly improve accuracy and response to selection, but requires that the pedigree is known. It is, therefore, advisable to genotype these CB animals as well. Finally, considering the breed-origin of alleles allows for modeling breed-specific effects in the CB, but this did not always lead to higher accuracies. Our review shows that the differences in accuracy and response to selection between strategies depend on several factors. One of the most important factors is rpc, and we, therefore, recommend to obtain accurate estimates of rpc of all breeding goal traits. Furthermore, knowledge about the importance of components of rpc (i.e., dominance, epistasis, and G × E) can help breeders to decide which model to use, and whether to collect data on animals in a CB environment. Future research should focus on the development of a tool that predicts accuracy and response to selection from scenario specific parameters.

Convergence behavior of single-step GBLUP and SNPBLUP for different termination criteria.

BACKGROUND: The preconditioned conjugate gradient (PCG) method is the current method of choice for iterative solving of genetic evaluations. The relative difference between two successive iterates and the relative residual of the system of equations are usually chosen as a termination criterion for the PCG method in animal breeding. However, our initial analyses showed that these two commonly used termination criteria may report that a PCG method applied to a single-step single nucleotide polymorphism best linear unbiased prediction (ssSNPBLUP) is not converged yet, whereas the solutions are accurate enough for practical use. Therefore, the aim of this study was to propose two termination criteria that have been (partly) developed in other fields, but are new in animal breeding, and to compare their behavior to that of the two termination criteria widely used in animal breeding for the PCG method applied to ssSNPBLUP. The convergence patterns of ssSNPBLUP were also compared to the convergence patterns of single-step genomic BLUP (ssGBLUP). RESULTS: Building upon previous work, we propose two termination criteria that take the properties of the system of equations into account. These two termination criteria are directly related to the relative error of the iterates with respect to the true solutions. Based on pig and dairy cattle datasets, we show that the preconditioned coefficient matrices of ssSNPBLUP and ssGBLUP have similar properties when using a second-level preconditioner for ssSNPBLUP. Therefore, the PCG method applied to ssSNPBLUP and ssGBLUP converged similarly based on the relative error of the iterates with respect to the true solutions. This similar convergence behavior between ssSNPBLUP and ssGBLUP was observed for both proposed termination criteria. This was, however, not the case for the termination criterion defined as the relative residual when applied to the dairy cattle evaluations. CONCLUSION: Our results showed that the PCG method can converge similarly when applied to ssSNPBLUP and to ssGBLUP. The two proposed termination criteria always depicted these similar convergence behaviors, and we recommend them for comparing convergence properties of different models and for routine evaluations.

Predicting the purebred-crossbred genetic correlation from the genetic variance components in the parental lines.

BACKGROUND: The genetic correlation between purebred and crossbred performance ([Formula: see text]) is an important parameter in pig and poultry breeding, because response to selection in crossbred performance depends on the value of [Formula: see text] when selection is based on purebred (PB) performance. The value of [Formula: see text] can be substantially lower than 1, which is partly due to differences in allele frequencies between parental lines when non-additive genetic effects are present. This relationship between [Formula: see text] and parental allele frequencies suggests that [Formula: see text] can be expressed as a function of genetic parameters for the trait in the parental lines. In this study, we derived expressions for [Formula: see text] based on genetic variances within, and the genetic covariance between parental lines. It is important to note that the variance components used in our expressions are not the components that are typically estimated in empirical data. The expressions were derived for a genetic model with additive and dominance effects (D), and additive and epistatic additive-by-additive effects (EAA). We validated our expressions using simulations of purebred parental lines and their crosses, where the parental lines were either selected or not. Finally, using these simulations, we investigated the value of [Formula: see text] for genetic models with both dominance and epistasis or with other types of epistasis, for which expressions could not be derived. RESULTS: Our simulations show that when non-additive effects are present, [Formula: see text] decreases with increasing differences in allele frequencies between the parental lines. Genetic models that involve dominance result in lower values of [Formula: see text] than genetic models that involve epistasis only. Using information of parental lines only, our expressions provide exact estimates of [Formula: see text] for models D and EAA, and accurate upper and lower bounds of [Formula: see text] for two other genetic models. CONCLUSION: This work lays the foundation to enable estimation of [Formula: see text] from information collected in PB parental lines only.

Technical note: Genetic groups in single-step single nucleotide polymorphism best linear unbiased predictor.

Genetic groups, also called unknown or phantom parents groups, are often used in dairy cattle genetic evaluations to account for selection that cannot be accounted for by known genetic relationships. With the advent of genomic evaluations, the theory of genetic groups was extended to the so-called single-step genomic BLUP (ssGBLUP). In short, genetic groups can be fitted in ssGBLUP through regression effects, or by including them in the pedigree and computing the adequate combined pedigree and genomic relationship matrix. In this study, we applied the so-called Quaas and Pollak transformation to a system of equations for single-step SNP BLUP (ssSNPBLUP), such that genetic groups can thereafter be included in the pedigree. The example in this study showed that including the genetic groups in the pedigree for ssSNPBLUP allowed reduced memory burden and computational costs in comparison to genetic groups fitted as covariates.