Domestic estimated breeding values and genomic enhanced breeding values of bulls in comparison with their foreign genomic enhanced breeding values.

Estimated breeding values (EBVs) and genomic enhanced breeding values (GEBVs) for milk production of young genotyped Holstein bulls were predicted using a conventional BLUP - Animal Model, a method fitting regression coefficients for loci (RRBLUP), a method utilizing the realized genomic relationship matrix (GBLUP), by a single-step procedure (ssGBLUP) and by a one-step blending procedure. Information sources for prediction were the nation-wide database of domestic Czech production records in the first lactation combined with deregressed proofs (DRP) from Interbull files (August 2013) and domestic test-day (TD) records for the first three lactations. Data from 2627 genotyped bulls were used, of which 2189 were already proven under domestic conditions. Analyses were run that used Interbull values for genotyped bulls only or that used Interbull values for all available sires. Resultant predictions were compared with GEBV of 96 young foreign bulls evaluated abroad and whose proofs were from Interbull method GMACE (August 2013) on the Czech scale. Correlations of predictions with GMACE values of foreign bulls ranged from 0.33 to 0.75. Combining domestic data with Interbull EBVs improved prediction of both EBV and GEBV. Predictions by Animal Model (traditional EBV) using only domestic first lactation records and GMACE values were correlated by only 0.33. Combining the nation-wide domestic database with all available DRP for genotyped and un-genotyped sires from Interbull resulted in an EBV correlation of 0.60, compared with 0.47 when only Interbull data were used. In all cases, GEBVs had higher correlations than traditional EBVs, and the highest correlations were for predictions from the ssGBLUP procedure using combined data (0.75), or with all available DRP from Interbull records only (one-step blending approach, 0.69). The ssGBLUP predictions using the first three domestic lactation records in the TD model were correlated with GMACE predictions by 0.69, 0.64 and 0.61 for milk yield, protein yield and fat yield, respectively.

Contribution of domestic production records, Interbull estimated breeding values, and single nucleotide polymorphism genetic markers to the single-step genomic evaluation of milk production.

Estimated breeding values (EBV) for first-lactation milk production of Holstein cattle in the Czech Republic were calculated using a conventional animal model and by single-step prediction of the genomic enhanced breeding value. Two overlapping data sets of milk production data were evaluated: (1) calving years 1991 to 2006, with 861,429 lactations and 1,918,901 animals in the pedigree and (2) calving years 1991 to 2010, with 1,097,319 lactations and 1,906,576 animals in the pedigree. Global Interbull (Uppsala, Sweden) deregressed proofs of 114,189 bulls were used in the analyses. Reliabilities of Interbull values were equivalent to an average of 8.53 effective records, which were used in a weighted analysis. A total of 1,341 bulls were genotyped using the Illumina BovineSNP50 BeadChip V2 (Illumina Inc., San Diego, CA). Among the genotyped bulls were 332 young bulls with no daughters in the first data set but more than 50 daughters (88.41, on average) with performance records in the second data set. For young bulls, correlations of EBV and genomic enhanced breeding value before and after progeny testing, corresponding average expected reliabilities, and effective daughter contributions (EDC) were calculated. The reliability of prediction pedigree EBV of young bulls was 0.41, corresponding to EDC=10.6. Including Interbull deregressed proofs improved the reliability of prediction by EDC=13.4 and including genotyping improved prediction reliability by EDC=6.2. Total average expected reliability of prediction reached 0.67, corresponding to EDC=30.2. The combination of domestic and Interbull sources for both genotyped and nongenotyped animals is valuable for improving the accuracy of genetic prediction in small populations of dairy cattle.

Effect of type traits on functional longevity of Czech Holstein cows estimated from a Cox proportional hazards model.

Relationships between conformation traits and functional longevity in Holstein cows were evaluated using survival analysis. Functional longevity was defined as the number of days between the first calving and culling; that is, length of productive life. The data set consisted of 116,369 Holstein cows that first calved from 2003 to 2008. All cows used in the analysis were scored for conformation between d 30 and d 210 of their first lactation. The data included 48% censored records. Analyses were done separately for 20 linear descriptive type traits, 6 composite traits, and height at sacrum measured in centimeters. Cox proportional hazard models were fitted to analyze data. The hazard function was described as the product of a baseline hazard function and the time-independent effects of age at first calving and sire (random), and the time-dependent effects of stage of lactation and lactation number, herd, year and season, herd size, and 305-d milk production. The strongest relationship between a composite trait and functional longevity was for dairy form, followed by udder and final score. Among the descriptive type traits, the strongest relationships with longevity were found for body condition score, angularity, traits related to udder attachment, and udder depth. Foot and leg traits showed substantially lower effect on functional longevity, and the effect of foot angle was minimal. Functional longevity declined with decreased body condition score of cows. Cows with deep udders had significantly lower functional survival compared with cows with shallow udders. In addition, weak central ligament was associated with significant reduction of cow longevity. For dairy form and angularity, cows classified as very good were the worst with respect to longevity, whereas cows classified as poor were the best. An intermediate optimum was evident for rear legs rear view and rear legs set (side view), whereas cows with sickled legs had lower longevity than cows with straighter legs.

Non-additive effects on milk production in Czech dairy cows.

Crossbreeding effects on milk production traits of Czech dual-purpose and dairy cattle breeds were estimated. Nearly 370,000 cows with known gene proportions from Czech Pied, Ayrshire or Holstein cattle were selected from the national milk recording data base. Single-trait animal models were calculated for milk, fat and protein yield, fat and protein content. The model of Dickerson including additive, additive maternal, heterotic and recombination effects was used for the part of the animal model describing the crossbreeding effects in all calculations. For milk yield, the additive genetic effect (defined as deviation from Czech Pied cattle) was 850-900 kg for Holstein and 240-480 kg for Ayrshire. The maternal effects were low and negative. Low significant positive heterotic effects were observed being up to approximately 100 kg for Czech Pied x Holstein. The recombination effects were negative and statistically significant for Czech Pied x Holstein. The results for fat and protein yield were similar to the results for milk yield. For fat and protein content, nearly no statistically significant crossbreeding effects were found.