标题

Transcription

1 动物营养学报 2017,29(8): Chinese Journal of Animal Nutrition doi: / j.issn x 云南水牛乳成分分析及乳能量预测模型的建立 李 1,2 清 1,2 毛华明 李 1,3 文 4 李建华 1 鲁琼芬 1 何鸿源 (1. 云南农业大学动物科学技术学院, 昆明 ;2. 云南省动物营养与饲料重点实验室, 昆明 ; 3. 湖南工业大学生命科学与化学学院, 株洲 ;4. 云南省安宁市畜牧兽医站, 昆明 ;5. 云南省昆明市奶牛生产性能测定中心, 昆明 ) 刘 5 琴 摘要 : 本试验旨在分析水牛乳成分中乳脂肪 乳蛋白 乳糖和乳尿素氮含量及其与乳能量的相关性, 建立产奶净能的预测模型, 用于奶水牛泌乳期科学的饲养管理 试验采取云南省主要奶水牛养殖场和养殖小区的 年 9 12 月份送达昆明市奶牛生产性能测定中心的 304 个原料乳样本, 水牛乳的乳成分由昆明市奶牛生产性能测定中心用 MilkoScan FT+FC 乳成分体细胞联用仪和 MilkoScan FT 120 乳成分分析仪测定, 鲜乳经真空干燥箱恒温烘干后用氧弹式热量计测定乳能量 结果表明 : 乳脂肪 乳蛋白 乳糖 乳尿素氮 乳总固形物含量和乳能量分别为 6.45% 4.55% 5.31% mg / dl 18.77% 和 4.02 MJ / kg; 乳能量分别与乳脂肪和乳蛋白含量呈现极显著正相关 ( r = r = ,P<0.01), 乳能量受乳脂肪和乳蛋白含量的影响较大 分别以乳脂肪 ( F), 乳脂肪和乳蛋白 ( F 和 P), 乳脂肪 乳蛋白和乳糖 ( F P 和 La), 乳脂肪 乳蛋白 乳糖和乳尿素氮含量 ( F P La 和 MUN) 为预测因子, 建立的预测乳能量 ( E) 的一元 二元 三元及四元回归方程的拟合度均在 0.90 以上, 方程分别为 :E = 0.388F ( R 2 = ,P< 0.01) ;E = F P ( R 2 = ,P < 0. 01) ;E = F P La ( R 2 = ,P<0.01) ;E = 0.397F+0.187P+0.106La+0.002MUN ( R 2 = ,P<0.01) 由此可见, 可通过水牛乳中的乳脂肪 乳蛋白 乳糖及乳尿素氮含量预测水牛乳产奶净能 关键词 : 水牛乳 ; 乳成分 ; 乳能量 ; 预测模型中图分类号 :S823 文献标识码 :A 文章编号 : X(2017) 能量在动物生长发育 繁殖 生产过程中是必不可少的, 对于泌乳牛, 乳的能量即产奶净能 ( NE L ) 目前, 采用净能体系的国家统一用产奶净能 4% 乳脂率和产奶量来计算泌乳牛的能量需求 [1-6] 乳的能量来源于乳脂肪 乳蛋白 乳糖等乳成分, 但是在实际生产中, 测定牛乳中的能量相对复杂, 而牛乳中的营养成分则是每个养牛场每月必测的指标, 较容易获得, 乳脂肪和乳蛋白含量不仅可预测乳中能量值, 还可预测产奶牛的能量是否平衡 [7-8] 美国奶牛饲养标准 英国奶牛饲养 标准及中国奶牛饲养标准中通过测定牛乳中的营养成分含量, 推算出产奶净能的预测模型, 从而指导生产 [4-6] 在这些预测模型中所采集的牛奶样本均为荷斯坦血系的奶牛, 它们的乳脂肪 乳蛋白 乳糖及乳尿素氮含量均低于水牛乳中相应含量, 二者乳成分差异较大, 尤其是水牛乳中尿素氮的含量较高, 其在乳能量中是否有贡献? 对于奶水牛的产奶净能能否套用荷斯坦奶牛的预测模型值得考量 据联合国粮食及农业组织 ( FAO) 2014 年统 收稿日期 : 基金项目 : 云南省生物重大专项 奶业专项 (2014ZA002) ; 湖南省科技计划项目 ( 2016NK2096) ; 中国博士后面上项目 ( 2016M592456) ; 长沙市科技局重点项目 ( kh ) 作者简介 : 李清 (1971 ), 女, 广东台山人, 副教授, 博士, 主要从事动物营养研究 E mail: @ qq.com 通信作者 : 毛华明, 教授, 博士生导师,E mail: maohm@ vip.sina.com; 李文, 副教授,E mail: @ qq.com

2 2876 动物营养学报 29 卷 计, 奶水牛全球存栏量为 万头, 鲜水牛乳占总奶量的比重已上升到 12.92% [9], 但水牛奶量单产低, 除了加强品种选育提高单产外, 合理的饲粮及科学的饲养管理则是提高单产的核心 [10] 但目前, 尚未有确定的奶水牛营养需要量指导生产, 因此获得奶水牛的营养需要量意义重大 本试验用 MilkoScan FT+FC 乳成分体细胞联用仪和 MilkoScan FT 120 乳成分分析仪测定水牛乳成分, 用氧弹式热量计测定乳能量, 探讨水牛乳能量与乳成分的相关性及通过回归分析建立乳能量的预测模型, 得到奶水牛的泌乳能量需要, 为奶水牛饲养标准的制定提供基础数据 1 材料与方法 1.1 试验材料为了使所测定样本具有代表性, 试验选择了云南省主要奶水牛养殖地不同的饲养模式及饲粮水平的 3 个奶水牛群体进行分析 3 个群体的试验动物每天饲喂 2 次, 自由饮水, 每日挤奶 2 次 槟榔江水牛乳来源于腾冲槟榔江水牛核心养殖场, 产犊胎次主要为 2 ~ 4 胎 泌乳早中期健康的 184 头槟榔江水牛, 饲粮为自配精料 全株玉米青贮 稻草, 饲喂量为每头牛每天精料 3 ~ 4 kg( 根据产奶量 ) 玉米青贮 20 ~ 25 kg 干稻草自由采食 ( 摄入量 4 ~ 5 kg), 蛋白质水平为 12% 大理杂交水牛乳来自于大理杂交水牛养殖小区, 产犊胎次主要为 1 ~ 5 胎 泌乳早中期健康的 56 头大理杂交水牛, 饲粮为市售精料 玉米秸秆青贮 稻草 苜蓿干草 羊草, 饲喂量为每头牛每天精料 3 ~ 4 kg( 根据产奶量 ) 玉米秸秆青贮 15 ~ 20 kg 干稻草自由采食 ( 摄入量 4 ~ 5 kg), 蛋白质水平为 14%; 泌乳中期增加苜蓿干草和羊草共 2 kg 德宏杂交水牛乳来源于德宏杂交水牛养殖小区, 产犊胎次主要为 1 ~ 5 胎 泌乳早中期健康的 64 头德宏杂交水牛, 饲粮为自配精料 甘蔗梢 稻草 啤酒糟 王草, 饲喂量为每头牛每天精料 1.5 ~ 2.0 kg( 根据产奶量 ) 甘蔗梢 20 kg 稻草和王草自由采食( 摄入量 4 ~ 5 kg), 蛋白质水平为 11%; 泌乳中期增加啤酒糟 2.0 kg 1.2 样品采集及测定分别采集云南省腾冲 大理和德宏 3 个主要奶水牛养殖场和养殖小区的 年 9 12 月份的 304 份牛奶样品进行乳成分含量及乳能量 的测定 采样时间均为每月的月底, 将每天 2 次 的牛奶样品按 1 1 充分混匀后倒入采样瓶摇匀, 确 保防腐剂溶于牛奶样品, 每头牛采样 40 ~ 50 ml, 分成 2 份, 将样品低温保存迅速运输至实验室 2 ~ 7 可保存 1 周 1.3 乳成分含量及乳能量的测定 乳成分含量的测定由昆明市奶牛生产性能测 定中心用 MilkoScan FT+FC 乳成分体细胞联用仪 和 MilkoScan FT 120 乳成分分析仪测定完成 乳能量测定在云南农业大学用 BH ⅢS 氧弹 式热量计测定 : 移取 5 ml 牛奶样品放入已知重量 的自封袋中称重 ( 自封袋置于古氏坩埚中 ), 放入 真空干燥箱 (60 ±5) 干燥后用氧弹式热量计测 定乳能量 1.4 数据处理 试验数据整理及作图用 Excel 2016 和 SPSS 22.0 软件进行, 差异显著性检验采用单因素方差 分析 ( one way ANOVA) 和 Duncan 氏多重比较 法, 由模型 Y i = μ+t i +ε i ( Y i 为观测值,μ 为总体平 均值,T i 为处理效应 ;ε i 为随机误差 ) 获得, 结果用 平均值 ± 标准差表示 相关性分析采用二元变量 相关分析中的 Pearson 相关系数判定 ; 相关分析模 型采用多元线性回归分析得到, 用模型 E = β 0 + β 1 F+β 2 P + β 3 La + β 4 MUN [ β 0 为常数 ; β 1 β 2 β 3 β 4 为偏回归系数 ;F 为乳脂肪含量 ( %) ;P 为乳蛋白 含量 ( %) ;La 为乳糖含量 ( %) ;MUN 为乳尿素氮 含量 ( mg / dl) ;E 为乳能量 ( MJ / kg) ] 表示 2 结果 2.1 水牛乳成分比较 由表 1 可知, 乳脂肪 乳蛋白 乳糖 乳尿素 氮 乳总固形物含量和乳能量的平均值分别为 6.45% 4. 55% 5. 31% mg / dl % 和 4.02 MJ / kg 乳脂肪含量在 3.30% ~ 11.10%, 且随 泌乳时间的变化而有所上升 ; 乳蛋白含量在 2.84% ~ 6.58%; 乳糖含量在 3.09% ~ 5.99%; 乳尿 素氮含量范围较广, 在 1.40 ~ mg / dl; 乳总固 形物含量在 11.42% ~ 25.93% 同一品种不同月份间水牛乳成分比较, 腾冲 槟榔江水牛乳的乳脂肪 乳蛋白 乳糖 乳尿素氮 乳总固形物含量和乳能量差异均不显著 ( P > 0.05) ; 大理杂交水牛乳和德宏杂交水牛乳中的乳 蛋白 乳糖 乳总固形物含量和乳能量差异均不显

3 8 期李清等 : 云南水牛乳成分分析及乳能量预测模型的建立 2877 著 ( P>0.05), 但二者 9 月份的乳脂肪 乳尿素氮含量均与 12 月份差异显著 ( P<0.05) Table 1 表 1 同一品种不同月份水牛乳成分比较 Comparison of buffalo milk composition in different months of the same breed 品种 Breeds 腾冲槟榔江水牛乳 Milk of Tengchong BLJ buffalo 大理杂交水牛乳 Milk of Dali hybrid buffalo 德宏杂交水牛乳 Milk of Dehong hybrid buffalo 月份 Months 总样本范围 Total sample range ( n = 304) 总样本平均值 Total sample mean ( n = 304) 样本数 Sample number 乳脂肪 Milk fat / % 乳蛋白 Milk protein / % 乳糖 Lactose / % 乳尿素氮 Milk urea nitrogen / ( mg / dl) 乳总固形物 Milk total solid / % 乳能量 Milk energy / ( MJ / kg) ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ±1.56 b 4.76± ± ±2.79 b 19.57± ± ±1.37 ab 4.78± ± ±2.51 ab 18.84± ± ±1.44 ab 4.22± ± ±2.20 ab 19.79± ± ±2.34 a 4.32± ± ±2.42 a 19.22± ± ±0.51 b 4.03± ± ±3.60 b 18.00± ± ±0.85 ab 3.92± ± ±3.40 ab 18.71± ± ±1.09 ab 3.98± ± ±2.80 ab 17.96± ± ±0.89 a 4.30± ± ±2.29 a 18.83± ± ~ ~ ~ ~ ~ ~ ± ± ± ± ± ±0.64 同一品种同列数据肩标不同小写字母表示差异显著 ( P<0.05), 无字母或相同小写字母表示差异不显著 ( 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 个品种水牛乳的乳脂肪 乳蛋白 乳糖 乳总固形物含量和乳能量差异均不显著 ( P>0.05) ; 大理杂交水牛乳的乳尿素氮含量分别与腾冲槟榔江 德宏杂交水牛乳的乳尿素氮含量 差异显著 ( P<0.05), 但腾冲槟榔江水牛乳的乳尿 素氮含量与德宏杂交水牛乳的乳尿素氮含量差异 不显著 ( P>0.05) Table 2 表 2 不同品种水牛乳成分比较 Comparison of different breeds buffalo milk composition 品种 Breeds 腾冲槟榔江水牛乳 Milk of Tengchong BLJ buffalo ( n = 184) 大理杂交水牛乳 Milk of Dali hybrid buffalo ( n = 56) 德宏杂交水牛乳 Milk of Dehong hybrid buffalo ( n = 64) 乳脂肪 Milk fat / % 乳蛋白 Milk protein / % 乳糖 Lactose / % 乳尿素氮 Milk urea nitrogen / ( mg / dl) 乳总固形物 Milk total solid / % 乳能量 Milk energy / ( MJ / kg) 6.40± ± ± ±4.97 bc 18.57± ± ± ± ± ±2.47 a 19.35± ± ± ± ± ±3.98 bc 18.38± ±0.63 同列数据肩标不同小写字母表示差异显著 ( P<0.05), 无字母或相同小写字母表示差异不显著 ( 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).

4 2878 动物营养学报 29 卷 数据无字母标注表示差异不显著 ( 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 个群体 304 个总样本的乳脂肪 乳蛋白 乳糖 乳尿素氮 乳总固形物含量和乳能量在月份间差异均不显著 ( P>0.05) 2.2 水牛乳能量与乳成分的相关性分析由表 3 可知, 乳能量分别与乳脂肪和乳蛋白含量呈现极显著正相关 ( r = r = , P<0.01) 由图 2 可知, 乳能量受乳脂肪及乳蛋白含量 的影响较大 2.3 回归方程预测乳能量的显著性和拟合度分析 由表 4 可知, 试验所得水牛乳能量方程计算 结果与中国奶牛饲养标准和中国奶牛饲养标准科 [3-4] 研协作组方程计算结果接近, 分别低于 Mus grave 等 [1] AFRC [5] 和 NRC [6] 的方程结果 和 MJ / kg, 但差异不大 由图 3 ~ 图 6 可知, 回归方程经统计学检验达 差异极显著水平 ( P < 0.01) 分别以乳脂肪 ( F), 乳脂肪和乳蛋白 ( F 和 P), 乳脂肪 乳蛋白和乳糖 ( F P 和 La), 乳脂肪 乳蛋白 乳糖和乳尿素氮含 量 ( F P La 和 MUN) 为预测因子, 建立的预测乳 能量 ( E) 的一元 二元 三元及四元回归方程的拟 合度均在 以上, 方程分别为 : E = F ( R 2 = , 残差范围 : ~ ,P < 0. 01, n = 304) ; E = F P ( R 2 = , 残差范围 : ~ , P<0. 01,n = 304) ;E = F P La (R 2 = , 残差范围 : ~ ,P< 0.01,n = 304 ) ; E = F P La MUN ( R 2 = , 残差范围 : ~ ,P<0.01,n = 304) Table 3 表 3 水牛乳乳能量与乳成分相关性分析 Analysis of correlation between energy and milk composition of buffalo milk 指标 乳脂肪 乳蛋白 乳糖 乳尿素氮 乳能量 Items Milk fat Milk protein Milk lactose Milk urea nitriogen Milk energy 乳脂肪 Milk fat 乳蛋白 Milk protein 乳糖 Milk lactose 乳尿素氮 Milk urea nitriogen 表示相关系数在 0.05 水平上显著 ( 双尾 ), 表示相关系数在 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). 从拟合度及乳能量的残差 ( 实测值 - 预测值 ) 可得, 分别以乳脂肪 乳蛋白和乳糖, 乳脂肪 乳蛋白 乳糖和乳尿素氮含量为预测因子的三元和四元回归方程的拟合度较好, 预测结果较为准确, 在只有乳脂肪和乳蛋白含量的基础上, 也可用回归方程预测产奶净能 3 讨论 3.1 乳成分比较及影响因素不同月份腾冲槟榔江水牛乳的乳脂肪 乳蛋白 乳糖 乳尿素氮 乳总固形物含量和乳能量差异均不显著, 这与整个样本采集期的饲粮组成及饲喂水平无改变有关 ; 大理杂交水牛乳和德宏杂交水牛乳中的乳蛋白 乳糖 乳总固形物含量和乳

5 ８期 李 清等 云南水牛乳成分分析及乳能量预测模型的建立 能量差异均 不 显 著 这 是 因 为 乳 蛋 白 和 乳 糖 含 量 受季节影响不 大 ５ １１ １５ 不 同 品 种 的 水 牛 乳 的 乳 ２８７９ 物 含 量 在 １１ ４２ ２５ ９３ 在 Ｎａｓｒ １１ Ｒａｆｉｑ 等 １８ 报 道 的 范 围 内 乳 能 量 在 ２ ０１ 脂肪 乳蛋白 乳糖 乳总固形物含量差异不显著 ６ １６ ＭＪ ｋｇ 均值为 ４ ０２ ＭＪ ｋｇ 高于荷斯坦牛奶 这是因 为 品 种 对 乳 糖 和 乳 蛋 白 含 量 的 影 响 非 常 的乳能量 ２ ８４ ＭＪ ｋｇ 低于邹彩霞等 ２３ 和 Ｃｌａｅｙｓ 小 ５ １１ １５ 乳脂肪含量的均值为 ６ ４５ 低于谢红 等 １６ Ｓｕｎ 等 １７ Ｒａｆｉｑ 等 １８ 的研究结果 高于 Ｉｓ ｌａｍ 等 １９ Ｊａｖｅｄ 等 ２０ 的研究结果 这可能与饲粮 中的粗脂肪 含 量 及 所 采 乳 样 为 放 乳 前 中 和 后 期 等 ２４ 报 道 的 ４ ２４ ４ ７８ ＭＪ ｋｇ 和 ４ ２０ ４ ７０ ＭＪ ｋｇ 这与乳成分受品种 牛的大小 温度 环境条件 管理和卫生的影响有关 １１ ２８ ３０ 或前 中和后期混合样密切相 关 ２１ ２２ 乳 脂 肪 含 量 在 ３ ３０ １１ １０ 在 前 人 研 究 报 道 １１ ２０ ２３ ２５ 的范围 内 乳 蛋 白 含 量 的 均 值 为 ４ ５５ 与 谢 红 等 １６ Ｊａｖｅｄ 等 ２０ 及 邹 彩 霞 等 ２３ 报 道 一 致 乳 蛋 白含量在 ２ ８４ ６ ５８ 在 Ｎａｓｒ １１ Ｊａｖｅｄ 等 ２０ 邹彩霞等 ２３ Ｃｌａｅｙｓ 等 ２４ 及 Ｚｏｔｏｓ 等 ２５ 报道的范 围内 乳糖含量在 ３ ０９ ５ ９９ 在谢红等 １６ 和 邹彩霞等 ２３ 的研究范围内 乳尿素氮含量范围较 广 在 １ ４０ ２６ ６０ ｍｇ ｄｌ 且大理杂交水牛乳与德 宏杂交水牛乳的乳尿素氮和乳脂肪含量呈现相同 规律 与 Ｒａｊａｌａ Ｓｃｈｕｌｔｚ 等 ２６ 研 究 的 结 果 一 致 乳 尿素氮含量 与 乳 脂 肪 含 量 呈 正 相 关 且 在 ９ 月 份 与 １２ 月份出现差异显著 以及大理杂交水牛乳中 图２ 乳尿素氮含 量 分 别 与 腾 冲 槟 榔 江 水 牛 乳 德 宏 杂 Ｆｉｇ ２ 交水牛乳中 乳 尿 素 氮 含 量 差 异 显 著 这 都 与 饲 粮 组成及动物个体氮的代谢密切相关 表４ Ｔａｂｌｅ ４ ２７ 乳脂肪与乳蛋白和乳能量的三维图 Ｔｈｅ ３Ｄ ｆｉｇｕｒｅ ｏｆ ｍｉｌｋ ｆａｔ ａｎｄ ｍｉｌｋ ｐｒｏｔｅｉｎ ｃｏｎｔｅｎｔｓ ｗｉｔｈ ｍｉｌｋ ｅｎｅｒｇｙ 乳总固形 不同乳能量预测公式在 ４ 校正乳下的乳能量比较 Ｃｏｍｐａｒｉｓｏｎ ｏｆ ｄｉｆｆｅｒｅｎｔ ｍｉｌｋ ｅｎｅｒｇｙ ｐｒｅｄｉｃｔｉｏｎ ｆｏｒｍｕｌａ ｉｎ ４ ｆａｔ ｃｏｒｒｅｃｔｅｄ ｍｉｌｋ 乳能量 Ｍｉｌｋ ｅｎｅｒｇｙ ＭＪ ｋｇ 公式 Ｆｏｒｍｕｌａ Ｅ ０ ３８８Ｆ １ ５４０ ３ ０９２ ０ Ｅ ０ ４３４ ４Ｆ １ ４３７ １ Ｅ ０ ４０６Ｆ １ ５０９ ６ ３ １７４ ６ Ｅ ０ ４０５ ４Ｆ １ ５０６ Ｅ ０ ４１５ ３Ｆ １ ４３３ ７ ４ Ｅ ０ ４１５ ３Ｆ １ ４３３ ６ ３ ３ １２７ ６ ３ １３３ ０ ５ ３ ０９４ ９ ３ ０９４ ８ Ｅ 乳能量 ｍｉｌｋ ｅｎｅｒｇｙ Ｆ ４ 校正乳 ４ ｆａｔ ｃｏｒｒｅｃｔｅｄ ｍｉｌｋ ３ ２ 乳成分与乳能量相关性分析 对乳能量有 一 定 的 贡 献 从 而 降 低 了 乳 脂 肪 含 量 乳能量分别与乳脂肪和乳蛋白含量呈现极显 著正相关 与 Ｍｕｓｇｒａｖｅ 等 １ 与乳能量的相关性 乳能量受乳脂肪及乳蛋白含 得到的乳脂肪含量与 量的影响较 大 乳 中 泌 乳 净 能 等 于 乳 中 各 组 分 燃 乳热量值存在显著直线关系一致 但乳脂肪含量 烧热值的总能 已报道的乳脂 乳蛋白及乳糖的燃 与乳能量的相关性低于 Ｔｙｒｒｅｌｌ 等 ２ 的 结 果 这 可 能是与品种有关 水牛乳有高的乳尿素氮含量 其 烧热分别为 ３８ ８７ ２３ ８９ 和 １６ ５３ ＭＪ ｋｇ 而 乳 糖 的变异程度较小 ５

6 ２８８０ 动 图３ 物 营 养 学 报 ２９ 卷 以乳脂肪含量为预测因子的乳能量预测值及残差的关系散点图 Ｆｉｇ ３ Ｒｅｌａｔｉｏｎｓｈｉｐ ｏｆ ｓｃａｔｔｅｒ ｄｉａｇｒａｍ ｂｅｔｗｅｅｎ ｐｒｅｄｉｃｔｅｄ ｖａｌｕｅ ｏｆ ｍｉｌｋ ｅｎｅｒｇｙ ａｎｄ ｒｅｓｉｄｕａｌ ｂｙ ｐｒｅｄｉｃｔｉｖｅ ｆａｃｔｏｒ ｏｆ ｍｉｌｋ ｆａｔ ｃｏｎｔｅｎｔ 图４ Ｆｉｇ ４ 以乳脂肪和乳蛋白含量为预测因子的乳能量预测值及残差的关系散点图 Ｒｅｌａｔｉｏｎｓｈｉｐ ｏｆ ｓｃａｔｔｅｒ ｄｉａｇｒａｍ ｂｅｔｗｅｅｎ ｐｒｅｄｉｃｔｅｄ ｖａｌｕｅ ｏｆ ｍｉｌｋ ｅｎｅｒｇｙ ａｎｄ ｒｅｓｉｄｕａｌ ｂｙ ｐｒｅｄｉｃｔｉｖｅ ｆａｃｔｏｒｓ ｏｆ ｍｉｌｋ ｆａｔ ａｎｄ ｍｉｌｋ ｐｒｏｔｅｉｎ ｃｏｎｔｅｎｔ 图５ Ｆｉｇ ５ 以乳脂肪 乳蛋白和乳糖含量为预测因子的乳能量预测值及残差的关系散点图 Ｒｅｌａｔｉｏｎｓｈｉｐ ｏｆ ｓｃａｔｔｅｒ ｄｉａｇｒａｍ ｂｅｔｗｅｅｎ ｐｒｅｄｉｃｔｅｄ ｖａｌｕｅ ｏｆ ｍｉｌｋ ｅｎｅｒｇｙ ａｎｄ ｒｅｓｉｄｕａｌ ｂｙ ｐｒｅｄｉｃｔｉｖｅ ｆａｃｔｏｒｓ ｏｆ ｍｉｌｋ ｆａｔ ｍｉｌｋ ｐｒｏｔｅｉｎ ａｎｄ ｍｉｌｋ ｌａｃｔｏｓｅ ｃｏｎｔｅｎｔ

7 ８期 李 图６ 清等 云南水牛乳成分分析及乳能量预测模型的建立 ２８８１ 以乳脂肪 乳蛋白 乳糖和乳尿素氮含量为预测因子的乳能量预测值及残差的关系散点图 Ｆｉｇ ６ Ｒｅｌａｔｉｏｎｓｈｉｐ ｏｆ ｓｃａｔｔｅｒ ｄｉａｇｒａｍ ｂｅｔｗｅｅｎ ｐｒｅｄｉｃｔｅｄ ｖａｌｕｅ ｏｆ ｍｉｌｋ ｅｎｅｒｇｙ ａｎｄ ｒｅｓｉｄｕａｌ ｂｙ ｐｒｅｄｉｃｔｉｖｅ ｆａｃｔｏｒｓ ｏｆ ｍｉｌｋ ｆａｔ ｍｉｌｋ ｐｒｏｔｅｉｎ ｍｉｌｋ ｌａｃｔｏｓｅ ａｎｄ ｍｉｌｋ ｕｒｅａ ｎｉｔｒｏｇｅｎ ｃｏｎｔｅｎｔ ３ ３ 回归方程预测乳能量的分析 元及四元回 归 方 程 的 残 差 范 围 较 小 其 模 型 假 设 早在 １９２８ 年 Ｇａｉｎｅｓ 提出用含脂 ４ 的校正乳 避免牛奶成分的改变而干扰乳成分与乳能量的正 常关系 Ｍｕｓｇｒａｖｅ 等 １ 认为 ４ 的校正乳方程将调 的合理性及数据的可靠性更高 ３１ ４ 结 整牛奶 产 量 的 恒 定 能 量 值 试 验 所 得 乳 脂 率 为 ４ 的水牛乳能量为 ３ ０９２ ０ ＭＪ ｋｇ 在报道的 ４ 乳脂的校 正 乳 能 量 为 ３ ０５ ３ １４ ＭＪ ｋｇ 内 ３ 与 中国奶牛饲养标准和中国奶牛饲养标准科研协作 组方程计算结果接近 ３ ４ 试验所得分别以乳脂 肪 乳蛋 白 乳 脂 肪 乳 蛋 白 乳 糖 乳 脂 肪 乳 蛋 白 乳糖 乳 尿 素 氮 含 量 为 预 测 因 子 的 二 元 三 元 及四元回归方程分别与英国 ＡＦＲＣ １９９３ 奶牛饲 养标准中乳能量 预 测 方 程 ６ ＮＲＣ ２００１ 奶 牛 饲 养标准中 乳 能 量 预 测 方 程 ５ 中 国 奶 牛 饲 养 标 准 ２００４ 中乳能量预测方程 ４ 及中国奶牛饲养标准 科研协作组 １９８６ 中乳能量预测方程 ３ 中的乳脂 乳糖的燃烧值分别为 ０ ２３８ 和 ０ １６５ ＭＪ ｋｇ 有所 ０ ９３３ ６ Ｐ ０ ０１ Ｅ ０ ３７３Ｆ ０ ２２１Ｐ ０ ４６０ Ｒ ２ ０ ９２６ ７ Ｐ ０ ０１ Ｅ ０ ３９６Ｆ ０ １８６Ｐ ０ １０５Ｌａ ０ １０４ Ｒ ２ ０ ９５４ ０ Ｐ ０ ０１ Ｅ ０ ３９７Ｆ ０ １８７Ｐ ０ １０６Ｌａ ０ ００２ＭＵＮ ０ １４６ Ｒ ２ ０ ９５８ ０ Ｐ ０ ０１ 参考文献 １ 乳糖的系数均与邹彩霞等 ２３ 乳能量预测方程中的 ２ ３ 乳脂肪 乳蛋白和乳糖的系数差异较大 其原因还 需进一步分析 预 测 方 程 式 的 残 差 范 围 分 别 为 １ １１３ ６ １ ４６６ ０ １ １２５ ０ １ ５７４ １ ０ ９９８ １ １ ５６０ ９ 和 １ ００７ ０ １ ５４７ ２ 以乳脂肪 乳蛋白 乳糖 乳脂 肪 乳蛋白 乳 糖 乳 尿 素 氮 含 量 为 预 测 因 子 的 三 ＭＵＳＧＲＡＶＥ Ｓ Ｄ ＳＡＬＩＳＢＵＲＹ Ｇ Ｗ Ｔｈｅ ｒｅｌａｔｉｏｎ ｓｈｉｐ ｏｆ ｍｉｌｋ ｅｎｅｒｇｙ ａｎｄ ｔｏｔａｌ ｐｒｏｔｅｉｎ ｔｏ ｐｅｒ ｃｅｎｔ ｆａｔ ｉｎ ｂｒｏｗｎ Ｓｗｉｓｓ ｈｅｒｄ ｍｉｌｋ Ｊ Ｊｏｕｒｎａｌ ｏｆ Ｄａｉｒｙ Ｓｃｉｅｎｃｅ 差异 这和乳能量受乳脂肪 乳蛋白和乳糖含量共 同影响有 关 ５ 而 其 方 程 中 的 乳 脂 肪 乳 蛋 白 和 可通过水牛 乳 中 的 乳 脂 肪 Ｆ 乳 蛋 白 Ｐ 乳糖 Ｌａ 及乳尿素氮 ＭＵＮ 含量预测 水 牛 乳 产 奶净能 Ｅ 方程分别为 Ｅ ０ ３８８Ｆ １ ５４０ Ｒ ２ 肪 乳蛋白和乳糖的系数相差不大 较接近于所报 道的乳脂肪的燃烧值 ０ ３８８ ＭＪ ｋｇ 而与乳蛋白及 论 ４ １９５２ ３５ ２ １７４ １７８ ＴＹＲＲＥＬＬ Ｈ Ｆ ＲＥＩＤ Ｊ Ｔ Ｐｒｅｄｉｃｔｉｏｎ ｏｆ ｔｈｅ ｅｎｅｒｇｙ ｖａｌｕｅ ｏｆ ｃｏｗ ｓ ｍｉｌｋ Ｊ Ｊｏｕｒｎａｌ ｏｆ Ｄａｉｒｙ Ｓｃｉｅｎｃｅ １９６５ ４８ ９ １２１５ １２２３ 中国奶牛饲养标准科研协作组 中国黑白花奶牛乳 中各种营养成分及能量间相关性的研究 Ｊ 中 国 畜牧杂志 １９８７ ２ ５ ７ 中华人民共和国农业部 ＮＹ Ｔ ３４ ２００４ 奶牛饲养 标准 Ｓ 北京 中国农业出版社 ２００４ ５ 美 国家科学研究 委 员 会 奶 牛 营 养 需 要 Ｍ 孟 ６ ＡＬＤＥＲＭＡＮ Ｇ Ｅｎｅｒｇｙ ａｎｄ ｐｒｏｔｅｉｎ ｒｅｑｕｉｒｅｍｅｎｔｓ ｏｆ 庆翔 译 北京 中国农业大学出版社 ２００２

8 2882 动物营养学报 29 卷 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 Sci ence,2007,90( 12) : [ 8 ] BRUNLAFLEUR L,DELABY L,HUSSON F.Predic ting energy protein interaction on milk yield and milk composition in dairy cows [ J]. Journal of Dairy Sci ence,2010,93( 9) : [ 9 ] 联合国粮农组织统计数据库 [ EB / OL].( ) [ ]. http: / / faostat3. fao. org / com pare / Q / QA / E. [10] 简保权, 秦学敏, 龚芳. 世界水牛奶业发展现状和典型模式分析 [ J]. 世界农业,2015(3) : [ 11] NASR M A. The impact of crossbreeding Egyptian and Italian buffalo on milk yield and composition un der subtropical environmental conditions [ J]. Journal of Dairy Research,2016,83( 2) : [ 12] JENKINS T C, MCGUIRE M A. Major advances in nutrition:impact on milk composition [ J]. Journal of Dairy Science,2006,89( 4) : [13] 塔娜, 吉日木图. 日粮中添加保护性脂肪对奶牛产奶量及乳成分的影响 [ J]. 中国奶牛,2008(9) : [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 buffa lo in China[ J]. Pakistan Journal of Zoology,2015,47 (5) : [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) : [16] 谢红, 刘琴, 周亚平. 云南省奶水牛乳成分与产奶量之间的关系 [ J]. 中国奶牛,2011(20) : [ 17] SUN Q, LV J P, LIU L, et al. Comparison of milk samples collected from some buffalo breeds and cross breeds in China [ J ]. Dairy Science & Technology, 2014,94( 4) : [ 18] RAFIQ S, HUMA N, PASHA I,et al.chemical com position,nitrogen fractions and amino acids profile of milk from different animal species[ J].Asian Australa sian Journal of Animal Sciences,2016,29( 7) : [19] ISLAM M A,ALAM M K,ISLAM M N,et al.princi pal milk components in buffalo,holstein cross,indige nous cattle and red chittagong cattle from Bangladesh [ J]. Asian Australasian Journal of Animal Sciences, 2014,27( 6) : [20] JAVED K,ABDULLAH M,KHALID M S, et al.in ter relationship of milk constituents with body and ud der measurements in NiLi Ravi buffaloes raised at commercial farms of Pakistan [ J]. Buffalo Bulletin, 2013,32( 2) : [21] 刘洪生. 影响乳成分的日粮调控技术 [ J]. 中国奶牛, 2008( 3) : [22] 王雅春, 杨雪丽, 张淑娟. 奶牛放乳前中后期乳成分变化的研究 [ J]. 中国奶牛,1997(4) : [23] 邹彩霞, 杨炳壮, 韦升菊, 等. 水牛乳能值与乳成分的相关性研究 [ J]. 动物营养学报,2010,22( 4) : [ 24] CLAEYS W L, VERRAES C, CARDOEN S, et al. Consumption of raw or heated milk from different spe cies: an evaluation of the nutritional and potential health benefits[ J].Food Control,2014,42: [ 25] ZOTOS A,BAMPIDIS V A.Milk fat quality of Greek buffalo ( Bubalus bubalis) [ J]. Journal of Food Com position and Analysis,2014,33( 2) : [ 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) : [27] 任春燕, 周凌云, 卜登攀, 等. 牛奶尿素氮的主要影响因素及其在奶牛生产中的应用 [ J]. 中国畜牧兽医, 2012,39( 9) : [28] FOODA T A,MOURAD K A,MOURAD G. Pheno typic and genetic trends for milk production in Egyp tian buffaloes[ J].Journal of American Science,2010, 6( 11) : [29] ARA U JO K B S,RANGEL A H N,FONSECA F C E,et al. Influence of the year and calving season on production, composition and mozzarella cheese yield of water buffalo in the State of Rio Grande Do Norte, Brazil[ J].Italian Journal of Animal Science,2012,11 (1) : [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) : [31] 冯岩松.SPSS 22.0 统计分析应用 [ M]. 北京 : 清华大学出版社,2015:6.

9 8 期李清等 : 云南水牛乳成分分析及乳能量预测模型的建立 2883 Establishment of Prediction Model of Milk Energy and Milk Composition Analysis of Yunnan Buffalo LI Qing 1,2 MAO Huaming 1,2 LI Wen 1,3 LI Jianhua 4 LU Qiongfen 1 HE Hongyuan 1 LIU Qin 5 ( 1. College of Animal Science and Technology, Yunnan Agricultural University, Kunming , China; 2. Yunnan Provincial Key Laboratory of Animal Nutrition and Feed Science, Kunming , China; 3. College of Life Science and Chemistry, Hunan University of Technology, Zhuzhou , China; 4. Animal Husbandry and Veterinary Station of Anning, Kunming , China; 5. Dairy Production Performance Measurement Center of Kunming, Kunming , China) Abstract: 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 in Yunnan province at 2014 to 2015 years 9 to 12 months, which were sent to Kunming Milk Production Center. The milk composi tion 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 tempera ture drying in vacuum drying oven. The results showed that milk fat, milk protein, milk lactose, milk urea ni trogen, milk total solids and milk energy were 6. 45%, 4. 55%, 5. 31%, mg / dl, % and 4.02 MJ / kg, respectively. And there was significant positive correlation between milk energy and milk fat, milk protein ( r = , r = , 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 ( F and P), milk fat, milk protein and milk lactose ( F, P and La), milk fat, milk protein, milk lactose and milk urea nitrogen contents ( F, P, La and MUN) as predictors, and regression equations were as followed: E = 0.388F ( R 2 = , P<0.01) ; E = F P ( R 2 = , P < 0. 01) ; E = F P La ( R 2 = , P<0.01) ; E = 0.397F+0.187P+0.106La+0.002MUN ( R 2 = , P<0.01). In conclusion, the milk production net energy can be predicted by the contents of milk fat, milk pro tein, milk lactose and milk urea nitrogen of buffalo milk. [ Chinese Journal of Animal Nutrition, 2017, 29 (8) : ] Key words: buffalo milk; milk composition; milk energy; prediction model Corresponding authors: MAO Huaming, professor, E mail: maohm@ vip.sina.com; LI Wen, associate professor, E mail: @ qq.com ( 责任编辑武海龙 )