吉林矮小蘆花雞早期生長曲線擬合及相關(guān)分析

柳儉強(qiáng)，劉臣，張芳毓，武斌，曹陽，肖成，劉艷光,2，劉洪亮，張立春*

(1.吉林省農(nóng)業(yè)科學(xué)院畜牧分院,吉林 公主嶺 136100；2.延邊大學(xué)農(nóng)學(xué)院,吉林 延吉 133002)

吉林地方品種以蘆花雞較為常見，通過多世代選育，培育出個體小、基礎(chǔ)代謝低、飼料報酬高、抗病抗應(yīng)激能力強(qiáng)的肉蛋兼用型品系吉林矮小蘆花雞，具有良好的市場前景。

家禽生長發(fā)育規(guī)律呈非線性，通過非線性模型進(jìn)行擬合分析，可較為動態(tài)了解其生長發(fā)育趨勢與規(guī)律，指導(dǎo)家禽的選育和飼養(yǎng)管理，使其生產(chǎn)潛力發(fā)揮到最大[1]。目前最為常用的是Logistic、Gompertz 和 Bertalanffy 3種非線性曲線模型，已經(jīng)廣泛用于家禽生長發(fā)育規(guī)律研究中[2-3]。本試驗通過體重和脛長指標(biāo)，并結(jié)合3種曲線模型對吉林矮小蘆花雞生長發(fā)育規(guī)律進(jìn)行擬合分析，旨在為其早期的選育提供科學(xué)依據(jù)。

1 材料與方法

1.1 試驗材料

選擇吉林省農(nóng)業(yè)科學(xué)院畜牧分院試驗雞場的吉林蘆花雞在人工選擇的基礎(chǔ)上，通過純繁獲得吉林矮小蘆花雞(dw)和吉林蘆花雞(Dw)各200只，公母各100只，共400只。

1.2 飼養(yǎng)管理

確定分組后將每只試驗雞穿翅號，自由采食和飲水，日常免疫接種與管理。飼養(yǎng)過程分為1～5周的育雛期和6～13周的育成期。

1.3 測定指標(biāo)

1.3.1 體重

用電子天平空腹稱量試驗雞的初生重以及0、1、3、5、7、9、11、13周齡的體重，并做好記錄。

1.3.2 脛長和脛骨長

用游標(biāo)卡尺測量試驗雞7、9、11、13周齡的脛長和脛骨長(在測量脛骨長時，若長度超出游標(biāo)卡尺最大范圍，可用皮尺替換)，并做好記錄。

1.4 生長曲線擬合

選擇Logistic、Gompertz 和 Bertalanffy 3種非線性曲線模型對0～13周生長曲線進(jìn)行擬合，具體表達(dá)式與參數(shù)見表1。

表1 3種非線性模型曲線

1.5 數(shù)據(jù)統(tǒng)計與分析

利用Excel進(jìn)行初步數(shù)據(jù)整理，結(jié)合SPSS 24.0單因素方差分析和Duncan’s多重比較，數(shù)值用“平均數(shù)±標(biāo)準(zhǔn)誤”表示，擬合3種非線性曲線，分別計算A、B、K的最優(yōu)估計值，建模后再計算出拐點體重、周齡以及最大周增重。然后對體重與脛長和脛骨長進(jìn)行相關(guān)性分析。

2 結(jié)果與分析

2.1 體重實際生長情況

由表2和圖1可知，dw母雞初生重與dw公雞、Dw母雞差異不顯著(P>0.05)，與Dw公雞差異顯著(P<0.05)；dw母雞從3周開始生長速度明顯加快，但與另外3組相比較體重保持最低，差異顯著(P<0.05)，在第7周時達(dá)到高峰，隨后生長速度放緩，與Dw母雞相比，從5周開始，生長趨勢大致保持一致。dw公雞5～9周增長迅速，并在9周時達(dá)到高峰，比Dw公雞高峰期提前了2周左右，隨后放緩在11周后又略有上升。在整個早期生長階段，dw公雞始終低于Dw公雞，差異顯著(P<0.05)。早期同一生長階段體重大致為：Dw公雞>Dw母雞≈dw公雞>dw母雞。

表2 0～13周齡體重測定值比較 g

圖1 0～13周齡吉林矮小蘆花雞周增重

2.2 生長曲線擬合分析

由表3可知，所有雞種生長曲線模型擬合度均超過了0.99，其中dw母雞Gompertz和Bertalanffy模型擬合度達(dá)到0.995，Gompertz模型拐點體重為607.924 g，拐點周齡為7.345，最大周增重達(dá)到113.682 g；dw公雞Logistic和Gompertz模型擬合度相同；Dw母雞Gompertz模型擬合度最高，達(dá)到0.995，拐點體重為808.599 g，拐點周齡7.342，最大周增重為164.146 g；Dw公雞Logistic和Gompertz模型擬合度均為0.995，拐點體重、周齡、最大周增重差異不明顯，均可作為適合模型，另外，Dw、dw公雞的拐點周齡均晚于Dw、dw母雞。

表3 3種生長曲線模型的擬合參數(shù)

2.3 生長曲線與3種擬合曲線的比較

如表4、表5所示，通過比較dw和Dw雞實際測定值和3種生長曲線模型的估計值發(fā)現(xiàn)，Logistic模型的初生重估計值遠(yuǎn)高于實際測定值，Bertalanffy模型的初生重估計值遠(yuǎn)低于實際測定值，而Gompertz模型的初生重估計值與實際測定值較為接近，整個早期生長階段的模擬值與實際測定值吻合度較高，因此Gompertz模型擬合效果優(yōu)于Logistic和Bertalanffy模型。

表4 dw雞實際測定值與3種擬合曲線的比較

表5 Dw雞實際測定值與3種擬合曲線的比較

2.4 脛長和脛骨長的發(fā)育情況

各組雞脛長和脛骨長發(fā)育情況見表6、表7。同一階段dw母雞脛長均顯著低于Dw母雞(P<0.05)，隨著周齡增加，dw母雞脛長分別占Dw母雞的82.7%、82.4%、75.5%和72.2%，差距逐漸變大，dw母雞在7～9周脛長增速明顯，9周之后增長速度放緩，Dw母雞在7～11周增速逐漸加快，11周以后增速放緩；同一階段dw母雞脛骨長均顯著低于Dw母雞(P<0.05)，在觀測期中，dw母雞脛骨長分別占Dw母雞的85.2%、80.0%、81.1%和68.4%，差距逐漸拉大，dw母雞脛骨長在7～11周平穩(wěn)增加，11周以后趨于平穩(wěn)，而Dw母雞脛骨長在觀測期一直逐漸增加；dw母雞脛長與脛骨長生長趨勢幾乎一致，Dw母雞脛長與脛骨長生長趨勢大致相同。

表6 7～13周齡脛長比較 mm

表7 7～13周齡脛骨長比較 mm

同一階段dw與Dw公雞脛長相比較差異顯著(P<0.05)，dw公雞脛長先增加在9周后逐漸穩(wěn)定，dw公雞脛長分別占Dw公雞的86.3%、84.9%、73.9%和73.7%，而Dw公雞脛長呈現(xiàn)7～11周逐漸增加，11周后增速放緩，脛長趨于穩(wěn)定，dw公雞脛骨長分別占Dw公雞的85.4%、81.9%、74.9%和73.0%；同一階段dw與Dw公雞脛骨長相比較差異顯著(P<0.05)，Dw公雞脛骨長逐漸增加，dw公雞脛骨的增長速度呈先增加后減少的趨勢，dw公雞脛長與脛骨長生長趨勢無明顯差別，Dw公雞脛長與脛骨長生長趨勢大致相同。

2.5 體重與脛長、脛骨長的相關(guān)性分析

體重與脛長和脛骨長的結(jié)果見表8和表9。不同周齡各組合的脛長與體重呈現(xiàn)出極顯著正相關(guān)(P<0.01)，其中dw與Dw母雞分別在7周齡和9周齡相關(guān)系數(shù)最大，為0.934和0.891；dw和Dw公雞均在9周齡相關(guān)系數(shù)最大；dw雞脛長與體重的相關(guān)程度高于正常的Dw雞。不同周齡各組合的脛骨長與體重除dw母雞13周齡呈現(xiàn)出顯著正相關(guān)(P<0.05)外，其余均呈現(xiàn)出極顯著正相關(guān)(P<0.01)，其中dw和Dw母雞脛骨長與體重的相關(guān)程度分別在7、9周齡時最大，而dw和Dw公雞則表現(xiàn)在9和11周齡。

表8 不同周齡脛長與體重相關(guān)性分析

表9 不同周齡脛骨長與體重相關(guān)性分析

脛長與脛骨長的相關(guān)性結(jié)果如表10所示，7～9周齡各組合脛長與脛骨長的相關(guān)性達(dá)到極顯著水平(P<0.01)；11周齡dw母雞和公雞脛長與脛骨長的相關(guān)性達(dá)到極顯著水平(P<0.01)，Dw母雞和Dw公雞脛長與脛骨長的相關(guān)性達(dá)到顯著水平(P<0.05)；13周齡各組合脛長與脛骨長的相關(guān)性達(dá)到顯著水平(P<0.05)。

表10 不同周齡脛長與脛骨長相關(guān)性分析

3 討論

3.1 生長曲線擬合特點

動物生長發(fā)育具有相對穩(wěn)定的基本特征,這是用數(shù)學(xué)模型描述動物體生長過程的基礎(chǔ)[4]。非線性生長模型是根據(jù)生物體生物量增加，描述其隨時間規(guī)律性的變化的過程[5]。動物生長發(fā)育在遺傳、環(huán)境及營養(yǎng)等多種因素的作用下，不同品種間差異大，因此不同品種要用不同的生長模型去擬合。目前常用的3種非線性生長模型是Logistic、Gompertz和Bertalanffy模型，Logistic模型適用于早期生長速度慢的品種，Gompertz則適用于早期生長速度快的品種，Bertalanffy適用于生長較為遲緩的動物[6]。本研究利用這3個非線性曲線模型對0～13周的吉林矮小蘆花雞和正常蘆花雞進(jìn)行了擬合分析。結(jié)果表明Gompertz模型擬合度(R2)最高達(dá)到0.995，擬合值與實際測定值最為接近，效果最佳，這一結(jié)果與Mohammed[7]、李青青等[8]、杜德英等[9]、付睿琦等[10]對羅斯308、獨龍雞、文昌雞和北京油雞的研究結(jié)果相一致。Manjula等[11]和包艷波等[12]則認(rèn)為Logistic模型更適用于韓國本地雞(KNC)和大骨雞的生長曲線擬合。姜自琴等[13]認(rèn)為，Bertallanffy 模型在擬合度和預(yù)測武定雞體重方面效果最好。吉林矮小蘆花雞和正常吉林蘆花雞中，公雞的生長周期始終要比母雞長，拐點周齡相對較晚,所以在飼養(yǎng)時要注意分群，并在此之前提供充足的養(yǎng)分，來滿足其生長發(fā)育的需求。因此，Gompertz模型對吉林矮小蘆花雞的擬合效果最佳，拐點周齡及體重與實測值最為接近。

3.2 脛長、脛骨長和體重的相關(guān)性分析

骨骼作為家禽機(jī)體的重要組成，在其生長發(fā)育過程中有著重要的意義，對于生產(chǎn)潛力挖掘產(chǎn)生了重要影響[14]。脛長和脛骨長則是衡量家禽早期生長發(fā)育的重要指標(biāo)，與體重結(jié)合進(jìn)行相關(guān)性分析，可為吉林矮小型蘆花雞的選育提供有價值的參考。本研究結(jié)果表明，dw和Dw公母雞在各自周齡脛長、脛骨長與體重呈極顯著正相關(guān)，這與陳銳等[15]認(rèn)為略陽烏雞脛長與體質(zhì)量呈極顯著正相關(guān)的結(jié)果相一致。葉偉慶等[16]的研究結(jié)果表明，信宜懷鄉(xiāng)雞體質(zhì)量與脛長極顯著正相關(guān)。成海榮等[17]提出，脛骨長度和重量是反映肉雞骨骼發(fā)育的重要指標(biāo)。本研究發(fā)現(xiàn)，dw和Dw公母雞脛骨相對生長在9～11周齡時達(dá)到最高，這與寧康健等[18]發(fā)現(xiàn)青腳麻公母雞脛骨的相對生長均在28日齡達(dá)到最高結(jié)果不一致。這可能是由于品種、營養(yǎng)、環(huán)境等因素造成的脛骨發(fā)育的差異。吉林矮小蘆花雞各周齡體重與脛長和脛骨長呈極顯著正相關(guān)，脛長和脛骨長可作為反映生長發(fā)育的主要指標(biāo)。