邢 磊 袁紅艷 張春華 趙樂樂 王文濤 焦 楠 陸雪林*
1.上海市動(dòng)物疫病預(yù)防控制中心,上海201103;2.上海欣灝珍禽育種有限公司,上海201408
申鴻七彩雉在2019年4月通過國(guó)家畜禽遺傳資源委員會(huì)審定,是我國(guó)雉雞行業(yè)第1 個(gè)具有獨(dú)立知識(shí)產(chǎn)權(quán)的國(guó)家審定品種,該品種體重較大、體型豐滿、屠宰率和全凈膛率高、肉質(zhì)優(yōu)、產(chǎn)蛋多、蛋品質(zhì)好、主要生產(chǎn)性狀遺傳穩(wěn)定性好。畜禽生長(zhǎng)一般均具有非線性的特點(diǎn),可以用非線性數(shù)學(xué)模型來(lái)描述動(dòng)物的生長(zhǎng)規(guī)律[1]。理想的動(dòng)物生長(zhǎng)曲線模型,有助于判斷和分析飼養(yǎng)、管理、防疫上的成功與失敗,可以用來(lái)比較和檢驗(yàn)畜禽不同品種類型、親緣關(guān)系和性別的遺傳品質(zhì)等[2]。目前,在畜禽中常用的擬合模型包括Logistic、Gompertz 和Von Berta?lanffy 等[3]。柳儉強(qiáng)等[4]采用Logistic、Gompertz 和Von Bertalanffy 3 種非線性生長(zhǎng)模型對(duì)所選育的吉林矮小蘆花雞0~13 周齡體重指標(biāo)的變化進(jìn)行擬合分析;苗小猛等[5]采用Gompertz、Logistic 和Von Ber?talanffy 模型對(duì)0~10 周齡的群體體重進(jìn)行曲線擬合分析,對(duì)10周齡三穗鴨公母之間的體尺性狀進(jìn)行比較,并分析體重與各體尺性狀間的相關(guān)性;牟騰慧等[6]采用Logistic 模型、Gompertz 模型以及Von Berta?lanffy 模型對(duì)威寧雞0~16 周齡體重變化進(jìn)行擬合分析。前期研究主要集中在申鴻七彩雉肉質(zhì)性狀[7]和蛋品質(zhì)[8]方面,缺少對(duì)申鴻七彩雉生長(zhǎng)曲線的研究。因此,本研究通過對(duì)申鴻七彩雉生長(zhǎng)測(cè)定記錄的分析,采用以上3種最常用的生長(zhǎng)模型,擬合其生長(zhǎng)曲線,從而獲得其生長(zhǎng)發(fā)育規(guī)律,為申鴻七彩雉生產(chǎn)與育種工作提供理論依據(jù)。
采用上海欣灝珍禽育種有限公司飼養(yǎng)剛出雛的雉雞公母各1 000 只,隨機(jī)抽取公母各50 只測(cè)定體重,以后每2 周測(cè)定1 次體重,直至24 周齡,自由采食及飲水,執(zhí)行常規(guī)免疫程序,并實(shí)施規(guī)范的飼養(yǎng)管理。采用Logistic、Gompertz、Von Bertalanffy 3種生長(zhǎng)曲線模型(表1),3 種曲線模型中待定參數(shù)A為極限生長(zhǎng)量(終極生長(zhǎng)量或成熟體重),k 為瞬時(shí)相對(duì)生長(zhǎng)率,B 為常數(shù)尺度[9]。應(yīng)用SPSS 19.0 軟件進(jìn)行統(tǒng)計(jì)分析,利用不同周齡體重?cái)?shù)據(jù)擬合算出模型參數(shù)的最優(yōu)估計(jì)值A(chǔ)、B、k,建立生長(zhǎng)模型,算出拐點(diǎn)體重、拐點(diǎn)周齡和最大周增重,根據(jù)擬合度(R2)評(píng)價(jià)生長(zhǎng)模型。
表1 3種常用動(dòng)物生長(zhǎng)曲線模型
由圖1可知,2 周齡開始公雉雞的體重大于母雉雞,體重增長(zhǎng)速度快。6 周齡后公雉雞生長(zhǎng)發(fā)育加快,公雉雞和母雉雞之間的體重差距逐漸拉大,14周齡以后相對(duì)減慢。
圖1 申鴻七彩雉生長(zhǎng)曲線
申鴻七彩雉公母的3種模型生長(zhǎng)曲線擬合參數(shù)估計(jì)值(A、B、k)、擬合度(R2)、拐點(diǎn)體重、拐點(diǎn)周齡、最大周增重見表2。由表2可知,申鴻七彩雉公雉雞和母雉雞的R2均在0.993以上,3個(gè)生長(zhǎng)曲線模型均適宜申鴻七彩雉公母雉雞體重生長(zhǎng)過程,公母雉雞用Gompertz 模型R2最高,分別為0.997 和0.996,其次為Von Bertalanffy,Logistic 模型R2最低,說(shuō)明Gompertz 更適合于公雉雞和母雉雞生長(zhǎng)曲線的擬合。Logistic 模型擬合的公雉雞和母雉雞的拐點(diǎn)體重最大,分別為741.293 g 和575.018 g;其次為Gom?pertz 模型;拐點(diǎn)體重最小的為Von Bertalanffy 模型。Logistic 模型擬合公母雉雞拐點(diǎn)周齡最大,分別為8.244 周齡和8.837 周齡;Logistic 模型擬合公母雉雞周增重最大,分別為117.124 g和83.234 g。
表2 3種生長(zhǎng)曲線擬合參數(shù)估計(jì)值和擬合度
3 種模型對(duì)申鴻七彩雉公雞和母雞擬合得到的相對(duì)生長(zhǎng)率均是隨日齡的增加而降低(表3),Logis?tic、Gompertz、Von Bertalanffy模型分別在4、6、8周齡前公雉雞的相對(duì)生長(zhǎng)率高于母雉雞。
表3 申鴻七彩雉的相對(duì)生長(zhǎng)率情況
由圖2和圖3可知,3 種模型均能較好地模擬申鴻七彩雉公雉雞和母雉雞的生長(zhǎng)曲線,3 種模型擬合的生長(zhǎng)曲線與申鴻七彩雉生長(zhǎng)發(fā)育的實(shí)測(cè)值曲線基本吻合。Logistic 模型的生長(zhǎng)體重?cái)M合值與實(shí)測(cè)值相比較,4 周齡前周齡明顯大于實(shí)測(cè)值,22~24 周齡略低于實(shí)測(cè)值。Gompertz 模型2周齡前略大于實(shí)測(cè)值,比Logistic 模型生長(zhǎng)體重的擬合值更接近于實(shí)測(cè)值。Von Bertalanffy 模型的生長(zhǎng)體重?cái)M合值與實(shí)測(cè)值最吻合。綜上所述,Gompertz 模型和Von Bertalanffy 模型用于估計(jì)申鴻七彩雉生長(zhǎng)發(fā)育階段體重均較好,Logistic 模型最差。
圖2 公雉雞體重3種模型擬合值和實(shí)測(cè)值生長(zhǎng)曲線
圖3 母雉雞體重3種模型擬合值和實(shí)測(cè)值生長(zhǎng)曲線
本研究發(fā)現(xiàn),申鴻七彩雉公、母雉雞在前2周體重?zé)o較大差異;而2周后,公雉雞體重長(zhǎng)勢(shì)明顯高于母雉雞;14周后,母雉雞生長(zhǎng)速度明顯降低,公雉雞依舊保持一定的生長(zhǎng)速度。張玲等[10]運(yùn)用Logistic曲線方程擬合不同性別余干烏骨雞的早期生長(zhǎng)發(fā)育規(guī)律時(shí)發(fā)現(xiàn),在2~4 周齡余干烏骨雞體尺生長(zhǎng)發(fā)育較快,4周齡以后變慢。因此,在體重快速增長(zhǎng)階段,應(yīng)提供豐富的營(yíng)養(yǎng),為后續(xù)生產(chǎn)做好準(zhǔn)備。
在本研究中,Logistic、Gompertz 和Von Berta?lanffy 3 種模型擬合申鴻七彩雉生長(zhǎng)曲線的R2均高于0.99。進(jìn)一步從統(tǒng)計(jì)和生物學(xué)意義角度評(píng)價(jià)3 種模型發(fā)現(xiàn),Gompertz 模型擬合效果最好,這與河北柴雞[11]、淮南麻黃雞[12]、北京油雞[13]的擬合結(jié)果一致。拐點(diǎn)周齡是畜禽生長(zhǎng)曲線的一個(gè)重要轉(zhuǎn)折點(diǎn),動(dòng)物生長(zhǎng)速度在拐點(diǎn)處由快變慢[14]。拐點(diǎn)體重隨拐點(diǎn)周齡越晚而越大,若體成熟越早,出欄也就越早,經(jīng)濟(jì)效益也會(huì)越大。本研究中Gompertz方程擬合公雉雞和母雉雞的生長(zhǎng)拐點(diǎn)周齡分別為6.43 周和6.91 周;公雉雞的拐點(diǎn)體重和最大周增重分別為571.68 g和112.91 g,母雉雞分別為448.51 g 和79.39 g。楊飛等[15]研究顯示,Gompertz 模型擬合的拜城油雞平均生長(zhǎng)拐點(diǎn)周齡為6.43周,拐點(diǎn)體重為451.326 2 g,這與本研究結(jié)果較為相似。
通過擬合畜禽生長(zhǎng)曲線,可根據(jù)重要生長(zhǎng)節(jié)點(diǎn)預(yù)測(cè)不同生長(zhǎng)階段的飼料需求,調(diào)整飼養(yǎng)方案,節(jié)省生產(chǎn)成本,增加企業(yè)效益。本研究發(fā)現(xiàn),擬合申鴻七彩雉公、母雉雞生長(zhǎng)曲線效果最好的是Gomp?ertz 模型,Gompertz 模型可預(yù)測(cè)和評(píng)估申鴻七彩雉生長(zhǎng)規(guī)律,為申鴻七彩雉生產(chǎn)與育種工作提供理論依據(jù)。