|
|
||||||||

* National Agricultural Research Center for Hokkaido Region, Hitsujigaoka 1, Toyohiraku, Sapporo, Japan 0628555
Dairy and Swine Research and Development Centre, Agriculture and Agri-Food Canada, University of Guelph, Guelph, Ontario, Canada N1G 2W1
1 Corresponding author: tkenji{at}naro.affrc.go.jp
| ABSTRACT |
|---|
|
|
|---|
Key Words: lactation curve lactation milk persistency eigenvector index
| INTRODUCTION |
|---|
|
|
|---|
Olori et al. (1999) reported that the first (leading) eigenvector of the additive genetic covariance function was related to increased mean milk yields at each stage of lactation; and the second eigenvector was related to the change in the shape of the lactation curve. Macciotta et al. (2004) used the eigenvectors of the correlation matrix of test-day records to derive 2 index scores, which were related to the rate of ascent to the lactation peak and the rate of decline after the peak (i.e., persistency). That study pointed out the possible relationship of the eigenvectors to lactation milk and persistency on a population basis, but did not show how to use them for genetic selection. The objectives of this paper are to 1) demonstrate the construction of various eigenvector indexes for improving lactation milk and persistency, 2) characterize the genetic response patterns across lactation of individual eigenvectors in response to selection, 3) study the selection effectiveness of eigenvector indexes as compared with conventional selection based on lactation EBV, and 4) determine how many eigenvectors are necessary to explain most of the variation in the breeding goal.
| MATERIALS AND METHODS |
|---|
|
|
|---|
= [
0
1 · · ·
k1]' be a (k x 1) vector of the additive genetic random regression (RR) coefficients due to the fitting of a Legendre polynomial of degree (k 1). The variance of vector
is a (k x k) additive genetic RR covariance matrix (K) with k pairs of eigen values (
i) and normalized (orthogonal) eigenvectors (ei, i = 1, 2, ..., k). Let E be a (k x k) matrix containing these orthogonal eigenvectors as columns. The genetic covariance matrix (G) of daily yields from DIM 5 to 305 is G =
K
' where
is a (301 x k) matrix of Legendre polynomial coefficients (i.e., covariates) evaluated from DIM 5 through 305. The index (IK) constructed based on the eigenvectors of K is defined as
![]() |
where b is a (k x 1) vector of index coefficients. By this definition, an "index trait" (
'ei) is a linear combination of additive genetic RR coefficients weighted by the elements of a given eigenvector, thus resulting in a total of k index traits. The first index trait corresponds to the first eigenvector with the largest eigen value, the second index trait corresponds to the second eigenvector with the second largest eigen value, and so on. In statistical terms, these "synthetic" index traits are principal components. The variance of index IK is
2IK = b'E'KEb = b'Db where D is a diagonal matrix with the eigen values of K as the diagonal elements (Searle, 1966), indicating that the index traits are uncorrelated. The total variance of these index traits is
with 1 being the summing vector of order k.
Construction of Individual Eigenvector Indexes
The K matrix used in this study was estimated using a quartic Legendre polynomial (k = 5) under a RR test-day animal model with the first-lactation milk of Japanese Holstein cows (Togashi et al., 2005). A quartic Legendre polynomial (k = 5) resulted in 5 eigenvectors and thus, 5 index traits. To study the characteristics of each eigenvector of K, the genetic response associated with each individual eigenvector was computed. Each index trait was used separately as a selection criterion (I(ith) =
'ei where ith = first, second, third, fourth, or fifth) to assess the impact of the ith eigenvector on lactation milk and persistency. Let g be a (301 x 1) vector of genetic values from DIM 5 to 305. The genetic responses from DIM 5 through 305 (
) to selection on I(ith) =
'ei is as follows:
![]() |
where
= [
G5
G6 · · ·
G305]', SD is selection differential,
is selection intensity, and
2I(ith) = ei'Kei =
i.
Sequential Eigenvector Indexes for Improving Lactation Milk
Let the breeding goal be to improve the genetic value of lactation milk (GL = 1'
). Then, the eigenvector index I can be derived as follows:
![]() |
Differentiating the function f with respect to b and setting the resulting derivatives equal to zero lead to the following:
![]() |
Thus,
![]() | [1] |
The genetic responses (
) of daily yields from DIM 5 to 305 are
![]() | [2] |
where
2I = b'Db.
The full eigenvector index denoted by I(5) consists of 5 index traits. The last index trait of the full eigenvector index was dropped sequentially to yield 4 reduced indexes I(4), I(3), I(2), and I(1) where the number in parenthesis indicates the number of index traits included in an index. For example, the last index trait of I(5) was dropped to yield I(4); the last index trait of I(4) was dropped to produce I(3); and so on. The full and reduced indexes were computed according to [1]. Because the index traits are orthogonal, the index coefficient for a given index trait is the same among these sequential indexes (e.g., b2 for the second index trait in I(2), I(3), I(4), and I(5) are identical). Genetic response to each of these sequential indexes was computed based on [2].
Sequential Eigenvector Index for Improving Lactation Milk and Persistency
Let g280 and g55 be the genetic values at DIM 280 and 55, respectively, and a1 and a2 be the economic weights of lactation milk and persistency. The eigenvector index (I*) designed to maximize a linear combination of lactation milk and persistency is defined as
![]() |
and net merit:
![]() |
where g* is a (299 x 1) vector of genetic values from DIM 5 to 305 excluding g280 and g55. The exclusion of g280 and g55 from g* is to avoid duplication because the second trait of the net merit (H) contains both g280 and g55. The measure of persistency is defined as the difference between g280 and g55. Togashi and Lin (2004b) compared the efficiency of 5 different selection criteria for persistency and found that selection on the difference in EBV between the peak and DIM 280 achieved the greatest persistency. First-lactation milk of the Japanese Holstein cows peaked at DIM 55 (Togashi et al., 2005). The eigenvector index designed to maximize net merit is derived as follows:
![]() |
where G* is a (299 x 299) submatrix of G by deleting the rows and columns of G corresponding to DIM 55 and 280; vectors G55 and G280 refer to the columns of G corresponding to DIM 55 and 280, respectively;
55 and
280 refer to the rows of
corresponding to DIM 55 and 280, respectively;
* is a (299 x k) matrix obtained by deleting the rows of
corresponding to DIM 55 and 280; and a1 and a2 are economic weights between lactation milk and persistency. Because the relative economic weights between lactation milk and persistency were not available in literature, their economic weights were assumed to be a1:a2 = 1:50 and a1:a2 = 1:100. Taking the derivative of the function f with respect to b and setting the partial derivatives to zero result in the following equations:
![]() |
Thus, the index for maximizing the response in net merit of lactation milk and persistency is,
![]() | [3] |
The genetic responses (
) of daily yields from DIM 5 to 305 are computed as follows:
![]() | [4] |
The fitting of quartic Legendre polynomial (k = 5) results in a full index, I*(5), that consists of 5 index traits resulting from the 5 eigenvectors. Dropping the last index trait sequentially produced the reduced indexes I*(4), I* (3), I*(2), and I*(1). Genetic responses to this series of sequential indexes were computed according to [4]. Selection intensity is set at 1.0 for all selection criteria compared in this study. The genetic responses to different indexes compared were computed and used for direct comparisons.
Numerical Example for Constructing Eigenvector Index
The construction of an eigenvector index for improving lactation milk and persistency with a1:a2 = 1:50 may be illustrated. The additive genetic RR covariance matrix (K) was estimated using a quartic Legendre polynomial (k = 5) under a test-day animal model (Togashi et al., 2005):
![]() |
Eigen values and eigenvectors of matrix K were shown in Table 1
. It follows that
|
![]() |
![]() |
is a (301 x 5) matrix of Legendre polynomial coefficients evaluated from DIM 5 to 305;
![]() |
55 and
280 refer to the rows of
corresponding to DIM 55 and 280, respectively:
![]() |
![]() |
and
* is a (299 x 5) matrix formed by deleting
55 and
280 from
.
Substituting the relevant matrix elements into [3] gives the following index coefficients,
![]() |
Therefore, I* = 213.45x1 + 71.40x2 + 60.44x3 24.91x4 + 20.40x5
where x1 = e'1
, x2 = e'2
, x3 = e'3
, x4 = e' 4
, and x5 = e'5
with ei being the ith eigenvector of K and
= (
0
1
2
3
4)' being the additive genetic RR coefficients of each animal.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
Genetic responses in lactation milk and persistency to the 5 individual eigenvector indexes (I(first), I(second), I(third), I(fourth), and I(fifth)) were compared with conventional selection on lactation milk in Table 2
. Genetic responses to the first eigenvector index (I(first)) and to selection on lactation milk (IGL ) were similar. The second eigenvector index decreased both lactation milk and peak yield but increased milk at DIM 280, thus increasing persistency. The third eigenvector index slightly increased both lactation milk and persistency whereas the fourth eigenvector index decreased both lactation milk and persistency. The genetic responses to the fifth eigenvector index in lactation milk and persistency were negligible.
|
|
Genetic responses in daily milk to the fourth and the fifth eigenvector indexes hovered around zero across lactation (Figure 1
), giving a combined gain of 2.3 kg in lactation milk. Furthermore, these 2 individual eigenvector indexes produced little response in persistency (Table 2
). The combined results indicate that both the fourth and fifth eigenvectors play a minimal role in genetic improvement of lactation milk and persistency, indicating that the fourth and fifth eigenvectors can be excluded from the analysis. Note that Legendre polynomial terms 4 and 5 do not correspond to eigenvectors 4 and 5. Therefore, it makes no difference in fitting quadratic (k = 3) or quartic (k = 5) Legendre polynomials to Japanese test-day data. Obviously, the use of a quadratic rather than a quartic Legendre polynomial would greatly reduce computational requirements under a test-day animal model. In this regard, Druet et al. (2003) pointed out that the model using the eigenvectors of the additive genetic covariance matrix would reduce the computational requirements and improve its convergence properties. The degree of the additive genetic RR Legendre polynomial fitted to a given data set should be determined not only by computational costs but genetic responses as well. It merits further study to utilize the genetic response pattern of individual eigenvectors to modify the shape of the lactation curve.
Sequential Eigenvector Indexes Designed to Maximize Lactation Milk
Genetic responses (kg) to a series of sequential eigenvector indexes designed to maximize lactation milk are in Table 3
. Genetic responses to each of the sequential eigenvector indexes constructed by dropping one eigenvector at a time (I(5), I(4), I(3), I(2), I(1)) are almost the same, suggesting that extra genetic gains by adding extra eigenvectors to the first (leading) eigenvector index are minimal when the breeding goal is to improve lactation milk alone. Index coefficient of 211.45 for the first eigenvector is much larger than any other 4 eigenvectors (Table 3
). Because the index traits are orthogonal, the genetic responses in the 5 index traits are additive and the variance of I(5) is equal to the sum of the variances for the 5 individual indexes.
|
H) to selection on I*(5) was 754.9 as compared with 749.7 for conventional selection on lactation EBV (IGL). When a1:a2 = 1:100, selection on I*(5) yielded
H = 780.5 vs. 759.9 for selection on IGL . Therefore, the advantage of I*(5) over IGL increases as the economic weight for persistency increases. Table 3
H increased as the number of eigenvectors used for index construction increased and approached the maximum when the first 3 eigenvectors were used to derive the index. The use of the first 2 eigenvectors explained 99.78% of the variation in the breeding goal of improving the lactation milk and persistency as opposed to 99.99% for the first 3 eigenvectors, suggesting that the use of the first 2 eigenvectors is sufficient. Although the contribution of the third eigenvector is small in this case, the importance of the third eigenvector increases with increasing economic weight on persistency; thus, it is worthwhile using the first 3 eigenvectors for the purpose of general application. The joint contribution of the fourth and fifth eigenvectors to the breeding goal is negligible. In this regard, Druet et al. (2003) suggested that using the reduced set of eigenvectors for genetic evaluation would reduce computational cost. Van der Werf et al. (1998) used the eigen values and eigenvectors to perform canonical transformation to reduce the number of RR coefficients needed and the number of equations in the mixed model equations for computational efficiency. | CONCLUSIONS |
|---|
|
|
|---|
H due to the inclusion of the second and third eigenvector in the index increases with increasing economic weight for persistency in relation to lactation milk. The first 3 eigenvectors are necessary to construct an index to maximize
H and the fourth and the fifth eigenvectors are negligible. This study demonstrates the construction of various eigenvector indexes based on a first-parity RR test-day model, but the same principle is applicable to a multiple-lactation, test-day model.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication September 20, 2005. Accepted for publication July 6, 2006.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
B. Albarran-Portillo and G. E. Pollott Genetic Parameters Derived From Using a Biological Model of Lactation on Records of Commercial Dairy Cows J Dairy Sci, September 1, 2008; 91(9): 3639 - 3648. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Togashi and C. Y. Lin Genetic Improvement of Total Milk Yield and Total Lactation Persistency of the First Three Lactations in Dairy Cattle J Dairy Sci, July 1, 2008; 91(7): 2836 - 2843. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Togashi and C. Y. Lin Genetic Modification of the Lactation Curve by Bending the Eigenvectors of the Additive Genetic Random Regression Coefficient Matrix J Dairy Sci, December 1, 2007; 90(12): 5753 - 5758. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. P. Berry, B. Horan, M. O'Donovan, F. Buckley, E. Kennedy, M. McEvoy, and P. Dillon Genetics of Grass Dry Matter Intake, Energy Balance, and Digestibility in Grazing Irish Dairy Cows J Dairy Sci, October 1, 2007; 90(10): 4835 - 4845. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |