Genetic characteristics of Japanese Holstein cows based on multiple-lactation random regression test-day animal models

The first three lactation curves of the Japanese Holstein cows were analyzed using a random regression (RR) test-day model with a cubic Legendre polynomial fitted to each of the three parities. The first three eigenvectors of the additive genetic RR covariance matrix explained 77.8, 10.9, and 4.2% of the total variance of the three parities and are associated mainly with the level of milk yield, the linear increase, and the concave curve, respectively. On a lactational basis, as the parity increases, the contribution of the first eigenvector to a lactational variation decreases whereas the contribution of the second eigenvector increases sharply. This means that the impact of the first eigenvector on the level of milk production decreases across parity whereas the effect of the second eigenvector on the shape of the lactation curve increases across parity. The first lactation curve was the most persistent, followed by the second and the third lactation. Persistency and days to reach peak yield decrease as the parity increases (45, 40, and 36 days for the first three parities). Daily heritabilities within lactation were lower for the first parity than for the second or the third parity. The first three lactation curves possess distinctive genetic characteristics that merit consideration when combining the proofs of the first three lactations to select for lifetime production. Within- and between-parity genetic correlations between the constant and the linear RR coefficients were all positive, suggesting that raising the level of milk production in one parity would increase the linear slope in all parities, thus improving persistency. Within- and between-parity genetic correlations between the constant and the quadratic RR coefficients were all negative, implying that increasing the level of production in one parity would deepen and/or widen the concave curve in all parities, thus decreasing persistency. The linear and quadratic RR coefficients were negatively correlated within or between parities and thus have antagonistic effects on persistency.