Using Sequence Variants in Linkage Disequilibrium with Causative Mutations to Improve Across-Breed Prediction in Dairy Cattle: A Simulation Study.

Sequence data are expected to increase the reliability of genomic prediction by containing causative mutations directly, especially in cases where low linkage disequilibrium between markers and causative mutations limits prediction reliability, such as across-breed prediction in dairy cattle. In practice, the causative mutations are unknown, and prediction with only variants in perfect linkage disequilibrium with the causative mutations is not realistic, leading to a reduced reliability compared to knowing the causative variants. Our objective was to use sequence data to investigate the potential benefits of sequence data for the prediction of genomic relationships, and consequently reliability of genomic breeding values. We used sequence data from five dairy cattle breeds, and a larger number of imputed sequences for two of the five breeds. We focused on the influence of linkage disequilibrium between markers and causative mutations, and assumed that a fraction of the causative mutations was shared across breeds and had the same effect across breeds. By comparing the loss in reliability of different scenarios, varying the distance between markers and causative mutations, using either all genome wide markers from commercial SNP chips, or only the markers closest to the causative mutations, we demonstrate the importance of using only variants very close to the causative mutations, especially for across-breed prediction. Rare variants improved prediction only if they were very close to rare causative mutations, and all causative mutations were rare. Our results show that sequence data can potentially improve genomic prediction, but careful selection of markers is essential.

van den Berg I ，Boichard D ，Guldbrandtsen B ，Lund MS ... -  《G3-Genes Genomes Genetics》

Sequence variants selected from a multi-breed GWAS can improve the reliability of genomic predictions in dairy cattle.

Sequence data can potentially increase the reliability of genomic predictions, because such data include causative mutations instead of relying on linkage disequilibrium (LD) between causative mutations and prediction variants. However, the location of the causative mutations is not known, and the presence of many variants that are in low LD with the causative mutations may reduce prediction reliability. Our objective was to investigate whether the use of variants at quantitative trait loci (QTL) that are identified in a multi-breed genome-wide association study (GWAS) for milk, fat and protein yield would increase the reliability of within- and multi-breed genomic predictions in Holstein, Jersey and Danish Red cattle. A wide range of scenarios that test different strategies to select prediction markers, for both within-breed and multi-breed prediction, were compared. For all breeds and traits, the use of variants selected from a multi-breed GWAS resulted in substantial increases in prediction reliabilities compared to within-breed prediction using a 50 K SNP array. Reliabilities depended highly on the choice of the prediction markers, and the scenario that led to the highest reliability varied between breeds and traits. While genomic correlations across breeds were low for genome-wide sequence variants, the effects of the QTL variants that yielded the highest reliabilities were highly correlated across breeds. Our results show that the use of sequence variants, which are located near peaks of QTL that are detected in a multi-breed GWAS, can increase reliability of genomic predictions.

van den Berg I ，Boichard D ，Lund MS 《-》

Accuracy of predicting genomic breeding values for residual feed intake in Angus and Charolais beef cattle.

In beef cattle, phenotypic data that are difficult and/or costly to measure, such as feed efficiency, and DNA marker genotypes are usually available on a small number of animals of different breeds or populations. To achieve a maximal accuracy of genomic prediction using the phenotype and genotype data, strategies for forming a training population to predict genomic breeding values (GEBV) of the selection candidates need to be evaluated. In this study, we examined the accuracy of predicting GEBV for residual feed intake (RFI) based on 522 Angus and 395 Charolais steers genotyped on SNP with the Illumina Bovine SNP50 Beadchip for 3 training population forming strategies: within breed, across breed, and by pooling data from the 2 breeds (i.e., combined). Two other scenarios with the training and validation data split by birth year and by sire family within a breed were also investigated to assess the impact of genetic relationships on the accuracy of genomic prediction. Three statistical methods including the best linear unbiased prediction with the relationship matrix defined based on the pedigree (PBLUP), based on the SNP genotypes (GBLUP), and a Bayesian method (BayesB) were used to predict the GEBV. The results showed that the accuracy of the GEBV prediction was the highest when the prediction was within breed and when the validation population had greater genetic relationships with the training population, with a maximum of 0.58 for Angus and 0.64 for Charolais. The within-breed prediction accuracies dropped to 0.29 and 0.38, respectively, when the validation populations had a minimal pedigree link with the training population. When the training population of a different breed was used to predict the GEBV of the validation population, that is, across-breed genomic prediction, the accuracies were further reduced to 0.10 to 0.22, depending on the prediction method used. Pooling data from the 2 breeds to form the training population resulted in accuracies increased to 0.31 and 0.43, respectively, for the Angus and Charolais validation populations. The results suggested that the genetic relationship of selection candidates with the training population has a greater impact on the accuracy of GEBV using the Illumina Bovine SNP50 Beadchip. Pooling data from different breeds to form the training population will improve the accuracy of across breed genomic prediction for RFI in beef cattle.

Chen L ，Schenkel F ，Vinsky M ，Crews DH Jr ，Li C ... -  《-》

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.

Steyn Y ，Lourenco DAL ，Misztal I 《-》

Use of meta-analyses and joint analyses to select variants in whole genome sequences for genomic evaluation: An application in milk production of French dairy cattle breeds.

As a result of the 1000 Bull Genome Project, it has become possible to impute millions of variants, with many of these potentially causative for traits of interest, for thousands of animals that have been genotyped with medium-density chips. This enormous source of data opens up very interesting possibilities for the inclusion of these variants in genomic evaluations. However, for computational reasons, it is not possible to include all variants in genomic evaluation procedures. One potential approach could be to select the most relevant variants based on the results of genome-wide association studies (GWAS); however, the identification of causative mutations is still difficult with this method, partly because of weak imputation accuracy for rare variants. To address this problem, this study assesses the ability of different approaches based on multi-breed GWAS (joint and meta-analyses) to identify single-nucleotide polymorphisms (SNP) for use in genomic evaluation in the 3 main French dairy cattle breeds. A total of 6,262 Holstein bulls, 2,434 Montbéliarde bulls, and 2,175 Normande bulls with daughter yield deviations for 5 milk production traits were imputed for 27 million variants. Within-breed and joint (including all 3 breeds) GWAS were performed and 3 models of meta-analysis were tested: fixed effect, random effect, and Z-score. Comparison of the results of within- and multi-breed GWAS showed that most of the quantitative trait loci identified using within-breed approaches were also found with multi-breed methods. However, the most significant variants identified in each region differed depending on the method used. To determine which approach highlighted the most predictive SNP for each trait, we used a marker-assisted best unbiased linear prediction model to evaluate lists of SNP generated by the different GWAS methods; each list contained between 25 and 2,000 candidate variants per trait, which were identified using a single within- or multi-breed GWAS approach. Among all the multi-breed methods tested in this study, variant selection based on meta-analysis (fixed effect) resulted in the most-accurate genomic evaluation (+1 to +3 points compared with other multi-breed approaches). However, the accuracies of genomic evaluation were always better when variants were selected using the results of within-breed GWAS. As has generally been found in studies of quantitative trait loci, these results suggest that part of the genetic variance of milk production traits is breed specific in Holstein, Montbéliarde, and Normande cattle.

Teissier M ，Sanchez MP ，Boussaha M ，Barbat A ，Hoze C ，Robert-Granie C ，Croiseau P ... -  《-》