任小麗，張震＊，閆躍飛，閆磊，白雪利，王笑笑，李靜茹，汪聰勇，耿繁軍

（1.國家奶牛產(chǎn)業(yè)體系鄭州綜合試驗站，河南鄭州450008；2.河南省奶牛生產(chǎn)性能測定中心）

河南省奶牛生產(chǎn)性能測定與育種進(jìn)展

任小麗1,2，張震1,2＊，閆躍飛1,2，閆磊1,2，白雪利1,2，王笑笑1,2，李靜茹1,2，汪聰勇1，耿繁軍1

（1.國家奶牛產(chǎn)業(yè)體系鄭州綜合試驗站，河南鄭州450008；2.河南省奶牛生產(chǎn)性能測定中心）

1 河南省奶牛生產(chǎn)性能測定情況

1.1 總體情況：2008-2016年參測個體、牧場和記錄數(shù)對比

2008年至2016年8月累積參測298個牧場1 748 254頭個體共2 459 532條記錄。2016年1月至8月共200個牧場參測累積數(shù)據(jù)1 073 93頭中國荷斯坦牛583 067條測定日記錄，其中有效記錄數(shù)463 123條（個體測定次數(shù)在6次以上的記錄）、1月至8月上報數(shù)據(jù)平均有效率91.8%；平均日產(chǎn)奶量26.19 kg、乳脂率4.04%、蛋白率3.25%、脂蛋比1.25、體細(xì)胞數(shù)352.25×1 000/ml、泌乳持續(xù)力108.13、高峰產(chǎn)奶量33.40 kg、305 d產(chǎn)奶量8 260.88 kg。（數(shù)據(jù)來源于中國奶牛數(shù)據(jù)中心和河南省奶牛數(shù)據(jù)綜合服務(wù)平臺）。

圖1 2008-2016年參測牧場數(shù)、個體數(shù)及記錄數(shù)趨勢

圖2 2008-2016年測定日產(chǎn)奶量、乳成分與全國DHI平均對比

1.2 2012-2016年8月奶牛群體泌乳曲線、乳成分對比

表1 2012-2016年8月測定日產(chǎn)奶量與乳成分對比

由表1可知，相比于2015年，2016年日奶量、乳脂率、蛋白率分別增加了9.38%、6.86%和3.57%。

表2 2012-2016年8月群體泌乳曲線變化

由表2可知，相比于2015年，2016年高峰奶和305 d奶量分別增加了13.25%和20.33%。

表3 2012-2016年同期測定記錄數(shù)、測定日產(chǎn)奶量與乳成分對比

由表3可知，2016年1月至8月共上報583 067條記錄，比上年同期增加了8.28%；測定日產(chǎn)奶量、乳脂率和乳蛋白率與上年同期相比分別增加了12.92%、11.60%與4.78%；平均體細(xì)胞數(shù)相比去年降低了1.71%，平均32.52萬/毫升。

圖3 2012-2016年8月每月上報數(shù)據(jù)量

1.3 2016年1月至8月上報數(shù)據(jù)有效率

圖4 2016年1月至8月上報有效數(shù)據(jù)比例

由圖4可知，2016年1月至8月上報數(shù)據(jù)有效率平均在91.87%。

1.4 參測場單產(chǎn)水平

2016年奶牛305 d產(chǎn)奶量平均為8 260.88 kg、高峰奶量33.40 kg、平均日產(chǎn)奶量26.19 kg，均高于2012年至2015年平均水平；低產(chǎn)牧場（單產(chǎn)小于7噸）比例下降（29%）、中高產(chǎn)牧場（中產(chǎn)：7～9噸，高產(chǎn)9～12噸）比例上升（14%和9%）、單產(chǎn)12噸以上牧場4個。

表4 2012-2016年牧場305 d產(chǎn)奶量及對不同305 d產(chǎn)奶量分段牧場分布

圖5 2012年至2016年各牧場平均305 d產(chǎn)奶量分布

2 河南省奶牛生產(chǎn)性能測定中心育種工作進(jìn)展

2.1 奶牛生產(chǎn)性能測定數(shù)據(jù)質(zhì)量分析

DHI數(shù)據(jù)質(zhì)量由高紅彬于2005年首次提出，張勝利和王雅春分別于2013年和2014年提出具體考核指標(biāo)，包括：系譜完整率、胎次、泌乳天數(shù)、測定間隔、首測泌乳天數(shù)、測定次數(shù)、首次產(chǎn)犢月齡和公牛后代女兒數(shù)及其分布場數(shù)。指標(biāo)中奶牛胎次對泌乳曲線擬合效果有極顯著影響、對女兒成年后產(chǎn)奶量有顯著影響且存在負(fù)相關(guān)，即：母牛胎次越高，女兒成年后產(chǎn)奶量越低；泌乳天數(shù)對于估計305 d產(chǎn)奶量至關(guān)重要；測定間隔、首測泌乳天數(shù)和測定次數(shù)等指標(biāo)，則可以反映個體泌乳曲線的擬合效果；首次產(chǎn)犢月齡是奶牛生產(chǎn)中非常重要的指標(biāo)，不僅影響母牛本身的繁殖性能，同時還影響其泌乳性能和生產(chǎn)壽命；系譜完整率指一代和三代系譜完整率，影響遺傳參數(shù)估計準(zhǔn)確性、后代選種選配和近交系數(shù)計算；女兒數(shù)及其分布場數(shù)是公牛育種值準(zhǔn)確性的主要衡量指標(biāo)。因此，以上指標(biāo)的數(shù)據(jù)質(zhì)量直接關(guān)系到后裔測定、遺傳評估、品種改良等工作能否順利開展，也是衡量上述工作質(zhì)量的綜合指標(biāo)。

根據(jù)以上指標(biāo)，建立河南省DHI育種數(shù)據(jù)篩選標(biāo)準(zhǔn)，具體如下：①刪除沒有父親的個體（保證系譜完整）；②刪除無胎次和胎次記錄不符的記錄（1～5胎次的測定記錄予以保留）；③刪除泌乳天數(shù)在5～305 d之外的記錄；④刪除性狀表型值缺失的記錄；⑤刪除產(chǎn)犢月齡不符合奶牛生理規(guī)律的記錄；⑥刪除在泌乳早期（DIM90 d前）無測定記錄的個體；⑦刪除個體在完整泌乳期內(nèi)兩次測定間隔大于50 d的個體；⑧刪除個體在完整泌乳期內(nèi)測定次數(shù)太少（6次以下）的個體；⑨刪除公牛后代女兒數(shù)和分布不夠廣泛的（后代數(shù)在20頭以下，分布場少于5場）個體。

表5 育種值估計數(shù)據(jù)篩選條件

2008年至2016年11月，河南省奶牛生產(chǎn)性能測定數(shù)據(jù)共包括：301個牧場183 640頭個體的2 718 347條測定記錄，經(jīng)過SAS 9.4軟件由上述指標(biāo)篩選后共223場56 522頭687 249條測定記錄（占原始記錄25.28%）可用于遺傳參數(shù)和育種值估計。由圖6可以看到，對河南省DHI遺傳評估數(shù)據(jù)篩選過程中，泌乳天數(shù)無效（在5～305 d之外）的記錄數(shù)（占總記錄28.73%），月齡不符合首測泌乳天數(shù)高于90 d、無父號不符的記錄數(shù)分別占原始記錄的10.89%和9.16%。

圖6 DHI數(shù)據(jù)質(zhì)量各指標(biāo)無效數(shù)據(jù)比例

經(jīng)過篩選，可用于育種值估計的數(shù)據(jù)中，個體數(shù)占總數(shù)據(jù)的比例（個體數(shù)占比）由2008年的1.38%逐年增加至2015年的39.24%；記錄數(shù)占總記錄數(shù)的比例（記錄數(shù)占比）由2008年的1.26%逐年增加至2015年的37.77%；2010年至2015年，牧場數(shù)占總數(shù)據(jù)的比例平均為93.78%，而2015年的牧場占比（87.30%）小于2014年（96.39%）。

各牧場育種數(shù)據(jù)篩選后，沒有育種數(shù)據(jù)的牧場有78個、育種數(shù)據(jù)占總數(shù)據(jù)比例在10%以下的牧場有50個、比例在10%～30%的有96個、比例在30%～50%的牧場有71個，比例在50%以上的牧場只有6個。

圖7 育種值數(shù)據(jù)占原始數(shù)據(jù)比例隨測定年的變化趨勢

表6 各牧場育種數(shù)據(jù)占總記錄比例分布

2.2 河南省奶牛群體遺傳改良體系建設(shè)

自開展奶牛生產(chǎn)性能測定工作以來，河南省奶牛生產(chǎn)水平快速提高（單產(chǎn)提升、奶質(zhì)提高、體細(xì)胞數(shù)降低）。與奶業(yè)發(fā)達(dá)國家、省市相比，河南省奶業(yè)存在的主要問題包括：奶牛養(yǎng)殖成本高和奶牛品種登記、體型鑒定和遺傳評估工作未全面開展，影響了奶牛群體遺傳水平改進(jìn)效率。

借鑒國外奶業(yè)發(fā)達(dá)國家發(fā)展經(jīng)驗，并結(jié)合河南省奶牛養(yǎng)殖的實際情況，我們已經(jīng)建立了河南省奶牛遺傳評估平臺。如圖8所示，河南省奶牛遺傳改良體系將從品種登記平臺搭建，到牧場核心群選擇，適配凍精選擇等逐步實現(xiàn)奶牛遺傳改良，提高河南省奶牛質(zhì)量。□

圖8 河南省奶牛遺傳改良圖示