主页 › 范文 › 作物学 › 本文

植物育种中品系亲缘系数的修正计算

摘要

文章揭示了传统亲缘系数（COP）计算方法在植物育种应用中的系统性偏差。通过数学证明和育种实例分析，指出个体COP公式直接用于品系（含遗传变异的群体）会导致亲缘关系高估。提出基于等位基因频率的修正方法：建立终端祖先等位基因贡献矩阵，通过正交投影计算品系间真实COP值。在麦芽大麦等案例中，修正后COP值较传统方法降低37.5%，准确反映杂交后代近交系数。

1. 引言

亲缘系数（COP）作为遗传相关性度量，在作物改良计划中被广泛用于预测杂交表现。传统COP计算基于动物育种模型，假设谱系节点为单一基因型个体。文章首次指出：植物品系实质为遗传变异群体，直接套用个体COP公式将导致近交程度高估。以大麦育种为例，"Morex"与"Robust"品系杂交的预期近交系数实测值为1/4，而传统方法计算结果为1/2，偏差达100%。

2. 核心偏差机制

2.1 个体与品系的本质差异

个体节点：特定基因型，等位基因传递概率为1/2
品系节点：多基因型群体，需考虑等位基因频率分布

如图3谱系所示：当Z、W为个体时COP=1/4；作为品系时（Z含AC/BC基因型，W含AD/CD/BD），随机抽样个体间等位基因一致性概率降至3/16。

2.2 近交系数的误导性

传统公式依赖祖先近交系数F_X：

\[ f_{XZ} = \frac{1}{4}(1+F_X) \]

该式对品系失效。纯合品系（如含AA/BB基因型）虽F_X=1，仍可传递双等位基因，杂交后代近交系数实测值仅为1/2，而非公式所得1。

3. 品系COP修正方法

3.1 矩阵算法框架

构建谱系关联矩阵P（p_ij=1表示品系i是j的亲本）：

\[ \mathbf{Q} = (\mathbf{I} - \frac{1}{2}\mathbf{P})^{-1} \]

提取终端祖先贡献矩阵EQ，品系间修正COP计算为：

\[ f^{*}_{XY} = \frac{1}{2}(\mathbf{EQ})^T(\mathbf{EQ}) \]

3.2 麦芽大麦案例验证

应用图4谱系：

终端祖先X,Y,Manker基因型分别为AA,BB,CC
"Morex"（XY杂交后代）等位基因频率：A:1/2, B:1/2
"Robust"（Morex×Manker）频率：A:1/4,B:1/4,C:1/2

修正COP计算：

\[ f^{*}_{(Morex,Robust)} = (\frac{1}{2}\times\frac{1}{4}) + (\frac{1}{2}\times\frac{1}{4}) = \frac{1}{4} \]

较传统方法结果1/2降低50%。

4. 算法优势与意义

数学证明$f_{XY} \geq f^{*}_{XY}$恒成立，修正方法避免近交高估：

图6谱系中21个品系对，1/3存在显著COP差异
开放授粉品系自交COP实测1/2，传统公式得1
矩阵算法兼容复杂回交谱系，仅需满足哈迪-温伯格平衡

5. 结论

文章建立植物品系专属COP计算范式，通过终端祖先等位基因频率投影准确量化遗传相关性。修正方法为杂交设计提供可靠理论依据，避免传统方法导致的育种选择偏差。