Integration of Interbull's multiple across-country evaluation approach breeding values into the multiple-trait single-step random regression test-day genetic evaluation for yield traits of Australian Red breeds.

Interbull's multiple across-country evaluaftion provides national breeding organizations with breeding values for internationally used bulls, which integrate performance data obtained in different breeding populations, environments, and production systems. However, breeding value-based selection decisions on domestic individuals born to foreign sires can only benefit from Interbull breeding values if they are integrated such that their information can contribute to the breeding values of all related domestic animals. For that purpose, several methods have been proposed which either model Interbull breeding values as prior information in a Bayesian approach, as additional pseudo data points, or as correlated traits, where these methods also differ in their software and parameterization requirements. Further, the complexity of integration also depends on the traits and genetic evaluation models. Especially random regression models require attention because of the dimensionality discrepancy between the number of Interbull breeding values and the number of modeled genetic effects. This paper presents the results from integrating 16,063 Interbull breeding values into the domestic single-step random regression test-day model for milk, fat, and protein yield for Australian Red dairy cattle breeds. Interbull breeding values were modeled as pseudo data points with data point-specific residual variances derived within animal across traits, ignoring relationships between integrated animals. Results suggest that the integration was successful with regard to alignment of Interbull breeding values with their domestic equivalent as well as with regard to the individual and population-wide increase in reliabilities. Depending on the relationship structure between integration candidates, further work is required to account for those relationships in a computationally feasible manner. Other traits with separate parity effects nationally could use a similar approach, even if not modeled with a test-day model.

Sharing of either phenotypes or genetic variants can increase the accuracy of genomic prediction of feed efficiency.

BACKGROUND: Sharing individual phenotype and genotype data between countries is complex and fraught with potential errors, while sharing summary statistics of genome-wide association studies (GWAS) is relatively straightforward, and thus would be especially useful for traits that are expensive or difficult-to-measure, such as feed efficiency. Here we examined: (1) the sharing of individual cow data from international partners; and (2) the use of sequence variants selected from GWAS of international cow data to evaluate the accuracy of genomic estimated breeding values (GEBV) for residual feed intake (RFI) in Australian cows. RESULTS: GEBV for RFI were estimated using genomic best linear unbiased prediction (GBLUP) with 50k or high-density single nucleotide polymorphisms (SNPs), from a training population of 3797 individuals in univariate to trivariate analyses where the three traits were RFI phenotypes calculated using 584 Australian lactating cows (AUSc), 824 growing heifers (AUSh), and 2526 international lactating cows (OVE). Accuracies of GEBV in AUSc were evaluated by either cohort-by-birth-year or fourfold random cross-validations. GEBV of AUSc were also predicted using only the AUS training population with a weighted genomic relationship matrix constructed with SNPs from the 50k array and sequence variants selected from a meta-GWAS that included only international datasets. The genomic heritabilities estimated using the AUSc, OVE and AUSh datasets were moderate, ranging from 0.20 to 0.36. The genetic correlations (rg) of traits between heifers and cows ranged from 0.30 to 0.95 but were associated with large standard errors. The mean accuracies of GEBV in Australian cows were up to 0.32 and almost doubled when either overseas cows, or both overseas cows and AUS heifers were included in the training population. They also increased when selected sequence variants were combined with 50k SNPs, but with a smaller relative increase. CONCLUSIONS: The accuracy of RFI GEBV increased when international data were used or when selected sequence variants were combined with 50k SNP array data. This suggests that if direct sharing of data is not feasible, a meta-analysis of summary GWAS statistics could provide selected SNPs for custom panels to use in genomic selection programs. However, since this finding is based on a small cross-validation study, confirmation through a larger study is recommended.

Short communication: Implementation of a breeding value for heat tolerance in Australian dairy cattle.

Excessive ambient temperature and humidity can impair milk production and fertility of dairy cows. Selection for heat-tolerant animals is one possible option to mitigate the effects of heat stress. To enable selection for this trait, we describe the development of a heat tolerance breeding value for Australian dairy cattle. We estimated the direct genomic values of decline in milk, fat, and protein yield per unit increase of temperature-humidity index (THI) using 46,726 single nucleotide polymorphisms and a reference population of 2,236 sires and 11,853 cows for Holsteins and 506 sires and 4,268 cows for Jerseys. This new direct genomic value is the Australian genomic breeding value for heat tolerance (HT ABVg). The components of the HT ABVg are the decline in milk, fat, and protein per unit increase in THI when THI increases above the threshold of 60. These components are weighted by their respective economic values, assumed to be equivalent to the weights applied to milk, fat, and protein yield in the Australian selection indices. Within each breed, the HT ABVg is then standardized to have a mean of 100 and standard deviation (SD) of 5, which is consistent with the presentation of breeding values for many other traits in Australia. The HT ABVg ranged from -4 to +3 SD in Holsteins and -3 to +4 SD in Jerseys. The mean reliabilities of HT ABVg among validation sires, calculated from the prediction error variance and additive genetic variance, were 38% in both breeds. The range in ABVg and their reliability suggests that HT can be improved using genomic selection. There has been a deterioration in the genetic trend of HT, and to moderate the decline it is suggested that the HT ABVg should be included in a multitrait economic index with other traits that contribute to farm profit.