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.

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.

In our previous study, we performed genome-wide association study (GWAS) to identify the genomic region associated with Fat area ratio to rib eye area (FAR) and detected a candidate in BTA7 at 10-30 Mbp. The present study aims to comprehensively detect all polymorphisms in the candidate region using whole-genome resequencing data. Based on whole-genome resequencing of eight animals, we detected 127,090 polymorphisms within the region. Of these, 31,945 were located within the genes. We further narrowed the polymorphisms to 6,044 with more than five allele differences between the high and low FAR groups that were located within 179 genes. We subsequently investigated the functions of these genes and selected 170 polymorphisms in eight genes as possible candidate polymorphisms. We focused on SLC27A6 K81M as a putative candidate polymorphism. We genotyped the SNP in a Japanese Black population (n = 904) to investigate the effect on FAR. Analysis of variance revealed that SLC27A6 K81M had a lower p-value (p = .0009) than the most significant SNP in GWAS (p = .0049). Although only SLC27A6 K81M was verified in the present study, subsequent verification of the remaining candidate genes and polymorphisms could lead to the identification of genes and polymorphisms responsible for FAR.