《关于用育种值提升一般配合力分析水平的研究》原文发表于《种子科技》2019年第01期，以下是摘录的部分内容，供育种者参考。

动物育种中育种值（Breedingvalue）广泛使用，甚至是动物育种学重要的教学内容之一，而植物育种尤其玉米育种，多采用配合力理论指导实际育种工作，育种值和配合力都在育种实践中取得了很好的育种结果。无论植物育种还是动物育种，都是对基因的选择，在解决多年多点配合力分析问题时，笔者发现一般配合力和育种值在模型构建上一致，进而提出了适用于植物育种使用的植物育种值模型来提升一般配合力的计算分析水平。

1、多年多点配合力分析问题

对一个群体进行过一次以上的配合力试验，如异地安排试验进行穿梭育种，如何解决多年多点配合力试验数据分析问题？很显然，多年多点配合力试验可以提供更丰富更准确的试验信息，这对亲本选择意义重大。

依据线性模型的效应可加性，多年多点配合力试验，可以分离的效应中包含年度效应、地点效应、一般配合力效应（GeneralCombining Ability, GCA）、特殊配合力效应（Special CombiningAbility, SCA）。单一年度单一地点的配合力试验是多年多点配合力试验的特例。在单一的配合力试验中，年度效应、地点效应因为无法有效分离，混杂在其余效应中。线性混合模型中固定效应可多个，随机效应多设为一个，为保持跟一般配合力定义一致，将模型中的GCA构建为随机效应（Σ效应=0）。在育种实践中，一个配制好的杂交组合，不种出来，无法知道其真实生产表现，一般认为这与SCA有关，因此玉米育种多采用GCA选择亲本自交系，采用SCA寻找优良组合，这也意味着SCA在育种选择方面事实上失去了计算意义。因而模型应侧重GCA，于是用年度效应、地点效应取代模型中的SCA。

得到了多年多点配合力分析模型：

Y =  μ + 年度效应 + 地点效应 + GCA + 误差

其中：年度效应、地点效应均为固定效应，GCA为随机效应。

跟动物育种中的动物模型进行比较，动物模型^[1]：

Y= μ +  + a+ e

模型中 b_j---第j个系统环境效应，为固定效应

a---加性遗传效应（育种值），为随机效应

e---随机残差

多年多点配合力分析模型中年度效应和地点效应均属于环境效应，于是，植物多年多点配合力模型和育种值动物模型两者模型一致。两个模型的区别仅在于同为随机效应的，一个是一般配合力，一个是育种值。

植物领域一般配合力使用，一直依据其定义（g_i =y_i. - y_..）进行计算，90年代，美国开始在植物育种上尝试采用BLUP分析方法，所用分析模型，包含了GCA和SCA效应，没有包含年度和地点效应。模型效应线性可加的情况下，应力求效应的生物学意义明晰，因此年度和地点效应有必要在模型中单独列出进行准确分离，这符合植物开放式生产环境的特点，也与生物统计方差分析的效应分离思路一致。考虑一般配合力（基因加性效应）对育种选择的重要作用，笔者提出了以育种值为随机效应的线性混合效应模型---植物育种值模型。

2、植物育种值模型

Y = μ + b_年度 + b_地点 +  a_j + e

Y：数据观测值，

μ：总体均值

b_年度：试验年度的效应,固定效应

b_地点：试验地点的效应,固定效应

a_i：第j个品种的育种值，随机效应

e：随机残差

考虑种群效应的育种值模型：

Y = μ + b_年度 + b_地点 + h_i + a_ij + e

Y：数据观测值，

μ：总体均值

b_年度：试验年度的效应,固定效应

b_地点：试验地点的效应,固定效应

h_i：第i个种群效应，固定效应

a_ij：第i个种群第j个品种的育种值，随机效应

e：随机残差

动物育种发展出动物模型、公畜模型、公畜母畜模型、外祖父模型等，植物育种可以套用动物模型。但是考虑植物育种的特点，植物育种值模型更能反映植物育种数据特点。每个观测值可以从效应的角度，剖分为总体均值，不同年度效应，不同地点效应，不同种群效应，品种的育种值效应。对照数量遗传学的基本模型，可以发现BLUP线性混合模型特别适合计算育种值。加性遗传模型下，G=M+ A，M为群体均值，A为加性效应，对应育种值。

表现型值P = 基因型值G + 环境效应E + 随机误差e

基因型值G=群体均值M+ 加性效应A

3、模型使用分析探讨

3.1育种值较配合力更符合数量遗传学模型

《植物育种的数量遗传学基础》在介绍1942年Sprague 和Tatum提出的配合力育种时，是这样表述的：“所谓‘一般配合力’指一个自交系或品种（纯合体）在一系列的杂交组合中的平均产量（或其他经济性状）的表现；所谓‘特殊配合力’指根据于一般配合力所预测的平均产量（或其他经济性状）的离差，或者说，和所有杂交组合的平均数比较而论，某一的杂交组合中所表现的产量（或其他经济性状）较之平均数为优或劣的结果”。孔繁玲《植物数量遗传学》则直接用“配合效应”来介绍配合力：“通常把亲本在实际组配中的表现称为该亲本或该对亲本的配合力。一个亲本在其所有杂交组合中表现的平均配合效应称之为该亲本的一般配合力效应（GCA）。一般配合力反映特定亲本交配效应的平均水平。”

配合力是从计算的角度而不是从效应的角度，给出的定义，找不到跟基因效应的直接字面联系，跟数量遗传学的效应模型也没有建立相应关联，这就造成了实际育种工作中对配合力的不同理解甚至滥用。依据一般配合力公式，获得的结果接近于亲代稳定遗传给子代的效应部分，亲代可稳定遗传效应部分，跟亲代基因有关，而加性效应是数量遗传学理论中可稳定遗传给子代的基因效应，就是育种值。数量遗传学对育种值分别给出了理论定义和实际定义，反映的就是要将数量遗传学效应模型和实际育种实践相结合的需要。

在植物分子育种领域，配合力概念*被育种值取代，因为直接面对基因的各种评价或者选择，就是各种基因效应，而跟育种选择密切有关的就是加性效应就是育种值。

3.2育种值模型对试验设计无特殊要求

配合力分析需要以配合力试验设计为基础，但是育种值模型却不受此限制。育种值模型可以用于分析专业的配合力试验，另外，无需试验设计，依据材料系谱，仍然可以进行育种材料的育种值评价。

3.3育种值模型在计算方法上优于配合力模型

相比一般配合力，结合材料间亲缘关系的BLUP育种值模型可以更准确计算材料的育种价值。

一般配合力的计算基于配合力的定义（g_i= y_i. - y_..），没有分离年度效应、地点效应，计算简单，不能分析多年多点配合力数据。育种值模型，将育种值设定为随机效应，其余，地点、年度、种群设为固定效应，更符合植物育种实际情况。充分考虑材料间的亲缘关系构建出了亲缘关系A矩阵，A矩阵相比配合力试验的父母本组配设计，可以更全面更准确反映材料亲缘关系，采用线性混合模型对育种值求解，从而可以更准确计算基因加性效应。

借助育种值模型，可以做到数据积累越多，对育种材料的育种价值把握越准确，这对育种工作尤为重要。

3.4建议用育种值提升一般配合力分析水平而不是*替代或否定配合力

从育种值模型看，模型主要改进的是一般配合力的分析水平，不涉及特殊配合力，因此主张用育种值提升一般配合力分析水平，而不是去*替代配合力，更不是否定配合力。

以玉米育种为例，当下育种者熟悉的育种路线是前期依据测验种决定的一般配合力选育自交系，后期利用za种优势，发现有特殊配合力效应的材料组合。育种值模型使用，不改变任何育种安排，仅是对测配试验，用育种值模型分析方法替代原有的一般配合力分析方法，来进行试验数据分析。

张勤详细介绍了BLUP相关知识，张沅给出了动物育种值具体计算示例，配合力试验数据直接用育种值模型分析及结果比对可以参照《育种实用统计学》。

......

农博士育种家软件提供BLUP育种值分析评价功能。