Significance testing and genomic inflation factor using high‐density genotypes or whole‐genome sequence data

Significance testing for genome‐wide association study (GWAS) with increasing SNP density up to whole‐genome sequence data (WGS) is not straightforward, because of strong LD between SNP and population stratification. Therefore, the objective of this study was to investigate genomic control and different significance testing procedures using data from a commercial pig breeding scheme. A GWAS was performed in GCTA with data of 4,964 Large White pigs using medium density, high density or imputed whole‐genome sequence data, fitting a genomic relationship matrix based on a leave‐one–chromosome‐out approach to account for population structure. Subsequently, genomic inflation factors were assessed on whole‐genome level and the chromosome level. To establish a significance threshold, permutation testing, Bonferroni corrections using either the total number of SNPs or the number of independent chromosome fragments, and false discovery rates (FDR) using either the Benjamini–Hochberg procedure or the Benjamini and Yekutieli procedure were evaluated. We found that genomic inflation factors did not differ between different density genotypes but do differ between chromosomes. Also, the leave‐one‐chromosome‐out approach for GWAS or using the pedigree relationships did not account appropriately for population stratification and gave strong genomic inflation. Regarding different procedures for significance testing, when the aim is to find QTL regions that are associated with a trait of interest, we recommend applying the FDR following the Benjamini and Yekutieli approach to establish a significance threshold that is adjusted for multiple testing. When the aim is to pinpoint a specific mutation, the more conservative Bonferroni correction based on the total number of SNPs is more appropriate, till an appropriate method is established to adjust for the number of independent tests.     

Improving accuracy of genomic prediction in Brangus cattle by adding animals with imputed low‐density SNP genotypes

Reliable genomic prediction of breeding values for quantitative traits requires the availability of sufficient number of animals with genotypes and phenotypes in the training set. As of 31 October 2016, there were 3,797 Brangus animals with genotypes and phenotypes. These Brangus animals were genotyped using different commercial SNP chips. Of them, the largest group consisted of 1,535 animals genotyped by the GGP‐LDV4 SNP chip. The remaining 2,262 genotypes were imputed to the SNP content of the GGP‐LDV4 chip, so that the number of animals available for training the genomic prediction models was more than doubled. The present study showed that the pooling of animals with both original or imputed 40K SNP genotypes substantially increased genomic prediction accuracies on the ten traits. By supplementing imputed genotypes, the relative gains in genomic prediction accuracies on estimated breeding values (EBV) were from 12.60% to 31.27%, and the relative gain in genomic prediction accuracies on de‐regressed EBV was slightly small (i.e. 0.87%–18.75%). The present study also compared the performance of five genomic prediction models and two cross‐validation methods. The five genomic models predicted EBV and de‐regressed EBV of the ten traits similarly well. Of the two cross‐validation methods, leave‐one‐out cross‐validation maximized the number of animals at the stage of training for genomic prediction. Genomic prediction accuracy (GPA) on the ten quantitative traits was validated in 1,106 newly genotyped Brangus animals based on the SNP effects estimated in the previous set of 3,797 Brangus animals, and they were slightly lower than GPA in the original data. The present study was the first to leverage currently available genotype and phenotype resources in order to harness genomic prediction in Brangus beef cattle.     

Assessment of bagging GBLUP for whole‐genome prediction of broiler chicken traits

Bootstrap aggregation (bagging) is a resampling method known to produce more accurate predictions when predictors are unstable or when the number of markers is much larger than sample size, because of variance reduction capabilities. The purpose of this study was to compare genomic best linear unbiased prediction (GBLUP) with bootstrap aggregated sampling GBLUP (Bagged GBLUP, or BGBLUP) in terms of prediction accuracy. We used a 600 K Affymetrix platform with 1351 birds genotyped and phenotyped for three traits in broiler chickens; body weight, ultrasound measurement of breast muscle and hen house egg production. The predictive performance of GBLUP versus BGBLUP was evaluated in different scenarios consisting of including or excluding the TOP 20 markers from a standard genome‐wide association study (GWAS) as fixed effects in the GBLUP model, and varying training sample sizes and allelic frequency bins. Predictive performance was assessed via five replications of a threefold cross‐validation using the correlation between observed and predicted values, and prediction mean‐squared error. GBLUP overfitted the training set data, and BGBLUP delivered a better predictive ability in testing sets. Treating the TOP 20 markers from the GWAS into the model as fixed effects improved prediction accuracy and added advantages to BGBLUP over GBLUP. The performance of GBLUP and BGBLUP at different allele frequency bins and training sample sizes was similar. In general, results of this study confirm that BGBLUP can be valuable for enhancing genome‐enabled prediction of complex traits.     

An association study using imputed whole‐genome sequence data identifies novel significant loci for growth‐related traits in a Duroc × Erhualian F2 population

The average daily gain (ADG) and body weight (BW) are very important traits for breeding programs and for the meat production industry, which have attracted many researchers to delineate the genetic architecture behind these traits. In the present study, single‐ and multi‐trait genome‐wide association studies (GWAS) were performed between imputed whole‐genome sequence data and the traits of the ADG and BW at different stages in a large‐scale White Duroc × Erhualian F₂ population. A bioinformatics annotation analysis was used to assist in the identification of candidate genes that are associated with these traits. Five and seven genome‐wide significant quantitative trait loci (QTLs) were identified by single‐ and multi‐trait GWAS, respectively. Furthermore, more than 40 genome‐wide suggestive loci were detected. On the basis of the whole‐genome sequence association study and the bioinformatics analysis, NDUFAF6, TNS1 and HMGA1 stood out as the strongest candidate genes. The presented single‐ and multi‐trait GWAS analysis using imputed whole‐genome sequence data identified several novel QTLs for pig growth‐related traits. Integrating the GWAS with bioinformatics analysis can facilitate the more accurate identification of candidate genes. Higher imputation accuracy, time‐saving algorithms, improved models and comprehensive databases will accelerate the identification of causal genes or mutations, which will contribute to genomic selection and pig breeding in the future.     

Estimation of variance and genomic prediction using genotypes imputed from low‐density marker subsets for carcass traits in Japanese black cattle

The influence of genotype imputation using low‐density single nucleotide polymorphism (SNP) marker subsets on the genomic relationship matrix (G matrix), genetic variance explained, and genomic prediction (GP) was investigated for carcass weight and marbling score in Japanese Black fattened steers, using genotype data of approximately 40,000 SNPs. Genotypes were imputed using equally spaced SNP subsets of different densities. Two different linear models were used. The first (model 1) incorporated one G matrix, while the second (model 2) used two different G matrices constructed using the selected and remaining SNPs. When using model 1, the estimated additive genetic variance was always larger when using all SNPs obtained via genotype imputation than when using only equally spaced SNP subsets. The correlations between the genomic estimated breeding values obtained using genotype imputation with at least 3,000 SNPs and those using all available SNPs without imputation were higher than 0.99 for both traits. While additive genetic variance was likely to be partitioned with model 2, it did not enhance the accuracy of GP compared with model 1. These results indicate that genotype imputation using an equally spaced low‐density panel of an appropriate size can be used to produce a cost‐effective, valid GP.     

Increasing imputation and prediction accuracy for Chinese Holsteins using joint Chinese–Nordic reference population

This study investigated the effect of including Nordic Holsteins in the reference population on the imputation accuracy and prediction accuracy for Chinese Holsteins. The data used in this study include 85 Chinese Holstein bulls genotyped with both 54K chip and 777K (HD) chip, 2862 Chinese cows genotyped with 54K chip, 510 Nordic Holstein bulls genotyped with HD chip, and 4398 Nordic Holstein bulls genotyped with 54K chip and with deregressed proofs for five milk production traits. Based on these data, the accuracy of imputation from 54K to HD marker data and the accuracy of genomic predictions in Chinese Holstein were assessed. The allele correct rate increased around 2.7 and 1.7% in imputation from the 54K to the HD marker data for Chinese Holstein bulls and cows, respectively, when the Nordic HD‐genotyped bulls were included in the reference data for imputation. However, the prediction accuracy was improved slightly when using the marker data imputed based on the combined HD reference data, compared with using the marker data imputed based on the Chinese HD reference data only. On the other hand, when using the combined reference population including 4398 Nordic Holstein bulls, the accuracy of genomic predictions increased 6.5 percentage points together with a reduction of prediction bias. The HD markers did not outperform the 54K markers in genomic prediction based on the present data. The results indicate that for Chinese Holsteins, it is necessary to genotype more individuals with 54K chip to increase reference population rather than increasing marker density.     

Accuracy of genomic prediction using different models and response variables in the Nordic Red cattle population

Breeding animals can be accurately evaluated using appropriate genomic prediction models, based on marker data and phenotype information. In this study, direct genomic values (DGV) were estimated for 16 traits of Nordic Total Merit (NTM) Index in Nordic Red cattle population using three models and two different response variables. The three models were as follows: a linear mixed model (GBLUP), a Bayesian variable selection model similar to BayesA (BayesA*) and a Bayesian least absolute shrinkage and selection operator model (Bayesian Lasso). The response variables were deregressed proofs (DRP) and conventional estimated breeding values (EBV). The reliability of genomic predictions was measured on bulls in the validation data set as the squared correlation between DGV and DRP divided by the reliability of DRP. Using DRP as response variable, the reliabilities of DGV among the 16 traits ranged from 0.151 to 0.569 (average 0.317) for GBLUP, from 0.152 to 0.576 (average 0.318) for BayesA* and from 0.150 to 0.570 (average 0.320) for Bayesian Lasso. Using EBV as response variable, the reliabilities ranged from 0.159 to 0.580 (average 0.322) for GBLUP, from 0.157 to 0.578 (average 0.319) for BayesA* and from 0.159 to 0.582 (average 0.325) for Bayesian Lasso. In summary, Bayesian Lasso performed slightly better than the other two models, and EBV performed slightly better than DRP as response variable, with regard to prediction reliability of DGV. However, these differences were not statistically significant. Moreover, using EBV as response variable would result in problems with the scale of the resulting DGV and potential problem due to double counting.     

Detection of candidate polymorphisms around the QTL for fat area ratio to rib eye area on BTA7 using whole‐genome resequencing in Japanese Black cattle

In our previous study, we performed genome‐wide association study (GWAS) to identify the genomic region associated with Fat area ratio to rib eye area (FAR) and detected a candidate in BTA7 at 10–30 Mbp. The present study aims to comprehensively detect all polymorphisms in the candidate region using whole‐genome resequencing data. Based on whole‐genome resequencing of eight animals, we detected 127,090 polymorphisms within the region. Of these, 31,945 were located within the genes. We further narrowed the polymorphisms to 6,044 with more than five allele differences between the high and low FAR groups that were located within 179 genes. We subsequently investigated the functions of these genes and selected 170 polymorphisms in eight genes as possible candidate polymorphisms. We focused on SLC27A6 K81M as a putative candidate polymorphism. We genotyped the SNP in a Japanese Black population (n = 904) to investigate the effect on FAR. Analysis of variance revealed that SLC27A6 K81M had a lower p‐value (p = .0009) than the most significant SNP in GWAS (p = .0049). Although only SLC27A6 K81M was verified in the present study, subsequent verification of the remaining candidate genes and polymorphisms could lead to the identification of genes and polymorphisms responsible for FAR.     

Effect of marker‐data editing on the accuracy of genomic prediction

Genomic selection is a method to predict breeding values using genome‐wide single‐nucleotide polymorphism (SNP) markers. High‐quality marker data are necessary for genomic selection. The aim of this study was to investigate the effect of marker‐editing criteria on the accuracy of genomic predictions in the Nordic Holstein and Jersey populations. Data included 4429 Holstein and 1071 Jersey bulls. In total, 48 222 SNP for Holstein and 44 305 SNP for Jersey were polymorphic. The SNP data were edited based on (i) minor allele frequencies (MAF) with thresholds of no limit, 0.001, 0.01, 0.02, 0.05 and 0.10, (ii) deviations from Hardy–Weinberg proportions (HWP) with thresholds of no limit, chi‐squared p‐values of 0.001, 0.02, 0.05 and 0.10, and (iii) GenCall (GC) scores with thresholds of 0.15, 0.55, 0.60, 0.65 and 0.70. The marker data sets edited with different criteria were used for genomic prediction of protein yield, fertility and mastitis using a Bayesian variable selection and a GBLUP model. De‐regressed EBV were used as response variables. The result showed little difference between prediction accuracies based on marker data sets edited with MAF and deviation from HWP. However, accuracy decreased with more stringent thresholds of GC score. According to the results of this study, it would be appropriate to edit data with restriction of MAF being between 0.01 and 0.02, a p‐value of deviation from HWP being 0.05, and keeping all individual SNP genotypes having a GC score over 0.15.     

Whole genome sequence analysis to detect signatures of positive selection for high fecundity in sheep

Ovulation rate and prolificacy are the most important reproductive traits that have major impact on the efficiency of lamb meat production. Here, we compared the whole genomes of the Romanov sheep, known as one of the high prolific breeds, and four other sheep breeds namely Assaf, Awassi, Cambridge and British du cher, to identify genetic mechanisms underlying prolificacy in sheep. Selection signature analysis revealed 637 and 477 protein‐coding genes under positive selection from F_ST and nucleotide diversity (Pi) statistics, respectively. Further analysis showed that several candidate genes including LEPR, PDGFRL and KLF5 genes are involved in sheep prolificacy. The identified candidate genes in the selected regions are novel and provide new insights into the genetic mechanisms underlying prolificacy in sheep and can be useful in sheep breeding programmes to develop improved breeds for high reproductive efficiency.     

Genomic prediction in a numerically small breed population using prioritized genetic markers from whole-genome sequence data

The objective of this study was to investigate the accuracy of genomic prediction of body weight and eating quality traits in a numerically small sheep population (Dorper sheep). Prediction was based on a large multi-breed/admixed reference population and using (a) 50k or 500k single nucleotide polymorphism (SNP) genotypes, (b) imputed whole-genome sequencing data (~31 million), (c) selected SNPs from whole genome sequence data and (d) 50k SNP genotypes plus selected SNPs from whole-genome sequence data. Furthermore, the impact of using a breed-adjusted genomic relationship matrix on accuracy of genomic breeding value was assessed. The selection of genetic variants was based on an association study performed on imputed whole-genome sequence data in an independent population, which was chosen either randomly from the base population or according to higher genetic proximity to the target population. Genomic prediction was based on genomic best linear unbiased prediction (GBLUP), and the accuracy of genomic prediction was assessed according to the correlation between genomic breeding value and corrected phenotypes divided by the square root of trait heritability. The accuracy of genomic prediction was between 0.20 and 0.30 across different traits based on common 50k SNP genotypes, which improved on average by 0.06 (absolute value) on average based on using prioritized genetic markers from whole-genome sequence data. Using prioritized genetic markers from a genetically more related GWAS population resulted in slightly higher prediction accuracy (0.02 absolute value) compared to genetic markers derived from a random GWAS population. Using high-density SNP genotypes or imputed whole-genome sequence data in GBLUP showed almost no improvement in genomic prediction accuracy however, accounting for different marker allele frequencies in reference population according to a breed-adjusted GRM resulted to on average 0.024 (absolute value) increase in accuracy of genomic prediction.     

Accuracy of genomic prediction of purebreds for cross bred performance in pigs

In pig breeding, as the final product is a cross bred (CB) animal, the goal is to increase the CB performance. This goal requires different strategies for the implementation of genomic selection from what is currently implemented in, for example dairy cattle breeding. A good strategy is to estimate marker effects on the basis of CB performance and subsequently use them to select pure bred (PB) breeding animals. The objective of our study was to assess empirically the predictive ability (accuracy) of direct genomic values of PB for CB performance across two traits using CB and PB genomic and phenotypic data. We studied three scenarios in which genetic merit was predicted within each population, and four scenarios where PB genetic merit for CB performance was predicted based on either CB or a PB training data. Accuracy of prediction of PB genetic merit for CB performance based on CB training data ranged from 0.23 to 0.27 for gestation length (GLE), whereas it ranged from 0.11 to 0.22 for total number of piglets born (TNB). When based on PB training data, it ranged from 0.35 to 0.55 for GLE and from 0.30 to 0.40 for TNB. Our results showed that it is possible to predict PB genetic merit for CB performance using CB training data, but predictive ability was lower than training using PB training data. This result is mainly due to the structure of our data, which had small‐to‐moderate size of the CB training data set, low relationship between the CB training and the PB validation populations, and a high genetic correlation (0.94 for GLE and 0.90 for TNB) between the studied traits in PB and CB individuals, thus favouring selection on the basis of PB data.     

Accuracy of imputation of single nucleotide polymorphism marker genotypes from low‐density panels in Japanese Black cattle

Using target and reference fattened steer populations, the performance of genotype imputation using lower‐density marker panels in Japanese Black cattle was evaluated. Population imputation was performed using BEAGLE software. Genotype information for approximately 40 000 single nucleotide polymorphism (SNP) markers by Illumina BovineSNP50 BeadChip was available, and imputation accuracy was assessed based on the average concordance rates of the genotypes, varying equally spaced SNP densities, and the number of individuals in the reference population. Two additional statistics were also calculated as indicators of imputation performance. The concordance rates tended to be lower for SNPs with greater minor allele frequencies, or those located near the ends of the chromosomes. Longer autosomes yielded greater imputation accuracies than shorter ones. When SNPs were selected based on linkage disequilibrium information, relative imputation accuracy was slightly improved. When 3000 and 10 000 equally spaced SNPs were used, the imputation accuracies were greater than 90% and approximately 97%, respectively. These results indicate that combining genotyping using a lower‐density SNP chip with genotype imputation based on a population of individuals genotyped using a higher‐density SNP chip is a cost‐effective and valid approach for genomic prediction.     

Molecular characterization of community‐associated methicillin‐resistant Staphylococcus aureus from pet dogs

Community‐associated methicillin‐resistant Staphylococcus aureus (MRSA) is a serious public health concern and in Australia, one that disproportionately affects Aboriginal people. Paralleling MRSA in human medicine, methicillin‐resistant S. pseudintermedius (MRSP) is an increasingly prevalent pathogen in veterinary medicine. We aimed to characterize the carriage of MRSA and MRSP in dogs and cats from predominantly Aboriginal communities in a very remote region of New South Wales (NSW), Australia. Pets (303 dogs and 80 cats) were recruited from six communities in western NSW. Three swabs were collected from each animal (anterior nares, oropharynx and perineum) and from skin lesions or wounds (if present) and cultured on selective media for methicillin‐resistant staphylococci. Human host‐adapted community‐associated MRSA representing four multilocus sequence types (ST1‐IV, ST5‐IV, ST72‐IV, ST93‐IV) were isolated from eight dogs (prevalence 2.6%, 95% confidence interval 1.3%–5.1%). Two ST5‐IV isolates from a single dog were phenotypically trimethoprim‐resistant, harbouring trimethoprim‐resistant gene dfrG within the SCCmec type IVo mobile genetic element. MRSA was not isolated from any cats and MRSP was not isolated from any dogs or cats. This study estimated a high prevalence of human host‐adapted community‐associated MRSA carriage in dogs despite an absence of MRSP. This suggests MRSA carried by dogs in remote NSW originate from human hosts. The cycle of transmission between people, dogs and common environmental sources warrants further investigation. To our knowledge, this is the first report of trimethoprim‐resistant ST5‐IV in eastern Australia and the first report of trimethoprim‐resistant ST5‐IV from a dog.     

A comparison of accuracy validation methods for genomic and pedigree‐based predictions of swine litter size traits using Large White and simulated data

The objective of this study was to compare and determine the optimal validation method when comparing accuracy from single‐step GBLUP (ssGBLUP) to traditional pedigree‐based BLUP. Field data included six litter size traits. Simulated data included ten replicates designed to mimic the field data in order to determine the method that was closest to the true accuracy. Data were split into training and validation sets. The methods used were as follows: (i) theoretical accuracy derived from the prediction error variance (PEV) of the direct inverse (iLHS), (ii) approximated accuracies from the accf90(GS) program in the BLUPF90 family of programs (Approx), (iii) correlation between predictions and the single‐step GEBVs from the full data set (GEBV_Full), (iv) correlation between predictions and the corrected phenotypes of females from the full data set (Y_c), (v) correlation from method iv divided by the square root of the heritability (Y_ch) and (vi) correlation between sire predictions and the average of their daughters' corrected phenotypes (Y_cs). Accuracies from iLHS increased from 0.27 to 0.37 (37%) in the Large White. Approximation accuracies were very consistent and close in absolute value (0.41 to 0.43). Both iLHS and Approx were much less variable than the corrected phenotype methods (ranging from 0.04 to 0.27). On average, simulated data showed an increase in accuracy from 0.34 to 0.44 (29%) using ssGBLUP. Both iLHS and Y_ch approximated the increase well, 0.30 to 0.46 and 0.36 to 0.45, respectively. GEBV_Full performed poorly in both data sets and is not recommended. Results suggest that for within‐breed selection, theoretical accuracy using PEV was consistent and accurate. When direct inversion is infeasible to get the PEV, correlating predictions to the corrected phenotypes divided by the square root of heritability is adequate given a large enough validation data set.     

Gaussian covariance graph models accounting for correlated marker effects in genome‐wide prediction

Several statistical models used in genome‐wide prediction assume uncorrelated marker allele substitution effects, but it is known that these effects may be correlated. In statistics, graphical models have been identified as a useful tool for covariance estimation in high‐dimensional problems and it is an area that has recently experienced a great expansion. In Gaussian covariance graph models (GCovGM), the joint distribution of a set of random variables is assumed to be Gaussian and the pattern of zeros of the covariance matrix is encoded in terms of an undirected graph G. In this study, methods adapting the theory of GCovGM to genome‐wide prediction were developed (Bayes GCov, Bayes GCov‐KR and Bayes GCov‐H). In simulated data sets, improvements in correlation between phenotypes and predicted breeding values and accuracies of predicted breeding values were found. Our models account for correlation of marker effects and permit to accommodate general structures as opposed to models proposed in previous studies, which consider spatial correlation only. In addition, they allow incorporation of biological information in the prediction process through its use when constructing graph G, and their extension to the multi‐allelic loci case is straightforward.     

Whole‐genome resequencing to identify candidate genes for the QTL for oleic acid percentage in Japanese Black cattle

In our previous study, we detected a QTL for the oleic acid percentage (C18:1) on BTA9 in Japanese Black cattle through a genome‐wide association study (GWAS). In this study, we performed whole‐genome resequencing on eight animals with higher and lower C18:1 to identify candidate polymorphisms for the QTL. A total of 39,658 polymorphisms were detected in the candidate region, which were narrowed to 1993 polymorphisms within 23 genes based on allele differences between the high and low C18:1 groups. We subsequently selected three candidate genes, that is, CYB5R4, MED23, and VNN1, among the 23 genes based on their function in fatty acid metabolism. In each candidate gene, three SNPs, that is, CYB5R4 c.*349G > T, MED23 c.3700G > A, and VNN1 c.197C > T, were selected as candidate SNPs to verify their effect on C18:1 in a Japanese Black cattle population (n = 889). The statistical analysis showed that these SNPs were significantly associated with C18:1 (p < 0.05), suggesting that they were candidates for the QTL. In conclusion, we successfully narrowed the candidates for the QTL by detecting possible polymorphisms located within the candidate region. It is expected that the responsible polymorphism can be identified by demonstrating their effect on the gene's function.     

Accounting for genetic architecture in single‐ and multipopulation genomic prediction using weights from genomewide association studies in pigs

We studied the effect of including GWAS results on the accuracy of single‐ and multipopulation genomic predictions. Phenotypes (backfat thickness) and genotypes of animals from two sire lines (SL1, n = 1146 and SL3, n = 1264) were used in the analyses. First, GWAS were conducted for each line and for a combined data set (both lines together) to estimate the genetic variance explained by each SNP. These estimates were used to build matrices of weights (D), which was incorporated into a GBLUP method. Single population evaluated with traditional GBLUP had accuracies of 0.30 for SL1 and 0.31 for SL3. When weights were employed in GBLUP, the accuracies for both lines increased (0.32 for SL1 and 0.34 for SL3). When a multipopulation reference set was used in GBLUP, the accuracies were higher (0.36 for SL1 and 0.32 for SL3) than in single‐population prediction. In addition, putting together the multipopulation reference set and the weights from the combined GWAS provided even higher accuracies (0.37 for SL1, and 0.34 for SL3). The use of multipopulation predictions and weights estimated from a combined GWAS increased the accuracy of genomic predictions.     

Additional considerations to the use of single‐step genomic predictions in a dominance setting

Recent publications indicate that single‐step models are suitable to estimate breeding values, dominance deviations and total genetic values with acceptable quality. Additive single‐step methods implicitly extend known number of allele information from genotyped to non‐genotyped animals. This theory is well derived in an additive setting. It was recently shown, at least empirically, that this basic strategy can be extended to dominance with reasonable prediction quality. Our study addressed two additional issues. It illustrated the theoretical basis for extension and validated genomic predictions to dominance based on single‐step genomic best linear unbiased prediction theory. This development was then extended to include inbreeding into dominance relationships, which is a currently not yet solved issue. Different parametrizations of dominance relationship matrices were proposed. Five dominance single‐step inverse matrices were tested and described as C¹ , C² , C³ , C⁴ and C⁵ . Genotypes were simulated for a real pig population (n = 11,943 animals). In order to avoid any confounding issues with additive effects, pseudo‐records including only dominance deviations and residuals were simulated. SNP effects of heterozygous genotypes were summed up to generate true dominance deviations. We added random noise to those values and used them as phenotypes. Accuracy was defined as correlation between true and predicted dominance deviations. We conducted five replicates and estimated accuracies in three sets: between all ( S₁ ), non‐genotyped ( S₂ ) and inbred non‐genotyped ( S₃ ) animals. Potential bias was assessed by regressing true dominance deviations on predicted values. Matrices accounting for inbreeding ( C³ , C⁴ and C⁵ ) best fit. Accuracies were on average 0.77, 0.40 and 0.46 in S₁ , S₂ and S₃ , respectively. In addition, C³ , C⁴ and C⁵ scenarios have shown better accuracies than C¹ and C² , and dominance deviations were less biased. Better matrix compatibility (accuracy and bias) was observed by re‐scaling diagonal elements to 1 minus the inbreeding coefficient ( C⁵ ).     

Comparison of heritabilities of dairy traits in Australian Holstein‐Friesian cattle from genomic and pedigree data and implications for genomic evaluations

The reliability of genomic evaluations depends on the proportion of genetic variation explained by the DNA markers. In this study, we have estimated the proportion of variance in daughter trait deviations (DTDs) of dairy bulls explained by 45 993 genome wide single‐nucleotide poly‐ morphism (SNP) markers for 29 traits in Australian Holstein‐Friesian dairy cattle. We compare these proportions to the proportion of variance in DTDs explained by the additive relationship matrix derived from the pedigree, as well as the sum of variance explained by both pedigree and marker information when these were fitted simultaneously. The propor‐ tion of genetic variance in DTDs relative to the total genetic variance (the total genetic variance explained by the genomic relationships and pedigree relationships when both were fitted simultaneously) varied from 32% for fertility to approximately 80% for milk yield traits. When fitting genomic and pedigree relationships simultaneously, the variance unexplained (i.e. the residual variance) in DTDs of the total variance for most traits was reduced compared to fitting either individually, suggesting that there is not complete overlap between the effects. The proportion of genetic variance accounted by the genomic relationships can be used to modify the blending equations used to calculate genomic estimated breeding value (GEBV) from direct genomic breeding value (DGV) and parent average. Our results, from a validation population of young dairy bulls with DTD, suggest that this modification can improve the reliability of GEBV by up to 5%.