特哈杂交羔羊消化道线虫的遗...参数估计及其与生长性状关系_贺鑫.pdf

2023年1月第1期（总第218期）草食家畜（双月刊）特哈杂交羔羊消化道线虫的遗传参数估计及其与生长性状关系贺鑫1，李佳莹1，何鹏飞1，巴音花1，梁龙2，贺三刚2*（1.新疆农业大学动物科学学院，乌鲁木齐830052；2.新疆畜牧科学院生物技术研究所/农业农村部草食家畜遗传育种与繁殖重点实验室/新疆动物生物技术重点实验室，乌鲁木齐830011）摘要：旨在分析特哈杂交羔羊的组别、性别、群别、母羊年龄、出生类型等固定效应对消化道线虫虫卵数和早期生长性状的影响，目标性状包括羔羊每克粪便中的虫卵数（FEC）、羔羊初生重（BW）、羔羊断奶重（WW）、羔羊6月龄体重（6WT）、羔羊出生到6月龄平均日增重（ADG）。对目标性状作出遗传参数估计（遗传力、方差组分、遗传相关和表型相关），为绵羊的遗传改良提供数据参考。使用ASReml软件进行分析，先利用混合线性模型分析固定效应对绵羊的消化道线虫虫卵数和早期生长性状的影响，再使用多性状动物模型对目标性状进行遗传参数估计。研究表明，羔羊的组别和群别对所有性状影响极显著（P0.001）；性别对虫卵数、初生重、断奶重和平均日增重影响极显著（P0.001）；母羊年龄对初生重影响极显著（P0.05）；出生类型对初生重、断奶重、6月龄重、平均日增重影响极显著（P0.001）；FEC、BW、WW、6WT、ADG的遗传力分别为0.142、0.156、0.365、0.294、0.024；FEC与BW（-0.889）之间均呈现出比较高的负遗传相关，FEC与6WT（0.096）、ADG（0.152）之间呈较低正遗传相关，与WW（-0.574）之间呈负遗传相关。因此，在遗传改良和选育过程中，不仅要提高对线虫感染的控制，还要重视并改善绵羊的管理方式、生长环境。本研究估计了特哈杂交羔羊消化道线虫虫卵数及早期生长性状的遗传参数，可为下一步的遗传改良提供理论参考。关键词：绵羊；消化道线虫；影响因素；遗传力中图分类号：S858.26文献标识码：A文章编号：10036377（2023）01001-08消化道线虫感染是影响中国和全世界绵羊或其它放牧反刍动物最重要的健康问题之一。据统计，澳大利亚绵羊业每年因线虫感染而造成的生产损失和治疗成本为4.36亿美元，预计在全球范围将达到数百亿美元1。寄生虫对放牧反刍动物的健康和生产力的影响已得到充分证明，包括体重减轻、腹泻、厌基金项目：新疆维吾尔自治区重点实验室开放课题“全基因组关联分析鉴定绵羊消化道线虫抗性基因”；国家重点研发政府间国际合作专项“气候变化下地方绵羊经济性状改良的组学研究”（2017YFE0117900）；新疆维吾尔自治区现代畜禽种业提升专项“牧区地方肉羊品种选育推广技术体系”。作者简介：贺鑫（1998-），男，陕西榆林人，硕士研究生，研究方向为动物遗传育种。E-mail：通讯作者：贺三刚（1982-），男，甘肃庆阳人，研究员，研究方向为动物遗传育种。E-mail：收稿日期：2022-10-07，修回日期：2022-10-2710.16863/ki.1003-6377.2023.01.00112023年1月第1期（总第218期）草食家畜（双月刊）食、贫血和死亡，在Singh等2的一项研究中，在每克粪便中虫卵数2 000的情况下，平均每只羊损失活重5.06 kg。在过去的几十年里，养羊业越来越依赖驱虫治疗作为寄生虫控制的方法。然而，驱虫治疗价格昂贵，治疗时的频繁用药也导致全世界绵羊的驱虫药抗性迅速增加3，使治疗效果大大降低4。因此，养殖户和育种人员需要寻找替代的方法来抵抗寄生虫的感染，以保证他们羊群的健康和生产力。有证据表明，一些绵羊种群对消化道线虫感染具有遗传抗性，随着时间的推移，这些选定的绵羊种群将对寄生虫感染产生抗性，用药将最小化甚至不再需要给药5。而对线虫抗性屠宰前的评估只能通过测量粪便中线虫虫卵的数量来进行，虫卵数（FEC）被证明是提高群体内对线虫抗性的有效选择标准，是一种简单且应用最广泛的性状6，7。因此，绵羊的选择性育种能增强对寄生虫感染的抵抗力8，9，提供了一个管理绵羊线虫感染的替代方法。准确的遗传参数估值可以为选种提供可靠的依据，以引进的特克塞尔羊为父本、本土品种哈萨克羊为母本生产的特哈杂交羔羊为试验动物，既有父系日增重和胴体重较高等特点，又有母系抗逆性和环境适应性较强等特点10。因此对其虫卵数的研究特别重要，选育好了既可提高对线虫的感染抗性，也可提高生长性状。本研究旨在评估特哈杂交羔羊虫卵数和早期生长性状遗传参数，为下一步合理制定羔羊的育种方案提供理论支撑。1材料与方法1.1试验动物本研究所有试验动物均来源于新疆维吾尔自治区伊犁哈萨克自治州昭苏县昭苏马场哈萨克羊繁育基地。为改良哈萨克羊的生产性能，课题组将哈萨克羊作母本、特克塞尔作父本进行杂交获得杂种F1代，再以F1代为母本与特克塞尔羊级进杂交获得子代即F2代。2021年3、4月份所有新生羔羊按4种杂交组合跟随母羊分为7个群体、封闭管理。在五月中旬，所有羔羊跟随母羊迁往牧区，每天放牧810 h。羔羊在断奶前驱虫，6月龄时进行粪便采样，直肠取粪法采集粪便，再通过改良的McMaster技术进行虫卵计数获得虫卵数数据，通过检测发现羔羊线虫感染率为100%。并收集有关2021年新生羔羊的配种产羔记录、各阶段体重记录和断奶鉴定记录。测定性状包括羔羊每克粪便中的虫卵数（Faecal Nematode Egg Count,FEC）、羔羊初生重（Birth Weight，BW）、羔羊断奶重（Weaning Weight,WW）、羔羊6月龄体重（6-Month Weight,6WT）、羔羊出生到6月龄平均日增重（Average Daily Gain At 6 Months Of Age,ADG）。1.2研究方法1.2.1数据处理共选取2021年生产的892只羔羊，使用Excel剔除表型记录不准确，筛选系谱错误及缺失的记录后，共获得554只羔羊，每只羔羊在感染期间的每克粪便中的虫卵数（FEC）取平均值。由于虫卵数不是正态分布，而是偏态分布，数据自然对数转换为LFEC=loge（FEC+100）。所有统计检验均应用转换后的数据进行分析11。对筛选后的数据进行描述性统计。1.2.2可行性分析使用ASReml统计软件绘制性状对应的残差直方图、残差和拟合值的散点图、Q-Q图以及残差和单元值的散点图，通过检测所用模型的准确性和数据做估计遗传参数的可行性。1.2.3固定效应处理影响绵羊线虫虫卵数和早期生长性状的固定效应包括组别效应、性别效应、群别效应、母羊年龄和出22023年1月第1期（总第218期）草食家畜（双月刊）生类型。组别效应：所有羔羊按不同杂交组合分为（F1F2、F2F1、F2哈萨克羊、纯特克塞尔F1）4个水平；性别效应：按照（公、母）划分为2个水平；群别效应：羔羊按不同群别划分为1、2、3、4、5、6、7个水平；母羊年龄：按照产羔母羊不同年龄划分为（2岁、3岁、4岁、5岁）4个水平；出生类型：按照（单羔、双羔）划分为2个水平。1.2.4固定效应的显著性检验通过ASReml的混合线性模型对五个性状进行分析，检验纳入模型的固定效应的显著性。采用的统计分析模型表达式为：Yjklmn=+zi+sj+ak+tm+wl+eijklmn其中，Yjklmn代表表型性状观察值；为总体均值；zi为羔羊组合群别第i个水平效应；sj为羔羊性别第j个水平效应；ak为羔羊群别第k个水平效应；tm为母羊年龄第m个水平效应；wl为出生类型第l个水平效应；eijklmn为随机残差效应。1.2.5遗传参数估计使用ASReml软件中的多性状动物模型对绵羊线虫虫卵数及早期生长性状进行遗传参数估计，其模型如下所示：yXbZae式中，y为个体的表型观测值；b是固定效应向量，包括羔羊的组别效应、性别效应、群别效应、母羊年龄和出生类型；X为b的关系矩阵；a为个体的加性遗传效应向量；Z为a的结构矩阵；e为随机残差效应向量。2结果与分析2.1各性状的描述性统计表1分别列出了羔羊线虫虫卵数及早期生长性状的观察值样本数、均值、标准差、最大值、最小值和变异系数等统计量。羔羊各性状变异系数范围在9.88%16.18%，羔羊消化道线虫虫卵数的变异系数与羔羊早期生长性状的变异系数相比较小。表1各性状的描述性统计性状样本数均值标准差最大值最小值变异系数FEC5547.0140.6939.425.019.88BW/kg5545.2170.8447.32.016.18WW/kg54428.7174.45541.5417.3815.516WT/kg45932.7254.88948.3220.5114.94ADG/g459184.04724.752258.649121.26313.452.2残差分析利用ASReml软件检验模型和数据做遗传参数估计的可行性。以羔羊每克粪便中的虫卵数（FEC）为例（图1），残差直方图呈正态分布趋势（图1a），Q-Q图近似45直线（图1c）,残差均匀的分布在直线的两方，与单元数之间没有明显的趋势（图1b和图1d）。结果说明，残差基本符合假设要求，且方差相同，为独立的正态分布，各个性状的表型数据均符合后续分析的条件。32023年1月第1期（总第218期）草食家畜（双月刊）残差样本数量残差-1012210-1-2210-1（a）（d）（c）-2026.46.87.27.602004006040200210-1注：（a）残差直方图；（b）残差和拟合值的散点图；（c）QQ图；（d）残差和单元值的散点图图1绵羊消化道线虫虫卵数残差图2.3固定效应显著性检验由表2的分析结果可知，羔羊的组别和群别对所有性状影响极显著（P0.001）；性别对虫卵数、初生重、断奶重和平均日增重影响均极显著（P0.001），对6月龄体重有影响（P0.05）；母羊年龄对初生重影响极显著（P0.001），对平均日增重影响显著（P0.01），对断奶重、6月龄重有影响（P0.05）；出生类型对初生重、断奶重、6月龄重和平均日增重影响极显著（P0.001）。表2固定效应显著性检验FECBWWW6WTADG组别DF33333P1.077e-06*2.2e-16*2.2e-16*2.2e-16*1.311e-12*性别DF11111P4.020e-07*1.256e-06*4.343e-14*0.02004*2.475e-06*群别DF66666P1.899e-12*2.121e-09*3.343e-10*1.998e-10*1.408e-8*母羊年龄DF44444P0.2259ns2.209e-06*0.0432*0.02161*0.001034*出生类型DF11111P0.5205ns 2.2e-16*2.2e-16*1.926e-13*4.020e-07*（b）注：ns表示P0.05；*表示P0.05；*表示P0.01；*表示P0.001。42023年1月第1期（总第218期）草食家畜（双月刊）2.4各性状的方差组分及遗传力使用ASReml软件中的动物模型对羔羊虫卵数、初生重、断奶重、6月龄重、平均日增重进行方差组分及遗传参数估计。表3为特哈杂交羔羊各性状方差组分和遗传参数结果。特哈杂交羔羊各性状的遗传力范围在0.020.37之 间，其 中FEC（0.1420.113）、BW（0.1560.110）、6WT（0.2940.186）、WW（0.3650.195）、ADG（0.0220.082）。表3方差组分及遗传力性状加性遗传方差残差方差表型方差遗传力标准误FEC0.0570.3460.4030.1420.113BW0.0680.3690.4370.1560.110WW5.1718.97614.1470.3650.1956WT4.50610.82415.330.2940.186ADG10.659440.442451.1010.0240.0772.5各性状的遗传相关及表型相关通过ASReml软件中多性状模型估计特哈杂交羔羊虫卵数、初生重、断奶重、6月龄重和平均日增重的遗传相关和表型相关，结果见表4。特哈杂交