生态因子对中国地方山羊品种繁殖性状的影响

胡紫意，王妍，李祥龙⋆，巩元芳

（1.河北科技师范学院动物科技学院，河北秦皇岛 066004；2.河北省秦皇岛水文勘测研究中心，河北秦皇岛 066000；3.河北科技师范学院动物科技学院，河北省特色动物种子资源挖掘与创新重点实验室，河北秦皇岛 066004）

近年来， 随着我国经济科技不断发展和全国人口的增加， 人们对优质肉食品需求越来越大，我国山羊业发展迅速。 山羊是家畜中地域分布最广，且其他家畜无法生存的环境，山羊都能生存和繁衍[1]。 山羊具有繁殖率高、适应性强、好管理等特点。 我国土地辽阔生态环境差异极大，不同的气候、海拔、气温、降水量的地域都有山羊存在。 中国山羊饲养历史久远，据记载，夏商时代就有饲养山羊记录。

人类定居后，活动范围相对缩小，随着人们的驯养、繁殖，自然环境的选择，山羊的地方品种开始出现[2]。 在不同的生态环境下，也形成了不同的体型、外貌、生产性能的品种。 按生产性能可分为毛用山羊、乳用山羊等多种用途的山羊品种。 因此，为了了解生态系统多样性对中国地方山羊品种繁殖性状的影响， 本文对我国58 个中国地方山羊品种繁殖性状与生态因子进行典型相关性分析[3]，得出繁殖性状与生态因子的内在联系，为研究中国地方山羊品种多样性与分类特征提供参考依据。

1 材料与方法

58 个中国地方山羊品种作为本次试验的对象，每组数据均来自山羊的各自原产地。 中国地方山羊品种包括2 种性状11 个变量，生态因子：平均海波（x1）、年均气温（x2）、降水量(x3)、年均日照时数(x4)、无霜期(x5），繁殖性状：羔羊断奶重（y1）、初生羔羊重（y2）、年平均产羔率（y3）、妊娠期（y4）、发情周期（y5）、初配年龄（y6）。数据来自国家畜禽遗传资源委员会组编的《中国畜禽遗传资源志：山羊志》[1]。

采用软件SPSS 进行分析，典型相关分析是利用综合变量对之间的相关关系来反映两组指标之间的整体相关性的多元统计分析方法[3-7]。

2 结果与分析

2.1 性状间的简单相关分析

58 个中国地方山羊品种的生态因子和繁殖性状11 个变量的测定值见表1。 由表1 和表2可知，11 个变量产生66 个相关关系数， 其中16个极显著水平相关 （P＜0.01），10 个显著水平相关（P＜0.05）。其中绝对值较大系数年均气温和年均降水量、年均日照时数与无霜期简单相关系数分别为0.573、0.589。

表1 地方山羊品种各性状平均值和标准差

表2 生态因子和繁殖性状间相关性分析结果

2.2 性状间典型相关分析

生态因子与繁殖性状间的典型相关分析见表3，由表3 可知，生态因子对中国地方山羊品种繁殖性状典型相关系数，第一对典型相关系数为0.916，第二对典型相关系数为0.843，第三对典型相关系数为0.402，第四对典型相关系数为0.291，第五对典型相关系数为0.196。 第一对贡献率是53.585%， 第二对贡献率是17.308%，第三对贡献率为11.489%， 第四对贡献率为9.866%，第五对贡献率为7.751%。其中第一对相关信息量占总相关信息量的53.585%，前两对相关信息量占总相关信息量的70.893%，前三对相关信息量占总相关信息量的82.383%，前四对相关信息量占总相关信息量的92.249%。说明前四对典型变量表达了两组性状间有较高的相关性，因此我们只需要讨论分析前四对变量。

表3 典型相关系数

2.3 性状间典型相关变量的构成与分析

生态因子和繁殖性状各对典型变量的构成为：

V1=-0.131x1+0.249x2-1.045x3+0.393x4+0.567x5

W1=-0.035y1+0.191y2+0.326y3+0.267y4-0.652y5+0.696y6

V2=-0.949x1+0.941x2-1.114x3-0.328x4+0.296x5

W2=0.499y1-0.034y2+0.963y3-0.014y4-0.201y5-0.12y5

V3=-0.04x1+0.325x2-0.467x3-0.204x4+1.097x5

W3=-0.657y1-0.331y2+0.532y3+0.145y4+0.2y5+0.607y6

V4=0.719x1+0.692x2-0.467x3-0.204x4-1.097x5

W4 =0.698y1 -0.48y2 +0.43y3 +0.697y4 +0.329y5+0.219y6

在生态因子和中国地方山羊品种繁殖性状间的第一对典型变量中， 年均降水量和初配年龄的相关系数最大且呈负相关， 年均降水量越大初配年龄越小。 第二对典型变量中，年均降水量与羔断奶重的相关系数最大且呈负相关，羔断奶重随着年均降水量的增加而减小。 在第三对典型变量中， 无霜期和羔断奶重的相关系数最大且呈负相关， 羔断奶重随这无霜期的增加而减小。 在第四对典型变量中，平均海拔和羔断奶重的相关系数最大且呈正相关，有正向趋势，羔断奶重随这平均海拔的增加而增加。 因此，生态因子对地方山羊品种繁殖性状间V1 和W1相关主要是由于年降水量和初配年龄引起的，两者是负向趋势，即初配年龄随年降水量的增加而减小。V2 和W2 的相关是由年降水量和羔断奶重之间密切相关引起的，两者呈负相关，即羔断奶重随着年降水量增加而减小。 V3 和W3相关是由无霜期和羔断奶重引起的，呈负相关。

3 结果与讨论

在简单的相关分析中，共有16 个极显著水平相关，10 个显著水平相关。其中绝大部分是性状内因子之间的显著， 并且典型相关分析中相关性达到显著水平的变量， 而在简单的相关分析中相关性不显著。 简单相关分析是分析变量和变量之间的关系， 相关关系是变量数值间的一种不确定的相互依存关系，就例如，山羊的生长快慢与饲喂量有关系， 但这并不能说饲喂越多山羊生长越快， 虽然它们之间有联系但存在多种影响因素，因此，得出的结果有可能也是假的，或者说结果反映的是两组性状的相关程度。而典型相关分析就是利用综合变量对之间的相关关系来反映两组指标之间的整体相关性的多元统计分析方法[8-10]。 利用这两个综合变量之间的相关关系来反映两组指标之间的整体相关性。 其分析结果是一个综合性结果，反应结果相对准确清晰。 本文对生态因子对中国地方山羊品种间的变量进行相关性分析结果表明， 年均降水量和初配年龄的相关系数最大且呈负相关，年均降水量越大初配年龄越小；年均降水量与羔断奶重的相关系数最大且呈负相关， 羔断奶重随着年均降水量的增加而减小； 无霜期和羔断奶重的相关系数最大且呈负相关， 羔断奶重随着无霜期的增加而减小。

典型相关分析研究生态因子对地方山羊品种繁殖性状的影响，揭示了其之间联系，在这么多信息中反映出之间的关联， 因此本研究所得出的结论为中国地方山羊品种遗传育种与繁殖提供科学的参考意见[11-12]。

4 结论

本文得出的最终结论是11 个变量产生66个相关关系数，其中16 个极显著水平相关（P＜0.01），10 个显著水平相关（P＜0.05）。其中绝对值较大系数年均气温和年均降水量、 年均日照时数与无霜期简单相关系数分别为0.573、0.589。典型相关分析结果表明， 两组性状间的第1 个典型相关系数为0.916，贡献率为53.585%。本文对生态因子对中国地方山羊品种间的变量进行相关性分析结果表明， 年均降水量和初配年龄的相关系数最大且呈负相关， 年均降水量越大初配年龄越小； 年均降水量与羔断奶重的相关系数最大且呈负相关， 羔断奶重随着年均降水量的增加而减小； 无霜期和羔断奶重的相关系数最大且呈负相关， 羔断奶重随着无霜期的增加而减小。