云南水牛乳成分分析及乳能量預(yù)測(cè)模型的建立

李 清 毛華明* 李 文 李建華 魯瓊芬 何鴻源 劉 琴

(1.云南農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)技術(shù)學(xué)院，昆明650201；2.云南省動(dòng)物營(yíng)養(yǎng)與飼料重點(diǎn)實(shí)驗(yàn)室，昆明650201；3.湖南工業(yè)大學(xué)生命科學(xué)與化學(xué)學(xué)院，株洲412007；4.云南省安寧市畜牧獸醫(yī)站，昆明650300；5.云南省昆明市奶牛生產(chǎn)性能測(cè)定中心，昆明650041)

云南水牛乳成分分析及乳能量預(yù)測(cè)模型的建立

李 清1，2毛華明1，2*李 文1，3*李建華4魯瓊芬1何鴻源1劉 琴5

(1.云南農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)技術(shù)學(xué)院，昆明650201；2.云南省動(dòng)物營(yíng)養(yǎng)與飼料重點(diǎn)實(shí)驗(yàn)室，昆明650201；3.湖南工業(yè)大學(xué)生命科學(xué)與化學(xué)學(xué)院，株洲412007；4.云南省安寧市畜牧獸醫(yī)站，昆明650300；5.云南省昆明市奶牛生產(chǎn)性能測(cè)定中心，昆明650041)

本試驗(yàn)旨在分析水牛乳成分中乳脂肪、乳蛋白、乳糖和乳尿素氮含量及其與乳能量的相關(guān)性，建立產(chǎn)奶凈能的預(yù)測(cè)模型，用于奶水牛泌乳期科學(xué)的飼養(yǎng)管理。試驗(yàn)采取云南省主要奶水牛養(yǎng)殖場(chǎng)和養(yǎng)殖小區(qū)的2014—2015年9—12月份送達(dá)昆明市奶牛生產(chǎn)性能測(cè)定中心的304個(gè)原料乳樣本，水牛乳的乳成分由昆明市奶牛生產(chǎn)性能測(cè)定中心用MilkoScan FT+FC乳成分體細(xì)胞聯(lián)用儀和MilkoScan FT 120乳成分分析儀測(cè)定，鮮乳經(jīng)真空干燥箱恒溫烘干后用氧彈式熱量計(jì)測(cè)定乳能量。結(jié)果表明：乳脂肪、乳蛋白、乳糖、乳尿素氮、乳總固形物含量和乳能量分別為6.45%、4.55%、5.31%、13.60 mg/dL、18.77%和4.02 MJ/kg；乳能量分別與乳脂肪和乳蛋白含量呈現(xiàn)極顯著正相關(guān)(r=0.896 0、r=0.563 0，P<0.01)，乳能量受乳脂肪和乳蛋白含量的影響較大。分別以乳脂肪(F)，乳脂肪和乳蛋白(F和P)，乳脂肪、乳蛋白和乳糖(F、P和La)，乳脂肪、乳蛋白、乳糖和乳尿素氮含量(F、P、La和MUN)為預(yù)測(cè)因子，建立的預(yù)測(cè)乳能量(E)的一元、二元、三元及四元回歸方程的擬合度均在0.90以上，方程分別為：E=0.388F+1.540 (R2=0.933 6，P<0.01)；E=0.373F+0.221P+0.460 (R2=0.926 7，P<0.01)；E=0.396F+0.186P+0.105La-0.104 (R2=0.954 0，P<0.01)；E=0.397F+0.187P+0.106La+0.002MUN-0.146 (R2=0.958 0，P<0.01)。由此可見(jiàn)，可通過(guò)水牛乳中的乳脂肪、乳蛋白、乳糖及乳尿素氮含量預(yù)測(cè)水牛乳產(chǎn)奶凈能。

水牛乳；乳成分；乳能量；預(yù)測(cè)模型

能量在動(dòng)物生長(zhǎng)發(fā)育、繁殖、生產(chǎn)過(guò)程中是必不可少的，對(duì)于泌乳牛，乳的能量即產(chǎn)奶凈能(NEL)。目前，采用凈能體系的國(guó)家統(tǒng)一用產(chǎn)奶凈能、4%乳脂率和產(chǎn)奶量來(lái)計(jì)算泌乳牛的能量需求[1-6]。乳的能量來(lái)源于乳脂肪、乳蛋白、乳糖等乳成分，但是在實(shí)際生產(chǎn)中，測(cè)定牛乳中的能量相對(duì)復(fù)雜，而牛乳中的營(yíng)養(yǎng)成分則是每個(gè)養(yǎng)牛場(chǎng)每月必測(cè)的指標(biāo)，較容易獲得，乳脂肪和乳蛋白含量不僅可預(yù)測(cè)乳中能量值，還可預(yù)測(cè)產(chǎn)奶牛的能量是否平衡[7-8]。美國(guó)奶牛飼養(yǎng)標(biāo)準(zhǔn)、英國(guó)奶牛飼養(yǎng)標(biāo)準(zhǔn)及中國(guó)奶牛飼養(yǎng)標(biāo)準(zhǔn)中通過(guò)測(cè)定牛乳中的營(yíng)養(yǎng)成分含量，推算出產(chǎn)奶凈能的預(yù)測(cè)模型，從而指導(dǎo)生產(chǎn)[4-6]。在這些預(yù)測(cè)模型中所采集的牛奶樣本均為荷斯坦血系的奶牛，它們的乳脂肪、乳蛋白、乳糖及乳尿素氮含量均低于水牛乳中相應(yīng)含量，二者乳成分差異較大，尤其是水牛乳中尿素氮的含量較高，其在乳能量中是否有貢獻(xiàn)？對(duì)于奶水牛的產(chǎn)奶凈能能否套用荷斯坦奶牛的預(yù)測(cè)模型值得考量。

據(jù)聯(lián)合國(guó)糧食及農(nóng)業(yè)組織(FAO)2014年統(tǒng)計(jì)，奶水牛全球存欄量為19 509.83萬(wàn)頭，鮮水牛乳占總奶量的比重已上升到12.92%[9]，但水牛奶量單產(chǎn)低，除了加強(qiáng)品種選育提高單產(chǎn)外，合理的飼糧及科學(xué)的飼養(yǎng)管理則是提高單產(chǎn)的核心[10]。但目前，尚未有確定的奶水牛營(yíng)養(yǎng)需要量指導(dǎo)生產(chǎn)，因此獲得奶水牛的營(yíng)養(yǎng)需要量意義重大。本試驗(yàn)用MilkoScan FT+FC乳成分體細(xì)胞聯(lián)用儀和MilkoScan FT 120乳成分分析儀測(cè)定水牛乳成分，用氧彈式熱量計(jì)測(cè)定乳能量，探討水牛乳能量與乳成分的相關(guān)性及通過(guò)回歸分析建立乳能量的預(yù)測(cè)模型，得到奶水牛的泌乳能量需要，為奶水牛飼養(yǎng)標(biāo)準(zhǔn)的制定提供基礎(chǔ)數(shù)據(jù)。

1 材料與方法

1.1 試驗(yàn)材料

為了使所測(cè)定樣本具有代表性，試驗(yàn)選擇了云南省主要奶水牛養(yǎng)殖地不同的飼養(yǎng)模式及飼糧水平的3個(gè)奶水牛群體進(jìn)行分析。3個(gè)群體的試驗(yàn)動(dòng)物每天飼喂2次，自由飲水，每日擠奶2次。

檳榔江水牛乳來(lái)源于騰沖檳榔江水牛核心養(yǎng)殖場(chǎng)，產(chǎn)犢胎次主要為2～4胎、泌乳早中期健康的184頭檳榔江水牛，飼糧為自配精料、全株玉米青貯、稻草，飼喂量為每頭牛每天精料3～4 kg(根據(jù)產(chǎn)奶量)、玉米青貯20～25 kg、干稻草自由采食(攝入量4～5 kg)，蛋白質(zhì)水平為12%。大理雜交水牛乳來(lái)自于大理雜交水牛養(yǎng)殖小區(qū)，產(chǎn)犢胎次主要為1～5胎、泌乳早中期健康的56頭大理雜交水牛，飼糧為市售精料、玉米秸稈青貯、稻草、苜蓿干草、羊草，飼喂量為每頭牛每天精料3～4 kg(根據(jù)產(chǎn)奶量)、玉米秸稈青貯15～20 kg、干稻草自由采食(攝入量4～5 kg)，蛋白質(zhì)水平為14%；泌乳中期增加苜蓿干草和羊草共2 kg。德宏雜交水牛乳來(lái)源于德宏雜交水牛養(yǎng)殖小區(qū)，產(chǎn)犢胎次主要為1～5胎、泌乳早中期健康的64頭德宏雜交水牛，飼糧為自配精料、甘蔗梢、稻草、啤酒糟、王草，飼喂量為每頭牛每天精料1.5～2.0 kg(根據(jù)產(chǎn)奶量)、甘蔗梢20 kg、稻草和王草自由采食(攝入量4～5 kg)，蛋白質(zhì)水平為11%；泌乳中期增加啤酒糟2.0 kg。

1.2 樣品采集及測(cè)定

分別采集云南省騰沖、大理和德宏3個(gè)主要奶水牛養(yǎng)殖場(chǎng)和養(yǎng)殖小區(qū)的2014—2015年9—12月份的304份牛奶樣品進(jìn)行乳成分含量及乳能量的測(cè)定。采樣時(shí)間均為每月的月底，將每天2次的牛奶樣品按1∶1充分混勻后倒入采樣瓶搖勻，確保防腐劑溶于牛奶樣品，每頭牛采樣40～50 mL，分成2份，將樣品低溫保存迅速運(yùn)輸至實(shí)驗(yàn)室2～7 ℃可保存1周。

1.3 乳成分含量及乳能量的測(cè)定

乳成分含量的測(cè)定由昆明市奶牛生產(chǎn)性能測(cè)定中心用MilkoScan FT+FC乳成分體細(xì)胞聯(lián)用儀和MilkoScan FT 120乳成分分析儀測(cè)定完成。

乳能量測(cè)定在云南農(nóng)業(yè)大學(xué)用BH-ⅢS氧彈式熱量計(jì)測(cè)定：移取5 mL牛奶樣品放入已知重量的自封袋中稱重(自封袋置于古氏坩堝中)，放入真空干燥箱(60±5) ℃干燥后用氧彈式熱量計(jì)測(cè)定乳能量。

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

試驗(yàn)數(shù)據(jù)整理及作圖用Excel 2016和SPSS 22.0軟件進(jìn)行，差異顯著性檢驗(yàn)采用單因素方差分析(one-way ANOVA)和Duncan氏多重比較法，由模型Yi=μ+Ti+εi(Yi為觀測(cè)值，μ為總體平均值，Ti為處理效應(yīng)；εi為隨機(jī)誤差)獲得，結(jié)果用平均值±標(biāo)準(zhǔn)差表示。相關(guān)性分析采用二元變量相關(guān)分析中的Pearson相關(guān)系數(shù)判定；相關(guān)分析模型采用多元線性回歸分析得到，用模型E=β0+β1F+β2P+β3La+β4MUN[β0為常數(shù)；β1、β2、β3、β4為偏回歸系數(shù)；F為乳脂肪含量(%)；P為乳蛋白含量(%)；La為乳糖含量(%)；MUN為乳尿素氮含量(mg/dL)；E為乳能量(MJ/kg)]表示。

2 結(jié) 果

2.1 水牛乳成分比較

由表1可知，乳脂肪、乳蛋白、乳糖、乳尿素氮、乳總固形物含量和乳能量的平均值分別為6.45%、4.55%、5.31%、13.60 mg/dL、18.77%和4.02 MJ/kg。乳脂肪含量在3.30%～11.10%，且隨泌乳時(shí)間的變化而有所上升；乳蛋白含量在2.84%～6.58%；乳糖含量在3.09%～5.99%；乳尿素氮含量范圍較廣，在1.40～26.60 mg/dL；乳總固形物含量在11.42%～25.93%。

同一品種不同月份間水牛乳成分比較，騰沖檳榔江水牛乳的乳脂肪、乳蛋白、乳糖、乳尿素氮、乳總固形物含量和乳能量差異均不顯著(P>0.05)；大理雜交水牛乳和德宏雜交水牛乳中的乳蛋白、乳糖、乳總固形物含量和乳能量差異均不顯著(P>0.05)，但二者9月份的乳脂肪、乳尿素氮含量均與12月份差異顯著(P<0.05)。

表1 同一品種不同月份水牛乳成分比較

同一品種同列數(shù)據(jù)肩標(biāo)不同小寫字母表示差異顯著(P<0.05)，無(wú)字母或相同小寫字母表示差異不顯著(P>0.05)。

In the same breed and column, values with different small letter superscripts mean significant difference (P<0.05), while with no letter or the same small letter superscripts mean no significant difference (P>0.05).

由表2可知，3個(gè)品種水牛乳的乳脂肪、乳蛋白、乳糖、乳總固形物含量和乳能量差異均不顯著(P>0.05)；大理雜交水牛乳的乳尿素氮含量分別與騰沖檳榔江、德宏雜交水牛乳的乳尿素氮含量差異顯著(P<0.05)，但騰沖檳榔江水牛乳的乳尿素氮含量與德宏雜交水牛乳的乳尿素氮含量差異不顯著(P>0.05)。

表2 不同品種水牛乳成分比較

同列數(shù)據(jù)肩標(biāo)不同小寫字母表示差異顯著(P<0.05)，無(wú)字母或相同小寫字母表示差異不顯著(P>0.05)。

In the same column, values with different small letter superscripts mean significant difference (P<0.05), while with no letter or the same small letter superscripts mean no significant difference (P>0.05).

數(shù)據(jù)無(wú)字母標(biāo)注表示差異不顯著(P>0.05)。

Data noted with no letter mean no significant difference (P>0.05).

圖1 不同月份總樣本的乳成分比較

Fig.1 Comparison of milk composition of total samples in different months

由圖1可知，3個(gè)群體304個(gè)總樣本的乳脂肪、乳蛋白、乳糖、乳尿素氮、乳總固形物含量和乳能量在月份間差異均不顯著(P>0.05)。

2.2 水牛乳能量與乳成分的相關(guān)性分析

由表3可知，乳能量分別與乳脂肪和乳蛋白含量呈現(xiàn)極顯著正相關(guān)(r=0.896 0、r=0.563 0，P<0.01)。

由圖2可知，乳能量受乳脂肪及乳蛋白含量的影響較大。

2.3 回歸方程預(yù)測(cè)乳能量的顯著性和擬合度分析

由表4可知，試驗(yàn)所得水牛乳能量方程計(jì)算結(jié)果與中國(guó)奶牛飼養(yǎng)標(biāo)準(zhǔn)和中國(guó)奶牛飼養(yǎng)標(biāo)準(zhǔn)科研協(xié)作組方程計(jì)算結(jié)果[3-4]接近，分別低于Musgrave等[1]、AFRC[5]和NRC[6]的方程結(jié)果0.082 6、0.041 0和0.035 6 MJ/kg，但差異不大。

由圖3～圖6可知，回歸方程經(jīng)統(tǒng)計(jì)學(xué)檢驗(yàn)達(dá)差異極顯著水平(P<0.01)。分別以乳脂肪(F)，乳脂肪和乳蛋白(F和P)，乳脂肪、乳蛋白和乳糖(F、P和La)，乳脂肪、乳蛋白、乳糖和乳尿素氮含量(F、P、La和MUN)為預(yù)測(cè)因子，建立的預(yù)測(cè)乳能量(E)的一元、二元、三元及四元回歸方程的擬合度均在0.90以上，方程分別為：E=0.388F+1.540 (R2=0.933 6，殘差范圍：-1.113 6～1.466 0，P<0.01，n=304)；E=0.373F+0.221P+0.460 (R2=0.926 7，殘差范圍：-1.125 0～1.574 1，P<0.01，n=304)；E=0.396F+0.186P+0.105La-0.104 (R2=0.954 0，殘差范圍：-0.998 1～1.560 9，P<0.01，n=304)；E=0.397F+0.187P+0.106La+0.002MUN-0.146 (R2=0.958 0，殘差范圍：-1.007 0～1.547 2，P<0.01，n=304)。

表3 水牛乳乳能量與乳成分相關(guān)性分析

*表示相關(guān)系數(shù)在0.05水平上顯著(雙尾)，**表示相關(guān)系數(shù)在0.01水平上極顯著(雙尾)。

* mean correlation is significant at the 0.05 level (2-tailed), and ** mean correlation is significant at the 0.01 level (2-tailed).

從擬合度及乳能量的殘差(實(shí)測(cè)值-預(yù)測(cè)值)可得，分別以乳脂肪、乳蛋白和乳糖，乳脂肪、乳蛋白、乳糖和乳尿素氮含量為預(yù)測(cè)因子的三元和四元回歸方程的擬合度較好，預(yù)測(cè)結(jié)果較為準(zhǔn)確，在只有乳脂肪和乳蛋白含量的基礎(chǔ)上，也可用回歸方程預(yù)測(cè)產(chǎn)奶凈能。

3 討 論

3.1 乳成分比較及影響因素

不同月份騰沖檳榔江水牛乳的乳脂肪、乳蛋白、乳糖、乳尿素氮、乳總固形物含量和乳能量差異均不顯著，這與整個(gè)樣本采集期的飼糧組成及飼喂水平無(wú)改變有關(guān)；大理雜交水牛乳和德宏雜交水牛乳中的乳蛋白、乳糖、乳總固形物含量和乳能量差異均不顯著，這是因?yàn)槿榈鞍缀腿樘呛渴芗竟?jié)影響不大[5，11-15]。不同品種的水牛乳的乳脂肪、乳蛋白、乳糖、乳總固形物含量差異不顯著，這是因?yàn)槠贩N對(duì)乳糖和乳蛋白含量的影響非常小[5，11-15]。乳脂肪含量的均值為6.45%，低于謝紅等[16]、Sun等[17]、Rafiq等[18]的研究結(jié)果，高于Islam等[19]、Javed等[20]的研究結(jié)果，這可能與飼糧中的粗脂肪含量及所采乳樣為放乳前、中和后期或前、中和后期混合樣密切相關(guān)[21-22]。乳脂肪含量在3.30%～11.10%，在前人研究報(bào)道[11，20，23-25]的范圍內(nèi)。乳蛋白含量的均值為4.55%，與謝紅等[16]、Javed等[20]及鄒彩霞等[23]報(bào)道一致；乳蛋白含量在2.84%～6.58%，在Nasr[11]、Javed等[20]、鄒彩霞等[23]、Claeys等[24]及Zotos等[25]報(bào)道的范圍內(nèi)。乳糖含量在3.09%～5.99%，在謝紅等[16]和鄒彩霞等[23]的研究范圍內(nèi)。乳尿素氮含量范圍較廣，在1.40～26.60 mg/dL，且大理雜交水牛乳與德宏雜交水牛乳的乳尿素氮和乳脂肪含量呈現(xiàn)相同規(guī)律，與Rajala-Schultz等[26]研究的結(jié)果一致，乳尿素氮含量與乳脂肪含量呈正相關(guān)，且在9月份與12月份出現(xiàn)差異顯著，以及大理雜交水牛乳中乳尿素氮含量分別與騰沖檳榔江水牛乳、德宏雜交水牛乳中乳尿素氮含量差異顯著，這都與飼糧組成及動(dòng)物個(gè)體氮的代謝密切相關(guān)[27]。乳總固形

物含量在11.42%～25.93%，在Nasr[11]、Rafiq等[18]報(bào)道的范圍內(nèi)。乳能量在2.01～6.16 MJ/kg，均值為4.02 MJ/kg，高于荷斯坦牛奶的乳能量2.84 MJ/kg，低于鄒彩霞等[23]和Claeys等[24]報(bào)道的4.24～4.78 MJ/kg和4.20～4.70 MJ/kg，這與乳成分受品種、牛的大小、溫度、環(huán)境條件、管理和衛(wèi)生的影響有關(guān)[11，28-30]。

圖2 乳脂肪與乳蛋白和乳能量的三維圖

公式Formula乳能量Milkenergy/(MJ/kg)E=0.388F+1.5403.0920E=0.4344F+1.437[1]3.1746E=0.406F+1.509[6]3.1330E=0.4054F+1.506[5]3.1276E=0.4153F+1.4337[4]3.0949E=0.4153F+1.4336[3]3.0948

E：乳能量 milk energy；F：4%校正乳 4% fat corrected milk。

3.2 乳成分與乳能量相關(guān)性分析

乳能量分別與乳脂肪和乳蛋白含量呈現(xiàn)極顯著正相關(guān)，與Musgrave等[1]得到的乳脂肪含量與乳熱量值存在顯著直線關(guān)系一致。但乳脂肪含量與乳能量的相關(guān)性低于Tyrrell等[2]的結(jié)果，這可能是與品種有關(guān)，水牛乳有高的乳尿素氮含量，其對(duì)乳能量有一定的貢獻(xiàn)，從而降低了乳脂肪含量與乳能量的相關(guān)性。乳能量受乳脂肪及乳蛋白含量的影響較大，乳中泌乳凈能等于乳中各組分燃燒熱值的總能，已報(bào)道的乳脂、乳蛋白及乳糖的燃燒熱分別為38.87、23.89和16.53 MJ/kg，而乳糖的變異程度較小[5]。

圖3 以乳脂肪含量為預(yù)測(cè)因子的乳能量預(yù)測(cè)值及殘差的關(guān)系散點(diǎn)圖

圖4 以乳脂肪和乳蛋白含量為預(yù)測(cè)因子的乳能量預(yù)測(cè)值及殘差的關(guān)系散點(diǎn)圖

圖5 以乳脂肪、乳蛋白和乳糖含量為預(yù)測(cè)因子的乳能量預(yù)測(cè)值及殘差的關(guān)系散點(diǎn)圖

圖6 以乳脂肪、乳蛋白、乳糖和乳尿素氮含量為預(yù)測(cè)因子的乳能量預(yù)測(cè)值及殘差的關(guān)系散點(diǎn)圖

3.3 回歸方程預(yù)測(cè)乳能量的分析

早在1928年Gaines提出用含脂4%的校正乳避免牛奶成分的改變而干擾乳成分與乳能量的正常關(guān)系，Musgrave等[1]認(rèn)為4%的校正乳方程將調(diào)整牛奶產(chǎn)量的恒定能量值。試驗(yàn)所得乳脂率為4%的水牛乳能量為3.092 0 MJ/kg，在報(bào)道的4%乳脂的校正乳能量為3.05～3.14 MJ/kg內(nèi)[3]，與中國(guó)奶牛飼養(yǎng)標(biāo)準(zhǔn)和中國(guó)奶牛飼養(yǎng)標(biāo)準(zhǔn)科研協(xié)作組方程計(jì)算結(jié)果接近[3-4]。試驗(yàn)所得分別以乳脂肪、乳蛋白，乳脂肪、乳蛋白、乳糖，乳脂肪、乳蛋白、乳糖、乳尿素氮含量為預(yù)測(cè)因子的二元、三元及四元回歸方程分別與英國(guó)AFRC(1993)奶牛飼養(yǎng)標(biāo)準(zhǔn)中乳能量預(yù)測(cè)方程[6]，NRC(2001)奶牛飼養(yǎng)標(biāo)準(zhǔn)中乳能量預(yù)測(cè)方程[5]，中國(guó)奶牛飼養(yǎng)標(biāo)準(zhǔn)(2004)中乳能量預(yù)測(cè)方程[4]及中國(guó)奶牛飼養(yǎng)標(biāo)準(zhǔn)科研協(xié)作組(1986)中乳能量預(yù)測(cè)方程[3]中的乳脂肪、乳蛋白和乳糖的系數(shù)相差不大，較接近于所報(bào)道的乳脂肪的燃燒值0.388 MJ/kg，而與乳蛋白及乳糖的燃燒值分別為(0.238和0.165 MJ/kg)有所差異，這和乳能量受乳脂肪、乳蛋白和乳糖含量共同影響有關(guān)[5]。而其方程中的乳脂肪、乳蛋白和乳糖的系數(shù)均與鄒彩霞等[23]乳能量預(yù)測(cè)方程中的乳脂肪、乳蛋白和乳糖的系數(shù)差異較大，其原因還需進(jìn)一步分析。

預(yù)測(cè)方程式的殘差范圍分別為-1.113 6～1.466 0、-1.125 0～1.574 1、-0.998 1～1.560 9和-1.007 0～1.547 2，以乳脂肪、乳蛋白、乳糖，乳脂肪、乳蛋白、乳糖、乳尿素氮含量為預(yù)測(cè)因子的三元及四元回歸方程的殘差范圍較小，其模型假設(shè)的合理性及數(shù)據(jù)的可靠性更高[31]。

4 結(jié) 論

可通過(guò)水牛乳中的乳脂肪(F)、乳蛋白(P)、乳糖(La)及乳尿素氮(MUN)含量預(yù)測(cè)水牛乳產(chǎn)奶凈能(E),方程分別為：E=0.388F+1.540 (R2=0.933 6，P<0.01)；E=0.373F+0.221P+0.460 (R2=0.926 7，P<0.01)；E=0.396F+0.186P+0.105La-0.104 (R2=0.954 0，P<0.01)；E=0.397F+0.187P+0.106La+0.002MUN-0.146 (R2=0.958 0，P<0.01)。

[1] MUSGRAVE S D,SALISBURY G W.The relationship of milk energy and total protein to per cent fat in brown Swiss herd milk[J].Journal of Dairy Science,1952,35(2):174-178.

[2] TYRRELL H F,REID J T.Prediction of the energy value of cow，s milk[J].Journal of Dairy Science,1965,48(9):1215-1223.

[3] 中國(guó)奶牛飼養(yǎng)標(biāo)準(zhǔn)科研協(xié)作組.中國(guó)黑白花奶牛乳中各種營(yíng)養(yǎng)成分及能量間相關(guān)性的研究[J].中國(guó)畜牧雜志,1987(2):5-7.

[4] 中華人民共和國(guó)農(nóng)業(yè)部.NY/T 34—2004奶牛飼養(yǎng)標(biāo)準(zhǔn)[S].北京:中國(guó)農(nóng)業(yè)出版社,2004.

[5] [美]國(guó)家科學(xué)研究委員會(huì).奶牛營(yíng)養(yǎng)需要[M].孟慶翔,譯.北京:中國(guó)農(nóng)業(yè)大學(xué)出版社,2002.

[6] ALDERMAN G.Energy and protein requirements of ruminants[M].Wallingford:CABI Publishing,1993.

[7] FRIGGENS N C,RIDDER C,L?VENDAHL P.On the use of milk composition measures to predict the energy balance of dairy cows[J].Journal of Dairy Science,2007,90(12):5453-5467.

[8] BRUNLAFLEUR L,DELABY L,HUSSON F.Predicting energy×protein interaction on milk yield and milk composition in dairy cows[J].Journal of Dairy Science,2010,93(9):4128-4143.

[9] 聯(lián)合國(guó)糧農(nóng)組織統(tǒng)計(jì)數(shù)據(jù)庫(kù)[EB/OL].(2001-12-19)[2017-01-01].http://faostat3.fao.org/compare/Q/QA/E.

[10] 簡(jiǎn)保權(quán),秦學(xué)敏,龔芳.世界水牛奶業(yè)發(fā)展現(xiàn)狀和典型模式分析[J].世界農(nóng)業(yè),2015(3):115-118.

[11] NASR M A.The impact of crossbreeding Egyptian and Italian buffalo on milk yield and composition under subtropical environmental conditions[J].Journal of Dairy Research ,2016,83(2):196-201.

[12] JENKINS T C,MCGUIRE M A.Major advances in nutrition:impact on milk composition[J].Journal of Dairy Science,2006,89(4):1302-1310.

[13] 塔娜,吉日木圖.日糧中添加保護(hù)性脂肪對(duì)奶牛產(chǎn)奶量及乳成分的影響[J].中國(guó)奶牛,2008(9):13-14.

[14] REN D X,ZOU C X,LIN B,et al.A comparison of milk protein,amino acid and fatty acid profiles of river buffalo and their F1 and F2 hybrids with swamp buffalo in China[J].Pakistan Journal of Zoology,2015,47(5):1459-1465.

[15] HECK J M L,VAN VALENBERG H J F,DIJKSTRA J,et al.Seasonal variation in the Dutch bovine raw milk composition[J].Journal of Dairy Science,2009,92(10):4745-4755.

[16] 謝紅,劉琴,周亞平.云南省奶水牛乳成分與產(chǎn)奶量之間的關(guān)系[J].中國(guó)奶牛,2011(20):54-57.

[17] SUN Q,LV J P,LIU L,et al.Comparison of milk samples collected from some buffalo breeds and crossbreeds in China[J].Dairy Science & Technology,2014,94(4):387-395.

[18] RAFIQ S, HUMA N, PASHA I,et al.Chemical composition,nitrogen fractions and amino acids profile of milk from different animal species[J].Asian-Australasian Journal of Animal Sciences,2016,29(7):1022-1028.

[19] ISLAM M A,ALAM M K,ISLAM M N,et al.Principal milk components in buffalo,holstein cross,indigenous cattle and red chittagong cattle from Bangladesh[J].Asian-Australasian Journal of Animal Sciences,2014,27(6):886-897.

[20] JAVED K,ABDULLAH M,KHALID M S, et al.Inter-relationship of milk constituents with body and udder measurements in NiLi-Ravi buffaloes raised at commercial farms of Pakistan[J].Buffalo Bulletin,2013,32(2):1170-1173.

[21] 劉洪生.影響乳成分的日糧調(diào)控技術(shù)[J].中國(guó)奶牛,2008(3):22-24.

[22] 王雅春,楊雪麗,張淑娟.奶牛放乳前中后期乳成分變化的研究[J].中國(guó)奶牛,1997(4):19-20.

[23] 鄒彩霞,楊炳壯,韋升菊,等.水牛乳能值與乳成分的相關(guān)性研究[J].動(dòng)物營(yíng)養(yǎng)學(xué)報(bào),2010,22(4):964-968.

[24] CLAEYS W L,VERRAES C,CARDOEN S,et al.Consumption of raw or heated milk from different species:an evaluation of the nutritional and potential health benefits[J].Food Control,2014,42:188-201.

[25] ZOTOS A,BAMPIDIS V A.Milk fat quality of Greek buffalo (Bubalusbubalis)[J].Journal of Food Composition and Analysis,2014,33(2):181-186.

[26] RAJALA-SCHULTZ P J,SAVILLE W J A.Sources of variation in milk urea nitrogen in Ohio dairy herds[J].Journal Dairy Science,2003,86(5):1653-1661.

[27] 任春燕,周凌云,卜登攀,等.牛奶尿素氮的主要影響因素及其在奶牛生產(chǎn)中的應(yīng)用[J].中國(guó)畜牧獸醫(yī),2012,39(9):129-133.

[28] FOODA T A,MOURAD K A,MOURAD G.Phenotypic and genetic trends for milk production in Egyptian buffaloes[J].Journal of American Science,2010,6(11):143-147.

[30] BHUTKAR S S,THOMBRE B M,BAINWAD D V.Effect of non-genetic factors on production traits in Deoni cows[J].Journal of Agriculture and Veterinary Science,2014,7(12):15-19.

[31] 馮巖松.SPSS 22.0統(tǒng)計(jì)分析應(yīng)用[M].北京:清華大學(xué)出版社,2015:6.

*Corresponding authors: MAO Huaming, professor, E-mail: maohm@vip.sina.com; LI Wen, associate professor, E-mail: 331631693@qq.com

(責(zé)任編輯 武海龍)

Establishment of Prediction Model of Milk Energy and Milk Composition Analysis ofYunnanBuffalo

LI Qing1,2MAO Huaming1,2*LI Wen1,3*LI Jianhua4LU Qiongfen1HE Hongyuan1LIU Qin5

(1.CollegeofAnimalScienceandTechnology,YunnanAgriculturalUniversity,Kunming650201,China; 2.YunnanProvincialKeyLaboratoryofAnimalNutritionandFeedScience,Kunming650201,China; 3.CollegeofLifeScienceandChemistry,HunanUniversityofTechnology,Zhuzhou412007,China; 4.AnimalHusbandryandVeterinaryStationofAnning,Kunming650300,China; 5.DairyProductionPerformanceMeasurementCenterofKunming,Kunming650041,China)

This experiment was conducted to analysis the contents of milk fat, milk protein, milk lactose and milk urea nitrogen of buffalo milk and those correlation with milk energy, and to established the prediction model of milk production net energy for scientific feeding management of lactating dairy buffalo. A total of 304 milk samples were collected from the main milk buffalo farms and farming communities inYunnanprovince at 2014 to 2015 years 9 to 12 months, which were sent toKunmingMilk Production Center. The milk composition of the buffalo milk was detected by the MilkoScan FT+FC and MilkoScan FT 120 analyzer in Kunming Dairy Production Center, and milk energy was determined by oxygen bomb calorimeter after constant temperature drying in vacuum drying oven. The results showed that milk fat, milk protein, milk lactose, milk urea nitrogen, milk total solids and milk energy were 6.45%, 4.55%, 5.31%, 13.60 mg/dL, 18.77% and 4.02 MJ/kg, respectively. And there was significant positive correlation between milk energy and milk fat, milk protein (r=0.896 0,r=0.563 0,P<0.01). The milk energy was greater impacted by milk fat and milk protein. Moreover, the fitting degrees of the one, two, three and four regression equations were above 0.90 to predict the net energy (E) of buffalo milk by using milk fat (F), milk fat and milk protein (FandP), milk fat, milk protein and milk lactose (F,PandLa), milk fat, milk protein, milk lactose and milk urea nitrogen contents (F,P,LaandMUN) as predictors, and regression equations were as followed:E=0.388F+1.540 (R2=0.933 6,P<0.01);E=0.373F+0.221P+0.460 (R2=0.926 7,P<0.01);E=0.396F+0.186P+0.105La-0.104 (R2=0.954 0,P<0.01);E=0.397F+0.187P+0.106La+0.002MUN-0.146 (R2=0.958 0,P<0.01). In conclusion, the milk production net energy can be predicted by the contents of milk fat, milk protein, milk lactose and milk urea nitrogen of buffalo milk.[ChineseJournalofAnimalNutrition, 2017, 29(8)：2875-2883]

buffalo milk; milk composition; milk energy; prediction model

10.3969/j.issn.1006-267x.2017.08.032

2017-01-01

云南省生物重大專項(xiàng)——奶業(yè)專項(xiàng)(2014ZA002)；湖南省科技計(jì)劃項(xiàng)目(2016NK2096)；中國(guó)博士后面上項(xiàng)目(2016M592456)；長(zhǎng)沙市科技局重點(diǎn)項(xiàng)目(kh16011067)

李 清(1971—)，女，廣東臺(tái)山人，副教授，博士，主要從事動(dòng)物營(yíng)養(yǎng)研究。E-mail: 814318012@qq.com

*通信作者:毛華明，教授，博士生導(dǎo)師，E-mail: maohm@vip.sina.com；李 文，副教授，E-mail: 331631693@qq.com

S823

A

1006-267X(2017)08-2875-09