常 瑤，朱瑩琳，施佳慶，安 濤，于亞冬，張海亮，徐 偉，張寶石，王志剛，王雅春*

(1. 中國農(nóng)業(yè)大學(xué)動物科學(xué)技術(shù)學(xué)院，北京 100193； 2. 大同市四方高科農(nóng)牧有限公司，大同 037001；3. 全國畜牧總站，北京 100125)

隨著奶牛生產(chǎn)水平的提高及長期的高度選育，高產(chǎn)荷斯坦牛的健康狀況、繁殖力等出現(xiàn)了下降的趨勢[1]。西門塔爾牛在世界上分布廣、數(shù)量多，具有抗逆性強(qiáng)、乳品質(zhì)優(yōu)良、繁殖率高等特點(diǎn)，是優(yōu)秀的乳肉兼用品種。因此在部分國家和地區(qū)，奶牛養(yǎng)殖者開始嘗試?yán)萌槿饧嬗眯臀鏖T塔爾牛和荷斯坦牛雜交(后代簡稱“西荷?！?繁育模式來提高牛群抗逆性和繁殖力以達(dá)到提高綜合經(jīng)濟(jì)效益的目的[2]。在泌乳性能方面，R.E.McDowell等[3]和G.Ahlborn-Breier等[4]研究表明，愛爾夏、瑞士褐牛、娟姍牛與荷斯坦牛雜種后代的乳成分含量高于荷斯坦牛；王志耕等[5]將西門塔爾牛和當(dāng)?shù)攸S牛雜交品種與荷斯坦牛進(jìn)行比較，結(jié)果表明雜交牛的乳脂量、乳蛋白量有顯著的提高；在繁殖力方面，B.J.Heins等[6]和孫少華等[7]研究發(fā)現(xiàn)，蒙貝利亞(法系西門塔爾牛)和荷斯坦牛的雜種后代空懷天數(shù)和初產(chǎn)月齡顯著低于荷斯坦牛。但目前國內(nèi)外對于西荷雜交效果的研究還不夠完善，只對泌乳、繁殖及生長性狀有所涉及，缺乏對其它重要經(jīng)濟(jì)性狀的綜合比較。

本研究就山西省某大規(guī)模奶牛場的法系西荷雜種牛、德系西荷雜種牛、自繁荷斯坦牛以及進(jìn)口澳系荷斯坦牛4個群體，除了比較前人研究中關(guān)注的泌乳、繁殖等重要經(jīng)濟(jì)性狀差異外，還就熱應(yīng)激相關(guān)性狀、生長、體況、步態(tài)、性情評分及皮膚褶皺厚度等性狀的測定結(jié)果進(jìn)行比較，以期對西荷雜交繁育體系的效果進(jìn)行更深入的探討與分析。

1 材料與方法

1.1 資料來源

本試驗(yàn)選取山西省某奶牛場為試驗(yàn)牛場。調(diào)查時該場存欄泌乳牛4 393頭，其中德系西荷雜種牛21頭，法系西荷雜種牛34頭，自繁荷斯坦牛2 428頭和進(jìn)口澳系荷斯坦牛1 910頭。試驗(yàn)測量數(shù)據(jù)由中國農(nóng)業(yè)大學(xué)牛百科團(tuán)隊(duì)6名隊(duì)員于2016年7月18日—7月28日集中測定，并收集了2012年—2016年間2 699頭奶牛的配種、產(chǎn)犢等繁殖記錄及2015年—2016年該場生產(chǎn)性能測定數(shù)據(jù)(Dairy herd improvement，DHI)。

1.2 測定項(xiàng)目

1.2.1 泌乳性狀 收集DHI數(shù)據(jù)中日產(chǎn)奶量、乳脂率、乳蛋白率、體細(xì)胞數(shù)(Somatic cell count，SCC)以及尿素氮指標(biāo)，依據(jù)公式SCS=log2(SCC/100)+3[8]，將SCC轉(zhuǎn)換成體細(xì)胞評分(Somatic cell score，SCS)，并計算出校正奶量(Herd test adjusted corrected milk，HTACM)，公式如下：

HTACM=((0.432×DMY)+(16.23×(DMY×F%))+((DIM-150)×0.002 9)×DMY)×PCF，式中，DMY為日產(chǎn)奶量，F(xiàn)%為乳脂率，DIM為泌乳天數(shù)，PCF為胎次校正系數(shù)，本研究泌乳性狀只選取頭胎牛，校正系數(shù)為1.064[9]。

1.2.2 繁殖性狀 本研究中針對6個繁殖性狀進(jìn)行分析：(1)青年牛首次配種日齡(Age at first service in heifer，AFSh)；(2)青年牛配妊次數(shù)(Number of services in heifer，NSh)；(3)青年牛妊娠期長度(Gestation length in heifer，GLh)；(4)成年母牛配妊次數(shù)(Number of services in cow，NSc)；(5)成年母牛產(chǎn)后首次配種時間(Calving to first service in cow，CTFSc)；(6)成年母牛空懷期(Days open in cow，DOc)。

1.2.3 生長發(fā)育及其他性狀 生長發(fā)育性狀測定項(xiàng)目包括體高、體斜長、胸圍和體深(最后一根肋骨處腰椎至腹底部的垂直距離,cm)，并采用約翰遜公式[10]估計體重，即：

體重=胸圍2×體斜長/10 800。

選取部分奶牛群體進(jìn)行體況、性情和步態(tài)評分，并測量皮膚褶皺厚度、直腸溫度及相應(yīng)環(huán)境溫濕度。

1.3 測定方法

1.3.1 體況評分(Body condition score，BCS) 依據(jù)E.E.Wildman等[11]和A.J.Edmonson等[12]的評分標(biāo)準(zhǔn)，通過視覺評估和觸摸判斷，結(jié)合整體印象進(jìn)行評分。每頭牛由兩人獨(dú)立進(jìn)行評分，并取兩人評分的平均值進(jìn)行結(jié)果分析。

1.3.2 直腸溫度 使用電子體溫計(OMRON)于上午07：00—11:00和下午14:30—19:00各測定一次直腸溫度。測量時體溫計在直腸內(nèi)停留時間10 s左右，插入深度10 cm左右，讀數(shù)精確到0.1 ℃。記錄測定直腸溫度前后的環(huán)境溫濕度、風(fēng)扇開關(guān)狀態(tài)及是否為擠奶前。溫濕度指數(shù)(Temperature humidity index，THI)計算參照董曉霞等[13]的方法。

1.3.3 性情評分 觀察并評價首次直腸溫度測定全過程奶牛的反應(yīng)進(jìn)行性情評分。評分原則：1分—一直很安靜，基本無反應(yīng)；2分—后軀左右擺動，但反應(yīng)較??；3分—有較大排斥反應(yīng)，甚至有踢人動作。

1.3.4 皮膚褶皺厚度 由一名測定人員使用電子游標(biāo)卡尺(型號：part#PD-151)分別在母牛右側(cè)頸部和肋部最后一節(jié)肋骨處測量得到。

1.3.5 步態(tài)評分(Locomotion score，LS) 步態(tài)評分標(biāo)準(zhǔn)主要根據(jù)自然站立時背部表現(xiàn)和行走表現(xiàn)綜合評定[14]，評分為1～5分。

1.4 數(shù)據(jù)處理

1.4.1 數(shù)據(jù)整理 結(jié)合在場牛只基本信息、DHI、繁殖記錄及各性狀測定數(shù)據(jù)，刪除：(1)信息無法匹配和記錄不完整的個體數(shù)據(jù)；(2)結(jié)合相關(guān)文獻(xiàn)[15]，不符合生理規(guī)律以及記錄或測量有誤的數(shù)據(jù)，其中繁殖數(shù)據(jù)的篩選標(biāo)準(zhǔn)見表1；(3)經(jīng)四分位數(shù)檢驗(yàn)法檢驗(yàn)判斷為異常值者；(4)德系西荷雜種牛產(chǎn)奶、繁殖性狀數(shù)據(jù)(頭數(shù)太少)；(5)泌乳天數(shù)大于330 d的牛只數(shù)據(jù)。

表1 繁殖數(shù)據(jù)篩選標(biāo)準(zhǔn)

Table 1 The standard for reproduction data filtering

項(xiàng)目Item篩選條件Filtercondition胎次Parity≤1青年牛首次配種日齡/dAFSh[270,900]青年牛或成年母牛妊娠期長度/dGLh[260,310]青年?；虺赡昴概Ｅ淙汛螖?shù)NSc[1,10]成年母牛產(chǎn)后首次配種天數(shù)/dCTFSc[20,250]成年母牛空懷天數(shù)/dDOc[20,400]

“[ ]”表示閉區(qū)間，即篩選條件為≥最小值，≤最大值

“[ ]”means closed interval，i.e.，the values≥ minimum and ≤ maximum are used as filter conditions

按照氣象信息季節(jié)劃分四季。為確保與雜種牛進(jìn)行對照，產(chǎn)奶性狀和步態(tài)評分選取頭胎牛，生長性狀和體況評分含1胎、2胎和3胎，其余性狀選取1胎和2胎兩個水平。泌乳階段劃分為4個階段：1～50 d(階段Ⅰ)、51～100 d(階段Ⅱ)、101～200 d(階段Ⅲ)及201～330 d(階段Ⅳ)。

1.4.2 統(tǒng)計方法 直腸溫度數(shù)據(jù)采用SAS9.2 Mixed過程，其余性狀均采用GLM過程分析，用Bonferronit檢驗(yàn)進(jìn)行多重比較，顯著性水平α=0.05。

對奶牛泌乳性狀分析采用固定模型1：

Yijkl=μ+YEARi+SEASONj+BREEDk+βDIM+eijkl，式中，Yijkl為頭胎奶牛的泌乳性狀(日產(chǎn)奶量、乳脂率、乳蛋白率、校正奶量、體細(xì)胞評分及尿素氮)，μ為總體平均值，YEARi為DHI采樣年份，SEASONj為采樣季節(jié)，BREEDk為品種效應(yīng)，DIM為泌乳天數(shù)，β為DIM的對應(yīng)回歸系數(shù)，eijkl為隨機(jī)殘差。

對青年牛繁殖性狀分析采用固定模型2：

Yijklm=μ+BREEDi+BYSj+SYS1k+SYSnl+eijklm，式中，Yijklm為各青年牛繁殖性狀的觀察值(AFSh，NSh，GLh)，μ為總體平均值，BREEDi為品種效應(yīng)，BYSj為出生年季，SYS1k和SYSnl為首次和末次配種年季，eijklm為隨機(jī)殘差。

對成母牛繁殖性狀分析采用固定模型3：

Yijklmn=μ+BREEDi+BYSj+SYS1k+SYSnl+CYSpm+eijklmn，式中，Yijklmn為各成母牛繁殖性狀的觀察值(NSc，CTFSc，DOc)，CYSpm為前次產(chǎn)犢年季，eijklmn為隨機(jī)殘差，其他效應(yīng)含義同前。

對奶牛直腸溫度分析采用混合模型4：

Yijkm=μ+BREEDi+PARITYj+STAGEk+β1BCS+β2THI+Zτ+eijkm，式中，Yijkm為奶牛上午或下午直腸溫度，PARITYj為胎次效應(yīng)，STAGEk為泌乳階段，BCS為協(xié)變量體況評分，THI為協(xié)變量溫濕度指數(shù)，β1、β2為回歸系數(shù)，Z為個體隨機(jī)效應(yīng)的設(shè)計矩陣，τ為個體隨機(jī)效應(yīng)，eijkm為隨機(jī)殘差，其他效應(yīng)含義同前。

其余性狀所選用的固定效應(yīng)見表2。

表2 各性狀分析模型中的固定效應(yīng)分布

Table 2 The fixed effects for analysis model of each trait

性狀Trait品種Breed胎次Parity泌乳階段Lactationstage泌乳天數(shù)Daysinmilk體況評分BSC體高、體斜長Bodyheightandslantinglength√√胸圍、體深、估算體重Chestcircumference,bodydepthandestimatedweight√√√體況評分、性情評分BCSandtemperamentscore√√√步態(tài)評分Locomotionscore√√皮膚褶皺厚度Skinfoldthickness√√√√

“√”表示該性狀的分析模型中包含對應(yīng)的固定效應(yīng)

“√” means the fixed effects for analysis model of each trait

2 結(jié) 果

2.1 西荷雜交牛與荷斯坦牛泌乳性狀對比

本研究中不同品種奶牛群體的泌乳性狀比較結(jié)果見表3。由表3可知，法系西荷牛的日產(chǎn)奶量和校正奶量均高于澳系和自繁荷斯坦牛，并與澳系荷斯坦牛有顯著差異(P<0.05)。在乳成分的比較上，法系西荷牛測定日乳脂率、乳蛋白率均為最高，分別高于自繁荷斯坦牛0.23和0.12個百分點(diǎn)，且乳蛋白率的差異達(dá)到顯著水平(P<0.05)。在體細(xì)胞評分上，法系西荷牛測定日SCS低于自繁荷斯坦牛，但兩品種差異不顯著(P>0.05)。3個品種群牛奶尿素氮含量無顯著差異(P>0.05)。

2.2 西荷雜交牛與荷斯坦牛繁殖性狀對比

由表4可知，法系西荷青年牛首次配種日齡顯著低于自繁荷斯坦牛和澳系荷斯坦牛(P<0.05)，分別提前18.09和24.43 d。成年母牛產(chǎn)后首次配種時間最短的為澳系荷斯坦牛，平均需要59.43 d的恢復(fù)期后即可進(jìn)行此泌乳期的第一次配種，顯著小于自繁荷斯坦牛(P<0.05)。法系西荷雜交牛產(chǎn)后到首次配種時間小于自繁荷斯坦牛6.66 d，但兩者無顯著差異(P>0.05)。在配妊次數(shù)上，澳系荷斯坦青年牛平均需要進(jìn)行1.31次配種妊娠，低于法系西荷牛和自繁荷斯坦牛，且與自繁荷斯坦牛有顯著差異(P<0.05)。經(jīng)產(chǎn)牛配妊次數(shù)整體多于青年牛，法系西荷成母牛平均所需次數(shù)最少(1.88次)。在空懷天數(shù)上，法系西荷牛均值比自繁和澳系荷斯坦牛分別少了38.43和20.93 d。

表3 不同品種奶牛群體泌乳性狀的比較(LSM±SE)

Table 3 Comparison on milk traits among different breeds of dairy groups(LSM±SE)

品種(數(shù)量)Breed(number)日產(chǎn)奶量/kgDailymilkyield校正奶量/kgHTACM乳蛋白率/%Milkproteinpercentage乳脂率/%Milkfatpercentage體細(xì)胞評分SCS尿素氮/(mg·dL-1)MUN法系西荷雜種牛MO×HO(N=86)29.91±0.93a32.86±1.47a3.50±0.03a3.64±0.11ab2.30±0.18ab14.56±0.35自繁荷斯坦牛HO(N=1921)28.54±0.20a30.23±0.24a3.38±0.01b3.41±0.02a2.42±0.04a14.80±0.08澳系荷斯坦牛AHO(N=2540)27.12±0.22b29.24±0.26b3.42±0.01a3.53±0.03b2.10±0.04b14.59±0.08

N. 數(shù)據(jù)量；同一因子的同列數(shù)據(jù)肩標(biāo)不同表示差異顯著(P<0.05)，相同表示差異不顯著(P>0.05)。下表同

N. Number; Different letters superscripts under the same factor within the same column means significant difference (P<0.05), same letter superscript means not significant difference(P>0.05). The same for the following tables

表4 不同品種奶牛群體繁殖性狀的比較

Table 4 Comparison on reproduction traits among different breeds of dairy groups

繁殖性狀Reproductiontrait法系西荷雜種牛MO×HO自繁荷斯坦牛HO澳系荷斯坦牛AHONLSM±SE/M±SDNLSM±SE/M±SDNLSM±SE/M±SD青年牛首次配種日齡/dAFSh27382.49±7.48b1129400.58±1.52c1174406.92±2.26a青年牛配妊次數(shù)NSh271.54±0.21ab10301.52±0.05b12651.31±0.07a青年牛妊娠期長度/dGLh26276.79±0.63733275.91±0.82786276.79±0.63成年母牛配妊次數(shù)NSc181.88±1.13ab4832.32±1.45b7682.26±1.53a成年母牛產(chǎn)后首次配種時間/dCTFSc2166.90±4.63ab90473.56±1.02b136259.43±2.05a成年母牛空懷天數(shù)/dDOc1895.94±42.75ab473134.37±79.20b756116.87±75.72a

M. 算數(shù)平均值；SD. 標(biāo)準(zhǔn)差。 NSc和DOc采用“算數(shù)平均數(shù)±標(biāo)準(zhǔn)差”表示，其余性狀采用“最小二乘均值±標(biāo)準(zhǔn)誤”

M. Means; SD. Standard deviation. The results of NSc and DOc are represented by “M±SD”, others are represented by “LSM±SE”

2.3 西荷雜交牛與荷斯坦牛生長發(fā)育及其他性狀對比

2.3.1 生長發(fā)育性狀及體況評分 對德系西荷牛、法系西荷牛和自繁荷斯坦牛3個群體的生長發(fā)育及體況評分指標(biāo)分析結(jié)果見表5。除體深外，其余性狀雜種牛均高于自繁荷斯坦牛。法系和德系西荷牛的胸圍分別大于自繁荷斯坦牛4.49和1.31 cm，且法系西荷牛與自繁荷斯坦牛差異顯著(P<0.05)。法系西荷牛體況評分最高，德系西荷牛其次，雜種牛的BCS顯著高于自繁荷斯坦牛(P<0.05)。

2.3.2 直腸溫度 圖1給出了試驗(yàn)測定期間牛場牛舍內(nèi)上午和下午溫濕度指數(shù)折線圖。在試驗(yàn)期間牛舍內(nèi)上午THI范圍為68.60～73.29，平均值為71.59，下午THI范圍為72.25～78.19，平均值為74.81。本研究中388頭奶牛上午直腸溫度平均為(38.83±0.45)℃，下午直腸溫度平均為(39.11±0.47)℃。

由表5可知，澳系荷斯坦牛上、下午的直腸溫度均為最高。雜種牛的上午直腸溫度整體低于荷斯坦牛，德系雜交牛直腸溫度最低(38.66 ℃)，在奶牛正常溫度38.3～38.7 ℃[16]范圍內(nèi)，且顯著低于澳系荷斯坦牛(P<0.05)。各品種的下午RT無顯著差異(P>0.05)。

2.3.3 性情評分 德系和法系西荷牛性情評分平均為1.14分，低于自繁和澳系荷斯坦牛，但各品種之間無顯著差異(P>0.05)。

表5 不同品種奶牛群體生長發(fā)育及其他性狀的比較(LSM±SE)

Table 5 Comparison on growth and other traits among different breeds of dairy groups (LSM±SE)

性狀Trait德系西荷雜種牛FL×HO法系西荷雜種牛MO×HO自繁荷斯坦牛HO澳系荷斯坦牛AHO體尺Bodysize數(shù)據(jù)量Number162498-體高/cmBodyheight148.39±0.85147.11±0.81147.10±0.45-體斜長/cmBodyslantinglength171.10±1.58171.35±1.50170.62±0.82-胸圍/cmChestcircumference208.33±1.95ab211.51±1.88a207.02±1.07b-體深/cmBodydepth85.52±0.7986.17±0.7687.23±0.43-估算體重/kgEstimatedweight689.86±15.12713.17±14.58681.90±8.27-體況Bodycodition數(shù)據(jù)量Number1014520-評分Score2.91±0.14a2.93±0.12a2.55±0.03b-直腸溫度/℃RT數(shù)據(jù)量Number111949047上午Morning38.66±0.14a38.72±0.11ab38.85±0.02ab38.97±0.07b數(shù)據(jù)量Number141859550下午Afternoon39.18±0.1439.16±0.1239.09±0.0239.25±0.08性情Temperament數(shù)據(jù)量Number81641833評分Score1.14±0.151.14±0.111.22±0.021.19±0.08皮厚/mmSkinfoldthickness數(shù)據(jù)量Number91340224肋部Thelastrib11.70±0.39a10.03±0.33bc10.61±0.07b9.69±0.25c頸部Neck9.10±0.42ab8.18±0.36ab8.73±0.08a7.97±0.27b步態(tài)數(shù)據(jù)量Number71620531Locomotion評分Score1.43±0.201.46±0.141.53±0.041.46±0.10

-. 代表該數(shù)據(jù)未測量

-. Means the data was not measured

圖1 牛場牛舍內(nèi)上午、下午溫濕度指數(shù)曲線圖Fig.1 The morning and afternoon THI curves in the barn in dairy farm

2.3.4 皮膚褶皺厚度 由表5可知，各品種肋部的皮膚褶皺厚度(簡稱皮褶厚度)明顯厚于頸部。本研究中，不論肋部或頸部，德系西荷牛皮褶厚度最厚，自繁荷斯坦牛和法系雜交牛次之，澳系荷斯坦牛最薄，且德系西荷牛的肋部皮褶厚度顯著高于其他品種(P<0.05)。

2.3.5 步態(tài)評分 由表5可知，德系西荷牛的步態(tài)評分均值最低，自繁荷斯坦牛的評分均值最高，但不同品種間差異不顯著(P>0.05)。各品種步態(tài)評分的比例分布見圖2。根據(jù)圖2可知，除自繁荷斯坦牛的評分有3分(0.98%)和4分(0.49%)外，其余均只有1分和2分。雜種牛的1分比例相對較高，2分中自繁荷斯坦牛的比例最高。

圖2 不同品種奶牛群體步態(tài)評分比例圖Fig.2 The LS distribution among different breeds of dairy groups

3 討 論

3.1 不同奶牛群體泌乳性狀比較

產(chǎn)奶量具有較高的經(jīng)濟(jì)價值，而隨著奶業(yè)市場的不斷發(fā)展，乳成分相關(guān)性狀也逐漸受到重視[2]。雜種后代在日產(chǎn)奶量、乳脂率、乳蛋白率和校正奶量等性狀具有一定的優(yōu)勢。其中荷斯坦牛與雜種牛比較的結(jié)果與B.J.Heins等[17-18]的研究結(jié)果不一致。本研究中的自繁荷斯坦牛生產(chǎn)性能較低，而法系西門塔爾牛其產(chǎn)奶量在法國僅次于荷斯坦[19]，自繁荷斯坦牛使用法系西門塔爾牛凍精產(chǎn)生的雜種后代，綜合遺傳互補(bǔ)和雜種優(yōu)勢效應(yīng)，產(chǎn)奶量水平較高?，F(xiàn)有研究結(jié)果表明，純種荷斯坦牛的產(chǎn)奶量與各類雜交牛相比具有明顯優(yōu)勢[2]，但少數(shù)特定組合雜種牛的產(chǎn)奶量表現(xiàn)亦不遜色，如A.R.Hazel等[20]在對處于第一個泌乳期的加州奶牛群體的研究中發(fā)現(xiàn)，法系西荷雜種牛的前150 d平均產(chǎn)奶量比同群純種荷斯坦牛略高(+2%)，但無顯著差異。

在乳成分上法系西荷牛優(yōu)勢明顯。L.B.Hansen等[21]對美國加利福尼亞7個牛場的系統(tǒng)研究中，法系西荷牛的日均乳成分產(chǎn)量(乳蛋白量和乳脂量)均顯著高于同群純種荷斯坦牛。何淑靜[22]研究表明，法系西荷牛乳脂率、乳蛋白率分別顯著和極顯著高于純種荷斯坦牛。另外，有研究[23]表明，澳系荷斯坦牛乳脂率和乳蛋白率平均為3.6%和3.4%，這與本研究結(jié)果一致。對于奶牛來講，飼料中的粗纖維經(jīng)瘤胃微生物發(fā)酵產(chǎn)生丁酸和乙酸，再由丁酸和乙酸作為原料進(jìn)一步合成乳脂，因此乳脂量與粗飼料利用率有關(guān)[24]。本研究中，法系西荷牛乳脂率較高，可能由于其粗飼料利用率較高[25]。遺傳是影響乳蛋白量的主要因素[19]，本研究中，雜種牛乳蛋白率高于自繁和澳系荷斯坦牛，可能是由于其繼承了西門塔爾牛乳蛋白含量較高的遺傳優(yōu)勢。

本研究中，法系西荷牛SCS較低，反映了雜交牛的乳房更為健康，這與F.Malchiod等[26]、B.J.Heins和L.B.Hansen[18]的研究結(jié)果相一致。S.M.Godden等[27]認(rèn)為，當(dāng)MUN在12～19 mg·dL-1時(乳蛋白率大于3%)能氮較為平衡。本研究中各品種MUN均處在正常含量范圍內(nèi)。

3.2 不同奶牛群體繁殖性狀比較

奶牛生產(chǎn)的經(jīng)濟(jì)效益很大程度上取決于牛群的繁殖效率，但繁殖性狀遺傳力低，而雜種優(yōu)勢較高，可期望在短時間內(nèi)通過雜交繁育提高繁殖效率[22]。

提前配種可使奶牛提早進(jìn)行生產(chǎn)，增加使用年限，同時可減少飼養(yǎng)青年牛所需的消耗，從而降低飼養(yǎng)成本[15]。何淑靜[22]研究表明，法系西荷母牛12月齡前的日增重顯著高于荷斯坦牛，因此雜種牛初配日齡較早可能是由于其生長速度快，能較早達(dá)到配種體重。毛永江等[28]研究得出，澳系進(jìn)口荷斯坦牛的初配月齡晚于中國荷斯坦牛2～3月左右，馬洪祥等[29]研究發(fā)現(xiàn)，澳系荷斯坦牛稍晚于西門塔爾牛，與本試驗(yàn)中澳系進(jìn)口荷斯坦牛初配日齡最大的結(jié)果相一致，推測這可能與其在原產(chǎn)地的飼養(yǎng)水平及長途運(yùn)輸有一定關(guān)系。

通常情況下，青年牛首次配種日齡稍晚有利于青年牛的體成熟及性成熟和成年母牛產(chǎn)犢后體能的恢復(fù)，因此青年牛首次配種日齡最晚的澳系荷斯坦牛產(chǎn)犢后恢復(fù)最快。但法系西荷雜交牛配種最早，但產(chǎn)后恢復(fù)時間仍短于自繁荷斯坦牛。D.A.Knob等[1]對比荷斯坦和西荷雜種牛的研究中，雜種牛產(chǎn)犢到首次配種的間隔極顯著低于荷斯坦牛(24 d)。結(jié)合本研究結(jié)果表明，與荷斯坦牛相比雜種牛產(chǎn)后恢復(fù)周期較短。L.G.D.Mendon?a等[30]認(rèn)為，分娩后雜種牛的黃體發(fā)育早于荷斯坦牛(28.4 dvs30.2 d)，同時首次排卵和產(chǎn)后第1個發(fā)情期也要早于荷斯坦牛，本研究結(jié)果也進(jìn)一步印證了這一理論。

由于初配日齡較晚的澳系牛性成熟及體成熟均發(fā)育較好，所以配妊次數(shù)少，配種妊娠效率高[14]。F.Malchiod等[31]對奶牛雜交效果的研究表明，不同群體間青年牛配妊次數(shù)差異不顯著，而經(jīng)產(chǎn)牛中法系西荷牛配妊次數(shù)(2.02次)顯著低于純種荷斯坦牛(2.53次)，這與本研究部分結(jié)果一致。

本研究中，雜種牛的空懷天數(shù)最小，與其他同類研究結(jié)果類似。B.J.Heins等[6]研究表明，法系西荷?？諔哑跒?31 d，顯著低于荷斯坦牛(150 d)。何淑靜[22]的研究同樣表明，法系西荷F1代母?？諔烟鞌?shù)(87.35 d)極顯著低于純種荷斯坦牛(20 d)。這可能與雜種母牛泌乳盛期飼料轉(zhuǎn)化率高，能更好的維持膘情，加快子宮復(fù)原有關(guān)[22]。

3.3 不同奶牛群體生長發(fā)育及其他性狀比較

3.3.1 生長發(fā)育性狀及體況評分 邢震全等[32]和何淑靜[22]報道了各生長階段德系和法系西荷牛體重均高于純種荷斯坦牛，其育肥產(chǎn)肉性能明顯提升。傅春泉等[33]對德系西荷牛與荷斯坦牛成年牛(36月齡)生長發(fā)育的研究表明，德系西荷牛體高顯著低于荷斯坦牛，體重、體斜長和胸圍均高于荷斯坦牛，除體高外的其他性狀結(jié)果與本研究一致。

體況評分是評價奶牛營養(yǎng)狀況和飼料好壞的指標(biāo)之一，可反映出奶牛的健康水平。雜種牛體況評分顯著高于自繁荷斯坦牛，說明西荷牛體脂沉積狀況優(yōu)于純種荷斯坦牛，這與相關(guān)文獻(xiàn)報道基本一致[1,30,34]。而體脂沉積狀況優(yōu)良，可更好地維持膘情，有利于減少能量負(fù)平衡，保持良好的泌乳持續(xù)力[22]。

本研究結(jié)果表明，雜種牛體型大、健壯，身軀具有粗而寬，膘情好的特點(diǎn)，與純種荷斯坦牛相比有更高的肉用價值，表現(xiàn)出乳肉兼用型西門塔爾牛與荷斯坦牛雜交在生長性能方面的良好效果。另外，綜合胸圍和體深2個性狀，可以推測雜種牛具有胸廓開張程度較大、肺活量大和血液循環(huán)旺盛等特點(diǎn)，需進(jìn)一步擴(kuò)大樣本量進(jìn)行驗(yàn)證。

3.3.2 直腸溫度 溫濕度指數(shù)是結(jié)合周圍氣體的溫度和相對濕度來綜合評價環(huán)境的一個指標(biāo)[35-36]。根據(jù)G.R.Monfore[37]結(jié)合THI判定熱應(yīng)激程度，7月21日和26日上午(THI<72)，奶牛未產(chǎn)生熱應(yīng)激；其余測量時間 (72≤THI<79)奶牛處于輕微熱應(yīng)激狀態(tài)。直腸溫度可作為評定牛耐熱性的生理指標(biāo)[38]。雜種牛的上午直腸溫度較低，說明耐熱性較好。傅春泉等[39]研究表明，在高溫環(huán)境下德系西荷牛夏季代謝率低，熱應(yīng)激不明顯。張林等[40]研究表明，西荷牛比荷斯坦牛表現(xiàn)出更好的耐熱性，這與本研究結(jié)果基本一致。澳系荷斯坦牛對高溫氣候具有較強(qiáng)的適應(yīng)性[41]，但本試驗(yàn)結(jié)果卻剛好相反，可能是由于該群體引入時間短、體況較差[18]，且C.B.Tucker等[42]研究認(rèn)為體況較差的奶牛體溫較高。

3.3.3 性情評分 本研究結(jié)果表明西荷牛性情較溫順。李英超等[43]研究表明，西荷F1代犢牛性情溫順，不善跑動，喜臥。澳牛生性膽小、易驚，警惕性高[25]，但本研究中澳牛比自繁荷斯坦牛性情評分低，需要進(jìn)一步擴(kuò)大樣本以證實(shí)。

3.3.4 皮膚褶皺厚度 各品種肋部的皮褶厚度厚于頸部，這與張海亮等[44]研究相一致。國外研究表明皮褶厚度受品種因素影響[45]。本研究中，德系西荷牛肋部和頸部的皮褶厚度最厚。這可能是由于德系西門塔爾是乳肉兼用品種，脂肪沉積能力強(qiáng)、較厚，與荷斯坦牛雜交的后代仍保持這一優(yōu)良性狀；文獻(xiàn)報道澳系牛具有皮薄的特點(diǎn)[23]，與本研究結(jié)果一致。

3.3.5 步態(tài)評分 步態(tài)評分可反映出奶牛肢蹄的健康狀況，步態(tài)評分過高意味患蹄病的風(fēng)險較大[46]。結(jié)合表5和圖2可知，與荷斯坦牛相比，雜種牛步態(tài)評分較低、肢蹄健康狀況較好。何淑靜等[22]對蹄病的研究初步說明雜種牛在肢蹄病上的優(yōu)勢。而步態(tài)評分與奶牛的產(chǎn)奶量、繁殖性能等呈負(fù)相關(guān)[47]，因此西荷牛步態(tài)評分較低，有利于降低肢蹄病造成的相關(guān)經(jīng)濟(jì)損失。

4 結(jié) 論

與荷斯坦牛相比，西荷雜種牛的乳成分具有顯著優(yōu)勢；體細(xì)胞評分低，抗病力強(qiáng)；在初配日齡和空懷期等繁殖性狀上有明顯優(yōu)勢；體型大、生長發(fā)育良好、皮膚褶皺厚度較厚，體脂沉積狀況優(yōu)于純種荷斯坦牛；抗逆性良好，具有更好的體溫調(diào)節(jié)能力；同時還具有性情溫順、肢蹄健康的優(yōu)點(diǎn)。綜合分析，初步表明西荷雜種一代牛綜合性能優(yōu)于純種荷斯坦牛，充分利用遺傳互補(bǔ)和雜種優(yōu)勢的奶牛雜交繁育體系具有廣闊的發(fā)展前景。

致謝：感謝羅漢鵬、李源韜、蘇義童、黃銘逸等在試驗(yàn)設(shè)計環(huán)節(jié)提供的幫助，同時也感謝山西省奶牛生產(chǎn)性能測定中心楊繼業(yè)主任對數(shù)據(jù)采集的大力支持。

[1] KNOB D A, ALESSIO D R M, NETO A T, et al. Reproductive performance and survival of Holstein and Holstein×Simmental crossbred cows[J].TropAnimHealthProd, 2016, 48(7): 1409-1413.

[2] 王 婷, 初 芹, 劉 林, 等. 雜交繁育在奶牛生產(chǎn)中的應(yīng)用與探討[J]. 中國畜牧雜志, 2010, 46(3): 53-57. WANG T, CHU Q, LIU L, et al. Research and application of crossbreeding to dairy cattle[J].ChineseJournalofAnimalScience, 2010, 46(3): 53-57. (in Chinese)

[3] MCDOWELL R E, MCDANIEL B T. Interbreed matings in dairy cattle. I. Yield traits, feed efficiency, type and rate of milking[J].JDairySci, 1968, 51(5): 767-777.

[4] AHLBORN-BREIER G, HOHENBOKEN W D. Additive and nonadditive genetic effects on milk production in dairy cattle: evidence for major individual heterosis[J].JDairySci, 1991, 74(2): 592-602.

[5] 王志耕, 梅 林, 錢東方, 等. 西門塔爾雜交黃牛奶品質(zhì)研究[J]. 黃牛雜志, 2002, 28(5): 3-6. WANG Z G, MEI L, QIAN D F, et al. Milk quality of crossbred progeny of simmental with indigenous yellow cattle[J].JournalofYellowCattleScience, 2002, 28(5): 3-6. (in Chinese)

[6] HEINS B J, HANSEN L B, SEYORA A J. Crossbreds of Normande-Holstein, Montbeliarde-Holstein, and Scandinavian Red-Holstein compared rate, days open, and survival[J].JAnimSci, 2004, 83(1): 92-93.

[7] 孫少華, 何淑靜, 孫志穎, 等. 雜交荷斯坦母牛與其純種在產(chǎn)奶、繁殖及抗病性狀上的綜合對比[J]. 中國獸醫(yī)學(xué)報, 2016, 36(8): 1416-1421. SUN S H, HE S J, SUN Z Y, et al. Comparative study on milk performance, fertility and mastitis resistance between Holstein crossing heifers and Holstein heifers[J].ChineseJournalofVeterinaryScience, 2016, 36(8): 1416-1421. (in Chinese)

[8] 羅漢鵬, 武俊達(dá), 邱文卿, 等. 奶牛乳頭長度與泌乳性能相關(guān)性的初步探究[J]. 畜牧與獸醫(yī), 2016, 48(6): 67-69. LUO H P, WU J D, QIU W Q, et al. Preliminary study on the correlationship between teat-length and milk performance in dairy cow[J].AnimalHusbandry&VeterinaryMedicine, 2016, 48(6): 67-69. (in Chinese)

[9] 呂文龍. DHI報告的正確解讀[J]. 今日畜牧獸醫(yī), 2014(8): 53-55. LV W L. The proper unscramble of DHI report[J].TodayAnimalHusbandryandVeterinaryMedicine, 2014(8): 53-55. (in Chinese)

[10] 李雪嬌, 曹 瓊, 王消消, 等. 應(yīng)用體尺指標(biāo)估測西雜后備母牛體重[J]. 中國奶牛, 2012(21): 19-23. LI X J, CAO Q, WANG X X, et al. Predicting body weight of simmental crossbred replacement heifers using body measurement[J].ChinaDairyCattle, 2012(21): 19-23. (in Chinese)

[11] WILDMAN E E, JONES G M, WAGNER P E, et al. A dairy cow body condition scoring system and its relationship to selected production characteristics[J].JDairySci, 1982, 65(3): 495-501.

[12] EDMONSON A J, LEAN I J, WEAVER L D, et al. A body condition scoring chart for Holstein dairy cows[J].JDairySci, 1989, 72(1): 68-78.

[13] 董曉霞, 劉浩淼, 張 超, 等. 北京市氣候變化對奶牛熱冷應(yīng)激的影響[J]. 農(nóng)業(yè)工程學(xué)報, 2013, 29(16): 198-205. DONG X X, LIU H M, ZHANG C, et al. Impact of climate change on heat and cold stress of cow breeding in Beijing[J].TransactionsoftheChineseSocietyofAgriculturalEngineering, 2013, 29(16): 198-205. (in Chinese)

[14] 趙 瑋. 關(guān)注奶牛肢蹄健康, 重視步態(tài)評分管理[J]. 黑龍江畜牧獸醫(yī), 2013(12): 84-85. ZHAO W. Pay attention to the health of limbs and feet and the management of gait score[J].HeilongjiangAnimalScienceandVeterinaryMedicine, 2013(12): 84-85. (in Chinese)

[15] 郭翔羽, 郭 剛, 李錫智, 等. 北京地區(qū)荷斯坦青年牛繁殖性狀影響因素分析[J]. 中國奶牛, 2011(24): 1-6. GUO X Y, GUO G, LI X Z, et al. Analysis of causal factors on reproduction traits for Holstein heifers in Beijing area[J].ChinaDairyCattle, 2011(24): 1-6. (in Chinese)

[16] 中華人民共和國農(nóng)業(yè)部. NY/T 2363-2013奶牛熱應(yīng)激評價技術(shù)規(guī)范[S]. 北京: 中國農(nóng)業(yè)出版社, 2013. The Ministry of Agriculture of the People's Republic of China. NY/T 2363-2013 Technical specification for heat of dairy cows[S]. Beijing: China Agriculture Press, 2013. (in Chinese)

[17] HEINS B J, HANSEN L B, SEYKORA A J. Production of pure Holsteins versus crossbreds of Holstein with Normande, Montbeliarde, and Scandinavian Red[J].JDairySci, 2006, 89(7): 2799-2804.

[18] HEINS B J, HANSEN L B. Short communication: Fertility, somatic cell score, and production of Normande×Holstein, Montbéliarde × Holstein, and Scandinavian Red × Holstein crossbreds versus pure Holsteins during their first 5 lactations[J].JDairySci, 2012, 95(2): 918-924.

[19] 米占鋒. 西雜牛乳用性能的開發(fā)利用[J]. 畜牧獸醫(yī)雜志, 2013, 32(3): 61-62. MI Z F. Development and utilization of milking performance of Simmental crossbred cows[J].JournalofAnimalScienceandVeterinaryMedicine, 2013, 32(3): 61-62. (in Chinese)

[20] HAZEL A R, HEINS B J, SEYKORA A J, et al. Montbéliarde-sired crossbreds compared with pure Holsteins for dry matter intake, production, and body traits during the first 150 days of first lactation[J].JDairySci, 2013, 96(3): 1915-1923.

[21] HANSEN L B, HEINS B J, SEYKORA A J. Is crossbreeding the answer for reproductive problems of dairy cattle?[J/OL].http://www.genusbreeding.co.uk/doeuments/upload/Is_Crossbreeding_theAnswer_for_Reproducfive_Preblems.pdf.

[22] 何淑靜. 蒙貝利亞牛雜交荷斯坦母牛效果的綜合比較研究[D]. 保定: 河北農(nóng)業(yè)大學(xué), 2014. HE S J. Comparative study on effects of Montbeliarde-sired crossbred Holstein cows[D]. Baoding: Hebei Agricultural University, 2014. (in Chinese)

[23] 王敬春, 李樹春, 李洪宇, 等. 澳牛生產(chǎn)性能與適應(yīng)性的觀察報告[J]. 中國奶牛, 2006(6): 25-27. WANG J C, LI S C, LI H Y, et al. Report on production performance and adaptability of Australia dairy cows in China[J].ChinaDairyCattle, 2006(6): 25-27. (in Chinese)

[24] 周振峰, 王 晶, 王加啟, 等. 裹包TMR飼喂對泌乳中期奶牛生產(chǎn)性能、養(yǎng)分表觀消化率及血液生化指標(biāo)的影響[J]. 草業(yè)學(xué)報, 2010, 19(5): 31-37. ZHOU Z F, WANG J, WANG J Q, et al. Effects of feeding baled TMR on lactation performance, apparent digestibility, and blood metabolites in lactating dairy cows[J].ActaPrataculturaeSinica, 2010, 19(5): 31-37. (in Chinese)

[25] 王林楓, 嚴(yán) 平. 德系西門塔爾改良牛的生產(chǎn)性能和生化指標(biāo)分析[C]//第十屆(2015)中國牛業(yè)發(fā)展大會論文匯編. 固原: 中國畜牧業(yè)協(xié)會, 2015. WANG L F, YAN P. Analysis of production performance and biochemical indexes of German Simmental crossbreed cattle[C]//Guyuan: China Animal Agriculture Association, 2015. (in Chinese)

[26] MALCHIOD F, PENASA M, TIEZZI F, et al. Milk yield traits, somatic cell score, milking time and age at calving of pure Holstein versus crossbred cow[J].AgricConspecSci, 2011, 76(3): 259-261.

[27] GODDEN S M, LISSEMORE K D, KELTON D F, et al. Factors associated with milk urea concentrations in Ontario dairy cows[J].JDairySci, 2001, 84(1): 107-114.

[28] 毛永江, 張亞琴, 張美榮, 等. 澳系進(jìn)口荷斯坦牛初配年齡對部分繁殖性能的影響[C]//中國奶業(yè)協(xié)會第26次繁殖學(xué)術(shù)年會暨國家肉牛牦牛/奶牛產(chǎn)業(yè)技術(shù)體系第3屆全國牛病防治學(xué)術(shù)研討會論文集. 蘭州: 中國奶業(yè)協(xié)會繁殖專業(yè)委員會, 國家肉牛產(chǎn)業(yè)技術(shù)體系疾病控制功能研究室, 國家奶牛產(chǎn)業(yè)技術(shù)體系疾病控制功能研究室, 華中農(nóng)業(yè)大學(xué), 家畜疫病病原微生物學(xué)國家重點(diǎn)實(shí)驗(yàn)室, 2011. MAO Y J, ZHANG Y Q, ZHANG M R, et al. Effects of the age of first service on some reproductive performance of Australian imported Holstein cows[C]//Lanzhou: Dairy Association of China, 2011. (in Chinese)

[29] 馬洪祥, 楊俊奎, 王清河, 等. 西門塔爾牛繁殖規(guī)律探討[J]. 中國畜牧雜志, 1997, 33(4): 18-19. MA H X, YANG J K, WANG Q H, et al. Study on reproduction of Simmental cattle[J].ChineseJournalofAnimalScience, 1997, 33(4): 18-19. (in Chinese)

[30] MENDON?A L G D, ABADE C C, DA SILVA E M, et al. Comparison of peripartum metabolic status and postpartum health of Holstein and Montbéliarde-sired crossbred dairy cows[J].JDairySci, 2014, 97(2): 805-818.

[31] MALCHIOD F, CECCHINATO A, BITTANTE G. Fertility traits of purebred Holsteins and 2- and 3-breed crossbred heifers and cows obtained from Swedish Red, Montbéliarde, and Brown Swiss sires[J].JDairySci, 2014, 97(12): 7916-7926.

[32] 邢震全, 郭光成. 德國弗萊維赫牛對荷斯坦奶牛的改良應(yīng)用[J]. 黑龍江動物繁殖, 2010, 18(2): 46-47. XING Z Q, GUO G C. The application of the German Fleckvieh cattle for improvement in Holstein cows[J].HeilongjiangJournalofAnimalReproduction, 2010, 18(2): 46-47. (in Chinese)

[33] 傅春泉, 姜俊芳, 周文仙, 等. 德系西門塔爾牛在散欄式牛場的雜交利用研究[J]. 中國奶牛, 2016(3): 21-24. FU C Q, JIANG J F, ZHOU W X, et al. Crossbreeding study of German Simmental cattle in free stall feeding farm[J].ChinaDairyCattle, 2016(3): 21-24. (in Chinese)

[34] HAZEL A R, HEINS B J, SEYKORA A J, et al. Production, fertility, survival, and body measurements of Montbéliarde-sired crossbreds compared with pure Holsteins during their first 5 lactations[J].JDairySci, 2014, 97(4): 2512-2525.

[35] 陸東林. 奶牛體況評分及其應(yīng)用[J]. 新疆畜牧業(yè), 2006(5): 19-21. LU D L. Body condition score of dairy cows and its application[J].AnimalHusbandryinXinjiang, 2006(5): 19-21. (in Chinese)

[36] 李朝明. 荷斯坦奶牛、娟姍牛和娟荷雜交牛熱適應(yīng)性的研究[D]. 重慶: 西南大學(xué), 2013. LI Z M. Study on heat tolerance of Jersey, Holstein and Jersey × Holstein dairy cows[D]. Chongqing: Southwest University, 2013. (in Chinese)

[37] MONFORE G R. Considerations for dairying in the tropics[J].JAnimSci, 1991, 69(6): 2463-2471.

[38] 史彬林, 閆素梅. 奶牛耐熱性評定指標(biāo)的研究[J]. 中國奶牛, 1996(2): 20-22. SHI B L, YAN S M. Study on the indexes of heat tolerance in dairy cows[J].ChinaDairyCattle, 1996(2): 20-22. (in Chinese)

[39] 傅春泉, 姜俊芳, 何英俊. 弗萊維赫牛改良荷斯坦牛抗逆性研究[J]. 畜牧與獸醫(yī), 2015, 47(1): 37-39. FU C Q, JIANG J F, HE Y J. Study on stress resistance of Holstein cattle improved by Fleckvieh cattle[J].AnimalHusbandry&VeterinaryMedicine, 2015, 47(1): 37-39. (in Chinese)

[40] 張 林, 劉 青, 李燁青, 等. 西荷雜交犢牛生長性能及耐熱性的研究[J]. 畜牧與獸醫(yī), 2017, 49(3): 1-4. ZHANG L, LIU Q, LI Y Q, et al. Growth performance and heat tolerance in dairy cows of Simmental×Holstein F1[J].AnimalHusbandry&VeterinaryMedicine, 2017, 49(3): 1-4. (in Chinese)

[41] 禹旺盛, 張佐忠, 考桂蘭, 等. 澳系荷斯坦奶牛在錫盟南部地區(qū)的適應(yīng)性觀察[J]. 中國奶牛, 2005(1): 40-42. YU W S, ZHANG Z Z, KAO G L, et al. Observation of adaptability of Australian Holstein cows in the southern region of Ximeng[J].ChinaDairyCattle, 2005(1): 40-42. (in Chinese)

[42] TUCKER C B, ROGERS A R, VERKERK G A, et al. Effects of shelter and body condition on the behaviour and physiology of dairy cattle in winter[J].ApplAnimBehavSci, 2007, 105(1-3): 1-13.

[43] 李英超, 侯秀賢, 和利民, 等. 德系西門塔爾牛與奶牛雜交效果觀察[J]. 河南畜牧獸醫(yī), 2015, 36(11): 3-4. LI Y C, HOU X X, HE L M, et al. Observation of efficiency of crossbred cows between German Simmental and Holstein[J].HenanJournalofAnimalHusbandryandVeterinaryMedicine, 2015, 36(11): 3-4. (in Chinese)

[44] 張海亮, 徐 偉, 李 想, 等. 北京地區(qū)荷斯坦牛皮膚褶皺厚度及其影響因素分析[J]. 中國畜牧雜志, 2016, 52(20): 57-62. ZHANG H L, XU W, LI X, et al. Measurement and influencing factors analysis of skinfold thickness in Holstein in Beijing[J].ChineseJournalofAnimalScience, 2016, 52(20): 57-62.

[45] DOWLING D F. The thickness of cattle skin[J].AustJAgricRes, 1955, 6(5): 776-785.

[46] 成 曉, 羅漢鵬, 阮含章, 等. 體況評分與步態(tài)評分對奶牛生產(chǎn)性能的影響[J]. 黑龍江畜牧獸醫(yī), 2015(21): 1-5. CHENG X, LUO H P, RUAN H Z, et al. Effects of body condition score and gait score on the productive performance of dairy cows[J].HeilongjiangAnimalScienceandVeterinaryMedicine, 2015(21): 1-5. (in Chinese)

[47] 李 想，鄢新義，曹 露，等. 新舊版本《中國荷斯坦牛體型鑒定規(guī)程》對比分析[J]. 畜牧獸醫(yī)學(xué)報, 2016, 47(12): 2414-2419. LI X, YAN X Y, CAO L, et al. Comparative analysis of two versions of code of practice of type classification in Chinese Holstein[J].ActaVeterinariaetZootechnicaSinica, 2016, 47(12): 2414-2419.(in Chinese)