DB65∕T 4460-2021 中国西门塔尔牛遗传评定技术规程(新疆维吾尔自治区).pdf

ICS 65. 020. 30 CCS B 43 0065 新疆维吾尔自治区地方标准DB 65/T 4460-2021 中国西门塔尔牛遗传评定技术规程The technical regulation of genetic assessment for Chinese simmental 2021 -11 -08发布2022 -01一01实施新疆维吾尔自治区市场监督管理局发布DB 65/T 4460-2021 目次前言.IIl 范围.12 规范性引用文件.13 术语和定义.14 遗传评估性状和筛选标准.14. 1 评估对象.14.2 评估性状和筛选标准.14.2.1 305 d产奶量.14.2.2 乳脂率.24.2.3 乳蛋白率.24.2.4 成母牛体重.24.2.5 外貌评分.25 评估方法.25. 1 遗传参数估计.25.2 育种值估计.25.2.1 305 d产奶量评估模型.25.2.2 乳脂率评估模型.35.2.3 乳蛋白评估模型.35.2.4 体重评估模型.35.2.5 外貌评分评估模型.35. 3 综合选择指数.I DB 65/T 4460-2021 目IJ吕本文件按照GB/T1. 1-2020 标准化工作导则第1部分:标准化文件的结构和起草规则的规定起草。本文件由新疆维吾尔自治区畜牧兽医局提出。本文件由新疆维吾尔自治区畜牧兽医局归口并组织实施。本文件起草单位:新疆农业大学、中国农业大学、新疆维吾尔自治区畜牧总站、新疆呼图壁种牛场有限公司、新疆天山畜牧生物工程股份有限公司。本文件主要起草人:黄锡霞、王雅春、程黎明、张晓雪、魏趁、张梦华、王丹、霄磊、葛建军、肖锐、帕尔哈提木铁力甫、王旭光、谭世新、李红燕、田晓阳、贾旭升、王爱芳、操礼军、董明明、种丽伟、刘婷婷、王彩云、苟穆荣。对本文件的修改意见建议，请反馈至新疆维吾尔自治区畜牧兽医局(乌鲁木齐市新华南路408号)、新疆农业大学(乌鲁木齐市农大东路311号)、新疆维吾尔自治区市场监督管理局(乌鲁木齐市新华南路167号)。新疆维吾尔自治区畜牧兽医局联系电话[phone];传真[phone];邮编830004。新疆农业大学联系电话[phone];传真[phone];邮编830052。中国农业大学联系电话[phone];传真[phone];邮编100193。新疆维吾尔自治区畜牧总站联系电话[phone];传真[phone];邮编830049。新疆呼图壁种牛场有限公司联系电话[phone];传真[phone];邮编8312030新疆天山畜牧生物工程股份有限公司联系电话:[phone];传真:0994一6566618;邮编:831100。新疆维吾尔自治区市场监督管理局联系电话:[phone];传真:[phone];邮编:8300040 II DB 65/T 4460-2021 中国西门塔尔牛遗传评定技术规程1 范围本文件规定了中国西门塔尔牛术语和定义、遗传评估性状和筛选标准、评估方法等的要求。本文件适用于新疆地区中国西门塔尔牛种牛遗传价值的评估。2 规范性引用文件下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款。其中，注日期的引用文件，仅该日期对应的版本适用于本文件;不注日期的引用文件，其最新版本(包括所有的修改单)适用于本文件。GB 19166-2003中国西门塔尔牛3 术语和定义下列术语和定义适用于本文件。3. 1 育种值breedingvalues 一个个体作为亲本(种畜)的价值，是衡量个体遗传价值高低的指标，利用各种可利用的表型或遗传信息和统计方法，将加性遗传效应从表型值中剖分出来。3.2 最佳线性无偏预测法bestlinear unbiased prediction 以线性混合模型为基础的育种值估计方法。可以在估计育种值的同时对系统环境误差进行估计和校正，以下简称BLUP。3.3 选择指数selectionindex 家畜选种工作中进行多性状选择时，根据育种目标中各性状的遗传基础和经济重要性，分别给予各性状一个适当的加权系数，综合到一个指数中，根据指数的大小进行选择。4 遗传评估性状和筛选标准4.1 评估对象新疆地区中国西门塔尔牛种牛。4.2 评估性状和筛选标准4. 2. 1 305 d产奶量自产棋后第一天开始到305d的总产奶量。泌乳天数不足305d或超过305d用中国西门塔尔牛305d 产奶量校正系数进行校正，校正系数见附录Ao305 d产奶量的筛选标准为2000kg13000 kg。DB 65/T 4460-2021 4.2.2 乳脂率乳中脂肪含量占全乳的百分比。由于产奶量数据是305d产奶量，应将测定日的乳脂率转化为305d 的数据，平均乳脂率按公式(1)计算。乳脂率的筛选标准为2.0%7. 0%。F=平(1)式中:F一一平均乳脂率(乳蛋白率)，单位为百分比例); IM-每一胎次测定日乳脂率(乳蛋白率)总和，单位为百分比例); N一一该胎次测定日乳脂率(乳蛋白率)测定月份数。4.2.3 乳蛋白率乳中蛋白含量占全乳的百分比。由于产奶量数据是305d产奶量，应将测定日的乳蛋白率转化为305d的数据，平均乳蛋白率按式(1)计算。乳蛋白率的筛选标准为2.0%5. 0%。4.2.4 成母牛体重成年母牛体重连续测定空腹重2 d，取其平均值。测定时，要求用分度值0.1kg的磅秤称量，单位为千克(kg)，结果保留2位小数。4.2.5 外貌评分应按照GB19166一2003中5.1条的规定进行外貌评定。5 评估方法5.1 遗传参数估计采用最小二乘法分析不同胎次、产棋年份、产棋季节、产按年龄对中国西门塔尔牛305d产奶量、乳脂率、乳蛋白率、体重、外貌评分的影响，利用多性状动物模型平均信息约束最大似然法对中国西门塔尔牛305d产奶量、乳脂率、军L蛋白率、体重、外貌评分进行遗传参数估计。评估模型见附录B。5.2 育种值估计采用动物模型BLUP进行个体育种值估计。5.2.1 305 d产奶量评估模型Yijklm =+bi+与+Yk + sl(bD +ijlkm + eijklm. 式中:Yijklm一一个体表型观察值:一一群体均值效应;bi一一胎次效应:马一一产按年份效应;2 DB 65/T 4460-2021 Yk一一产棋季节效应;sl(bD一一产按年龄嵌套于胎次效应:ijlkm一一随机加性遗传效应:eijklm一一随机残差效应。5.2.2 乳脂率评估模型Yijkm =+bi十马+Yk+ijkm + eijkm. (3) 式中:Y批m一一个体表型观察值:一一群体均值效应:bi一胎次效应:马一一产按年份效应;Yk一一产棋季节效应:aijkm一一随机加性遗传效应;eijkm一一随机残差效应。5.2.3 乳蛋白评估模型式中:Yjkm一一个体表型观察值:一一群体均值效应:马一产按年份效应:Yk一一产棋季节效应:jkm一一随机加性遗传效应:ejkm一随机残差效应。5.2.4 体重评估模型Yjkm =+与十Yk+jkm+ ejkm. . (引Yij1m =+bi+马+sl(bD +ijlm十eijlm.(5) 式中:Yij1m一一个体表型观察值;一一群体均值效应:bi一一胎次效应:马一一产按年份效应;sl(bD一-产按年龄嵌套于胎次效应:ijlm一一随机加性遗传效应:E仙n一一随机残差效应。5.2.5 外貌评分评估模型3 DB 65/T 4460-2021 Yijklm =+与十Yk+ sl(ba +ijklm十eijklm. . . . (6) 式中:Yijklm一一个体表型观察值:一一群体均值效应:马一一产按年份效应:Yk一一产棋季节效应:sl(bi)一一产棋年龄嵌套于胎次效应:ijlkm一一随机加性遗传效应:句jklm一一随机残差效应。5.3 综合选择指数利用本文件5.1所得遗传方差计算各性状遗传标准差，将本文件5.2所得305d产奶量、乳脂率、乳蛋白率、体重、外貌评分的育种值依据相对经济重要性加权构成一个总性能指数TPI。TPI = 13 x乒!L+乒立+乒+乒旦旦+乒主1x 50 . (7) MYB FPB PPB MCWB TCBJ 式中:MYB一-305d产奶量的育种值;SMYB一-305d产奶量遗传标准差:FPB一一乳脂率的育种值:SFPB一一乳脂率遗传标准差:PPB一一乳蛋白率的育种值:SPPB一一乳蛋白率遗传标准差:MCWB-体重的育种值:SMCWB一一体重遗传标准差:TCB一一外貌评分的育种值:SrCB一一外貌评分遗传标准差:50一一指数系数。4 附录A(资料性产奶量校正系数人1中国西门塔尔牛泌乳不足305d校正系数DB 65/T 4460-2021 泌乳天数不足305d用中国西门塔尔牛泌乳不足305d校正系数进行校正，校正系数见表A.1。表A.1中国西门塔尔牛泌乳不足305d校正系数泌乳第一胎第二胎第三胎第四胎四胎泌乳第一胎第二胎第三胎第四胎四胎天数以上天数以上90 3. 109 2.847 2. 742 2. 764 2. 781 200 1. 411 1. 387 1. 332 1. 350 1. 359 95 2.903 2. 665 2. 578 2.585 2. 605 205 1. 385 1. 363 1. 307 1. 327 1. 332 100 2.742 2. 514 2.427 2.437 2.456 210 1. 364 1. 342 1. 287 1. 305 1. 316 105 2.626 2. 408 2.322 2.334 2. 352 215 1. 328 1. 309 1. 259 1. 272 1. 284 110 2.491 2. 325 2. 249 2. 263 2.271 220 1. 298 1. 279 1. 230 1. 243 1. 255 115 2. 397 2.241 2. 159 2. 182 2. 180 225 1. 276 1. 259 1. 210 1. 224 1. 236 120 2.324 2. 168 2. 091 2.107 2.119 230 1. 248 1. 243 1. 199 1. 211 1. 221 125 2.210 2. 068 1. 999 2.007 2. 022 235 1. 229 1. 225 1. 181 1. 195 1. 202 130 2.118 1. 982 1. 912 1. 923 1. 936 240 1. 214 1. 210 1. 168 1. 179 1. 190 135 2.051 1. 921 1. 851 1. 863 1. 877 245 1. 187 1. 184 1. 147 1. 154 1. 166 140 1. 970 1. 873 1. 808 1. 822 1. 830 250 1. 165 1. 161 1. 126 1. 132 1. 144 145 1. 913 1. 822 1. 753 1. 774 1. 774 255 1. 149 1. 145 1. 113 1. 118 1. 129 150 1. 868 1. 778 1. 710 1. 727 1. 737 260 1. 127 1. 133 1. 105 1. 109 1. 118 155 1. 795 1. 713 1. 651 1. 663 1. 675 265 1. 113 1. 119 1. 092 1. 097 1. 102 160 1. 736 1. 657 1. 594 1. 608 1. 619 270 1. 102 1. 107 1. 084 1. 085 1. 094 165 1. 692 1. 617 1. 553 1. 569 1. 581 275 1. 081 1. 085 1. 068 1. 066 1. 074 170 1. 638 1. 586 1. 525 1. 543 1. 550 280 1. 063 1. 066 1. 051 1. 049 1. 056 175 1. 600 1. 551 1. 488 1. 510 1. 513 285 1. 051 1. 053 1. 040 1. 038 1. 044 180 1. 570 1. 522 1. 460 1. 479 1. 489 290 1. 034 1. 042 1. 034 1. 032 1. 035 185 1. 520 1. 477 1. 419 1. 434 1. 446 295 1. 024 1. 030 1. 023 1. 024 1. 023 190 1. 479 1. 437 1. 379 1. 395 1. 406 300 1. 016 1. 019 1. 015 1. 015 1. 016 195 1. 449 1. 409 1. 351 1. 368 1. 379 305 1. 000 1. 000 1. 000 1. 000 1. 000 5 DB 65/T 4460-2021 人2中国西门塔尔牛泌乳超过305d校正系数泌乳天数超过305d用中国西门塔尔牛泌乳超过305d校正系数进行校正，校正系数见表A.2。表A.2中国西门塔尔牛泊乳超过305d校正系数泌军L第1台第二胎第二胎第四胎四胎泌乳第一胎第二胎第三胎第四胎四胎| 天数以上天数以上305 1. 000 1. 000 1. 000 1. 000 1. 000 370 0.976 0.840 0.801 0.974 0.947 310 0.987 0.983 O. 984 0.987 O. 985 375 O. 998 0.872 O. 781 O. 989 0.961 315 0.979 0.970 0.972 O. 980 0.976 380 1. 022 0.827 O. 764 1. 014 O. 982 320 0.968 0.960 O. 965 0.977 0.970 385 1. 061 0.822 0.743 1. 043 1. 007 325 。.9620.949 O. 952 O. 973 0.961 390 1. 115 0.819 O. 724 1. 081 1. 047 330 0.960 0.938 O. 941 0.969 0.958 395 1. 175 0.813 O. 701 1. 121 1. 088 335 0.950 0.921 0.924 0.959 0.947 400 1. 261 0.810 0.677 1. 176 1. 146 340 0.945 0.905 0.905 0.952 0.938 405 1. 391 0.812 0.657 1. 260 1. 234 345 0.945 0.894 O. 890 0.951 O. 935 410 1. 571 0.818 0.638 1. 385 1. 363 350 O. 943 0.885 。.877O. 956 O. 935 415 1. 870 0.828 0.617 1. 557 1. 548 355 0.948 0.875 O. 860 0.959 0.935 420 2.414 0.842 0.598 1. 819 1. 857 360 O. 958 0.866 0.844 0.965 O. 941 425 3. 516 0.857 O. 578 2.224 2. 367 365 O. 964 0.852 0.823 0.967 O. 942 430 7.596 0.879 0.557 3.016 3.494 6 B. 1 非遗传因素影响分析附录B(规范性)育种值估计DB 65/T 4460-2021 利用最小二乘法分析不同胎次、产棋年份、产棋季节、产按年龄中国西门塔尔牛305d产奶量、乳脂率、乳蛋白率、体重、外貌评分的影响，将对各性状有显著影响的效应作为B.2遗传参数和育种值估计模型中的固定效应。模型为:Yklm =十bi+马+Yk十sl(bD+ eijklm(B.1) 式中:Yijklm一一个体表型观察值:一一群体均值效应:bi一一胎次效应:与一一产棋年份效应:Yk一一产棋季节效应;sl(bD一一产较年龄嵌套于胎次效应:eijklm一一随机残差效应。B.2 遗传参数估计利用多性状动物模型平均信息约束最大似然法对中国西门塔尔牛305d产奶量、乳脂率、乳蛋白率、体重、外貌评分进行遗传参数和育种值估计，模型为:式中:Y一一性状的观察值向量:b一一固定效应向量:一一动物个体加性遗传效应向量;X一一固定效应关联矩阵;Z一一加性遗传效应关联矩阵:e一一随机残差向量。模型假设前提条件:Y=Xb+Z+ e . . . . . . . . (B. 2) E(y) = X卢，E() = 0, Vr() = A(T， Var(m) = A马，Cov(，m) =Am， Var(e) = A(T 根据计算得到的方差组分按式(B.3), (B.4) , (B.5)计算遗传力、遗传相关、表型相关:h2 = (T/. (B. 3) 仄习:.(B. 4) 一苦?(B5) DB 65/T 4460-2021 式中:h2一一遗传力:一一加性遗传方差:3一-表型方差:rA一一性状间的遗传相关:rp一一性状间的表型相关:Cov(AVA2)一一性状间的遗传协方差:COV(Pl，P2)一一性状间的表型协方差。8