Genetic parameters of carcass and meat quality traits in different muscles (longissimus dorsi and semimembranosus) of Hanwoo (Korean cattle)

Estimation of genetic correlations and genomic prediction accuracy for reproductive and carcass traits in Hanwoo cows

This study estimated the heritabilities (h2) and genetic and phenotypic correlations between reproductive traits, including calving interval (CI), age at first calving (AFC), gestation length (GL), number of artificial inseminations per conception (NAIPC), and carcass traits, including carcass weight (CWT), eye muscle area (EMA), backfat thickness (BF), and marbling score (MS) in Korean Hanwoo cows. In addition, the accuracy of genomic predictions of breeding values was evaluated by applying the genomic best linear unbiased prediction (GBLUP) and the weighted GBLUP (WGBLUP) method. The phenotypic data for reproductive and carcass traits were collected from 1,544 Hanwoo cows, and all animals were genotyped using Illumina Bovine 50K single nucleotide polymorphism (SNP) chip. The genetic parameters were estimated using a multi-trait animal model using the MTG2 program. The estimated h2 for CI, AFC, GL, NAIPC, CWT, EMA, BF, and MS were 0.10, 0.13, 0.17, 0.11, 0.37, 0.35, 0.27, and 0.45, respectively, according to the GBLUP model. The GBLUP accuracy estimates ranged from 0.51 to 0.74, while the WGBLUP accuracy estimates for the traits under study ranged from 0.51 to 0.79. Strong and favorable genetic correlations were observed between GL and NAIPC (0.61), CWT and EMA (0.60), NAIPC and CWT (0.49), AFC and CWT (0.48), CI and GL (0.36), BF and MS (0.35), NAIPC and EMA (0.35), CI and BF (0.30), EMA and MS (0.28), CI and AFC (0.26), AFC and EMA (0.24), and AFC and BF (0.21). The present study identified low to moderate positive genetic correlations between reproductive and CWT traits, suggesting that a heavier body weight may lead to a longer CI, AFC, GL, and NAIPC. The moderately positive genetic correlation between CWT and AFC, and NAIPC, with a phenotypic correlation of nearly zero, suggesting that the genotype-environment interactions are more likely to be responsible for the phenotypic manifestation of these traits. As a result, the inclusion of these traits by breeders as selection criteria may present a good opportunity for developing a selection index to increase the response to the selection and identification of candidate animals, which can result in significantly increased profitability of production systems.

Parameter estimation and assessment of bias in genetic evaluation of carcass traits in Hanwoo cattle using real and simulated data

Most carcass and meat quality traits are moderate to highly heritable, indicating that they can be improved through selection. Genetic evaluation for these types of traits is performed using performance data obtained from commercial and progeny testing evaluation. The performance data from commercial farms are available in large volume, however, some drawbacks have been observed. The drawback of the commercial data is mainly due to sorting of animals based on live weight prior to slaughter, and this could lead to bias in the genetic evaluation of later measured traits such as carcass traits. The current study has two components to address the drawback of the commercial data. The first component of the study aimed to estimate genetic parameters for carcass and meat quality traits in Korean Hanwoo cattle using a large sample size of industry-based carcass performance records (n = 469,002). The second component of the study aimed to describe the impact of sorting animals into different contemporary groups based on an early measured trait and then examine the effect on the genetic evaluation of subsequently measured traits. To demonstrate our objectives, we used real performance data to estimate genetic parameters and simulated data was used to assess the bias in genetic evaluation. The results of our first study showed that commercial data obtained from slaughterhouses is a potential source of carcass performance data and useful for genetic evaluation of carcass traits to improve beef cattle performance. However, we observed some harvesting effect which leads to bias in genetic evaluation of carcass traits. This is mainly due to the selection of animal based on their body weight before arrival to slaughterhouse. Overall, the non-random allocation of animals into a contemporary group leads to a biased estimated breeding value in genetic evaluation, the severity of which increases when the evaluation traits are highly correlated.

Development of nutrigenomic based precision management model for Hanwoo steers

Focusing high marble deposition, Hanwoo feedlot system uses high-energy diet over the prolonged fattening period. However, due to the individual genetic variation, around 40% of them are graded into inferior quality grades (QG), despite they utilized the same resources. Therefore, focusing on development of a nutrigenomic based precision management model, this study was to evaluated the response to the divergent selection on genetic merit for marbling score (MS), under different dietary total digestible nutrient (TDN) levels. Total of 111 calves were genotyped and initially grouped according to estimated breeding value (high and low) for marbling score (MS-EBV). Subsequently, managed under two levels of feed TDN%, over the calf period, early, middle, and final fattening periods following 2 × 2 factorial arrangement. Carcasses were evaluated for MS, Back fat thickness (BFT) and Korean beef quality grading standard. As the direct response to the selection was significant, the results confirmed the importance of initial genetic grouping of Hanwoo steers for MS-EBV. However, dietary TDN level did not show an effect (p > 0.05) on the MS. Furthermore, no genetic-by-nutrition interaction for MS (p > 0.05) was also observed. The present results showed no correlation response on BFT (p > 0.05), which indicates that the selection based on MS-EBV can be used to enhance the MS without undesirable effect on BFT. Ultimate turnover of the Hanwoo feedlot operation is primarily determined by the QGs. The present model shows that the initial grouping for MS-EBV increased the proportion of carcasses graded for higher QGs (QG1⁺⁺ and QG1⁺) by approximately 20%. Moreover, there appear to be a potential to increase the proportion of QG 1⁺⁺ animals among the high-genetic group by further increasing the dietary energy content. Overall, this precision management strategy suggests the importance of adopting an MS based initial genetic grouping system for Hanwoo steers with a subsequent divergent management based on dietary energy level.

The study of "muscle eye" in bulls of Ukrainian black-spotted dairy-meat breed as a factor in improving the properties of meat products

The impact of age, live weight, and growth rate of the bulls of Ukrainian breeds on the area of “muscle eye” (cross-section of m. longissimus dorsi when the carcass is divided into front and rear between the 12th and 13th ribs) was studied. The correlation between the size of the “muscle eye” and the carcass's characteristics and the meat's qualitative indicators was also determined. The research was conducted on the bulls of Ukrainian black-and-white dairy (UBWDB) and Ukrainian meat (UMB) breeds. Living animals “muscle eye” area was determined with the ultrasonic analyser Emperor 860, after slaughter. It was found that UMB bulls have the area of “muscle eye” twice as big as their UBWDB peers. The “muscle eye” area increases when growing the cattle to 400 – 450 kg. In the future, it will be practically independent of the age and weight of the animals and remains stable. An increase in the average daily gains within the breed leads to an increase in the “muscle eye” area. The area of “muscle eye” has a weak negative connection (r = -0.193) with meat tenderness and dry matter content (r = -0.345) and a positive one with slaughter weight (r = 0.614) and slaughter yield (r = 0.653). Of the three parameters (length, depth, and area) of “muscle eye”, the greatest impact on the technological properties of meat has depth. Its increase has a negative connection with meat tenderness (r = -0.810) and moisture (r = -0.474), but it has a positive impact on the moisture retention capacity (r = 0.338) and weight of weighed portion after heat treatment. The obtained results can be used to clarify the optimal growing parameters of the bulls of Ukrainian black-and-white dairy and meat breeds for meat and determine the optimal age and live weight of the cattle slaughter.

Estimation of Genetic Correlations of Primal Cut Yields with Carcass Traits in Hanwoo Beef Cattle

Simple Summary

Recently, there is a growing interest in the genetic improvement of carcass merit in the Korean beef industry. Primal cut yields have been proposed to characterize the meat quality and quantity in beef cattle and are known to be genetically correlated with the carcass merit and premium prices. Hence, knowledge of genetic parameters is required to select the weight of primal cuts that may be used as the selection criteria for designing future breeding programs. This study aimed to estimate the heritability and genetic and phenotypic correlations of primal cut yields and carcass traits in Hanwoo cattle. All traits presented a medium to high heritability, which indicates a probable increase in their response to selection. In addition, moderate to highly favorable genetic correlations were found between the primal cut yields and carcass traits, such as the carcass weight and eye muscle area. Therefore, results suggested inclusion of primal cut traits as a selection objective in Hanwoo breeding programs to meet the growing demand for high quality products.

Abstract

This study was carried out to estimate the variance components, heritability, and genetic correlations between the carcass traits and primal cut yields in Hanwoo cattle. Carcass traits comprising 5622 records included back fat thickness (BFT), carcass weight (CW), eye muscle area (EMA), and marbling score (MS). The 10 primal cut yields from 3467 Hanwoo steers included the tenderloin (TLN), sirloin (SLN), striploin (STLN), chuck (CHK), brisket (BSK), top round (TRD), bottom round (BRD), rib (RB), shank (SK), and flank (FK). In addition, three composite traits were formed by combining primal cut yields as novel traits according to consumer preferences and market price: high-value cuts (HVC), medium-value cuts (MVC), and low-value cuts (LVC). Heritability estimates for the interest of traits were moderate to high, ranging from 0.21 ± 0.04 for CHK to 0.59 ± 0.05 for MS. Except genetic correlations between RB and other primal cut traits, favorable and moderate to high correlations were observed among the yields of primal cut that ranged from 0.38 ± 0.14 (CHK and FK) to 0.93 ± 0.01 (TRD and BRD). Moreover, the estimated genetic correlations of CW and EMA with primal cut yields and three composite traits were positive and moderate to strong, except for BFT, which was negative. These results indicate that genetic progress can be achieved for all traits, and selection to increase the yields of primal cuts can lead to considerable profitability in the Hanwoo beef industry.

Transcriptome profiling analysis of muscle tissue reveals potential candidate genes affecting water holding capacity in Chinese Simmental beef cattle

Water holding capacity (WHC) is an important sensory attribute that greatly influences meat quality. However, the molecular mechanism that regulates the beef WHC remains to be elucidated. In this study, the longissimus dorsi (LD) muscles of 49 Chinese Simmental beef cattle were measured for meat quality traits and subjected to RNA sequencing. WHC had significant correlation with 35 kg water loss (r = − 0.99, p < 0.01) and IMF content (r = 0.31, p < 0.05), but not with SF (r = − 0.20, p = 0.18) and pH (r = 0.11, p = 0.44). Eight individuals with the highest WHC (H-WHC) and the lowest WHC (L-WHC) were selected for transcriptome analysis. A total of 865 genes were identified as differentially expressed genes (DEGs) between two groups, of which 633 genes were up-regulated and 232 genes were down-regulated. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment revealed that DEGs were significantly enriched in 15 GO terms and 96 pathways. Additionally, based on protein–protein interaction (PPI) network, animal QTL database (QTLdb), and relevant literature, the study not only confirmed seven genes (HSPA12A, HSPA13, PPARγ, MYL2, MYPN, TPI, and ATP2A1) influenced WHC in accordance with previous studies, but also identified ATP2B4, ACTN1, ITGAV, TGFBR1, THBS1, and TEK as the most promising novel candidate genes affecting the WHC. These findings could offer important insight for exploring the molecular mechanism underlying the WHC trait and facilitate the improvement of beef quality.

Estimation of Genetic Parameters and Correlation between Yearling Ultrasound Measurements and Carcass Traits in Hanwoo Cattle

Simple Summary

Knowledge of genetic parameters is essential to obtain breeding values in order to increase the response to selection and incorporate novel traits in designing a breeding program. There is a growing demand for the genetic improvement of carcass traits in the Korean beef industry. The use of yearling ultrasound measurements as indicator traits can be an efficient way to evaluate carcass traits. To date, the assessment of genetic parameters for ultrasound measurements in Hanwoo cattle is still limited. Therefore, this study was conducted to estimate the heritability, and the genetic and phenotypic correlations of yearling ultrasonic and carcass traits in Hanwoo cattle. The results revealed moderate to high heritability estimates for the traits of interest, which indicate a probable increase in the response to selection for these traits. Moreover, high and favorable genetic correlations were observed between carcass traits and their corresponding ultrasound measurements. Our findings suggest that the inclusion of yearling ultrasound data on potential replacements would be suitable as a selection tool for genetic improvement of carcass traits in Hanwoo breeding programs.

Abstract

Genetic parameters have a significant role in designing a breeding program and are required to evaluate economically important traits. The objective of this study was to estimate heritability and genetic correlation between yearling ultrasound measurements, such as backfat thickness (UBFT), eye muscle area (UEMA), intramuscular fat content (UIMF), and carcass traits, such as backfat thickness (BFT), carcass weight (CW), eye muscle area (EMA), marbling score (MS) at approximately 24 months of age, as well as yearling weight (YW) in Hanwoo bulls (15,796) and steers (5682). The (co) variance components were estimated using a multi-trait animal model. Moderate to high heritability estimates were obtained and were 0.42, 0.50, 0.56, and 0.59 for CW, EMA, BFT, and MS, respectively. Heritability estimates for yearling measurements of YW, UEMA, UBFT, and UIMF were 0.31, 0.32, 0.30, and 0.19, respectively. Favorable and strong genetic correlations were observed between UIMF and MS (0.78), UBFT and BFT (0.63), and UEMA and EMA (0.65). Moreover, the estimated genetic correlation between YW and CW was high (0.84) and relatively moderate between YW and EMA (0.43). These results suggest that genetic improvement can be achieved for carcass traits when using yearling ultrasound measurements as selection criteria in ongoing Hanwoo breeding programs.

RNA-seq analysis identifies cytoskeletal structural genes and pathways for meat quality in beef

RNA sequencing (RNA-seq) has allowed for transcriptional profiling of biological systems through the identification of differentially expressed (DE) genes and pathways. A total of 80 steers with extreme phenotypes were selected from the University of Florida multibreed Angus-Brahman herd. The average slaughter age was 12.91±8.69 months. Tenderness, juiciness and connective tissue assessed by sensory panel, along with marbling, Warner-Bratzler Shear Force (WBSF) and cooking loss, were measured in longissimus dorsi muscle. Total RNA was extracted from muscle and one RNA-seq library per sample was constructed, multiplexed, and sequenced based on protocols by Illumina HiSeq-3000 platform to generate 2×101 bp paired-end reads. The overall read mapping rate using the Btau_4.6.1 reference genome was 63%. A total of 8,799 genes were analyzed using two different methodologies, an expression association and a DE analysis. A gene and exon expression association analysis was carried out using a meat quality index on all 80 samples as a continuous response variable. The expression of 208 genes and 3,280 exons from 1,565 genes was associated with the meat quality index (p-value ≤ 0.05). A gene and isoform DE evaluation was performed analyzing two groups with extreme WBSF, tenderness and marbling. A total of 676 (adjusted p-value≤0.05), 70 (adjusted p-value≤0.1) and 198 (adjusted p-value≤0.1) genes were DE for WBSF, tenderness and marbling, respectively. A total of 106 isoforms from 98 genes for WBSF, 13 isoforms from 13 genes for tenderness and 43 isoforms from 42 genes for marbling (FDR≤0.1) were DE. Cytoskeletal and transmembrane anchoring genes and pathways were identified in the expression association, DE and the gene enrichment analyses; these proteins can have a direct effect on meat quality. Cytoskeletal proteins and transmembrane anchoring molecules can influence meat quality by allowing cytoskeletal interaction with myocyte and organelle membranes, contributing to cytoskeletal structure and architecture maintenance postmortem.

Genetic parameters and correlations of related feed efficiency, growth, and carcass traits in Hanwoo beef cattle

Objective

This study aimed to estimate the genetic parameters and genetic correlations for related feed efficiency, growth, and carcass traits in Hanwoo cattle.

Methods

Phenotypic data from 15,279 animals born between 1989 and 2015 were considered. The related feed efficiency traits considered were Kleiber ratio (KR) and relative growth rate (RGR). Carcass traits analyzed were backfat thickness (BT), carcass weight, eye muscle area, and marbling score. Growth traits were assessed by the average daily gain (ADG), metabolic body weight (MBW) at mid-test age from 6 to 24 months, and yearling weight (YW). Variance and covariance components were estimated using restricted maximum likelihood using nine multi-trait animal models.

Results

The heritability estimates for related feed efficiency (0.28±0.04 for KR and RGR) and growth traits (0.26±0.02 to 0.33±0.04) were moderate, but the carcass traits tended to be higher (0.38±0.04 to 0.61±0.06). The related feed efficiency traits were positively genetically correlated with all the carcass traits (0.37±0.09 to 0.47±0.07 for KR, and 0.14±0.09 to 0.37±0.09 for RGR), except for BT, which showed null to weak correlation. Conversely, the genetic correlations of RGR with MBW (−0.36±0.08) and YW (−0.30±0.08) were negative, and those of KR with MBW and YW were close to zero, whereas the genetic correlations of ADG with RGR (0.40±0.08) and KR (0.70±0.05) were positive and relatively moderate to high. The genetic (0.92±0.02) correlations between KR and RGR were very high.

Conclusion

Sufficient genetic variability and heritability were observed for traits of interest. Moreover, the inclusion of KR and/or RGR in Hanwoo cattle breeding programs could improve the feed efficiency without producing any unfavorable effects on the carcass traits.

Haplotype-Based Genome-Wide Association Study and Identification of Candidate Genes Associated with Carcass Traits in Hanwoo Cattle

Hanwoo, is the most popular native beef cattle in South Korea. Due to its extensive popularity, research is ongoing to enhance its carcass quality and marbling traits. In this study we conducted a haplotype-based genome-wide association study (GWAS) by constructing haplotype blocks by three methods: number of single nucleotide polymorphisms (SNPs) in a haplotype block (nsnp), length of genomic region in kb (Len) and linkage disequilibrium (LD). Significant haplotype blocks and genes associated with them were identified for carcass traits such as BFT (back fat thickness), EMA (eye Muscle area), CWT (carcass weight) and MS (marbling score). Gene-set enrichment analysis and functional annotation of genes in the significantly-associated loci revealed candidate genes, including PLCB1 and PLCB4 present on BTA13, coding for phospholipases, which might be important candidates for increasing fat deposition due to their role in lipid metabolism and adipogenesis. CEL (carboxyl ester lipase), a bile-salt activated lipase, responsible for lipid catabolic process was also identified within the significantly-associated haplotype block on BTA11. The results were validated in a different Hanwoo population. The genes and pathways identified in this study may serve as good candidates for improving carcass traits in Hanwoo cattle.

A Gene-Set Enrichment and Protein-Protein Interaction Network-Based GWAS with Regulatory SNPs Identifies Candidate Genes and Pathways Associated with Carcass Traits in Hanwoo Cattle

Non-synonymous SNPs and protein coding SNPs within the promoter region of genes (regulatory SNPs) might have a significant effect on carcass traits. Imputed sequence level data of 10,215 Hanwoo bulls, annotated and filtered to include only regulatory SNPs (450,062 SNPs), were used in a genome-wide association study (GWAS) to identify loci associated with backfat thickness (BFT), carcass weight (CWT), eye muscle area (EMA), and marbling score (MS). A total of 15, 176, and 1 SNPs were found to be significantly associated (p < 1.11 × 10-7) with BFT, CWT, and EMA, respectively. The significant loci were BTA4 (CWT), BTA6 (CWT), BTA14 (CWT and EMA), and BTA19 (BFT). BayesR estimated that 1.1%~1.9% of the SNPs contributed to more than 0.01% of the phenotypic variance. So, the GWAS was complemented by a gene-set enrichment (GSEA) and protein-protein interaction network (PPIN) analysis in identifying the pathways affecting carcass traits. At p < 0.005 (~2,261 SNPs), 25 GO and 18 KEGG categories, including calcium signaling, cell proliferation, and folate biosynthesis, were found to be enriched through GSEA. The PPIN analysis showed enrichment for 81 candidate genes involved in various pathways, including the PI3K-AKT, calcium, and FoxO signaling pathways. Our finding provides insight into the effects of regulatory SNPs on carcass traits.

Estimation of Genetic Parameters by Single-Trait and Multi-Trait Models for Carcass Traits in Hanwoo Cattle

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.

Genome-Wide Association Study of Meat Quality Traits in Hanwoo Beef Cattle Using Imputed Whole-Genome Sequence Data

The discovery of single nucleotide polymorphisms (SNP) and the subsequent genotyping of large numbers of animals have enabled large-scale analyses to begin to understand the biological processes that underpin variation in animal populations. In beef cattle, genome-wide association studies using genotype arrays have revealed many quantitative trait loci (QTL) for various production traits such as growth, efficiency and meat quality. Most studies regarding meat quality have focused on marbling, which is a key trait associated with meat eating quality. However, other important traits like meat color, texture and fat color have not commonly been studied. Developments in genome sequencing technologies provide new opportunities to identify regions associated with these traits more precisely. The objective of this study was to estimate variance components and identify significant variants underpinning variation in meat quality traits using imputed whole genome sequence data. Phenotypic and genomic data from 2,110 Hanwoo cattle were used. The estimated heritabilities for the studied traits were 0.01, 0.16, 0.31, and 0.49 for fat color, meat color, meat texture and marbling score, respectively. Marbling score and meat texture were highly correlated. The genome-wide association study revealed 107 significant SNPs located on 14 selected chromosomes (one QTL region per selected chromosome). Four QTL regions were identified on BTA2, 12, 16, and 24 for marbling score and two QTL regions were found for meat texture trait on BTA12 and 29. Similarly, three QTL regions were identified for meat color on BTA2, 14 and 24 and five QTL regions for fat color on BTA7, 10, 12, 16, and 21. Candidate genes were identified for all traits, and their potential influence on the given trait was discussed. The significant SNP will be an important inclusion into commercial genotyping arrays to select new breeding animals more accurately.

Comparison of conventional BLUP and single-step genomic BLUP evaluations for yearling weight and carcass traits in Hanwoo beef cattle using single trait and multi-trait models

Hanwoo, an important indigenous and popular breed of beef cattle in Korea, shows rapid growth and has high meat quality. Its yearling weight (YW) and carcass traits (backfat thickness, carcass weight- CW, eye muscle area, and marbling score) are economically important for selection of young and proven bulls. However, measuring carcass traits is difficult and expensive, and can only be performed postmortem. Genomic selection has become an appealing procedure for genetic evaluation of these traits (by inclusion of the genomic data) along with the possibility of multi-trait analysis. The aim of this study was to compare conventional best linear unbiased prediction (BLUP) and single-step genomic BLUP (ssGBLUP) methods, using both single-trait (ST-BLUP, ST-ssGBLUP) and multi-trait (MT-BLUP, MT-ssGBLUP) models to investigate the improvement of breeding-value accuracy for carcass traits and YW. The data comprised of 15,279 phenotypic records for YW and 5,824 records for carcass traits, and 1,541 genotyped animals for 34,479 single-nucleotide polymorphisms. Accuracy for each trait and model was estimated only for genotyped animals by five-fold cross-validation. ssGBLUP models (ST-ssGBLUP and MT-ssGBLUP) showed ~19% and ~36% greater accuracy than conventional BLUP models (ST-BLUP and MT-BLUP) for YW and carcass traits, respectively. Within ssGBLUP models, the accuracy of the genomically estimated breeding value for CW increased (19%) when ST-ssGBLUP was replaced with the MT-ssGBLUP model, as the inclusion of YW in the analysis led to a strong genetic correlation with CW (0.76). For backfat thickness, eye muscle area, and marbling score, ST- and MT-ssGBLUP models yielded similar accuracy. Thus, combining pedigree and genomic data via the ssGBLUP model may be a promising way to ensure acceptable accuracy of predictions, especially among young animals, for ongoing Hanwoo cattle breeding programs. MT-ssGBLUP is highly recommended when phenotypic records are limited for one of the two highly correlated genetic traits.

Functional Partitioning of Genomic Variance and Genome-Wide Association Study for Carcass Traits in Korean Hanwoo Cattle Using Imputed Sequence Level SNP Data

Quantitative traits are usually controlled by numerous genomic variants with small individual effects, and variances associated with those traits are explained in a continuous manner. However, the relative contributions of genomic regions to observed genetic variations have not been well explored using sequence level single nucleotide polymorphism (SNP) information. Here, imputed sequence level SNP data (11,278,153 SNPs) of 2109 Hanwoo steers (Korean native cattle) were partitioned according to functional annotation, chromosome, and minor allele frequency (MAF). Genomic relationship matrices (GRMs) were constructed for each classified region and fitted in the model both separately and together for carcass weight (CWT), eye muscle area (EMA), backfat thickness (BFT), and marbling score (MS) traits. A genome-wide association study (GWAS) was performed to identify significantly associated variants in genic and exon regions using a linear mixed model, and the genetic contribution of each exonic SNP was determined using a Bayesian mixture model. Considering all SNPs together, the heritability estimates for CWT, EMA, BFT, and MS were 0.57 ± 0.05, 0.46 ± 0.05, 0.45 ± 0.05, and 0.49 ± 0.05, respectively, which reflected substantial genomic contributions. Joint analysis revealed that the variance explained by each chromosome was proportional to its physical length with weak linear relationships for all traits. Moreover, genomic variances explained by functional category and MAF class differed greatly among the traits studied in joint analysis. For example, exon regions had larger contributions for BFT (0.13 ± 0.08) and MS (0.22 ± 0.08), whereas intron and intergenic regions explained most of the total genomic variances for CWT and EMA (0.22 ± 0.09–0.32 ± 0.11). Considering different functional classes of exon regions and the per SNP contribution revealed the largest proportion of genetic variance was attributable to synonymous variants. GWAS detected 206 and 27 SNPs in genic and exon regions, respectively, on BTA4, BTA6, and BTA14 that were significantly associated with CWT and EMA. These SNPs were harbored by 31 candidate genes, among which TOX, FAM184B, PPARGC1A, PRKDC, LCORL, and COL1A2 were noteworthy. BayesR analysis found that most SNPs (>93%) had very small effects and the 4.02–6.92% that had larger effects (10^-4 × σA2, 10^-3 × σA2, and 10^-2 × σA2) explained most of the total genetic variance, confirming polygenic components of the traits studied.