Identification of potential biomarkers associated with meat tenderness in Hanwoo (Korean cattle): An expression quantitative trait loci analysis.

Meat tenderness is considered the most important trait contributing to beef quality, level of consumer satisfaction, willingness to pay premium prices and industry profit. Genomic selection method would be helpful for genetic improvement of traits with low heritability and that are difficult to measure. The identification of core genes can aid genomic selection for complex traits with low heritability that are difficult to measure. We performed statistical analysis of associations between longissimus dorsi muscle tenderness and gene expression in 20 Hanwoo cattle, using Warner-Bratzler shear force and RNAseq data, respectively. We found a total of 166 core genes, from which six (ASAP1, CAPN5, ELN, SUMF2, TTC8 and MGAT4A) were regulated by 16 cis-expression quantitative trait loci (eQTL) SNPs. Notably, we found that a cis-eQTL SNP of the ELN gene contained an MFZ-1 binding site in its putative promoter region. These findings provide useful information for genomic prediction of beef tenderness in Hanwoo cattle.

Genomic partitioning of growth traits using a high-density single nucleotide polymorphism array in Hanwoo (Korean cattle).

OBJECTIVE: The objective of this study was to characterize the number of loci affecting growth traits and the distribution of single nucleotide polymorphism (SNP) effects on growth traits, and to understand the genetic architecture for growth traits in Hanwoo (Korean cattle) using genome-wide association study (GWAS), genomic partitioning, and hierarchical Bayesian mixture models. METHODS: GWAS: A single-marker regression-based mixed model was used to test the association between SNPs and causal variants. A genotype relationship matrix was fitted as a random effect in this linear mixed model to correct the genetic structure of a sire family. Genomic restricted maximum likelihood and BayesR: A priori information included setting the fixed additive genetic variance to a pre-specified value; the first mixture component was set to zero, the second to 0.0001×σ_g^2, the third 0.001 × σ_g^2, d the fourth to 0.01 × σ_g^2. BayesR fixed a priori information was not more than 1% of the genetic variance for each of the SNPs affecting the mixed distribution. RESULTS: The GWAS revealed common genomic regions of 2 Mb on bovine chromosome 14 (BTA14) and 3 had a moderate effect that may contain causal variants for body weight at 6, 12, 18, and 24 months. This genomic region explained approximately 10% of the variance against total additive genetic variance and body weight heritability at 12, 18, and 24 months. BayesR identified the exact genomic region containing causal SNPs on BTA14, 3, and 22. However, the genetic variance explained by each chromosome or SNP was estimated to be very small compared to the total additive genetic variance. Causal SNPs for growth trait on BTA14 explained only 0.04% to 0.5% of the genetic variance. CONCLUSION: Segregating mutations have a moderate effect on BTA14, 3, and 19; many other loci with small effects on growth traits at different ages were also identified.

Major histocompatibility complex genes exhibit a potential immunological role in mixed Eimeria-infected broiler cecum analyzed using RNA sequencing.

OBJECTIVE: This study was conducted to investigate the differential expression of the major histocompatibility complex (MHC) gene region in Eimeria-infected broiler. METHODS: We profiled gene expression of Eimeria-infected and uninfected ceca of broilers sampled at 4, 7, and 21 days post-infection (dpi) using RNA sequencing. Differentially expressed genes (DEGs) between two sample groups were identified at each time point. DEGs located on chicken chromosome 16 were used for further analysis. Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis was conducted for the functional annotation of DEGs. RESULTS: Fourteen significant (false discovery rate <0.1) DEGs were identified at 4 and 7 dpi and categorized into three groups: MHC-Y class I genes, MHC-B region genes, and non-MHC genes. In Eimeria-infected broilers, MHC-Y class I genes were upregulated at 4 dpi but downregulated at 7 dpi. This result implies that MHC-Y class I genes initially activated an immune response, which was then suppressed by Eimeria. Of the MHC-B region genes, the DMB1 gene was upregulated, and TAP-related genes significantly implemented antigen processing for MHC class I at 4 dpi, which was supported by KEGG pathway analysis. CONCLUSION: This study is the first to investigate MHC gene responses to coccidia infection in chickens using RNA sequencing. MHC-B and MHC-Y genes showed their immune responses in reaction to Eimeria infection. These findings are valuable for understanding chicken MHC gene function.

Utilization of Electrical Conductivity to Improve Prediction Accuracy of Cooking Loss of Pork Loin.

This study investigated the predictability of cooking loss of pork loin through relatively easy and quick measurable quality properties. The pH, color, moisture, protein content, and cooking loss of 100 pork loins were measured. The explanatory variables included in all linear regression models with an adjust-r2 value of ≥0.5 were pH and the protein content. In the linear regression model predicting cooking loss, the highest adjust-r2 value was 0.7, with pH, CIE L*, CIE b*, moisture, and protein content as the explanatory variables. In 30 pork loins, electrical conductivity was additionally measured, and as a result of linear regression analysis for predicting cooking loss, the highest adjust-r2 value was 0.646 with electrical conductivity measured at 40 Hz, with pH and color as the explanatory variables. Ordinal logistic regression analysis was performed to predict the three grades (low, middle, and high) of loin cooking loss using pH, color, and 40 Hz electrical conductivity as the explanatory variables, and the percent concordance was 93.8%. In conclusion, the addition of electrical conductivity as an explanatory variable did not increase the prediction accuracy of the linear regression model for predicting cooking loss; however, it was demonstrated that it is possible to predict and classify the cooking loss grade of pork loin through quality properties that can be measured quickly and easily.

Identification of Missense Variants Affecting Carcass Traits for Hanwoo Precision Breeding.

This study aimed to identify causal variants associated with important carcass traits such as weight and meat quality in Hanwoo cattle. We analyzed missense mutations extracted from imputed sequence data (ARS-UCD1.2) and performed an exon-specific association test on the carcass traits of 16,970 commercial Hanwoo. We found 33, 2, 1, and 3 significant SNPs associated with carcass weight (CW), backfat thickness (BFT), eye muscle area (EMA), and marbling score (MS), respectively. In CW and EMA, the most significant missense SNP was identified at 19,524,263 on BTA14 and involved the PRKDC. A missense SNP in the ZFAND2B, located at 107,160,304 on BTA2 was identified as being involved in BFT. For MS, missense SNP in the ACVR2B gene, located at 11,849,704 in BTA22 was identified as the most significant marker. The contribution of the most significant missense SNPs to genetic variance was confirmed to be 8.47%, 2.08%, 1.73%, and 1.19% in CW, BFT, EMA, and MS, respectively. We generated favorable and unfavorable haplotype combinations based on the significant SNPs for CW. Significant differences in GEBV (Genomic Estimated Breeding Values) were observed between groups with each favorable and unfavorable haplotype combination. In particular, the missense SNPs in PRKDC, MRPL9, and ANKFN1 appear to significantly affect the protein's function and structure, making them strong candidates as causal mutations. These missense SNPs have the potential to serve as valuable markers for improving carcass traits in Hanwoo commercial farms.

Use of a graph neural network to the weighted gene co-expression network analysis of Korean native cattle.

In the general framework of the weighted gene co-expression network analysis (WGCNA), a hierarchical clustering algorithm is commonly used to module definition. However, hierarchical clustering depends strongly on the topological overlap measure. In other words, this algorithm may assign two genes with low topological overlap to different modules even though their expression patterns are similar. Here, a novel gene module clustering algorithm for WGCNA is proposed. We develop a gene module clustering network (gmcNet), which simultaneously addresses single-level expression and topological overlap measure. The proposed gmcNet includes a "co-expression pattern recognizer" (CEPR) and "module classifier". The CEPR incorporates expression features of single genes into the topological features of co-expressed ones. Given this CEPR-embedded feature, the module classifier computes module assignment probabilities. We validated gmcNet performance using 4,976 genes from 20 native Korean cattle. We observed that the CEPR generates more robust features than single-level expression or topological overlap measure. Given the CEPR-embedded feature, gmcNet achieved the best performance in terms of modularity (0.261) and the differentially expressed signal (27.739) compared with other clustering methods tested. Furthermore, gmcNet detected some interesting biological functionalities for carcass weight, backfat thickness, intramuscular fat, and beef tenderness of Korean native cattle. Therefore, gmcNet is a useful framework for WGCNA module clustering.

Time-series transcriptome analysis identified differentially expressed genes in broiler chicken infected with mixed Eimeria species.

Coccidiosis caused by the Eimeria species is a highly problematic disease in the chicken industry. Here, we used RNA sequencing to observe the time-dependent host responses of Eimeria-infected chickens to examine the genes and biological functions associated with immunity to the parasite. Transcriptome analysis was performed at three time points: 4, 7, and 21 days post-infection (dpi). Based on the changes in gene expression patterns, we defined three groups of genes that showed differential expression. This enabled us to capture evidence of endoplasmic reticulum stress at the initial stage of Eimeria infection. Furthermore, we found that innate immune responses against the parasite were activated at the first exposure; they then showed gradual normalization. Although the cytokine-cytokine receptor interaction pathway was significantly operative at 4 dpi, its downregulation led to an anti-inflammatory effect. Additionally, the construction of gene co-expression networks enabled identification of immunoregulation hub genes and critical pattern recognition receptors after Eimeria infection. Our results provide a detailed understanding of the host-pathogen interaction between chicken and Eimeria. The clusters of genes defined in this study can be utilized to improve chickens for coccidiosis control.

Accuracy of genotype imputation based on reference population size and marker density in Hanwoo cattle.

Recently, the cattle genome sequence has been completed, followed by developing a commercial single nucleotide polymorphism (SNP) chip panel in the animal genome industry. In order to increase statistical power for detecting quantitative trait locus (QTL), a number of animals should be genotyped. However, a high-density chip for many animals would be increasing the genotyping cost. Therefore, statistical inference of genotype imputation (low-density chip to high-density) will be useful in the animal industry. The purpose of this study is to investigate the effect of the reference population size and marker density on the imputation accuracy and to suggest the appropriate number of reference population sets for the imputation in Hanwoo cattle. A total of 3,821 Hanwoo cattle were divided into reference and validation populations. The reference sets consisted of 50k (38,916) marker data and different population sizes (500, 1,000, 1,500, 2,000, and 3,600). The validation sets consisted of four validation sets (Total 889) and the different marker density (5k [5,000], 10k [10,000], and 15k [15,000]). The accuracy of imputation was calculated by direct comparison of the true genotype and the imputed genotype. In conclusion, when the lowest marker density (5k) was used in the validation set, according to the reference population size, the imputation accuracy was 0.793 to 0.929. On the other hand, when the highest marker density (15k), according to the reference population size, the imputation accuracy was 0.904 to 0.967. Moreover, the reference population size should be more than 1,000 to obtain at least 88% imputation accuracy in Hanwoo cattle.

THI Modulation of Genetic and Non-genetic Variance Components for Carcass Traits in Hanwoo Cattle.

The phenotype of carcass traits in beef cattle are affected by random genetic and non-genetic effects, which both can be modulated by an environmental variable such as Temperature-Humidity Index (THI), a key environmental factor in cattle production. In this study, a multivariate reaction norm model (MRNM) was used to assess if the random genetic and non-genetic (i.e., residual) effects of carcass weight (CW), back fat thickness (BFT), eye muscle area (EMA), and marbling score (MS) were modulated by THI, using 9,318 Hanwoo steers (N = 8,964) and cows (N = 354) that were genotyped on the Illumina Bovine SNP50 BeadChip (50K). THI was measured based on the period of 15-45 days before slaughter. Both the correlation and the interaction between THI and random genetic and non-genetic effects were accounted for in the model. In the analyses, it was shown that the genetic effects of EMA and the non-genetic effects of CW and MS were significantly modulated by THI. No significant THI modulation of such effects was found for BFT. These results highlight the relevance of THI changes for the genetic and non-genetic variation of CW, EMA, and MS in Hanwoo beef cattle. Importantly, heritability estimates for CW, EMA, and MS from additive models without considering THI interactions were underestimated. Moreover, the significance of interaction can be biased if not properly accounting for the correlation between THI and genetic and non-genetic effects. Thus, we argue that the estimation of genetic parameters should be based on appropriate models to avoid any potential bias of estimates. Our finding should serve as a basis for future studies aiming at revealing genotype by environment interaction in estimation and genomic prediction of breeding values.