Genetic parameters for udder conformation traits derived from Cartesian coordinates generated by robotic milking systems in North American Holstein cattle

Udder conformation is directly related to milk yield, cow health, workability, and welfare. Automatic milking systems (AMS, also known as milking robots) have become popular worldwide, and the number of dairy farms adopting these systems have increased considerably over the past years. In each milking visit, AMS record the location of the 4 teats as Cartesian coordinates in a xyz plan, which can then be used to derive udder conformation traits. AMS generate a large amount of per milking visit data for individual cows, which contribute to an accurate assessment of important traits such as udder conformation without the addition of human classifier errors (in subjective scoring systems). Therefore, the primary objectives of this study were to estimate genomic-based genetic parameters for udder conformation traits derived from AMS records in North American Holstein cattle and to assess the genetic correlation between the derived traits for evaluating the feasibility of multi-trait genomic selection for breeding cows that are more suitable for milking in AMS. The Cartesian teat coordinates measured during each milking visit were collected by 36 milking robots in 4,480 Holstein cows from 2017 to 2021, resulting in 5,317,488 records. A total of 4,118 of these Holstein cows were also genotyped for 57,600 single nucleotide polymorphisms. Five udder conformation traits were derived: udder balance (UB, mm), udder depth (UD, mm), front teat distance (FTD, mm), rear teat distance (RTD, mm), and distance front-rear (DFR, mm). In addition, 2 traits directly related to cow productivity in the system were added to the study: daily milk yield (DY) and milk electroconductivity (EC; as an indicator of mastitis). Variance components and genetic parameters for UB, UD, FTD, RTD, DFR, DY, and EC were estimated based on repeatability animal models. The estimates of heritability (±standard error, SE) for UB, UD, FTD, RTD, DFR, DY, and EC were 0.41 ± 0.02, 0.79 ± 0.01, 0.53 ± 0.02, 0.40 ± 0.02, 0.65 ± 0.02, 0.20 ± 0.02, and 0.46 ± 0.02, respectively. The repeatability estimates (±SE) for UB, UD, FTD, RTD, and DFR were 0.82 ± 0.01, 0.93 ± 0.01, 0.87 ± 0.01, 0.83 ± 0.01, and 0.88 ± 0.01, respectively. The strongest genetic correlations were observed between the FTD and RTD (0.54 ± 0.03), UD and DFR (−0.47 ± 0.03), DFR and FTD (0.32 ± 0.03), and UD and FTD (−0.31 ± 0.03). These results suggest that udder conformation traits derived from Cartesian coordinates from AMS are moderately to highly heritable. Furthermore, the moderate genetic correlations between these traits should be considered when developing selection sub-indexes. The most relevant genetic correlations between traits related to cow milk productivity and udder conformation traits were between UD and EC (−0.25 ± 0.03) and between DFR and DY (0.30 ± 0.04), in which both genetic correlations are favorable. These findings will contribute to the design of genomic selection schemes for improving udder conformation in North American Holstein cattle, especially in precision dairy farms.