|
|
||||||||
ANIMAL GENETICS |


,
,3
* Department of Animal Science, University of Nebraska, Lincoln 68583-0908; and ARS, USDA, Roman L. Hruska U.S. Meat Animal Research Center,
and
Clay Center, NE 68933 and
and
Lincoln, NE 68583-0908
| Abstract |
|---|
|
|
|---|
Key Words: Beef Cattle Longevity Selection
| Introduction |
|---|
|
|
|---|
Selection for length of productive life would be practical only with indicators of length of productive life which can be obtained early in life and which show genetic variation (Burnside and Wilton, 1970
). The importance of length of productive life for economic value was outlined by Rendel and Robertson (1950)
. Longer productive life may increase profits by 1) decreasing annual cost of replacements of cows, 2) increasing herd production through an increase in proportion of cows in higher-producing age groups, 3) decreasing the number of replacements to be reared, and therefore allowing an increase in size of the productive herd, and 4) increasing voluntary culling.
Weight of weaned calves per female accumulated over a lifetime is a comprehensive measure of fertility (pregnancy and calving rate), maternal ability (weaning rate), milking capacity (weaning weight), and cow survival (ability of the cow to delay culling or death), which has been proposed as an integral measurement of lifetime production (Tanida et al., 1988
).
The objectives of this study were to 1) estimate genetic parameters for six measures of length of productive life and three measures of lifetime production evaluated to 6 y after calving, and 2) estimate genetic and environmental trends for these measures of length of productive life and lifetime production for a control and three selected lines of Hereford cows.
| Materials and Methods |
|---|
|
|
|---|
The Population
The three selection lines were established in 1960 by randomly allocating 325 cows from 14 Hereford herds to the weaning weight line (WWL), the yearling weight line (YWL), and the index line based on yearling weight and muscle score (IXL). The 42 foundation sires used from 1957 to 1963 were from 11 of the same 14 herds as the cows and from two other herds (Koch et al., 1974a
,b
; 1994
).
The cattle were maintained at Fort Robinson Beef Cattle Research Station, Crawford, NE, until 1971, when they were moved to USMARC, near Clay Center, NE. The projected herd size of 150 breeding females per line was reached by 1964 and maintained until the end of the experiment. Approximately 225 foundation cows and other cows from the three selection lines that had been replaced were artificially inseminated with semen from seven of the foundation sires from 1968 through 1971 to provide the basis for a control line. In 1971, the control line (CTL) was established from 20 representative sons and heifer calves from matings with the 225 cows (Koch et al., 1974a
,b
; 1994
).
Selection in WWL was based on weight adjusted to 200 d of age. Selection in YWL was based on weight at 452 d (approximately 15 mo of age) for bulls and at 550 d (approximately 18 mo of age) for heifers. Selection in IXL was based on an index giving equal emphasis to muscle score and yearling weight when both were expressed in standard measure. Selection of heifers in IXL from birth years 1960 through 1965 was based on yearling weight alone. Originally, only bulls were evaluated for muscle score, but beginning in the 1966 birth year, heifers were also evaluated for muscle score and were selected for an index of muscle score and yearling weight until the end of the experiment (Koch et al., 1974a
,b
; 1994
). Selection continued through matings to produce the 1982 calf crop.
Bulls were selected at 2 yr of age. Through 1970, two bulls were chosen from each year of birth and used to sire calves when 3, 4, and 5 yr old. After 1970, three 2-yr-old bulls were selected each year to be used for 2 yr (i.e., they sired calves when 3 and 4 yr of age. Bulls were removed from service early only because of breeding unsoundness. To minimize inbreeding, no more than two sons of a given sire or dam were selected (Koch et al., 1974a
,b
; Koch et al., 1994
). Heifers born in 1964 and later were bred to calve at 2 yr of age. All heifers were exposed to bulls. Selection was practiced only among those that were pregnant. Through 1970, the top 25 heifers from each line were selected based on their line. After 1970, the top 35 heifers were selected. After 1973, essentially all pregnant heifers were kept in the herd. Cows were removed without regard to progeny performance based on the following criteria (Koch et al., 1974a
,b
; 1994
): 1) not pregnant at weaning time; 2) serious unsoundness (e.g., cancer eye, chronic bloat, broken teeth); 3) failure to raise a live calf for two consecutive years; and 4) older cows were removed if, after Criteria 1 through 3 were applied, too many cows were still left.
Exceptions for Criteria 1 and 3 were made only when additional cows were needed to maintain the herd size.
For each breeding season, mating sires were randomly assigned to females within each age and line, except that half-sib or more closely related matings were avoided. The 60-d breeding season was generally during the months of June and July.
All lines were maintained as one herd except during the breeding season. The herd calved only in the spring. Calves were weaned together each year when the average age was about 200 d (Koch et al., 1974a
,b
; 1994
).
Analysis of Length of Productive Life and Lifetime Production
Analysis of length of productive life requires that only observations from complete cohorts can be used (Schons et al., 1985
; Tanida et al., 1988
). Although a female receives a length of productive life observation when she leaves the herd, that observation cannot be used until all her contemporaries have left the herd. However, some of the drawbacks can be overcome using the concept of opportunity groups (Hudson and Van Vleck, 1981
), which offer considerable flexibility in specification of length of productive life.
Length of productive life (LPL|O) was defined as the number of days a cow survived in the next N years given that she was alive at first calving and had the opportunity to stay in the herd at least N years. Opportunity, in this case, is determined by when the selection experiment was terminated. A cow born 3 yr before termination of data collection and calving at 2 yr of age would have the opportunity to live only 1 yr. Length of productive life was defined for six different opportunity groups: L1, L2, L3, L4, L5, and L6, where, for example, L1 is the length of productive life in days in the next year of a cow that was alive at first calving and had the opportunity to live at least 1 yr longer. Lifetime production (LP) was defined in three ways (all corresponding to opportunity group, L6): 1) total number of calves born by 6 yr after first calving (NB); 2) total number of calves weaned by 6 yr after first calving (NW); and 3) total weaning weight in kilograms of all calves (adjusted to 200-d, and for sex and age of dam), by 6 yr after first calving (CW).
Records of a cow qualified for inclusion only if the cow calved as a 2 yr old. The NB excluded abortions and twins, and therefore is equivalent to total number of single calves born alive or dead. A total of 1,886 cows had records for longevity analysis (WWL, 507; YWL, 506; IXL, 485; and CTL, 388).
Genetic parameters were estimated using a two-trait animal model with data pooled from all lines. For all traits, the model included year the cow was born as a fixed effect and animal direct genetic as a random effect. The two-trait animal model was as follows:
![]() |
where y1 = an n1 x 1 vector of observations for Trait 1, y2 = an n2 x 1 vector of observations for Trait 2, b1 = a vector of fixed effects for Trait 1, b2 = a vector of fixed effects for Trait 2, u1 = a vector of random animal direct genetic values for Trait 1, u2 = vector of random animal direct genetic values for Trait 2, e1 = an n1 x 1 vector of random residual effects for Trait 1, e2 = an n2 x 1 vector of random residual effects for Trait 2, and X, Z = known incidence matrices relating the observations to fixed and random effects, respectively.
The first and second moments of the model are assumed to be:
![]() |
where A = Wrights numerator relationship matrix:
![]() |
,
,
:
![]() |
,
,
, I1 = an identity matrix with order, number of animals with records only for Trait 1, I2 = an identity matrix with order, number of animals with records for both traits, I3 = an identity matrix with order, number of animals with records only for Trait 2.
For pairwise analyses with no missing values (e.g., L6 with NB, NW, and CW [n1 = n2]), I1 and I3 will collapse into I2 (i.e., Re = R0
I2).
Estimates of genetic parameters were obtained using a multiple trait derivative-free algorithm to obtain restricted maximum likelihood estimates with MTDFREML (Boldman et al., 1995
) for combinations of traits with equal number of observations and, an AIREML algorithm (Kachman, 2001
) with MATVEC (Wang, et al., 2002
) for combinations of traits with unequal number of observations. Starting values for the estimates of (co)variance components were from the literature (Tanida et al., 1988
, Morris et al., 1993
, Arthur et al., 1994
). The simplex algorithm was stopped when the variance of the function values (i.e., 2logL with L = likelihood given y) in the simplex was less than 1 x 106. Once the convergence criterion was reached, fresh restarts from those estimates were continued until 2logL differed by less than 1 x 102 between successive restarts. The (co)variance components attained from the last restart were used for estimating breeding values.
Standard errors were obtained directly from MTDFREML program for the two-trait analyses for NB, NW, CW, and L6 because they had the same number of observations using the average information matrix. Due to a restriction in the MTDFREML program, standard errors for the estimates of parameters for multiple-trait analyses cannot be calculated when there are missing observations. Therefore, an AIREML algorithm (Kachman, 2001
) was used with MATVEC (Wang, et al., 2002
) to obtain the information matrix at convergence from the two-trait analyses for the six measures of length of productive life. Standard errors were estimated using the "delta method" and the average information matrix at convergence (Searle et al., 1992
).
Genetic trends were estimated by regressing averages of estimated breeding values by year of birth and line on birth year. The environmental trend was estimated by regressing the solutions for birth year on birth year. Genetic and environmental trends were estimated only for L6, NB, NW, and CW.
| Results and Discussion |
|---|
|
|
|---|
|
|
|
Estimates of genetic correlations were in the range of 0.74 to 1.00 for L1 with L6, for L2 with L3, and L4 with L5, respectively, as would be expected from the part-whole relationships. The estimates of genetic correlations tended to decrease with each added 1-yr increments between the opportunity groups. Thus, for analyses of L1 with L2, L3, L4, L5, and L6, estimates of genetic correlations were 0.89, 0.83, 0.83, 0.76, and 0.74, respectively. The pattern was similar for the other combinations of LPL|O. Similarly, estimates of environmental correlations were high with a range between 0.67 to 0.98, for L1 with L6, and for L5 with L6, respectively. Thus, for analyses of L1 with L2, L3, L4, L5, and L6, estimates of environmental correlations were 0.91, 0.81, 0.71, 0.70, and 0.67, respectively. As with the estimates of genetic correlations, estimates of environmental correlations decreased with each additional 1-yr increment between opportunity groups. The pattern was similar for the other combinations of LPL|O.
Table 4
summarizes results from two-trait analyses for estimates of genetic parameters for length of productive life to 6 y after first calving (L6) and lifetime production defined as NB, NW, or CW by 6 y after first calving. Estimates of heritability for L6, NB, NW, and CW ranged from low to medium. The estimates of heritability averaged 0.15, 0.17, 0.21, and 0.18 for L6, NB, NW, and CW, respectively. These estimates are similar to those reported by Arthur and Makarechian (1992)
, Arthur et al. (1994)
, and Morris et al. (1993)
.
|
Estimates of the genetic and environmental correlations are high, ranging from 0.96 to 1.00 and from 0.56 to 0.99, respectively. These estimates are within the range found in the literature (Tanida et al., 1988
; Arthur and Makarechian, 1992
). For Hereford cows, Tanida et al. (1988)
reported estimates of genetic and phenotypic correlations between length of productive life measured as time between first calving and disposal date of 0.98 and 0.96, respectively. Similarly, estimates of genetic and phenotypic correlations for Angus cows were 1.05 (estimated from Hendersons Method III) and 0.96, respectively. Arthur and Makarechian (1992)
reported estimates of genetic and phenotypic correlations between length of productive life measured as time between first calving and disposal date and number of calves weaned of 0.81 (0.13) and 0.90, respectively. Estimates of genetic and phenotypic correlations between length of productive life and cumulative weaning weight were 0.82 (0.13) and 0.89, respectively (Arthur and Makarechian, 1992
). Estimates of genetic and phenotypic correlations between number of calves weaned and cumulative weaning weight were 1.00 (0.01) and 0.99 (Arthur and Makarechian, 1992
).
Estimates of Genetic and Environmental Changes
Table 5
summarizes estimates of genetic and environmental changes per year for L6, NB, NW, and CW. The environmental trends were significantly negative for all traits. The genetic trends were not different from zero for NB and NW. Although none was statistically significant, estimates of genetic change for CW were 1.24, 2.81, 0.66, and 3.29 kg/yr for WWL, YWL, IXL, and CTL, respectively. The estimates of genetic change for L6 were 1.34, 3.53, 2.05, and 4.63 d/yr for WWL, YWL, IXL, and CTL, respectively, but only the change in CTL was significantly different from zero. No statistically significant differences were found between the selected lines (WWL, YWL, and IXL) and CTL (P < 0.05) for any of the traits. Unfortunately, a lack of research exists about previous estimates of genetic and environmental correlations between measures of length of productive life with different opportunity groups and lifetime production.
|
| Implications |
|---|
|
|
|---|
| Footnotes |
|---|
2 Current address: Instituto de Producción Animal, Facultad de Agronomia, Universidad Central de Venezuela, Apartado Postal 4579, Maracay 2101, Aragua, Venezuela. ![]()
3 Correspondence: A218 Animal Science (phone: 402-472-6010; fax: 402-472-6362; e-mail: lvanvleck{at}unlnotes.unl.edu).
Received for publication September 10, 2003. Accepted for publication April 12, 2004.
| Literature Cited |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. H. Van Melis, J. P. Eler, H. N. Oliveira, G. J. M. Rosa, J. A. V. Silva II, J. B. S. Ferraz, and E. Pereira Study of stayability in Nellore cows using a threshold model J Anim Sci, July 1, 2007; 85(7): 1780 - 1786. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |