This study investigated the effects of heat stress on milk production in Korean Holstein cows using large-scale data. Heat stress was assessed using the temperature-humidity index (THI). Weather records (2016 to 2020) were collected from 70 regional weather stations using an installed automated surface observing system (ASOS). A dataset of 2,094,436 milk production records from 215,276 Holstein cows obtained from the Dairy Cattle Genetic Improvement Center was analyzed. Stepwise selection was used to select the input variables, including the daily maximum THI (THI_max). Least-squares means were calculated for milk yield, fat and protein corrected milk (FPCM), fat and protein yield, fat-to-protein ratio, solids not fat, and lactation persistency. Segmented linear regression analysis determined the break points (BPs) of the THI_max. Over the five years, heat stress exposure increased, particularly from May to September. This study identified BPs around THI_max of 80–82 for milk yield and FPCM. Similar patterns for other milk traits were observed, which significantly decreased beyond their respective BPs. These findings indicate that THI variations adversely affect milk yield and composition in dairy cows, highlighting the importance of appropriate feeding management strategies to ensure the optimal productivity of Holstein cows under varying climatic conditions.

Objective: This study was conducted to test the efficiency of genomic selection for milk production traits in a Korean Holstein cattle population.Methods: A total of 506,481 milk production records from 293,855 animals (2,090 heads with single nucleotide polymorphism information) were used to estimate breeding value by single step best linear unbiased prediction.Results: The heritability estimates for milk, fat, and protein yields in the first parity were 0.28, 0.26, and 0.23, respectively. As the parity increased, the heritability decreased for all milk production traits. The estimated generation intervals of sire for the production of bulls (LSB) and that for the production of cows (LSC) were 7.9 and 8.1 years, respectively, and the estimated generation intervals of dams for the production of bulls (LDB) and cows (LDC) were 4.9 and 4.2 years, respectively. In the overall data set, the reliability of genomic estimated breeding value (GEBV) increased by 9% on average over that of estimated breeding value (EBV), and increased by 7% in cows with test records, about 4% in bulls with progeny records, and 13% in heifers without test records. The difference in the reliability between GEBV and EBV was especially significant for the data from young bulls, i.e. 17% on average for milk (39% vs 22%), fat (39% vs 22%), and protein (37% vs 22%) yields, respectively. When selected for the milk yield using GEBV, the genetic gain increased about 7.1% over the gain with the EBV in the cows with test records, and by 2.9% in bulls with progeny records, while the genetic gain increased by about 24.2% in heifers without test records and by 35% in young bulls without progeny records.Conclusion: More genetic gains can be expected through the use of GEBV than EBV, and genomic selection was more effective in the selection of young bulls and heifers without test records.

Holsteins are known as the world's highest-milk producing dairy cattle. The purpose of this study was to identify genetic regions strongly associated with milk traits (milk production, fat, and protein) using Korean Holstein data. This study was performed using single nucleotide polymorphism (SNP) chip data (Illumina BovineSNP50 Beadchip) of 911 Korean Holstein individuals. We inferred each genomic estimated breeding values based on best linear unbiased prediction (BLUP) and ridge regression using BLUPF90 and R. We then performed a genome-wide association study and identified genetic regions related to milk traits. We identified 9, 6, and 17 significant genetic regions related to milk production, fat and protein, respectively. These genes are newly reported in the genetic association with milk traits of Holstein. This study complements a recent Holstein genome-wide association studies that identified other SNPs and genes as the most significant variants. These results will help to expand the knowledge of the polygenic nature of milk production in Holsteins.

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.

Objective: The objective of this study was to investigate the genetic components of daily milk yield and to re-rank bulls in South Korea by estimated breeding value (EBV) under heat stress using the temperature-humidity index (THI). Methods: This study was conducted using 125,312 monthly test-day records, collected from January 2000 to February 2017 for 19,889 Holstein cows from 647 farms in South Korea. Milk production data were collected from two agencies, the Dairy Cattle Genetic Improvement Center and the Korea Animal Improvement Association, and meteorological data were obtained from 41 regional weather stations using the Automated Surface Observing System (ASOS) installed throughout South Korea. A random regression model using the THI was applied to estimate genetic parameters of heat tolerance based on the test-day records. The model included herd-year-season, calving age, and days-in-milk as fixed effects, as well as heat tolerance as an additive genetic effect, permanent environmental effect, and direct additive and permanent environmental effect. Results: Below the THI threshold (<= 72; no heat stress), the variance in heat tolerance was zero. However, the heat tolerance variance began to increase as THI exceeded the threshold. The covariance between the genetic additive effect and the heat tolerance effect was -0.33. Heritability estimates of milk yield ranged from 0.111 to 0.176 (average: 0.128). Heritability decreased slightly as THI increased, and began to increase at a THI of 79. The predicted bull EBV ranking varied with THI. Conclusion: We conclude that genetic evaluation using the THI function could be useful for selecting bulls for heat tolerance in South Korea.

The objectives of this study were to compare identified informative regions through two genome-wide association study (GWAS) approaches and determine the accuracy and bias of the direct genomic value (DGV) for milk production traits in Korean Holstein cattle, using two genomic prediction approaches: single-step genomic best linear unbiased prediction (ss-GBLUP) and Bayesian Bayes-B. Records on production traits such as adjusted 305-day milk (MY305), fat (FY305), and protein (PY305) yields were collected from 265,271 first parity cows. After quality control, 50,765 single-nucleotide polymorphic genotypes were available for analysis. In GWAS for ss-GBLUP (ssGWAS) and Bayes-B (BayesGWAS), the proportion of genetic variance for each 1-Mb genomic window was calculated and used to identify informative genomic regions. Accuracy of the DGV was estimated by a five-fold cross-validation with random clustering. As a measure of accuracy for DGV, we also assessed the correlation between DGV and deregressed-estimated breeding value (DEBV). The bias of DGV for each method was obtained by determining regression coefficients. A total of nine and five significant windows (1 Mb) were identified for MY305 using ssGWAS and BayesGWAS, respectively. Using ssGWAS and BayesGWAS, we also detected multiple significant regions for FY305 (12 and 7) and PY305 (14 and 2), respectively. Both single-step DGV and Bayes DGV also showed somewhat moderate accuracy ranges for MY305 (0.32 to 0.34), FY305 (0.37 to 0.39), and PY305 (0.35 to 0.36) traits, respectively. The mean biases of DGVs determined using the single-step and Bayesian methods were 1.50±0.21 and 1.18±0.26 for MY305, 1.75±0.33 and 1.14±0.20 for FY305, and 1.59±0.20 and 1.14±0.15 for PY305, respectively. From the bias perspective, we believe that genomic selection based on the application of Bayesian approaches would be more suitable than application of ss-GBLUP in Korean Holstein populations.