任小麗，栗敏杰，白雪利，劉 歡，王雅春，閆 磊,4，閆躍飛,4，蘇銀池，張 震*

（1.河南省奶牛生產(chǎn)性能測定中心，河南鄭州 450000；2.河南農(nóng)業(yè)大學牧醫(yī)工程學院，河南鄭州 450000；3.中國農(nóng)業(yè)大學動物科學技術學院，北京 100193；4.河南乾祥農(nóng)牧科技有限公司，河南鄭州 450000；5.河南省種牛遺傳性能測定中心，河南鄭州 450000）

奶牛高產(chǎn)性狀的高強度選育工作取得了顯著效果，同時造成奶牛的健康水平、長壽性和繁殖性能下降，生產(chǎn)成本逐步增加[1]。2003年，低繁殖力成為英國奶牛群體非自然淘汰的主要原因[2]。繁殖性狀與臨床乳房炎[3]、飼料轉化效率[4]、體況評分[5]、產(chǎn)奶量和尿素氮[6]存在不同程度的遺傳相關，是影響奶牛養(yǎng)殖效率的重要性狀。近年來，隨著平衡育種理念的發(fā)展，越來越多的國家逐漸將繁殖性狀納入到育種目標中[7-8]。

奶牛繁殖性狀中，青年牛初次產(chǎn)犢年齡能反映母牛個體的性成熟早晚及其受孕能力[9]，最佳的初次產(chǎn)犢年齡保證個體能更早地獲得生產(chǎn)效益并適當提高利用年限[10]；研究表明，青年牛初次產(chǎn)犢年齡與成年母牛產(chǎn)犢間隔存在相關關系[10]，且高產(chǎn)奶牛存在產(chǎn)犢間隔增加、繁殖性能衰退的趨勢[11]；娟珊牛和荷斯坦牛群體的產(chǎn)犢間隔與產(chǎn)后初配天數(shù)、懷孕、泌乳期存在正相關，在荷斯坦牛群體中與使用年限有負相關關系[12]。

奶牛重要經(jīng)濟性狀遺傳參數(shù)具有群體特異性。本研究收集河南省17個牧場2008—2017年奶牛生產(chǎn)性能記錄獲得14 406頭青年牛初次產(chǎn)犢年齡記錄與8 912頭成年母牛產(chǎn)犢間隔記錄，使用多性狀免求導（Derivative Free Multivariate，DMU）軟件DMU_AI模塊結合動物模型估計遺傳參數(shù)，為區(qū)域性、精細化育種目標的制定提供理論參考。

1 材料與方法

1.1 數(shù)據(jù)收集與整理 收集2008—2017年河南省飼養(yǎng)規(guī)模500頭左右、參加奶牛生產(chǎn)性能測定的17個牧場的19 228頭共341 582條生產(chǎn)性能測定記錄，包括牛號、父號、母號、出生日期、產(chǎn)犢日期、胎次、牧場信息。

青年牛初次產(chǎn)犢年齡：頭胎牛只產(chǎn)犢日期與出生日期相差天數(shù)（d），質(zhì)控標準為550～1 100 d。獲得14 406頭青年牛初次產(chǎn)犢年齡記錄。

成年母牛產(chǎn)犢間隔：連續(xù)參測牛只（胎次連續(xù)）相鄰胎次間產(chǎn)犢日期間相差天數(shù)（d），質(zhì)控標準為280～600 d。獲得8 912頭成年母牛產(chǎn)犢間隔記錄。

表1 各牧場青年牛初次產(chǎn)犢年齡和成年母牛產(chǎn)犢間隔數(shù)據(jù)

表2 成年母牛各胎次間隔中的成年母牛產(chǎn)犢間隔變化

質(zhì)控后數(shù)據(jù)記錄數(shù)分布見表1和表2。

1.2 遺傳評估模型 本研究使用動物模型，借助DMU 軟件DMU_AI模塊采用AI-REML（AI：Average Inforlihood，約束最大似然法）結合期望最大化（Expectation Maximization，EM）算法估計初次產(chǎn)犢年齡和產(chǎn)犢間隔遺傳參數(shù)。模型通式：

其中，y為觀察值向量（初次產(chǎn)犢年齡和產(chǎn)犢間隔）；β為固定效應向量，初次產(chǎn)犢年齡模型中固定效應包括出生年季、產(chǎn)犢年季和牧場效應，產(chǎn)犢間隔模型中固定效應包括出生年季、產(chǎn)犢年季、牧場和胎次效應；α為隨機效應向量，在初次產(chǎn)犢年齡和產(chǎn)犢間隔模型中為加性效應；e為隨機殘差向量。X和Z分別為結構矩陣。

模型前提假設：

其中，σα表示加性效應方差，σe表示殘差效應方差，σp表示表型方差， h2表示遺傳力。

遺傳評估采用的系譜為大群系譜共5 962頭公牛、21 6947母牛，包含本研究中16 443頭母牛一代、二代、三代母系系譜比例為71.91%、45.82%、16.27%。個體育種值為DMU軟件計算得到的加性效應最佳線性無偏預 測（Best Linear Unbiased Prediction，E_BLUP） 估計值。

根據(jù)河南省溫濕度變化劃分為以下4個季度：第一季度為3、4、5月，第二季度為6、7、8月，第三季度為9、10、11月，第四季度為12、1、2月[13]。

2 結果與分析

2.1 性狀描述性統(tǒng)計量 整理后共有14 406頭（14 406條記錄）青年牛初次產(chǎn)犢年齡數(shù)據(jù)，平均值為794.65 d（標準差100.60 d），最大值為1 100 d，最小值為551 d；8 912頭（14 787條記錄）成年母牛14 787條產(chǎn)犢間隔記錄，平均值為415.17 d（標準差73.20 d），最大值為600 d，最小值為280 d。

2.2 青年牛初次產(chǎn)犢年齡隨出生年和產(chǎn)犢年的變化趨勢 由圖1可知，2004年出生的青年牛初次產(chǎn)犢年齡為906.00 d，2015年降至710.10 d，2008—2017年青年牛初次產(chǎn)犢年齡由857.28 d降至767.68 d。

圖1 青年牛初次產(chǎn)犢年齡隨出生年份（上）和產(chǎn)犢年份（下）的變化趨勢

2.3 成年母牛產(chǎn)犢間隔隨出生年和產(chǎn)犢年的變化趨勢由圖2可知，2000年出生的成年母牛多胎次間的產(chǎn)犢間隔為400.86 d，2014年降至367.82 d；成年母牛產(chǎn)犢間隔隨著產(chǎn)犢年份沒有很明顯變化。

圖2 成年母牛產(chǎn)犢間隔隨出生年份（上）和產(chǎn)犢年份（下）的變化趨勢

2.4 青年牛初次產(chǎn)犢年齡和成年母牛產(chǎn)犢間隔的方差組分和遺傳力估計結果 青年牛初次產(chǎn)犢年齡加性效應方差為1 152.90 d，殘差效應為4 330.07 d，遺傳力為0.21；成年母牛產(chǎn)犢間隔加性效應方差為229.73 d，殘差效應方差為4 280.99 d，遺傳力為0.05。

2.5 青年牛初次產(chǎn)犢年齡和成年母牛產(chǎn)犢間隔的遺傳趨勢變化 由圖3可知，2002—2015年青年牛初次產(chǎn)犢年齡和成年母牛產(chǎn)犢間隔育種值均呈下降趨勢，青年牛初次產(chǎn)犢年齡育種值下降趨勢更大，由0.30下降到-4.84，成年母牛產(chǎn)犢間隔育種值由0.77下降到-0.01。

圖3 青年牛初次產(chǎn)犢年齡和成年母牛產(chǎn)犢間隔的遺傳趨勢

3 討 論

適當降低初次產(chǎn)犢年齡可以增加荷斯坦牛的繁殖力和產(chǎn)奶量，但同時也會增加死胎率。美國荷斯坦牛和瑞士褐牛群體中關于初次產(chǎn)犢年齡的研究建議最佳初次產(chǎn)犢年齡為21～22月齡（約630～660 d）[14]。本研究中青年牛初次產(chǎn)犢年齡為794.65 d ，與北京地區(qū)中國荷斯坦牛（798.47 d）[1]、墨西哥荷斯坦牛群體（798 d）[15]相近，低于新疆地區(qū)（814.4 d）[16]、中國荷斯坦牛的803.17 d[17]；且隨著出生年份增加，青年牛初次產(chǎn)犢年齡逐漸降至2105年的710.10 d（約24月齡），達到奶牛較適宜的初次產(chǎn)犢年齡范圍。

在初次產(chǎn)犢年齡和產(chǎn)犢間隔遺傳參數(shù)估計模型研究中，Sun等[18]比較公畜模型、公母畜模型和動物模型估計丹麥荷斯坦牛群體繁殖性狀遺傳參數(shù)效果，表明動物模型更適合用于繁殖性狀遺傳參數(shù)估計。Brito等[19]研究表明，雙閾值-線性不考慮缺失信息的模型更適合用于肉用牛初次產(chǎn)犢年齡性狀的遺傳評估。de Araujo Neto等[20]采用多性狀互作標準模型（Multi-trait Reaction Norm Models）中加入基因與環(huán)境互作效應，估計得到內(nèi)絡爾牛初次產(chǎn)犢年齡遺傳參數(shù)為0.222～0.316。另有研究表明，同一繁殖性狀不同胎次的遺傳力不同，如Liu等[21]關于中國荷斯坦牛、Wasike等[22]關于博蘭牛的研究表明，在第1、2、3、4胎次，產(chǎn)犢間隔遺傳力分別為0.00、0.15、0.00、0.00。本研究對初次產(chǎn)犢年齡和產(chǎn)犢間隔均采用動物模型進行遺傳參數(shù)估計，不同胎次間的產(chǎn)犢間隔按照同一性狀對待，遺傳力結果與Wasike等[22]研究結果相近。

本研究中，青年牛初次產(chǎn)犢月齡遺傳力計算結果為0.21，這與荷斯坦牛相關研究的0.28[15]、0.25[17]、相近，高于北京地區(qū)中國荷斯坦牛0.066 4[1]、新疆地區(qū)荷斯坦牛群體0.07[16]。初次產(chǎn)犢年齡在不同品種間遺傳力估計結果不同，但均屬于低遺傳力性狀，如采用多性狀動物模型估計博蘭牛為0.04[22]、3 955頭瑞士褐牛為0.08[23]、用REML的動物模型估計355頭娟珊牛雜交后代為0.220[24]、乳肉兼用牛為0.03[25]、內(nèi)絡爾牛為0.13[26]。本研究中，青年牛初次產(chǎn)犢年齡性狀表型值是根據(jù)個體出生日期與頭胎產(chǎn)犢日期的間隔天數(shù)計算得來，這或許是造成遺傳力結果有差異的原因之一。

本研究中，8 912頭成母牛產(chǎn)犢間隔平均值為415.17 d，不同于其他研究中荷斯坦奶牛的397.32 d[1]、417.8 d[16]。與初次產(chǎn)犢年齡性狀趨勢一致，隨著生產(chǎn)管理水平的提高成年母牛產(chǎn)犢間隔降低至2014年的367.82 d，也接近奶牛適宜的產(chǎn)犢間隔（365 d）[27]。遺傳參數(shù)估計結果表明，成年母牛產(chǎn)犢間隔遺傳力結果為0.05，這與相關研究的0.054 9[1]、0.04[28]相近。而Wall等[2]用多性狀公畜-外祖父模型估計英國荷斯坦牛群體產(chǎn)犢間隔遺傳力為0.033。

4 結 論

本研究使用DMU軟件DMU_AI模塊采用AIREML結合EM算法的動物模型估計得到青年牛初次產(chǎn)犢年齡的遺傳力為0.21、成年母牛產(chǎn)犢間隔性狀遺傳力為0.05，為區(qū)域性、精細化育種目標的制定奠定理論基礎。