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

* Department of Animal Sciences, Colorado State University, Fort Collins 80523-1171; and
and
Institute of Veterinary, Animal and Biomedical Sciences, Massey University, Palmerston North, New Zealand
| Abstract |
|---|
|
|
|---|
Key Words: beef cattle heritability weaning weight environment
| INTRODUCTION |
|---|
|
|
|---|
The partitioning of variation in weaning weights recorded under differing environmental conditions based on individual cow circumstances across years has not been investigated. If estimates of heritability differ for good vs. poor cow environments, and genetic correlations are not unity between environmental conditions, national cattle evaluations would not appropriately account for these differences, and as a result, may inappropriately rank sires.
The purpose of this study was to evaluate the influence of cow environment on expression of genetic differences in weaning weight. The goal was to categorize individual cows as willing to erode their own BW to support lactation for the benefit of their calf or, in contrast, to achieve above average increases in their own body reserves. Heritability and genetic correlations for weaning weights recorded in these alternative cow nutritional environments were estimated to identify whether a single value is appropriate for the differing environmental conditions experienced by cows from year to year.
| MATERIALS AND METHODS |
|---|
|
|
|---|
The general methodology for this study was to classify weaning weights in the RAAA database as occurring in a good or poor cow environment and to subsequently estimate direct, maternal, and permanent environmental variance components from both environments and to compare those with estimates from a model typically used in national cattle evaluation. Good and bad cow environments were determined based on changes in individual cow BW from year to year. The specifics of these procedures are presented in more detail below.
Cow BW Data
Cow BW observations were collected at calf weaning and spanned the years 1981 to 2004. Weight observations of the same cows as yearlings were pooled with these numbers and used to calculate BW change from 1 to 2 yr of age as described later.
Contemporary groups for cow and yearling BW were formed in a similar manner as used in the national cattle evaluation, in accordance with RAAA guidelines. These included percent Red Angus (%RA) of the cow (50
%RA < 87.5; %RA
87.5) and the weaning contemporary group of her calf (at the time at which the cow BW was recorded). Weaning contemporary groups for the calf were formed using work and management group of the calf, calf sex (female, male, and castrated male), breed group of the calf (50
%RA < 87.5; %RA
87.5), weaning date, weaning age group (130 to 159 d, 160 to 250 d, or 251 to 280 d), and Brahman influence (yes/no). Yearling contemporary groups were formed from the weaning contemporary group designation outlined above, in combination with yearling work and management groups. Weaning and yearling management groups are breeder-defined cohorts that have experienced similar management practices and opportunities to perform. Weaning and yearling workgroups identified individuals belonging to the same ranch or herd.
Cow BW observations must have been associated with a valid contemporary group designation (9,912 groups). Some 2,683 observations were removed because their contemporary group exhibited no variation. The final cow BW data set contained 96,064 observations representing 34,223 individual cows.
A 3-generation pedigree consisting of animal ID, sire ID, and dam ID was built based on animals included in the final data set. This pedigree consisted of 69,049 individuals, including 8,071 sires, 39,949 dams, and 6,043 maternal grandsires.
Weaning Weight Data
Age-adjusted weaning weight data from the RAAA national database was used in the analysis. Weaning weight observations were discarded if
A calfs weaning observation, to be included, must have been associated with a corresponding BW observation on its dam. That dam must have had a BW observation taken in a previous year also so that we could determine whether the weaning weight was recorded in a good or a poor environment. These requirements resulted in a data set including 27,534 weaning weight observations in 1,418 contemporary groups.
A 3-generation pedigree was constructed for the calves in this final data file. Foundation dams (dams whose sire and dam were unknown) with only 1 calf and who were missing their own weaning weight observation (n = 3,703) were removed from the pedigree. This resulted in a final pedigree consisting of 47,185 individual animals, including 3,620 sires, 19,765 dams, and 2,818 maternal grandsires.
Weaning Weight Environmental Classification
Yearling and older cow BW were analyzed using a single trait, repeatability animal model to classify each weaning weight observation as occurring in a good or poor cow environment. Fixed effects included in the analysis were cow BW contemporary group and age of cow (1 to 18 yr). Random effects included direct animal and permanent environmental effects. The cow BW model used in the analysis was
![]() |
where X, Z, and W were incidence matrices relating cow BW observations in y to fixed (b), random genetic (u), and permanent environmental effects (p), with e defining a vector of random residual errors. Random effects in the model were assumed to have zero means and variances, represented as follows:
![]() |
where A was Wrights numerator relationship matrix, ID and IN were identity matrices with order equivalent to the number of females with mature BW observations and the total number of BW observations, respectively. Additive genetic (
2u), permanent environmental (
2p), and residual (
2e) variances of cow BW were assumed to be 1,872, 533, and 427 kg2, respectively (Evans, 2001
). This linear model was implemented and analyzed using ASREML (Gilmour et al., 2002
).
Best linear unbiased estimates of cow age solutions were then used to classify the calfs weaning weight as occurring in a good or a poor cow environment. If a cows actual change in BW from one year to the next met or exceeded the expected BW change as determined by the differences between age solutions for cow BW obtained above, her calfs weaning weight, corresponding to the end of this time frame, was classified as having occurred in a good environment. If the observed BW change failed to meet this expected BW change, her calfs weaning weight, corresponding to the end of this time period, was classified as having occurred in a poor environment. For example, according to the best linear unbiased estimates of cow BW per year of age, cow BW was expected to increase 41.6 kg as the cow aged from 2 to 3 yr of age. If the cows actual BW change from 2 to 3 yr of age met or exceeded 41.6 kg, her calf weaned at 3 yr of age was classified as having a weaning weight in a good environment.
Estimation of Weaning Weight Variance Components
Variance components for weaning weight were estimated using 2 models. In model 1, weaning weight was analyzed using a single trait, multiple component animal model represented by
![]() |
where y was a vector of weaning weight observations; X was a known incidence matrix relating the fixed effects in b to the observations in y; and Zw, Zm, and Zp were known incidence matrices relating the random animal effects uw, um, up to the observations in y for direct, maternal, and maternal permanent environmental effects, respectively. The vector e represented random residuals unique to each weaning weight observation in y. The random effects were assumed to have means of zero and variances, shown as
![]() |
where ID and IN were identity matrices with order equal to the number of dams with weaning weights on their calves and the number of weaning weight observations, respectively. Fixed effects included in the model were weaning contemporary group of the calf and a sex x age of dam interaction.
In model 2, weaning weight observations classified as having been observed in good or poor environments were analyzed as correlated traits using a bivariate model. The model used in the analysis is represented in matrix form as
![]() |
where y was a vector of weaning weight observations classified by subscript as good (g) or below average (b), respectively; X was an incidence matrix relating good and poor weaning weight observations to their respective fixed effects in b; Z was a known incidence matrix relating good and poor weaning weight observations to their random effects in u, a vector representing random additive effects for direct, maternal, and maternal permanent environmental effects for good and poor weaning weight, respectively; and e was a vector of good and poor weaning weight random residuals unique to each observation. The random effects were assumed to have means of zero and genetic variances represented by
![]() |
maternal permanent environmental variances represented by
![]() |
and residual variances as
![]() |
In the above equations, subscripts wg, mg, pg, wb, mb, and pb refer to weaning "good" direct, maternal, maternal permanent environment; and weaning "poor" direct, maternal, and maternal permanent environment, respectively. The IDg, IDb, Ig, and Ib were identity matrices with order equal to the number of dams whose offspring have good weaning observations, have poor weaning observations, the total number of good observations and the total number of poor observations, respectively. The matrix E contained all zeros, except for locations corresponding to particular dams with calves in both environments, in which case that element of E was unity. Residual covariances were zero because of the nature of this analysis, in which no individual calf could have both a good and a poor observation recorded.
Variance component estimates were obtained for the above 2 weaning weight models using ASREML (Gilmour et al., 2002
). The numerator relationship matrix was constructed for 47,185 individuals with an average inbreeding level of 0.026.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
The best linear unbiased estimates for cow BW changes are in Table 1
. These differences were used to classify the calfs weaning weight as coming from a good or poor cow environment. Cows were expected to gain 118.9 kg from 1 to 2 yr of age and 41.6 kg from 2 to 3 yr of age based on the fixed effects solutions from the cow BW analysis. This large discrepancy is a result of a varying amount of time elapsing between these measurements. Yearling BW was recorded at approximately 1 yr of age, and 2-yr-old cow BW was taken approximately 1.5 yr later at the time of weaning her first calf. During the first year of these 1.5 yr, the animals are still growing and likely have adequate to abundant fat stores. A trend for gaining BW was evident up to 7 yr of age at which point cows experienced BW loss until 14 yr of age. Few cows were represented above 14 yr of age, and the corresponding age change effects had large standard errors. Weaning weights of cows older than 14 yr were discarded for all subsequent results reported.
|
|
Table 2
shows the number of calves attributed to cows in the data set. About one-third of cows had a single calf in the analysis. About one-sixth of the cows had 2 calves, and these were approximately equally distributed across the 4 sequences poor-poor, poor-good, good-poor, and good-good. There were a total of 13,582 dams that represented 11 different good vs. poor combinations. There were 7,288 dams with only 1 observation, whereas only 2 dams had 11 different observations. The largest percentage (58%) of good observations was represented where 10 observations were recorded. However, there were approximately equal numbers of good and poor observations across all combinations.
|
Weaning Weight Analysis
Weaning weight summary statistics are reported in Table 3
. The mean weaning weight of calves in the poor environment (240.3 kg) was approximately 6 kg lighter than in the good environment (246.9 kg). The poor environment did not contain the lowest weaning weight observation, but the highest weaning weight (375.1 kg) observation in this class was approximately 26 kg less than the maximum weaning weight observation in the good environment (401.0 kg). These simple summary statistics show no practical difference in weaning weight observations between the overall, good, and poor weaning weight observations.
|
|
|
Maternal heritability estimates in the poor environment were marginally higher than in the good environment (0.13 ± 0.02 in the good vs. 0.14 ± 0.02 in the poor). The difference between the maternal estimates of heritability in the good and poor environments was not significant when considering the standard error of the estimates. These results conflict with those reported by Brown et al. (1993)
, where it was found that maternal effects for growth tend to vary with nutritional environment.
Correlations between the direct effects (0.96) in good and poor cow environments and between maternal effects (0.89) in good and poor environments were extremely high. On a wider scale, De Mattos et al. (2000)
found genetic correlations in the Hereford breed for direct and maternal effects across Canada, Brazil, and the United States to be greater than 0.80, indicating sires ranked similarly in all environments. The high genetic correlations found here for direct and maternal effects across good and poor environments and the findings reported by De Mattos et al., (2000)
above indicate weaning weight observations in these 2 environments may not be different traits.
The negative genetic correlations between the direct and maternal components found in model 1 were confirmed in the results from fitting model 2 (0.47 in the good environment and 0.20 in the poor environment). However, weaning weight maternal was more strongly correlated to its direct counterpart in the good environment than in the poor environment. This agrees with our interpretation of the results reported by Snelling et al. (1996)
where the direct maternal genetic correlation decreased from an average of 0.03 to 0.21 with improved environment. Direct comparisons of results are difficult because that study did not attempt to identify cow environment within herd, but appear to support the finding reported herein. There was little difference in estimated genetic correlations between direct genetic effects in one environment and maternal genetic effects in the other environment (0.25 and 0.29).
The maternal permanent environmental effects for these 2 traits are in Table 5
. Maternal permanent environment effects accounted for a similar proportion of variation in weaning weight in good (0.17) and in poor environments (0.15). The correlation between these 2 effects was 0.72.
It might be argued that classification of calves according to the BW gain of their dams imposes a form of indirect selection on weaning weight. Should that be the case, estimates of genetic and perhaps phenotypic variation from fitting model 2 (
2PGood: 539 kg2,
2PBad: 478 kg2) would be eroded in relation to the estimates obtained from model 1 (
2P: 421 kg2). The results do not provide any evidence of reduced direct, maternal genetic, or phenotypic variation.
It was hypothesized that anabolic and catabolic systems are under different genetic control, and therefore the inheritance of maternal effects may differ in cows that have above- vs. below-average BW gains during lactation. However, unexpected findings were the near-perfect genetic correlation and high permanent environmental correlation between maternal effects in the 2 environments. These results may argue that the classification of cows was poorly achieved using the single mature BW observed at weaning each year. Further studies with cows classified using a better measure of cow energy status (e.g., monthly BW, condition scores, or ultrasonic measures) would be desirable.
Additional analysis prompted by the comments of an anonymous referee were undertaken using model 2 but with stricter classification of cows into the poor environment. A cows yearly BW change had to be in the bottom 10% of BW gain for their age group for their calfs corresponding weaning weight to be classified into the poor environment. All other BW were classified as having occurred in a good environment. Heritabilities and genetic correlations for this good/poor division are in Table 6
. No differences were seen between heritability estimates in both good and poor environments for both weaning weight direct (0.24 ± 0.03 in the good environment vs. 0.25 ± 0.09 in the poor environment) and weaning weight maternal (0.13 ± 0.02 in the good environment vs. 0.13 ± 0.08 in the poor environment).
|
Alternative methods for classifying individual, yearly, cow environment were considered, but all introduced confounding of genetic potential and environmental classification or altered the question addressed by classifying entire contemporary groups as performing in a single environment. For instance, grouping calves based on cow milk EPD would reduce genetic variance in milk production.
In an effort to assure that there was no confounding of growth and milk production potential in this study, an additional analysis was performed. Cows with a single calf weaned were eliminated from the data because these seemed the most likely candidates for confounding because cows with multiple observations readily reclassified as previously discussed. All heritability and genetic correlation estimates were within 1 standard error as reported in Table 5
with the exception of maternal heritability in the good environment which increased to 0.30 ± 0.03 (as did maternal heritability in the poor environment) and C2, which was 0.14 ± 0.02 and 0.11 ± 0.02 in the good and poor environments, respectively.
| IMPLICATIONS |
|---|
|
|
|---|
| Footnotes |
|---|
2 Corresponding author: mark.enns{at}colostate.edu
Received for publication February 16, 2006. Accepted for publication October 26, 2006.
| LITERATURE CITED |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
K. L. Shirley, D. W. Beckman, and D. J. Garrick Inheritance of pulmonary arterial pressure in Angus cattle and its correlation with growth J Anim Sci, April 1, 2008; 86(4): 815 - 819. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |