Value of sharing cow reference population between countries on reliability of genomic prediction for milk yield traits.

Increasing the reliability of genomic prediction (GP) of economic traits in the pasture-based dairy production systems of New Zealand (NZ) and Australia (AU) is important to both countries. This study assessed if sharing cow phenotype and genotype data of NZ and AU improves the reliability of GP for NZ bulls. Data from approximately 32,000 NZ genotyped cows and their contemporaries were included in the May 2018 routine genetic evaluation of the Australian Dairy cattle in an attempt to provide consistent phenotypes for both countries. After the genetic evaluation, deregressed proofs of cows were calculated for milk yield traits. The April 2018 multiple across-country evaluation of Interbull was also used to calculate deregressed proofs for bulls on the NZ scale. Approximately 1,178 Jersey (Jer) and 6,422 Holstein (Hol) bulls had genotype and phenotype data. In addition to NZ cows, phenotype data of close to 60,000 genotyped Australian (AU) cows from the same genetic evaluation run as NZ cows were used. All AU and NZ females were genotyped using low-density SNP chips (<10K SNP) and were imputed first to 50K and then to ∼600K (referred to as high density; HD). We used up to 98,000 animals in the reference populations, both by expanding the NZ reference set (cow, bull, single breed to multi-breed set) and by adding AU cows. Reliabilities of GP were calculated for 508 Jer and 1,251 Hol bulls whose sires are not included in the reference set (RS) to ensure that real differences are not masked by close relationships. The GP was tested using 50K or high-density SNP chip using genomic BLUP in bivariate (considering country as a trait) or single trait models. The RS that gave the highest reliability for each breed were also tested using a hybrid GP method that combines expectation maximization with Bayes R. The addition of the AU cows to an NZ RS that included either NZ cows only, or cows and bulls, improved the reliability of GP for both NZ Hol and Jer validation bulls for all traits. Using single breed reference populations also increased reliability when NZ crossbred cows were added to reference populations that included only purebred NZ bulls and cows and AU cows. The full multi-breed RS (all NZ cows and bulls and AU cows) provided similar reliabilities in NZ Hol bulls, when compared with the single breed reference with crossbred NZ cows. For Jer validation bulls, the RS that included Jer cows and bulls and crossbred cows from NZ and Jer cows from AU was marginally better than the all-breed, all-country RS. In terms of reliability, the advantage of the HD SNP chip was small but captured more of the genomic variance than the 50K, particularly for Hol. The expectation maximization Bayes R GP method was slightly (up to 3 percentage points) better than genomic BLUP. We conclude that GP of milk production traits in NZ bulls improves by up to 7 percentage points in reliability by expanding the NZ reference population to include AU cows.

Haile-Mariam M ，MacLeod IM ，Bolormaa S ，Schrooten C ，O'Connor E ，de Jong G ，Daetwyler HD ，Pryce JE ... -  《-》

Including overseas performance information in genomic evaluations of Australian dairy cattle.

In dairy cattle, the rate of genetic gain from genomic selection depends on reliability of direct genomic values (DGV). One option to increase reliabilities could be to increase the size of the reference set used for prediction, by using genotyped bulls with daughter information in countries other than the evaluating country. The increase in reliabilities of DGV from using this information will depend on the extent of genotype by environment interaction between the evaluating country and countries contributing information, and whether this is correctly accounted for in the prediction method. As the genotype by environment interaction between Australia and Europe or North America is greater than between Europe and North America for most dairy traits, ways of including information from other countries in Australian genomic evaluations were examined. Thus, alternative approaches for including information from other countries and their effect on the reliability and bias of DGV of selection candidates were assessed. We also investigated the effect of including overseas (OS) information on reliabilities of DGV for selection candidates that had weaker relationships to the current Australian reference set. The DGV were predicted either using daughter trait deviations (DTD) for the bulls with daughters in Australia, or using this information as well as OS information by including deregressed proofs (DRP) from Interbull for bulls with only OS daughters in either single trait or bivariate models. In the bivariate models, DTD and DRP were considered as different traits. Analyses were performed for Holstein and Jersey bulls for milk yield traits, fertility, cell count, survival, and some type traits. For Holsteins, the data used included up to 3,580 bulls with DTD and up to 5,720 bulls with only DRP. For Jersey, about 900 bulls with DTD and 1,820 bulls with DRP were used. Bulls born after 2003 and genotyped cows that were not dams of genotyped bulls were used for validation. The results showed that the combined use of DRP on bulls with OS daughters only and DTD for Australian bulls in either the single trait or bivariate model increased the coefficient of determination [(R(2)) (DGV,DTD)] in the validation set, averaged across 6 main traits, by 3% in Holstein and by 5% in Jersey validation bulls relative to the use of DTD only. Gains in reliability and unbiasedness of DGV were similar for the single trait and bivariate models for production traits, whereas the bivariate model performed slightly better for somatic cell count in Holstein. The increase in R(2) (DGV,DTD) as a result of using bulls with OS daughters was relatively higher for those bulls and cows in the validation sets that were less related to the current reference set. For example, in Holstein, the average increase in R(2) for milk yield traits when DTD and DRP were used in a single trait model was 23% in the least-related cow group, but only 3% in the most-related cow group. In general, for both breeds the use of DTD from domestic sources and DRP from Interbull in a single trait or bivariate model can increase reliability of DGV for selection candidates.

Haile-Mariam M ，Pryce JE ，Schrooten C ，Hayes BJ ... -  《-》

Optimizing genomic prediction for Australian Red dairy cattle.

The reliability of genomic prediction is influenced by several factors, including the size of the reference population, which makes genomic prediction for breeds with a relatively small population size challenging, such as Australian Red dairy cattle. Including other breeds in the reference population may help to increase the size of the reference population, but the reliability of genomic prediction is also influenced by the relatedness between the reference and validation population. Our objective was to optimize the reference population for genomic prediction of Australian Red dairy cattle. A reference population comprising up to 3,248 Holstein bulls, 48,386 Holstein cows, 807 Jersey bulls, 8,734 Jersey cows, and 3,041 Australian Red cows and a validation population with between 208 and 224 Australian Red Bulls were used, with records for milk, fat, and protein yield, somatic cell count, fertility, and survival. Three different analyses were implemented: single-trait genomic best linear unbiased predictor (GBLUP), multi-trait GBLUP, and single-trait Bayes R, using 2 different medium-density SNP panels: the standard 50K chip and a custom array of variants that were expected to be enriched for causative mutations. Various reference populations were constructed containing the Australian Red cows and all Holstein and Jersey bulls and cows, all Holstein and Jersey bulls, all Holstein bulls and cows, all Holstein bulls, and a subset of the Holstein individuals varying the relatedness between Holsteins and Australian Reds and the number of Holsteins. Varying the relatedness between reference and validation populations only led to small changes in reliability. Whereas adding a limited number of closely related Holsteins increased reliabilities compared with within-breed prediction, increasing the number of Holsteins decreased the reliability. The multi-trait GBLUP, which considered the same trait in different breeds as correlated traits, yielded higher reliabilities than the single-trait GBLUP. Bayes R yielded lower reliabilities than multi-trait GBLUP and outperformed single-trait GBLUP for larger reference populations. Our results show that increasing the size of a multi-breed reference population may result in a reference population dominated by one breed and reduce the reliability to predict in other breeds.

van den Berg I ，MacLeod IM ，Reich CM ，Breen EJ ，Pryce JE ... -  《-》

Accuracy of genomic predictions in Gyr (Bos indicus) dairy cattle.

Genomic selection may accelerate genetic progress in breeding programs of indicine breeds when compared with traditional selection methods. We present results of genomic predictions in Gyr (Bos indicus) dairy cattle of Brazil for milk yield (MY), fat yield (FY), protein yield (PY), and age at first calving using information from bulls and cows. Four different single nucleotide polymorphism (SNP) chips were studied. Additionally, the effect of the use of imputed data on genomic prediction accuracy was studied. A total of 474 bulls and 1,688 cows were genotyped with the Illumina BovineHD (HD; San Diego, CA) and BovineSNP50 (50K) chip, respectively. Genotypes of cows were imputed to HD using FImpute v2.2. After quality check of data, 496,606 markers remained. The HD markers present on the GeneSeek SGGP-20Ki (15,727; Lincoln, NE), 50K (22,152), and GeneSeek GGP-75Ki (65,018) were subset and used to assess the effect of lower SNP density on accuracy of prediction. Deregressed breeding values were used as pseudophenotypes for model training. Data were split into reference and validation to mimic a forward prediction scheme. The reference population consisted of animals whose birth year was ≤2004 and consisted of either only bulls (TR1) or a combination of bulls and dams (TR2), whereas the validation set consisted of younger bulls (born after 2004). Genomic BLUP was used to estimate genomic breeding values (GEBV) and reliability of GEBV (R) was based on the prediction error variance approach. Reliability of GEBV ranged from ∼0.46 (FY and PY) to 0.56 (MY) with TR1 and from 0.51 (PY) to 0.65 (MY) with TR2. When averaged across all traits, R were substantially higher (R of TR1 = 0.50 and TR2 = 0.57) compared with reliabilities of parent averages (0.35) computed from pedigree data and based on diagonals of the coefficient matrix (prediction error variance approach). Reliability was similar for all the 4 marker panels using either TR1 or TR2, except that imputed HD cow data set led to an inflation of reliability. Reliability of GEBV could be increased by enlarging the limited bull reference population with cow information. A reduced panel of ∼15K markers resulted in reliabilities similar to using HD markers. Reliability of GEBV could be increased by enlarging the limited bull reference population with cow information.

Boison SA ，Utsunomiya ATH ，Santos DJA ，Neves HHR ，Carvalheiro R ，Mészáros G ，Utsunomiya YT ，do Carmo AS ，Verneque RS ，Machado MA ，Panetto JCC ，Garcia JF ，Sölkner J ，da Silva MVGB ... -  《-》

Efficiency of multi-breed genomic selection for dairy cattle breeds with different sizes of reference population.

Single-breed genomic selection (GS) based on medium single nucleotide polymorphism (SNP) density (~50,000; 50K) is now routinely implemented in several large cattle breeds. However, building large enough reference populations remains a challenge for many medium or small breeds. The high-density BovineHD BeadChip (HD chip; Illumina Inc., San Diego, CA) containing 777,609 SNP developed in 2010 is characterized by short-distance linkage disequilibrium expected to be maintained across breeds. Therefore, combining reference populations can be envisioned. A population of 1,869 influential ancestors from 3 dairy breeds (Holstein, Montbéliarde, and Normande) was genotyped with the HD chip. Using this sample, 50K genotypes were imputed within breed to high-density genotypes, leading to a large HD reference population. This population was used to develop a multi-breed genomic evaluation. The goal of this paper was to investigate the gain of multi-breed genomic evaluation for a small breed. The advantage of using a large breed (Normande in the present study) to mimic a small breed is the large potential validation population to compare alternative genomic selection approaches more reliably. In the Normande breed, 3 training sets were defined with 1,597, 404, and 198 bulls, and a unique validation set included the 394 youngest bulls. For each training set, estimated breeding values (EBV) were computed using pedigree-based BLUP, single-breed BayesC, or multi-breed BayesC for which the reference population was formed by any of the Normande training data sets and 4,989 Holstein and 1,788 Montbéliarde bulls. Phenotypes were standardized by within-breed genetic standard deviation, the proportion of polygenic variance was set to 30%, and the estimated number of SNP with a nonzero effect was about 7,000. The 2 genomic selection (GS) approaches were performed using either the 50K or HD genotypes. The correlations between EBV and observed daughter yield deviations (DYD) were computed for 6 traits and using the different prediction approaches. Compared with pedigree-based BLUP, the average gain in accuracy with GS in small populations was 0.057 for the single-breed and 0.086 for multi-breed approach. This gain was up to 0.193 and 0.209, respectively, with the large reference population. Improvement of EBV prediction due to the multi-breed evaluation was higher for animals not closely related to the reference population. In the case of a breed with a small reference population size, the increase in correlation due to multi-breed GS was 0.141 for bulls without their sire in reference population compared with 0.016 for bulls with their sire in reference population. These results demonstrate that multi-breed GS can contribute to increase genomic evaluation accuracy in small breeds.

Hozé C ，Fritz S ，Phocas F ，Boichard D ，Ducrocq V ，Croiseau P ... -  《-》