王 潔，梁 艷，張 強(qiáng)，王海洋，郭夢(mèng)玲，楊章平，毛永江*

(1. 揚(yáng)州大學(xué) 動(dòng)物科學(xué)與技術(shù)學(xué)院，江蘇 揚(yáng)州225009；2. 教育部農(nóng)業(yè)與農(nóng)產(chǎn)品安全國(guó)際合作聯(lián)合實(shí)驗(yàn)室，江蘇 揚(yáng)州225009)

乳蛋白是由乳腺上皮細(xì)胞合成并分泌的營(yíng)養(yǎng)物質(zhì)，在牛奶與乳制品中均起著重要作用。泌乳期指的是奶牛分娩后從開(kāi)始泌乳到泌乳結(jié)束這一階段。泌乳曲線是用來(lái)描述泌乳期內(nèi)奶牛群體或個(gè)體泌乳性狀隨時(shí)間變化而變化的一個(gè)數(shù)學(xué)模型，在奶牛生產(chǎn)和育種方面的有較強(qiáng)的指導(dǎo)作用，而乳蛋白率是奶牛泌乳性狀的重要指標(biāo)之一，也是衡量牛奶品質(zhì)的重要指標(biāo)。相關(guān)研究表明，乳蛋白率在泌乳初期逐漸下降，一定時(shí)間后達(dá)到最低點(diǎn)，后又隨時(shí)間推移而逐漸回升，在不同泌乳階段呈規(guī)律性變化，許多學(xué)者通過(guò)建立或借鑒各種數(shù)學(xué)模型來(lái)分析、描述、量化這一變化規(guī)律，發(fā)現(xiàn)Wood不完全伽瑪函數(shù)模型(Wood模型)、逆多項(xiàng)式模型(Nelder模型)、Wilmink模型、Ali-Schaeffer模型和Dijkstra模型等對(duì)泌乳曲線有較好的擬合效果。孫先枝等研究發(fā)現(xiàn)，Wood模型適用于長(zhǎng)三角地區(qū)荷斯坦牛乳蛋白率變化曲線的擬合。毛永江等研究發(fā)現(xiàn)，Wood模型對(duì)南方中國(guó)荷斯坦牛乳蛋白率變化曲線擬合度均在0.99以上。

相關(guān)研究表明，奶牛泌乳曲線模型參數(shù)受場(chǎng)效應(yīng)、胎次、產(chǎn)犢季節(jié)等多種非遺傳因素影響。近年來(lái)，隨著牧場(chǎng)硬件條件的改善和飼養(yǎng)管理水平的提升，奶牛生產(chǎn)性能得到大幅度提高，從而對(duì)奶牛乳蛋白率泌乳曲線產(chǎn)生一定影響。因此，本研究擬用Wood模型對(duì)江蘇地區(qū)12個(gè)牧場(chǎng)2017-2019年荷斯坦牛乳蛋白率的變化進(jìn)行曲線擬合，定量的描述其變化趨勢(shì)，并分析不同牧場(chǎng)規(guī)模、胎次、產(chǎn)犢季節(jié)、產(chǎn)犢間隔和305 d產(chǎn)奶量等非遺傳因素對(duì)荷斯坦牛乳蛋白率Wood泌乳曲線參數(shù)的影響，從而為江蘇地區(qū)荷斯坦?？茖W(xué)飼養(yǎng)等方面提供參考依據(jù)。

1 材料與方法

1.1 數(shù)據(jù)來(lái)源

供試數(shù)據(jù)來(lái)自江蘇地區(qū)12個(gè)奶牛場(chǎng)2017～2019年荷斯坦牛DHI測(cè)定日記錄共580 025條，主要包括：牛號(hào)、測(cè)定日期、胎次、泌乳天數(shù)、乳蛋白率、305 d產(chǎn)奶量等。為保證最終結(jié)果的有效性，數(shù)據(jù)質(zhì)控條件如下：胎次為1～5胎，泌乳天數(shù)≤305 d，乳蛋白率為1%～7%。最后合格記錄數(shù)467 396條。根據(jù)江蘇省氣候特點(diǎn)，將產(chǎn)犢季節(jié)劃分為春季(3-5月)、夏季(6-8月)、秋季(9-11月)、冬季(12月到次年2月)；每30 d劃分為1個(gè)泌乳月，最后一個(gè)泌乳月為271～305 d，共10個(gè)泌乳月。表1為數(shù)據(jù)的基本情況。

表1 泌乳數(shù)據(jù)的描述性統(tǒng)計(jì)Table 1 Descriptive statistils of lactating data

1.2 泌乳曲線擬合

用Wood不完全伽瑪函數(shù)模型對(duì)測(cè)定日乳蛋白率變化進(jìn)行擬合，基本模型如下：

Y=ate

式中：t為泌乳月；Y為時(shí)間t的乳蛋白率；a、b、c是模型參數(shù)。參數(shù)a是泌乳潛力；c是達(dá)到乳蛋白率最低點(diǎn)的速度；b是從乳蛋白率最低點(diǎn)上升的速度。a、b、c均由平均值±標(biāo)準(zhǔn)誤形式表示。在進(jìn)行曲線擬合時(shí)，各參數(shù)初始值的確定采用Olori等計(jì)算所得的結(jié)果。再使用如下公式換算出以下三個(gè)二級(jí)參數(shù)：

T= b/c

Y= a(b/c)e

Per= -(b+1)lnc

式中：T為達(dá)到乳蛋白率最低點(diǎn)所需的時(shí)間；Y為乳蛋白率的最低值；Per為乳蛋白率的持久力。先用Excel整合乳蛋白率數(shù)據(jù)，然后用SPSS(Ver 26.0)非線性回歸子程序(Nonlinear Regression)進(jìn)行擬合。擬合效果用擬合度R和誤差均方進(jìn)行評(píng)價(jià)。

采用最小二乘模型分析非遺傳因素對(duì)奶牛乳蛋白率的影響，模型如下：

Y=μ+F+P+S+D+N+e

式中：Y為乳蛋白率的觀察值；F為牛場(chǎng)規(guī)模的固定效應(yīng)；P為胎次的固定效應(yīng)；S為產(chǎn)犢季節(jié)的固定效應(yīng)；D為產(chǎn)犢間隔的固定效應(yīng)；N為305 d產(chǎn)奶量的固定效應(yīng)；e為隨機(jī)殘差。各因素不同水平間的多重比較使用Duncan's法，顯著性水平定義如下：≤0.01為極顯著，≤0.05為顯著，>0.05為不顯著。根據(jù)上述影響測(cè)定日乳蛋白率的因素和水平，分別選取相應(yīng)數(shù)據(jù)，利用如前所述方法分別進(jìn)行Wood泌乳曲線擬合，并計(jì)算相應(yīng)參數(shù)。

2 結(jié)果與分析

2.1 不同因素對(duì)荷斯坦牛乳蛋白率的影響

由表2可知，牧場(chǎng)、胎次、產(chǎn)犢季節(jié)、產(chǎn)犢間隔和305 d產(chǎn)奶量對(duì)荷斯坦牛乳蛋白率有極顯著影響(<0.01)。其中，規(guī)模在5 000頭以上的牧場(chǎng)荷斯坦牛乳蛋白率最高(3.31%)，規(guī)模在2 001～5 000頭的牧場(chǎng)荷斯坦牛乳蛋白率最低(3.23%)；第2胎荷斯坦牛乳蛋白率為3.32%，極顯著高于3～5胎(<0.01)，第5胎乳蛋白率為3.25%，極顯著低于其它胎次(<0.01)；冬季產(chǎn)犢的荷斯坦牛乳蛋白率為3.31%，極顯著高于其他產(chǎn)犢季節(jié)(<0.01)，春季產(chǎn)犢乳蛋白率為3.28%，極顯著低于其它產(chǎn)犢季節(jié)(<0.01)；荷斯坦牛乳蛋白率隨著產(chǎn)犢間隔的增加而極顯著下降，分別為3.32%(300～365 d)、3.30%(366～400 d)、3.29%(401～420 d)、3.28%(421～440 d)和3.27%(>441 d)；305 d產(chǎn)奶量為3 000～5 000 kg和5 001～7 000 kg的荷斯坦牛乳蛋白率極顯著高于其他305 d產(chǎn)奶量范圍奶牛的乳蛋白率(<0.01)。

表2 不同因素對(duì)乳蛋白率的影響Table 2 Effects of different factors on milk protein percentage

2.2 不同因素對(duì)日乳蛋白率曲線擬合的影響

由圖1和表3可知，規(guī)模在5 000頭以上的奶牛乳蛋白率曲線擬合度最大(0.9932)且誤差均方最小(0.010)，規(guī)模在2 001～5 000頭的奶牛乳蛋白率曲線擬合度最小(0.9921)；規(guī)模在5 000頭以上的奶牛乳蛋白率最低值最大(3.11%)，規(guī)模在1 000頭以下的奶牛乳蛋白率最低值最小(3.02%)；規(guī)模在1 000頭以下的奶牛達(dá)到乳蛋白率最低值的時(shí)間最大(87.43 d)，規(guī)模在5 000頭以上的奶牛達(dá)到乳蛋白率最低值的時(shí)間最小(77.00 d)；規(guī)模在1 000～2 000頭的奶牛乳蛋白率持久力最大(3.81)，規(guī)模在1 000頭以下的奶牛乳蛋白率持久力最小(3.69)；規(guī)模在1 000～2 000頭的奶牛乳蛋白率潛力最大(3.106)，規(guī)模在1 000頭以下的奶牛乳蛋白率潛力最小(3.038)；規(guī)模在1 000頭以下的奶牛達(dá)到乳蛋白率最低點(diǎn)的速度和最低點(diǎn)上升的速度最大，分別為0.035和0.102。

圖1 不同牧場(chǎng)規(guī)模下Wood模型擬合泌乳曲線Fig.1 Lactation curves of Wood model for different farm sizes

結(jié)合圖2和表3可知，荷斯坦牛第1胎的乳蛋白率曲線擬合度最大(0.9937)且誤差均方最小(0.009)，第5胎擬合度最小(0.9918)；荷斯坦牛第2胎乳蛋白率最低值最大(3.10%)，第5胎乳蛋白率最低值最小(3.06%)；荷斯坦牛第4胎達(dá)到乳蛋白率最低值的時(shí)間最大(90.00 d)，第1胎達(dá)到乳蛋白率最低值的時(shí)間最小(63.75 d)；荷斯坦牛第1胎乳蛋白率持久力最大(3.92)，第3胎乳蛋白率持久力最小(3.69)；荷斯坦牛第3胎乳蛋白率潛力最大(3.108)，第1胎乳蛋白率潛力最小(3.050)；荷斯坦牛第3、4胎達(dá)到乳蛋白率最低點(diǎn)的速度最大(0.035)，第4胎由最低點(diǎn)上升的速度最大(0.105)。

圖2 不同胎次Wood模型擬合泌乳曲線Fig.2 Lactation curves of Wood model for different parities

結(jié)合圖3和表3可知，荷斯坦牛冬季產(chǎn)犢乳蛋白率曲線擬合度最大(0.9933)，夏季產(chǎn)犢乳蛋白率曲線擬合度最小(0.9925)；冬季產(chǎn)犢乳蛋白率最低值最大(3.14%)，春季產(chǎn)犢乳蛋白率最低值最小(2.99%)；冬季產(chǎn)犢達(dá)到乳蛋白率最低值的時(shí)間最晚(99.33 d)，夏季產(chǎn)犢達(dá)到乳蛋白率最低值的時(shí)間最早(28.57 d)；秋季產(chǎn)犢乳蛋白率持久力最大(4.08)，春季產(chǎn)犢乳蛋白率持久力最小(3.43)；春季產(chǎn)犢達(dá)到乳蛋白率最低點(diǎn)的速度最大(0.047)，冬季產(chǎn)犢乳蛋白率潛力和乳蛋白率由最低點(diǎn)上升的速度最大，分別為3.159和0.149。

圖3 不同產(chǎn)犢季節(jié)Wood模型擬合泌乳曲線Fig.3 Lactation curves of Wood model for different calving seasons

由圖4和表3可知，荷斯坦牛產(chǎn)犢間隔為401～420 d時(shí)乳蛋白率曲線擬合度最大(0.9928)且誤差均方最小(0.010)，產(chǎn)犢間隔441 d以上時(shí)乳蛋白率曲線擬合度最小(0.9919)；產(chǎn)犢間隔為300～365 d時(shí)乳蛋白率最低值最大(3.13%)，產(chǎn)犢間隔為441 d以上時(shí)乳蛋白率最低值最小(3.02%)；產(chǎn)犢間

圖4 不同產(chǎn)犢間隔Wood模型擬合泌乳曲線Fig.4 Lactation curves of Wood model for different calving interval

表3 不同因素日乳蛋白率Wood模型擬合參數(shù)Table 3 Wood model fitting parameters of daily milk protein ratio by different factors

隔為421～440 d時(shí)達(dá)到乳蛋白率最低值的時(shí)間最晚(92.00 d)，產(chǎn)犢間隔441 d以上時(shí)達(dá)到乳蛋白率最低值的時(shí)間最早(75.38 d)；產(chǎn)犢間隔為300～365 d時(shí)乳蛋白率持久力最大(3.77)，產(chǎn)犢間隔為421～440 d時(shí)乳蛋白率持久力最小(3.53)；產(chǎn)犢間隔為300～365 d時(shí)乳蛋白率潛力最大(3.156)，441 d以上時(shí)乳蛋白率潛力最?。划a(chǎn)犢間隔為421～440 d時(shí)達(dá)到乳蛋白率最低點(diǎn)的速度和由最低點(diǎn)上升的速度最大，分別為0.045和0.138。

結(jié)合圖5和表3可知，305 d產(chǎn)奶量為11 001～13 000 kg的荷斯坦牛乳蛋白率曲線擬合度最大(0.9944)且誤差均方最小(0.008)，3 000～5 000 kg的荷斯坦牛乳蛋白率曲線擬合度最小(0.9919)；305 d產(chǎn)奶量為7 001～9 000 kg的荷斯坦牛乳蛋白率最低值最大(4.17%)，13 001～15 000 kg的荷斯坦牛乳蛋白率最低值最小(3.04%)；305 d產(chǎn)奶量為7 001～9 000 kg的荷斯坦牛達(dá)到乳蛋白率最低值的時(shí)間最大(3420.00 d)，13 001～15 000 kg的荷斯坦牛達(dá)到乳蛋白率最低值的時(shí)間最小(80.00 d)；305 d產(chǎn)奶量為9 001～11 000 kg的荷斯坦牛乳蛋白率持久力最大(4.55)，13 001～15 000 kg的荷斯坦牛乳蛋白率持久力最小(3.87)；305 d產(chǎn)奶量為13 001～15 000 kg的荷斯坦牛乳蛋白率潛力、達(dá)到乳蛋白率最低點(diǎn)的速度和由最低點(diǎn)上升的速度最大，分別為3.038、0.027和0.072。

圖5 不同305 d產(chǎn)奶量Wood模型擬合泌乳曲線Fig.5 Lactation curves of Wood model for different 305-day milk production

3 討 論

3.1 牧場(chǎng)規(guī)模對(duì)Wood泌乳曲線模型的影響

王瑞軍等和張佳蘭等研究表明，場(chǎng)效應(yīng)對(duì)荷斯坦牛泌乳曲線模型參數(shù)有顯著影響，不同牧場(chǎng)泌乳曲線模型參數(shù)的差異顯著。本研究發(fā)現(xiàn)，不同規(guī)模的牧場(chǎng)對(duì)荷斯坦牛乳蛋白率Wood泌乳曲線模型參數(shù)有一定影響，其中規(guī)模在5000頭以上的牧場(chǎng)荷斯坦牛乳蛋白率曲線擬合度和乳蛋白率最低值最大，達(dá)到乳蛋白率最低值的時(shí)間最快。相關(guān)研究表明，規(guī)模在5 000頭以上的牧場(chǎng)屬于大規(guī)模牧場(chǎng)，場(chǎng)中設(shè)備機(jī)械化程度高，操作技術(shù)力量強(qiáng)，并且對(duì)不同泌乳期和產(chǎn)奶量水平的奶牛使用不同飼料配方的全混合日糧(Total mixed rations, TMR)進(jìn)行分群飼養(yǎng)，進(jìn)而提高了奶牛泌乳性能和泌乳效率，因此平均乳蛋白率和乳蛋白潛力較高。此外，大規(guī)模牛場(chǎng)較之中小型牛場(chǎng)具有更完善的飼養(yǎng)管理措施，并經(jīng)過(guò)規(guī)范且科學(xué)的選育提高，因此泌乳效率高于小型牧場(chǎng)。

3.2 胎次對(duì)Wood泌乳曲線模型的影響

本研究發(fā)現(xiàn)，第5胎奶牛乳蛋白率曲線擬合度最低，乳蛋白率最低值最?。坏?胎牛乳蛋白率持久力最大，潛力最小；第3胎牛乳蛋白率潛力最大，持久力最小。Roberto等通過(guò)其研究發(fā)現(xiàn)，Wood模型擬合群體泌乳曲線參數(shù)極顯著受到胎次影響，頭胎牛泌乳潛力較低。孫先枝等研究發(fā)現(xiàn)，1胎奶牛乳蛋白率持久力最高，3胎牛乳蛋白率潛力最大，持久力最小。上述研究與本研究基本一致。對(duì)于頭胎牛，奶牛機(jī)體發(fā)育尚不完善，乳蛋白率潛力較低，泌乳后期乳蛋白率上升較慢，乳蛋白率持久力強(qiáng)。奶牛泌乳性能、乳蛋白潛力隨著胎次的增加逐漸提升，第3胎時(shí)乳蛋白率潛力達(dá)到最大。就乳蛋白率而言，持久力越大，說(shuō)明其維持較低值的能力越強(qiáng)，乳蛋白率越低，第3胎乳蛋白率持久力最小，泌乳性能最佳，乳蛋白率潛力最大。對(duì)于荷斯坦牛第5胎乳蛋白率曲線擬合度最低和乳蛋白率最低值最小，Knaus等和Oltenacu等研究表明奧地利和美國(guó)奶牛的平均淘汰胎次為3.3胎，鄢新義等研究發(fā)現(xiàn)中國(guó)荷斯坦牛平均淘汰胎次為2.86胎；高胎次奶牛機(jī)體免疫力降低，泌乳機(jī)能下降，但也不排除部分5胎以上奶牛生產(chǎn)性能較好，因此造成5胎奶牛乳蛋白率變化差異較大，從而出現(xiàn)泌乳曲線擬合度最低和乳蛋白率最低值最小。但這一結(jié)果還有待擴(kuò)大樣本量，進(jìn)一步深入研究。

3.3 產(chǎn)犢季節(jié)對(duì)Wood泌乳曲線模型的影響

亓建剛等研究發(fā)現(xiàn)，產(chǎn)犢季節(jié)對(duì)Wood泌乳曲線模型參數(shù)均有一定影響，其中冬季產(chǎn)犢擬合度最高。本研究中，不同產(chǎn)犢季節(jié)對(duì)荷斯坦牛乳蛋白率和Wood泌乳曲線模型參數(shù)均有相應(yīng)程度的影響，其中冬季產(chǎn)犢乳蛋白率變化曲線擬合度、乳蛋白率最低值和乳蛋白率潛力最大；夏季產(chǎn)犢擬合度最小，達(dá)到乳蛋白率最低值的時(shí)間也最快。相關(guān)研究表明，奶牛在夏季產(chǎn)犢機(jī)體免疫力下降，同時(shí)最易受外界環(huán)境影響，在熱應(yīng)激情況下奶牛的各項(xiàng)生理機(jī)能出現(xiàn)紊亂，為補(bǔ)足能量和糖的需求，糖原分解和糖異生增加，這使得蛋白質(zhì)代謝增強(qiáng)，曲線擬合度下降。同時(shí)，奶牛夏季產(chǎn)犢受熱應(yīng)激影響，降低了奶牛干物質(zhì)采食量和脂肪組織動(dòng)員能力，并通過(guò)激素信號(hào)通路促進(jìn)了葡萄糖的利用，從而使維持需求增加，用于泌乳的能量減少，乳蛋白率快速降低，達(dá)到乳蛋白率最低值的時(shí)間也最小。因此，從維持奶牛身體健康和提高泌乳性能方面而言，牧場(chǎng)應(yīng)做好有針對(duì)性的防暑降溫工作，同時(shí)在可能的范圍內(nèi)合理安排牛群配種時(shí)間，盡量避免奶牛在酷暑季節(jié)產(chǎn)犢，以獲得更高的養(yǎng)殖效益。

3.4 產(chǎn)犢間隔對(duì)Wood泌乳曲線模型的影響

本研究發(fā)現(xiàn)，奶牛產(chǎn)犢間隔對(duì)乳蛋白率和Wood泌乳曲線模型參數(shù)有一定影響，其中產(chǎn)犢間隔為401～420 d時(shí)泌乳曲線擬合度最大，產(chǎn)犢間隔為441 d以上時(shí)乳蛋白率變化曲線擬合度、乳蛋白率潛力、乳蛋白率最低值和達(dá)到乳蛋白率最低值的時(shí)間最小。研究表明，奶牛下丘腦-垂體-卵巢軸受到生長(zhǎng)激素(GH)、胰島素樣生長(zhǎng)因子-I(IGF-I)、胰島素、甲狀腺激素和瘦蛋白等多種代謝物和激素調(diào)節(jié)。對(duì)于產(chǎn)犢間隔為441 d以上的荷斯坦牛，體內(nèi)分泌信號(hào)因子IGF-I和瘦蛋白處于負(fù)平衡狀態(tài)，機(jī)體的新陳代謝能力和生殖激素的分泌能力下降，繁殖疾病增多，從而使產(chǎn)妊間隔延長(zhǎng)，泌乳天數(shù)也隨之延長(zhǎng)，從而導(dǎo)致產(chǎn)犢間隔較長(zhǎng)，因此產(chǎn)犢間隔441 d以上的奶牛乳蛋白率變化曲線擬合度最低，乳蛋白率潛力和乳蛋白率最低值最小。

3.5 305 d產(chǎn)奶量對(duì)Wood泌乳曲線模型的影響

白音塔拉研究表明，不同305 d產(chǎn)奶量對(duì)乳蛋白率Wood泌乳曲線模型有一定影響，其中經(jīng)產(chǎn)牛平均305 d產(chǎn)奶量為11 328.8 kg時(shí)Wood泌乳曲線模型擬合效果最好。本研究發(fā)現(xiàn)，305 d產(chǎn)奶量為11 001～13 000 kg的荷斯坦牛乳蛋白率變化曲線擬合度最大，3 000～5 000 kg的荷斯坦牛乳蛋白率變化曲線擬合度最??；305 d產(chǎn)奶量為13 001～15 000 kg的荷斯坦牛泌乳潛力、達(dá)到乳蛋白率最低點(diǎn)的速度和由最低點(diǎn)上升的速度最大，且乳蛋白率持久力最小。此外，305 d產(chǎn)奶量為3 000～5 000 kg的荷斯坦牛部分因?yàn)樵绠a(chǎn)，部分可能為遺傳和疾病等原因，導(dǎo)致泌乳性能下降，乳蛋白率曲線擬合度降低。對(duì)于305 d產(chǎn)奶量為13 001～15 000 kg時(shí)Wood泌乳曲線模型參數(shù)變化，其原因可能是由于本研究中荷斯坦牛305 d產(chǎn)奶量多集中在7001～13 000 kg之間，而產(chǎn)奶量在13 001～15 000 kg的樣本量過(guò)少所致，因此還有待擴(kuò)大樣本量，以進(jìn)行更深入的研究。

4 結(jié) 論

Wood模型適用于江蘇地區(qū)荷斯坦牛測(cè)定日乳蛋白率變化曲線的擬合分析，且牧場(chǎng)規(guī)模、胎次、產(chǎn)犢季節(jié)、產(chǎn)犢間隔、305 d產(chǎn)奶量等非遺傳因素對(duì)乳蛋白率Wood泌乳曲線模型參數(shù)均有顯著影響。該結(jié)果為江蘇地區(qū)荷斯坦牛精準(zhǔn)養(yǎng)殖、選種選配等方面提供了參考依據(jù)。