李佳莹,梁 龙,贺 鑫,巴音花,何鹏飞,张 宁,贺三刚*
(1.新疆农业大学动物科学学院,乌鲁木齐 830052;2.新疆畜牧科学院生物技术研究所/农业农村部草食家畜遗传育种与繁殖重点实验室/新疆动物生物技术重点实验室,乌鲁木齐 830011)
在绵羊的选种选育中,为了便于生产管理,减少人力资源的经济成本,管理者通常采取统一管理的方式进行统一配种、统一注射疫苗、统一饲喂、统一培养、统一接羔,统一断奶,统一测量体重和体尺体况。这些指标检测都要在同一天进行,但是我们不能保证母羊同时发情,也不能保证羔羊统一出生,更不能计算羔羊的断奶天数让它们刚好满100 d进行断奶,苏保全[1]指出:断乳羔羊育肥是羊肉生产的主要形式。断奶重作为生产发育和评定羔羊重要的指标,同一天断奶的数据并不能作为评价羔羊断奶重的指标。根据出生类型和性别进行划分羔羊类型进行断奶重的校正,使断奶体重校正到同一水平,可以为接下来的研究提供理论基础,使羔羊的选种选育更加科学。
绵羊不仅是人类需要的主要肉类食品和辅助畜产品的来源,也是人类农业文明传播的主要组成部分[2]。断奶重是反映羔羊早期生长发育健康状况、羔羊生产性能的重要性状,受到遗传、环境和母体效应等多种因素的影响[3]。羔羊的断奶鉴定是绵羊育种工作中早期选种的重要环节,尤其是在种公羊的选种时,科研人员和农场养殖户会根据断奶时期羔羊的体重作为评判羔羊好坏的唯一指标。断奶体重作为一项重要的中间经济性状,是影响家畜生产能力的决定指标,倍受畜牧工作者的重视[4]。此时便突出了羔羊断奶体重的重要性,成年母羊配种周期长,通常为45 d左右,在日后的产羔过程中,难免会造成羔羊出生日期的差异,而牧场和农户为了节约时间、经济、人力的成本通常会选择统一断奶,彼时测定的断奶体重严重缺乏说服力,不能作为评判羔羊优劣的标准。将不同日龄断奶羔羊体重校正到标准化的同一水平上,可将优秀的羔羊集中在优秀的羊群中饲喂培养,为社会挑选出优良的羊种。
1.1.1 数据来源
采用新疆维吾尔自治区伊犁哈萨克自治州昭苏马场萨福克种羊场2021年7月中旬断奶鉴定数据共1 030条,羔羊品种为特哈杂交F2代羔羊,即以特克塞尔羊为父本,以特克塞尔×哈萨克F1代羊为母本,杂交繁育获得级进F2代羊。
1.1.2 试验数据
采集的数据包括羔羊个体号、性别、出生日期(5月出生羔羊数量较少,出生日期太晚,不在人工配种预期的产羔周期内,故将其剔除)、母羊年龄、单双羔(由于3胎数量极少<10只,故将3胎数据划归到双胎数据之中)、初生重、断奶重、死亡记录等。将数据录入Excel表格中进行整理,为保证数据的准确性,剔除死亡、缺失数据。羔羊断奶重数据记为平均数加减3倍标准差来确定,筛选删除异常值,其中包括羔羊死胎数据101条;出生于2月份羔羊数据3条;出生于5月份羔羊数据36条;不明羔羊性别的数据1条;不明羔羊母亲年龄数据17条;共删除异常值158条,最终有效数据872条。
1.2.1 水平划分
将试验类别按照出生类型和性别进行划分为四个类型,分别为:单胎公羊、单胎母羊、双胎公羊、双胎母羊。根据其各个类型的日龄和断奶重进行校正。
1.2.2 统计分析
运用SPSS26.0统计软件进行描点绘制散点图,运用双变量相关程序求得断奶日龄与断奶重的相关系数,并进行线性回归分析,最终得到理想校正值。
2.1.1 断奶体重与断奶日龄的散点图
为了简洁直观的观察断奶日龄与断奶重之间的变化趋势,特制定平面散点图将每一对观测值描点,形成散点图,见图1,横轴x轴代表断奶日龄(d),断奶重(kg)在y轴上表示。初步判断断奶日龄与断奶重
图1 断奶日龄与断奶重的散点图
[4]呈直线回归关系。
2.1.2 相关系数与回归方程
为了将影响断奶校正的因素降至最低,故将出生类型与性别分类校正。在断奶日龄中显示,四组类别的断奶日龄差距跨度较大,范围较广,增加了不同日龄断奶羔羊同一时间断奶造成羔羊断奶重的误差,体现了校正断奶重的必要性。
表1、2、3、4的结果显示,单胎公羊、单胎母羊、双胎公羊、双胎母羊的断奶日龄与断奶重的回归方程分别为y=0.205x+12.404;y=0.195x+11.472;y=0.212x+3.218;y=0.216x+1.203。四组类别的相关系数分别为0.468、0.571、0.503、0.516,均有极显著的影响(P<0.01),并大于总体相关系数。按照系数显著,回归方程一定显著;系数不显著,回归方程一定不显著的原则,表中的结果数据可信度高,为接下来的分析提供了理论基础。
表1 羔羊断奶日龄和断奶体重相关系数和回归方程
2.1.3 断奶体重理论值和校正系数
断奶日龄每间隔5 d划分为一组,取中位数和目标断奶日龄100 d,将其代入单胎公羊表2、单胎母羊表3、双胎公羊表4、双胎母羊表5的回归方程,分别得出断奶重的理论值和校正日龄体重值。校正系数为校正日龄体重值除以断奶重的理论值。最后,用校正系数乘以实际断奶重得出校正体重值。
表2 单胎公羊不同日龄断奶重理论值和校正系数
表3 单胎母羊不同日龄断奶重理论值和校正系数
表4 双胎公羊不同日龄断奶重理论值和校正系数
表5 双胎母羊不同日龄断奶重理论值和校正系数
本文所采用“日龄校正方法”结果显示回归方程达到极显著的水平(P<0.01),说明用本文的四组回归方程校正系数是可靠的。而散点图显示断奶日龄与断奶体重呈正相关,说明了校正断奶体重的必要性,能够进行更好的选种选育。王瑞金[5]郑浩[6]等人同样运用了求得相关系数与回归方程最后通过日龄来校正体重。
在早些年前,当时的计算机系统还不算发达的时候,李蒙召[7]张延龄[8]通过影响羔羊断奶重的因素如单双羔、胎次、哺乳期等进行羔羊分类,提出校正羔羊断奶重的方法,运用日增重法、回归法、多项式的方法校正羔羊的体重。齐侃虎等[9]运用最小二乘法制定校正初生重与断奶重的系数。为后续新兴的计算机统计软件算法奠定了基础。刘望宏[10]将三种不同品种的猪进行了早期和全期生长性状的体重测定,运用线性回归及非线性回归的模型,根据数据的特点进行了抽样、分阶段、连续测定的方法,选取了最适合的模型。结果显示:Logistic、Gompertz、Saturation等模型适合体重性状与日龄的关系,相关拟合平均指数达到0.9,达到极显著(P<0.01)其中Gompertz模型更具有生物学意义,Saturation的模型回归更接近一条直线,最适合做校正公式。吕凤林[11]对内江仔猪运用不同的校正方法并进行比较分析,结果以平均日增重和Tayl or公式对内江仔猪体重校正较理想。邵学芬等[12]采用“全样本”、“等日龄平均”、“期间5日平均”、“升序10羊平均”四种模型断对奶日龄的分割方法进行断奶重的校正,结果显示,四种方法的拟合度均很高,t检验的结果均为极显著(P<0.01),通过相关系数的比较和与原始样本断奶重的差异分析,“升序10羊平均”的模型最适合。此种方法是近年来国内鲜见的仅用断奶日龄划分模型的方法,具有创新性,对断奶重校正提供了诸多灵感。
对于断奶重校正的方法,大多数的科研技术人员会采用最简便的“日增重校正法”,通过羔羊断奶前的平均日增重乘以校正天数。但此种方法的弊端在于不能保证羔羊每天的日增重相同,误差较大。而建立数学模型进行断奶重校正,保证其相关系数显著、拟合值接近1,t检验显著,产生的误差远远小于日增重法。而建立数学模型需要在养殖场和牧民中的相关科技人员掌握一定的理论知识和学科基础。但目前看来大部分人员实际动手操作能力强,但理论知识都达不到相关要求,很难建立起数学模型。
通过划分生长类型和性别不同分别建立四种断奶天数与断奶重双变量相关系数,得出断奶天数和断奶重均有极显著相关(P<0.01)的结果,建立散点回归图,通过趋势线初步判断断奶天数与断奶重呈线性回归,再建立断奶重与断奶天数的线性回归方程,按照5日间隔将断奶天数划分14组,计算各组中数,分别代入目标校正天数100 d求出校正系数,与实际断奶重相乘得出校正断奶重。为后续的选种选育建立特哈杂交F2代种羊遗传评估体系提供了一定的参考标准。