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

* Graduate School of Agricultural Science, Faculty of Agriculture, Tohoku University, Sendai, Miyagi Prefecture, 981-8555 Japan; and
and
Miyagi Prefecture Animal Industry Experiment Station, Iwadeyama-cho, Miyagi Prefecture, 989-6445 Japan
| Abstract |
|---|
|
|
|---|
Key Words: Genetic Effects Growth Japanese Black Cattle Random Regression
| Introduction |
|---|
|
|
|---|
Japanese Black cattle are famous for the high marbling score of their beef. Genetic studies for carcass traits are numerous for this breed (Mukai, 1994
; Kitamura et al., 1999
; Uchida et al., 2001
). However, there are comparatively few studies dealing with parameter estimates for growth traits (Shimada et al., 1995
; Mukai et al., 2000
). Estimates of genetic parameters for weights in a sequence of small age intervals are not available. Of interest are genetic parameters that describe the change of such traits in time. From an animal breeding point of view, these genetic parameters reflect to what extent and how genetic changes in growth patterns over time can be achieved by selection.
The objectives of this study were 1) to estimate coefficients of covariance functions for direct and maternal genetic and permanent environmental effects, and 2) to use random regression to estimate genetic, environmental, and phenotypic parameters of the aforementioned effects for weights taken from birth to 356 d of age in a herd of Japanese Black calves.
| Materials and Methods |
|---|
|
|
|---|
Japanese Black is the most common cattle breed in Japan. Japanese Black cattle were established in 1944 by crossing Japanese native cattle with various breeds such as Shorthorn, Simmental, and Brown Swiss. The breeds used and the structure of the crossbreeding program varied among prefectural populations. Even after the establishment of the breed, most breeding programs continued to be carried out separately within each prefecture, although there has been some genetic flow among the various subpopulations (Nomura and Sasaki, 1988
). In fact, the breed is genetically divided into several subpopulations, each with a different genetic background. The current breed is characterized by high meat quality, with prominent i.m. fat deposition (i.e., marbling).
The initial data set comprised 17,368 weight records on 3,272 Japanese Black calves. These data were collected from 1997 to 2003 in a farm located in Shiroshi city, Myiagi prefecture, Japan. This farm has an operational genetic improvement program that considers growth, marbling score, fertility, and maternal ability traits. Calves are weighed at birth and subsequently at monthly intervals. They are weaned at approximately 6 mo of age. Calves are given free access to concentrate mixed with 10% chopped rice straw by weight. The concentrate comprised 20% ground barley, 35% ground yellow corn, 20% wheat bran, 17% defatted rice bran, 6% soybean meal, 1% NaCl, and 1% calcium carbonate (as-fed basis). In addition, roughage (in the form of green forage, silage, or hay) and water are supplied ad libitum. A unique number assigned by the Wagyu Registry Association, which maintains pedigree records, identifies animals. Information available is animal name and number, parental identifications, and date and place of birth. Traits included in this study were birth weight and nine consecutive monthly weights up to 356 d of age. There were few records after 356 d, and these were eliminated to avoid computational problems due to small numbers of records per subclass. Records with implausible dates or weights were also eliminated. The edited data set included 11,815 weights from 2,394 calves, offspring of 45 sires, and 631 dams. Calves were born year round. All animals had records for birth weight. The number of age subclasses including birth weight was 269. Contemporary groups (CG) were defined as year and season of birth-sex of calf. The number of animals in the relationship matrix was 3,082. All calves and dams had their parents identified. The percentage of dams with at least two progeny was 90. There was some overlapping of ages for different measurements. Figure 1
shows the distribution of weight recordings over age in days. Table 1
contains number of records, unadjusted means, standard deviations, and CV for the 10 consecutive weights and their corresponding age subclasses. Figure 2
depicts BW means throughout the complete test period.
|
|
|
The model considered for fitting repeated weights of calves included parity of dam as a fixed effect and fixed regressions of weights on fourth-order orthogonal Legendre polynomials of actual day of age for CG effects. Random regressions on fourth-order Legendre polynomials were included for animal additive genetic and animal permanent environmental effects, whereas third-order Legendre polynomials were used for maternal genetic and maternal permanent environmental effects. The model was as follows:
![]() |
where yij is the jth weight of the ith calf, bn(j) is the fixed linear regression coefficient of the nth parity for the jth measurement, Pn is the nth parity of dam effect, tij is the ijth actual day of age at recording for which Legendre polynomials are defined, standardized to the range 1 to +1, and
m(tij) is the corresponding mth Legendre polynomials, ßm is the mth fixed regression coefficient for CG,
m,
m,
m, and
m are the mth-order random regression coefficients for animal additive genetic (G), maternal genetic (M), animal permanent environmental (P), and maternal permanent environmental (C) effects, respectively, NG and NP denote the orders of fit (G = 4 and P = 3), and
ij is the ijth temporary random environmental effect (residual) that was assumed to be independently distributed with heterogeneous variances. Heterogeneity of variances was modeled as a step function of the covariable, defining e = 10 classes, including a separate class for birth weight, (0, 1 to 90, 91 to 120, 121 to 150, 151 to 180, 181 to 210, 211 to 240, 241 to 270, 271 to 300, and 301 to 356 d). The ranges specified assumed that a separate error variance needed to be fitted for records taken at birth (0 d), and subsequently every month (30 d), except for the second and the last classes, which span 90 and 56 d, respectively.
The model written in matrix notation was:
![]() |
where y is the vector of records; ß is the vector of fixed effects, including parity of dam and CG;
,
,
, and
are vectors of animal additive genetic, maternal genetic, animal permanent environmental, and maternal permanent environmental effects, respectively; X is an incidence matrix for fixed effects; Z1, Z2, Z3, and Z4 are incidence matrices for G, M, P, and C effects, respectively; and e is the residual vector.
The set of orthogonal polynomials for age can be arranged in the following matrix (
):
![]() |
where M is a matrix containing the polynomials of the corresponding standardized ages and
is a matrix containing the coefficients of the Legendre polynomials of a given order. The elements of the additive genetic covariance matrix G can be expressed as follows:
![]() |
where K is a matrix of coefficients that define the covariance function. Covariance matrices for the other random effects (M, P, and C) can be expressed in similar fashion.
The first and second moments of the model are:
![]() |
where A and I are the numerator relationship matrix and the identity matrix, respectively; K
, K
, K
, and K
are the matrices of coefficients of the covariance functions for animal additive genetic, maternal genetic, animal permanent environmental, and maternal permanent environmental effects, with dimensions equal to the order of the covariances matrices in the analysis; and R is the variance of the residual effects, which is a diagonal matrix allowing for different variances for each age t.
The data allowed permanent environmental and maternal genetic effects to be separated, as there was a large number of dams with more than one progeny with records. The random regression model (RRM) was fitted on Legendre polynomials of age. A fourth-order polynomial was fitted for G and P effects, whereas a third-order polynomial was fitted to M and C effects. A preliminary analysis indicated that higher order polynomials were unnecessary. Covariances between random regression coefficients pertaining to different random factors were assumed to be zero. Matrices K
, K
, K
, and K
corresponding to G, M, P, and C effects were estimated as covariance matrices for the corresponding random regression coefficients. Elements of K
, K
, K
, or K
are the coefficients of the covariance functions defining covariances between any two ages in the data for the corresponding random effect (Meyer and Hill, 1997
). Thus, the estimated covariance for the rth random regression (r represents
,
,
, or
) for animal i at ages tij and tij'was as follows:
![]() |
where
rmn is the mnth element of Kr (Kirkpatrick, et al. 1990
). According to Kirkpatrick and Heckman (1989)
, fewer dimensions of Kr often produce a good approximation of
r, and eigenvalues and eigenvectors of Kr can be used as approximations of eigenvalues and eigenfunctions of
r.
Statistical Analyses
Covariances between random regression coefficients of K
, K
, K
, or K
were estimated by Bayesian procedures. Gibbs sampler was used to obtain marginal posterior distributions for all parameters of the model.
Estimates were obtained by RRGIBBS program (Meyer, 2002
). This program generates a Markov chain for the parameters in the model, sampling from conditional posterior distributions as described by Jamrozik and Schaeffer (1997)
and Rekaya et al. (1999)
. It allows the Bayesian estimation of (co)variance components for random regression models using Gibbs sampling. Convergence of the Gibbs chain was monitored by visual inspection of plots of samples for selected components. Flat priors were assumed for fixed location parameters. Prior inverted Wishart distributions were assumed for G, M, P, and C effects, and prior independent scaled inverted
2 distributions were assumed for residual variances. The hyperparameters of the inverted Wishart distributions (
and S) were 5 and 0 for G and P, and 4 and 0 for M and C effects, whereas those of the scaled inverted
2 distributions (
and s2) were equal to 2 and 0, making the distributions flat (i.e., uninformative). After a burn-in period of 20,000 samples, 270,000 samples were used to estimate means for all (co)variances of the RRM. From the estimates of the random regression parameters, variance and covariance components and genetic and phenotypic parameters were calculated for the ages of interest (i.e., at birth, 90, 120, 150, 180, 210, 240, 270, 300, 330, and 356 d). These ages were chosen because they almost cover the age trajectory. Heritability estimate at age i was calculated as follows:
![]() |
where
2a is the additive genetic variance and
2p is the phenotypic variance, which is the sum of additive genetic, maternal genetic, direct permanent environmental, maternal permanent environmental, and error variances at age i. Ratios of maternal genetic, direct permanent environmental, and maternal permanent environmental variances to phenotypic variances were expressed in similar fashion. Posterior standard deviations of the estimated genetic and phenotypic parameters were estimated by the Delta method as described by Lynch and Walsh (1997)
.
| Results and Discussion |
|---|
|
|
|---|
Table 2
presents the estimated coefficient matrices (K
, K
, K
, and K
) of the covariance functions, their respective eigenvalues, their proportions out of the total variation and the correlations between the coefficients for G, P, M, and C effects. There were strong positive correlations between the intercept and the linear coefficients of the different effects. The correlations between the intercept and the linear coefficients of G and P were higher than those of M and C. The correlations between the intercept and the quadratic regression coefficients of all effects were negative, ranging from low (P and C) to moderate (G and M). Albuquerque and Meyer (2001a)
described similar results. They reported moderate to high positive correlations between intercept and linear coefficients and low negative correlations between intercept and quadratic regression coefficients, when fitting random regression for the growth curve from birth to 630 d of age in Nelore cattle.
|
Animal Parameters
Estimates of heritability (h2), ratio of animal permanent environmental variance to phenotypic variance (p2), and repeatability (t = h2 + p2) for the whole period are shown in Figure 3
. Estimates of h2, genetic (rG) and phenotypic (rPH) correlations for weights at 0, 90, 120, 150, 180, 210, 240, 270, 300, 330, and 356 d of age are presented in Table 3
. Estimates of h2 increased from 0.38 at birth up to 0.66 at 96 d, fluctuated slightly around this value up to 120 d, then decreased to 0.38 at 300 d and increased again up to 0.47 at 356 d. The estimates of h2 were consistently higher due to higher estimates of
2A (data not shown). Residual variances for weights at all calf ages are shown at the bottom of Table 3
. Residual variances were heterogeneous across ages. Their contributions to the total variation for weight ranged from 4 to 24%.
|
|
|
Information on parameter estimates for growth of Japanese black cattle is scarce. All estimates found in the literature were calculated by methods other than random regression. Estimates were variable, reflecting the different genetic backgrounds of Japanese Black cattle subpopulations. Sasaki et al. (1982)
and Kitamura et al. (1999)
reported lower h2 values for BW at 210 and 252 d of age, accounting for 0.20 and 0.22, respectively, compared with the estimates of 0.32 and 0.38 obtained here at comparable ages (210 and 240 d, respectively). Mukai et al. (2000)
estimated h2 values for data collected from three prefectures on BW of Japanese Black calves at birth and 300 d of age. Their estimates ranged between 0.19 and 0.39 and between 0.25 and 0.30 for both traits, respectively. These estimates are in partial agreement with the h2 value of 0.38 reported here at the same ages (Table 3
). On the other hand, Ujiie et al. (1990)
reported markedly higher h2 estimate (0.82) for final BW at 240 d of Japanese Black calves compared with the corresponding estimate of 0.48 obtained here. Treating BW at different ages as different traits, Shimada et al. (1995)
reported h2 estimates accounting for 0.61, 0.29, 0.28, 0.30, and 0.33 for BW of Japanese Black calves at birth, 90, 120, 150, and 180 d of age. Although their h2 estimate of birth weight is higher than the corresponding value (0.38) reported here, the other h2 estimates were lower. The discrepancy observed between the estimates reported in the literature, concerning growth of the Japanese Black cattle, and those of the present study could be attributed to the different methods of analyses, different sizes of data, and/or the genetic diversity of the breed among subpopulations.
In general, genetic, phenotypic, and environmental correlations decreased as the interval between ages increased. Genetic correlations between birth weight and other weights were lower than the corresponding values observed between weights at adjacent ages, suggesting that birth weight is not under the same genetic control as weights at older ages. Albuquerque and Meyer (2001a)
reported similar findings; correlations of birth weight with weights taken at other ages were lower than those between weights measured at adjacent ages. Our results are also consistent with the results obtained by Meyer (2001)
, Huisman et al. (2002)
, and Schenkel et al. (2002)
, who found that correlations between adjacent ages were high and decreased as the interval between ages increased. Genetic correlations between weights at 180 and 210 d, between weights at 210 and 240 d, and between weights at 240 and 270 d approached unity. The genetic correlations were all positive, suggesting that selection for weight at an early age would have a positive effect on weights at later ages. Working on Japanese Black cattle, Mukai et al. (2000)
reported genetic correlations between weights at birth and at 300 d of age ranging between 0.40 and 0.87. Our rG estimate for this trait (0.78) falls within this range.
Phenotypic correlations were lower than those due to additive genetic, with a minimum of 0.47 between weights at birth and at 90 d, and a maximum of 0.93 between weights at 240 d and 270 d of age. Their pattern, however, was similar (i.e., correlations between adjacent ages were high and decreased as the intervals increased).
Animal permanent environmental correlations between weights at birth and other weights were lower than those among weights at other ages. They were slightly lower than those due to additive genetic, but they were higher than the phenotypic ones, with a minimum of 0.40 between weights at birth and 356 d, and a maximum approached unity between weights at 90 d and 120 d, between 120 d and 150 d, between 150 d and 180 d and between 180 d and 210 d of age.
Maternal Parameters
Maternal genetic effects of growth traits are known to be important in beef cattle. Partitioning maternal effects into maternal genetic and maternal permanent environmental effects presents some difficulties, according to Willham (1980), Gerstmayr (1992)
, and Meyer (1992)
. There is evidence that most of the maternal variation would be accounted for by including only one of the maternal effects (maternal genetic ro maternal permanent environmental), Meyer (1992)
.
Maternal heritability (m2) and maternal permanent environmental estimates (c2) for the whole period are presented in Figure 4
. Estimates of m2 and maternal genetic (rM) correlations for weights at 0, 90, 120, 150, 180, 210, 240, 270, 300, 330, and 356 d of age are presented in Table 5
. Estimates of c2 and maternal permanent environmental (rM) correlations for the aforementioned ages are shown in Table 6
. Maternal heritabilities increased from 0.04 at birth up to 0.09 at 120 d, and were around 0.06 thereafter. The c2 estimates were approximately equal over the age ranges, accounting for 0.03 of variance. The results of Meyer et al. (1993)
indicated that maternal effect might affect postweaning growth as a carryover effect from weaning weight. If this is the case, maternal effect might have some influence on postweaning growth. Meyer (1992)
presented an extensive list of literature estimates of maternal heritability and showed that estimates for postweaning gain and final weight are generally low, ranging between 0.01 and 0.09. The m2 estimates for postweaning growth, reported here, are in agreement with the values reviewed by Meyer (1992)
. Waldron et al. (1993)
estimated m2 and c2 for weight traits in one Hereford and two Angus experimental herds in Australia, using a univariate animal model. The m2 values reported in their study were 0.07, 0.13, 0.04, and 0.04 for weights at birth, weaning, yearling, and 18 mo of age, respectively, and the respective c2 estimates were 0.05, 0.15, 0.03, and 0.01. The corresponding m2 and c2 estimates reported here are in partial agreement with their findings. Furthermore, they concluded that m2 progressively decreased in importance with time after weaning. Our results indicated that maternal effect decreased from 0.09 at 120 d of age to 0.06 at 180 d of age (weaning age) and remained constant around this value up to the end of the age trajectory. Furthermore, Albuquerque and Meyer (2001b)
reported results similar to those obtained here. They indicated that m2 increased from 0.01 at birth to 0.07 and 0.08 at 180 d and 210 d, after which estimates decreased slowly with age. Conversely, Meyer (2001)
reported c2 estimates higher than those of m2, noting that the later estimates varied slightly over the range of ages. In contrast, our m2 estimates were slightly larger than those of c2.
|
|
|
Similar patterns to the estimates of rG and rP were also observed for the maternal genetic (rM) and maternal permanent environment (rC) correlations (i.e., estimates decreased with increasing time between measurements). Estimates of rM were lower than those of rC. Although rM estimates were lower than those due to additive genetic, rC estimates were approximately equal to those of rP. The estimates of rM and rC of birth weight with the other weights were lower than those of weights taken at older ages. The lowest estimates of rM and rC were between birth weight and weight at 356 d of age, accounting for 0.19 and 0.35, respectively. Estimates of rM between weights at 210 d and 240 d and between 240 d and 270 d were the highest (0.98), whereas rC estimates between weights at 90 d and 120 d and between 120 d and 150 d approached unity. There were positive and strong estimates of rM and rC between weights at adjacent ages, starting from 3 mo, indicating that maternal effects were essentially identical.
General Discussion
To improve growth performance of Japanese Black cattle, it is necessary to appropriately model the genetic and environmental factors contributing to the performance over the test period. By using random regression, it is possible to fit a more complete model than with traditional methods. Random regression allows the trajectory of genetic and environmental parameters over the test period to be investigated. Random regression models accommodate time dependent changes in the genetic and environmental factors; however, the key issues in application of RRM to growth traits are the number of times individuals need to be measured, at what times in their lives, and what will be the upper age range. Here, covariance matrices for G, P, M, and C were estimated by Bayesian procedures via Gibbs sampling. As indicated by the eigenvalues, the intercept and the linear coefficients accounted for most of the variance in weights of calves. The quadratic and cubic terms seem to still be relevant for additive genetic and direct permanent environmental effects, as well as the quadratic term for the maternal genetic and maternal permanent environmental effects.
Information obtained here is useful for determining the optimal length of test period, to assess the role of environmental effects along the test period, and to reflect to what extent and how genetic changes in performance over time can be achieved by selection. The random regression model allowed the comparison of different testing ages and the results showed that h2 and p2 of weights taken at the different ages were moderate to high, whereas those due to maternal genetic and maternal permanent environment were low, but significant. In beef cattle, maternal effects are important for growth traits until weaning although significant effects remain thereafter.
Direct and maternal heritabilities were highest at 120 d. Direct heritability estimates of weights at all ages were moderate to high, indicating that selection on weight at any age is feasible. Genetic correlations across ages were positive. Consequently, selection for increased weight at any age would tend to increase weights throughout the entire growth period. Another important consequence is that preweaning weights can be used to select for weights at later ages.
Genetic evaluation of growth of Japanese Black cattle often considers three different traits, namely birth weight, starting weight, and final weight at fattening period. Obviously, the last two weights are taken at arbitrary ages and estimating their genetic parameters will have little effect on growth improvement in the breed. However, random regression models using records taken at different ages would be more appropriate because they would allow the use of more weights and, hence, more accuracy will be achieved.
Differences between estimated genetic parameters reported in the literature on Japanese Black cattle and those estimated here could be attributed to different methods of estimation used and to the use of data sets from different subpopulations of Japanese Black cattle with distinct genetic backgrounds. In particular, inferences from the results of this study cannot be extrapolated to the complete population of Japanese Black cattle because the data used represent a single herd from only one prefecture.
| Implications |
|---|
|
|
|---|
1 Correspondencephone: +81 22 717 8697; e-mail: maaziz{at}bios.tohoku.ac.jp or fassad9{at}yahoo.com.
Received for publication August 18, 2004. Accepted for publication December 8, 2004.
| Literature Cited |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |