|
|
||||||||
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 |
|---|
|
|
|---|
The ability of cows to have a long productive life is important for commercial beef producers because a longer productive life means lower costs for rearing replacements, fewer young cows, and thus more and heavier calves available for sale (Rendel and Robertson, 1950
). For purebred cattle breeders, increased lifetime production would permit greater selection intensity for other important traits. The weight of weaned calves per initial replacement female accumulated over a lifetime is a measure of the cows contribution to the genotype of its calves for growth, its fertility (pregnancy and calving rate), its maternal ability (weaning rate), milking capacity (maternal weaning weight), and its survival (ability of the cow to delay culling or death), and has been proposed as a comprehensive measure of lifetime production (Tanida et al., 1988
; Nuñez-Dominguez et al., 1991
).
A previous study (Martinez et al., 2004) on six measures of length of productive life and three measures of lifetime production by 6 yr after first calving, reported high genetic and environmental correlations among the three measures of lifetime production.
The objectives of this study were 1) to estimate genetic parameters for three measures of lifetime production conditional on the opportunity to produce 2 through 7 yr of age, and 2) to estimate genetic and environmental trends for these measures of lifetime production for a control and three selected lines of Hereford cows.
| Materials and Methods |
|---|
|
|
|---|
Data were from the Nebraska Agric. Exp. Stn. Project 40-002 entitled "Effect of selection for weaning weight, yearling weight, and muscling in beef cattle," in cooperation with the Roman L. Hruska U.S. Meat Animal Research Center (USMARC) (Koch et al., 1974a
,b
; 1994
). Data used were from animals born from 1964 through 1980.
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 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 (CTL). In 1971, the CTL was established from 20 representative sons and heifer calves from matings with the 225 cows (Koch et al., 1974a
,b
; 1994
).
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
; 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. From 1964 through 1970, the top 25 heifers from each line were selected based on their line data. From 1971 through 1973, the top 35 heifers were selected. From 1974 through 1984, all pregnant heifers were kept in the herd. No heifers that were pregnant were eliminated from the analysis. 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, which was rare (e.g., cancer eye, chronic bloat, bad mouth); 3) failure to raise a live calf for two consecutive years; and 4) older cows were removed if cows needed to be culled to maintain 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.
These criteria were the basis of the management protocol to maintain herd size, although exceptions were occasionally made if more cows were needed (e.g., nonpregnant cows would be carried over to maintain herd size). The number of cows in these categories was not counted.
All lines were maintained as one herd except during the 60-d breeding season. The herd calved only in the spring, mostly in March and April. Cows were palpated for pregnancy once at weaning. Heifers were palpated at about 18 mo of age. Calves were weaned together each year when the average age was about 200 d (Koch et al., 1974a
,b
; 1994
).
Selection Objectives in the Original Project
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 also 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.
Analysis of Measures of Lifetime Production
Records qualified for analysis only if the cow were available for breeding as a yearling. Thus, lifetime production traits were defined only for heifers that entered the breeding program as yearlings. For each of the three measures of lifetime production, six opportunity groups were defined, based on whether the cow was born early enough in the experiment to have the opportunity to live to calve at 2, 3, 4, 5, 6, and 7 y of age.
The three general measures of lifetime production (LP) were defined as 1) total number of calves born (NB) by 2 (NB2), 3 (NB3), 4 (NB4), 5 (NB5), 6 (NB6), and 7 yr of age (NB7); 2) total number of calves weaned (NW) by 2 (NW2), 3 (NW3), 4 (NW4), 5 (NW5), 6 (NW6), and 7 yr of age (NW7); and 3) total weaning weight in kg of all calves (CW; adjusted to 200 d of age and for sex and age of dam) by 2 (CW2), 3 (CW3), 4 (CW4), 5 (CW5), 6 (CW6), and 7 yr of age (CW7).
These definitions imply part-whole correlations within a type of measure. Cows that were not pregnant in any year were given values of zero for NB, NW, and CW, as were pregnant cows that were kept in the herd after weaning but were culled before calving.
Number of calves born did not include abortions at any stage of pregnancy and did not include twins, and therefore is equivalent to total single calves born alive or dead. Records of a cow that had twins were deleted from all analyses.
Genetic parameters were estimated using two-trait animal models within type of measure (e.g., NB2 with NB3, NB4, NB5, NB6, and NB7) and similarly for NW and CW. In addition, two-trait animal models were used for analyses of NB2 with all opportunity groups for NW and CW and of NW2 with all opportunity groups for CW.
For all traits, the model included birth year of the cow as a fixed effect and random direct genetic 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 = vector of fixed effects for Trait 1; b2 = vector of fixed effects for Trait 2; u1 = vector of random animal direct genetic effects for Trait 1; u2 = vector of random animal direct genetic effects 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, with u1, u2, Z1, and Z2 augmented for relatives without records.
The first and second moments of the model are assumed to be:
![]() |
where A = Wrights numerator relationship matrix and
![]() |
= additive genetic variance for Trait 1;
= additive genetic variance for Trait 2;
u1u2 = additive genetic covariance between Traits 1 and 2:
![]() |
= residual variance for Trait 1;
= residual variance for Trait 2;
e1e2 = residual covariance between Traits 1 and 2; 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; and I3 = an identity matrix with order number of animals with records only for Trait 2.
For pairwise analyses with no missing values (e.g., NB2 with NW2, NB2 with CW2, and NW2 with CW2 [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
). Starting values for the estimates of (co)variance components were from the literature (Morris et al., 1987
, 1993
; Tanida et al., 1988
). 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 x106. Once that convergence criterion was reached, fresh restarts from those estimates were continued until 2logL differed 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 for two-trait analyses with the same number of observations were obtained directly from MTDFREML using the average information matrix. The MTDFREML program does not calculate standard errors for estimates of genetic parameters for multiple-trait analyses when there are missing observations. Therefore, an AIREML algorithm (Kachman, 2001
) was used with MATVEC (Wang et al., 2002
) to obtain the average information matrix at convergence. Standard errors were estimated using the "delta method" and the average information matrix at convergence (Searle et al., 1992
).
Estimates of genetic change were obtained by regressing the average of estimated breeding values by year of birth for each line on birth year. Estimates of environmental change were estimated by regressing the solutions for year of birth on birth year. Estimates of genetic and environmental change were estimated only for NB2, NB6, NW2, NW6, CW2, and CW6.
| Results and Discussion |
|---|
|
|
|---|
|
Table 2
summarizes estimates of heritability and genetic and environmental correlations for NB2 to NB7. Estimates of heritability for NB were low, ranging from 0.08 for NB2 to 0.16 for NB7. Estimates of heritability tended to increase with added length of the opportunity group with averages of 0.08, 0.10, 0.11, 0.12, 0.14, and 0.16 for NB2, NB3, NB4, NB5, NB6, and NB7, respectively. The increase in estimates of heritability may be due to increased opportunity for genetic differences to be expressed. The estimates of heritability are similar to previous reports. Morris et al. (1987)
reported an estimate of heritability for number born (single calves alive or dead) of 0.05 (0.05) for crossbred cows. Morris et al. (1993)
reported an estimate of heritability of 0.06 (0.22) for Hereford cows, 0.11 (0.09) for crossbred cows, and 0.11 (0.08) for combined data of Hereford and crossbred cows.
|
Estimates of environmental correlations were slightly larger than estimates of genetic correlations, but as with estimates of genetic correlations the estimate decreased with each added increment of one year for the opportunity group. Estimates of environmental correlations for NB2 with NB3, NB4, NB5, NB6, and NB7, were 0.87, 0.78, 0.73, 0.69, and 0.66, respectively.
Table 3
summarizes estimates of heritability and genetic and environmental correlations for NW2 to NW7. Estimates of heritability for NW were low, ranging from 0.05 for NW2 to 0.16 for NW7. Estimates of heritability tended to increase with each increment in length for the opportunity group, with average estimates of 0.05, 0.09, 0.11, 0.13, 0.15, and 0.16 for NW2, NW3, NW4, NW5, NW6, and NW7, respectively. Estimates of heritability were similar to those found in the literature. Morris et al. (1987)
reported an estimate of heritability for number weaned of 0.02 (0.04) for crossbred cows. Tanida et al. (1988)
reported an estimate of heritability of 0.16 (0.08) from daughter-dam regression, and 0.22 (0.08) from paternal half-sib analysis for Hereford cows, and 0.03 (0.14) from paternal half-sib analysis for Angus cows. Arthur and Makarechian (1992)
indicated an estimate of heritability of 0.28 (0.14), and Arthur et al. (1994)
reported an estimate of 0.24 (0.02) for Hereford cows. Morris et al. (1993)
reported an estimate of heritability of 0.13 (0.23) for Hereford cows, 0.03 (0.08) for crossbred cows, and 0.15 (0.08) for the combined analysis of Hereford and crossbred cows.
|
Table 4
summarizes estimates of heritability and genetic and environmental correlations for CW2 to CW7. Estimates of heritability for CW were low, ranging from 0.06 for CW2 to 0.16 for CW7. Estimates of heritability tended to increase with increased length of opportunity group with average estimates of 0.06, 0.10, 0.11, 0.12, 0.14, and 0.16 for CW2, CW3, CW4, CW5, CW6, and CW7, respectively. Estimates of heritability are smaller than those reported by Arthur and Makarechian (1992)
, who reported an estimate of 0.27 (0.14) with a sire model and Arthur et al. (1994)
, who reported an estimate of 0.30 (0.01) for Hereford cows with an animal model.
|
Table 5
summarizes estimates of heritability and genetic and environmental correlations for NB2 with NW2 through NW7. In general, estimates of heritability were similar to those from two-trait analyses within type of measurement. Estimates of genetic and environmental correlations were moderate to high. Estimates of genetic correlations ranged from 0.47 between NB2 and NW4 to 0.71 between NB2 and NW2. The general tendency was for estimates of correlations to decrease with increase in length of time between opportunity groups. Estimates of environmental correlations ranged from 0.61 for NB2 with NW7 to 0.68 for NB2 with both NW2 and NW3. As with estimates of genetic correlations, estimates of environmental correlations tended to decline with added length of time between opportunity groups for NW.
|
|
|
Estimates of Genetic and Environmental Changes
Table 8
summarizes estimates of genetic and environmental change per year for NB2, NB6, NW2, NW6, CW2, and CW6. The estimates of genetic trend were near to zero for NB2, NB6, NW2, and NW6. Although estimates of genetic trend for CW2 (P < 0.05) were less than 1 kg/year, all were different from zero for all lines. However, estimates of 2 to 5 kg/year did not differ from zero for CW6 (P > 0.05). Genetic trends were not significantly different between selected lines and the control line. The environmental trend was slightly positive for traits measured early in life such as NB2, NW2, and CW2, but the trends were negative for NB6, NW6, and CW6, which are measured later in life.
|
| 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:
![]() |
K. M. Cammack, M. G. Thomas, and R. M. Enns Review: Reproductive Traits and Their Heritabilities in Beef Cattle Professional Animal Scientist, October 1, 2009; 25(5): 517 - 528. [Abstract] [PDF] |
||||
![]() |
R. A. Cushman, M. F. Allan, L. A. Kuehn, W. M. Snelling, A. S. Cupp, and H. C. Freetly Evaluation of antral follicle count and ovarian morphology in crossbred beef cows: Investigation of influence of stage of the estrous cycle, age, and birth weight J Anim Sci, June 1, 2009; 87(6): 1971 - 1980. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |