■趙江濤 艾 琴，2 左建軍 馮定遠(yuǎn)* 譚會澤

（1.廣東溫氏食品集團(tuán)股份有限公司，廣東云浮 527439；2.廣東溫氏大華農(nóng)生物科技有限公司，廣東云浮 527439；3.華南農(nóng)業(yè)大學(xué)動物科學(xué)院，廣東廣州 510642）

從20世紀(jì)80年代飼用非淀粉多糖酶制劑引入中國養(yǎng)殖業(yè)以來，木聚糖酶的功效一直受到爭議，主要存在如下觀點(diǎn)：①木聚糖酶對生產(chǎn)性能無顯著影響；②木聚糖酶可提高肉雞日增重，降低日采食量和料重比；③飼糧中添加木聚糖酶后生產(chǎn)性能反而出現(xiàn)降低趨勢；④木聚糖酶在仔雞階段的作用效果優(yōu)于生長雞階段。為了尋求一個(gè)統(tǒng)一、穩(wěn)定的答案，國內(nèi)外眾多科研機(jī)構(gòu)或生產(chǎn)單位每年都投入大量資金和人力展開研究。主要的試驗(yàn)設(shè)計(jì)思路包括：①設(shè)計(jì)木聚糖酶的不同添加梯度；②設(shè)計(jì)低能日糧對照組；③采用不同的日糧結(jié)構(gòu)；④針對不同的飼養(yǎng)階段。由于酶活力的定義和測定標(biāo)準(zhǔn)以及原料營養(yǎng)水平和飼養(yǎng)標(biāo)準(zhǔn)不同，因此，酶的不同添加梯度和飼糧能量水平不能作為評價(jià)酶制劑的指標(biāo)。那么，日糧結(jié)構(gòu)和日齡是否會影響木聚糖酶對肉雞生產(chǎn)性能的作用效果呢？針對以上問題，本文采用Meta分析統(tǒng)計(jì)了2002～2012年間公開發(fā)表的數(shù)據(jù)，系統(tǒng)評價(jià)了木聚糖酶對肉雞生產(chǎn)性能指標(biāo)的影響，并通過Meta回歸對其影響因子進(jìn)行了分析，為木聚糖酶作用效果的進(jìn)一步研究提供理論和數(shù)據(jù)依據(jù)。

1 材料與方法

1.1 分析方法

本文采用Meta分析的方法對54篇相關(guān)的隨機(jī)對照試驗(yàn)進(jìn)行了綜合評價(jià)。Meta分析是一種對多個(gè)同類研究結(jié)果進(jìn)行合并匯總的分析方法，從統(tǒng)計(jì)的角度來講，該分析達(dá)到了增大樣本含量，提高檢驗(yàn)效能的目的，尤其是當(dāng)多個(gè)研究結(jié)果不一致或都沒有統(tǒng)計(jì)學(xué)意義時(shí)，采用Meta分析可得到更加接近真實(shí)情況的綜合分析結(jié)果。

處理方法包括：①異質(zhì)性檢驗(yàn)（齊性檢驗(yàn)）；②統(tǒng)計(jì)合并效應(yīng)量（加權(quán)合并，計(jì)算效應(yīng)尺度及95%的置信區(qū)間）并進(jìn)行統(tǒng)計(jì)推斷；③圖示單個(gè)試驗(yàn)的結(jié)果和合并后的結(jié)果；④敏感性分析；⑤通過“倒漏斗圖”了解潛在的發(fā)表偏倚。

1.2 文獻(xiàn)收集

相關(guān)文獻(xiàn)利用華南農(nóng)業(yè)大學(xué)圖書館（http://202.116.174.99/find/find_journal.html）SCI（ISTP,BP,MEDLIN）平臺和Google Scholar search engines搜索并由人工篩選整理。檢索條件為：“標(biāo)題”xylanase and broiler and performance and“出版年”2002～2012。

1.3 納入標(biāo)準(zhǔn)

選用的研究需滿足以下標(biāo)準(zhǔn)：①2002～2012年間發(fā)表的英文文獻(xiàn)；②研究對象為肉雞；③研究中包含內(nèi)切木聚糖酶單酶（endo-xylanase）空白組和試驗(yàn)組；④雞生產(chǎn)性能指標(biāo)包括1～21 d和22～42 d平均日增重（ADG）、平均日采食量（ADFI）和料重比（F∶G）；⑤研究中給出指標(biāo)的數(shù)值或平均值（mean）、標(biāo)準(zhǔn)差（SD）或標(biāo)準(zhǔn)誤（SE）。SD重在描述變量值的分布，SE重在反映抽樣誤差的大小，二者轉(zhuǎn)化公式為SE=SD/，n表示樣本量；⑥試驗(yàn)設(shè)計(jì)為完全隨機(jī)對照試驗(yàn)。

1.4 排除標(biāo)準(zhǔn)

①文獻(xiàn)摘要、會議論文和未發(fā)表的報(bào)道；②木聚糖酶與其他NSP或消化酶復(fù)合使用的研究；③數(shù)據(jù)不全或無法利用的文獻(xiàn)。

1.5 數(shù)據(jù)提取和整理

提取的數(shù)據(jù)包括每個(gè)處理的效應(yīng)因子平均值SD，試驗(yàn)動物個(gè)數(shù)的數(shù)據(jù)。分析的參數(shù)主要有平均日增重、平均日采食量、料重比。評價(jià)的精確性以文獻(xiàn)中報(bào)道的SD為依據(jù)，如果SD沒有提供，則利用SD=SE×計(jì)算得到。

1.6 統(tǒng)計(jì)分析

本研究利用 STATA（Intercooled Stata V.12.0，College Station,TX）軟件的固定效應(yīng)模型計(jì)算了每項(xiàng)研究的效應(yīng)值、權(quán)重、方差及綜合效應(yīng)量。采用標(biāo)準(zhǔn)差作為權(quán)重的衡量依據(jù)，標(biāo)準(zhǔn)差越小，權(quán)重越大；綜合效應(yīng)采用標(biāo)準(zhǔn)均數(shù)差（standardized mean difference，SMD）表示，SMD可簡單地理解為兩均數(shù)的差值再除以合并標(biāo)準(zhǔn)差的商，它不僅消除了多個(gè)研究間絕對值大小的影響，還消除了多個(gè)研究測量單位不同的影響。本次Meta分析報(bào)告了固定效應(yīng)模型的結(jié)果，用Q檢驗(yàn)（服從卡方分布）評價(jià)各研究之間的異質(zhì)性，用漏斗圖、Begg秩相關(guān)法、Egger回歸法檢驗(yàn)了發(fā)表偏倚，并對納入的研究結(jié)果進(jìn)行敏感性分析。

2 結(jié)果與分析

2.1 納入文獻(xiàn)的特征

數(shù)據(jù)庫初檢出文獻(xiàn)54篇相關(guān)的隨機(jī)對照試驗(yàn)，通過上述納入標(biāo)準(zhǔn)逐步篩選后，最終有9篇研究納入本次Meta分析（表1～表2）。被納入的文獻(xiàn)發(fā)表于2002～2012年間，互相之間沒有重復(fù)。這些研究觀察了添加木聚糖酶后，對照組和試驗(yàn)中平均日增重、平均日采食量和料重比的對比情況。個(gè)別研究中有設(shè)計(jì)木聚糖酶不同的添加梯度或能量梯度，但本次研究中只引入飼養(yǎng)階段和日糧結(jié)構(gòu)兩個(gè)因素，因此單篇研究可能有多個(gè)數(shù)據(jù)納入，其中平均日增重、平均日采食量和料重比的納入數(shù)據(jù)分別為28、27組和25組。

2.2 Meta分析結(jié)果

飼糧中添加木聚糖酶對肉雞生產(chǎn)性能影響的效應(yīng)因子Meta分析結(jié)果見表3。

表1 Meta分析肉雞生產(chǎn)性能數(shù)據(jù)

表2 Meta分析肉雞生產(chǎn)性能數(shù)據(jù)總匯

表3 Meta分析中木聚糖酶對生產(chǎn)性能影響的效應(yīng)因子表

由表3可見，通過對所有研究的分析，飼糧中添加木聚糖酶可使肉雞的平均日增重提高0.46（0.19～0.73）g/(d·只)[ES（效應(yīng)量）=0.20，P2=0.027，n（引用文獻(xiàn)中重復(fù)數(shù)）=244]，可使料重比降低0.02（ES=-0.277，P2=0.006，n=243），但對平均日采食量的影響不大（ES=0.115，P2=0.216，n=214）。通過I2值計(jì)算發(fā)現(xiàn)平均日增重、平均日采食量和料重比的異質(zhì)性均不顯著（分別為I2=0%、12.6%、14.1%；P1=0.808、0.276、0.262）。因此，沒有必要進(jìn)行異質(zhì)性檢驗(yàn)分析，并且所有模型均選用固定效應(yīng)模型計(jì)算。

ADG、ADFI和F∶G影響的Meta分析森林圖見圖1～圖3。圖中X軸表示SMD值（利用標(biāo)準(zhǔn)化的Z統(tǒng)計(jì)計(jì)算得到），水平線的長度代表每個(gè)研究的95%置信區(qū)間，中間方形面積大小代表該研究在所有研究中所占比重，具體數(shù)值列于右側(cè)。圖1中虛線上的菱形代表飼糧中添加木聚糖酶可顯著提高肉雞平均日增重（P2=0.027），所有研究異質(zhì)性指數(shù)的I2=0%、P1=0.808，說明研究異質(zhì)性不顯著。圖2中圖中虛線上的菱形代表飼糧中添加木聚糖酶可提高肉雞的平均日采食量，但提高的程度不顯著（P2=0.216），所有研究異質(zhì)性指數(shù)的I2=12.6%、P1=0.276，說明研究異質(zhì)性不顯著。圖3中虛線上的菱形代表添加木聚糖酶明顯降低肉雞料重比（P2=0.006），所有研究異質(zhì)性指數(shù)的I2=14.1%、P1=0.262，說明研究異質(zhì)性不顯著。

圖1 日糧中添加木聚糖酶對肉雞平均日增重影響的Meta分析森林圖

圖2 日糧中添加木聚糖酶對肉雞平均日采食量影響的Meta分析森林圖

圖3 日糧中添加木聚糖酶對肉雞料重比影響的Meta分析森林圖

2.3 發(fā)表偏倚的檢驗(yàn)

為了對納入的研究進(jìn)行發(fā)表偏倚的檢驗(yàn)，對納入文獻(xiàn)的研究數(shù)據(jù)進(jìn)行了Egger's檢驗(yàn)（見表4），結(jié)果發(fā)現(xiàn)平均日增重、平均日采食量和料重比的P值分別為0.002、＜0.001和＜0.001，即所有試驗(yàn)因子均存在發(fā)表偏倚。

表4 變量發(fā)表偏倚Egger's檢測

同時(shí)利用Begg's檢測其發(fā)表偏倚并以圖形（見圖4～圖6）形式直觀的顯示發(fā)表偏倚狀況，其中X軸代表SMD的SE，Y軸代表SMD，水平線代表全部的效應(yīng)因子預(yù)測值。兩條斜線代表效應(yīng)因子預(yù)測值的95%置信區(qū)間。水平線兩側(cè)的研究個(gè)數(shù)代表研究間的發(fā)表偏倚。由圖4～圖6可見，平均日增重和料重比均存在顯著的發(fā)表偏倚（P＜0.001，＜0.001），但平均日采食量的發(fā)表偏倚不顯著（P=0.05）。

圖4 飼糧中添加木聚糖酶對肉雞平均日增重影響的發(fā)表偏倚漏斗圖

2.4 敏感性分析

各因子的敏感性分析見圖7～圖9，圖中X軸代表SMD，Y軸代表單個(gè)研究。水平線的長度代表95%置信區(qū)間，其中圓圈代表去除左側(cè)對應(yīng)的研究后，剩余所有研究的SMD值。所有圓圈均處于所有研究的95%置信區(qū)間內(nèi)，表明研究結(jié)論非常穩(wěn)定。由圖7～圖9可見，本試驗(yàn)敏感性分析表明不存在敏感研究，去掉任何一個(gè)單獨(dú)研究都不會改變整體趨勢，對于全部研究結(jié)果影響不大。

圖5 飼糧中添加木聚糖酶對肉雞平均日采食量影響的發(fā)表偏倚漏斗圖

圖6 飼糧中添加木聚糖酶對肉雞料重比影響的發(fā)表偏倚漏斗圖

圖7 日糧中添加木聚糖酶對肉雞平均日增重影響的敏感性分析

3 小結(jié)

3.1 Meta分析9篇文獻(xiàn)28組研究數(shù)據(jù)顯示肉雞日糧中添加內(nèi)切木聚糖酶可使ADG提高0.46（0.19～0.73）g/(d·只)（P＜0.05），F(xiàn)∶G 降低 0.02（0.01～0.03）（P＜0.05），對ADFI影響不顯著（P ＞0.05）。

3.2 所有研究存在的顯著的發(fā)表偏倚（P＜0.05），添加木聚糖酶后提高日增重、降低采食量和料重比的研究發(fā)表數(shù)較多。

圖8 日糧中添加木聚糖酶對肉雞平均日采食量影響的敏感性分析

圖9 日糧中添加木聚糖酶對肉雞料重比影響的敏感性分析

（參考文獻(xiàn)刊略，需者可函索zhaojt186@sina.com）