KIRREL3内含子多态性与沃金牛胴体性状、肉品质及肌肉脂肪酸组成的关联研究

《Animals》：Association of KIRREL3 Intronic Polymorphisms with Carcass Traits, Meat Quality, and Muscle Fatty Acid Composition in Woking Cattle

　　

肌内脂肪（IMF）沉积是决定高端牛肉价值的关键因素，但表型选择成本高昂且周期漫长。尽管全基因组研究已将KIRREL3区域与肉品质相关联，但内含子变异是否直接影响胴体性状并在不同生理阶段保持预测效力仍未得到证实。研究人员对中国吉林省培育的高品质肉牛品种沃金牛（n = 473；其中36月龄和48月龄母牛分别为226头和206头，另设独立验证群n = 41）的两个内含子单核苷酸多态性（SNP）位点（g.29756010C>T和g.29756149C>T）进行基因分型，并评估了脂肪酸组成、肉品质及多组织基因表达。两个位点均具有多态性且符合哈迪-温伯格平衡。在36月龄时，两个位点的TT个体在屠宰率、肋厚和肉品等级方面均高于CC个体。至48月龄时，63-bp位点与胴体性状的关联范围更广，而203-bp位点的效应减弱。在验证群中，CC个体的IMF含量和肌肉KIRREL3表达量高于TT个体，但肉品等级却较低。经Benjamini-Hochberg校正后，有5种脂肪酸在基因型间存在显著差异（q < 0.05），而主要脂肪酸（C16:0、C18:1n9c、C18:0、C18:2n6c）未见显著比例差异。在检测的8种组织中，KIRREL3在肾脏中的表达量最高。这些发现确定63-bp位点为中国地方黄牛肉品质的候选标记，具有标记辅助选择的潜力，但需在更大规模的群体和功能研究中进一步验证。

沃金牛（Woking cattle）KIRREL3基因内含子多态性与胴体性状、肉品质及肌肉脂肪酸组成的关联研究——论文解读

一、研究背景与科学问题

肌内脂肪（Intramuscular fat, IMF）沉积是决定牛肉市场价值的最重要因素，直接影响牛肉的嫩度、多汁性和风味。然而，IMF受多基因控制且对环境高度敏感，传统的表型选择需要将动物饲养至屠宰月龄。在中国优质牛肉群体中，这一周期通常超过36个月，生产者缺乏对大理石花纹潜力的早期遗传筛选手段。这种在数年昂贵饲喂之前无法识别优良个体的局限性，构成了肉牛遗传改良中的重大经济制约。虽然安格斯牛（Angus）和和牛（Wagyu）等品种通过数十年的选择育种建立了大理石花纹基准，但包括沃金牛在内的中国地方品种在分子水平上仍 largely 未被探索， validated 的候选基因标记十分稀缺。

沃金牛是吉林省培育的高品质肉牛品种，以大理石花纹牛肉著称，达到市场级大理石花纹通常需要超过36个月的育肥期。然而，肉品质只能在屠宰后才能评估，生产者往往在可能最终评级不佳的动物上投入多年的饲料和劳动力。候选基因关联分析为实施大理石花纹潜力的标记辅助选择提供了切实可行的第一步。

KIRREL3（kin of IRRE-like 3, NEPH2）编码免疫球蛋白超家族的嗜同性细胞粘附分子，最初因其在神经元突触形成中的作用而被表征。在功能层面，KIRREL3与Na⁺,K⁺-ATP酶α3亚基共表达，提示其参与细胞能量代谢和信号转导。值得注意的是，包含KIRREL3的BTA29区间已在家牛肉牛和日本和牛中鉴定出选择信号，但这些基因组扫描发现尚未转化为经过验证的标记-性状关联。此外，尚无研究检验KIRREL3标记效应是否在育肥期内持续存在——大多数候选基因研究依赖单一时间点，而IMF持续积累且代谢通路活性在36月龄左右发生显著变化，同一标记在36月龄和48月龄可能表现不同。

二、研究目的与研究内容

本研究旨在通过整合多层面分析，系统评估KIRREL3基因两个内含子SNP（g.29756010C>T和g.29756149C>T）与沃金牛胴体性状、肉品质和肌肉脂肪酸组成的关联。研究人员采用直接测序法对沃金牛进行基因分型，在36月龄和48月龄两个时间点开展关联分析，并纳入独立验证群、肌肉脂肪酸谱分析和多组织表达分析，以期鉴定出足以支持早期辅助选择的稳健标记。

三、关键技术方法

本研究采用了以下主要技术方法：（1）PCR扩增与Sanger测序：针对KIRREL3基因第一内含子设计特异性引物，扩增603 bp片段覆盖两个SNP位点，通过双向测序进行基因分型；（2）脂肪酸气相色谱-质谱联用（GC-MS/MS）分析：采用酸水解法提取肌肉脂质，甲酯化后经DB-FastFame毛细管柱分离，以37组分脂肪酸甲酯标准品进行外标定性和峰识别；（3）实时荧光定量PCR（RT-qPCR）：以β-actin为内参基因，采用2^?ΔΔCt法计算KIRREL3相对表达量，以心脏组织为校准组织；（4）群体遗传学分析：计算基因型频率、等位基因频率、期望杂合度（He）、多态信息含量（PIC）和有效等位基因数（Ne），采用皮尔逊卡方检验评估哈迪-温伯格平衡；（5）生物信息学虚拟敲除（in silico virtual knockout）预测：整合群体基因组选择信号、功能注释和网络扰动预测，模拟KIRREL3基因敲除对成肌和成脂调控网络的影响。

样本来源方面，主要关联群体为吉林省长春皓月有限公司（Changchun Haoyue Co., Ltd.）的432头健康沃金母牛，其中36月龄226头（体重798.71 ± 28.93 kg），48月龄206头（体重818.66 ± 30.52 kg）。独立验证群从986头已知KIRREL3基因型的群体中筛选41头健康母牛，按三种基因型（CC、TC、TT）分层平衡选择。

四、研究结果

（一）测序与多态性分析结果

直接测序603 bp KIRREL3内含子片段揭示了两个SNP位点：g.29756010C>T（63-bp位点）和g.29756149C>T（203-bp位点）。两个位点均呈中等多态性（0.25 < PIC < 0.50）。在63-bp位点，CT杂合子频率最高（51.77%），其次为CC（29.20%）和TT（19.03%）；等位基因频率C为0.5509，T为0.4491。在203-bp位点，TC同样占主导（52.94%），其次为CC（31.86%）和TT（15.20%）。两个位点均未偏离哈迪-温伯格平衡（63-bp: χ² = 0.48, p = 0.487; 203-bp: χ² = 1.62, p = 0.203），表明群体遗传结构稳定。

（二）36月龄胴体性状关联分析

在63-bp位点，基因型显著影响屠宰率、肋厚、肌间脂肪厚度和肉品等级（均p < 0.05，单因素方差分析）。Tukey HSD事后检验显示，TT个体在屠宰率、肌间脂肪厚度和肉品等级方面显著高于CC个体；肋厚方面TT显著厚于CC，TC居中且与两者均无显著差异。在203-bp位点，TT在屠宰率、肋厚和肉品等级方面再次显著高于CC，但肌间脂肪厚度在该位点无显著关联。

（三）48月龄胴体性状关联分析

至48月龄，63-bp位点对屠宰率、肉品等级、肋厚和皮下脂肪厚度均显示显著基因型效应。事后分析表明TC和TT在肉品等级和肋厚方面均显著高于CC（均p < 0.01），TT在屠宰率方面仍显著高于CC（p < 0.05）。值得注意的是，皮下脂肪厚度在TC中最高，显著高于TT（p < 0.05），CC居中，提示存在杂合效应。在203-bp位点，TT维持较高的屠宰率（61.68 ± 1.77%）和肋厚（7.81 ± 0.66 cm），但关联范围明显收窄。63-bp位点在48月龄时影响范围扩大，而203-bp位点仅保留有限关联。

（四）脂肪酸组成分析

共检测到33种脂肪酸，其中32种在全部三组基因型中可定量。总体脂肪酸谱以油酸（C18:1n9c，约22–25%）、棕榈酸（C16:0，约17–18%）、硬脂酸（C18:0，约15–16%）和亚油酸（C18:2n6c，约10–11%）为主。经Benjamini-Hochberg（BH）多重检验校正后，5种脂肪酸在基因型间差异显著（q < 0.05）：顺-10-十五碳烯酸（C15:1, q = 0.0194）、γ-亚麻酸（C18:3n6, q = 0.0282）、顺-11,14,17-二十碳三烯酸（C20:3n3, q = 0.0008）、芥酸（C22:1, q = 0.0067）和神经酸（C24:1, q = 0.0000）。所有显著差异均涉及低丰度脂肪酸（占总鉴定FA的<1%），主要脂肪酸（C16:0、C18:0、C18:1n9c、C18:2n6c）在基因型间无显著差异（均q > 0.05）。

（五）验证群肉品质性状分析

验证群41头母牛按基因型分层（CC, n = 12; TC, n = 20; TT, n = 9）进行肉品质分析。IMF含量在基因型间差异显著（p = 0.0267），CC显著高于TT（p < 0.05），TC居中。CC的L\*和b\*值显著高于TT（分别p < 0.001和p < 0.05）。持水力方面，TT的离心损失较高（p < 0.05），而CC的蒸煮损失较低（p < 0.01 vs. TC）。TT的粗蛋白含量显著高于CC（p < 0.01）。加压损失、滴水损失、嫩度和a\*值无显著差异。

（六）草原红牛和延边牛的多态性验证

为明确KIRREL3标记在中国东北牛群体中的遗传背景，研究人员在草原红牛（n = 50）和延边牛（n = 50）中对两个SNP进行基因分型。两个品种均表现中等多态性（PIC = 0.37），有效等位基因数接近2，符合哈迪-温伯格平衡。品种间基因型分布（χ² = 2.496, p = 0.287）和等位基因频率（χ² = 2.010, p = 0.156）均无显著差异，遗传分化系数F_ST为0.0122，表明分化程度极弱。

（七）KIRREL3组织表达谱

RT-qPCR检测发现KIRREL3 mRNA在全部8种组织中均有表达。肾脏表达量最高（为心脏的36.25倍，组a），其次为背脂、肝脏、最长肌和腹脂（11.67至17.13倍，组b），均显著高于心脏（p < 0.05）。脾脏和肺脏表达量分别为心脏的2.98倍和3.37倍，与心脏同属最低表达层（组c）。

（八）肌肉中KIRREL3基因型特异性表达

在验证群肌肉样本中，CC基因型的KIRREL3相对表达量显著高于TC和TT（p < 0.05），TC与TT之间无显著差异。值得注意的是，CC表现出最高的肌肉表达量却对应最低的肉品等级，而TT表达量最低却获得最高等级，呈现反向模式，提示KIRREL3表达与大理石花纹评分之间存在复杂关系。

（九）虚拟敲除预测

虚拟敲除预测表明，KIRREL3缺失将导致核心成肌因子（MYOD、MYOG、MEF2C）下调和脂肪生成因子（PPARG、FABP4、FASN、SCD）及细胞外基质重塑因子（MMP2）上调。在"粘附屏障假说"框架下，KIRREL3缺失将松弛成肌细胞间的粘附，允许前脂肪细胞浸润至肌纤维间空隙，同时阻断成肌细胞融合，导致肌纤维变细和肌肉密度增加。尽管总IMF含量预计下降，但脂肪细胞在可见肌束间区域的空间再分布预计将提高客观大理石花纹评分。

五、讨论与结论

本研究表明，KIRREL3基因两个内含子SNP与沃金牛胴体品质相关，细化了此前通过全基因组重测序检测到的BTA29选择信号。63-bp位点在36月龄和48月龄均与屠宰率、肋厚和肉品等级显著关联，且48月龄时呈现T等位基因完全显性模式（TC ≈ TT > CC），这对于育种实践具有直接优势——单个T等位基因即可提供完整的表型效益，育种者可以保留杂合性并避免近交衰退。值得注意的是，两个位点均不影响活重、胴体重和眼肌面积，表明选择T等位基因可在不牺牲生长性能的前提下改善肉品质。

CC基因型动物积累更多IMF并表现出更高的肌肉KIRREL3表达，却获得较低的肉品等级，这一看似矛盾的现象可通过将大理石花纹视为空间过程而非单纯脂质储存来解释。研究提出KIRREL3通过结构性"粘附屏障"机制影响大理石花纹可见性：高表达封闭成肌细胞膜并限制前脂肪细胞浸润，产生高化学IMF但可见大理石花纹差；低表达（或虚拟敲除）则松弛该屏障，允许脂肪细胞再分布并改善肉品等级。然而，这些预测仍属计算性质，需通过CRISPR-Cas9敲除和等位基因特异性表达分析进行功能验证。

研究结论：在沃金牛中鉴定出两个KIRREL3内含子SNP（g.29756010C>T和g.29756149C>T），两个位点均呈中等多态性且符合哈迪-温伯格平衡。63-bp位点在36月龄和48月龄均与屠宰率、肋厚和肉品等级相关，且遗传模式随时间从TT优越性转变为T显性。203-bp位点在36月龄显示相似关联但至48月龄减弱。在验证群中，CC动物积累更多IMF并表现出更高的肌肉KIRREL3表达，但肉品等级较低。经BH校正后5种低丰度脂肪酸在基因型间差异显著，主要脂肪酸无显著差异。KIRREL3在肾脏中表达量最高。这些发现支持63-bp SNP作为中国东北地方牛种肉品质的实用候选标记，鉴于其跨阶段持续性和在草原红牛及延边牛中的存在性。但本研究确立的是统计关联而非因果关系，内含子变异可能标记连锁的调控位点而非直接发挥效应。验证样本量较小（n = 41），且群体在不同年份的商业条件下饲养。在将该标记整合到标记辅助选择程序之前，必须在更大规模的多农场群体中进行确认，并结合荧光素酶报告基因检测、EMSA和肌肉中KIRREL3表达的单细胞分辨率等功能实验。