The objectives of the study were to estimate genetic parameters for milk production traits of Holstein cattle using random regression models (RRMs), and to compare the goodness of fit of various RRMs with homogeneous and heterogeneous residual variances. A total of 126,980 test-day milk production records of the first parity Holstein cows between 2007 and 2014 from the Dairy Cattle Improvement Center of National Agricultural Cooperative Federation in South Korea were used. These records included milk yield (MILK), fat yield (FAT), protein yield (PROT), and solids-not-fat yield (SNF). The statistical models included random effects of genetic and permanent environments using Legendre polynomials (LP) of the third to fifth order (L3–L5), fixed effects of herd-test day, year-season at calving, and a fixed regression for the test-day record (third to fifth order). The residual variances in the models were either homogeneous (HOM) or heterogeneous (15 classes, HET15; 60 classes, HET60). A total of nine models (3 orders of polynomials×3 types of residual variance) including L3-HOM, L3-HET15, L3-HET60, L4-HOM, L4-HET15, L4-HET60, L5-HOM, L5-HET15, and L5-HET60 were compared using Akaike information criteria (AIC) and/or Schwarz Bayesian information criteria (BIC) statistics to identify the model(s) of best fit for their respective traits. The lowest BIC value was observed for the models L5-HET15 (MILK; PROT; SNF) and L4-HET15 (FAT), which fit the best. In general, the BIC values of HET15 models for a particular polynomial order was lower than that of the HET60 model in most cases. This implies that the orders of LP and types of residual variances affect the goodness of models. Also, the heterogeneity of residual variances should be considered for the test-day analysis. The heritability estimates of from the best fitted models ranged from 0.08 to 0.15 for MILK, 0.06 to 0.14 for FAT, 0.08 to 0.12 for PROT, and 0.07 to 0.13 for SNF according to days in milk of first lactation. Genetic variances for studied traits tended to decrease during the earlier stages of lactation, which were followed by increases in the middle and decreases further at the end of lactation. With regards to the fitness of the models and the differential genetic parameters across the lactation stages, we could estimate genetic parameters more accurately from RRMs than from lactation models. Therefore, we suggest using RRMs in place of lactation models to make national dairy cattle genetic evaluations for milk production traits in Korea.

The study was conducted to analyze the genetic parameters of somatic cell score (SCS) of Holstein cows, which is an important indicator to udder health. Test-day records of somatic cell counts (SCC) of 305-day lactation design from first to fifth lactations were collected on Holsteins in Korea during 2000 to 2012. Records of animals within 18 to 42 months, 30 to 54 months, 42 to 66 months, 54 to 78 months, and 66 to 90 months of age at the first, second, third, fourth and fifth parities were analyzed, respectively. Somatic cell scores were calculated, and adjusted for lactation production stages by Wilmink’s function. Lactation averages of SCS (LSCS1 through LSCS5) were derived by further adjustments of each test-day SCS for five age groups in particular lactations. Two datasets were prepared through restrictions on number of sires/herd and dams/herd, progenies/sire, and number of parities/cow to reduce data size and attain better relationships among animals. All LSCS traits were treated as individual trait and, analyzed through multiple-trait sire models and single trait animal models via VCE 6.0 software package. Herd-year was fitted as a random effect. Age at calving was regressed as a fixed covariate. The mean LSCS of five lactations were between 3.507 and 4.322 that corresponded to a SCC range between 71,000 and 125,000 cells/mL; with coefficient of variation from 28.2% to 29.9%. Heritability estimates from sire models were within the range of 0.10 to 0.16 for all LSCS. Heritability was the highest at lactation 2 from both datasets (0.14/0.16) and lowest at lactation 5 (0.11/0.10) using sire model. Heritabilities from single trait animal model analyses were slightly higher than sire models. Genetic correlations between LSCS traits were strong (0.62 to 0.99). Very strong associations (0.96 to 0.99) were present between successive records of later lactations. Phenotypic correlations were relatively weaker (<0.55). All correlations became weaker at distant lactations. The estimated breeding values (EBVs) of LSCS traits were somewhat similar over the years for a particular lactation, but increased with lactation number increment. The lowest EBV in first lactation indicated that selection for SCS (mastitis resistance) might be better with later lactation records. It is expected that results obtained from these multi-trait lactation model analyses, being the first large scale SCS data analysis in Korea, would create a good starting step for application of advanced statistical tools for future genomic studies focusing on selection for mastitis resistance in Holsteins of Korea.