The importance of historically popular sires on the accuracy of genomic predictions of young animals in the US Holstein population.

The dairy industry is known for its extensive use of artificial insemination, which has resulted in a population where most animals can be traced back to only a few sires. Due to their relatedness to the population, old influential sires could still contribute to the accuracy of genomic predictions. The objective of the study was to identify the impact of historically influential sires on the recent population. This was tested by constructing a genomic relationship matrix using recursion with different sets of sires. Differences in prediction accuracies with different sets are indicative of how important each set is. Recursion coefficients linking young animals to those sets reveal the relative importance of specific sires to the prediction accuracy of recent animals. The data included ∼10 million scores for stature and fore udder attachment (FUA) measured from 1983. Genotypes of 569,404 animals were available. Sire sets included the 100 most popular sires born within different time periods. Computations were with single-step genomic BLUP. In general, the younger sires had higher prediction accuracies than the oldest sires, even though they generally have fewer progeny. The accuracy of evaluation for stature was increased from 0.54 with the most popular sires born before 1981 to 0.69 with sires born from 2001 to 2010, while the accuracy for FUA increased from 0.47 to 0.61. The accuracy achieved using the overall 100 most used sires was 0.66 for stature and 0.58 for FUA. All 100 sires from each period were combined in a subset to determine the importance of each sire relative to all 400 animals in the combined subset. The highest relative impact of a sire that was born within the different time sets was 1.97 for Valiant (before 1981), 1.94 for Blackstar (1981 to 1990), 4.38 for Shottle (1991 to 2000), and 3.09 for Planet (2001 to 2010). The 3 sires among the 400 with the greatest impact were Shottle, Goldwyn (3.73), and Planet. The relative impact of a sire was not strongly related to the number of progeny. For instance, the relative impact of Durham with 34K progeny was 2.29, whereas the impact of O Man with 15K progeny was 3.13. The impact of a sire is also influenced by whether it was used as a sire of sires. Results show that younger sires are more relevant to the accuracy of breeding value prediction in the recent population.

Genome-Wide Association Study for Body Conformation Traits and Fitness in Czech Holsteins.

The aim of this study was a genome-wide association study (GWAS) on conformation traits using 25,486 genotyped Czech Holsteins, with 35,227 common SNPs for each genotype. Linear trait records were collected between 1995 and 2020. The Interbull information from Multiple Across Country Evaluation (MACE) was included for bulls that mostly had daughter records in a foreign country. When using the Bonferroni correction, the number of SNPs that were either significant or approached the significance threshold was low-dairy capacity composite on BTA4, feet and legs composite BTA21, total score BTA10, stature BTA24, body depth BTA6, angularity BTA20, fore udder attachment BTA10. Without the Bonferroni correction, the total number of significant or near of significance SNPs was 32. The SNPs were localized on BTA1,2,4,5,6,7,8,18,22,25,26,28 for dairy capacity composite, BTA15,21 for feet and legs composite, BTA10 for total score, BTA24 stature, BTA6,23 body depth, BTA20 angularity, BTA2 rump angle, BTA9,10 rear legs rear view, BTA2,19 rear legs side view, BTA10 fore udder attachment, BTA2 udder depth, BTA10 rear udder height, BTA12 central alignment, BTA24 rear teat placement, BTA8,29 rear udder width. The results provide biological information for the improvement of body conformation and fitness in the Holstein population.

The prediction accuracies of linear-type traits in Czech Holstein cattle when using ssGBLUP or wssGBLUP.

The aim of this study was to assess the contribution of the weighted single-step genomic best linear unbiased prediction (wssGBLUP) method compared to the single-step genomic best linear unbiased prediction (ssGBLUP) method for genomic evaluation of 25 linear-type traits in the Czech Holstein cattle population. The nationwide database of linear-type traits with 6,99,681 records combined with deregressed proofs from Interbull (MACE method) was used as the input data. Genomic breeding values (GEBVs) were predicted based on these phenotypes using ssGBLUP and wssGBLUP methods using the BLUPF90 software. The bull validation test was employed which was based on comparing GEBVs of young bulls (N = 334) with no progeny in 2016. A minimum of 50 daughters with their own performance in 2020 was chosen to verify the contribution to the GEBV prediction, GEBV reliability, validation reliabilities (R2), and regression coefficients (b1). The results showed that the differences between the two methods were negligible. The low benefit of wssGBLUP may be due to the inclusion of a small number of SNPs; therefore, most predictions rely on polygenic relationships between animals. Nevertheless, the benefits of wssGBLUP analysis should be assessed with respect to specific population structures and given traits.

Animal breeding is based on statistical and mathematical approaches. With the development of computer technology, these procedures have become more efficient and have shed light upon an increasing amount of information, particularly in the field of molecular genetics. This results in a more accurate prediction of the breeding value. The single-step approach is the most popular genomic breeding value-prediction method. Single-step genomic BLUP (ssGBLUP) assumes that all single nucleotide polymorphisms (SNPs) explain the same fraction of genetic variance. In contrast, unequal variance and SNP weights were considered in weighted ssGBLUP (wssGBLUP). The aim of this study was to assess the contribution of wssGBLUP compared to ssGBLUP for the genomic evaluation of 25 linear-type traits in Czech Holstein cattle. We can conclude that no significant differences were found between these methods for the evaluated traits, probably because of the low number of included SNPs, which may not cover all significant SNPs in the genome.

Optimal definition of contemporary groups for crossbred pigs in a joint purebred and crossbred genetic evaluation.

In the pig industry, purebred animals are raised in nucleus herds and selected to produce crossbred progeny to perform in commercial environments. Crossbred and purebred performances are different, correlated traits. All purebreds in a pen have their performance assessed together at the end of a performance test. However, only selected crossbreds are removed (based on visual inspection) and measured at different times creating many small contemporary groups (CGs). This may reduce estimated breeding value (EBV) prediction accuracies. Considering this sequential recording of crossbreds, the objective was to investigate the impact of different CG definitions on genetic parameters and EBV prediction accuracy for crossbred traits. Growth rate (GP) and ultrasound backfat (BFP) records were available for purebreds. Lifetime growth (GX) and backfat (BFX) were recorded on crossbreds. Different CGs were tested: CG_all included farm, sex, birth year, and birth week; CG_week added slaughter week; and CG_day used slaughter day instead of week. Data of 124,709 crossbreds were used. The purebred phenotypes (62,274 animals) included three generations of purebred ancestors of these crossbreds and their CG mates. Variance components for four-trait models with different CG definitions were estimated with average information restricted maximum likelihood. Purebred traits' variance components remained stable across CG definitions and varied slightly for BFX. Additive genetic variances (and heritabilities) for GX fluctuated more: 812 ± 36 (0.28 ± 0.01), 257 ± 15 (0.17 ± 0.01), and 204 ± 13 (0.15 ± 0.01) for CG_all, CG_week, and CG_day, respectively. Age at slaughter (AAS) and hot carcass weight (HCW) adjusted for age were investigated as alternatives for GX. Both have potential for selection but lower heritabilities compared with GX: 0.21 ± 0.01 (0.18 ± 0.01), 0.16 ± 0.02 (0.16 + 0.01), and 0.10 ± 0.01 (0.14 ± 0.01) for AAS (HCW) using CG_all, CG_week, and CG_day, respectively. The predictive ability, linear regression (LR) accuracy, bias, and dispersion of crossbred traits in crossbreds favored CG_day, but correlations with unadjusted phenotypes favored CG_all. In purebreds, CG_all showed the best LR accuracy, while showing small relative differences in bias and dispersion. Different CG scenarios showed no relevant impact on BFX EBV. This study shows that different CG definitions may affect evaluation stability and animal ranking. Results suggest that ignoring slaughter dates in CG is more appropriate for estimating crossbred trait EBV for purebred animals.

Genomic predictions in purebreds with a multibreed genomic relationship matrix1.

Combining breeds in a multibreed evaluation can have a negative impact on prediction accuracy, especially if single nucleotide polymorphism (SNP) effects differ among breeds. The aim of this study was to evaluate the use of a multibreed genomic relationship matrix (G), where SNP effects are considered to be unique to each breed, that is, nonshared. This multibreed G was created by treating SNP of different breeds as if they were on nonoverlapping positions on the chromosome, although, in reality, they were not. This simple setup may avoid spurious Identity by state (IBS) relationships between breeds and automatically considers breed-specific allele frequencies. This scenario was contrasted to a regular multibreed evaluation where all SNPs were shared, that is, the same position, and to single-breed evaluations. Different SNP densities (9k and 45k) and different effective population sizes (Ne) were tested. Five breeds mimicking recent beef cattle populations that diverged from the same historical population were simulated using different selection criteria. It was assumed that quantitative trait locus (QTL) effects were the same over all breeds. For the recent population, generations 1-9 had approximately half of the animals genotyped, whereas all animals in generation 10 were genotyped. Generation 10 animals were set for validation; therefore, each breed had a validation group. Analyses were performed using single-step genomic best linear unbiased prediction. Prediction accuracy was calculated as the correlation between true (T) and genomic estimated breeding values (GEBV). Accuracies of GEBV were lower for the larger Ne and low SNP density. All three evaluation scenarios using 45k resulted in similar accuracies, suggesting that the marker density is high enough to account for relationships and linkage disequilibrium with QTL. A shared multibreed evaluation using 9k resulted in a decrease of accuracy of 0.08 for a smaller Ne and 0.12 for a larger Ne. This loss was mostly avoided when markers were treated as nonshared within the same G matrix. A G matrix with nonshared SNP enables multibreed evaluations without considerably changing accuracy, especially with limited information per breed.