Genetic parameters of milk and lactation curve traits of dairy cattle from research farms in Thailand

Objective: This study was aimed to estimate the genetic parameters, including genetic and phenotypic correlations, of milk yield, lactation curve traits and milk composition of Thai dairy cattle from three government research farms. Methods: The data of 25,789 test-day milk yield and milk composition records of 1,468 cattle from lactation 1 to 3 of Holstein Friesian (HF) and crossbred HF dairy cattle calved between 1990 and 2015 from three government research farms in Thailand were analysed. 305-day milk yield was estimated by the Wood model and a test interval method. The Wood model was used for estimating cumulative 305-day milk yield, peak milk yield, days to peak milk yield and persistency. Genetic parameters were estimated using linear mixed models with herd, breed group, year and season of calving as fixed effects, and animals linked to a pedigree as random effects, together with a residual error. Univariate models were used to estimate variance components, heritability, estimated breeding values (EBVs) and repeatability of each trait, while pairwise bivariate models were used to estimate covariance components and correlations between traits in the same lactation and in the same trait across lactations. Results: The heritability of 305-day milk yield, peak milk yield and protein percentage have moderate to high estimates ranging from 0.19 to 0.45 while days to peak milk yield, persistency and fat percentage have low heritability ranging from 0.08 to 0.14 in lactation 1 cows. Further, heritability of most traits considered was higher in lactation 1 compared with lactations 2 and 3. For cows in lactation 1, high genetic correlations were found between 305-day milk yield and peak milk yield (0.86±0.07) and days to peak milk yield and persistency (0.99±0.02) while estimates of genetic correlations between the remaining traits were imprecise due to the high standard errors. The genetic correlations within the traits across lactation were high. There was no consistent trend of EBVs for most traits in the first lactation over the study period. Conclusion: Both the Wood model and test interval method can be used for milk yield estimates in these herds. However, the Wood model has advantages over the test interval method as it can be fitted using fewer test-day records and the estimated model parameters can be used to derive estimates of other lactation curve parameters. Milk yield, peak milk yield and protein percentage can be improved by a selection and mating program while days to peak milk yield, persistency and fat percentage can be improved by including into a selection index.