Reinhard Reents

A common reference population from four European Holstein populations increases reliability of genomic predictions.

BackgroundSize of the reference population and reliability of phenotypes are crucial factors influencing the reliability of genomic predictions. It is therefore useful to combine closely related populations. Increased accuracies of genomic predictions depend on the number of individuals added to the reference population, the reliability of their phenotypes, and the relatedness of the populations that are combined.MethodsThis paper assesses the increase in reliability achieved when combining four Holstein reference populations of 4000 bulls each, from European breeding organizations, i.e. UNCEIA (France), VikingGenetics (Denmark, Sweden, Finland), DHV-VIT (Germany) and CRV (The Netherlands, Flanders). Each partner validated its own bulls using their national reference data and the combined data, respectively.ResultsCombining the data significantly increased the reliability of genomic predictions for bulls in all four populations. Reliabilities increased by 10%, compared to reliabilities obtained with national reference populations alone, when they were averaged over countries and the traits evaluated. For different traits and countries, the increase in reliability ranged from 2% to 19%.ConclusionsGenomic selection programs benefit greatly from combining data from several closely related populations into a single large reference population.

Impacts of both reference population size and inclusion of a residual polygenic effect on the accuracy of genomic prediction

BackgroundThe purpose of this work was to study the impact of both the size of genomic reference populations and the inclusion of a residual polygenic effect on dairy cattle genetic evaluations enhanced with genomic information.MethodsDirect genomic values were estimated for German Holstein cattle with a genomic BLUP model including a residual polygenic effect. A total of 17,429 genotyped Holstein bulls were evaluated using the phenotypes of 44 traits. The Interbull genomic validation test was implemented to investigate how the inclusion of a residual polygenic effect impacted genomic estimated breeding values.ResultsAs the number of reference bulls increased, both the variance of the estimates of single nucleotide polymorphism effects and the reliability of the direct genomic values of selection candidates increased. Fitting a residual polygenic effect in the model resulted in less biased genome-enhanced breeding values and decreased the correlation between direct genomic values and estimated breeding values of sires in the reference population.ConclusionsGenetic evaluation of dairy cattle enhanced with genomic information is highly effective in increasing reliability, as well as using large genomic reference populations. We found that fitting a residual polygenic effect reduced the bias in genome-enhanced breeding values, decreased the correlation between direct genomic values and sires estimated breeding values and made genome-enhanced breeding values more consistent in mean and variance as is the case for pedigree-based estimated breeding values.

Combined analysis of data from two granddaughter designs: A simple strategy for QTL confirmation and increasing experimental power in dairy cattle

A joint analysis of five paternal half-sib Holstein families that were part of two different granddaughter designs (ADR- or Inra-design) was carried out for five milk production traits and somatic cell score in order to conduct a QTL confirmation study and to increase the experimental power. Data were exchanged in a coded and standardised form. The combined data set (JOINT-design) consisted of on average 231 sires per grandsire. Genetic maps were calculated for 133 markers distributed over nine chromosomes. QTL analyses were performed separately for each design and each trait. The results revealed QTL for milk production on chromosome 14, for milk yield on chromosome 5, and for fat content on chromosome 19 in both the ADR- and the Inra-design (confirmed within this study). Some QTL could only be mapped in either the ADR- or in the Inra-design (not confirmed within this study). Additional QTL previously undetected in the single designs were mapped in the JOINT-design for fat yield (chromosome 19 and 26), protein yield (chromosome 26), protein content (chromosome 5), and somatic cell score (chromosome 2 and 19) with genomewide significance. This study demonstrated the potential benefits of a combined analysis of data from different granddaughter designs.

A single-step genomic model with direct estimation of marker effects

Compared with the currently widely used multi-step genomic models for genomic evaluation, single-step genomic models can provide more accurate genomic evaluation by jointly analyzing phenotypes and genotypes of all animals and can properly correct for the effect of genomic preselection on genetic evaluations. The objectives of this study were to introduce a single-step genomic model, allowing a direct estimation of single nucleotide polymorphism (SNP) effects, and to develop efficient computing algorithms for solving equations of the single-step SNP model. We proposed an alternative to the current single-step genomic model based on the genomic relationship matrix by including an additional step for estimating the effects of SNP markers. Our single-step SNP model allowed flexible modeling of SNP effects in terms of the number and variance of SNP markers. Moreover, our single-step SNP model included a residual polygenic effect with trait-specific variance for reducing inflation in genomic prediction. A kernel calculation of the SNP model involved repeated multiplications of the inverse of the pedigree relationship matrix of genotyped animals with a vector, for which numerical methods such as preconditioned conjugate gradients can be used. For estimating SNP effects, a special updating algorithm was proposed to separate residual polygenic effects from the SNP effects. We extended our single-step SNP model to general multiple-trait cases. By taking advantage of a block-diagonal (co)variance matrix of SNP effects, we showed how to estimate multivariate SNP effects in an efficient way. A general prediction formula was derived for candidates without phenotypes, which can be used for frequent, interim genomic evaluations without running the whole genomic evaluation process. We discussed various issues related to implementation of the single-step SNP model in Holstein populations with an across-country genomic reference population.

A whole genome scan for differences in recombination rates among three Bos taurus breeds

Abstract. Twenty paternal half-sib families of a granddaughter design were genotyped for 265 genetic markers, most of them microsatellites. These were 16 Holstein families, 3 Simmental families, and 1 Brown Swiss family. The number of sires per breed was 872, 170, and 32, respectively. Two-point recombination rates were estimated both jointly for all breeds and each single breed separately. Of 1168 marker intervals, 865 provided estimates for at least two breeds. Differences between breeds were tested by likelihood ratio tests. Four marker intervals, representing three genomic regions on BTA19, BTA24, and BTA27, show a significant impact of the breed at a false discovery rate of 0.23 and indicate a genetic component of observed heterogeneity of recombination. The variability of recombination rates between cattle breeds might not be a common feature of the whole genome, but rather might be restricted to certain chromosomal segments. Thus, attention should be paid to heterogeneities when pooling data of such regions from different breeds.

Data transformation for rank reduction in multi-trait MACE model for international bull comparison

Since many countries use multiple lactation random regression test day models in national evaluations for milk production traits, a random regression multiple across-country evaluation (MACE) model permitting a variable number of correlated traits per country should be used in international dairy evaluations. In order to reduce the number of within country traits for international comparison, three different MACE models were implemented based on German daughter yield deviation data and compared to the random regression MACE. The multiple lactation MACE model analysed daughter yield deviations on a lactation basis reducing the rank from nine random regression coefficients to three lactations. The lactation breeding values were very accurate for old bulls, but not for the youngest bulls with daughters with short lactations. The other two models applied principal component analysis as the dimension reduction technique: one based on eigenvalues of a genetic correlation matrix and the other on eigenvalues of a combined lactation matrix. The first one showed that German data can be transformed from nine traits to five eigenfunctions without losing much accuracy in any of the estimated random regression coefficients. The second one allowed performing rank reductions to three eigenfunctions without having the problem of young bulls with daughters with short lactations.

Technical note: Equivalent genomic models with a residual polygenic effect

Routine genomic evaluations in animal breeding are usually based on either a BLUP with genomic relationship matrix (GBLUP) or single nucleotide polymorphism (SNP) BLUP model. For a multi-step genomic evaluation, these 2 alternative genomic models were proven to give equivalent predictions for genomic reference animals. The model equivalence was verified also for young genotyped animals without phenotypes. Due to incomplete linkage disequilibrium of SNP markers to genes or causal mutations responsible for genetic inheritance of quantitative traits, SNP markers cannot explain all the genetic variance. A residual polygenic effect is normally fitted in the genomic model to account for the incomplete linkage disequilibrium. In this study, we start by showing the proof that the multi-step GBLUP and SNP BLUP models are equivalent for the reference animals, when they have a residual polygenic effect included. Second, the equivalence of both multi-step genomic models with a residual polygenic effect was also verified for young genotyped animals without phenotypes. Additionally, we derived formulas to convert genomic estimated breeding values of the GBLUP model to its components, direct genomic values and residual polygenic effect. Third, we made a proof that the equivalence of these 2 genomic models with a residual polygenic effect holds also for single-step genomic evaluation. Both the single-step GBLUP and SNP BLUP models lead to equal prediction for genotyped animals with phenotypes (e.g., reference animals), as well as for (young) genotyped animals without phenotypes. Finally, these 2 single-step genomic models with a residual polygenic effect were proven to be equivalent for estimation of SNP effects, too.