JARB Journal of Animal Reproduction and Biotehnology

OPEN ACCESS pISSN: 2671-4639
eISSN: 2671-4663

Article Search

Original Article

Article Original Article
Split Viewer

Journal of Animal Reproduction and Biotechnology 2023; 38(4): 268-274

Published online December 31, 2023

https://doi.org/10.12750/JARB.38.4.268

Copyright © The Korean Society of Animal Reproduction and Biotechnology.

Genetic evaluation for economic traits of commercial Hanwoo population using single-step GBLUP

Gwang Hyeon Lee1 , Khaliunaa Tseveen1 , Yoon Seok Lee1,2,3 and Hong Sik Kong1,2,3,*

1Department of Applied Biotechnology, The Graduate School of Hankyong National University, Anseong 17579, Korea
2Genomic Informatics Center, Hankyong National University, Anseong 17579, Korea
3Gyeonggi Regional Research Center, Hankyong National University, Anseong 17579, Korea

Correspondence to: Hong Sik Kong
E-mail: kebinkhs@hknu.ac.kr

Received: December 12, 2023; Revised: December 18, 2023; Accepted: December 19, 2023

This is an open-access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/), which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Background: Recently, the single-step genomic best linear unbiased prediction (ssGBLUP) method, which incorporates not only genomic information but also phenotypic information of pedigree, is under study. In this study, we performed a ssGBLUP analysis on a commercial Hanwoo population using phenotypic, genotypic, and pedigree data.
Methods: The test population comprised Hanwoo 1,740 heads raised in four regions of Korea, while the reference population used Hanwoo 18,499 heads raised across the country and two-generation pedigree data. Analysis was performed using genotype data generated by the Hanwoo 50 K SNP beadchip.
Results: The mean Genome estimated breeding values (GEBVs) estimated using the ssGBLUP methods for carcass weight (CWT), eye muscle area (EMA), back fat thickness (BFT), and marbling score (MS) were 7.348, 1.515, -0.355, and 0.040, respectively, while the accuracy of each trait was 0.749, 0.733, 0.769, and 0.768, respectively. When the correlation analysis between the GEBVs as a result of this study and the actual slaughter performance was confirmed, CWT, EMA, BFT, and MS were reported to be 0.519, 0.435, 0.444, and 0.543, respectively.
Conclusions: Our results suggest that the ssGBLUP method enables a more accurate evaluation because it conducts a genetic evaluation of an individual using not only genotype information but also phenotypic information of the pedigree. Individual evaluation using the ssGBLUP method is considered effective for enhancing the genetic ability of farms and enabling accurate and rapid improvements. It is considered that if more pedigree information of reference population is collected for analysis, genetic ability can be evaluated more accurately.

Keywords: accuracy, carcass traits, genetic evaluation, Hanwoo, ssGBLUP

Hanwoo (also known as Korean Native Cattle; Bos taurus coronae) is a native breed of beef cattle in South Korea, well known (highly sought) for its excellent meat quality. Excessive marbling has emerged as a significant meat product quality determinant in commercial beef markets (Choy et al., 2012; Choi et al., 2015). Research on Hanwoo breeds has been conducted for many years and has expanded in recent years. In breeding, selection is commonly based on the breeding value calculated using the best linear unbiased prediction (BLUP) approach, which utilizes an individual’s phenotypic and pedigree data. BLUP is arguably the most frequently used selection technique for determining estimated breeding values (EBVs) (Henderson, 1984; Meuwissen et al., 2001). Genomic selection (GS) is a widely recognized tool for genetic improvement. GS implementation has evolved rapidly in various livestock and plant species (Meuwissen et al., 2001; VanRaden, 2012; Misztal et al., 2020). The statistical method for GS is the genomic best linear unbiased prediction (GBLUP) method, which is based on the infinitesimal model, and the breeding value is evaluated using the genomic relationship matrix (GRM) rather than the Numerator Relationship Matrix (Garrick, 2007; VanRaden, 2007; Zhang, 2007). The GBLUP method calculates single nucleotide polymorphism (SNP) effects based on individual SNP and phenotypic information using SNP data evenly distributed over the whole genome using SNP chips (Meuwissen et al., 2001). The breeding value of individuals lacking phenotypic data can be estimated using their SNP information and pre-calculated SNP effect values. Additionally, even with an identical pedigree, the breeding value of an individual can be estimated precisely because of the ability of the GBLUP method to create a GRM.

The single-step genomic best linear unbiased prediction (ssGBLUP) integrates all phenotypic, pedigree, and genomic information (genotypic data) in a single-step analysis, which is available simultaneously to predict genomic merit values for genotyped and non-genotyped animals through a combined matrix (Gao et al., 2019; Song et al., 2019). Furthermore, it enables more precise and effective estimates of genetic parameters. Numerous genomic selection studies in beef cattle (Gordo et al., 2016; Lee et al., 2017), dairy cattle (VanRaden, 2020), pigs (Forni et al., 2011; Xiang et al., 2016), and chickens (Chen et al., 2011) have demonstrated the advantages of ssGBLUP in various domesticated animals. Therefore, this study investigated the genomic estimated breeding value (GEBV) of a commercial Hanwoo population using the ssGBLUP method and phenotypic, genotypic, and pedigree information. Additionally, the accuracy of the estimated GEBV was determined by analyzing its accuracy.

Test population

We performed a ssGBLUP analysis on a commercial Hanwoo population using phenotypic, genotypic, and pedigree data (Table 1). The analyzed population comprised a reference population and a test population. To determine the genotype of the test population, DNA was collected from 1,740 tail hair root samples from 1,221 Hanwoo cows and 519 steers raised in the Gyeonggi, Gangwon, Chungbuk, and Gyeongbuk provinces. The extracted DNA was analyzed by determining the genotype using a Hanwoo 50 K SNP BeadChip (Illumina, San Diego, CA, USA). For the test population, information on 1,740 heads was collected from three generations of pedigree.

Table 1 . Total number of references and test population

SexTest populationReference population

PedigreeGenotype
Cow1,22136,327355
Steer51938,40318,144
Total1,74074,73018,499


Reference population

The reference population for the analysis used genotype information from 18,499 cattle raised and transported across the country. Pedigree information was collected from the Korean Animal Improvement Association and used for analysis. For the reference population, information from two generations of pedigree on 18,499 heads was collected, and a total of 74,730 heads information were analyzed. The Livestock Product Quality Evaluation Institute provided phenotypic information on 62,567 of the 74,730 heads, including carcass weight (CWT), eye muscle area (EMA), back fat thickness (BFT), and marbling score (MS).

Genomic information quality control (QC)

To increase the accuracy of the analysis, QC was performed on the SNP information, and unnecessary SNPs were removed. Genotype information was generated from raw data using the GenomeStudio 2.0 software (Illumina Inc., 2016). Plink version 1.9 (Purcell et al., 2007; Chang et al., 2015) was used for QC after the final report file was converted into an input file format. The QC criteria were minor allele frequency (MAF) < 1% (4,721), SNP call rate < 90% (730), and Hardy–Weinberg equilibrium p-value < 1 × 10-7 (1,898). After removing 6,647 SNPs that did not meet the QC criteria, The final dataset included 45,548 SNPs from 52,195 SNPs.

Statistical analysis

The program used for object evaluation was the BLUPF90 family program (Misztal et al., 2014). Genetic parameters were estimated to assess GEBV. Genetic parameters were estimated using the AIREMLF90 program (Misztal et al., 2014), with year and month of birth, year and month of slaughter, and sex as fixed effects. The mixed linear model used to estimate the dielectric parameters of the additive genetic effect is as follows (Henderson, 1975):

YP=Xβ+Zμ+e

where YP is the actual observed values of the four major economic traits (CWT, EMA, BFT, and MS), where Z is the random effect vector of individuals, β is a vector of estimates for unknown fixed effects (year and month of birth, year and month of slaughter, and sex), μ is the additive genetic effect, e is the random environmental effect (Henderson and Quaas, 1976). The breeding value estimated using the ssGBLUP method was BLUPF90 (Misztal et al., 2014), and H-1 was generated using PREGSF90 (Misztal et al., 2014). The analysis was performed single trait analysis and the estimation of H-1 was performed using the following formula (Misztal et al., 2014):

H1= A1+ 000αG 1 +βA22 1

H is a relationship matrix constructed from the pedigree and genotype information. A is an NRM based on pedigree, whereas G is a GRM based on genotype. α and β are the correction values for G and A, respectively. α (0.95) and β (0.05) were analyzed using default parameter settings.

The heritability was calculated using the following formula using the estimated genetic parameters.

h2 = σα2σP2              σp2= σα2+ σe2

Where h2 is heritability for each trait, σα2 is the represents genetic variance, σe2 is the represents residual variance and σp2 is phenotypic variance. Heritability can be calculated as the percentage of genetic variance to phenotypic variance, which is the sum of genetic and residual variances.

Accuracy analysis

The accuracy of the GEBV estimated using the ssGLBUP method was estimated. Accuracy was calculated using the standard error of the breeding value estimated for each individual among the solution values generated using the BLUPF90 program (Misztal et al., 2014) and the additive genetic variance estimated using the AIREMLF90 program (Misztal et al., 2014) analysis. Accuracy was calculated as follows:

Accuracy = 1PEVσα2

where accuracy is the estimated breeding value accuracy, PEV is the standard error of the estimated breeding value for each individual, and σα2 is the estimated genetic variance. Accuracy was analyzed by determining the correlation coefficient between the observed phenotype of the test population and GEBV using the ssGBLUP method, and the proximity between carcass characteristics and GEBV was confirmed.

Basic statistical analysis of the reference population

To estimate the GEBV using the ssGBLUP method, phenotypic data from 62,567 heads of the Hanwoo population slaughtered throughout the country were used for analysis, and genotype data from 18,499 heads of the Hanwoo population slaughtered throughout the country served as reference population genotype data. Table 2 presents the basic statistics. The basic statistics for the entire pedigree information were 349.37 ± 54.13 kg, 83.86 ± 12.67 cm2, 13.80 ± 6.09 mm, and 4.01 ± 1.93 score for Hanwoo cows CWT, EMA, BFT, and MS, respectively. However, for steers, the statistics were 442.02 ± 55.20 kg, 95.15 ± 12.33 cm2, 13.80 ± 5.10 mm, and 5.82 ± 1.95 score, respectively. Kim et al. (2022) compiled the genomic and phenotypic information of Hanwoo cows and steers shipped across the country as the reference population for ssGBLUP analysis after researching the basic statistics of the reference population analyzed in previous studies. The basic statistics including cows without genomic information 352.92 ± 46.77 kg, 83.73 ± 11.49 cm2, 13.32 ± 5.37 mm, and 4.36 ± 1.9 score, respectively. Furthermore, for steers without genomic information, the basic statistics were 443.51 ± 49.80 kg, 93.72 ± 11.53 cm2, 13.62 ± 5.05 mm, and 5.94 ± 1.84 score, respectively.

Table 2 . Basic statistics of reference population according to the ssGBLUP method

SubsetSexNo.TraitsMeanS.D.MaxMin
Pedigree + genotypeCow27,264CWT (kg)349.3754.13667110
EMA (cm2)83.8612.6714722
BFT (mm)13.806.09771
MS (score)4.011.9391
Steer35,303CWT (kg)442.0255.20745148
EMA (cm2)95.1512.3317120
BFT (mm)13.805.10511
MS (score)5.821.9591
GenotypeCow355CWT (kg)356.1647.94498159
EMA (cm2)87.8712.1711938
BFT (mm)13.895.03374
MS (score)4.951.9991
Steer18,144CWT (kg)444.0149.57692160
EMA (cm2)96.7012.2516020
BFT (mm)14.274.93471
MS (score)5.981.8791

CWT, carcass weight; EMA, eye muscle area; BFT, back fat thickness; MS, marbling score; SD, standard deviation.



Park et al. (2020) conducted a genetic evaluation on Korean Hanwoo steer information raised by the National Institute of Animal Science. They included Hanwoo steers without genomic information and computed basic statistics for CWT, EMA, BFT, and MS. The reported scores were 371.63 ± 42.07 kg, 81.65 ± 8.82 cm2, 9.71 ± 3.82 mm, and 3.58 ± 1.62 score, respectively. When compared with the results of this study, the reference livestock basic statistics value of Kim et al. (2022) showed similar values to the study’s results, although it was higher than the findings of Park et al. (2020). Park et al. (2020) used nationally managed steering data from 1997 to 2018, which was thought to be the result of including relatively young individuals.

To test the normality of the complete reference population, the findings of quantile-quantile plots (Q-Q plots) for each trait (CWT, EMA, BFT, and MS) are shown in Fig. 1. As a result of the Q-Q plots, outlier values that strayed from the line were identified, but they were included in the analysis because they represented the diversity of the individual.

Figure 1. Quantile-quantile plot (Q-Q plot) for the four carcass traits in the ssGBLUP reference population.

Estimation of genetic parameter

Genetic parameters were estimated by the ssGBLUP method. Table 3 presents the estimation results of the genetic parameters for genetic variance, residual variance, and heritability. The genetic variance was estimated to be 980.84, 50.594, 15.072, and 1.609 for CWT, EMA, BFT, and MS, respectively. The residual variances were estimated to be 1530.4, 97.642, 17.660, and 1.894, and the heritability values using genetic and residual variances were estimated to be 0.391, 0.341, 0.460, and 0.459, respectively. After reviewing previous studies, Kim et al. (2022) estimated the ssGBLUP accuracy based on the amount of genomic information available in the pedigree. They reported heritability values of 0.429, 0.372, 0.357, and 0.471 for the CWT, EMA, BFT, and MS treatments, respectively. Lopez et al. (2020) found heritability values of 0.38, 0.34, 0.38, and 0.48 in an imputation accuracy study using the Bovine HD BeadChip (777 k) and Hanwoo 50 K SNP beadchip. Compared to the results of this study, the analysis included the phenotype of cows, which is thought to be the reason why the BFT heritability was estimated to be low. According to Park et al. (2020), heritability estimates for CWT, EMA, BFT, and MS were 0.35, 0.44, 0.46, and 0.56, respectively. In contrast to the results of this study, the heritability of BFT and CWT was low because the population used in the previous study differed from the population used in this study by using nationally managed individuals.

Table 3 . Genetic parameters estimated using the ssGBLUP method based on genotype and pedigree information

Analysis methodTypeTrait

CWTEMABFTMS
ssGBLUP methodGenetic
Variance
980.8450.59415.0721.609
Residual
Variance
1,530.497.64217.6601.894
Heritability0.3910.3410.4600.459

ssGBLUP, single-step genomic best linear unbiased prediction; CWT, carcass weight; EMA, eye muscle area; BFT, back fat thickness; MS, marbling score.



Accuracy analysis of breeding value

Table 4 lists the GEBV, GEBV accuracy, and realized accuracy of the subjects, estimated using the ssGBLUP method. The mean GEBV estimated using the ssGBLUP methods for CWT, EMA, BFT, and MS were 7.348, 1.515, -0.355, and 0.040, respectively and the accuracy of each trait was 0.749, 0.733, 0.769, and 0.768, respectively. The accuracy of GEBV was highest for BFT and lowest for EMA. Because the breeding value accuracy was affected by the estimated heritability, it showed the same pattern as that of the previous heritability. Kim et al. (2022) conducted a comparative study of the accuracy of ssGBLUP based on the amount of genomic information in the pedigree. They observed that the highest accuracies among the GEBV were 0.736, 0.718, 0.713, and 0.747 for CWT, EMA, BFT, and MS, respectively. Park et al. (2020) reported that for CWT, EMA, BFT, and MS, the accuracies of the estimated breeding values using the ssGBLUP method were 0.68, 0.69, 0.69, and 0.73, respectively. Lopez et al. (2020) reported that the breeding value accuracy estimated using the ssGBLUP method through single-trait analysis was 0.47, 0.42, 0.33, and 0.37 for CWT, EMA, BFT, and MS, respectively. Previous studies have reported a lower accuracy of the mean GEBV than the results obtained in this study. This difference is attributed to the larger size of the reference population used in this study compared to previous studies. Because the ssGBLUP method estimates GEBV using individual genotype information and pedigree phenotypic information, the larger the reference population, the more diverse the information can be explained.

Table 4 . Breeding value, accuracy, and realized accuracy were estimated using the ssGBLUP method

Analysis methodReference population sizeTypeTrait

CWT (kg)EMA (cm2)BFT (mm)MS (score)
ssGBLUP method18,499 + 74,730GEBV7.3481.515-0.3550.040
Standard
Error
20.6894.8242.4750.809
Accuracy0.7490.7330.7690.768
Realized accuracy0.5190.4350.4440.543

ssGBLUP, single-step genomic best linear unbiased prediction; CWT, carcass weight; EMA, eye muscle area; BFT, back fat thickness; MS, marbling score; GEBV, genomic estimated breeding value.



After confirming the correlation analysis (realized accuracy) between the GEBV and the actual slaughter performance, CWT, EMA, BFT, and MS were reported to be 0.519, 0.435, 0.444, and 0.543, respectively. The realized accuracy was the highest for MS and the lowest for EMA. In a previous study, Lee et al. (2023) demonstrated that the realized accuracy values in the genetic ability evaluation of a commercial Hanwoo population using the GBLUP method were 0.522, 0.404, 0.444, and 0.539 for CWT, EMA, BFT, and MS. Kim (2021) reported the realized accuracy of GEBV as 0.509, 0.446, 0.415, and 0.555, respectively, when using the ssGBLUP method based on the phenotypic information of 202,808 heads and the genotype information of 14,000 heads. When comparing the results of this study, the realized accuracy values demonstrated an increase compared with the GBLUP method. Collecting more phenotypic data related to pedigree will enable more accurate genetic evaluation.

The results of this study show that the ssGBLUP method allows for a more accurate evaluation of individuals because it incorporates into its analysis not only the genotype and phenotype information used in the current GBLUP method but also the pedigree and phenotype information of the reference population. Currently, the Korean Hanwoo industry is experiencing difficulties, such as rising feed prices and inflation, owing to changes in the international situation. To ensure a stable income, it is necessary to enhance the average genetic ability of Hanwoo cattle raised on farms to produce stable, high-grade Hanwoo cattle. Individual evaluation using the ssGBLUP method is considered effective for enhancing the genetic ability of farms and enabling accurate and rapid improvements. This study used the two generations of pedigree information on the reference population for analysis. Based on the results of this study, it is considered that obtaining further pedigree information of the reference population will enable high-accuracy individual selection.

Conceptualization, H.S.K.; methodology, G.H.L., K.T., Y.S.L., H.S.K.; investigation, G.H.L., and K.T.; data curation, G.H.L., K.T.; writing—original draft preparation, G.H.L., K.T.; writing—review and editing, G.H.L., K.T., Y.S.L., H.S.K.; supervision, H.S.K.; project administration, H.S.K.; funding acquisition, H.S.K.

  1. Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. 2015. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience 4:7.
    Pubmed KoreaMed CrossRef
  2. Chen CY, Misztal I, Aguilar I, Legarra A, Muir WM. 2011. Effect of different genomic relationship matrices on accuracy and scale. J. Anim. Sci. 89:2673-2679.
    Pubmed CrossRef
  3. Choi TJ, Alam M, Cho CI, Lee JG, Park B, Kim S, Koo Y, Roh SH. 2015. Genetic parameters for yearling weight, carcass traits, and primal-cut yields of Hanwoo cattle. J. Anim. Sci. 93:1511-1521.
    Pubmed CrossRef
  4. Choy YH, Park BH, Choi TJ, Choi JG, Cho KH, Lee SS, Choi YL, Koh KC, Kim HS. 2012. Estimation of relative economic weights of hanwoo carcass traits based on carcass market price. Asian-Australas. J. Anim. Sci. 25:1667-1673.
    Pubmed KoreaMed CrossRef
  5. Forni S, Aguilar I, Misztal I. 2011. Different genomic relationship matrices for single-step analysis using phenotypic, pedigree and genomic information. Genet. Sel. Evol. 43:1.
    Pubmed KoreaMed CrossRef
  6. Gao N, Teng J, Pan R, Li X, Ye S, Li J, Zhang H, Zhang X, Zhang Z. 2019. Accuracy of whole genome prediction with single-step GBLUP in a Chinese yellow-feathered chicken population. Livest. Sci. 230:103817.
    CrossRef
  7. Garrick DJ. 2007. Equivalent mixed model equations for genomic selection. J. Anim. Sci. 85:376.
  8. Gordo DG, Espigolan R, Tonussi RL, Júnior GA, Bresolin T, Magalhães AF, Feitosa FL, Baldi F, Carvalheiro R, Tonhati H, de Oliveira HN, Chardulo LA, de Albuquerque LG. 2016. Genetic parameter estimates for carcass traits and visual scores including or not genomic information. J. Anim. Sci. 94:1821-1826.
    Pubmed CrossRef
  9. Henderson CR. 1975. Best linear unbiased estimation and prediction under a selection model. Biometrics 31:423-447.
    Pubmed CrossRef
  10. Henderson CR. 1984. Applications of Linear Models in Animal Breeding, University of Guelph, Guelph.
  11. Henderson CR and Quaas RL. 1976. Multiple trait evaluation using relatives' records. J. Anim. Sci. 43:1188-1197.
    CrossRef
  12. Kim DH. 2021. Studies on the genetic evaluation of Hanwoo using the genomic information. [Doctoral dissertation, Jeonbuk National University]. RISS. http://www.riss.kr/link?id=T15777082
  13. Kim EH, Sun DW, Kang HC, Myung CH, Kim JY, Lee DH, Lee SH, Lim HT. 2022. Estimated of genomic estimated breeding value and accuracy analysis according to the amount of genotypes in the full-sib family. J. Agric. Life Sci. 56:171-178.
    CrossRef
  14. Kim SJ, Choi TJ, Son JI, Lee DM, Lee JJ, Lee JG, Lim HT, Koo YM. 2022. Estimation of genomic estimated breeding value(GEBV) and reliability for Hanwoo carcass traits using ssGBLUP. J. Anim. Breed. Genom. 6:57-72.
    CrossRef
  15. Lee GH, Lee YH, Kong HS. 2023. Genetic evaluation and accuracy analysis of commercial Hanwoo population using genomic data. J. Anim. Reprod. Biotechnol. 38:32-37.
    CrossRef
  16. Lee J, Cheng H, Garrick D, Golden B, Dekkers J, Park K, Lee D, Fernando R. 2017. Comparison of alternative approaches to single-trait genomic prediction using genotyped and non-genotyped Hanwoo beef cattle. Genet. Sel. Evol. 49:2.
    Pubmed KoreaMed CrossRef
  17. Lopez BI, Lee SH, Shin DH, Oh JD, Chai HH, Park W, Park JE, Lim D. 2020. Accuracy of genomic evaluation using imputed high-density genotypes for carcass traits in commercial Hanwoo population. Livest. Sci. 241:104256.
    CrossRef
  18. Meuwissen TH, Hayes BJ, Goddard ME. 2001. Prediction of total genetic value using genome-wide dense marker maps. Genetics 157:1819-1829.
    Pubmed KoreaMed CrossRef
  19. Misztal I, Lourenco D, Legarra A. 2020. Current status of genomic evaluation. J. Anim. Sci. 98:skaa101.
    Pubmed KoreaMed CrossRef
  20. Misztal I, Tsuruta S, Lourenco D, Masuda Y, Aguilar I, Legarra A, Vitezica Z. 2014. Manual for BLUPF90 Family of Programs, University of Georgia, Athens.
  21. Park MN, Alam M, Kim S, Park B, Lee SH, Lee SS. 2020. Genomic selection through single-step genomic best linear unbiased prediction improves the accuracy of evaluation in Hanwoo cattle. Asian-Australas. J. Anim. Sci. 33:1544-1557.
    Pubmed KoreaMed CrossRef
  22. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC. 2007. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81:559-575.
    Pubmed KoreaMed CrossRef
  23. Song H, Zhang J, Zhang Q, Ding X. 2019. Using different single-step strategies to improve the efficiency of genomic prediction on body measurement traits in pig. Front. Genet. 9:730.
    Pubmed KoreaMed CrossRef
  24. VanRaden PM. 2007. Genomic measures of relationship and inbreeding. Interbull Bull. 37:33-36.
  25. VanRaden PM. 2012. Avoiding bias from genomic pre-selection in converting daughter information across countries. Interbull Bull. 45.
  26. VanRaden PM. 2020. Symposium review: how to implement genomic selection. J. Dairy Sci. 103:5291-5301.
    Pubmed CrossRef
  27. Xiang T, Nielsen B, Su G, Legarra A, Christensen OF. 2016. Application of single-step genomic evaluation for crossbred performance in pig. J. Anim. Sci. 94:936-948.
    Pubmed CrossRef
  28. Zhang Z, Todhunter RJ, Buckler ES, Van Vleck LD. 2007. Technical note: use of marker-based relationships with multiple-trait derivative-free restricted maximal likelihood. J. Anim. Sci. 85:881-885.
    Pubmed CrossRef

Article

Original Article

Journal of Animal Reproduction and Biotechnology 2023; 38(4): 268-274

Published online December 31, 2023 https://doi.org/10.12750/JARB.38.4.268

Copyright © The Korean Society of Animal Reproduction and Biotechnology.

Genetic evaluation for economic traits of commercial Hanwoo population using single-step GBLUP

Gwang Hyeon Lee1 , Khaliunaa Tseveen1 , Yoon Seok Lee1,2,3 and Hong Sik Kong1,2,3,*

1Department of Applied Biotechnology, The Graduate School of Hankyong National University, Anseong 17579, Korea
2Genomic Informatics Center, Hankyong National University, Anseong 17579, Korea
3Gyeonggi Regional Research Center, Hankyong National University, Anseong 17579, Korea

Correspondence to:Hong Sik Kong
E-mail: kebinkhs@hknu.ac.kr

Received: December 12, 2023; Revised: December 18, 2023; Accepted: December 19, 2023

This is an open-access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/), which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background: Recently, the single-step genomic best linear unbiased prediction (ssGBLUP) method, which incorporates not only genomic information but also phenotypic information of pedigree, is under study. In this study, we performed a ssGBLUP analysis on a commercial Hanwoo population using phenotypic, genotypic, and pedigree data.
Methods: The test population comprised Hanwoo 1,740 heads raised in four regions of Korea, while the reference population used Hanwoo 18,499 heads raised across the country and two-generation pedigree data. Analysis was performed using genotype data generated by the Hanwoo 50 K SNP beadchip.
Results: The mean Genome estimated breeding values (GEBVs) estimated using the ssGBLUP methods for carcass weight (CWT), eye muscle area (EMA), back fat thickness (BFT), and marbling score (MS) were 7.348, 1.515, -0.355, and 0.040, respectively, while the accuracy of each trait was 0.749, 0.733, 0.769, and 0.768, respectively. When the correlation analysis between the GEBVs as a result of this study and the actual slaughter performance was confirmed, CWT, EMA, BFT, and MS were reported to be 0.519, 0.435, 0.444, and 0.543, respectively.
Conclusions: Our results suggest that the ssGBLUP method enables a more accurate evaluation because it conducts a genetic evaluation of an individual using not only genotype information but also phenotypic information of the pedigree. Individual evaluation using the ssGBLUP method is considered effective for enhancing the genetic ability of farms and enabling accurate and rapid improvements. It is considered that if more pedigree information of reference population is collected for analysis, genetic ability can be evaluated more accurately.

Keywords: accuracy, carcass traits, genetic evaluation, Hanwoo, ssGBLUP

INTRODUCTION

Hanwoo (also known as Korean Native Cattle; Bos taurus coronae) is a native breed of beef cattle in South Korea, well known (highly sought) for its excellent meat quality. Excessive marbling has emerged as a significant meat product quality determinant in commercial beef markets (Choy et al., 2012; Choi et al., 2015). Research on Hanwoo breeds has been conducted for many years and has expanded in recent years. In breeding, selection is commonly based on the breeding value calculated using the best linear unbiased prediction (BLUP) approach, which utilizes an individual’s phenotypic and pedigree data. BLUP is arguably the most frequently used selection technique for determining estimated breeding values (EBVs) (Henderson, 1984; Meuwissen et al., 2001). Genomic selection (GS) is a widely recognized tool for genetic improvement. GS implementation has evolved rapidly in various livestock and plant species (Meuwissen et al., 2001; VanRaden, 2012; Misztal et al., 2020). The statistical method for GS is the genomic best linear unbiased prediction (GBLUP) method, which is based on the infinitesimal model, and the breeding value is evaluated using the genomic relationship matrix (GRM) rather than the Numerator Relationship Matrix (Garrick, 2007; VanRaden, 2007; Zhang, 2007). The GBLUP method calculates single nucleotide polymorphism (SNP) effects based on individual SNP and phenotypic information using SNP data evenly distributed over the whole genome using SNP chips (Meuwissen et al., 2001). The breeding value of individuals lacking phenotypic data can be estimated using their SNP information and pre-calculated SNP effect values. Additionally, even with an identical pedigree, the breeding value of an individual can be estimated precisely because of the ability of the GBLUP method to create a GRM.

The single-step genomic best linear unbiased prediction (ssGBLUP) integrates all phenotypic, pedigree, and genomic information (genotypic data) in a single-step analysis, which is available simultaneously to predict genomic merit values for genotyped and non-genotyped animals through a combined matrix (Gao et al., 2019; Song et al., 2019). Furthermore, it enables more precise and effective estimates of genetic parameters. Numerous genomic selection studies in beef cattle (Gordo et al., 2016; Lee et al., 2017), dairy cattle (VanRaden, 2020), pigs (Forni et al., 2011; Xiang et al., 2016), and chickens (Chen et al., 2011) have demonstrated the advantages of ssGBLUP in various domesticated animals. Therefore, this study investigated the genomic estimated breeding value (GEBV) of a commercial Hanwoo population using the ssGBLUP method and phenotypic, genotypic, and pedigree information. Additionally, the accuracy of the estimated GEBV was determined by analyzing its accuracy.

MATERIALS AND METHODS

Test population

We performed a ssGBLUP analysis on a commercial Hanwoo population using phenotypic, genotypic, and pedigree data (Table 1). The analyzed population comprised a reference population and a test population. To determine the genotype of the test population, DNA was collected from 1,740 tail hair root samples from 1,221 Hanwoo cows and 519 steers raised in the Gyeonggi, Gangwon, Chungbuk, and Gyeongbuk provinces. The extracted DNA was analyzed by determining the genotype using a Hanwoo 50 K SNP BeadChip (Illumina, San Diego, CA, USA). For the test population, information on 1,740 heads was collected from three generations of pedigree.

Table 1. Total number of references and test population.

SexTest populationReference population

PedigreeGenotype
Cow1,22136,327355
Steer51938,40318,144
Total1,74074,73018,499


Reference population

The reference population for the analysis used genotype information from 18,499 cattle raised and transported across the country. Pedigree information was collected from the Korean Animal Improvement Association and used for analysis. For the reference population, information from two generations of pedigree on 18,499 heads was collected, and a total of 74,730 heads information were analyzed. The Livestock Product Quality Evaluation Institute provided phenotypic information on 62,567 of the 74,730 heads, including carcass weight (CWT), eye muscle area (EMA), back fat thickness (BFT), and marbling score (MS).

Genomic information quality control (QC)

To increase the accuracy of the analysis, QC was performed on the SNP information, and unnecessary SNPs were removed. Genotype information was generated from raw data using the GenomeStudio 2.0 software (Illumina Inc., 2016). Plink version 1.9 (Purcell et al., 2007; Chang et al., 2015) was used for QC after the final report file was converted into an input file format. The QC criteria were minor allele frequency (MAF) < 1% (4,721), SNP call rate < 90% (730), and Hardy–Weinberg equilibrium p-value < 1 × 10-7 (1,898). After removing 6,647 SNPs that did not meet the QC criteria, The final dataset included 45,548 SNPs from 52,195 SNPs.

Statistical analysis

The program used for object evaluation was the BLUPF90 family program (Misztal et al., 2014). Genetic parameters were estimated to assess GEBV. Genetic parameters were estimated using the AIREMLF90 program (Misztal et al., 2014), with year and month of birth, year and month of slaughter, and sex as fixed effects. The mixed linear model used to estimate the dielectric parameters of the additive genetic effect is as follows (Henderson, 1975):

YP=Xβ+Zμ+e

where YP is the actual observed values of the four major economic traits (CWT, EMA, BFT, and MS), where Z is the random effect vector of individuals, β is a vector of estimates for unknown fixed effects (year and month of birth, year and month of slaughter, and sex), μ is the additive genetic effect, e is the random environmental effect (Henderson and Quaas, 1976). The breeding value estimated using the ssGBLUP method was BLUPF90 (Misztal et al., 2014), and H-1 was generated using PREGSF90 (Misztal et al., 2014). The analysis was performed single trait analysis and the estimation of H-1 was performed using the following formula (Misztal et al., 2014):

H1= A1+ 000αG 1 +βA22 1

H is a relationship matrix constructed from the pedigree and genotype information. A is an NRM based on pedigree, whereas G is a GRM based on genotype. α and β are the correction values for G and A, respectively. α (0.95) and β (0.05) were analyzed using default parameter settings.

The heritability was calculated using the following formula using the estimated genetic parameters.

h2 = σα2σP2              σp2= σα2+ σe2

Where h2 is heritability for each trait, σα2 is the represents genetic variance, σe2 is the represents residual variance and σp2 is phenotypic variance. Heritability can be calculated as the percentage of genetic variance to phenotypic variance, which is the sum of genetic and residual variances.

Accuracy analysis

The accuracy of the GEBV estimated using the ssGLBUP method was estimated. Accuracy was calculated using the standard error of the breeding value estimated for each individual among the solution values generated using the BLUPF90 program (Misztal et al., 2014) and the additive genetic variance estimated using the AIREMLF90 program (Misztal et al., 2014) analysis. Accuracy was calculated as follows:

Accuracy = 1PEVσα2

where accuracy is the estimated breeding value accuracy, PEV is the standard error of the estimated breeding value for each individual, and σα2 is the estimated genetic variance. Accuracy was analyzed by determining the correlation coefficient between the observed phenotype of the test population and GEBV using the ssGBLUP method, and the proximity between carcass characteristics and GEBV was confirmed.

RESULTS

Basic statistical analysis of the reference population

To estimate the GEBV using the ssGBLUP method, phenotypic data from 62,567 heads of the Hanwoo population slaughtered throughout the country were used for analysis, and genotype data from 18,499 heads of the Hanwoo population slaughtered throughout the country served as reference population genotype data. Table 2 presents the basic statistics. The basic statistics for the entire pedigree information were 349.37 ± 54.13 kg, 83.86 ± 12.67 cm2, 13.80 ± 6.09 mm, and 4.01 ± 1.93 score for Hanwoo cows CWT, EMA, BFT, and MS, respectively. However, for steers, the statistics were 442.02 ± 55.20 kg, 95.15 ± 12.33 cm2, 13.80 ± 5.10 mm, and 5.82 ± 1.95 score, respectively. Kim et al. (2022) compiled the genomic and phenotypic information of Hanwoo cows and steers shipped across the country as the reference population for ssGBLUP analysis after researching the basic statistics of the reference population analyzed in previous studies. The basic statistics including cows without genomic information 352.92 ± 46.77 kg, 83.73 ± 11.49 cm2, 13.32 ± 5.37 mm, and 4.36 ± 1.9 score, respectively. Furthermore, for steers without genomic information, the basic statistics were 443.51 ± 49.80 kg, 93.72 ± 11.53 cm2, 13.62 ± 5.05 mm, and 5.94 ± 1.84 score, respectively.

Table 2. Basic statistics of reference population according to the ssGBLUP method.

SubsetSexNo.TraitsMeanS.D.MaxMin
Pedigree + genotypeCow27,264CWT (kg)349.3754.13667110
EMA (cm2)83.8612.6714722
BFT (mm)13.806.09771
MS (score)4.011.9391
Steer35,303CWT (kg)442.0255.20745148
EMA (cm2)95.1512.3317120
BFT (mm)13.805.10511
MS (score)5.821.9591
GenotypeCow355CWT (kg)356.1647.94498159
EMA (cm2)87.8712.1711938
BFT (mm)13.895.03374
MS (score)4.951.9991
Steer18,144CWT (kg)444.0149.57692160
EMA (cm2)96.7012.2516020
BFT (mm)14.274.93471
MS (score)5.981.8791

CWT, carcass weight; EMA, eye muscle area; BFT, back fat thickness; MS, marbling score; SD, standard deviation..



Park et al. (2020) conducted a genetic evaluation on Korean Hanwoo steer information raised by the National Institute of Animal Science. They included Hanwoo steers without genomic information and computed basic statistics for CWT, EMA, BFT, and MS. The reported scores were 371.63 ± 42.07 kg, 81.65 ± 8.82 cm2, 9.71 ± 3.82 mm, and 3.58 ± 1.62 score, respectively. When compared with the results of this study, the reference livestock basic statistics value of Kim et al. (2022) showed similar values to the study’s results, although it was higher than the findings of Park et al. (2020). Park et al. (2020) used nationally managed steering data from 1997 to 2018, which was thought to be the result of including relatively young individuals.

To test the normality of the complete reference population, the findings of quantile-quantile plots (Q-Q plots) for each trait (CWT, EMA, BFT, and MS) are shown in Fig. 1. As a result of the Q-Q plots, outlier values that strayed from the line were identified, but they were included in the analysis because they represented the diversity of the individual.

Figure 1.Quantile-quantile plot (Q-Q plot) for the four carcass traits in the ssGBLUP reference population.

Estimation of genetic parameter

Genetic parameters were estimated by the ssGBLUP method. Table 3 presents the estimation results of the genetic parameters for genetic variance, residual variance, and heritability. The genetic variance was estimated to be 980.84, 50.594, 15.072, and 1.609 for CWT, EMA, BFT, and MS, respectively. The residual variances were estimated to be 1530.4, 97.642, 17.660, and 1.894, and the heritability values using genetic and residual variances were estimated to be 0.391, 0.341, 0.460, and 0.459, respectively. After reviewing previous studies, Kim et al. (2022) estimated the ssGBLUP accuracy based on the amount of genomic information available in the pedigree. They reported heritability values of 0.429, 0.372, 0.357, and 0.471 for the CWT, EMA, BFT, and MS treatments, respectively. Lopez et al. (2020) found heritability values of 0.38, 0.34, 0.38, and 0.48 in an imputation accuracy study using the Bovine HD BeadChip (777 k) and Hanwoo 50 K SNP beadchip. Compared to the results of this study, the analysis included the phenotype of cows, which is thought to be the reason why the BFT heritability was estimated to be low. According to Park et al. (2020), heritability estimates for CWT, EMA, BFT, and MS were 0.35, 0.44, 0.46, and 0.56, respectively. In contrast to the results of this study, the heritability of BFT and CWT was low because the population used in the previous study differed from the population used in this study by using nationally managed individuals.

Table 3. Genetic parameters estimated using the ssGBLUP method based on genotype and pedigree information.

Analysis methodTypeTrait

CWTEMABFTMS
ssGBLUP methodGenetic
Variance
980.8450.59415.0721.609
Residual
Variance
1,530.497.64217.6601.894
Heritability0.3910.3410.4600.459

ssGBLUP, single-step genomic best linear unbiased prediction; CWT, carcass weight; EMA, eye muscle area; BFT, back fat thickness; MS, marbling score..



Accuracy analysis of breeding value

Table 4 lists the GEBV, GEBV accuracy, and realized accuracy of the subjects, estimated using the ssGBLUP method. The mean GEBV estimated using the ssGBLUP methods for CWT, EMA, BFT, and MS were 7.348, 1.515, -0.355, and 0.040, respectively and the accuracy of each trait was 0.749, 0.733, 0.769, and 0.768, respectively. The accuracy of GEBV was highest for BFT and lowest for EMA. Because the breeding value accuracy was affected by the estimated heritability, it showed the same pattern as that of the previous heritability. Kim et al. (2022) conducted a comparative study of the accuracy of ssGBLUP based on the amount of genomic information in the pedigree. They observed that the highest accuracies among the GEBV were 0.736, 0.718, 0.713, and 0.747 for CWT, EMA, BFT, and MS, respectively. Park et al. (2020) reported that for CWT, EMA, BFT, and MS, the accuracies of the estimated breeding values using the ssGBLUP method were 0.68, 0.69, 0.69, and 0.73, respectively. Lopez et al. (2020) reported that the breeding value accuracy estimated using the ssGBLUP method through single-trait analysis was 0.47, 0.42, 0.33, and 0.37 for CWT, EMA, BFT, and MS, respectively. Previous studies have reported a lower accuracy of the mean GEBV than the results obtained in this study. This difference is attributed to the larger size of the reference population used in this study compared to previous studies. Because the ssGBLUP method estimates GEBV using individual genotype information and pedigree phenotypic information, the larger the reference population, the more diverse the information can be explained.

Table 4. Breeding value, accuracy, and realized accuracy were estimated using the ssGBLUP method.

Analysis methodReference population sizeTypeTrait

CWT (kg)EMA (cm2)BFT (mm)MS (score)
ssGBLUP method18,499 + 74,730GEBV7.3481.515-0.3550.040
Standard
Error
20.6894.8242.4750.809
Accuracy0.7490.7330.7690.768
Realized accuracy0.5190.4350.4440.543

ssGBLUP, single-step genomic best linear unbiased prediction; CWT, carcass weight; EMA, eye muscle area; BFT, back fat thickness; MS, marbling score; GEBV, genomic estimated breeding value..



After confirming the correlation analysis (realized accuracy) between the GEBV and the actual slaughter performance, CWT, EMA, BFT, and MS were reported to be 0.519, 0.435, 0.444, and 0.543, respectively. The realized accuracy was the highest for MS and the lowest for EMA. In a previous study, Lee et al. (2023) demonstrated that the realized accuracy values in the genetic ability evaluation of a commercial Hanwoo population using the GBLUP method were 0.522, 0.404, 0.444, and 0.539 for CWT, EMA, BFT, and MS. Kim (2021) reported the realized accuracy of GEBV as 0.509, 0.446, 0.415, and 0.555, respectively, when using the ssGBLUP method based on the phenotypic information of 202,808 heads and the genotype information of 14,000 heads. When comparing the results of this study, the realized accuracy values demonstrated an increase compared with the GBLUP method. Collecting more phenotypic data related to pedigree will enable more accurate genetic evaluation.

DISCUSSION

The results of this study show that the ssGBLUP method allows for a more accurate evaluation of individuals because it incorporates into its analysis not only the genotype and phenotype information used in the current GBLUP method but also the pedigree and phenotype information of the reference population. Currently, the Korean Hanwoo industry is experiencing difficulties, such as rising feed prices and inflation, owing to changes in the international situation. To ensure a stable income, it is necessary to enhance the average genetic ability of Hanwoo cattle raised on farms to produce stable, high-grade Hanwoo cattle. Individual evaluation using the ssGBLUP method is considered effective for enhancing the genetic ability of farms and enabling accurate and rapid improvements. This study used the two generations of pedigree information on the reference population for analysis. Based on the results of this study, it is considered that obtaining further pedigree information of the reference population will enable high-accuracy individual selection.

Acknowledgements

None.

Author Contributions

Conceptualization, H.S.K.; methodology, G.H.L., K.T., Y.S.L., H.S.K.; investigation, G.H.L., and K.T.; data curation, G.H.L., K.T.; writing—original draft preparation, G.H.L., K.T.; writing—review and editing, G.H.L., K.T., Y.S.L., H.S.K.; supervision, H.S.K.; project administration, H.S.K.; funding acquisition, H.S.K.

Funding

None.

Ethical Approval

The study was approved by the Hankyong National University Animal Ethics Committee (No.2021-3).

Consent to Participate

Not applicable.

Consent to Publish

Not applicable.

Availability of Data and Materials

Not applicable.

Conflicts of Interest

No potential conflict of interest relevant to this article was reported.

Fig 1.

Figure 1.Quantile-quantile plot (Q-Q plot) for the four carcass traits in the ssGBLUP reference population.
Journal of Animal Reproduction and Biotechnology 2023; 38: 268-274https://doi.org/10.12750/JARB.38.4.268

Table 1 . Total number of references and test population.

SexTest populationReference population

PedigreeGenotype
Cow1,22136,327355
Steer51938,40318,144
Total1,74074,73018,499

Table 2 . Basic statistics of reference population according to the ssGBLUP method.

SubsetSexNo.TraitsMeanS.D.MaxMin
Pedigree + genotypeCow27,264CWT (kg)349.3754.13667110
EMA (cm2)83.8612.6714722
BFT (mm)13.806.09771
MS (score)4.011.9391
Steer35,303CWT (kg)442.0255.20745148
EMA (cm2)95.1512.3317120
BFT (mm)13.805.10511
MS (score)5.821.9591
GenotypeCow355CWT (kg)356.1647.94498159
EMA (cm2)87.8712.1711938
BFT (mm)13.895.03374
MS (score)4.951.9991
Steer18,144CWT (kg)444.0149.57692160
EMA (cm2)96.7012.2516020
BFT (mm)14.274.93471
MS (score)5.981.8791

CWT, carcass weight; EMA, eye muscle area; BFT, back fat thickness; MS, marbling score; SD, standard deviation..


Table 3 . Genetic parameters estimated using the ssGBLUP method based on genotype and pedigree information.

Analysis methodTypeTrait

CWTEMABFTMS
ssGBLUP methodGenetic
Variance
980.8450.59415.0721.609
Residual
Variance
1,530.497.64217.6601.894
Heritability0.3910.3410.4600.459

ssGBLUP, single-step genomic best linear unbiased prediction; CWT, carcass weight; EMA, eye muscle area; BFT, back fat thickness; MS, marbling score..


Table 4 . Breeding value, accuracy, and realized accuracy were estimated using the ssGBLUP method.

Analysis methodReference population sizeTypeTrait

CWT (kg)EMA (cm2)BFT (mm)MS (score)
ssGBLUP method18,499 + 74,730GEBV7.3481.515-0.3550.040
Standard
Error
20.6894.8242.4750.809
Accuracy0.7490.7330.7690.768
Realized accuracy0.5190.4350.4440.543

ssGBLUP, single-step genomic best linear unbiased prediction; CWT, carcass weight; EMA, eye muscle area; BFT, back fat thickness; MS, marbling score; GEBV, genomic estimated breeding value..


References

  1. Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. 2015. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience 4:7.
    Pubmed KoreaMed CrossRef
  2. Chen CY, Misztal I, Aguilar I, Legarra A, Muir WM. 2011. Effect of different genomic relationship matrices on accuracy and scale. J. Anim. Sci. 89:2673-2679.
    Pubmed CrossRef
  3. Choi TJ, Alam M, Cho CI, Lee JG, Park B, Kim S, Koo Y, Roh SH. 2015. Genetic parameters for yearling weight, carcass traits, and primal-cut yields of Hanwoo cattle. J. Anim. Sci. 93:1511-1521.
    Pubmed CrossRef
  4. Choy YH, Park BH, Choi TJ, Choi JG, Cho KH, Lee SS, Choi YL, Koh KC, Kim HS. 2012. Estimation of relative economic weights of hanwoo carcass traits based on carcass market price. Asian-Australas. J. Anim. Sci. 25:1667-1673.
    Pubmed KoreaMed CrossRef
  5. Forni S, Aguilar I, Misztal I. 2011. Different genomic relationship matrices for single-step analysis using phenotypic, pedigree and genomic information. Genet. Sel. Evol. 43:1.
    Pubmed KoreaMed CrossRef
  6. Gao N, Teng J, Pan R, Li X, Ye S, Li J, Zhang H, Zhang X, Zhang Z. 2019. Accuracy of whole genome prediction with single-step GBLUP in a Chinese yellow-feathered chicken population. Livest. Sci. 230:103817.
    CrossRef
  7. Garrick DJ. 2007. Equivalent mixed model equations for genomic selection. J. Anim. Sci. 85:376.
  8. Gordo DG, Espigolan R, Tonussi RL, Júnior GA, Bresolin T, Magalhães AF, Feitosa FL, Baldi F, Carvalheiro R, Tonhati H, de Oliveira HN, Chardulo LA, de Albuquerque LG. 2016. Genetic parameter estimates for carcass traits and visual scores including or not genomic information. J. Anim. Sci. 94:1821-1826.
    Pubmed CrossRef
  9. Henderson CR. 1975. Best linear unbiased estimation and prediction under a selection model. Biometrics 31:423-447.
    Pubmed CrossRef
  10. Henderson CR. 1984. Applications of Linear Models in Animal Breeding, University of Guelph, Guelph.
  11. Henderson CR and Quaas RL. 1976. Multiple trait evaluation using relatives' records. J. Anim. Sci. 43:1188-1197.
    CrossRef
  12. Kim DH. 2021. Studies on the genetic evaluation of Hanwoo using the genomic information. [Doctoral dissertation, Jeonbuk National University]. RISS. http://www.riss.kr/link?id=T15777082
  13. Kim EH, Sun DW, Kang HC, Myung CH, Kim JY, Lee DH, Lee SH, Lim HT. 2022. Estimated of genomic estimated breeding value and accuracy analysis according to the amount of genotypes in the full-sib family. J. Agric. Life Sci. 56:171-178.
    CrossRef
  14. Kim SJ, Choi TJ, Son JI, Lee DM, Lee JJ, Lee JG, Lim HT, Koo YM. 2022. Estimation of genomic estimated breeding value(GEBV) and reliability for Hanwoo carcass traits using ssGBLUP. J. Anim. Breed. Genom. 6:57-72.
    CrossRef
  15. Lee GH, Lee YH, Kong HS. 2023. Genetic evaluation and accuracy analysis of commercial Hanwoo population using genomic data. J. Anim. Reprod. Biotechnol. 38:32-37.
    CrossRef
  16. Lee J, Cheng H, Garrick D, Golden B, Dekkers J, Park K, Lee D, Fernando R. 2017. Comparison of alternative approaches to single-trait genomic prediction using genotyped and non-genotyped Hanwoo beef cattle. Genet. Sel. Evol. 49:2.
    Pubmed KoreaMed CrossRef
  17. Lopez BI, Lee SH, Shin DH, Oh JD, Chai HH, Park W, Park JE, Lim D. 2020. Accuracy of genomic evaluation using imputed high-density genotypes for carcass traits in commercial Hanwoo population. Livest. Sci. 241:104256.
    CrossRef
  18. Meuwissen TH, Hayes BJ, Goddard ME. 2001. Prediction of total genetic value using genome-wide dense marker maps. Genetics 157:1819-1829.
    Pubmed KoreaMed CrossRef
  19. Misztal I, Lourenco D, Legarra A. 2020. Current status of genomic evaluation. J. Anim. Sci. 98:skaa101.
    Pubmed KoreaMed CrossRef
  20. Misztal I, Tsuruta S, Lourenco D, Masuda Y, Aguilar I, Legarra A, Vitezica Z. 2014. Manual for BLUPF90 Family of Programs, University of Georgia, Athens.
  21. Park MN, Alam M, Kim S, Park B, Lee SH, Lee SS. 2020. Genomic selection through single-step genomic best linear unbiased prediction improves the accuracy of evaluation in Hanwoo cattle. Asian-Australas. J. Anim. Sci. 33:1544-1557.
    Pubmed KoreaMed CrossRef
  22. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC. 2007. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81:559-575.
    Pubmed KoreaMed CrossRef
  23. Song H, Zhang J, Zhang Q, Ding X. 2019. Using different single-step strategies to improve the efficiency of genomic prediction on body measurement traits in pig. Front. Genet. 9:730.
    Pubmed KoreaMed CrossRef
  24. VanRaden PM. 2007. Genomic measures of relationship and inbreeding. Interbull Bull. 37:33-36.
  25. VanRaden PM. 2012. Avoiding bias from genomic pre-selection in converting daughter information across countries. Interbull Bull. 45.
  26. VanRaden PM. 2020. Symposium review: how to implement genomic selection. J. Dairy Sci. 103:5291-5301.
    Pubmed CrossRef
  27. Xiang T, Nielsen B, Su G, Legarra A, Christensen OF. 2016. Application of single-step genomic evaluation for crossbred performance in pig. J. Anim. Sci. 94:936-948.
    Pubmed CrossRef
  28. Zhang Z, Todhunter RJ, Buckler ES, Van Vleck LD. 2007. Technical note: use of marker-based relationships with multiple-trait derivative-free restricted maximal likelihood. J. Anim. Sci. 85:881-885.
    Pubmed CrossRef

JARB Journal of Animal Reproduction and Biotehnology

qr code

OPEN ACCESS pISSN: 2671-4639
eISSN: 2671-4663