Investigating the persistence of accuracy of genomic predictions over time in broilers

Accuracy of genomic predictions is an important component of the selection response. The objectives of this research were: 1) to investigate trends for prediction accuracies over time in a broiler population of accumulated phenotypes, genotypes, and pedigrees and 2) to test if data from distant generations are useful to maintain prediction accuracies in selection candidates. The data contained 820K phenotypes for a growth trait (GT), 200K for two feed efficiency traits (FE1 and FE2), and 42K for a carcass yield trait (CY). The pedigree included 1,252,619 birds hatched over 7 years, of which 154,318 from the last 4 years were genotyped. Training populations were constructed adding 1 year of data sequentially, persistency of accuracy over time was evaluated using predictions from birds hatched in the three generations following or in the years after the training populations. In the first generation, before genotypes became available for the training populations (first 3 years of data), accuracies remained almost stable with successive additions of phenotypes and pedigree to the accumulated dataset. The inclusion of 1 year of genotypes in addition to 4 years of phenotypes and pedigree in the training population led to increases in accuracy of 54% for GT, 76% for FE1, 110% for CY, and 38% for FE2; on average, 74% of the increase was due to genomics. Prediction accuracies declined faster without than with genomic information in the training populations. When genotypes were unavailable, the average decline in prediction accuracy across traits was 41% from the first to the second generation of validation, and 51% from the second to the third generation of validation. When genotypes were available, the average decline across traits was 14% from the first to the second generation of validation, and 3% from the second to the third generation of validation. Prediction accuracies in the last three generations were the same when the training population included 5 or 2 years of data, and a decrease of ~7% was observed when the training population included only 1 year of data. Training sets including genomic information provided an increase in accuracy and persistence of genomic predictions compared with training sets without genomic data. The two most recent years of pedigree, phenotypic, and genomic data were sufficient to maintain prediction accuracies in selection candidates. Similar conclusions were obtained using validation populations per year.     

Validation of single‐step GBLUP genomic predictions from threshold models using the linear regression method: An application in chicken mortality

The objective of this study was to determine whether the linear regression (LR) method could be used to validate genomic threshold models. Statistics for the LR method were computed from estimated breeding values (EBVs) using the whole and truncated data sets with variances from the reference and validation populations. The method was tested using simulated and real chicken data sets. The simulated data set included 10 generations of 4,500 birds each; genotypes were available for the last three generations. Each animal was assigned a continuous trait, which was converted to a binary score assuming an incidence of failure of 7%. The real data set included the survival status of 186,596 broilers (mortality rate equal to 7.2%) and genotypes of 18,047 birds. Both data sets were analysed using best linear unbiased predictor (BLUP) or single‐step GBLUP (ssGBLUP). The whole data set included all phenotypes available, whereas in the partial data set, phenotypes of the most recent generation were removed. In the simulated data set, the accuracies based on the LR formulas were 0.45 for BLUP and 0.76 for ssGBLUP, whereas the correlations between true breeding values and EBVs (i.e. true accuracies) were 0.37 and 0.65, respectively. The gain in accuracy by adding genomic information was overestimated by 0.09 when using the LR method compared to the true increase in accuracy. However, when the estimated ratio between the additive variance computed based on pedigree only and on pedigree and genomic information was considered, the difference between true and estimated gain was <0.02. Accuracies of BLUP and ssGBLUP with the real data set were 0.41 and 0.47, respectively. This small improvement in accuracy when using ssGBLUP with the real data set was due to population structure and lower heritability. The LR method is a useful tool for estimating improvements in accuracy of EBVs due to the inclusion of genomic information when traditional validation methods as k‐fold validation and predictive ability are not applicable.     

Selection of core animals in the Algorithm for Proven and Young using a simulation model

The Algorithm for Proven and Young (APY) enables the implementation of single‐step genomic BLUP (ssGBLUP) in large, genotyped populations by separating genotyped animals into core and non‐core subsets and creating a computationally efficient inverse for the genomic relationship matrix ( G ). As APY became the choice for large‐scale genomic evaluations in BLUP‐based methods, a common question is how to choose the animals in the core subset. We compared several core definitions to answer this question. Simulations comprised a moderately heritable trait for 95,010 animals and 50,000 genotypes for animals across five generations. Genotypes consisted of 25,500 SNP distributed across 15 chromosomes. Genotyping errors and missing pedigree were also mimicked. Core animals were defined based on individual generations, equal representation across generations, and at random. For a sufficiently large core size, core definitions had the same accuracies and biases, even if the core animals had imperfect genotypes. When genotyped animals had unknown parents, accuracy and bias were significantly better (p ≤ .05) for random and across generation core definitions.     

Genome-wide marker-assisted selection combining all pedigree phenotypic information with genotypic data in one step: An example using broiler chickens

Data of broiler chickens for 2 pure lines across 3 generations were used for genomic evaluation. A complete population (full data set; FDS) consisted of 183,784 and 164,246 broilers for the 2 lines. The genotyped subsets (SUB) consisted of 3,284 and 3,098 broilers with 57,636 SNP. Genotyped animals were preselected based on more than 20 traits with different index applied to each line. Three traits were analyzed: BW at 6 wk (BW6), ultrasound measurement of breast meat (BM), and leg score (LS) coded 1 = no and 2 = yes for leg defect. Some phenotypes were missing for BM. The training population consisted of the first 2 generations including all animals in FDS or only genotyped animals in SUB. The validation data set contained only genotyped animals in the third generation. Genetic evaluations were performed using 3 approaches: 1) phenotypic BLUP, 2) extending BLUP methodologies to utilize pedigree and genomic information in a single step (ssGBLUP), and 3) Bayes A. Whereas BLUP and ssGBLUP utilized all phenotypic data, Bayes A could use only those of the genotyped subset. Heritabilities were 0.17 to 0.20 for BW6, 0.30 to 0.35 for BM, and 0.09 to 0.11 for LS. The average accuracies of the validation population with BLUP for BW6, BM, and LS were 0.46, 0.30, and <0 with SUB and 0.51, 0.34, and 0.28 with FDS. With ssGBLUP, those accuracies were 0.60, 0.34, and 0.06 with SUB and 0.61, 0.40, and 0.37 with FDS, respectively. With Bayes A, the accuracies were 0.60, 0.36, and 0.09 with SUB. With SUB, Bayes A and ssGBLUP had similar accuracies. For traits of high heritability, the accuracy of Bayes A/SUB and ssGBLUP/FDS were similar, and up to 50% better than BLUP/FDS. However, with low heritability, ssGBLUP/FDS was 4 to 6 times more accurate than Bayes A/SUB and 50% better than BLUP/FDS. An optimal genomic evaluation would be multi-trait and involve all traits and records on which selection is based.     

Implications of SNP weighting on single‐step genomic predictions for different reference population sizes

We investigated the importance of SNP weighting in populations with 2,000 to 25,000 genotyped animals. Populations were simulated with two effective sizes (20 or 100) and three numbers of QTL (10, 50 or 500). Pedigree information was available for six generations; phenotypes were recorded for the four middle generations. Animals from the last three generations were genotyped for 45,000 SNP. Single‐step genomic BLUP (ssGBLUP) and weighted ssGBLUP (WssGBLUP) were used to estimate genomic EBV using a genomic relationship matrix ( G ). The WssGBLUP performed better in small genotyped populations; however, any advantage for WssGBLUP was reduced or eliminated when more animals were genotyped. WssGBLUP had greater resolution for genome‐wide association (GWA) as did increasing the number of genotyped animals. For few QTL, accuracy was greater for WssGBLUP than ssGBLUP; however, for many QTL, accuracy was the same for both methods. The largest genotyped set was used to assess the dimensionality of genomic information (number of effective SNP). The number of effective SNP was considerably less in weighted G than in unweighted G . Once the number of independent SNP is well represented in the genotyped population, the impact of SNP weighting becomes less important.     

Parameters affecting genome simulation for evaluating genomic selection method

The present study investigated the parameter settings for obtaining a simulated genome at steady state of allele frequency (mutation–drift equilibrium) and linkage disequilibrium (LD), and evaluated the impact of whether or not the simulated genome reached steady state of allele frequency and LD on the accuracy of genomic estimated breeding values (GEBVs). After 500 to 50 000 historical generations, the base population and subsequent seven generations were generated as recent populations. The allele frequency distribution of the last generations of the historical population and LD in the base population were calculated when varying the values of five parameters: initial minor allele frequency, mutation rate, effective population size, number of markers and chromosome length. The accuracies of GEBVs in the last generation of the recent population were calculated by genomic best linear unbiased prediction. The number of historical generations required to reach mutation–drift equilibrium depended on the initial allele frequency and mutation rate. Regardless of the parameters, LD reached a steady state before allele frequency distribution reached mutation–drift equilibrium. The accuracies of GEBVs largely reflect the extent of linkage disequilibrium with the exception of varying chromosome length, although there were no associations between the accuracies of GEBVs and allele frequency distribution.     

Impact of fitting dominance and additive effects on accuracy of genomic prediction of breeding values in layers

Most genomic prediction studies fit only additive effects in models to estimate genomic breeding values (GEBV). However, if dominance genetic effects are an important source of variation for complex traits, accounting for them may improve the accuracy of GEBV. We investigated the effect of fitting dominance and additive effects on the accuracy of GEBV for eight egg production and quality traits in a purebred line of brown layers using pedigree or genomic information (42K single‐nucleotide polymorphism (SNP) panel). Phenotypes were corrected for the effect of hatch date. Additive and dominance genetic variances were estimated using genomic‐based [genomic best linear unbiased prediction (GBLUP)‐REML and BayesC] and pedigree‐based (PBLUP‐REML) methods. Breeding values were predicted using a model that included both additive and dominance effects and a model that included only additive effects. The reference population consisted of approximately 1800 animals hatched between 2004 and 2009, while approximately 300 young animals hatched in 2010 were used for validation. Accuracy of prediction was computed as the correlation between phenotypes and estimated breeding values of the validation animals divided by the square root of the estimate of heritability in the whole population. The proportion of dominance variance to total phenotypic variance ranged from 0.03 to 0.22 with PBLUP‐REML across traits, from 0 to 0.03 with GBLUP‐REML and from 0.01 to 0.05 with BayesC. Accuracies of GEBV ranged from 0.28 to 0.60 across traits. Inclusion of dominance effects did not improve the accuracy of GEBV, and differences in their accuracies between genomic‐based methods were small (0.01–0.05), with GBLUP‐REML yielding higher prediction accuracies than BayesC for egg production, egg colour and yolk weight, while BayesC yielded higher accuracies than GBLUP‐REML for the other traits. In conclusion, fitting dominance effects did not impact accuracy of genomic prediction of breeding values in this population.     

Partial-genome evaluation of postweaning feed intake and efficiency of crossbred beef cattle

The effects of individual SNP and the variation explained by sets of SNP associated with DMI, metabolic midtest BW, BW gain, and feed efficiency, expressed as phenotypic and genetic residual feed intake, were estimated from BW and the individual feed intake of 1,159 steers on dry lot offered a 3.0 Mcal/kg ration for at least 119 d before slaughter. Parents of these F(1) × F(1) (F(1)(2)) steers were AI-sired F(1) progeny of Angus, Charolais, Gelbvieh, Hereford, Limousin, Red Angus, and Simmental bulls mated to US Meat Animal Research Center Angus, Hereford, and MARC III composite females. Steers were genotyped with the BovineSNP50 BeadChip assay (Illumina Inc., San Diego, CA). Effects of 44,163 SNP having minor allele frequencies >0.05 in the F(1)(2) generation were estimated with a mixed model that included genotype, breed composition, heterosis, age of dam, and slaughter date contemporary groups as fixed effects, and a random additive genetic effect with recorded pedigree relationships among animals. Variance in this population attributable to sets of SNP was estimated with models that partitioned the additive genetic effect into a polygenic component attributable to pedigree relationships and a genotypic component attributable to genotypic relationships. The sets of SNP evaluated were the full set of 44,163 SNP and subsets containing 6 to 40,000 SNP selected according to association with phenotype. Ninety SNP were strongly associated (P < 0.0001) with at least 1 efficiency or component trait; these 90 accounted for 28 to 46% of the total additive genetic variance of each trait. Trait-specific sets containing 96 SNP having the strongest associations with each trait explained 50 to 87% of additive variance for that trait. Expected accuracy of steer breeding values predicted with pedigree and genotypic relationships exceeded the accuracy of their sires predicted without genotypic information, although gains in accuracy were not sufficient to encourage that performance testing be replaced by genotyping and genomic evaluations.     

一种适用于单胎动物多品种选择的基因组关系矩阵

旨在提出一种新型基因组关系矩阵并验证其在多品种联合群体中的模拟应用效果。本研究利用QMsim软件模拟牛的表型数据和基因型数据;利用Gmatrix软件构建常规G阵;利用R语言构建新型G阵,新型G阵在常规G阵的基础上,将多品种联合群体的非哈代-温伯格平衡位点考虑在内;利用DMU软件使用“一步”法模型计算基因组估计育种值（estimated genomic breeding value,GEBV）;比较不同情况下使用两种G阵的GEBV预测准确性。结果表明,在不同遗传力及QTL数下,不对新型G阵使用A₂₂阵加权就能达到常规G阵使用A₂₂阵加权时的GEBV预测准确性。在系谱部分缺失时,新型G阵不加权较常规G阵加权时GEBV预测准确性高。证明,在系谱有部分缺失时,新型G阵对多品种GEBV的预测有一定优势。     

Effect of genotyped bulls with different numbers of phenotyped progenies on quantitative trait loci detection and genomic evaluation in a simulated cattle population

The objective of this study was to assess the effect of genotyped bulls with different numbers of phenotyped progenies on quantitative trait loci (QTL) detection and genomic evaluation using a simulated cattle population. Twelve generations (G1–G12) were simulated from the base generation (G0). The recent population had different effective population sizes, heritability, and number of QTL. G0–G4 were used for pedigree information. A total of 300 genotyped bulls from G5–G10 were randomly selected. Their progenies were generated in G6–G11 with different numbers of progeny per bull. Scenarios were considered according to the number of progenies and whether the genotypes were possessed by the bulls or the progenies. A genome‐wide association study and genomic evaluation were performed with a single‐step genomic best linear unbiased prediction method to calculate the power of QTL detection and the genomic estimated breeding value (GEBV). We found that genotyped bulls could be available for QTL detection depending on conditions. Additionally, using a reference population, including genotyped bulls, which had more progeny phenotypes, enabled a more accurate prediction of GEBV. However, it is desirable to have more than 4,500 individuals consisting of both genotypes and phenotypes for practical genomic evaluation.     

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.     

Persistency of accuracy of genomic breeding values for different simulated pig breeding programs in developing countries

Genetic improvement of pigs in tropical developing countries has focused on imported exotic populations which have been subjected to intensive selection with attendant high population‐wide linkage disequilibrium (LD). Presently, indigenous pig population with limited selection and low LD are being considered for improvement. Given that the infrastructure for genetic improvement using the conventional BLUP selection methods are lacking, a genome‐wide selection (GS) program was proposed for developing countries. A simulation study was conducted to evaluate the option of using 60 K SNP panel and observed amount of LD in the exotic and indigenous pig populations. Several scenarios were evaluated including different size and structure of training and validation populations, different selection methods and long‐term accuracy of GS in different population/breeding structures and traits. The training set included previously selected exotic population, unselected indigenous population and their crossbreds. Traits studied included number born alive (NBA), average daily gain (ADG) and back fat thickness (BFT). The ridge regression method was used to train the prediction model. The results showed that accuracies of genomic breeding values (GBVs) in the range of 0.30 (NBA) to 0.86 (BFT) in the validation population are expected if high density marker panels are utilized. The GS method improved accuracy of breeding values better than pedigree‐based approach for traits with low heritability and in young animals with no performance data. Crossbred training population performed better than purebreds when validation was in populations with similar or a different structure as in the training set. Genome‐wide selection holds promise for genetic improvement of pigs in the tropics.     

Genomic predictions across Nordic Holstein and Nordic Red using the genomic best linear unbiased prediction model with different genomic relationship matrices

This study investigated genomic predictions across Nordic Holstein and Nordic Red using various genomic relationship matrices. Different sources of information, such as consistencies of linkage disequilibrium (LD) phase and marker effects, were used to construct the genomic relationship matrices (G‐matrices) across these two breeds. Single‐trait genomic best linear unbiased prediction (GBLUP) model and two‐trait GBLUP model were used for single‐breed and two‐breed genomic predictions. The data included 5215 Nordic Holstein bulls and 4361 Nordic Red bulls, which was composed of three populations: Danish Red, Swedish Red and Finnish Ayrshire. The bulls were genotyped with 50 000 SNP chip. Using the two‐breed predictions with a joint Nordic Holstein and Nordic Red reference population, accuracies increased slightly for all traits in Nordic Red, but only for some traits in Nordic Holstein. Among the three subpopulations of Nordic Red, accuracies increased more for Danish Red than for Swedish Red and Finnish Ayrshire. This is because closer genetic relationships exist between Danish Red and Nordic Holstein. Among Danish Red, individuals with higher genomic relationship coefficients with Nordic Holstein showed more increased accuracies in the two‐breed predictions. Weighting the two‐breed G‐matrices by LD phase consistencies, marker effects or both did not further improve accuracies of the two‐breed predictions.     

结合GWAS先验标记信息的肉鸡RFI性状全基因组选择研究

旨在比较结合全基因组关联分析（genome-wide association study,GWAS）先验标记信息的基因组育种值（genomic estimated breeding value,GEBV）估计与基因组最佳线性无偏预测（genomic best linear unbiased prediction,GBLUP）方法对鸡剩余采食量性状育种值估计的准确性,为提高基因组选择准确性提供理论与技术支持。本研究选用广西金陵花鸡3个世代共2 510个个体作为素材,其中公鸡1 648只,母鸡862只,以42～56日龄期间的剩余采食量（residual feed intake,RFI）为目标性状,将试验群体随机分为两组,其中一组作为先验标记信息发现群体,用于GWAS分析并筛选最显著的top5%、top10%、top15%和top20%的位点作为先验标记信息;另外一组分别结合不同的先验标记信息进行遗传参数估计并比较基因组育种值的预测准确性,使用重复10次的五倍交叉验证法获取准确性,随后两组群体再进行交叉验证。研究结果表明,GBLUP计算RFI的遗传力为0.153,预测准确性为0.387～0.429,结合GWAS先验标记信息的基因组选择方法计算RFI的遗传力为0.139～0.157,预测准确性为0.401～0.448。将GWAS结果中P值最显著的top10%～top15%的SNPs作为先验信息整合至基因组选择模型中可以将RFI的预测准确性提升2.10%～5.17%。     

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.     

Performance of genomic selection under a single-step approach in autochthonous Spanish beef cattle populations

This study evaluated different strategies for implementing a single-step genomic selection programme in two autochthonous Spanish beef cattle populations (Pirenaica—Pi and Rubia Gallega—RG). The strategies were compared in terms of accuracy attained under different scenarios by simulating genomic data over the known genealogy. Several genotyping approaches were tested, as well as, other factors like marker density, effective population size, mutation rate and heritability of the trait. The results obtained showed gains in accuracy with respect to pedigree BLUP evaluation in all cases. The greatest benefit was obtained when the candidates to selection had their genotypes included in the evaluation. Moreover, genotyping the individuals with the most accurate predictions maximized the gains but other suboptimal strategies also yielded satisfactory results. Furthermore, the gains in accuracy increased with the marker density reaching a plateau at around 50,000 markers. Likewise, the effective population size and the mutation rate have also shown an effect, both increasing the accuracy with decreasing values of these population parameters. Finally, the results obtained for the RG population showed greater gains compared to the Pi population, probably attributed to the wider implantation of artificial insemination.     

Estimates of genetic variance and variance of predicted genetic merits using pedigree or genomic relationship matrices in six Brown Swiss cattle populations for different traits

The amount of variance captured in genetic estimations may depend on whether a pedigree‐based or genomic relationship matrix is used. The purpose of this study was to investigate the genetic variance as well as the variance of predicted genetic merits (PGM) using pedigree‐based or genomic relationship matrices in Brown Swiss cattle. We examined a range of traits in six populations amounting to 173 population‐trait combinations. A main aim was to determine how using different relationship matrices affect variance estimation. We calculated ratios between different types of estimates and analysed the impact of trait heritability and population size. The genetic variances estimated by REML using a genomic relationship matrix were always smaller than the variances that were similarly estimated using a pedigree‐based relationship matrix. The variances from the genomic relationship matrix became closer to estimates from a pedigree relationship matrix as heritability and population size increased. In contrast, variances of predicted genetic merits obtained using a genomic relationship matrix were mostly larger than variances of genetic merit predicted using pedigree‐based relationship matrix. The ratio of the genomic to pedigree‐based PGM variances decreased as heritability and population size rose. The increased variance among predicted genetic merits is important for animal breeding because this is one of the factors influencing genetic progress.     

Genomic evaluation of regional dairy cattle breeds in single‐breed and multibreed contexts

An important prerequisite for high prediction accuracy in genomic prediction is the availability of a large training population, which allows accurate marker effect estimation. This requirement is not fulfilled in case of regional breeds with a limited number of breeding animals. We assessed the efficiency of the current French routine genomic evaluation procedure in four regional breeds (Abondance, Tarentaise, French Simmental and Vosgienne) as well as the potential benefits when the training populations consisting of males and females of these breeds are merged to form a multibreed training population. Genomic evaluation was 5–11% more accurate than a pedigree‐based BLUP in three of the four breeds, while the numerically smallest breed showed a < 1% increase in accuracy. Multibreed genomic evaluation was beneficial for two breeds (Abondance and French Simmental) with maximum gains of 5 and 8% in correlation coefficients between yield deviations and genomic estimated breeding values, when compared to the single‐breed genomic evaluation results. Inflation of genomic evaluation of young candidates was also reduced. Our results indicate that genomic selection can be effective in regional breeds as well. Here, we provide empirical evidence proving that genetic distance between breeds is only one of the factors affecting the efficiency of multibreed genomic evaluation.     

Is the use of formulae a reliable way to predict the accuracy of genomic selection?

We studied four formulae used to predict the accuracy of genomic selection prior to genotyping. The objectives of our study were to investigate the impact of the parameters of each formula on the values of accuracy calculated using these formulae, and to check whether the accuracies reported in the literature are in agreement with the formulae. First, we computed the marginal distribution of accuracy (by integration) for each parameter of all four formulae: heritability h², reference population size T, number of markers M and number of effective segments in the genome M_e. Then, we collected 145 accuracies and corresponding parameters reported in 13 publications on genomic selection (mainly in dairy cattle), and performed analysis of variance to test the differences between observed and predicted accuracy with effects of formulae and parameters. The variation of accuracy for different values of each parameter indicated that two parameters, T and M_e, had a significant impact and that considerable differences existed between the formulae (mean accuracies differed by up to 0.20 point). The results of our meta‐analysis showed a big formula effect on the accuracies predicted using each formula, and also a significant effect of the value obtained for M_e calculated from N_e (effective population size). Each formula can therefore be demonstrated to be optimal depending on the assumption used for M_e. In conclusion, no rules can be applied to predict the reliability of genomic selection using these formulae.     

Differences between genomic‐based and pedigree‐based relationships in a chicken population,as a function of quality control and pedigree links among individuals

This work studied differences between expected (calculated from pedigree) and realized (genomic, from markers) relationships in a real population, the influence of quality control on these differences, and their fit to current theory. Data included 4940 pure line chickens across five generations genotyped for 57 636 SNP. Pedigrees (5762 animals) were available for the five generations, pedigree starting on the first one. Three levels of quality control were used. With no quality control, mean difference between realized and expected relationships for different type of relationships was ≤ 0.04 with standard deviation ≤ 0.10. With strong quality control (call rate ≥ 0.9, parent‐progeny conflicts, minor allele frequency and use of only autosomal chromosomes), these numbers reduced to ≤ 0.02 and ≤ 0.04, respectively. While the maximum difference was 1.02 with the complete data, it was only 0.18 with the latest three generations of genotypes (but including all pedigrees). Variation of expected minus realized relationships agreed with theoretical developments and suggests an effective number of loci of 70 for this population. When the pedigree is complete and as deep as the genotypes, the standard deviation of difference between the expected and realized relationships is around 0.04, all categories confounded. Standard deviation of differences larger than 0.10 suggests bad quality control, mistakes in pedigree recording or genotype labelling, or insufficient depth of pedigree.