Genomic Selection Using Single-Step Genomic BLUP on the Number of Services per Conception Trait in Thai-Holstein Crossbreeds.

Reproductive traits are important traits that directly affect a farmer's income and are difficult to improve upon using traditional genetic methods. Therefore, there is a need to consider new options for increasing the accuracy of the genetic selection of dairy cows. The objective of this study was to compare the genetic methods of the traditional BLUP and ssGBLUP techniques in terms of the estimated genetic parameters and accuracy of the estimated breeding values. The data comprised 101,331 services per conception (NSPC) records from 54,027 Thai-Holstein crossbred cows, 109,233 pedigree data, and 770 genotyped animals. A Bayesian analysis via threshold Gibbs sampling was used to analyze the estimated variance components and genetic parameters. The results showed that the means of the NSPC data were 2.21, 2.31, and 2.42 for less than 87.5% for Holstein genetics (breed group; BG1), 87.5 to 93.6% for Holstein genetics (BG2), and greater than 93.7% for Holstein genetics (BG3), respectively. The estimated heritability values were 0.038 and 0.051, and the repeatability values were 0.149 and 0.157 for the traditional BLUP and ssGBLUP methods, respectively. The accuracy of the estimated breeding values from the ssGBLUP method was higher than that from the traditional BLUP method, ranging from 6.05 to 17.69%, depending on the dataset, especially in the top 20% of the bull dataset had the highest values. In conclusion, the ssGBLUP method could improve the heritability value and accuracy of the breeding values compared with the traditional BLUP method. Therefore, switching from traditional methods to the ssGBLUP method for the Thai dairy cattle breeding program is a viable option.

Determining Heat Stress Effects of Multiple Genetic Traits in Tropical Dairy Cattle Using Single-Step Genomic BLUP.

Heat stress is becoming a significant problem in dairy farming, especially in tropical countries, making accurate genetic selection for heat tolerance a priority. This study investigated the effect of heat stress manifestation on genetics for milk yield, milk quality, and dairy health traits with and without genomic information using single-step genomic best linear unbiased prediction (ssGBLUP) and BLUP in Thai-Holstein crossbred cows. The dataset contained 104,150 test-day records from the first lactation of 15,380 Thai-Holstein crossbred cows. A multiple-trait random regression test-day model on a temperature-humidity index (THI) function was used to estimate the genetic parameters and genetic values. Heat stress started at a THI of 76, and the heritability estimates ranged from moderate to low. The genetic correlation between those traits and heat stress in both BLUP methods was negative. The accuracy of genomic predictions in the ssGBLUP method was higher than the BLUP method. In conclusion, heat stress negatively impacted milk production, increased the somatic cell score, and disrupted the energy balance. Therefore, in dairy cattle genetic improvement programs, heat tolerance is an important trait. The new genetic evaluation method (ssGBLUP) should replace the traditional method (BLUP) for more accurate genetic selection.

Validation of single-step genomic predictions using the linear regression method for milk yield and heat tolerance in a Thai-Holstein population.

BACKGROUND AND AIM: Genomic selection improves accuracy and decreases the generation interval, increasing the selection response. This study was conducted to assess the benefits of using single-step genomic best linear unbiased prediction (ssGBLUP) for genomic evaluations of milk yield and heat tolerance in Thai-Holstein cows and to test the value of old phenotypic data to maintain the accuracy of predictions. MATERIALS AND METHODS: The dataset included 104,150 milk yield records collected from 1999 to 2018 from 15,380 cows. The pedigree contained 33,799 animals born between 1944 and 2016, of which 882 were genotyped. Analyses were performed with and without genomic information using ssGBLUP and BLUP, respectively. Statistics for bias, dispersion, the ratio of accuracies, and the accuracy of estimated breeding values were calculated using the linear regression (LR) method. A partial dataset excluded the phenotypes of the last generation, and 66 bulls were identified as validation individuals. RESULTS: Bias was considerable for BLUP (0.44) but negligible (-0.04) for ssGBLUP; dispersion was similar for both techniques (0.84 vs. 1.06 for BLUP and ssGBLUP, respectively). The ratio of accuracies was 0.33 for BLUP and 0.97 for ssGBLUP, indicating more stable predictions for ssGBLUP. The accuracy of predictions was 0.18 for BLUP and 0.36 for ssGBLUP. Excluding the first 10 years of phenotypic data (i.e., 1999-2008) decreased the accuracy to 0.09 for BLUP and 0.32 for ssGBLUP. Genomic information doubled the accuracy and increased the persistence of genomic estimated breeding values when old phenotypes were removed. CONCLUSION: The LR method is useful for estimating accuracies and bias in complex models. When the population size is small, old data are useful, and even a small amount of genomic information can substantially improve the accuracy. The effect of heat stress on first parity milk yield is small.

Estimation of genetic parameters and trends for production traits of dairy cattle in Thailand using a multiple-trait multiple-lactation test day model.

OBJECTIVE: The objective of this study was to estimate the genetic parameters and trends for milk, fat, and protein yields in the first three lactations of Thai dairy cattle using a 3-trait,- 3-lactation random regression test-day model. METHODS: Data included 168,996, 63,388, and 27,145 test-day records from the first, second, and third lactations, respectively. Records were from 19,068 cows calving from 1993 to 2013 in 124 herds. (Co) variance components were estimated by Bayesian methods. Gibbs sampling was used to obtain posterior distributions. The model included herd-year-month of testing, breed group-season of calving-month in tested milk group, linear and quadratic age at calving as fixed effects, and random regression coefficients for additive genetic and permanent environmental effects, which were defined as modified constant, linear, quadratic, cubic and quartic Legendre coefficients. RESULTS: Average daily heritabilities ranged from 0.36 to 0.48 for milk, 0.33 to 0.44 for fat and 0.37 to 0.48 for protein yields; they were higher in the third lactation for all traits. Heritabilities of test-day milk and protein yields for selected days in milk were higher in the middle than at the beginning or end of lactation, whereas those for test-day fat yields were high at the beginning and end of lactation. Genetics correlations (305-d yield) among production yields within lactations (0.44 to 0.69) were higher than those across lactations (0.36 to 0.68). The largest genetic correlation was observed between the first and second lactation. The genetic trends of 305-d milk, fat and protein yields were 230 to 250, 25 to 29, and 30 to 35 kg per year, respectively. CONCLUSION: A random regression model seems to be a flexible and reliable procedure for the genetic evaluation of production yields. It can be used to perform breeding value estimation for national genetic evaluation in the Thai dairy cattle population.

Genetic associations between milk fat-to-protein ratio, milk production and fertility in the first two lactations of Thai Holsteins dairy cattle.

The aims of this study were to estimate, simultaneously, the genetic parameters of test-day milk fat-to-protein ratio (FPR), test-day milk yield (MY), and days-open (DO) in the first two lactations of Thai Holsteins. A total of 76 194 test-day production records collected from 8874 cows with 8674 DO records between 2001 and 2011 from different lactations were treated as separated traits. The estimates of heritability for test-day FPR in the first lactation showed an increasing trend, whereas the estimates in the second lactation showed a U-shape trend. Genetic correlations for FPR-DO and MY-DO showed a decreasing trend along days in milk (DIM) in both lactations, whereas genetic correlations for FPR-MY increased along DIM in the first lactation but decreased in the second lactation. Genetic correlations of FPR between consecutive DIM were moderate to high, which showed the effectiveness of simultaneous analyses. Selection of FPR in the early stage has no adverse effect on MY and DO for the first lactation but has a negative effect on MY and positive effect on DO for the second lactation. This study showed that genetic improvement of the energy balance using FPR, MY and DO with multi-trait test day model could be applied in a Thailand dairy cattle breeding program.

Genetic analysis of the rates of conception using a longitudinal threshold model with random regression in dairy crossbreeding within a tropical environment.

This study was designed to: (i) estimate genetic parameters and breeding values for conception rates (CR) using the repeatability threshold model (RP-THM) and random regression threshold models (RR-THM); and (ii) compare covariance functions for modeling the additive genetic (AG) and permanent environmental (PE) effects in the RR-THM. The CR was defined as the outcome of an insemination. A data set of 130 592 first-lactation insemination records of 55 789 Thai dairy cows, calving between 1996 and 2011, was used in the analyses. All models included fixed effects of year × month of insemination, breed × day in milk to insemination class and age at calving. The random effects consisted of herd × year interaction, service sire, PE, AG and residual. Variance components were estimated using a Bayesian method via Gibbs sampling. Heritability estimates of CR ranged from 0.032 to 0.067, 0.037 to 0.165 and 0.045 to 0.218 for RR-THM with the second, third and fourth-order of Legendre polynomials, respectively. The heritability estimated from RP-THM was 0.056. Model comparisons based on goodness of fit, predictive abilities, predicted service results of animal, and pattern of genetic parameter estimates, indicated that the model which fit the desired outcome of insemination was the RR-THM with two regression coefficients.

Genetic relationships of fertility traits with test-day milk yield and fat-to-protein ratio in tropical smallholder dairy farms.

The test-day milk fat-to-protein ratio (TD-FPR) could serve as a measure of energy balance status and might be used as a criterion to improve metabolic stability and fertility through genetic selection. Therefore, genetic parameters for fertility traits, test-day milk yield (TD-MY) and TD-FPR, as well as, their relationships during different stages of lactation, were estimated on data collected from 25 968 primiparous Thai dairy crossbred cows. Gibbs sampling algorithms were implemented to obtain (co)variance components using both univariate linear and threshold animal models and bivariate linear-linear and linear-threshold animal models with random regression. Average TD-MY and TD-FPR were 12.60 and 1.15. Heritability estimates for TD-MY, TD-FPR and selected fertility traits ranged from 0.31 to 0.58, 0.17 to 0.19 and 0.02 to 0.05, respectively. Genetic correlations among TD-FPR and TD-MY, TD-FPR and fertility traits, and TD-MY and fertility traits ranged from 0.05 to -0.44, from -0.98 to 0.98 and -0.22 to 0.79, respectively. Selection for lower TD-FPR would decrease numbers of inseminations per conception and increase conception at first service and pregnancy within 90 days. In addition, cow selection based only on high milk production has strong effects to prolong days to first service, days open and calving interval.