黑白花兔繁殖性状的相关及通径分析

陈岩锋，陈冬金，孙世坤，桑 雷，王锦祥，谢喜平，许建书

（1.福建省农业科学院畜牧兽医研究所，福州 350013；2.龙岩市万家兔业发展有限公司，福建 龙岩 364100）

黑白花兔是福建省龙岩市养殖的地方兔群中出现的毛色黑白花个体，其被毛黑白花，头型清秀、眼珠黑色、耳短小、颈稍短、无肉髯、结构匀称、肌肉紧凑结实，属于小型肉用兔，当地群众认为其肉质鲜美，深受消费者青睐［1］。繁殖性状是家兔的重要经济性状之一，其性能高低不仅影响母兔的繁殖力，而且直接影响养兔生产的经济效益。繁殖性状主要包括产仔数、产活仔数、初生窝重、21日龄窝重（泌乳力）、断奶仔数、断奶窝重等。本文根据龙岩市万家兔业发展有限公司种兔场饲养的黑白花兔繁育群母兔繁殖记录资料，对黑白花兔繁殖性状的多个指标进行了相关分析、回归分析和通径分析，研究了有关繁殖性状之间的相互关系，剖析了影响断奶窝重的直接效应和间接效应，以便为黑白花兔繁殖性状的选育工作提供理论基础。

1 材料与方法

1.1 材料

试验所用数据材料来自龙岩市万家兔业发展有限公司种兔场黑白花兔繁育群母兔2016年4月至9月繁殖记录共计194窝。繁殖性状包括产仔数（X1）、产活仔数（X2）、初生窝重（X3）、21 日龄只数（X4）、21 日龄窝重（泌乳力）（X5）、28日龄断奶只数（X6）和 28日龄断奶窝重（Y）。

1.2 饲养管理

公兔4.5～5月龄、母兔4～4.5月龄开始配种，采用自然发情、人工辅助配种，第2胎以后在产后18 d配种，夏季（7—8月份）停配1个月，避开全同胞、半同胞进行随机交配；母兔每天饲喂2次全价颗粒饲料，其日粮营养水平为 DE 10.8～11.3 MJ/kg，CP 16.5%～17.5%，CF 12%～16%，Ca 0.5%～1.5%，P 0.4%～0.7%。自由饮水；仔兔断奶前实行母仔分养、定时哺乳，18日龄开始补料，按兔舍常规饲养管理程序进行。

1.3 统计分析方法

利用Excel数据分析工具计算出各繁殖性状的平均数、标准差、变异系数及繁殖性状间的简单相关系数，并以产仔数（X1）、产活仔数（X2）、初生窝重（X3）、21日龄只数（X4）、21日龄窝重（泌乳力）（X5）、28日龄只数（X6）为自变量，28日龄断奶窝重（Y）为依变量计算出偏回归系数，再运用有关数量遗传学理论及通径系数原理对其进行通径分析［2］。

2 结果与分析

2.1 繁殖性状的表型参数

黑白花兔繁殖性状的平均数、标准差、变异系数统计结果见表1，其平均产仔数、产活仔数、断奶仔数分别为7.15、6.76和5.70只，28日龄断奶窝重为2 770.1 g，表明其繁殖性能较高；21日龄窝重变异系数最大，为31.41%，28日龄断奶窝重次之，为20.71%，表明这些性状具有较大的选育提高潜力。

表1黑白花兔繁殖性能表型参数

2.2 繁殖性状间的简单相关系数

由表2可知，两繁殖性状之间的简单相关均呈正相关，且表型相关系数均达极显著水平（P＜0.01）。21日龄只数与28日龄只数的相关系数最大，为0.876 9，其次是产仔数与产活仔数的相关系数，为0.832 2。

表2繁殖性状间的相关系数

2.3 多元回归分析

以产仔数（X1）、产活仔数（X2）、初生窝重（X3）、21日龄只数（X4）、21日龄窝重（泌乳力）（X5）、28日龄断奶只数（X6）为自变量，28日龄断奶窝重（Y）为依变量，运用Excel数据分析工具进行繁殖性状间的回归分析，建立回归方程：Y=343.295 6+1.350 3X1-61.411X2+1.192 2X3+1.989 0X4+0.385 8X5+301.185 2X6，总回归F=63.4，回归关系达极显著水平（P＜0.01），而其复相关系数R2=0.6705，回归方程的误差e对y也产生作用，其大小以剩余通径系数Pye表示。剩余通径系数Pye=0.574 0，其值较大，说明还有一些对28日龄断奶窝重有影响的因素没有被考虑进去，对28日龄断奶窝重影响因素的全面分析还有待于进一步深入研究。

2.4 繁殖性状间的通径分析

通径系数为标准化的偏回归系数，其计算公式为：Piy＝bi×Sxi/Sy，bi为偏回归系数、Sxi为Xi的标准差、Sy为Y的标准差，根据以上偏回归系数和标准差计算出各自变量X对依变量Y的通径系数。

通径分析结果（表3）表明，对28日龄断奶窝重（Y）直接影响最大的是28日龄断奶只数（X6），其通径系数为0.548 5，其次是 21日龄窝重（泌乳力）（X5），其通径系数为 0.373 0；产仔数（X1）和 21 日龄只数（X4）对断奶窝重的直接影响最小（通径系数分别为0.003 0和0.003 9），但通过断奶只数（X6）对断奶窝重（Y）产生较大的间接影响，其间接效应分别为0.3814和0.4810。可见，产仔数、21日龄只数对断奶窝重的影响主要是通过断奶只数来实现的。

表3繁殖性状间的直接效应和间接效应

2.5 决定系数分析

在实践中，除了研究各个变量（Xi）对依变量（Y）的相对影响程度外，还要研究各个变量（Xi）对依变量（Y）的绝对影响程度，绝对影响程度用决定系数表示。它是通径系数的平方（直接决定效应）或两个相关原因的通径系数和相关系数乘积的两倍（间接决定效应）。计算公式为：dy.xi=p2y.xidyxi.xj=2rxi.xjPy.xiPy.xj，式中dy.xi为Xi对Y的决定系数，dyxi.xj为Xi与Xj性状对Y的决定系数，其计算结果（表4）表明：对断奶窝重（Y）的影响最大的是断奶只数（X6）（0.300 85），其次是 21 日龄窝重（X5）与断奶只数（X6）的互作效应（0.230 37）、21 日龄窝重（X5）（0.139 13）。决定系数总和∑d＝0.670 6，误差项为1-0.670 6=0.329 4，表明影响断奶窝重的主要因素都已分析到，但还可能存在观测值误差或者可能还有对断乳窝重有影响的性状如胎次、季节等因素未被考虑到，有待进一步研究。

表4各繁殖性状的决定系数剖析

3 讨论与结论

从繁殖性状表型参数可见，黑白花兔21日龄窝重变异系数较大，为31.41%，28日龄断奶窝重次之，为20.71%，繁殖性能属于低遗传力性状，受环境因素影响较大。黑白花兔各两繁殖性状间的简单相关均呈正相关，且表型相关系数均达极显著水平。其中28日龄断奶窝重与28日龄断奶只数、21日龄只数、21日龄窝重密切相关，相关系数分别为0.748 2、0.700 5和0.687 6，该结果与桑雷等［3］和马旭平等［4-5］分别对福建黄兔和塞北兔断奶窝重与断奶头数、泌乳力密切相关的研究结果相似。

断奶窝重是母兔繁殖性能各项指标的总体现，是反映母兔繁殖力的一项重要经济性状。本试验结果表明，断奶只数对断奶窝重的决定系数最大（0.300 85），直接影响也最大，其通径系数为0.548 5，断奶只数、21日龄窝重共同对28日龄断奶窝重的相对决定程度为0.230 37，其绝对值在各决定系数中居第2位；直接影响位居第2位的也是21日龄窝重，其通径系数为0.373 0，说明在注意断奶只数的同时，还应注意21日龄窝重这一性状。

回归方程中剩余通径系数值（0.574 0）较大和决定系数分析中误差项对断乳窝重的相对决定程度为0.329 4，表明测定的有关性状资料影响断乳窝重的主要性状已包括在内，但还可能存在观测值误差或者可能还有对断乳窝重有影响的性状和因素未被考虑到，有待进一步研究。

本研究采用了相关分析和通径分析方法，通过计算各个繁殖性状间的直接影响和间接影响，揭示了各个繁殖性状间的真实关系，为选育工作提供了理论依据。从相关分析中了解到，断奶窝重是一个复合性状，它与其他繁殖性状关系密切，可重点选择该指标。在生产实际中断奶窝重也是母兔繁殖性能各项指标的总体体现。由通径分析可知，断奶只数对断奶窝重的直接影响最大，其次是21日龄窝重。因此，在生产实际中，提高黑白花兔断奶窝重最重要的是要提高断奶只数，断奶仔数也是决定母兔年生产能力的重要指标，所以要通过加强饲养管理，健全仔兔保健等措施来提高窝成活数；其次通过选育和加强母兔饲养管理来提高母兔泌乳量，以保证仔兔在哺乳阶段获得充足的营养，从而从整体上提高黑白花兔的繁殖性能。