Multivariate analysis to evaluate genetic groups and production traits of crossbred Holstein × Zebu cows

The objective of this research was to try to unveil the relationship between production traits and genotypic proportions of crossbred dairy cattle using principal component analysis (PCA) and cluster analysis. The herd consists of crossbred animals of Holstein (H) and Zebu (Z) (Gir and Guzerat) in different genotypic proportions; the composition of which varies from 12.5 to 100.0 % of the genetic group H. For this study, 834 milk production records from 257 cows from the years 1997 to 2014 were analyzed. The animals were all managed at a farm located in northeastern Brazil. The variables in the PCA were total milk yield per lactation (MY), milk yield adjusted to 305 days (MY305), lactation length (LL), and proportion of H and Z breeding. This analysis reduced the size of the sample space from the original five variables to two principal components (PCs) that together explained 89.4 % of the total variation. MY, MY305, LL, and genotypic proportion of H all contributed positively to PC1. The genotypic proportion of Z contributed negatively, which established a contrast between H and Z. Further cluster analysis identified two distinct groups when considering production performance and genotype of the animals. The high-performance group was predominantly Holstein breeding, while the lower performing group consisted mostly of Zebu. Under the environmental and management conditions in which this research was conducted, the best performances for the traits considered were achieved from cows whose genotypic proportion was between 38.0 and 94.0 % Holstein breeding.