不同羽色及性別烏蒙鳳雞育成期生長(zhǎng)曲線擬合研究

王天松 ，張 勇*，莫先艇，楊 紅，黃明捷，白 優(yōu)，陸 靜，郭 勇

（1.貴州大學(xué)高原山地動(dòng)物遺傳育種與繁殖教育部重點(diǎn)實(shí)驗(yàn)室，貴州大學(xué)動(dòng)物科學(xué)學(xué)院，貴州省動(dòng)物遺傳育種與繁殖重點(diǎn)實(shí)驗(yàn)室，貴州貴陽(yáng) 550025；2.畢節(jié)市畜牧遺傳資源管理站，貴州畢節(jié) 551700）

動(dòng)物生長(zhǎng)取決遺傳基因、飼糧營(yíng)養(yǎng)成分、環(huán)境等多種因素，建立數(shù)學(xué)生長(zhǎng)模型的目的是通過(guò)生長(zhǎng)曲線對(duì)生長(zhǎng)數(shù)據(jù)進(jìn)行簡(jiǎn)化描述，家禽常見(jiàn)的非線性生長(zhǎng)曲線模型有Brody[1]、EVF[2]、Gompertz[3]、Logistic[4]、Richards[5]、von Bertalanffy[1]和Weibull[6]等。生長(zhǎng)曲線模型能描述畜禽的生長(zhǎng)發(fā)育規(guī)律，其中的參數(shù)估計(jì)值也能反映不同種群在最大體重及脛長(zhǎng)、生長(zhǎng)速率等方面的差異。

烏蒙鳳雞原產(chǎn)于貴州省六盤(pán)水市水城縣，因頭部有明顯鳳頭特征，故命名為“烏蒙鳳雞”，該雞中等體型、體態(tài)勻稱、耐粗飼、抗逆性強(qiáng)、基因純合、遺傳性能穩(wěn)定，屬于肉蛋兼用型雞種。目前有關(guān)不同羽色及性別烏蒙鳳雞育成期生長(zhǎng)曲線預(yù)測(cè)模型尚未有研究報(bào)道。本研究用Logistic、Gompertz、von Bertalanffy 3 種非線性生長(zhǎng)曲線數(shù)學(xué)模型對(duì)育成期不同羽色及性別烏蒙鳳雞體重和脛長(zhǎng)進(jìn)行生長(zhǎng)曲線擬合，制定非線性曲線回歸模型，并用最佳生長(zhǎng)曲線描述不同羽色及性別烏蒙鳳雞的生長(zhǎng)特點(diǎn)及發(fā)育規(guī)律，為烏蒙鳳雞后期的科學(xué)飼養(yǎng)、選種選育提供數(shù)據(jù)參考。

1 材料與方法

1.1 實(shí)驗(yàn)動(dòng)物及方法 從六盤(pán)水市金鳳棲自然牧業(yè)有限公司烏蒙鳳雞保種擴(kuò)繁場(chǎng)選取同批次出殼體型適中、健康無(wú)疾病且有鳳冠的黃羽和麻羽烏蒙鳳雞公母各120 只（共480 只），同一飼養(yǎng)管理?xiàng)l件下科學(xué)飼養(yǎng)。參照《家禽生產(chǎn)性能名詞術(shù)語(yǔ)和度量統(tǒng)計(jì)方法》（NY/T 823-004）測(cè)定7～18 周齡的體重和脛長(zhǎng)。

1.2 擬合數(shù)據(jù)模型 利用Logistic、Gompertz 和von Bertalanffy 3 種生長(zhǎng)曲線模型分析不同羽色及性別烏蒙鳳雞育成期體重和脛長(zhǎng)發(fā)育特征，3 種模型的表達(dá)式及參數(shù)見(jiàn)表1，其中A 為最大體重（脛長(zhǎng)），B 為參數(shù)，k 為瞬時(shí)生長(zhǎng)速度，t 為周齡。

表1 3 種生長(zhǎng)曲線模型

1.3 統(tǒng)計(jì)分析 使用Excel 2010 和 SPSS 22.0 對(duì)體重和脛長(zhǎng)數(shù)據(jù)進(jìn)行分析、擬合，對(duì)函數(shù)參數(shù)A、B 和k 進(jìn)行估計(jì)，建立生長(zhǎng)曲線模型，根據(jù)擬合度（R2）評(píng)價(jià)生長(zhǎng)曲線模型。根據(jù)模型參數(shù)估計(jì)值計(jì)算7～18 周齡體重和脛長(zhǎng)模型預(yù)測(cè)值、拐點(diǎn)周齡、拐點(diǎn)體重和拐點(diǎn)脛長(zhǎng)。

2 結(jié)果與分析

2.1 不同羽色烏蒙鳳雞育成期生長(zhǎng)發(fā)育規(guī)律 由表2、3 可知，烏蒙鳳雞麻羽和黃羽公雞在12 周齡時(shí)體重與脛長(zhǎng)增速最快，麻羽公雞育成期的平均周增重大于黃羽公雞，但黃羽公雞的脛長(zhǎng)平均周增長(zhǎng)大于麻羽公雞，12周齡后平均周增重及脛長(zhǎng)增長(zhǎng)減少。黃羽母雞和麻羽母雞在10 周齡時(shí)平均周增重最大，為149.88 g 和163.32 g，11 周齡后平均周增重緩慢。麻羽母雞、黃羽母雞脛長(zhǎng)在7～11 周齡增長(zhǎng)較快，平均周增長(zhǎng)0.37 cm。由圖1、2 可知，烏蒙鳳雞黃羽和麻羽公母雞體重和脛長(zhǎng)生長(zhǎng)曲線基本重合，呈現(xiàn)出公雞7～11 周齡體重增長(zhǎng)緩慢，12～15 周齡增長(zhǎng)較快，12 周時(shí)達(dá)到增重最大；脛長(zhǎng)0～12 周齡增長(zhǎng)最快，12 周齡時(shí)達(dá)到增長(zhǎng)頂峰。母雞體重7～9 周齡增長(zhǎng)緩慢，10～15 周齡增長(zhǎng)較快，10 周齡增重最大；脛長(zhǎng)7～12 周齡增長(zhǎng)較快，11～18 周齡平緩增長(zhǎng)。

2.2 3 種生長(zhǎng)曲線模型的擬合分析 由表4、5 可以看出，3 種模型均能很好地?cái)M合烏蒙鳳雞黃羽公雞、麻羽公雞、黃羽母雞和麻羽母雞育成期體重和脛長(zhǎng)的生長(zhǎng)曲線，其R2都達(dá)到0.9 以上，3 種模型對(duì)體重和脛長(zhǎng)的擬合曲線與實(shí)際曲線基本一致。

由表4～7 可知，Logistic 模型對(duì)黃羽公雞、麻羽公雞、黃羽母雞和麻羽母雞體重和脛長(zhǎng)的R2大于Gompertz 模型和von Bertalanffy 模型，R2為L(zhǎng)ogistic>Gompertz>von Bertalanffy。Logistic 模型對(duì)體重的R2分別為0.980、0.980、0.990、0.986，對(duì)脛長(zhǎng)的R2分別為0.973、0.981、0.976、0.988。黃羽公雞體重拐點(diǎn)周齡為10.839 周，拐點(diǎn)體重為823.531 g，脛長(zhǎng)拐點(diǎn)周齡為6.82 周，拐點(diǎn)脛長(zhǎng)為5.62 cm；麻羽公雞體重拐點(diǎn)周齡為10.695 周，拐點(diǎn)體重為802.678 g，脛長(zhǎng)拐點(diǎn)周齡為6.31 周，拐點(diǎn)脛長(zhǎng)為5.61 cm；黃羽母雞體重拐點(diǎn)周齡為10.134 周，拐點(diǎn)體重為590.079 g，脛長(zhǎng)拐點(diǎn)周齡為4.73 周，拐點(diǎn)脛長(zhǎng)為4.05 cm；麻羽母雞拐點(diǎn)周齡為9.786 周，拐點(diǎn)體重為570.913 g，脛長(zhǎng)拐點(diǎn)周齡為4.53 周，拐點(diǎn)脛長(zhǎng)為4.26 cm。Logistic 模型對(duì)黃羽公雞、麻羽公雞、黃羽母雞和麻羽母雞育成期體重和脛長(zhǎng)估測(cè)值更接近于體重觀測(cè)值，且擬合效果和精度優(yōu)于Gompertz 模型和von Bertalanffy 模型，所以Logistic 模型能更好地?cái)M合不同羽色和性別烏蒙鳳雞的體重和脛長(zhǎng)生長(zhǎng)規(guī)律。

表2 不同羽色及性別烏蒙鳳雞體重和平均周增重 g

表3 不同羽色及性別烏蒙鳳雞脛長(zhǎng)和平均周增長(zhǎng) cm

圖1 烏蒙鳳雞育成期體重生長(zhǎng)曲線

圖2 烏蒙鳳雞育成期脛長(zhǎng)生長(zhǎng)曲線

3 討 論

3.1 生長(zhǎng)發(fā)育規(guī)律分析 體尺性狀是反映動(dòng)物生長(zhǎng)發(fā)育快慢的重要指標(biāo)，受遺傳與環(huán)境等多種因素影響，而羽色作為禽類選種選育過(guò)程中的重要參考指標(biāo)，對(duì)禽類生產(chǎn)及其市場(chǎng)需求具有很大影響。本研究發(fā)現(xiàn)育成期烏蒙鳳雞黃羽公雞與麻羽公雞的體重和脛長(zhǎng)在12 周齡時(shí)增長(zhǎng)最快；黃羽母雞和麻羽母雞體重在10 周齡時(shí)增長(zhǎng)最快，脛長(zhǎng)分別在9 周齡及11 周齡增長(zhǎng)最快。因此，在烏蒙鳳雞育成期相應(yīng)周齡，應(yīng)提供足夠的能量、蛋白質(zhì)、鈣等營(yíng)養(yǎng)物質(zhì)以維持烏蒙鳳雞體重、脛長(zhǎng)指標(biāo)的生長(zhǎng)發(fā)育需要。本研究結(jié)果顯示黃羽和麻羽公母雞體重和脛長(zhǎng)生長(zhǎng)曲線基本吻合，說(shuō)明不同羽色烏蒙鳳雞生長(zhǎng)發(fā)育不存在差異，與黃炎坤[7]、趙書(shū)雨[8]、王曉薇[9]等研究中關(guān)于羽色對(duì)雞群的生長(zhǎng)發(fā)育無(wú)明顯影響的結(jié)果一致。本實(shí)驗(yàn)中不同羽色及性別烏蒙鳳雞體重的變化與李巖等[10]和李義書(shū)等[11]研究結(jié)果一致。

3.2 體重和脛長(zhǎng)生長(zhǎng)曲線擬合結(jié)果分析 理想的生長(zhǎng)曲線模型可描述畜禽的生長(zhǎng)發(fā)育情況，確定其生長(zhǎng)發(fā)育基本模式，同時(shí)預(yù)測(cè)其實(shí)際生產(chǎn)條件下所需要的營(yíng)養(yǎng)并采取相應(yīng)的飼養(yǎng)方案，對(duì)其飼養(yǎng)管理具有十分重要的指導(dǎo)意義。本研究中Logistic、Gompertz 和von Bertalanffy 3 種生長(zhǎng)曲線模均能很好地?cái)M合不同羽色和不同性別烏蒙鳳雞育成期體重和脛長(zhǎng)生長(zhǎng)曲線，R2都達(dá)到0.9 以上，3 種模型對(duì)黃羽公雞、麻羽公雞、黃羽母雞和麻羽母雞體重和脛長(zhǎng)的擬合精度為L(zhǎng)ogistic>Gompertz>von Bertalanffy，Logistic 模型對(duì)體重和脛長(zhǎng)估測(cè)值更接近于體重觀測(cè)值。白優(yōu)等[12]對(duì)烏蒙鳳雞0、4、8、12、16、20、24 周齡的體重和脛長(zhǎng)進(jìn)行擬合分析，得出Gompertz 模型對(duì)體重生長(zhǎng)規(guī)律擬合效果更好，與本研究Logistic 模型對(duì)體重?cái)M合最優(yōu)結(jié)果不同，可能與2 個(gè)研究的數(shù)據(jù)統(tǒng)計(jì)間隔時(shí)間存在差異有關(guān)，但2 個(gè)研究在脛長(zhǎng)的擬合分析結(jié)果上是一致的。劉小輝等[13]對(duì)壩上長(zhǎng)尾雞育雛育成期生長(zhǎng)曲線進(jìn)行擬合研究，得到公雞拐點(diǎn)周齡（9.859 周）、拐點(diǎn)體重（741.470 g），母雞拐點(diǎn)周齡（8.696 周）、拐點(diǎn)體重（513.608 g），而本研究中烏蒙鳳雞黃羽公雞體重拐點(diǎn)周齡為10.839 周，拐點(diǎn)體重為823.531 g，麻羽公雞體重拐點(diǎn)周齡為10.695 周，拐點(diǎn)體重為802.678 g，黃羽母雞體重拐點(diǎn)周齡為10.134周，拐點(diǎn)體重為590.079 g，麻羽母雞拐點(diǎn)周齡為9.786 周，拐點(diǎn)體重為570.913 g，均大于壩上長(zhǎng)尾雞[13]，說(shuō)明雞種、飼養(yǎng)管理、飼養(yǎng)環(huán)境等條件不同，擬合結(jié)果存在差異。李洪林等[14]研究表明Logistic 模型對(duì)貴州黃雞慢羽系體重?cái)M合效果最佳，與本研究結(jié)果相同。宋丹等[15]對(duì)京紅蛋雞、劉聰?shù)萚16]對(duì)麻雞、鄒劍敏等[17]對(duì)清遠(yuǎn)麻雞脛長(zhǎng)擬合研究結(jié)果均表明Logistic 模型的擬合效果最好，與本研究結(jié)果一致，上述研究中脛長(zhǎng)的拐點(diǎn)周齡、拐點(diǎn)脛長(zhǎng)值與本研究存在差異。

表4 不同羽色及性別烏蒙鳳雞育成期體重參數(shù)估計(jì)值和R2

表5 不同羽色及性別烏蒙鳳雞育成期脛長(zhǎng)參數(shù)估計(jì)值和R2

4 結(jié) 論

本研究運(yùn)用Logistic、Gompertz 和von Bertalanffy 3 種生長(zhǎng)曲線模對(duì)育成期不同羽色及性別烏蒙鳳雞的體重及脛長(zhǎng)進(jìn)行了生長(zhǎng)曲線擬合和分析，發(fā)現(xiàn)3 種模型的擬合效果均較好，但以Logistic 模型的擬合效果最佳，且公雞擬合效果比母雞好。

表6 不同羽色及性別烏蒙鳳雞育成期體重觀測(cè)值和擬合曲線估計(jì)值

表7 不同羽色及性別烏蒙鳳雞育成期脛長(zhǎng)觀測(cè)值和擬合曲線估計(jì)值