Genetic and phenotypic variation and consistency in cow preference and circadian use of robotic milking units

Many dairy herds use automatic milking stations (AMS), with cows in large herds often having access to 2 or more AMS, and must choose between them when they go for milking. Individual cows acquire routines of either consistently using a specific milking box or consistently using any available milking box. Here, we hypothesized that the degree of use of the same milking box was an expression of preference, and quantified it as preference consistency score (PCS). The PCS was calculated as a ratio between the excess frequencies of the first choice over the base frequency of "not first choice" over 15-d segments of lactation. This ratio was 0 if all choices were taken equally, and became 1.0 if only the first choice was taken in all events. We investigated the consistency of milking box preference in 2 cohorts (one Holstein and one Jersey) across 6 commercial dairy herds in Denmark (n = 4,665 cows total). In addition to PCS, we recorded and analyzed associated milking and behavior traits, including a time profile index showing use of specific clock hours when cows were milked (Time_profile, based on excess use of specific clock hours), milking frequency, time spent in the milking box, and milk yield. Records from each milking event were condensed into 15-d segments based on days in milk. The data were analyzed using a linear mixed model, with random genetic and individual cow effects, to estimate heritability (h²), repeatability (t), and individual level correlations (r_i) between traits. The average PCS was 0.43 and 0.41 in Holstein and Jersey, respectively, showing that cows developed routines for consistently using the same milking box; however, some cows had lower preference (i.e., greater flexibility in use). The Time_profile indicated that some cows were milked in a few hour-bins, whereas others were more flexible. The PCS and Time_profile traits had low heritability (h², PCS/Time_profile = 0.07 ± 0.02/0.11 ± 0.02 Holstein, 0.13 ± 0.03/0.04 ± 0.02 Jersey) and moderate repeatability (t, PCS/Time_profile = 0.47/0.40 Holstein, 0.50/0.42 Jersey). The 2 traits were weakly correlated with each other (r_i = 0.18 and 0.17), and were weakly correlated with milk yield (r_i range: 0.0 to -0.10). However, the time profile was strongly correlated with milking frequency (r_i range: -0.81 to -0.73), and was moderately correlated with daily box time (r_i range: -0.43 to -0.35). In general, Holstein and Jersey parameter estimates were of similar size, and thus in good agreement. Overall, individual cows covered a broad spectrum of preference consistency, both regarding the use of specific milking boxes and time profiles, with these 2 traits representing different aspects or dimensions of milking behavior. The findings that some cows have strong preferences for specific AMS may be most useful in herd management and farm design. The weak correlation to milk yield indicated that yield minimally affected these 2 milking associated behavior traits. In conclusion, although the traits were repeatable, heritability was low; thus, genetic selection for milk yield might minimally affect these 2 traits.