Sharing reference data and including cows in the reference population improve genomic predictions in Danish Jersey

A breed-of-origin of alleles model that includes crossbred data improves predictive ability for crossbred animals in a multi-breed population

BACKGROUND

Recently, crossbred animals have begun to be used as parents in the next generations of dairy and beef cattle systems, which has increased the interest in predicting the genetic merit of those animals. The primary objective of this study was to investigate three available methods for genomic prediction of crossbred animals. In the first two methods, SNP effects from within-breed evaluations are used by weighting them by the average breed proportions across the genome (BPM method) or by their breed-of-origin (BOM method). The third method differs from the BOM in that it estimates breed-specific SNP effects using purebred and crossbred data, considering the breed-of-origin of alleles (BOA method). For within-breed evaluations, and thus for BPM and BOM, 5948 Charolais, 6771 Limousin and 7552 Others (a combined population of other breeds) were used to estimate SNP effects separately within each breed. For the BOA, the purebreds' data were enhanced with data from ~ 4K, ~ 8K or ~ 18K crossbred animals. For each animal, its predictor of genetic merit (PGM) was estimated by considering the breed-specific SNP effects. Predictive ability and absence of bias were estimated for crossbreds and the Limousin and Charolais animals. Predictive ability was measured as the correlation between PGM and the adjusted phenotype, while the regression of the adjusted phenotype on PGM was estimated as a measure of bias.

RESULTS

With BPM and BOM, the predictive abilities for crossbreds were 0.468 and 0.472, respectively, and with the BOA method, they ranged from 0.490 to 0.510. The performance of the BOA method improved as the number of crossbred animals in the reference increased and with the use of the correlated approach, in which the correlation of SNP effects across the genome of the different breeds was considered. The slopes of regression for PGM on adjusted phenotypes for crossbreds showed overdispersion of the genetic merits for all methods but this bias tended to be reduced by the use of the BOA method and by increasing the number of crossbred animals.

CONCLUSIONS

For the estimation of the genetic merit of crossbred animals, the results from this study suggest that the BOA method that accommodates crossbred data can yield more accurate predictions than the methods that use SNP effects from separate within-breed evaluations.

Genomic prediction of service sire effect on female reproductive performance in Holstein cattle: A comparison between different methods, validation population and marker densities

Reproductive traits of dairy cattle are bound to the actual efficiency of farm operation, which therefore show great economic importance. Among them, some traits were deemed to be simultaneously affected by service sire and mating cow. Service sires are proved to play an important role in reproduction process of cows. However, limited study explored the genetic effect of service sire (GESS), let alone the genomic prediction of this effect. In the present study, 2244 genotyped bulls together with phenotypic records were used to predict the GESS on conception rate, 56-day non-return rate, calving ease, stillbirth and gestation length. The feasibilities of multi-step genomic best linear unbiased predictor (msGBLUP) and single-step genomic best linear unbiased predictor (ssGBLUP) were investigated under different scenarios, that is, different marker densities and validation population. The predictive accuracies and unbiasedness for GESS ranged from 0.159 to 0.647 and from 0.202 to 2.018, respectively, when validated on young bulls, while the accuracies and unbiasedness ranged from 0.409 to 0.802 and 0.333 to 1.146 when validated on random split data sets. It is feasible to predict GESS on reproductive traits by using a linear mixed model and genomic data, and high-density marker panel had limited contribution to the prediction. This research investigated the potential factors that influence the genomic prediction of GESS on reproductive traits and indicated the possibility of genomic selection on GESS, both in ideal and practical circumstances.

Improving Genomic Prediction Accuracy in the Chinese Holstein Population by Combining with the Nordic Holstein Reference Population

The size of the reference population is critical in order to improve the accuracy of genomic prediction. Indeed, improving genomic prediction accuracy by combining multinational reference populations has proven to be effective. In this study, we investigated the improvement of genomic prediction accuracy in seven complex traits (i.e., milk yield; fat yield; protein yield; somatic cell count; body conformation; feet and legs; and mammary system conformation) by combining the Chinese and Nordic Holstein reference populations. The estimated genetic correlations between the Chinese and Nordic Holstein populations are high with respect to protein yield, fat yield, and milk yield-whereby these correlations range from 0.621 to 0.720-and are moderate with respect to somatic cell count (0.449), but low for the three conformation traits (which range from 0.144 to 0.236). When utilizing the joint reference data and a two-trait GBLUP model, the genomic prediction accuracy in the Chinese Holsteins improves considerably with respect to the traits with moderate-to-high genetic correlations, whereas the improvement in Nordic Holsteins is small. When compared with the single population analysis, using the joint reference population for genomic prediction in younger animals, results in a 2.3 to 8.1 percent improvement in accuracy. Meanwhile, 10 replications of five-fold cross-validation were also implemented in order to evaluate the performance of joint genomic prediction, thereby resulting in a 1.6 to 5.2 percent increase in accuracy. With respect to joint genomic prediction, the bias was found to be quite low. However, for traits with low genetic correlations, the joint reference data do not improve the prediction accuracy substantially for either population.

Genomic Selection for Milk Production Traits in Xinjiang Brown Cattle

One-step genomic selection is a method for improving the reliability of the breeding value estimation. This study aimed to compare the reliability of pedigree-based best linear unbiased prediction (PBLUP) and single-step genomic best linear unbiased prediction (ssGBLUP), single-trait and multitrait models, and the restricted maximum likelihood (REML) and Bayesian methods. Data were collected from the production performance records of 2207 Xinjiang Brown cattle in Xinjiang from 1983 to 2018. A cross test was designed to calculate the genetic parameters and reliability of the breeding value of 305 daily milk yield (305 dMY), milk fat yield (MFY), milk protein yield (MPY), and somatic cell score (SCS) of Xinjiang Brown cattle. The heritability of 305 dMY, MFY, MPY, and SCS estimated using the REML and Bayesian multitrait models was approximately 0.39 (0.02), 0.40 (0.03), 0.49 (0.02), and 0.07 (0.02), respectively. The heritability and estimated breeding value (EBV) and the reliability of milk production traits of these cattle calculated based on PBLUP and ssGBLUP using the multitrait model REML and Bayesian methods were higher than those of the single-trait model REML method; the ssGBLUP method was significantly better than the PBLUP method. The reliability of the estimated breeding value can be improved from 0.9% to 3.6%, and the reliability of the genomic estimated breeding value (GEBV) for the genotyped population can reach 83%. Therefore, the genetic evaluation of the multitrait model is better than that of the single-trait model. Thus, genomic selection can be applied to small population varieties such as Xinjiang Brown cattle, in improving the reliability of the genomic estimated breeding value.

Genomic Estimated Breeding Value of Milk Performance and Fertility Traits in the Russian Black-and-White Cattle Population

A breakthrough in cattle breeding was achieved with the incorporation of animal genomic data into breeding programs. The introduction of genomic selection has a major impact on traditional genetic assessment systems and animal genetic improvement programs. Since 2010, genomic selection has been officially introduced in the evaluation of the breeding and genetic potential of cattle in Europe, the U.S., Canada, and many other developed countries. The purpose of this study is to develop a system for a genomic evaluation of the breeding value of the domestic livestock of Black-and-White and Russian Holstein cattle based on 3 milk performance traits: daily milk yield (kg), daily milk fat (%), and daily milk protein content (%) and 6 fertility traits: age at first calving (AFC), calving interval (CI), calving to first insemination interval (CFI), interval between first and last insemination (IFL), days open (DO), and number of services (NS). We built a unified database of breeding animals from 523 breeding farms in the Russian Federation. The database included pedigree information on 2,551,529 cows and 69,131 bulls of the Russian Holstein and Black-and-White cattle breeds, as well as information on the milk performance of 1,597,426 cows with 4,771,366 completed lactations. The date of birth of the animals included in the database was between 1975 and 2017. Genotyping was performed in 672 animals using a BovineSNP50 v3 DNA Analysis BeadChip microarray (Illumina, USA). The genomic estimated breeding value (GEBV) was evaluated only for 644 animals (427 bulls and 217 cows) using the single-step genomic best linear unbiased prediction - animal model (ssGBLUP-AM). The mean genetic potential was +0.88 and +1.03 kg for the daily milk yield, -0.002% for the milk fat content, and -0.003 and 0.001% for the milk protein content in the cows and bulls, respectively. There was negative genetic progress in the fertility traits in the studied population between 1975 and 2017. The reliability of the estimated breeding value (EBV) for genotyped bulls ranged from 89 to 93% for the milk performance traits and 85 to 90% for the fertility traits, whereas the reliability of the genomic estimated breeding value (GEBV) varied 54 to 64% for the milk traits and 23 to 60% for the fertility traits. This result shows that it is possible to use the genomic estimated breeding value with rather high reliability to evaluate the domestic livestock of Russian Holstein and Black-and-White cattle breeds for fertility and milk performance traits. This system of genomic evaluation may help bring domestic breeding in line with modern competitive practices and estimate the breeding value of cattle at birth based on information on the animal's genome.

Across-country genomic prediction of bull fertility in Jersey dairy cattle

The use of information across populations is an attractive approach to increase the accuracy of genomic predictions for numerically small breeds and traits that are time-consuming and difficult to measure, such as male fertility in cattle. This study was conducted to evaluate genomic prediction of Jersey bull fertility using an across-country reference population combining records from the United States and Australia. The data set consisted of 1,570 US Jersey bulls with sire conception rate (SCR) records, 603 Australian Jersey bulls with semen fertility value (SFV) records and SNP genotypes for roughly 90,000 loci. Both SCR and SFV are evaluations of service sire fertility based on cow field data, and both are intended as phenotypic evaluations because the estimates include genetic and nongenetic effects. Within- and across-country genomic predictions were evaluated using univariate and bivariate genomic best linear unbiased prediction models. Predictive ability was assessed in 5-fold cross-validation using the correlation between observed and predicted fertility values and mean squared error of prediction. Within-country genomic predictions exhibited predictive correlations of around 0.28 and 0.02 for the United States and Australia, respectively. The Australian Jersey population is genetically diverse and small in size, so careful selection of the reference population by including only closely related animals (e.g., excluding New Zealand bulls, which is a less-related population) increased the predictive correlations up to 0.20. Notably, the use of bivariate models fitting all US Jersey records and the optimized Australian population resulted in predictive correlations around of 0.24 for SFV values, which is a relative increase in predictive ability of 20%. Conversely, for predicting SCR values, the use of an across-country reference population did not outperform the standard approach using pure US Jersey reference data set. Our findings indicate that genomic prediction of male fertility in dairy cattle is feasible, and the use of an across-country reference population would be beneficial when local populations are small and genetically diverse.

Maternal effect on body measurement and meat production traits in purebred Duroc pigs

We investigated maternal effect on nine body measurement traits (body height, body length, front width (FW), chest width (CW), hind width (HW), chest depth, chest girth (CHG), front cannon circumference (FCC) and rear cannon circumference (RCC)) measured at the end of performance testing and five meat production traits (ages at the start and end of performance testing (D30 and D105), average daily gain (ADG), backfat thickness and loin muscle area) in purebred Duroc pigs. Genetic parameters for each trait were estimated by using six single-trait models with and without common litter environmental effect, maternal genetic effect and direct-maternal genetic correlation. The value of Akaike's information criterion was lowest with the model including direct additive genetic and common litter environmental effects for 10 traits. The estimated proportion of common litter environmental variance to phenotypic variance was approximately ≥0.1 for D30, D105, ADG, FW, CW, HW, CHG, FCC and RCC. Using a model without common litter environmental effect would overestimate the direct heritability of most traits. Standard errors of estimated genetic parameters tended to be larger in models including maternal genetic effect. The results indicate that a compromise could be made for accurate genetic parameter estimation for body measurement traits, as well as meat production traits, in pigs by considering common litter environmental effect.

Imputation for sequencing variants preselected to a customized low-density chip

The sequencing variants preselected from association analyses and bioinformatics analyses could improve genomic prediction. In this study, the imputation of sequencing SNPs preselected from major dairy breeds in Denmark-Finland-Sweden (DFS) and France (FRA) was investigated for both contemporary animals and old bulls in Danish Jersey. For contemporary animals, a two-step imputation which first imputed to 54 K and then to 54 K + DFS + FRA SNPs achieved highest accuracy. Correlations between observed and imputed genotypes were 91.6% for DFS SNPs and 87.6% for FRA SNPs, while concordance rates were 96.6% for DFS SNPs and 93.5% for FRA SNPs. The SNPs with lower minor allele frequency (MAF) tended to have lower correlations but higher concordance rates. For old bulls, imputation for DFS and FRA SNPs were relatively accurate even for bulls without progenies (correlations higher than 97.2% and concordance rates higher than 98.4%). For contemporary animals, given limited imputation accuracy of preselected sequencing SNPs especially for SNPs with low MAF, it would be a good strategy to directly genotype preselected sequencing SNPs with a customized SNP chip. For old bulls, given high imputation accuracy for preselected sequencing SNPs with all MAF ranges, it would be unnecessary to re-genotype preselected sequencing SNPs.

Symposium review: How to implement genomic selection

Genomic selection was adopted very quickly in the 10 yr after first implementation, and breeders continue to find new uses for genomic testing. Breeding values with higher reliability earlier in life are estimated by combining DNA genotypes for many thousands of loci using existing identification, pedigree, and phenotype databases for millions of animals. Quality control for both new and previous data is greatly improved by comparing genomic and pedigree relationships to correct parent-progeny conflicts and discover many additional ancestors. Many quantitative trait loci and gene tests have been added to previous assays that used only evenly spaced, highly polymorphic markers. Imputation now combines genotypes from many assays of differing marker densities. Prediction models have gradually advanced from normal or Bayesian distributions within trait and breed to single-step, multitrait, or other more complex models, such as multibreed models that may be needed for crossbred prediction. Genomic selection was initially applied to males to predict progeny performance but is now widely applied to females or even embryos to predict their own later performance. The initial focus on additive merit has expanded to include mating programs, genomic inbreeding, and recessive alleles. Many producers now use DNA testing to decide which heifers should be inseminated with elite dairy, beef, or sex-sorted semen, which should be embryo donors or recipients, or which should be sold or kept for breeding. Because some of these decisions are expensive to delay, predictions are now provided weekly instead of every few months. Predictions from international genomic databases are often more accurate and cost-effective than those from within-country databases that were previously designed for progeny testing unless local breeds, conditions, or traits differ greatly from the larger database. Selection indexes include many new traits, often with lower heritability or requiring large initial investments to obtain phenotypes, which provide further incentive to cooperate internationally. The genomic prediction methods developed for dairy cattle are now applied widely to many animal, human, and plant populations and could be applied to many more.

High-Throughput Sequencing With the Preselection of Markers Is a Good Alternative to SNP Chips for Genomic Prediction in Broilers

The choice of a genetic marker genotyping platform is important for genomic prediction in livestock and poultry. High-throughput sequencing can produce more genetic markers, but the genotype quality is lower than that obtained with single nucleotide polymorphism (SNP) chips. The aim of this study was to compare the accuracy of genomic prediction between high-throughput sequencing and SNP chips in broilers. In this study, we developed a new SNP marker screening method, the pre-marker-selection (PMS) method, to determine whether an SNP marker can be used for genomic prediction. We also compared a method which preselection marker based results from genome-wide association studies (GWAS). With the two methods, we analysed body weight at the12^th week (BW) and feed conversion ratio (FCR) in a local broiler population. A total of 395 birds were selected from the F2 generation of the population, and 10X specific-locus amplified fragment sequencing (SLAF-seq) and the Illumina Chicken 60K SNP Beadchip were used for genotyping. The genomic best linear unbiased prediction method (GBLUP) was used to predict the genomic breeding values. The accuracy of genomic prediction was validated by the leave-one-out cross-validation method. Without SNP marker screening, the accuracies of the genomic estimated breeding value (GEBV) of BW and FCR were 0.509 and 0.249, respectively, when using SLAF-seq, and the accuracies were 0.516 and 0.232, respectively, when using the SNP chip. With SNP marker screening by the PMS method, the accuracies of GEBV of the two traits were 0.671 and 0.499, respectively, when using SLAF-seq, and 0.605 and 0.422, respectively, when using the SNP chip. Our SNP marker screening method led to an increase of prediction accuracy by 0.089-0.250. With SNP marker screening by the GWAS method, the accuracies of genomic prediction for the two traits were also improved, but the gains of accuracy were less than the gains with PMS method for all traits. The results from this study indicate that our PMS method can improve the accuracy of GEBV, and that more accurate genomic prediction can be obtained from an increased number of genomic markers when using high-throughput sequencing in local broiler populations. Due to its lower genotyping cost, high-throughput sequencing could be a good alternative to SNP chips for genomic prediction in breeding programmes of local broiler populations.

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.

Improvement of genomic prediction by integrating additional single nucleotide polymorphisms selected from imputed whole genome sequencing data

The availability of whole genome sequencing (WGS) data enables the discovery of causative single nucleotide polymorphisms (SNPs) or SNPs in high linkage disequilibrium with causative SNPs. This study investigated effects of integrating SNPs selected from imputed WGS data into the data of 54K chip on genomic prediction in Danish Jersey. The WGS SNPs, mainly including peaks of quantitative trait loci, structure variants, regulatory regions of genes, and SNPs within genes with strong effects predicted with variant effect predictor, were selected in previous analyses for dairy breeds in Denmark–Finland–Sweden (DFS) and France (FRA). Animals genotyped with 54K chip, standard LD chip, and customized LD chip which covered selected WGS SNPs and SNPs in the standard LD chip, were imputed to 54K together with DFS and FRA SNPs. Genomic best linear unbiased prediction (GBLUP) and Bayesian four-distribution mixture models considering 54K and selected WGS SNPs as one (a one-component model) or two separate genetic components (a two-component model) were used to predict breeding values. For milk production traits and mastitis, both DFS (0.025) and FRA (0.029) sets of additional WGS SNPs improved reliabilities, and inclusions of all selected WGS SNPs generally achieved highest improvements of reliabilities (0.034). A Bayesian four-distribution model yielded higher reliabilities than a GBLUP model for milk and protein, but extra gains in reliabilities from using selected WGS SNPs were smaller for a Bayesian four-distribution model than a GBLUP model. Generally, no significant difference was observed between one-component and two-component models, except for using GBLUP models for milk.

Effects of preselection of genotyped animals on reliability and bias of genomic prediction in dairy cattle

OBJECTIVE

Models for genomic selection assume that the reference population is an unselected population. However, in practice, genotyped individuals, such as progeny-tested bulls, are highly selected, and the reference population is created after preselection. In dairy cattle, the intensity of selection is higher in males than in females, suggesting that cows can be added to the reference population with less bias and loss of accuracy. The objective is to develop formulas applied to any genomic prediction studies or practice with preselected animals as reference population.

METHODS

We developed formulas for calculating the reliability and bias of genomically enhanced breeding values (GEBV) in the reference population where individuals are preselected on estimated breeding values. Based on the formulas presented, deterministic simulation was conducted by varying heritability, preselection percentage, and the reference population size.

RESULTS

The number of bulls equal to a cow regarding the reliability of GEBV was expressed through a simple formula for the reference population consisting of preselected animals. The bull population was vastly superior to the cow population regarding the reliability of GEBV for low-heritability traits. However, the superiority of reliability from the bull reference population over the cow population decreased as heritability increased. Bias was greater for bulls than cows. Bias and reduction in reliability of GEBV due to preselection was alleviated by expanding reference population.

CONCLUSION

Cows are easier in expanding reference population size compared with bulls and alleviate bias and reduction in reliability of GEBV of bulls which are highly preselected than cows by expanding the cow reference population.

Genomic Prediction Using Multi-trait Weighted GBLUP Accounting for Heterogeneous Variances and Covariances Across the Genome

Implicit assumption of common (co)variance for all loci in multi-trait Genomic Best Linear Unbiased Prediction (GBLUP) results in a genomic relationship matrix (G) that is common to all traits. When this assumption is violated, Bayesian whole genome regression methods may be superior to GBLUP by accounting for unequal (co)variance for all loci or genome regions. This study aimed to develop a strategy to improve the accuracy of GBLUP for multi-trait genomic prediction, using (co)variance estimates of SNP effects from Bayesian whole genome regression methods. Five generations (G1-G5, test populations) of genotype data were available by simulations based on data of 2,200 Danish Holstein cows (G0, reference population). Two correlated traits with heritabilities of 0.1 or 0.4, and a genetic correlation of 0.45 were generated. First, SNP effects and breeding values were estimated using BayesAS method, assuming (co)variance was the same for SNPs within a genome region, and different between regions. Region size was set as one SNP, 100 SNPs, a whole chromosome or whole genome. Second, posterior (co)variances of SNP effects were used to weight SNPs in construction of G matrices. In general, region size of 100 SNPs led to highest prediction accuracies using BayesAS, and wGBLUP outperformed GBLUP at this region size. Our results suggest that when genetic architectures of traits favor Bayesian methods, the accuracy of multi-trait GBLUP can be as high as the Bayesian method if SNPs are weighted by the Bayesian posterior (co)variances.

Using a very low-density SNP panel for genomic selection in a breeding program for sheep

Background

Building an efficient reference population for genomic selection is an issue when the recorded population is small and phenotypes are poorly informed, which is often the case in sheep breeding programs. Using stochastic simulation, we evaluated a genomic design based on a reference population with medium-density genotypes [around 45 K single nucleotide polymorphisms (SNPs)] of dams that were imputed from very low-density genotypes (≤ 1000 SNPs).

Methods

A population under selection for a maternal trait was simulated using real genotypes. Genetic gains realized from classical selection and genomic selection designs were compared. Genomic selection scenarios that differed in reference population structure (whether or not dams were included in the reference) and genotype quality (medium-density or imputed to medium-density from very low-density) were evaluated.

Results

The genomic design increased genetic gain by 26% when the reference population was based on sire medium-density genotypes and by 54% when the reference population included both sire and dam medium-density genotypes. When medium-density genotypes of male candidates and dams were replaced by imputed genotypes from very low-density SNP genotypes (1000 SNPs), the increase in gain was 22% for the sire reference population and 42% for the sire and dam reference population. The rate of increase in inbreeding was lower (from − 20 to − 34%) for the genomic design than for the classical design regardless of the genomic scenario.

Conclusions

We show that very low-density genotypes of male candidates and dams combined with an imputation process result in a substantial increase in genetic gain for small sheep breeding programs.

Genomic evaluation using SNP- and haplotype-based genomic relationship matrices in a closed line of Duroc pigs

A simulation analysis and real phenotype analysis were performed to evaluate the impact of three different relationship matrices on heritability estimation and prediction accuracy in a closed-line breeding of Duroc pigs. The numerator relationship matrix (NRM), single nucleotide polymorphism (SNP)-based genomic relationship matrix (GRM) (G_S ), and haplotype-based GRM (G_H ) were applied in this study. We used PorcineSNP60 genotype array data (38 114 SNPs) of 831 Duroc pigs with four selection traits. In both heritability estimation and prediction accuracy, the accuracy depended on the number of animals with records. For heritability estimation, a large difference in the results among three relationship matrices was not shown, but the trend of the estimated heritabilities between GRMs, that is G_S  < G_H , was shown in this population. For the accuracy of prediction values in test animals, the accuracies of prediction values obtained by two GRMs were higher than that by the NRM in this population. The accuracies obtained by GRMs using animals with no records were lower than that by the NRM using animals with their performance records, but were close to that by the NRM using animals with full-sib testing records.

Systematic genotyping of groups of cows to improve genomic estimated breeding values of selection candidates

Background

Extending the reference set for genomic predictions in dairy cattle by adding large numbers of cows with genotypes and phenotypes has been proposed as a means to increase reliability of selection decisions for candidates.

Methods

In this study, we explored the potential of increasing the reliability of breeding values of young selection candidates by genotyping a fixed number of first-crop daughters of each sire from one or two generations in a balanced and regular system of genotyping. Using stochastic simulation, we developed a basic population scenario that mimics the situation in dual-purpose Fleckvieh cattle with respect to important key parameters. Starting with a reference set consisting of only genotyped bulls, we extended this reference set by including increasing numbers of daughter genotypes and phenotypes. We studied the effects on model-derived reliabilities, validation reliabilities and unbiasedness of predicted values for selection candidates. We also illustrate and discuss the effects of a selected sample and an unbalanced sampling of daughters. Furthermore, we quantified the role of selection with respect to the influence on validation reliabilities and contrasted these to model-derived reliabilities.

Results

In the most extended design, with 200 daughters per sire genotyped from two generations, single nucleotide polymorphism (SNP) effects were estimated from a reference set of 420,000 cows and 4200 bulls. For this design, the validation reliabilities for candidates reached 80 % or more, thereby exceeding the reliabilities that were achieved in traditional progeny-testing designs for a trait with moderate to high heritability. We demonstrate that even a moderate number of 25 genotyped daughters per sire will lead to considerable improvement in the reliability of predicted breeding values for selection candidates. Our results illustrate that the strategy applied to sample females for genotyping has a large impact on the benefits that can be achieved.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0250-9) contains supplementary material, which is available to authorized users.