Genetic parameters for milk fatty acids,milk yield and quality traits of a Holstein cattle population reared under tropical conditions

Information about genetic parameters is essential for selection decisions and genetic evaluation. These estimates are population specific; however, there are few studies with dairy cattle populations reared under tropical and sub‐tropical conditions. Thus, the aim was to obtain estimates of heritability and genetic correlations for milk yield and quality traits using pedigree and genomic information from a Holstein population maintained in a tropical environment. Phenotypic records (n = 36 457) of 4203 cows as well as the genotypes for 57 368 single nucleotide polymorphisms from 755 of these cows were used. Covariance components were estimated using the restricted maximum likelihood method under a mixed animal model, considering a pedigree‐based relationship matrix or a combined pedigree‐genomic matrix. High heritabilities (around 0.30) were estimated for lactose and protein content in milk whereas moderate values (between 0.19 and 0.26) were obtained for percentages of fat, saturated fatty acids and palmitic acid in milk. Genetic correlations ranging from −0.38 to −0.13 were determined between milk yield and composition traits. The smaller estimates compared to other similar studies can be due to poor environmental conditions, which may reduce genetic variability. These results highlight the importance in using genetic parameters estimated in the population under evaluation for selection decisions.     

闭锁群体BLUP选择下遗传方差和近交系数的变化

采用计算机随机模拟方法模拟了在一个闭锁群体内连续对单个性状进行 1 5个世代选择的情况。选择过程中世代不重叠 ,每个世代的种畜根据动物模型最佳线性无偏预测 (BLUP)法估计的育种值进行选留 ,并在此基础上系统地比较了不同群体规模、公母比例和性状遗传力对群体遗传方差和近交系数变化的影响。结果表明 ,扩大育种群规模、增加公畜比例以及对低遗传力性状进行选择时 ,群体遗传方差降低的速度和近交系数上升的速度会更慢 ,在长期选择时可望获得更大的持续进展和适宜的近交增量     

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.     

Using the genomic relationship matrix to predict the accuracy of genomic selection

Estimated breeding values (EBVs) using data from genetic markers can be predicted using a genomic relationship matrix, derived from animal's genotypes, and best linear unbiased prediction. However, if the accuracy of the EBVs is calculated in the usual manner (from the inverse element of the coefficient matrix), it is likely to be overestimated owing to sampling errors in elements of the genomic relationship matrix. We show here that the correct accuracy can be obtained by regressing the relationship matrix towards the pedigree relationship matrix so that it is an unbiased estimate of the relationships at the QTL controlling the trait. This method shows how the accuracy increases as the number of markers used increases because the regression coefficient (of genomic relationship towards pedigree relationship) increases. We also present a deterministic method for predicting the accuracy of such genomic EBVs before data on individual animals are collected. This method estimates the proportion of genetic variance explained by the markers, which is equal to the regression coefficient described above, and the accuracy with which marker effects are estimated. The latter depends on the variance in relationship between pairs of animals, which equals the mean linkage disequilibrium over all pairs of loci. The theory was validated using simulated data and data on fat concentration in the milk of Holstein cattle.     

Adjusting for heterogeneity of experimental data in genetic evaluation of dry matter intake in dairy cattle

The objectives of the present study were (i) to find the best fitted model for repeatedly measured daily dry matter intake (DMI) data obtained from different herds and experiments across lactations and (ii) to get better estimates of the genetic parameters and better genetic evaluations. After editing, there were 572,512 daily DMI records of 3,495 animals (Holstein cows) from 11 different herds across 13 lactations and the animals were under 110 different nutritional experiments. The fitted model for this data set was a univariate repeated‐measure animal model (called model 1) in which additive genetic and permanent environmental (within and across lactations) effects were fitted as random. Model 1 was fitted as two distinct models (called models 2 and 3) based on alternative fixed effect corrections. For unscaled data, each model (models 2 and 3) was fitted as a homoscedastic (HOM) model first and then as a heteroscedastic (HET) model. Then, data were scaled by multiplying with particular herd‐scaling factors, which were calculated by accounting for heterogeneity of phenotypic within‐herd variances. Models were selected based on cross‐validation and prediction accuracy results. Scaling factors were re‐estimated to determine the effectiveness of accounting for herd heterogeneity. Variance components and respective heritability and repeatability were estimated based on a pedigree‐based relationship matrix. Results indicated that the model fitted for scaled data showed better fit than the models (HOM or HET) fitted for unscaled data. The heritability estimates of the models 2 and 3 fitted for scaled data were 0.30 and 0.08, respectively. The repeatability estimates of the model fitted for scaled data ranged from 0.51 to 0.63. The re‐estimated scaling factor after accounting for heterogeneity of residual variances was close to 1.0, indicating the stabilization of residual variances and herd accounted for most of the heterogeneity. The rank correlation of EBVs between scaled and unscaled data ranged from 0.96 to 0.97.     

Bayesian inference strategies for the prediction of genetic merit using threshold models with an application to calving ease scores in Italian Piemontese cattle

First parity calving difficulty scores from Italian Piemontese cattle were analysed using a threshold mixed effects model. The model included the fixed effects of age of dam and sex of calf and their interaction and the random effects of sire, maternal grandsire, and herd‐year‐season. Covariances between sire and maternal grandsire effects were modelled using a numerator relationship matrix based on male ancestors. Field data consisted of 23 953 records collected between 1989 and 1998 from 4741 herd‐year‐seasons. Variance and covariance components were estimated using two alternative approximate marginal maximum likelihood (MML) methods, one based on expectation‐maximization (EM) and the other based on Laplacian integration. Inferences were compared to those based on three separate runs or sequences of Markov Chain Monte Carlo (MCMC) sampling in order to assess the validity of approximate MML estimates derived from data with similar size and design structure. Point estimates of direct heritability were 0.24, 0.25 and 0.26 for EM, Laplacian and MCMC (posterior mean), respectively, whereas corresponding maternal heritability estimates were 0.10, 0.11 and 0.12, respectively. The covariance between additive direct and maternal effects was found to be not different from zero based on MCMC‐derived confidence sets. The conventional joint modal estimates of sire effects and associated standard errors based on MML estimates of variance and covariance components differed little from the respective posterior means and standard deviations derived from MCMC. Therefore, there may be little need to pursue computation‐intensive MCMC methods for inference on genetic parameters and genetic merits using conventional threshold sire and maternal grandsire models for large datasets on calving ease.     

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%.     

Variational bayesian method of estimating variance components

We developed a Bayesian analysis approach by using a variational inference method, a so‐called variational Bayesian method, to determine the posterior distributions of variance components. This variational Bayesian method and an alternative Bayesian method using Gibbs sampling were compared in estimating genetic and residual variance components from both simulated data and publically available real pig data. In the simulated data set, we observed strong bias toward overestimation of genetic variance for the variational Bayesian method in the case of low heritability and low population size, and less bias was detected with larger population sizes in both methods examined. The differences in the estimates of variance components between the variational Bayesian and the Gibbs sampling were not found in the real pig data. However, the posterior distributions of the variance components obtained with the variational Bayesian method had shorter tails than those obtained with the Gibbs sampling. Consequently, the posterior standard deviations of the genetic and residual variances of the variational Bayesian method were lower than those of the method using Gibbs sampling. The computing time required was much shorter with the variational Bayesian method than with the method using Gibbs sampling.     

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.     

Cluster and meta‐analyses of genetic parameters for feed intake traits in growing beef cattle

A data set based on 50 studies including feed intake and utilization traits was used to perform a meta‐analysis to obtain pooled estimates using the variance between studies of genetic parameters for average daily gain (ADG); residual feed intake (RFI); metabolic body weight (MBW); feed conversion ratio (FCR); and daily dry matter intake (DMI) in beef cattle. The total data set included 128 heritability and 122 genetic correlation estimates published in the literature from 1961 to 2012. The meta‐analysis was performed using a random effects model where the restricted maximum likelihood estimator was used to evaluate variances among clusters. Also, a meta‐analysis using the method of cluster analysis was used to group the heritability estimates. Two clusters were obtained for each trait by different variables. It was observed, for all traits, that the heterogeneity of variance was significant between clusters and studies for genetic correlation estimates. The pooled estimates, adding the variance between clusters, for direct heritability estimates for ADG, DMI, RFI, MBW and FCR were 0.32 ± 0.04, 0.39 ± 0.03, 0.31 ± 0.02, 0.31 ± 0.03 and 0.26 ± 0.03, respectively. Pooled genetic correlation estimates ranged from ?0.15 to 0.67 among ADG, DMI, RFI, MBW and FCR. These pooled estimates of genetic parameters could be used to solve genetic prediction equations in populations where data is insufficient for variance component estimation. Cluster analysis is recommended as a statistical procedure to combine results from different studies to account for heterogeneity.     

Use of molecular markers to improve relationship information in the genetic evaluation of beef cattle tick resistance under pedigree‐based models

The selection of genetically superior individuals is conditional upon accurate breeding value predictions which, in turn, are highly depend on how precisely relationship is represented by pedigree. For that purpose, the numerator relationship matrix is essential as a priori information in mixed model equations. The presence of pedigree errors and/or the lack of relationship information affect the genetic gain because it reduces the correlation between the true and estimated breeding values. Thus, this study aimed to evaluate the effects of correcting the pedigree relationships using single‐nucleotide polymorphism (SNP) markers on genetic evaluation accuracies for resistance of beef cattle to ticks. Tick count data from Hereford and Braford cattle breeds were used as phenotype. Genotyping was carried out using a high‐density panel (BovineHD ‐ Illumina^® bead chip with 777 962 SNPs) for sires and the Illumina BovineSNP50 panel (54 609 SNPs) for their progenies. The relationship between the parents and progenies of genotyped animals was evaluated, and mismatches were based on the Mendelian conflicts counts. Variance components and genetic parameters estimates were obtained using a Bayesian approach via Gibbs sampling, and the breeding values were predicted assuming a repeatability model. A total of 460 corrections in relationship definitions were made (Table 1) corresponding to 1018 (9.5%) tick count records. Among these changes, 97.17% (447) were related to the sire's information, and 2.8% (13) were related to the dam's information. We observed 27.2% (236/868) of Mendelian conflicts for sire–progeny genotyped pairs and 14.3% (13/91) for dam–progeny genotyped pairs. We performed 2174 new definitions of half‐siblings according to the correlation coefficient between the coancestry and molecular coancestry matrices. It was observed that higher‐quality genetic relationships did not result in significant differences of variance components estimates; however, they resulted in more accurate breeding values predictions. Using SNPs to assess conflicts between parents and progenies increases certainty in relationships and consequently the accuracy of breeding value predictions of candidate animals for selection. Thus, higher genetic gains are expected when compared to the traditional non‐corrected relationship matrix.     

On estimating components of genetic variance in diallel matings 1

1. The relative importance of additive and non‐additive genetic effects on body weight, egg weight, maturity and rate of egg production were studied from diallel matings in a Leghorn population.

2. From analyses of variance, heritability estimates of the additive fractions, based on half‐sib variances, and the non‐additive or dominance fractions, based on the sire x dam interaction component were obtained.

3. Non‐additive genetic effects were not statistically significant for any of the traits, though for rate of egg production at 32 and 62 weeks, the non‐additive effects as proportion of total variances were 0.29 (P<0.10) and 0.20 (P<0.16), respectively, compared additive effects of 0.08 (NS) and 0.11 (P<0.05).

4. The ratios of non‐additive to additive variances, 1.89 and 3.62 respectively, give support to inbreeding and hybridisation or reciprocal recurrent selection as methods of genetic improvement of egg production.     

Estimation of heritability for Tying-up syndrome in the Thoroughbred racehorse by Gibbs sampling

Tying‐up is a condition that primarily affects the muscles of horses. In this study, the heritability of the Tying‐up syndrome in the Thoroughbred racehorse was estimated by Bayesian analysis with Gibbs sampling based on the threshold model for binary traits. The data used were the clinical data in racehorses diagnosed by veterinarians of the Racehorse Clinics of Japan Racing Association from 2000 to 2003. The health status of the Tying‐up was treated as a binary trait. In the genetic analysis, the effect of changing the amount of the pedigree or inbreeding information on the estimation of heritability was investigated, too. The heritability estimates with non‐zero probability in the posterior densities were approximately 0.16–0.18 in minimum, suggesting that the heritability of the Tying‐up is not zero at least. The posterior density distributions of the heritability estimates were generally more pointed and sharp with using inbreeding coefficients than without using it, suggesting that more stable estimations were obtained when inbreeding coefficients were used. Among the different amounts of pedigree and inbreeding information, the heritabilities obtained with three or four generations of pedigree using inbreeding coefficients seems to be preferable, i.e. heritability of 0.42 or 0.43 for Tying‐up.     

Estimation of additive genetic variance when base populations are selected

A population of size 40 was simulated 1,000 times for 10 generations. Five out of twenty males were selected each generation, and each male was mated to four females to have two progeny. The additive genetic variance (sigma 2a) before selection was 10, and the initial heritability was .5. Due to covariances among animals, inbreeding and gametic disequilibrium, the genetic variance was reduced to 6.72 after 10 generations of selection. Reduction of variance was lower in another population simulated with size 400 and 10% of the males selected. Restricted Maximum Likelihood was used to estimate sigma 2a using an animal model. The estimate of sigma 2a was empirically unbiased when all data and all relationships were used. Omitting data from selected ancestors caused biased estimates of sigma 2a due to not accounting for all gametic disequilibrium. Including additional relationships between assumed base animals adjusted for inbreeding and for covariances. Bias from gametic disequilibrium decreased slightly with the use of more relationship information, and it was smaller in the small population and(or) when selection had been practiced for just a few generations.     

Estimation of genetic parameters and response to selection for a continuous trait subject to culling before testing

The consequences of assuming a zero environmental covariance between a binary trait 'test-status' and a continuous trait on the estimates of genetic parameters by restricted maximum likelihood and Gibbs sampling and on response from genetic selection when the true environmental covariance deviates from zero were studied. Data were simulated for two traits (one that culling was based on and a continuous trait) using the following true parameters, on the underlying scale: h2 = 0.4; r(A) = 0.5; r(E) = 0.5, 0.0 or -0.5. The selection on the continuous trait was applied to five subsequent generations where 25 sires and 500 dams produced 1500 offspring per generation. Mass selection was applied in the analysis of the effect on estimation of genetic parameters. Estimated breeding values were used in the study of the effect of genetic selection on response and accuracy. The culling frequency was either 0.5 or 0.8 within each generation. Each of 10 replicates included 7500 records on 'test-status' and 9600 animals in the pedigree file. Results from bivariate analysis showed unbiased estimates of variance components and genetic parameters when true r(E) = 0.0. For r(E) = 0.5, variance components (13-19% bias) and especially (50-80%) were underestimated for the continuous trait, while heritability estimates were unbiased. For r(E) = -0.5, heritability estimates of test-status were unbiased, while genetic variance and heritability of the continuous trait together with were overestimated (25-50%). The bias was larger for the higher culling frequency. Culling always reduced genetic progress from selection, but the genetic progress was found to be robust to the use of wrong parameter values of the true environmental correlation between test-status and the continuous trait. Use of a bivariate linear-linear model reduced bias in genetic evaluations, when data were subject to culling.     

The efficiency of genome-wide selection for genetic improvement of net merit

Four methods of selection for net merit comprising 2 correlated traits were compared in this study: 1) EBV-only index (I?), which consists of the EBV of both traits (i.e., traditional 2-trait BLUP selection); 2) GEBV-only index (I?), which comprises the genomic EBV (GEBV) of both traits; 3) GEBV-assisted index (I?), which combines both the EBV and the GEBV of both traits; and 4) GBV-assisted index (I?), which combines both the EBV and the true genomic breeding value (GBV) of both traits. Comparisons of these indices were based on 3 evaluation criteria [selection accuracy, genetic response (ΔH), and relative efficiency] under 64 scenarios that arise from combining 2 levels of genetic correlation (r(G)), 2 ratios of genetic variances between traits, 2 ratios of the genomic variance to total genetic variances for trait 1, 4 accuracies of EBV, and 2 proportions of r(G) explained by the GBV. Both selection accuracy and genetic responses of the indices I?, I?, and I? increased as the accuracy of EBV increased, but the efficiency of the indices I? and I? relative to I? decreased as the accuracy of EBV increased. The relative efficiency of both I? and I? was generally greater when the accuracy of EBV was 0.6 than when it was 0.9, suggesting that the genomic markers are most useful to assist selection when the accuracy of EBV is low. The GBV-assisted index I? was superior to the GEBV-assisted I? in all 64 cases examined, indicating the importance of improving the accuracy of prediction of genomic breeding values. Other parameters being identical, increasing the genetic variance of a high heritability trait would increase the genetic response of the genomic indices (I?, I?, and I?). The genetic responses to I?, I?, and I(4) was greater when the genetic correlation between traits was positive (r(G) = 0.5) than when it was negative (r(G) = -0.5). The results of this study indicate that the effectiveness of the GEBV-assisted index I? is affected by heritability of and genetic correlation between traits, the ratio of genetic variances between traits, the genomic-genetic variance ratio of each index trait, the proportion of genetic correlation accounted for by the genomic markers, and the accuracy of predictions of both EBV and GBV. However, most of these affecting factors are genetic characteristics of a population that is beyond the control of the breeders. The key factor subject to manipulation is to maximize both the proportion of the genetic variance explained by GEBV and the accuracy of both GEBV and EBV. The developed procedures provide means to investigate the efficiency of various genomic indices for any given combination of the genetic factors studied.     

The use of coancestry based on shared segments for maintaining genetic diversity

We have evaluated the use of genomic coancestry coefficients based on shared segments for the maintenance of genetic diversity through optimal contributions methodology for populations of three different Austrian cattle breeds. This coancestry measure has been compared with the genomic coancestry coefficient calculated on a SNP‐by‐SNP basis and with pedigree‐based coancestry. The regressions of the shared segments coancestry on the other two coefficients suggest that the former mainly reflect Identity By Descent but with the advantage over pedigree‐based coancestry of providing the realized Identity By Descent rather than an expectation. The effective population size estimated from the rate of coancestry based on shared segments was very similar to those obtained with the other coefficients and of small magnitude (from 26.24 to 111.90). This result highlights the importance of implementing active management strategies to control the increase of inbreeding and the loss of genetic diversity in livestock breeds, even when the population size is reasonably large. One problem for the implementation of coancestry based on shared segments is the need of estimating the gametic phases of the SNPs which, given the techniques used to obtain the genotypes, are a priori unknown. This study shows, through computer simulations, that using estimates of gametic phases for computing coancestry based on shared segments does not lead to a significant loss in the diversity maintained. This has been shown to be true even when the size of the population is very small as it is usually the case in populations subjected to conservation programmes.     

Impact of variation in length of individual testing periods on estimation of (co)variance components of a random regression model for feed intake of growing pigs

A simulation study was conducted to assess the influence of differences in the length of individual testing periods on estimates of (co)variance components of a random regression model for daily feed intake of growing pigs performance tested between 30 and 100 kg live weight. A quadratic polynomial in days on test with fixed regressions for sex, random regressions for additive genetic and permanent environmental effects and a constant residual variance was used for a bivariate simulation of feed intake and daily gain. (Co)variance components were estimated for feed intake only by means of a Bayesian analysis using Gibbs sampling and restricted maximum likelihood (REML). A single trait random regression model analogous to the one used for data simulation was used to analyse two versions of the data: full data sets with 18 weekly means of feed intake per animal and reduced data sets with the individual length of testing periods determined when tested animals reached 100 kg live weight. Only one significant difference between estimates from full and reduced data (REML estimate of genetic covariance between linear and quadratic regression parameters) and two significant differences from expected values (Gibbs estimates of permanent environmental variance of quadratic regression parameters) occurred. These differences are believed to be negligible, as the number lies within the expected range of type I error when testing at the 5% level. The course of test day variances calculated from estimates of additive genetic and permanent environmental covariance matrices also supports the conclusion that no bias in estimates of (co)variance components occurs due to the individual length of testing periods of performance‐tested growing pigs. A lower number of records per tested animal only results in more variation among estimates of (co)variance components from reduced compared with full data sets. Compared with the full data, the effective sample size of Gibbs samples from the reduced data decreased to 18% for residual variance and increased up to five times for other (co)variances. The data structure seems to influence the mixing of Gibbs chains.     

The effect of including genomic relationships in the estimation of genetic parameters of functional traits in pigs

The term functionality in animal breeding is used for traits that increase the efficiency of production by lowering the input cost, such as animal health and leg weakness related to longevity. The main objective of the study was to investigate the impact of genomic information, in a multivariate variance component analysis, on some of these traits. In addition, the effect of the inclusion was studied by testing the model's prediction ability based on best linear unbiased estimates for fixed and random effects. The material in this study consists of phenotypes from 76 683 animals, of which 4933 animals are genotyped. The heritabilities for front leg conformation, stayability, osteochondrosis and arched back, estimated using the traditional pedigree, were found to be between 0.12 and 0.29. When using the combined genomic and pedigree relationship matrix, the heritabilities were between 0.14 and 0.36. The results show that the combined relationship matrix can be used for the estimation of (co)variance components, and that the predictive ability of the model in this study marginally increases with the inclusion of genomic information.     

基于元共祖的基因组联合育种模拟研究

旨在将整合元共祖的一步法(single-step genomic best linear unbiased prediction with metafounders,MF-SSGBLUP)应用到基因组联合育种中,并与其他经典基因组选择方法进行比较分析。本研究使用QMSim软件模拟3个系谱相互独立的奶牛群体;分别使用广义最小二乘法(generalized least squares,GLS)和原始方法(naïve,NAI)估计不同群体间的祖先关系矩阵Γ;将MF-SSGBLUP、SSGBLUP和BLUP用于3个模拟群体的联合育种,评估各方法在遗传参数和育种值估计方面的差异。在不同遗传力下,GLS所得的Γ矩阵在对角线元素上略低于NAI法,在非对角线元素上没有明显差异,且基因组关系矩阵与基于元共祖构建的亲缘关系矩阵对角线元素相关系数(0.750~0.775)高于基因组关系矩阵与传统的亲缘关系矩阵相关系数(0.508~0.572)。MF-SSGBLUP遗传力估计值(0.138、0.140、0.297和0.298)与当代群体遗传力(0.107和0.296)的偏差小于其余两种方法(0.145、0.173、0.273和0.340),且MF-SSGBLUP估计育种值准确性(0.888~0.908)高于SSGBLUP法(0.863~0.876)和BLUP法(0.854~0.871)。表明,MF-SSGBLUP的遗传参数估计值无偏性更好,估计育种值准确性更高。根据上述模拟数据结果表明,在联合育种中,整合元共祖的基因组选择方法优于其他经典基因组选择方法。