|
|
||||||||
ANIMAL GENETICS |
ARS, USDA, Fort Keogh Livestock and Range Research Laboratory, Miles City, MT 59301
| Abstract |
|---|
|
|
|---|
Key Words: Beef Cattle Cow Indexes Genetic Gain Selection Responses
| Introduction |
|---|
|
|
|---|
To justify using a ratio to express information about two quantities, three conditions should be satisfied: 1) the relationship between components is linear; 2) the intercept of the linear regression of numerator on denominator should be the origin; and 3) the variance of the ratio should increase with increasing values of numerator and denominator (Weil, 1962
). In addition, a ratio of two normally distributed variables is not normally distributed (Fieller, 1932
; Hinkley, 1969
). Thus, concern was raised for the traditional concept of heritability of a ratio (Gunsett, 1987
).
Despite theoretical issues with the ratio of calf weaning weight to cow weight, continuing interest in its use in selection and a need to improve the efficiency of beef production points to a need for genetic evaluation of this ratio and its components. Thus, the objectives of this research were to estimate genetic parameters for the ratio, its components, and other growth traits, and to evaluate responses to selection based on the ratio.
| Materials and Methods |
|---|
|
|
|---|
Red Angus,
Charolais, and
Tarentaise. Animals of this breed composition were inter se mated in all subsequent generations. No additional germplasm has been introduced into CGC since the foundation matings. Thus, direct breed effects were stabilized in the F2 generation, direct heterosis and maternal breed effects were stabilized in the F3 generation, and maternal heterosis effects were stabilized in the F4 generation. Additional details describing the formation of the CGC population can be found in Newman et al. (1993a)
In 1988, the population was divided randomly into three lines. Bulls for the first line were selected based on the index I = 365-d weight 3.2(birth weight) (Dickerson et al., 1974
). Response to selection for the index proposed by Dickerson et al. (1974)
has been previously reported (MacNeil, 2003
). Bulls for the second line were selected based on 100 times the ratio of their age-of-dam-adjusted 200-d weight to the coincident mature equivalent weight of their dam. Selection decisions were based on phenotypic performance within year. Virtually all selection pressure was applied to males, and most females were exposed for breeding as yearlings. Bulls for the third line were selected at random. All bulls were required to pass a breeding soundness examination as yearlings before being used for breeding.
In 1988 and 1989, the process of stabilizing the genetic composition of CGC was ongoing. Only cows belonging to generations subsequent to the F2 contributed to the selection experiment. Hence, the numbers of calves in the control and ratio selection lines were less than the numbers of calves in subsequent years (Table 1
). In 1989 and 1990, yearling heifers and bulls were sent to the U.S. Sheep Exp. Stn., Dubois, ID, for use with sheep in a mixed-species grazing trial. Mating occurred during the trial and all animals were returned to Miles City in late August. Thus, yearling heifers were mated to different bulls than the remainder of the cows in those years. Except during the mixed-species grazing trial, mating assignments were subject to the constraint that approximately equal numbers of like-aged females were assigned to each sire. In 1995 to 1997, a sample of cows was removed from the experiment (primarily from the control line) to initiate a project for mapping QTL (Mac-Neil and Grosz, 2002
). Except as described below, four yearling bulls were used as sires in each line during each breeding season, with one sire from each line, selected at random, used the next year. In the control line: nine yearling bulls were used in the 1989 and 1990 breeding seasons and the same six bulls were used in the 1996 to 1998 breeding seasons; in the 1999 breeding season, four of the six bulls used in 1998 were used again, and the other two bulls used in 1998 were replaced by sons. In the ratio selection line, nine yearling bulls were used in the 1989 and 1990 breeding seasons, and in the 1992 and 1996 to 1999 breeding seasons, five yearling bulls were used each year. Increasing the number of bulls used per year, while decreasing the number of cows, allowed the effective population size of the control line to remain approximately equal to that of the ratio selection line during the latter years of this research. Within line, bulls selected for breeding were assigned to cows to avoid inbreeding of the progeny produced in the next generation.
|
A 45-d breeding season began on approximately June 15 of each year. Cows were exposed for breeding in single-sire pastures until approximately August 1, with the assignment to breeding pastures at random with respect to line. After the breeding season, cows were grouped by geographic proximity of the breeding herds into two larger herds, and grazed native rangeland until vegetation was covered by snow. Cows were weighed and pregnancy tested when their calves were weaned in early October. Nonpregnant females were culled. Cows and heifers were managed separately through winter, with approximately 1.4 kg of alfalfa pellets per animal provided daily when forage was available and hay fed when the forage was covered by snow. In preparation for calving, cows were moved to small pastures and fed approximately 9 kg of alfalfa hay per cow daily. As calving approached, first-calf heifers were observed periodically throughout the day. Older cows were observed only during daylight. Cow weights were adjusted to a mature (5-yr-old) equivalent basis using multiplicative adjustment factors derived from the data first in 1987 and subsequently in 1990. The adjustment factors were 1.197, 1.110, and 1.045 for 2-, 3-, and 4-yr-old cows, respectively.
Most calves were born during April and May of each year. Male calves were not castrated, and no creep feed was offered to the calves. Calves were weaned at an average age of approximately 180 d. After weaning, the calves were moved to a feedlot for a 140-d evaluation of postweaning growth. Before the test period began, the calves were allowed a pretest adjustment period of 2 to 4 wk following weaning. Bull calves were fed a diet that was formulated to support an average growth rate of 1.4 kg/d. The energy density of the diet was approximately 2.7 Mcal of ME/kg of DM, and the CP content was approximately 12% of DM. Primary ingredients in the diet were corn silage, barley, and a protein and mineral supplement. Heifer calves were fed to gain 0.8 kg/d. The approximate energy density and CP content of the diet were 2.4 Mcal of ME/kg of DM and 9% of DM, respectively. The calves were weighed on Monday and Wednesday of the weeks at the beginning and end of the postweaning test period, and the pairs of weights were averaged to calculate initial and final weights, respectively. Birth weight and gain from birth to weaning were adjusted for differences in age of dam following the results of Newman et al. (1993b)
. Weaning and yearling weights were then adjusted to age-constant bases as recommended by the BIF (1996)
.
Selection applied and turnover of generations were monitored for each line, separately. Individual selection differentials (ID) were computed for the ratio and component traits within year-sex-line subclasses based on the adjusted phenotypes. Cumulative selection differentials were calculated following the methodology of Newman et al. (1973)
, with differences in selection applied between sires and dams summarized by year-line subclass following Koch et al. (1994)
. Performance records were not available for the original purchased females and their contemporaries or for the purebred bulls used by AI and their contemporaries. Therefore, in calculating the cumulative selection differentials, the deviation of a phenotype of an individual from the contemporary group mean was assumed to be zero for all animals born before 1980. For an animal, the cumulative selection differential (CS) was calculated recursively as:
![]() |
where CSs = cumulative sire differential and CSd = cumulative dam differential. The generation number (GN) of each animal was calculated similarly as:
![]() |
where GNs and GNd = generation numbers of the sire and dam, respectively (Brinks et al., 1964
). Generation number also was assumed to be zero for all animals born before 1980.
Preliminary analyses to derive estimates of parameters of the prior distribution for a Bayesian analysis of these data were conducted using multiple-trait, derivative-free REML (Smith and Graser, 1986
; Graser et al., 1987
) methods as implemented by Boldman et al. (1995)
. Genetic parameters were estimated from phenotypes recorded from all animals born between 1980 and 2000. Pedigree information was extended to the founding pure-breds, which were assumed to be unrelated. Earlier results indicated that year of birth, sex of calf, age of dam, and a variety of two-factor interactions among these traits affected birth, 200-d, and 365-d weights. Hence, contemporary groups for these traits were formed as year, sex, and age of dam subclasses for further analyses. Linear models for birth weight (kg; n = 5,083), 200-d weight (kg; n = 4,902), 365-d weight (kg; n = 4,626), and the ratio (100ªkg/kg; n = 4,184) were similar and included fixed contemporary groups and correlated random direct and maternal additive effects and uncorrelated random maternal permanent environmental effects of the dams. Cows ranged in age from 2 to 11 yr, with an average age of 3.6 yr. For cow weight (kg), contemporary groups were formed as year-age subclasses, with cows at least 5 yr of age being grouped together. A total of 1,433 cows had at least one recorded weight, and there were 4,375 observations over time. The linear model for cow weight included fixed contemporary group effects, random additive direct effects, and uncorrelated random permanent environmental effects to allow for repeated records of the cow. The same contemporary group effects were fitted in models for the ratio of 200-d weight to cow weight and 200-d weight. Multiple trait analyses, constructed using the models described above, allowed for correlations among the additive effects, the uncorrelated random effects, and the residual errors. This multiple-trait, derivative-free REML method involves maximizing the likelihood function (
) given the data and is equivalent to minimizing 2 log
. Each analysis was assumed to have converged to a global maximum when the variance of 2 log
in the simplex was less than 1010, and the scaled parameter estimates changed by less than 0.01 in a reanalysis of the data using an updated starting value.
Final data analyses were conducted using a multiple-trait Gibbs sampler for animal models (VanTassell and Van Vleck, 1996
). The linear models used for each trait were as described above. Inverted Wishart distributions, with shape parameter = 9, were used as prior distributions of the (co)variance components. For a univariate analysis of the ratio of weaning weight to cow weight and for a bivariate analysis of weaning weight and cow weight, initial Gibbs samples were obtained for chains run 20,000 rounds, saving every 20 rounds, and without discarding earlier rounds for burn-in. The Gibbsit program of Raftery and Lewis (1996)
was used to evaluate the length of the burn-in period and thinning interval needed to obtain stationary chains of independent samples and to estimate the total number of Gibbs samples required for the cumulative distribution function of the 0.025 quantile to be estimated within ± 0.0125, with a probability of 0.95 for all genetic parameters. The recommended burn-in periods were from 40 to 160 rounds. The thinning interval recommended by Gibbsit was consistently less than or equal to 20 rounds; however, it has been shown that the thinning intervals recommended by Gibbsit are not sufficiently large (Van Tassell and Van Vleck, 1996
). Thus, final analyses were conducted using a conservative burn-in of 5,000 rounds, 50,000 rounds of post-burn-in Gibbs sampling, and a thinning interval of 100 rounds. In summarizing the Gibbs samples, results are reported as mean ± SD of the posterior Gibbs samples.
| Results and Discussion |
|---|
|
|
|---|
Selection differentials for calves born in the early years of this experiment resulted from selection applied during formation of the CGC composite population and stabilization of its breed composition. Initial differences in selection differentials between the ratio and control lines primarily reflect differences in individual deviations of the selected sires (Table 1
). In the line selected for increased values of the ratio of calf weaning weight to weight of its dam at weaning, the total (sire + dam) indirect cumulative selection differential for 200-d weight increased by 153 kg (6.4 SD) from 1989 to 2000, whereas the direct cumulative selection differential increased by 28.67 ratio units (4.7 SD) over the same time period. However, the pattern of sire cumulative selection differentials for cow weight is remarkable and documents inconsistency in the indirect selection applied to cow weight through direct selection on the ratio of calf weaning weight to weight of its dam at weaning. Trends in the cumulative selection differentials for 200-d weight, cow weight, and their ratio were all significantly positive across years in the control line. Rogers et al. (2004)
documented the increased risks that lighter cows and those that experience dystocia have of not being pregnant, and therefore being culled from this population. Thus, these positive trends of the selection differentials in the control may be compromised by small females failing to breed and/or calve successfully and from the relationship between the numerator and denominator of the ratio.
The mean estimate of phenotypic variance of the ratio of calf 200-d weight to cow weight in the posterior distribution was 37.6. Mean estimates of direct and maternal heritability for the ratio were 0.20 ± 0.03 and 0.58 ± 0.05, respectively. The estimated posterior mean of the genetic correlation between direct and maternal effects was 0.77 ± 0.04. Permanent environmental effects due to females and residual effects, on average, accounted for 23 ± 3% and 25 ± 2% of phenotypic variance, respectively. Distributions of the Gibbs sample estimates for each of these parameters are shown in Figure 1
. No comparable partitioning of the variance of the ratio of calf weaning weight to cow weight was found in the literature. However, the strong genetic antagonism between direct and maternal genetic effects found here and the average posterior estimate of total heritability (Willham, 1972
) of 0.10 ± 0.02 foretell of challenges in the genetic improvement of this trait through selection based on phenotype.
|
|
The regression of 200-d weight on cow weight after adjusting both weights for contemporary group effects was nonlinear (P = 0.05) and the y-intercept of the regression line was positive (P < 0.01), complicating interpretation of the ratio. Figure 3
illustrates the relationship between cow weight and 200-d weight of her calf as implied by the ratio vis a vis the relationship derived from the data via regression. Estimates of genetic trends in 200-d weight and cow weight for the control line and in response to selection on the ratio of 200-d weight to cow weight are shown in Figures 4
and 5
, respectively. In the control line, the regressions of breeding value on generation number for direct and maternal 200-d weight were not significant (1.28 ± 0.91 and 0.62 ± 0.92 kg/ generation, respectively), and the regression of direct breeding value for cow weight on generation number was negative 5.72 ± 2.80 kg/generation. In the line selected for the ratio of calf 200-d weight to cow weight at weaning, the regressions of breeding value on generation number for direct and maternal effects on 200-d weight were 1.43 ± 0.79 and 2.90 ± 0.80 kg/generation, respectively, and the genetic trend in cow weight was not different from zero (2.79 ± 2.43 kg/generation).
|
|
|
To further examine the relationships between the ratio and it components, additional bivariate analyses were conducted with the ratio and its numerator and denominator. When the relationship between the ratio and its 200-d weight was examined in these data, correlations between direct, permanent environmental, and residual effects were all large (0.95 ± 0.01, 0.82 ± 0.03, and 0.79 ± 0.02, respectively). The correlation between maternal effects was smaller, 0.53 ± 0.6, perhaps due in part to the negative covariance of maternal genetic effects on 200-d weight and direct effects on cow weight. The posterior mean phenotypic correlation between the ratio and its numerator was 0.62. Correlations of direct effects on cow weight with both direct and maternal effects on the ratio were 0.87 ± 0.02 and 0.80 ± 0.03, respectively. Likewise, the correlations between permanent and temporary environmental effects on cow weight and the ratio were also large and negative: 0.81 ± 0.04 and 0.55 ± 0.02, respectively. The correlation between direct effects on cow weight and the ratio, largely reflective of the corresponding correlation between direct effects on cow weight and 200-d weight, was 0.82 ± 0.02. Thus, the posterior mean phenotypic correlation between cow weight and the ratio was 0.62, which is approximately equivalent to that predicted following Sutherland (1965)
.
The mean estimate of phenotypic variance for birth weight in the posterior distribution was 32.0 kg2. Estimates of heritability for direct and maternal genetic effects on birth weight were 0.46 ± 0.04 and 0.10 ± 0.02, somewhat greater and less than the corresponding average estimates from the literature reported by Koots et al. (1994a)
, and similar to those reported for this population by MacNeil et al. (2003)
. Genetic correlations between direct effects on birth weight and both direct and maternal effects on the ratio were large: 0.72 ± 0.05 and 0.60 ± 0.05, respectively. Genetic correlations of maternal genetic effects on birth weight with direct and maternal genetic effects on the ratio and with direct genetic effects on birth weight were small, 0.15 ± 0.11, 0.03 ± 0.10, and 0.05 ± 0.11, respectively. As in Mac-Neil et al. (2003)
, permanent environmental effects on birth weight due to dams were essentially nil, and temporary environmental effects accounted for 45 ± 3% of total variance. Correlations of permanent and temporary environmental effects on the ratio and birth weight were 0.19 ± 0.24 and 0.18 ± 0.04, respectively. Genetic trends in birth weight for the control line and in response to selection on the ratio of 200-d weight of the calf to the weight of its dam at weaning are shown in Figure 6
. Only the regression of direct genetic effects on generation number in the selected line (0.13 ± 0.06 kg/generation) differed from zero.
|
|
The ratio of calf weaning weight to cow weight is likely not the contemporary phenotype of choice on which to base a selection program. Its use arose during an era when genetic improvement programs were based on a culture of only measuring weights and trying to extract the most information from them. The ratio undoubtedly had and continues to have some intuitive appeal due to its ease of measurement and perceived relationship with efficiency. In this research, phenotypic selection for the ratio of a calfs 200-d weight to the mature equivalent weight of its dam resulted in a significant genetic trend being established for maternal effects on the ratio. In addition, significant correlated responses were observed in direct effects on birth weight, direct and maternal effects on 200-d weight, and maternal effects on 365-d weight. No significant genetic trend in mature equivalent cow weight was observed. However, a significant negative genetic trend in direct effects on cow weight was observed in the control line, and if this trend resulted from natural selection, rather than drift, it might be inferred that selection for the ratio would affect cow weight positively.
| Implications |
|---|
|
|
|---|
| Footnotes |
|---|
2 Mention of a proprietary product does not constitute a guarantee or warranty of the product by USDA, Montana Agric. Exp. Sta. or the authors and does not imply its approval to the exclusion of other products that may also be suitable. ![]()
3 Several people, including W. L. Reynolds, J. J. Urick, R. A. Bellows, R. E. Short, R. B. Staigmiller, B. W. Knapp, and S. Newman, made significant contributions to the design and/or initial implementation of this research. Their contributions are greatly appreciated. ![]()
4 Correspondence: 243 Fort Keogh Rd. (phone: 406-232-8213; fax: 406-232-8209; e-mail: mike{at}larrl.ars.usda.gov).
Received for publication November 12, 2004. Accepted for publication January 5, 2005.
| Literature Cited |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
J. L. Williams, D. J. Garrick, and S. E. Speidel Reducing bias in maintenance energy expected progeny difference by accounting for selection on weaning and yearling weights J Anim Sci, May 1, 2009; 87(5): 1628 - 1637. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. L. Vallet and B. A. Freking Differences in placental structure during gestation associated with large and small pig fetuses J Anim Sci, December 1, 2007; 85(12): 3267 - 3275. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. G. Riley, S. W. Coleman, C. C. Chase Jr., T. A. Olson, and A. C. Hammond Genetic parameters for body weight, hip height, and the ratio of weight to hip height from random regression analyses of Brahman feedlot cattle J Anim Sci, January 1, 2007; 85(1): 42 - 52. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |