模拟复合肉牛品种非名义后代风险的基因组架构分析

《Genetics Selection Evolution》：Genomic architecture of non-nominal offspring risk in a simulated composite cattle breed

　　

背景：复合肉牛品种旨在结合不同创始品种的互补特性，同时保留杂交优势的部分效应。婆罗格斯牛名义上定义为5/8安格斯和3/8婆罗门，但相同的全局品种组成可通过不同的基因组架构产生。这一区别可能很重要，因为具有相同婆罗门祖先比例的动物可能在祖先杂合度、纯合祖先组成以及作为亲本时的分离行为上存在差异。本研究的目的是利用受控模拟评估基因组架构是否会影响产生名义婆罗格斯牛祖先范围之外后代的概率。 结果：研究人员模拟了两种达到相同预期婆罗门祖先比例0.375的替代途径。尽管两种途径产生了相似的平均婆罗门祖先比例，但直接途径比间接途径产生了更宽的祖先分布，且位于名义区间之外的个体比例更高。将候选种公牛过滤至名义区间后，依赖于形成途径的基因组架构差异仍然明显。在一项扩展至100代和三种群体规模情景的混合互交分析中，途径特异性的祖先状态差异逐渐被侵蚀，婆罗门祖先比例的标准差在早期达到最大值后下降，但在模拟的时间范围内，非名义频率仍不可忽视。家系水平模拟显示，名义亲本之间的交配在其产生非名义女儿的概率上存在显著差异。这种风险取决于预期家系均值与目标祖先之间的距离，以及亲本基因组的区块加权分离潜力。一项探索性准二项式模型支持非名义风险与偏离目标的程度以及高方差亲本区块结构之间存在正相关关系。最后，一种启发式辅助的、基于模拟信息的种公牛分配策略相对于随机种公牛分配提高了名义女儿的比例，且从更广泛的种公牛目录中获得的增益最大。 结论：名义品种组成并不定义唯一的基因组状态。在模拟的婆罗格斯牛中，全局祖先等效的动物可能在基因组架构以及特定交配产生的后代祖先分布上存在差异。这些结果表明，复合品种管理不仅应考虑个体动物的祖先比例，还应考虑公母配对组合的预期结果。

研究背景

复合肉牛品种的设计初衷在于整合不同创始品种的互补优势，同时保留杂交优势的部分效应。例如，婆罗格斯牛名义上定义为5/8安格斯和3/8婆罗门，然而，相同的全局品种构成可以通过截然不同的基因组架构实现。以往的研究表明，即使两个体拥有相同的全球祖先比例，它们在祖先杂合度、纯合祖先组成以及作为亲本时的分离行为上也可能存在显著差异。因此，一个核心的管理问题是，仅凭个体的名义祖先比例是否足以预测其后代的遗传组成。当前研究的主要不足在于，大多数工作集中于描述个体或群体的祖先组成，而未能深入探讨不同基因组架构如何影响特定交配组合的后代祖先分布风险。为此，研究人员开展了此项受控模拟研究，旨在评估基因组架构是否会影响产生超出名义婆罗格斯牛祖先范围后代的概率，从而为复合品种的育种管理提供新的视角。

研究方法

研究人员使用AlphaSimR软件包进行模拟。首先，创建了两个完全分化的固定创始品种，分别代表安格斯（基因型0/0）和婆罗门（基因型1/1）。随后设计了两种达到相同预期婆罗门祖先比例（0.375）的形成途径：直接途径和间接途径。研究中计算了婆罗门祖先比例、基因组三态组成（安格斯纯合、杂合、婆罗门纯合）、预期家系均值、以及基于亲本联合状态区块的区块加权分离潜力指数等指标。通过家系水平模拟和混合互交分析，并结合探索性准二项式模型，最终提出了一种启发式辅助的、基于模拟信息的种公牛分配策略。

研究结果

分布特征与基因组异质性

直接途径产生了比间接途径更宽的祖先分布和更高的标准差（SD=0.040 vs. 0.026），导致更大比例的个体偏离名义区间。间接途径则表现出显著更高的祖先杂合度（0.746 vs. 0.507）。这表明，即使平均祖先比例相同，两种途径生成的群体在个体水平和基因组架构上存在本质差异。

名义种公牛的基因组约束

将候选种公牛筛选至名义祖先区间后，路径依赖性的基因组架构差异依然存在。直接途径的种公牛分布在较宽的三态组成空间内，而间接途径的种公牛则集中在高杂合度区域，缺乏婆罗门纯合位点。这说明名义祖先过滤并不能消除形成途径所施加的结构性差异。

商业牛群动态

使用名义种公牛更新商业母牛群10代后，虽然群体平均祖先比例保持稳定，但非名义个体比例在各场景中仍居高不下（43.4%-49.6%）。特别是，间接途径形成的母牛群在更新过程中其高杂合度架构迅速被稀释。这表明仅靠使用名义种公牛不足以维持整个群体处于名义区间内。

混合互交分析

在长达100代的混合互交分析中，无论群体规模大小，婆罗门祖先比例的标准差均呈现先快速上升后缓慢下降的趋势，但到第100代时仍高于独立位点的理论期望值。非名义后代的频率在早期达到峰值后虽有所下降，但在所有规模场景下仍维持在13%-18%的不可忽视水平，表明该风险是长期存在的。

家系水平风险分解

在家系水平分析中，所有亲本均为名义个体，但不同途径组合的交配场景显示出显著不同的非名义后代风险：直接×直接最低（38.1%），间接×间接最高（51.9%）。研究发现，这种风险主要由两个因素决定：一是预期家系均值与目标祖先的距离，二是由亲本基因组中高变异区块（尤其是双杂合子区块）的空间分布所决定的区块加权分离潜力。大的高变异区块对后代方差有不成比例的巨大贡献。

种公牛分配优化

在概念验证分析中，研究人员提出的启发式辅助模拟信息种公牛分配策略，相较于随机分配，在所有三种种公牛目录（严格、中等、宽松）下均能提高名义女儿的比例。特别是在最宽松的种公牛目录下，该策略带来的增益最为显著，证明了利用配对特异性基因组信息进行管理决策的潜在价值。

讨论与结论

本研究的结论明确指出，名义品种组成并不定义一个唯一的基因组状态。在模拟的婆罗格斯牛中，全球祖先等效的动物在基因组架构以及特定交配产生的后代祖先分布上可以存在显著差异。这些发现有力地证明，复合品种的管理不应仅仅关注单个动物的祖先比例，而应将其扩展到对特定公母配对组合预期结果的评估。研究人员认为，将公母配对视为控制品种组成和降低非名义后代风险的关键管理单元，将为未来的复合品种育种实践提供重要的理论依据和实用工具。