Genetic Evaluation of Resilience Indicators in Holstein Cows

The analysis of resilience indicators was based on daily milk yields recorded from 3347 lactations of 3080 Holstein cows located on 10 farms between 2022 and 2024. Six farms used an automatic milking system. A random regression function with a fourth-degree Legendre polynomial was used to predict the lactation curve. The indicators were the natural log-transformed variance (LnVar), lag-1 autocorrelation (r-auto), and skewness (skew) of daily milk yield (DMY) deviations from the predicted lactation curve, as well as the log-transformed variance of DMY (Var). The single-step genomic prediction method (ssGBLUP) was used for genomic evaluation. A total of 9845 genotyped animals and 36,839 SNPs were included. Heritability estimates were low (0.02-0.13). The strongest genetic correlation (0.87) was found between LnVar and Var. The genetic correlation between r-auto and skew was also strong but negative (-0.73). Resilience indicators showed a negative correlation with milk yield per lactation and a positive correlation with fat and protein contents. The negative correlation between fertility and two resilience indicators may be due to the evaluation period (50th-150th day of lactation) being when cows are most often bred after calving, and a decrease in production may accompany a significant oestrus. The associations between resilience indicators and health traits (clinical mastitis, claw health) were weak but mostly favourable.