鉴定与日本黑牛BTA17染色体脂肪面积比相关候选因果变异：GGT1与SLC5A1基因的功能关联

《Genes》：Identification of Candidate Causal Polymorphisms in GGT1 and SLC5A1 Associated with Fat Area Ratio on BTA17 in Japanese Black Cattle Shinji Sasazaki, Hikari Ito, Ryoto Adachi, Eiji Iwamoto, Emi Yoshida, Fuki Kawaguchi, Kenji Oyama and Hideyuki Mannen

本研究聚焦于和牛大理石纹肉品质遗传改良的关键问题，通过整合全基因组关联研究(GWAS)、高密度连锁不平衡(LD)精细定位与功能注释分析，在兵库县日本黑牛群体中成功鉴定了位于BTA17染色体上与肋眼肌脂肪面积比(FAR)显著相关的候选因果多态性，特别是GGT1基因的错义突变(c.589G>A)和5′UTR变异，以及SLC5A1基因的错义突变。该研究为深入理解肌肉脂肪沉积的分子机制提供了新线索，并为该品种大理石纹性状的精准分子育种提供了潜在的有效标记。

　　

大理石纹，即肌肉脂肪，是决定牛肉口感、多汁性和风味的关键因素。尤其在日本黑牛（一种重要的和牛品种）中，卓越的肌肉脂肪沉积能力是其高品质牛肉的标志，经过数十年对肉质性状的强化选育，形成了独特的遗传背景。然而，随着消费者健康意识的提升和口味的多样化，市场对脂肪沉积的精准遗传调控提出了更高要求——不再是单纯增加大理石纹，而是将其调整至适中的理想水平。当前，尽管已知大理石纹受遗传控制且遗传力较高，但控制该性状的主要致病基因和功能性多态性在很大程度上仍是未知的，这为实现精准的遗传改良带来了挑战。牛肉大理石纹标准(BMS)是传统的视觉评分体系，而肋眼肌脂肪面积比(FAR)作为一种通过图像分析得出的定量指标，与BMS评分高度相关，为遗传研究提供了更客观、可量化的表型。在此背景下，本研究旨在利用先进的基因组学方法，在特定日本黑牛群体中深入挖掘影响FAR的关键遗传变异。

本研究主要采用了以下几个关键技术方法：1. 样本与表型：研究使用了432头来自日本兵库县的日本黑牛（343头阉牛，89头母牛），通过图像分析精确测量了其肋眼肌的脂肪面积比(FAR)。2. 全基因组关联研究(GWAS)：使用Illumina BovineSNP50v3芯片对样本进行基因分型，通过线性混合模型校正固定效应和随机效应后，利用残差进行关联分析，鉴定与FAR显著相关的基因座。3. 精细定位与变异筛选：对42头兵库县日本黑牛的全基因组重测序数据进行分析，在GWAS信号区域（BTA17: 71-74 Mbp）内构建高密度连锁不平衡(LD)图谱，并结合变异功能预测（如错义变异、UTR变异等），从与top SNP存在LD(r²≥ 0.1)的基因内变异中筛选出96个候选多态性。4. 候选变异验证：在原始群体中对候选变异进行基因分型，并通过方差分析(ANOVA)评估其与FAR的关联性。

3. 结果与讨论

• •
  GWAS鉴定出BTA17上的显著位点：GWAS分析在BTA17染色体上检测到一个与FAR相关的显著信号，其top SNP位于17:72,329,662 (p = 3.60 × 10^-6)。该数量性状位点(QTL)可能反映了兵库县日本黑牛群体因封闭育种而形成的独特遗传背景。
• •
  高密度LD分析界定候选区域：利用全基因组重测序数据对top SNP周围进行高密度LD分析，发现了一个明确的LD区块，范围在71-74 Mbp之间，从而将该区域界定为候选区间。
• •
  筛选与验证候选功能性变异：在候选区域内，从与top SNP存在LD(r²≥ 0.1)的基因内4292个变异中，依据LD强度和功能预测注释筛选出96个候选变异进行基因分型验证。关联分析显示，位于γ-谷氨酰转移酶1(GGT1)基因的一个错义突变c.589G>A (p.Asp197Asn)与FAR的关联最为显著(p = 2.86 × 10^-6)。同时，GGT1基因的一个5′UTR变异(c. -256G>T)也显示出类似的强关联性。这两个变异与top SNP均存在强LD(r²= 0.748)。
• •
  发现另一个生物学合理候选基因SLC5A1：基于基因功能分析，溶质载体家族5成员1(SLC5A1)被确定为另一个有生物学意义的候选基因。该基因编码钠-葡萄糖协同转运蛋白1(SGLT1)，与营养感知和代谢调节有关。研究发现该基因的一个错义突变c.32C>T (p.Thr11Met)与FAR存在显著关联(p = 2.03 × 10^-5)，且与top SNP的LD值很高(r²= 0.775)。

4. 结论

本研究通过GWAS结合基于全基因组测序的精细定位，在兵库县日本黑牛中成功鉴定出BTA17染色体上一个与肋眼肌脂肪面积比(FAR)相关的新QTL。该研究进一步将候选范围缩小至71-74 Mbp的LD区块，并通过功能筛选和验证，首次将GGT1基因的错义突变(c.589G>A)和5′UTR变异(c. -256G>T)以及SLC5A1基因的错义突变(c.32C>T)确定为与该性状相关的强候选因果多态性。GGT1参与谷胱甘肽代谢和氧化应激调节，可能通过影响胰岛素敏感性来调控脂肪代谢；SLC5A1则涉及葡萄糖吸收和代谢调节，可能与脂肪细胞分化相关。尽管这些变异的功能影响尚需后续实验验证，但它们为理解肌肉脂肪沉积的遗传机制提供了重要的位置和生物学候选，并有望成为用于兵库县日本黑牛大理石纹相关性状遗传改良的实用分子标记。这项研究展示了一种高效整合多种基因组学数据鉴定复杂性状候选因果变异的方法，为家畜重要经济性状的精准育种提供了有价值的案例。