Variance components and breeding values for growth traits from different statistical models.

Estimates of genetic parameters resulting from various analytical models for birth weight (BWT, n = 4,155), 205-d weight (WWT, n = 3,884), and 365-d weight (YWT, n = 3,476) were compared. Data consisted of records for Line 1 Hereford cattle selected for postweaning growth from 1934 to 1989 at ARS-USDA, Miles City, MT. Twelve models were compared. Model 1 included fixed effects of year, sex, age of dam; covariates for birth day and inbreeding coefficients of animal and of dam; and random animal genetic and residual effects. Model 2 was the same as Model 1 but ignored inbreeding coefficients. Model 3 was the same as Model 1 and included random maternal genetic effects with covariance between direct and maternal genetic effects, and maternal permanent environmental effects. Model 4 was the same as Model 3 but ignored inbreeding. Model 5 was the same as Model 1 but with a random sire effect instead of animal genetic effect. Model 6 was the same as Model 5 but ignored inbreeding. Model 7 was a sire model that considered relationships among males. Model 8 was a sire model, assuming sires to be unrelated, but with dam effects as uncorrelated random effects to account for maternal effects. Model 9 was a sire and dam model but with relationships to account for direct and maternal genetic effects; dams also were included as uncorrelated random effects to account for maternal permanent environmental effects. Model 10 was a sire model with maternal grandsire and dam effects all as uncorrelated random effects. Model 11 was a sire and maternal grandsire model, with dams as uncorrelated random effects but with sires and maternal grandsires assumed to be related using male relationships. Model 12 was the same as Model 11 but with all pedigree relationships from the full animal model for sires and maternal grandsires. Rankings on predictions of breeding values were the same regardless of whether inbreeding coefficients for animal and dam were included in the models. Heritability estimates were similar regardless of whether inbreeding effects were in the model. Models 3 and 9 best fit the data for estimation of variances and covariances for direct, maternal genetic, and permanent environmental effects. Other models resulted in changes in ranking for predicted breeding values and for estimates of direct and maternal heritability. Heritability estimates of direct effects were smallest with sire and sire-maternal grandsire models.     

Comparison between genomic predictions using daughter yield deviation and conventional estimated breeding value as response variables

This study compared genomic predictions using conventional estimated breeding values (EBV) and daughter yield deviations (DYD) as response variables based on simulated data. Eight scenarios were simulated in regard to heritability (0.05 and 0.30), number of daughters per sire (30, 100, and unequal numbers with an average of 100 per sire) and numbers of genotyped sires (all or half of sires were genotyped). The simulated genome had a length of 1200 cM with 15,000 equally spaced Single-nucleotide polymorphism (SNP) markers and 500 randomly distributed Quantitative trait locus (QTL). In the simulated scenarios, the EBV approach was as effective as or slightly better than the DYD approach at predicting breeding value, dependent on simulated scenarios and statistical models. Applying a Bayesian common prior model (the same prior distribution of marker effect variance) and a linear mixed model (GBLUP), the EBV and DYD approaches provided similar genomic estimated breeding value (GEBV) reliabilities, except for scenarios with unequal numbers of daughters and half of sires without genotype, for which the EBV approach was superior to the DYD approach (by 1.2 and 2.4%). Using a Bayesian mixture prior model (mixture prior distribution of marker effect variance), the EBV approach resulted in slightly higher reliabilities of GEBV than the DYD approach (by 0.3-3.6% with an average of 1.9%), and more obvious in scenarios with low heritability, small or unequal numbers of daughters, and half of sires without genotype. Moreover, the results showed that the correlation between GEBV and conventional parent average (PA) was lower (corresponding to a relatively larger gain by including PA) when using the DYD approach than when using the EBV approach. Consequently, the two approaches led to similar reliability of an index combining GEBV and PA in most scenarios. These results indicate that EBV can be used as an alternative response variable for genomic prediction.     

Comparison of methods for handling censored records in beef fertility data: simulation study

A simulation study was conducted to compare methods for handling censored records for days to calving in beef cattle data. Days to calving was defined as the time, in days, between when a bull is turned out in the pasture and the subsequent parturition. Simulated data were generated to have data structure and genetic relationships similar to an available field data set. Records were simulated for 33,176 daughters of 4,238 sires. Data were simulated using a mixed linear model that included the fixed effects of contemporary group and sex of calf, linear and quadratic covariates for age at mating, and random effects of animal and residual error. Two methods for handling censored records were evaluated, and two censoring rates of 12 and 20% were applied to assess the influence of higher censoring rates on inferences. Censored records were assigned penalty values on a within-contemporary group basis under the first method (DCPEN). Under the second method (DCSIM), censored records were drawn from their respective predictive distributions. A Bayesian approach via Gibbs sampling was used to estimate variance components and predict breeding values. Posterior means (PM) and standard deviations (SD) of additive genetic variance for DCPEN at 12 and 20% censoring were 23.2 (3.7) and 21.0 (3.6), respectively, whereas the same estimates for DCSIM at 12 and 20% censoring were 23.7(3.3) and 21.9 (3.4), respectively. In all cases, the true value of the genetic variance was within the high posterior density (HPD) interval (95%). The PM (SD) of residual variance for DCPEN at 12 and 20% censoring were 415.7 (4.7) and 440.0 (4.8) respectively, whereas the same estimates for DCSIM at 12 and 20% censoring were 371.0 (4.3) and 365.4 (4.4), respectively. The true value of the residual variance was within the HPD (95%) for DCSIM, but it was outside this interval for DCPEN at both censoring rates, indicating a systematic bias for this parameter. Bayes Factor and Deviance Information Criteria were used for model comparisons, and both criteria indicated the superiority of the DCSIM method. However, little difference was observed between the two methods for correlations between true breeding values and posterior means of animal effects for sires, indicating that no major reranking of sires would be expected. This finding suggests that either censored data handling technique can be successfully used in a genetic evaluation for days to calving.     

Progeny testing sires selected by independent culling levels for below-average birth weight and high yearling weight or by mass selection for high yearling weight.

Breeding values of sires resulting from selection either for reduced birth weight and increased yearling weight (YB, n = 8) or for increased yearling weight alone (YW, n = 9) were compared with each other and with sires representative of the population before selection began (BS, n = 12) using progeny testing. Reference sires (n = 6) connected these Line 1 sires with the Hereford international genetic evaluation. Thirty-five sires produced 525 progeny that were evaluated through weaning. After weaning, 225 steer progeny were individually fed, slaughtered, and carcass data collected. Data were analyzed using restricted maximum likelihood procedures for multiple traits to estimate breeding values for traits measured on the top-cross progeny while simultaneously accounting for selection of the sires. Results of the progeny test substantiate within-line results for traits upon which sires were selected. Breeding values for gestation length were greater in YB sires than in YW sires and were unchanged relative to BS sires. Breeding values for growth rate and feed intake for the YB and YW sires were greater than for BS sires. Predicted breeding values for indicators of fat deposition tended to be greater in YB sires and less in YW sires relative to BS sires, although YB and YW sires had similar breeding values for marbling score. Selection based on easily and routinely measured growth traits, although achieving the intended direct responses, may not favorably affect all components of production efficiency. Further, divergence of selection lines may not be easily anticipated from preexisting parameter estimates, particularly when selection is based on more than one trait.     

Relationship between litters per sow per year sire breeding values and sire progeny means for farrowing rate,removal parity and lifetime born alive

The objective of this study was to determine the relationship between individual sire estimated breeding values (EBV) for litters/sow/year (LSY) and sire progeny means for farrowing rate (FR), removal parity and lifetime born alive (LTBA). Genetic parameters and breeding values were estimated using ASREML. The heritability estimate for LSY was 0.11. When all sires with 10 or more daughters with records were included in the analysis, Spearman rank correlations between the sire's LSY EBV and the sires' daughter means for FR, removal parity and LTBA were 0.49, 0.23 and 0.25 (p < 0.01). The sire EBV for LSY was favourably correlated with sires' daughter means for all three traits. This provides evidence that selecting sires with high EBV for LSY could improve herd FR, removal parity and LTBA. By including LSY as part of the selection criterion, the LTBA may be indirectly improved. The positive genetic correlation between LTBA and LSY may be a result of the improved longevity of sows with greater LSY compared with sows with lower LSY. The relationships between LSY and FR, removal parity and LTBA are strongly supported by the correlations between the sire progeny means for each trait and the sire LSY EBV.     

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.     

Comparisons of Angus, Charolais, Galloway, Hereford, Longhorn, Nellore, Piedmontese, Salers, and Shorthorn breeds for weight, weight adjusted for condition score, height, and condition score of cows

Breed differences for weight (CW), height (CH), and condition score (CS) were estimated from records (n = 12,188) of 2- to 6-yr-old cows (n = 744) from Cycle IV of the U.S. Meat Animal Research Center's Germplasm Evaluation (GPE) Program. Cows were produced from mating Angus and Hereford dams to Angus, Hereford, Charolais, Shorthorn, Galloway, Longhorn, Nellore, Piedmontese, and Salers sires. Samples of Angus and Hereford sires were 1) reference sires born from 1962 through 1970 and 2) 1980s sires born in 1980 through 1987. The mixed model included cow age, season of measurement and their interactions, year of birth, pregnancy-lactation code (PL), and breedgroup as fixed effects for CW and CS. Analyses of weight adjusted for condition score included CS as a linear covariate. The model for CH excluded PL. Random effects were additive genetic and permanent environmental effects associated with the cow. Differences among breed groups were significant (P < 0.05) for all traits and were maintained through maturity with few interchanges in ranking. The order of F1 cows for weight was as follows: Charolais (506 to 635 kg for different ages), Shorthorn and Salers, reciprocal Hereford-Angus (HA) with 1980s sires, Nellore, HA with reference sires, Galloway, Piedmontese, and Longhorn (412 to 525 kg for different ages). Order for height was as follows: Nellore (136 to 140 cm), Charolais, Shorthorn, Salers, HA with 1980s sires, Piedmontese, Longhorn, Galloway and HA with reference sires (126 to 128 cm). Hereford and Angus cows with reference sires were generally lighter than those with 1980s sires. In general, breed differences for height followed those for weight except that F1 Nellore cows were tallest, which may in part be due to Bos taurus-Bos indicus heterosis for size.     

Comparison of methods for handling censored records in beef fertility data: field data

The purpose of this study was to compare methods for handling censored days to calving records in beef cattle data, and verify results of an earlier simulation study. Data were records from natural service matings of 33,176 first-calf females in Australian Angus herds. Three methods for handling censored records were evaluated. Censored records (records on noncalving females) were assigned penalty values on a within-contemporary group basis under the first method (DCPEN). Under the second method (DCSIM), censored records were drawn from their respective predictive truncated normal distributions, whereas censored records were deleted under the third method (DCMISS). Data were analyzed using a mixed linear model that included the fixed effects of contemporary group and sex of calf, linear and quadratic covariates for age at mating, and random effects of animal and residual error. A Bayesian approach via Gibbs sampling was used to estimate variance components and predict breeding values. Posterior means (PM) (SD) of additive genetic variance for DCPEN, DCSIM, and DCMISS were 22.6d2 (4.2d2), 26.1d2 (3.6d2), and 13.5d2 (2.9d2), respectively. The PM (SD) of residual variance for DCPEN, DCSIM, and DCMISS were 431.4d2 (5.0d2), 371.4d2 (4.5d2), and 262.2d2 (3.4d2), respectively. The PM (SD) of heritability for DCPEN, DCSIM, and DCMISS were 0.05 (0.01), 0.07 (0.01), and 0.05 (0.01), respectively. Simulating trait records for noncalving females resulted in similar heritability to the penalty method but lower residual variance. Pearson correlations between posterior means of animal effects for sires with more than 20 daughters with records were 0.99 between DCPEN and DCSIM, 0.77 between DCPEN and DCMISS, and 0.81 between DCSIM and DCMISS. Of the 424 sires ranked in the top 10% and bottom 10% of sires in DCPEN, 91% and 89%, respectively, were also ranked in the top 10% and bottom 10% in DCSIM. Little difference was observed between DCPEN and DCSIM for correlations between posterior means of animal effects for sires, indicating that no major reranking of sires would be expected. This finding suggests little difference between these two censored data handling techniques for use in genetic evaluation of days to calving.     

Sire x herd interactions for weaning weight in beef cattle.

Weaning weight records of 44,357 Australian Angus calves produced by 1,020 sires in 90 herds were used to evaluate the importance of sire x herd interactions. Models fitted fixed effects of contemporary group (herd-year-date of weighing subclass), sex, calf age, and dam age and random effects of sire or of sire and sire x herd interaction using REML. Effects of standardizing the data, including sire relationships and including dam maternal breeding values (MBV) as a covariate were also investigated. Sire x herd interactions were found (P less than .05) in all cases and, in the most complete model, accounted for 3.3% of phenotypic variance. Across-herd heritabilities ranged from .19 to .28. Differential nonrandom mating among herds seemed to occur in the data. Significant sire x herd effects were observed for dam MBV, and adjustment for dam MBV yielded the smallest estimates of interaction variance and across-herd heritability. If sire x herd interactions were due only to genotype x environment interaction, within-herd heritabilities would range from .33 to .49. These estimates are larger than previously reported estimates. Thus, unreported environmental effects common to progeny of individual sires may also be involved in the observed interaction but could not be disentangled from true genotype x environment interaction effects using these data. Results of these analyses suggest that some accommodation of sire x herd interaction effects on weaning weight may be needed in beef cattle genetic evaluations, but a compelling case for development of herd-specific breeding value prediction cannot be made.     

Evaluation of expected response to selection for orthopedic health and performance traits in Hanoverian Warmblood horses

OBJECTIVE: To determine whether selection schemes accounting for orthopedic health traits were compatible with breeding progress in performance parameters in Hanoverian Warmblood horses. ANIMALS: 5,928 horses. PROCEDURE: Relative breeding values (RBVs) were predicted for osseous fragments in fetlock (metacarpo- and metatarsophalangeal) and tarsal joints, deforming arthropathy in tarsal joints, and pathologic changes in distal sesamoid bones. Selection schemes were developed on the basis of total indices for radiographic findings (TIR), dressage (TID), and jumping (TIJ). Response to selection was traced over 2 generations of horses for dressage and jumping ability and all-purpose breeding. Development of mean RBVs and mean total indices in sires and prevalences of orthopedic health traits in their offspring were used to assess response to selection. RESULTS: Giving equal weight toTIR andTID, TIJ, or a combined index of 60% TID and 40% TIJ, 43% to 53% of paternal grandsires and 70% to 82% of descending sires passed selection. In each case, RBVs and total indices increased by as much as 9% in selected sires, when compared with all sires, and prevalences of orthopedic health traits in offspring of selected sires decreased relatively by as much as 16%. When selection was exclusively based on TID, TIJ, or TID and TIJ, percentages of selected sires were 44% to 66% in the first and 73% to 84% in the second generation and TID and TIJ increased by 9% to 10% and 19% to 23%, respectively. CONCLUSIONS AND CLINICAL RELEVANCE: Compared with exclusively performance-based selection, percentages of selected sires changed slightly and breeding progress in TID, TIJ, or TID and TIJ was only slightly decreased; however, prevalences of orthopedic health traits decreased in offspring of TIR-selected sires.     

Estimated breeding values for meat characteristics of crossbred cattle with an animal model.

Longissimus muscle area, shear force measure, and sensory panel scores for flavor, juiciness, and tenderness, and marbling score were obtained from 682 steer carcasses, resulting from crosses among five Bos taurus and Bos indicus breeds. The single-trait model used included birth year and as covariates breed fractions, weaning age, and days on feed. The numerator relationship matrix was for 1,350 animals (682 steers, 74 pure breed and 52 F1-cross sires and 542 dams). The coefficient matrix was inverted to examine standard errors of prediction. Estimated breeding value is the sum of the estimate of genetic deviation and the weighted (fractions) sum of estimates of breed effects. Heritabilities used in estimating breeding values were .62, .06, .05, .11, .05, and .43 for longissimus muscle area, shear force, flavor, juiciness, tenderness, and marbling score. Sires within a breed or crossbred group tended to rank similarly due to large differences among breed effects (e.g., the six Sahiwal sires ranked in the highest six places for shear force). These results illustrate that for traits with large breed differences, selection of the proper breed should be done before selection within that breed.     

Pregnancy rate and first-service conception rate in Angus heifers

The objective of this project was to determine the genetic control of conception rate, or pregnancy percentage in Angus beef heifers. Producers from 6 herds in 5 states provided 3,144 heifer records that included breeding dates, breeding contemporary groups, service sires, and pregnancy check information. Two hundred fourteen sires of the heifers were represented; with 104 sires having less than 5 progeny, and 14 sires having greater than 50 progeny. These data were combined with performance and pedigree information, including actual and adjusted birth weights, weaning weights, and yearling weights, from the American Angus Association database. Heifer pregnancy rate varied from 75 to 95% between herds, and from 65 to 100% between sires, with an overall pregnancy rate of 93%, measured as the percentage of heifers pregnant at pregnancy check after the breeding season. Pregnancy was analyzed as a threshold trait with an underlying continuous distribution. A generalized linear animal model, using a relationship matrix, was fitted. This model included the fixed effects of contemporary group, age of dam, and first AI service sire, and the covariates of heifer age at the beginning of breeding, adjusted birth weight, adjusted weaning weight, and adjusted yearling weight. The relationship matrix included 4 generations of pedigree. The heritability of pregnancy and first-service conception rates on the underlying scale was 0.13 +/- 0.07 and 0.03 +/- 0.03, respectively. Estimated breeding values for pregnancy rate on the observed scale ranged from -0.02 to 0.05 for sires of heifers. Including growth traits with pregnancy rate as 2-trait analyses did not change the heritability of pregnancy rate. As expected for a reproductive trait, the heritability of pregnancy rate was low. Because of its low heritability, genetic improvement in fertility by selection on heifer pregnancy rate would be expected to be slow.     

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.     

Genetic and environmental trends for litter size in swine

Best linear unbiased predictors (BLUP) of breeding values for additive direct and additive maternal genetic effects were estimated from 3,944 purebred Yorkshire and Landrace first-parity litters recorded on the Quebec Record of Performance Sow Productivity Program and born between 1977 and 1987. Breeding values for gilts, dams, and sires were estimated using an individual animal model for measures of litter size of total number born (NOBN), number born alive (NOBA), and number weaned (NOWN). Environmental trends were estimated from average herd-year solutions, and genetic trends were estimated by regression of estimated breeding value on year of birth. Environmental trends were positive for all traits in both breeds but were significant only for NOWN in Landrace (.051 +/- .021 pigs/yr). Genetic trends were very small but were mainly negative for direct breeding value and combined direct and maternal breeding value. Significant estimates of genetic trends (P less than .05) were observed only within the Yorkshire breed, and these ranged from -.012 +/- .004 to .004 +/- .002 pigs/yr.     

Three statistical approaches for genetic analysis of disease resistance to vibriosis in Atlantic cod (Gadus morhua L.)

The objective of this study was to estimate and compare variance components and sire breeding values for disease resistance to vibriosis in Atlantic cod (Gadus morhua L.) using 3 statistical approaches. A total of 3,576 individually tagged juvenile cod from 50 full-sib families were infected with Vibrio anguillarum, which causes vibriosis, a frequently reported disease in cod aquaculture. The experimental fish were progeny of captured wild cod from populations of southern coastal cod (POP1), and northern coastal cod and northeast Arctic cod (combined as POP2 in the genetic analyses). Fish were randomly assigned to 1 of 3 test tanks, and daily mortality was recorded until the termination of the experiment at d 31 postinfection. Variance components were estimated separately for the 2 populations using a Cox regression model, univariate linear model, and a linear model that accounted for censoring. With all approaches, the additive genetic sire variance estimated from POP1 was greater than for POP2. Heritability estimates across populations varied from 0.08 to 0.17 depending on the method used. The Cox regression model and univariate linear model resulted in greater heritability estimates for POP1 (0.10 and 0.16) than for POP2 (0.08 and 0.13), whereas the contrary was true with a linear model that accounted for censoring (0.17 vs. 0.14). The predicted breeding values for the sires from the 3 approaches were highly correlated (0.97 to 0.99). This is likely due to the fact that censoring only occurred at the end of the test; i.e., observations of the most resistant fish were censored. The considerable genetic variation found in this study suggests that vibriosis resistance may be improved through selective breeding. The univariate linear model, even without censoring of the data, was robust for the estimation of breeding values using the present data. Therefore, inclusion of vibriosis resistance in the multivariate linear estimation of breeding values for the traits of economic importance in Atlantic cod seems appropriate.     

Approaches to the management of inbreeding and relationship in the German Holstein dairy cattle population

The aim of this study was to estimate the current level of inbreeding in the German cow population and for bull dams born in Germany, to find out sires most related to different subsets of their breed and to demonstrate the negative effect of homozygosity in the case of complex vertebral malformation (CVM). Further on, the application of optimum genetic contribution (OGC) theory for the selection of bull dams and bull sires in different breeding scenarios was investigated. Levels of inbreeding for the cow population were in a low range from 0.97% to 1.70% evaluating birth years from 1996 to 1999 in a total dataset of 244,427 registered Holstein cows. The inbreeding coefficient of 8030 bull dams was much higher, i.e. 3.71%, for the birth year 1999. Increases in inbreeding of 0.19% per year indicated an effective population size of only 52 animals. Individual sires like R.O.R.A. Elevation and Hannoverhill Starbuck were highly related to potential bull dams with coefficients of relationship of 13.4% and 12.9%, respectively, whereas P.F. Arlinda Chief (16.3%) and Carlin-M Ivanhoe Bell (16.1%) were highest related to the best available AI sires. Coefficients of relationship were calculated by classes of estimated breeding values (EBV) for production traits showing highest values above 7% in the two highest EBV-classes. The optimum genetic contribution theory using official EBVs and approximative, for zero inbreeding corrected EBVs, was applied for elite matings in a breeding program embracing 30 young bulls per year to find the optimal allocations of bull sires and bull dams. Compared with the actual breeding program applied in practice, OGC-theory has the potential to increase genetic gain under the same constraint for the increase of average relationship by 13.1%. A more relaxed constraint on increase in inbreeding allowed even higher expected genetic gain whereas a more severe constraint resulted in more equal contributions of selected bull sires. Contributions from 21 selected bull sires and 30 selected bull dams for a scenario at 5% constrained relationship were used to develop a specific mating plan to minimise inbreeding in the short term in the following generation applying a simulated annealing algorithm. The expected coefficient of inbreeding of progeny was 66.3% less then the one resulting from random mating. Mating programs can address inbreeding concerns on the farm, at least in the short term, but long-term control of inbreeding in a dairy population requires consideration of relationships between young bulls entering AI progeny test programs. Significantly better EBVs of CVM-free bulls compared with CVM-carriers for the paternal fertility justify the application of OGC for elite matings.     

Optimum contribution selection using traditional best linear unbiased prediction and genomic breeding values in aquaculture breeding schemes

The aim of this study was to compare genetic gain for a traditional aquaculture sib breeding scheme with breeding values based on phenotypic data (TBLUP) with a breeding scheme with genome-wide (GW) breeding values. Both breeding schemes were closed nuclei with discrete generations modeled by stochastic simulation. Optimum contribution selection was applied to restrict pedigree-based inbreeding to either 0.5 or 1% per generation. There were 1,000 selection candidates and a sib test group of either 4,000 or 8,000 fish. The number of selected dams and sires to create full sib families in each generation was determined from the optimum contribution selection method. True breeding values for a trait were simulated by summing the number of each QTL allele and the true effect of each of the 1,000 simulated QTL. Breeding values in TBLUP were predicted from phenotypic and pedigree information, whereas genomic breeding values were computed from genetic markers whose effects were estimated using a genomic BLUP model. In generation 5, genetic gain was 70 and 74% greater for the GW scheme than for the TBLUP scheme for inbreeding rates of 0.5 and 1%. The reduction in genetic variance was, however, greater for the GW scheme than for the TBLUP scheme due to fixation of some QTL. As expected, accuracy of selection increased with increasing heritability (e.g., from 0.77 with a heritability of 0.2 to 0.87 with a heritability of 0.6 for GW, and from 0.53 and 0.58 for TBLUP in generation 5 with sib information only). When the trait was measured on the selection candidate compared with only on sibs and the heritability was 0.4, accuracy increased from 0.55 to 0.69 for TBLUP and from 0.83 to 0.86 for GW. The number of selected sires to get the desired rate of inbreeding was in general less in GW than in TBLUP and was 33 for GW and 83 for TBLUP (rate of inbreeding 1% and heritability 0.4). With truncation selection, genetic gain for the scheme with GW breeding values was nearly twice as large as a scheme with traditional BLUP breeding values. The results indicate that the benefits of applying GW breeding values compared with TBLUP are reduced when contributions are optimized. In conclusion, genetic gain in aquaculture breeding schemes with optimized contributions can increase by as much as 81% by applying genome-wide breeding values compared with traditional BLUP breeding values.     

Effectiveness of adjusting for heterogeneity of variance in genetic evaluation of Japanese Black cattle

Heterogeneity of variance among subclasses of an effect is a potential source of bias in genetic evaluation. The objectives of this study were to quantify the heterogeneity of variance in carcass weight in Japanese Black cattle, to develop an adjustment method to account for the heterogeneity, and to evaluate the effectiveness of the method. A total of 96,950 records were collected from steers and heifers slaughtered from 1997 to 2005. These records were grouped into 2,767 farm-market-year-sex subclasses. Fourteen log-linear models for the variances were set up to estimate the heterogeneous phenotypic variances within subclasses. Schwarz's Bayesian information criterion was used for model selection. The preadjustment of records to a baseline variance was based on maximum likelihood estimates obtained from the selected model. As a result of adjustment, the SD, the CV, and the Gini coefficient for the phenotypic variance decreased by 68.6, 69.8, and 70.1%, respectively. When the top 5% of sires and top 1% of dams were selected, Spearman's rank correlation coefficients between the adjusted and unadjusted data were 0.95 for the selected sires and 0.78 for the selected dams. The effectiveness of the adjustment was evaluated in terms of the ability to predict breeding values, using the results of the successive genetic evaluations. Mean squared error between the parent averages and actual predicted values of the genetic merit for the sires whose progeny had a carcass record only from 2003 to 2005 was significantly reduced by the adjustment (P < 0.05). The results suggest that the genetic evaluation becomes more accurate by adjusting the data using the procedure developed in this study.     

Breeding performance of beef bulls assigned to two-sire cow groups on pasture

Angus, Polled Hereford and Santa Gertrudis bulls from ages 1 through 5 and 7 yr were assigned to 26 two-sire breeding groups. Each year, straightbred and crossbred cows of these breeds were allotted at random within breed composition, age of dam and calving date to breeding groups on pasture. Sires within each breeding group or pair were the same age at breeding and were two of the three breeds of sires. Neither calving rate nor the proportion of calves born by one vs the other sire in the two-sire breeding groups was affected by sire age among breeding groups. For a given breed, there was no uniformity among the sires in the proportion of calves they sired in their two-sire breeding groups. The proportion of calves born for the 26 sire pairs averaged .64 vs .36 (SE = 0.4 for either high or low value) for one vs the other sire in a sire pair with no indication that calving rate was affected by unequal proportions of calves by sires within sire pairs. Cows calved significantly earlier in the calving period (b = -.775 +/- .127) as calving rate increased among sire pairs. The number of days from the start of the breeding period to calf birth was affected by differences between sires in sire pairs for 8 of the 26 pairs, but there were no significant differences due to sire pair or breed of sire because of interaction between these two variables.     

Estimates of genetic parameters for live animal ultrasound, actual carcass data, and growth traits in beef cattle

Growth and carcass measurements were made on 2,411 Hereford steers slaughtered at a constant weight from a designed reference sire program involving 137 sires. A second data set consisted of ultrasound measures of backfat (USFAT) and longissimus muscle area (USREA) from 3,482 yearling Hereford cattle representing 441 sires. Restricted maximum likelihood procedures were used to estimate genetic parameters among carcass traits and live animal weight traits from these two separate data sets. Heritability estimates for the slaughter weight constant steer carcass backfat (FAT) and longissimus muscle area (REA) were .49 and .46, respectively. In addition, FAT had a negative genetic correlation with REA (-.37), weaning weight (-.28), and yearling weight (-.13) but positive with marbling (.19) and carcass weight (.36). Marbling was moderately heritable (.35) and highly correlated with total postweaning average daily gain (.54) and feedlot relative growth rate (.62). Heritability estimates for weight constant USFAT and USREA were .26 and .25, respectively. The genetic correlation between weight constant USFAT and USREA was positive (.39), indicating that in these young animals USFAT does not seem to be an indication of maturity. Mean USFAT measures and variability were small (.48 +/- .17 cm, n = 3,482). Results indicate that carcass fat on slaughter steers and ultrasound measures of backfat on young breeding animals may have different relationships with growth and muscling. These relationships need to be explored before wide scale selection based on ultrasound is implemented.