Estimation of 305-day milk yield from test-day records of Chinese Holstein cattle

Exploring milk loss and variability during environmental perturbations across lactation stages as resilience indicators in Holstein cattle

Genetic selection for resilience is essential to improve the long-term sustainability of the dairy cattle industry, especially the ability of cows to maintain their level of production when exposed to environmental disturbances. Recording of daily milk yield provides an opportunity to develop resilience indicators based on milk losses and fluctuations in daily milk yield caused by environmental disturbances. In this context, our study aimed to explore milk loss traits and measures of variability in daily milk yield, including log-transformed standard deviation of milk deviations (Lnsd), lag-1 autocorrelation (Ra), and skewness of the deviations (Ske), as indicators of general resilience in dairy cows. The unperturbed dynamics of milk yield as well as milk loss were predicted using an iterative procedure of lactation curve modeling. Milk fluctuations were defined as a period of at least 10 successive days of negative deviations in which milk yield dropped at least once below 90% of the expected values. Genetic parameters of these indicators and their genetic correlation with economically important traits were estimated using single-trait and bivariate animal models and 8,935 lactations (after quality control) from 6,816 Chinese Holstein cows. In general, cows experienced an average of 3.73 environmental disturbances with a milk loss of 267 kg of milk per lactation. Each fluctuation lasted for 19.80 ± 11.46 days. Milk loss traits are heritable with heritability estimates ranging from 0.004 to 0.061. The heritabilities differed between Lnsd (0.135-0.250), Ra (0.008-0.058), and Ske (0.001-0.075), with the highest heritability estimate of 0.250 ± 0.020 for Lnsd when removing the first and last 10 days in milk in a lactation (Lnsd2). Based on moderate to high genetic correlations, lower Lnsd2 is associated with less milk losses, better reproductive performance, and lower disease incidence. These findings indicate that among the variables evaluated, Lnsd2 is the most promising indicator for breeding for improved resilience in Holstein cattle.

Comparison and estimation of different linear and nonlinear lactation curve submodels in random regression analyses on dairy cattle

In the random regression model (RRM) for milk yield, by replacing empirical lactation curves with the five-order Legendre polynomial to fit fixed groups, the RRM can be transformed to a hierarchical model that consisted of a RRM in the first hierarchy with Legendre polynomials as individuals’ lactation curves resolved by restricted maximum likelihood (REML) software, and a multivariate animal model for phenotypic regression coefficients in the second hierarchy resolved by DMU software. Some empirical lactation functions can be embedded into the RRM at the first hierarchy to well fit phenotypic lactation curve of the average observations across all animals. The functional relationship between each parameter and time can be described by a Legendre polynomial or an empirical curve usually called submodel, and according to three commonly used criteria, the optimal submodels were picked from linear and nonlinear submodels except for polynomials. The so-called hierarchical estimation for the RRMs in dairy cattle indicated that more biologically meaningful models were available to fit the lactation curves; moreover, with the same number of parameters, the empirical lactation curves (MIL1, MIL5, and MK1 for 3, 4, and 5 parameters, respectively) performed higher goodness of fit than Legendre polynomial when modelling individuals’ phenotypic lactation curves.