季索菲,侯冠軍,趙秀俠,崔 凱,王利軍,芮 葉
(1. 安徽省農(nóng)業(yè)科學(xué)院水產(chǎn)研究所,水產(chǎn)增養(yǎng)殖安徽省重點(diǎn)實(shí)驗(yàn)室,合肥 230031;2.滁州市隆財(cái)漁業(yè)有限公司, 滁州沙河鎮(zhèn) 239060)
翹嘴鱖(Sinipercachuatsi)也稱桂魚、桂花魚,廣泛分布在我國的江河湖泊中,南至廣東,北至黑龍江水域,以長江中下游地區(qū)最多。翹嘴鱖高蛋白低脂肪,且富含n-3和n-6系列的多不飽和脂肪酸,是久負(fù)盛名的淡水名貴魚類。但是由于近交頻繁,造成鱖魚苗種種質(zhì)退化,抗病性降低,畸形率上升,養(yǎng)殖成本增加等問題日益突出,尤其暴發(fā)性的疾病給養(yǎng)殖生產(chǎn)帶來巨大的經(jīng)濟(jì)損失,嚴(yán)重阻礙了鱖魚養(yǎng)殖業(yè)的正常發(fā)展[1,2]。質(zhì)量和形態(tài)特征等指標(biāo)是選擇育種最常見的數(shù)據(jù)[1],體長、體質(zhì)量和體高之間的關(guān)系在魚類研究中應(yīng)用廣泛,不僅反映了個(gè)體生長情況及種群發(fā)展變化,也是反映魚類生活史的指標(biāo)之一[3,4]。體質(zhì)量增加快是優(yōu)質(zhì)良種的目標(biāo)之一,是決定生產(chǎn)性能的主要指標(biāo),也是養(yǎng)殖者的目標(biāo)[5~7]。外部形態(tài)特征觀察作為種質(zhì)資源鑒定的方法之一,具有簡單、直觀的優(yōu)勢[8]。Gompertz、Logistic和Von Bertalanffy 三種非線性模型常被用于魚類生長曲線研究[9,10],本研究采用三種非線性模型對翹嘴鱖形態(tài)特征以及生長性狀進(jìn)行擬合分析,探討池塘翹嘴鱖生長規(guī)律變化,可為翹嘴鱖的育種工作提供理論指導(dǎo)。
翹嘴鱖取自滁州市隆財(cái)漁業(yè)有限公司池塘。養(yǎng)殖池塘為東西走向長方形池塘(面積為2 000m2),水源充足、無污染,池底淤泥較少,進(jìn)排水方便,水位1.8~2.0m,池塘水體各項(xiàng)指標(biāo)符合淡水養(yǎng)殖用水標(biāo)準(zhǔn)(NY 5051),初始放養(yǎng)密度為1.5尾/m2,實(shí)驗(yàn)魚為1月齡翹嘴鱖,起始體質(zhì)量0.648a±0.309 g,體長為3.176±0.309 cm,體高為1.136±0.352 cm.試驗(yàn)池塘主養(yǎng)魚種為安徽本地繁育的翹嘴鱖,鱖魚與配套餌料魚培育池塘面積按1∶3配比,餌料魚為團(tuán)頭魴魚、鯉魚、鯽魚、鯪魚、鰱魚、草魚。每天巡塘,及時(shí)管理。
養(yǎng)殖時(shí)間為2019年5~11月,每月20日上午8點(diǎn)從池塘中(2 000m2)抬網(wǎng)隨機(jī)采樣30尾,暫養(yǎng)于暫存桶中,桶內(nèi)放置經(jīng)清水洗凈后并經(jīng)5% NaCl溶液消毒后的喜旱蓮子草,同時(shí)采取充氧、降溫措施,待翹嘴鱖應(yīng)激反應(yīng)過后進(jìn)行稱重測量記錄。長度使用游標(biāo)卡尺進(jìn)行測量,精確度0.01mm,質(zhì)量使用電子天平,精確度0.001g.測定的數(shù)據(jù)為翹嘴鱖體質(zhì)量、體長和體高,以平均值±標(biāo)準(zhǔn)偏差表示。
塘口試驗(yàn)結(jié)束后采用Excel 2010進(jìn)行統(tǒng)計(jì)分析,繪制翹嘴鱖體質(zhì)量、體長和體高1~6月齡的累計(jì)生長曲線、絕對生長曲線、相對生長曲線,以及體質(zhì)量和體長的關(guān)系曲線。
絕對生長率:G絕w=(W2-W1)/(t2-t1)
G絕L=(L2-L1)/(t2-t1)
G絕H=(H2-H1)/(t2-t1)
相對生長率:G相W=(W2-W1)/[W1(t2-t1)]×100%
G相L= (L2-L1)/[L1(t2-t1)]×100%
G相H= (H2-H1)/[H1(t2-t1)]×100%
其中,G為生長率;W為體質(zhì)量;t為時(shí)間;L為體長;H為體高。運(yùn)用SPSS 20.0軟件Regression中的Nonlinear進(jìn)行分析,采用表1中Gompertz、Logistic和Von Bertalanffy曲線模型對翹嘴鱖體質(zhì)量、體長和體高進(jìn)行生長擬合,計(jì)算模型參數(shù)的最優(yōu)估計(jì)值,建立生長模型,再依據(jù)擬合精度(R2)評價(jià)生長模型。
表1 用于擬合的三種非線性模型
翹嘴鱖1~6月的生長數(shù)據(jù)見表2.經(jīng)過6個(gè)月的池塘養(yǎng)殖,收獲1.38尾/m2,體質(zhì)量,體長和體高每月顯著性增加,體質(zhì)量從0.648±0.309g增長至387.3±50.425g,體長從3.176±0.309cm增長至25.056±1.079cm,體高從3.561±0.545cm增長至8.996±0.371cm.
表2 翹嘴鱖1~6月齡體質(zhì)量和體長的測量值(平均值±標(biāo)準(zhǔn)差)
圖1~9是1~6月齡翹嘴鱖質(zhì)量、體長和體高的累計(jì)生長曲線,月絕對生長曲線和月相對生長曲線。由圖1~3可見,翹嘴鱖體質(zhì)量、體長和體高隨月齡的增加而增加,其中4月齡之后體質(zhì)量增長較快,體長和體高在孵化后到4月齡之間生長較快。由圖4~6可見,體質(zhì)量和體高的絕對生長率趨勢一致,先上升后下降,而體長則呈現(xiàn)出下降上升后下降的趨勢,三者均在4月齡時(shí)達(dá)到峰值。圖7~9可見,體質(zhì)量和體長的相對生長率一直呈下降趨勢,體高相對生長率呈先下降后小幅度增加的趨勢。
圖1 翹嘴鱖1~6月齡體質(zhì)量累計(jì)生長曲線 圖2 翹嘴鱖1~6月齡體長累計(jì)生長曲線
圖3 翹嘴鱖1~6月齡體高累計(jì)生長曲線 圖4 翹嘴鱖1~6月齡體質(zhì)量月絕對生長曲線
圖7 翹嘴鱖1~6月齡體質(zhì)量月相對生長曲線 圖8 翹嘴鱖1~6月齡體長月相對生長曲線
圖9 翹嘴鱖1~6月齡體高月相對生長曲線 圖10 翹嘴鱖體質(zhì)量和體長的生長關(guān)系
根據(jù)測量的體質(zhì)量和體長數(shù)值進(jìn)行曲線擬合,運(yùn)用Excel選擇冪函數(shù)模型W=aLb對生長進(jìn)行擬合,W代表體質(zhì)量,L代表體長,擬合關(guān)系式為W=0.0174L3.108(R2=0.9992)。如圖10所示。
采用Gompertz、Logistic和Von Bertalanffy 3種非線性模型對翹嘴鱖的體質(zhì)量、體長和體高變化規(guī)律進(jìn)行擬合,結(jié)果見表3,擬合精度R2均達(dá)到0.921以上。安徽本地繁育翹嘴鱖體質(zhì)量、體長和體高的隨機(jī)觀測值和模型估計(jì)值見表4.用Logistic模型對1月齡和2月齡的翹嘴鱖體質(zhì)量的估計(jì)值與實(shí)際測定值存在較大的誤差,用Von Bertalanffy模型對1~3月齡的翹嘴鱖體質(zhì)量的估計(jì)值與實(shí)際測定值存在較大的誤差,而Gompertz模型估計(jì)結(jié)果較好。對于體長和體高的估計(jì),三種模型的預(yù)測效果均較好,只有Von Bertalanffy模型在1月齡與其他兩個(gè)模型估計(jì)值存在較大不同,結(jié)合估計(jì)值和擬合精度R2綜合考慮,決定采用Gompertz模型對體長進(jìn)行擬合,采用Logistic模型對體高進(jìn)行擬合。
表3 用3種模型擬合翹嘴鱖體質(zhì)量、體長和體高生長曲線的參數(shù)估計(jì)值和擬合度
表4 翹嘴鱖體質(zhì)量、體長和體高實(shí)際觀測值和擬合曲線估計(jì)值的比較
劉賢德等人[11,12]研究表明體長、頭長、全長、體高等形態(tài)性狀是影響體質(zhì)量的主要因素,其中體長和體高對翹嘴鱖早期體質(zhì)量影響最大[13],因此本試驗(yàn)以1~6月翹嘴鱖體長和體高為主要研究形態(tài)性狀。試驗(yàn)研究發(fā)現(xiàn)對翹嘴鱖體質(zhì)量有顯著影響的表型性狀為體長和體高,與張進(jìn)等人研究結(jié)果一致,在今后的生產(chǎn)實(shí)踐中可通過優(yōu)化這兩種表型性狀達(dá)到增加體質(zhì)量的目的,以促進(jìn)翹嘴鱖優(yōu)良品種選育和增強(qiáng)日??茖W(xué)養(yǎng)殖管理。
通過圖1~9可得出,翹嘴鱖的體質(zhì)量、體長、體高都隨著年齡的增長不斷增加,絕對生長率在4個(gè)月之前達(dá)到最大值,5~6個(gè)月呈現(xiàn)下降趨勢,相對生長率隨著年齡的增長不斷下降,說明體質(zhì)量、體長、體高會(huì)隨著年齡的增長而增長,但相對生長率會(huì)逐漸趨于零,這也與許多其他魚類的生長速度變化趨勢相仿,如對海南絲蝴蝶魚的研究表明體長第一年增長最快,以后增加會(huì)逐步放慢,體重增重率也隨年齡的增長逐步下降[14];對美洲紅點(diǎn)鮭(Salvelinus fontinalis)和白斑點(diǎn)鮭(Salvelinus leucomaenis)的研究表明體質(zhì)量的生長速度有一個(gè)峰值,峰值后生長速度開始減慢[9]。
魚類的體質(zhì)量和體長一般存在冪指數(shù)關(guān)系:W=aLb,一般指數(shù)b的值在2.5~4.0,根據(jù)本試驗(yàn)測量的數(shù)值得出的關(guān)系式為W=0.0174L3.108(R2=0.9992),其中b為3.108,也是在2.5~4.0這個(gè)范圍內(nèi),這與相關(guān)報(bào)道相符[15~17]。本試驗(yàn)中,冪指數(shù)b是反映翹嘴鱖生長情況的一個(gè)特征數(shù),有研究認(rèn)為冪指數(shù)是否接近3是判斷該魚類是否勻速生長的指標(biāo)[18,19],本指數(shù)關(guān)系式中b值為3.108,接近于3,說明該測量群體為近似勻速生長的魚類。翹嘴鱖在由小到大的生長過程中,三個(gè)線度(長、寬、高)方向生長的速度相等,體形隨個(gè)體生長按等比例放大[4],在1~6個(gè)月具有體形不變和比重不變的特征。R2大于0.99,說明擬合度效果好,但該生長關(guān)系式是由1~6月生長性狀測量值得出,是否適合翹嘴鱖其他生長時(shí)期,還有待進(jìn)一步研究。
擬合精度R2越接近于1,說明回歸模型的擬合效果較好。本試驗(yàn)采用Gompertz、Logistic和VonBertalanffy 3種生長模型對安徽本地翹嘴鱖體質(zhì)量、體長、體高的生長進(jìn)行擬合,R2最小值達(dá)到0.921,三種生長模型均有較好的擬合度。但從實(shí)際擬合效果看,Gompertz模型對體質(zhì)量和體長擬合效果較好,體質(zhì)量和體長的拐點(diǎn)月齡分別為3.52月齡和2.09月齡,拐點(diǎn)處體質(zhì)量和體長分別為160.51g和9.91cm.這與郭媛等對美洲紅點(diǎn)鮭體質(zhì)量和體長擬合結(jié)果一致,說明翹嘴鱖體長發(fā)育較體質(zhì)量發(fā)育時(shí)間更快。2~3月翹嘴鱖體長絕對生長率一直維持在較高水平,可能與其身體內(nèi)臟發(fā)育有關(guān)。3~4月是翹嘴鱖體重增長的關(guān)鍵時(shí)期,在此時(shí)期加強(qiáng)餌料魚的投喂和日常塘口管理對其體重增加尤為重要。Logistic模型對體高擬合較好,體高的拐點(diǎn)月齡為2.44月齡,拐點(diǎn)處體高為4.52cm.體高和體長拐點(diǎn)均出現(xiàn)在2月齡左右,這可能是翹嘴鱖早期骨骼和內(nèi)臟發(fā)育所致。從體質(zhì)量和體長的生長曲線中,根據(jù)實(shí)際觀測值和模型值進(jìn)行比照,5~6月齡生長逐漸減緩,體質(zhì)量隨著年齡的增長而增加,近似勻速增長,這也再次驗(yàn)證了冪指數(shù)b的判斷[20,21]。
體質(zhì)量是眾多水產(chǎn)動(dòng)物遺傳育種改良的重要目標(biāo)[22],翹嘴鱖體質(zhì)量除了遺傳因素外,養(yǎng)殖環(huán)境和管理水平的影響也很大,難以直接確定品種優(yōu)劣。本研究確立了體質(zhì)量與體長冪指數(shù)關(guān)系,在對翹嘴鱖進(jìn)行育種時(shí),除了對體質(zhì)量進(jìn)行直接選育外,還可以對體長進(jìn)行輔助選育,提高育種效率。此外,通過對體長、體質(zhì)量進(jìn)行預(yù)測,可對選育品系進(jìn)行提前篩選,加速育種進(jìn)程,減低育種成本。
冪指數(shù)生長關(guān)系式可反映魚體生長的內(nèi)在聯(lián)系及外界環(huán)境因素對生長的關(guān)系[10,23],利用本試驗(yàn)得出的生長關(guān)系式,在漁業(yè)生產(chǎn)過程中,若在相同養(yǎng)殖條件下月齡相同的翹嘴鱖生長發(fā)現(xiàn)異常,可查找原因,調(diào)整養(yǎng)殖措施,這在一定程度上也給翹嘴鱖養(yǎng)殖管理提供參考。
總之,通過對安徽本地翹嘴鱖體質(zhì)量、體長和體高進(jìn)行測量和分析,找出對翹嘴鱖體質(zhì)量、體長、體高較合適的非線性模型,得到其冪指數(shù)方程W=0.0174L3.108(R2=0.9992)。本研究的結(jié)果為研究翹嘴鱖早期生長性狀提供一定的參考依據(jù),同時(shí)為翹嘴鱖早期養(yǎng)殖管理中合理投喂餌料魚提供參考,為科學(xué)人工養(yǎng)殖和選育翹嘴鱖提供技術(shù)指導(dǎo),但對翹嘴鱖肌肉營養(yǎng)和肉質(zhì)等品質(zhì)有待進(jìn)一步研究。