Background In this study, we assessed the accuracy of genomic prediction for carcass weight (CWT), marbling score (MS), eye muscle area (EMA) and back fat thickness (BFT) in Hanwoo cattle when using genomic best linear unbiased prediction (GBLUP), weighted GBLUP (wGBLUP), and a BayesR model. For these models, we investigated the potential gain from using pre-selected single nucleotide polymorphisms (SNPs) from a genome-wide association study (GWAS) on imputed sequence data and from gene expression information. We used data on 13,717 animals with carcass phenotypes and imputed sequence genotypes that were split in an independent GWAS discovery set of varying size and a remaining set for validation of prediction. Expression data were used from a Hanwoo gene expression experiment based on 45 animals. Results Using a larger number of animals in the reference set increased the accuracy of genomic prediction whereas a larger independent GWAS discovery dataset improved identification of predictive SNPs. Using pre-selected SNPs from GWAS in GBLUP improved accuracy of prediction by 0.02 for EMA and up to 0.05 for BFT, CWT, and MS, compared to a 50 k standard SNP array that gave accuracies of 0.50, 0.47, 0.58, and 0.47, respectively. Accuracy of prediction of BFT and CWT increased when BayesR was applied with the 50 k SNP array (0.02 and 0.03, respectively) and was further improved by combining the 50 k array with the top-SNPs (0.06 and 0.04, respectively). By contrast, using BayesR resulted in limited improvement for EMA and MS. wGBLUP did not improve accuracy but increased prediction bias. Based on the RNA-seq experiment, we identified informative expression quantitative trait loci, which, when used in GBLUP, improved the accuracy of prediction slightly, i.e. between 0.01 and 0.02. SNPs that were located in genes, the expression of which was associated with differences in trait phenotype, did not contribute to a higher prediction accuracy. Conclusions Our results show that, in Hanwoo beef cattle, when SNPs are pre-selected from GWAS on imputed sequence data, the accuracy of prediction improves only slightly whereas the contribution of SNPs that are selected based on gene expression is not significant. The benefit of statistical models to prioritize selected SNPs for estimating genomic breeding values is trait-specific and depends on the genetic architecture of each trait.

Data collected from 690 purebred Duroc pigs from 2009 to 2012 were used to estimate the heritability, and genetic and phenotypic correlations between production and meat quality traits. Variance components were obtained through the restricted maximum likelihood procedure using Wombat and SAS version 9.0. Animals were raised under the same management in five different breeding farms. The average daily gain, loin muscle area (LMA), backfat thickness (BF), and lean percent (LP) were measured as production traits. Meat quality traits included pH, cooking loss, lightness (L*), redness (a*), yellowness (b*), marbling score (MS), moisture content (MC), water holding capacity (WHC), and shear force. The results showed that the heritability estimates for meat quality traits varied largely from 0.19 to 0.79. Production traits were moderate to highly heritable from 0.41 to 0.73. Genotypically, the BF was positively correlated (p<0.05) with MC (0.786), WHC (0.904), and pH (0.328) but negatively correlated with shear force (−0.533). The results of genetic correlations indicated that selection for less BF could decrease pH, moisture content, and WHC and increase the shear force of meat. Additionally, a significant positive correlation was recorded between average daily gain and WHC, which indicates pork from faster-growing animals has higher WHC. Furthermore, selection for larger LMA and LP could increase MS and lightness color of meat. The meat quality and production traits could be improved simultaneously if desired. Hence, to avoid further deterioration of pork characteristics, appropriate selection of traits should be considered.

Whole-genome sequence (WGS) data are increasingly being applied into genomic predictions, offering a higher predictive ability by including causal mutations or single-nucleotide polymorphisms (SNPs) putatively in strong linkage disequilibrium with causal mutations affecting the trait. This study aimed to improve the predictive performance of the customized Hanwoo 50 k SNP panel for four carcass traits in commercial Hanwoo population by adding highly predictive variants from sequence data. A total of 16,892 Hanwoo cattle with phenotypes (i.e., backfat thickness, carcass weight, longissimus muscle area, and marbling score), 50 k genotypes, and WGS imputed genotypes were used. We partitioned imputed WGS data according to functional annotation [intergenic (IGR), intron (ITR), regulatory (REG), synonymous (SYN), and non-synonymous (NSY)] to characterize the genomic regions that will deliver higher predictive power for the traits investigated. Animals were assigned into two groups, the discovery set (7324 animals) used for predictive variant detection and the cross-validation set for genomic prediction. Genome-wide association studies were performed by trait to every genomic region and entire WGS data for the pre-selection of variants. Each set of pre-selected SNPs with different density (1000, 3000, 5000, or 10,000) were added to the 50 k genotypes separately and the predictive performance of each set of genotypes was assessed using the genomic best linear unbiased prediction (GBLUP). Results showed that the predictive performance of the customized Hanwoo 50 k SNP panel can be improved by the addition of pre-selected variants from the WGS data, particularly 3000 variants from each trait, which is then sufficient to improve the prediction accuracy for all traits. When 12,000 pre-selected variants (3000 variants from each trait) were added to the 50 k genotypes, the prediction accuracies increased by 9.9, 9.2, 6.4, and 4.7% for backfat thickness, carcass weight, longissimus muscle area, and marbling score compared to the regular 50 k SNP panel, respectively. In terms of prediction bias, regression coefficients for all sets of genotypes in all traits were close to 1, indicating an unbiased prediction. The strategy used to select variants based on functional annotation did not show a clear advantage compared to using whole-genome. Nonetheless, such pre-selected SNPs from the IGR region gave the highest improvement in prediction accuracy among genomic regions and the values were close to those obtained using the WGS data for all traits. We concluded that additional gain in prediction accuracy when using pre-selected variants appears to be trait-dependent, and using WGS data remained more accurate compared to using a specific genomic region.

Simple SummaryDue to the extensive marbling of its beef, Hanwoo (Korean native cattle) has continuously gained popularity and has become a mainstay in South Korea’s animal industry. In any beef cattle production system, reproductive performance is one of the main economic aspects taken into consideration. Therefore, genetic parameter estimates are necessary to obtain indices in order to maximize the response to selection, which in turn could lead to higher profitability. To date, knowledge on the genetic parameters for reproductive traits in Hanwoo cattle is still limited. Therefore, this study estimated the variance components, heritability, phenotypic, and genetic correlations of age at first calving (AFC), calving interval (CI), days open (DO), and gestation length (GL) of Hanwoo cattle. This was done using single-trait and multi-trait animal models. Results revealed the low heritability estimates for AFC, CI, DO, and GL in both single-trait and multi-trait models, which indicated the probable slow response of these traits due to direct selection. Moreover, phenotypic and genetic correlations varied from low to high among the reproductive traits of interest. Nevertheless, heritability estimates and genetic correlations shown in this study will prove to be vital as initial estimates are considered in the genetic improvement program of Hanwoo cattle.AbstractGenetic parameters for the reproductive traits of Hanwoo cattle were estimated using data obtained from 15,355 cows in 92 herds across South Korea, which were inseminated from May 1997 to July 2016. An “average information” restricted maximum likelihood (REML) procedure that fit in single-trait and multi-trait animal models was used to estimate the variance components of age at first calving (AFC), calving interval (CI), days open (DO), and gestation length (GL). Results showed the low estimates of heritability for all reproductive traits from both single-trait and multi-trait models. Estimates of heritability for AFC were 0.08 and 0.10 with single-trait and multi-trait models, respectively, while the estimates of heritability using the same animal models ranged from 0.01 to 0.07, 0.01 to 0.09, and 0.10 to 0.16 for CI, DO, and GL, accordingly. While AFC showed positive genetic correlations of 0.52 and 0.46 with CI and DO, respectively, the estimates of genetic and phenotypic correlations of GL with AFC were close to zero. Moreover, phenotypic correlations of GL with CI and DO were also close to zero; however, the corresponding genetic correlations were 0.13 and –0.38 for CI and DO, respectively. These estimated variance components and genetic correlations for reproductive traits can be utilized for genetic improvement programs of Hanwoo cattle.

The genomic best linear unbiased prediction (GBLUP) method has been widely used in routine genomic evaluation as it assumes a common variance for all single nucleotide polymorphism (SNP). However, this is unlikely in the case of traits influenced by major SNP. Hence, the present study aimed to improve the accuracy of GBLUP by using the weighted GBLUP (WGBLUP), which gives more weight to important markers for various carcass traits of Hanwoo cattle, such as backfat thickness (BFT), carcass weight (CWT), eye muscle area (EMA), and marbling score (MS). Linear and different nonlinearA SNP weighting procedures under WGBLUP were evaluated and compared with unweighted GBLUP and traditional pedigree-based methods (PBLUP). WGBLUP methods were assessed over ten iterations. Phenotypic data from 10,215 animals from different commercial herds that were slaughtered at approximately 30-month-old of age were used. All these animals were genotyped using Illumina Bovine 50k SNP chip and were divided into a training and a validation population by birth date on 1 November 2015. Genomic prediction accuracies obtained in the nonlinearA weighting methods were higher than those of the linear weighting for all traits. Moreover, unlike with linear methods, no sudden drops in the accuracy were noted after the peak was reached in nonlinearA methods. The average accuracies using PBLUP were 0.37, 0.49, 0.40, and 0.37, and 0.62, 0.74, 0.67, and 0.65 using GBLUP for BFT, CWT, EMA, and MS, respectively. Moreover, these accuracies of genomic prediction were further increased to 4.84% and 2.70% for BFT and CWT, respectively by using the nonlinearA method under the WGBLUP model. For EMA and MS, WGBLUP was as accurate as GBLUP. Our results indicate that the WGBLUP using a nonlinearA weighting method provides improved predictions for CWT and BFT, suggesting that the ability of WGBLUP over the other models by weighting selected SNPs appears to be trait-dependent.

Hanwoo breed is preferred in South Korea because of the high standards in marbling and the palatability of its meat. Numerous studies have been conducted and are ongoing to increase the meat production and quality in this beef population. The aim of this study was to estimate and compare genetic parameters for carcass traits using BLUPF90 software. Four models were constructed, single trait pedigree model (STPM), single-trait genomic model (STGM), multi-trait pedigree model (MTPM), and multi-trait genomic model (MTGM), using the pedigree, phenotype, and genomic information of 7991 Hanwoo cattle. Four carcass traits were evaluated: Back fat thickness (BFT), carcass weight (CWT), eye muscle area (EMA), and marbling score (MS). Heritability estimates of 0.40 and 0.41 for BFT, 0.33 and 0.34 for CWT, 0.36 and 0.37 for EMA, and 0.35 and 0.38 for MS were obtained for the single-trait pedigree model and the multi-trait pedigree model, respectively, in Hanwoo. Further, the genomic model showed more improved results compared to the pedigree model, with heritability of 0.39 (CWT), 0.39 (EMA), and 0.46 (MS), except for 0.39 (BFT), which may be due to random events. Utilization of genomic information in the form of single nucleotide polymorphisms (SNPs) has allowed more capturing of the variance from the traits improving the variance components.

Genetic parameters and trends in the average daily gain (ADG), backfat thickness (BF), loin muscle area (LMA), lean percentage (LP), and age at 90 kg (D90) were estimated for populations of Landrace and Yorkshire pigs. Additionally, the correlations between these production traits and litter traits were estimated. Litter traits included total born (TB) and number born alive (NBA). The data used for this study were obtained from eight farms during 1999 to 2016. Analyses were carried out with a multivariate animal model to estimate genetic parameters for production traits while bivariate analyses were performed to estimate the correlations between production and litter traits. The heritability estimates were 0.52 and 0.43 for ADG; 0.54 and 0.45 for BF; 0.25 and 0.26 for LMA; 0.54 and 0.48 for LP; and 0.56 and 0.46 for D90 in the Landrace and Yorkshire breeds, respectively. The ADG and D90 showed low genetic correlation with BF and LP. The LMA had -0.40, -0.32, 0.49, and 0.39 genetic correlations with ADG, BF, LP, and D90, respectively. Genetic correlations between production and litter traits were generally low, except for the correlations between LMA and TB (-0.23) in Landrace and ADG and TB (-0.16), ADG and NBA (-0.18), D90 and TB (0.19), and D90 and NBA (0.20) in Yorkshire. Genetic trends in production traits were all favorable except for LMA.

ObjectiveThis study was conducted to estimate breed-specific variance components for total number born (TNB), number born alive (NBA) and mortality rate from birth through weaning including stillbirths (MORT) of three main swine breeds in Korea. In addition, the importance of including maternal genetic and service sire effects in estimation models was evaluated.MethodsRecords of farrowing traits from 6,412 Duroc, 18,020 Landrace, and 54,254 Yorkshire sows collected from January 2001 to September 2016 from different farms in Korea were used in the analysis. Animal models and the restricted maximum likelihood method were used to estimate variances in animal genetic, permanent environmental, maternal genetic, service sire and residuals.ResultsThe heritability estimates ranged from 0.072 to 0.102, 0.090 to 0.099, and 0.109 to 0.121 for TNB; 0.087 to 0.110, 0.088 to 0.100, and 0.099 to 0.107 for NBA; and 0.027 to 0.031, 0.050 to 0.053, and 0.073 to 0.081 for MORT in the Duroc, Landrace and Yorkshire breeds, respectively. The proportion of the total variation due to permanent environmental effects, maternal genetic effects, and service sire effects ranged from 0.042 to 0.088, 0.001 to 0.031, and 0.001 to 0.021, respectively. Spearman rank correlations among models ranged from 0.98 to 0.99, demonstrating that the maternal genetic and service sire effects have small effects on the precision of the breeding value.ConclusionModels that include additive genetic and permanent environmental effects are suitable for farrowing traits in Duroc, Landrace, and Yorkshire populations in Korea. This breed-specific variance components estimates for litter traits can be utilized for pig improvement programs in Korea.