Detection of candidate polymorphisms around the QTL for fat area ratio to rib eye area on BTA7 using whole-genome resequencing in Japanese Black cattle

Comprehensive genotyping analysis of single nucleotide polymorphisms responsible for beef marbling in Japanese Black cattle

BACKGROUND

Beef marbling is considered a desirable trait in the meat industry. Therefore, understanding the genetic factors that cause marbling is important. Previously, we performed a genome-wide association study to examine genetic factors associated with beef marbling in Japanese Black cattle and identified a candidate region between 10-30 Mbp on chromosome 7. We verified the effect of the SNPs in this region on beef marbling using linkage disequilibrium block analysis. We narrowed down the candidate region to a range of 15.8-16.1 Mbp. In this study, we comprehensively detected all of the SNPs in this region and verified their effects on beef marbling.

RESULTS

Genome resequencing using four animals exhibiting high beef marbling standard (BMS) and four with low BMS revealed a total of 1,846 polymorphisms within the candidate region. Based on the annotation, we selected 13 SNPs exhibiting a moderate impact, as no high-impact SNPs were detected. All of the SNPs represented missense polymorphisms and were located in the following seven genes: RDH8, ANGPTL6, DNMT1, MRPL4, ICAM1, ICAM3, and ICAM5. Finally, we determined the effects of these SNPs on the BMS of a Japanese Black cattle population (n = 529). Analysis of variance revealed that the five SNPs were located in genes encoding the intercellular adhesion molecules (ICAM1, ICAM3, and ICAM5), and showed a highly significant association compared with the remainder (p < 0.01). The lowest p-value was observed for ICAM3_c.739G > A (p = 1.18E-04). Previous studies have suggested that intercellular adhesion molecules (ICAM) may be an upstream factor that regulates adipocyte differentiation. Therefore, considering the polymorphism and putative gene function, we suggest that ICAM1 is potentially responsible for beef marbling. c.470C > G and/or c.994G > A on ICAM1 may be responsible for this quantitative trait locus.

CONCLUSIONS

Promising SNP candidates responsible for beef marbling were identified using extensive polymorphism verification in a previously reported QTL region. We aim to elucidate the mechanism of beef marbling in future studies by investigating how these polymorphisms alter protein structure and function.

Verification of Candidate SNP Effects Reveals Two QTLs on BTA7 for Beef Marbling in Two Japanese Black Cattle Populations

In our previous study, we used genome resequencing to detect all candidate polymorphisms within a quantitative trait loci (QTL) region for beef marbling reported previously at 10–30 Mbp on bovine chromosome 7, and we selected 6044 polymorphisms as candidate quantitative trait nucleotides (QTNs). In the present study, we aimed to identify quantitative trait genes (QTGs) and QTNs in this QTL region by verifying the effect of SNPs on beef marbling in two Japanese Black cattle populations using a Dynamic Array integrated fluidic circuit. In total, 96 selected SNPs were genotyped in 441 and 529 animals in Hyogo and Miyazaki cattle populations, respectively. The most significant p-values were detected in a SNP in a splice region of ALDH7A1 (SNP93_ALDH7A1; p = 3.46 × 10−5) in Hyogo cattle and a missense polymorphism of intercellular adhesion molecule-1 (ICAM1) (SNP37_ICAM1; p = 3.33 × 10−4) in Miyazaki cattle. Interestingly, SNP93_ALDH7A1 was not significant (p = 0.459) in Miyazaki cattle, and SNP37_ICAM1 showed a weakly significant association (p = 0.043) in Hyogo cattle. Thus, each population would likely have different QTGs and QTNs for beef marbling in the QTL region. In the Hyogo population, it was not possible to determine the accurate range of the linkage disequilibrium (LD) block in LD block analysis because of a strong LD structure throughout the assessed region. In Miyazaki cattle, however, an LD block containing SNP37_ICAM1 had a range of 15.8–16.1 Mbp, suggesting that QTNs would be located within this region. The functions of 19 genes in the LD block were investigated. ICAM1 is known to play an important role in adipocyte differentiation; given this function and the effect of amino acid substitution, SNP37_ICAM1 was identified as a promising candidate QTN for beef marbling. Further research on the effect of SNP37_ICAM1 on adipocyte differentiation is expected to provide insights into the mechanism underlying beef marbling formation.

Development of DNA markers for improvement of meat quality in a Japanese Black cattle population in Hyogo Prefecture

The polymorphisms associated with economic traits in livestock animals provide useful information as genetic indicators for breeding improvement. Over the last two decades, several DNA markers have been developed in Japanese Black cattle; however, the effect of these markers differs across populations due to differences in their genetic structures and backgrounds. As such, there is a need to verify the effectiveness of these markers in each population. This review summarizes the effectiveness of previously reported markers on carcass traits and the development of novel DNA markers in a Japanese Black cattle population in Hyogo Prefecture. As result of genome wide association studies and resequencing analyses, two novel significant markers associated with meat quality-related traits (beef marbling and fatty acid composition) were developed. These findings will lead to the identification of responsible genes and polymorphisms and contribute to the development of novel DNA markers for numerous traits in various cattle populations.

Estimation of genetic parameters for carcass grading traits, image analysis traits, and monounsaturated fatty acids in Japanese Black cattle from Hyogo Prefecture

Genetic parameters for carcass grading traits, image analysis traits, and monounsaturated fatty acid (MUFA) percentages were estimated in 29,942 Japanese Black cattle from Hyogo Prefecture. The analyzed traits included five carcass grading traits, two image analysis traits, fat area ratio and fineness index, and two MUFA traits, one measured in intermuscular fat using near-infrared spectroscopy (NIRS) and the other in intramuscular fat using gas chromatography (GC). The heritability estimates of image analysis traits and MUFA were moderate to high, ranging from 0.395 to 0.740, and it was considered that they could be improved simultaneously with carcass grading traits because no severe genetic antagonism was observed. Although the heritability of the NIRS-based intermuscular MUFA was slightly lower than that of the GC-based intramuscular MUFA, the genetic correlation between the two methods was as high as 0.804. These results indicate that the NIRS method can be used as an alternative evaluation procedure to predict MUFA in intramuscular fat in the longissimus muscle.

Gene Expression Analysis Provides New Insights into the Mechanism of Intramuscular Fat Formation in Japanese Black Cattle

Japanese Black cattle (Japanese Wagyu) have a unique phenotype in which ectopic intramuscular fat accumulates in skeletal muscle, producing finely marbled beef. However, the mechanism of intramuscular fat formation in Japanese Black cattle remains unclear. To investigate the key genes involved in intramuscular fat accumulation, we comprehensively analyzed mRNA levels in subcutaneous and intramuscular fat tissues using RNA sequence (RNA-seq) analysis, which detected 27,606 genes. We identified eight key genes, namely carboxypeptidase E, tenascin C, transgelin, collagen type IV alpha 5 (COL4A5), cysteine and glycine-rich protein 2, PDZ, and LIM domain 3, phosphatase 1 regulatory inhibitor subunit 14A, and regulator of calcineurin 2. These genes were highly and specifically expressed in intramuscular fat tissue. Immunohistochemical analysis revealed a collagen network, including COL4A5, in the basement membrane around the intramuscular fat tissue. Moreover, pathway analysis revealed that, in intramuscular fat tissue, differentially expressed genes are related to cell adhesion, proliferation, and cancer pathways. Furthermore, pathway analysis showed that the transforming growth factor-β (TGF-β) and small GTPases regulators RASGRP3, ARHGEF26, ARHGAP10, ARHGAP24, and DLC were upregulated in intramuscular fat. Our study suggests that these genes are involved in intramuscular fat formation in Japanese Black cattle.