開口餌料對泥鰍苗種存活率的影響及其生長曲線擬合分析

陸專靈 趙忠添 侯樹鑒 盧智發 李青 李登明 張益峰 黃黎明

摘要：【目的】確定最適生長曲線擬合模型及篩選出最適開口餌料，為泥鰍人工育苗提供科學依據。【方法】分別以豐年蟲、輪蟲、水蚯蚓、蛋黃、豆漿和鰻魚料為開口餌料，培育泥鰍苗種20 d后比較其對泥鰍苗種存活率及生長性能的影響，并采用Logistic、Gompertz和Bertalanffy 3種非線性模型對培育20 d泥鰍苗種的體長數據進行生長曲線擬合與回歸分析。【結果】以不同開口餌料培育20 d后，泥鰍苗種的存活率排序為：豐年蟲組>輪蟲組>水蚯蚓組>蛋黃組>鰻魚料組>豆漿組，各處理組間存在顯著差異（P<0.05）。其中，以豐年蟲組的存活率最高，為52.7%，豆漿組的最低，僅為10.3%。在生長性能方面，也以豐年蟲組泥鰍苗種的生長速率最快，豆漿組的生長速率最慢。Logistic、Gompertz和Bertalanffy模型均能較好地擬合泥鰍苗種的生長曲線，但Logistic模型擬合獲得的生長拐點日齡比Gompertz和Bertallanffy模型的擬合結果略晚，且從擬合度（R2）來看，Bertallanffy模型的R2更接近于1.000，即Bertallanffy模型更適合用于擬合泥鰍苗種的生長曲線。【結論】以豐年蟲為開口餌料能夠顯著提高泥鰍苗種培育存活率和生長速度，若兼顧養殖成本，則以輪蟲更適合作為泥鰍苗種的開口餌料。Bertallanffy模型更適合用于擬合泥鰍苗種的生長曲線，擬合得出的泥鰍苗種生長拐點在7～10日齡。

關鍵詞： 泥鰍苗種；開口餌料；存活率；生長曲線；Bertallanffy模型

中圖分類號： S966.4 文獻標志碼：A 文章編號：2095-1191（2016）08-1411-05

Abstract：【Objective】The present experiment was conducted to find out optimal growth curve fitting model and select the optimal initial bait， in order to provide scientific basis for artificial breeding of Misgurnus anguillicaudatus larvae. 【Method】Artemia salina， rotifera， egg yolk， soybean milk， tubificidae and eel feed were used， as initial feeds to breed M. anguillicaudatus larvae for 20 d， respectively. Then the survival rate and growth performance were determined， so as to compare effects of initial feeds on survival rate and growth performance. Moreover， the data about body length of M. anguillicaudatus larvae were determined for growth curve fitting and regression analysis using 3 nonlinear models viz.， Logistic， Gompertz and Bertalanffy. 【Result】After M. anguillicaudatus larvae were bred with different initial feeds for 20 days， the order of their survival rates was A. salina group>rotifera group>tubificidae group>egg yolk group>eel feed group>soybean milk group， there were significant differences in survival rate between different groups（P<0.05）. The survival rate of M. anguillicaudatus larvae in A. salina group was the highest（52.7%）， while that of M. anguillicaudatus larvae in soybean milk group was the lowest（10.3%）. In term of growth performance， M. anguillicaudatus larvae in fairy shrimp group had the fastest growth rate， those in soybean milk group had the slowest one. The Logistic， Gompertz and Bertalanffy nonlinear model all could well fit growth curve of M. anguillicaudatus， but the day old at inflection point of growth curve fitted by Logistic model was slightly later than those of growth curves fitted by Gompertz and Bertalanffy models. Furthermore， the fitting degree（R2） of Bertalanffy models was more closer to 1.000， so the Bertalanffy model was more suitable for fitting growth curve of M. anguillicaudatus larvae. 【Conclusion】A. salina used as initial bait can significantly improve survival rate and growth rate of M. anguillicaudatus larvae， while considering production cost， rotifera was more suitable as initial bait for M. anguillicaudatus larvae. The Bertalanffy model is more suitable for fitting growth curve of M. anguillicaudatus larvae， and the inflection point of fitting growth were at 7-10 days old.

Key words： Misgurnus anguillicaudatus larvae； initial feed； survival rate； growth curve； Bertalanffy model

0 引言

【研究意義】動物的生長一般都具有非線性特性，因此可用非線性數學模型來描述其生長規律，常用的非線性數學模型有Logistic、Gompertz和Bertalanffy 3種模型。在養殖過程中，可利用不同時間序列的動物體重、體長或其他指標擬合目的性狀的生長曲線，建立數學模型并估測其生長參數（Knízetová et al.，1991；楊運清和繆堯源，1992）。泥鰍（Misgurnus anguillicaudatus）隸屬于鯉形目（Cypriniformes）鰍科（Cobitidae），在我國及東南亞地區廣泛分布，是具有較高經濟價值的淡水魚類之一，其當前的研究熱點是開口餌料選擇及生長發育規律分析。因此，研究開口餌料對泥鰍存活率的影響及運用非線性數學模型對其生長曲線進行擬合分析，對提高泥鰍培育效果具有重要意義。【前人研究進展】近年來，泥鰍養殖產業得到快速發展。徐如衛等（2014）總結了泥鰍循環水集約化養殖的技術要點，且認為其養殖模式已由傳統池塘圍網養殖模式向高密度集約化養殖模式轉變。朱光來和趙子明（2014）對泥鰍繁殖的最佳催產劑配比、催產溫度進行比較分析，總結出泥鰍的人工繁育技術。在實際生產中，泥鰍苗種培育存活率低一直制約著產業的發展，而開口餌料是影響泥鰍苗種存活率和后期生長發育的最主要因素。邵力和童勇義（2002）研究發現，采用蛋黃、輪蟲、藻類作為泥鰍仔魚的開口餌料較適宜；蘇應兵和楊代勤（2010）通過比較不同開口餌料對泥鰍仔魚存活率和生長的影響，發現輪蟲組的存活率最高。由于動物生長發育一般都具有非線性特征，因此常用非線性模型來描述其生長發育規律，如Logistic、Gompertz和Bertalanffy等經典的非線性生長模型已廣泛應用于畜禽及水產經濟動物（Kwakkel et al.，1993；Sales et al.，1997；李俊營等，2014）。何小燕等（2011）分別用Logistic、Gompertz和Bertalanffy模型分析大口黑鱸早期的生長發育規律；肖俊等（2012）采用Logistic、Gompertz和Bertallanffy等非線性模型對尼羅羅非魚的生長發育規律進行研究。【本研究切入點】至今，除了何國森等（2015）曾采用Bertalanffy、Gompertz和Logistic模型擬合了泥鰍苗種的生長曲線外，鮮見其他有關開口餌料對泥鰍生長曲線擬合分析的研究報道。【擬解決的關鍵問題】以輪蟲、水蚯蚓、豐年蟲、蛋黃、豆漿和鰻魚料為開口餌料，比較其對泥鰍苗種存活率及生長性能的影響，再運用Logistic、Gompertz和Bertalanffy模型對泥鰍苗種體長生長曲線進行擬合分析，確定最適生長曲線擬合模型及篩選出最適開口餌料，為泥鰍人工育苗提供科學依據。

1 材料與方法

1. 1 試驗用泥鰍及餌料

泥鰍苗種由廣西水產科學研究院人工繁殖提供，即經72 h自行發育出膜3 d的魚苗，其卵黃囊接近消失，健康活潑、規格一致，初始全長4.30±0.15 mm。開口餌料分別為豐年蟲、輪蟲、蛋黃、豆漿、水蚯蚓和鰻魚料。試驗用水為曝氣自來水，水體指標：水溫（26±1）℃，pH 7.6±0.5，溶解氧>5 mg/L，氨氮<0.1 mg/L。

1. 2 試驗方法

根據開口餌料共設6個處理組，每組300尾泥鰍苗種。每天投喂3次，投喂量以吃飽為準。每次投喂1 h后用虹吸方法清除剩余的餌料殘渣及糞便，并加水至原水位高度（10 cm）。試驗周期20 d，每天隨機抽取30尾泥鰍進行全長測量及存活率統計。

1. 3 數據統計及曲線擬合模型

試驗數據以Excel 2007進行整理，采用SPSS 17.0進行方差分析，用GraphPad Prism 5.0進行生長曲線繪制，再根據非線性回歸計算各模型中的A、B、K參數值，建立生長曲線模型，計算拐點體長、拐點日齡和擬合度（R2），各計算公式見表1。式中，Lt為t時的體長估計值，A為最大理論生長點，B為常數尺度，K為瞬時相對生長率，L為拐點體長；At為拐點日齡。

2 結果與分析

2. 1 開口餌料對泥鰍苗種存活率的影響

由表2可知，在試驗前5 d，豐年蟲組的泥鰍苗種存活率最高，為83.2%，其次是蛋黃組，鰻魚料組的存活率最低（69.5%）。至試驗中期（10～15 d），各處理組泥鰍苗種的存活率均呈下降趨勢，且各處理組間的差異也越來越明顯。至試驗結束時（第20 d），泥鰍苗種的培育存活率排序為：豐年蟲組>輪蟲組>水蚯蚓組>蛋黃組>鰻魚料組>豆漿組。其中，豐年蟲組的泥鰍苗種存活率（52.7%）顯著高于其他處理組（P<0.05，下同） ，而豆漿組的存活率最低，僅為10.3%。

2. 2 開口餌料對泥鰍苗種生長的影響

由表3可知，在試驗前期5 d，豐年蟲組和輪蟲組泥鰍苗種的體長分別為6.80和6.54 mm，均顯著高于其他處理組。至試驗中期（10～15 d），仍以豐年蟲組泥鰍苗種的生長最快，其體長顯著高于其他處理組，其次是輪蟲組和水蚯蚓組，二者間差異不顯著（P>0.05，下同），豆漿組的生長最慢。至試驗結束時（第20 d），豐年蟲組泥鰍苗種的全長為14.86 mm，其次是輪蟲組和水蚯蚓組，分別為13.89和13.56 mm，與蛋黃組、鰻魚料組、豆漿組間存在顯著差異。綜上所述，豐年蟲組泥鰍苗種的生長速率最快，而豆漿組的生長速率最慢。

2. 3 開口餌料對泥鰍全長生長曲線擬合效果的影響

不同開口餌料對泥鰍苗種全長的Logistic、Gompertz和Bertalanffy模型擬合生長曲線參數估計值及R2見表4。由表4可以看出，3種模型均能較好地模擬出泥鰍的全長生長曲線，其中Bertalanffy模型的R2最高，其次是Gompertz模型，Logistic模型的略低。R2是衡量模型擬合效果的重要指標之一，R2越接近于1.000，曲線擬合效果越佳。因此，確定Bertalanffy模型的擬合效果最佳。在Bertalanffy模型擬合生長曲線中，不同開口餌料處理組的R2也存在明顯差異，其排序為：輪蟲組（1.000）>豐年蟲組（0.999）=水蚯蚓組（0.999）>豆漿組（0.998）>蛋黃組（0.990）>鰻魚料組（0.989）。在拐點體長和拐點日齡參數中，Logistic模型的拐點體長和拐點日齡最高，其次為Gompertz模型，Bertalanffy模型最低。以Bertalanffy模型進行分析可知，豐年蟲組、蛋黃組和豆漿組泥鰍苗種的生長拐點日齡為7 d，輪蟲組的生長拐點日齡為8 d，水蚯蚓和鰻魚料組的生長拐點日齡為10 d。

3 討論

魚苗的初孵苗種依靠內源性餌料卵黃生存，在從內源性營養轉化到混合性營養或外源性營養的過渡階段，通常會出現大批量死亡（殷名稱，1995；龔全等，2015）。魚苗初次攝食外源性餌料的種類、顆粒大小和成分對其存活與生長發育有密切關系。泥鰍苗種的植物性開口餌料包括多種單細胞藻類、光合細菌等，動物性開口餌料包括輪蟲、豐年蟲、水蚯蚓及格足類和枝角類等，人工開口餌料有蛋黃、豆漿、魚粉、鰻魚料等。陳舒泛（2001）研究表明，輪蟲適口性較好，游泳能力弱，粗蛋白含量為58.20%，粗脂肪含量為14.20%，且含有苗種生長發育所必需的不飽和脂肪酸，是目前使用最廣泛的魚苗開口餌料之一。潘茜等（2009）研究表明，豐年蟲無節幼體因個體較小、運動能力弱、適口性好，且具有增強免疫力和活力等特點，現已廣泛應用在魚、蝦、蟹繁殖生產中，但豐年蟲成本較高，通常只作為經濟價值較高的魚苗開口餌料。蛋黃及鰻魚料投喂后會沉入水中，靜止的餌料無法促使魚苗產生攝食刺激，餌料利用很低，極易造成水體污染；豆漿最容易惡化水質，魚苗培育存活和生長效果最差。水蚯蚓碎段無游泳能力，顆粒大小不一，前期開口效果較差，導致苗種生長不勻，但其蛋白質含量高達70%，粗蛋白中氨基酸齊全，對泥鰍苗種后期生長有較好的促進作用。本研究結果表明，使用動物性開口餌料豐年蟲、輪蟲和水蚯蚓培育泥鰍苗種的存活率顯著高于蛋黃、鰻魚料和豆漿等人工開口餌料。雖然豐年蟲組的泥鰍苗種存活率及生長速度均最高，輪蟲組的次之，但兼顧規模養殖過程中的生產成本，發現輪蟲更適合作為泥鰍苗種的開口餌料（蘇應兵和楊代勤，2010）。此外，目前市場上尚缺少質量可靠的泥鰍苗種培育階段人工開口餌料，因此還需深入研究泥鰍苗種消化道中各種消化酶種類、活性及其食性轉換機理，不斷探索泥鰍苗種開口餌料中動、植性餌料的比例及蛋白含量、餌料大小、形態、沉降特點等，才能提高泥鰍苗的培育效果。

運用Logistic、Gompertz和Bertalanffy非線性數學模型對動物的生長發育進行擬合時，因研究對象的生長環境、管理水平等因素差異，其擬合結果也不同。王新安等（2011）研究表明，在大菱鲆選育中不同組別的生長曲線以Gompertz和Bertallanffy模型擬合效果較好，而對照組以Logistic模型擬合效果較優。肖俊等（2012）對尼羅羅非魚生長曲線進行擬合分析，結果發現以Bertallanffy模型的擬合效果較好。李軍等（2015）在擬合唇■的生長方程擬合過程中，發現Bertallanffy模型的R2最高。本研究選用泥鰍苗種體長作為指標，其擬合生長曲線表明，Logistic、Gompertz和Bertalanffy模型均能較好地擬合生長曲線，但Logistic模型擬合獲得的生長拐點日齡比Gompertz和Bertallanffy模型的擬合結果略晚，且從R2來看，Bertallanffy模型的R2更接近于1.000，即Bertallanffy模型更適合用于擬合泥鰍苗種的生長曲線。生長拐點是機體生長速率最快的時間，其出現越早表明生長速率越高（項云和王偉，2000；肖俊等，2012）。根據生長拐點可篩選出適宜的營養需求（開口餌料），進而提高機體生長速度及培育存活率。本研通過Bertallanffy模型擬合的泥鰍苗種生長曲線，結果發現豐年蟲組泥鰍苗種的生長拐點日齡為7 d，明顯早于其他開口餌料組，表明以豐年蟲為開口餌料時泥鰍苗種的生長速率最快，與何國森等（2015）的研究結果相似。

4 結論

以豐年蟲為開口餌料能夠顯著提高泥鰍苗種培育成活率和生長速度，若兼顧養殖成本，則以輪蟲更適宜作為泥鰍苗種的開口餌料。Bertallanffy模型更適合用于擬合泥鰍苗種的生長曲線，擬合得出的泥鰍苗種生長拐點在7～10日齡。

參考文獻：

陳舒泛. 2001. 淡水漁業開口餌料輪蟲產業化生產[J]. 南京曉莊學院學報，17（4）：51-56.

Chen S F. 2001. Mass production of rotifer as first feeding for freshwater fish larviculture[J]. Journal of Nanjing Xiaozhuang College，17（4）：51-56.

龔全，劉亞，賴見生，趙剛，杜軍，李麗鵑. 2015. 不同開口餌料對達氏鱘魚苗生長的影響[J]. 西南農業學報，28（5）：2297-2300.

Gong Q，Liu Y，Lai J S，Zhao G，Du J，Li L J. 2015. Effects of initial feeding on growth of Acipenser dabryanus Dumeril[J]. Southwest China Journal of Agricultural Sciences， 28（5）：2297-2300.

何國森，陳度煌，曾占壯，王劍鋒，王壽昆. 2015. 不同餌料對泥鰍苗種前期培育效果的研究[J]. 福建水產，37（1）：62-67.

He G S，Chen D H，Zeng Z Z，Wang J F，Wang S K. 2015. Effects of different initial feed on rearing of larva Misgurnus anguillicaudatus[J]. Journal of Fujian Fisheries，37（1）：62-67.

何小燕，白俊杰，樊佳佳，李勝杰，劉小林. 2011. 大口黑鱸早期生長發育規律的研究[J]. 大連海洋大學學報，26（1）：23-29.

He X Y，Bai J J，Fan J J，Li S J，Liu X L. 2011. Analysis of fitting growth model in early stages of largemouth bass Microperus salmoides[J]. Journal of Dalian Ocean University，26（1）：23-29.

李軍，駱小年，金廣海，魏宏祥. 2015. 網箱養殖唇■的生長性能與生長模型[J]. 廣東海洋大學學報，35（1）：99-103.

Li J，Luo X N，Jin G H，Wei H X. 2015. Growth performance and models of Hemibarbus labeo cultrued in net cages[J]. Journal of Guangdong Ocean University，35（1）：99-103.

李俊營，詹凱，楊秀娟，張偉，劉偉，羅聯輝. 2014. 黃山黑雞體重生長發育規律研究[J]. 西南農業學報，27（3）：1295-1299.

Li J Y，Zhan K，Yang X J，Zhang W，Liu W，Luo L H. 2014. Research of growth and development rules of body weight of Huangshan Black Chicken[J]. Southwest China Journal of Agricultural Sciences，27（3）：1295-1299.

潘茜，葉金云，陳建明，王友慧，沈斌乾，郭建林. 2009. 古根豐年蟲營養成分分析[J]. 淡水漁業，39（5）：43-46.

Pan Q，Ye J Y，Chen J M，Wang Y H，Shen B Q，Guo J L. 2009. Nutritional composition analysis of Branchinella kugenumaesis[J]. Freshwater Fisheries，39（5）：43-46.

邵力，童勇義. 2002. 泥鰍仔魚開口餌料與生長的初步研究[J]. 浙江海洋學院學報（自然科學版），21（3）：216-219.

Shao L，Tong Y Y. 2002. Study on the initial bait and growth of Misgurnus anguillicaudatus[J]. Journal of Zhejiang Ocean University（Natural Science），21（3）：216-219.

蘇應兵，楊代勤. 2010. 不同開口餌料對泥鰍仔魚成活率和生長的影響[J]. 長江大學學報（自然科學版），7（3）：37-39.

Su Y B，Yang D Q. 2010. Effects of different initial feeds on survival and growth of larva Misgurnus anguillicaudatus[J]. Journal of Yangtze University（Natural Science Edition），7（3）：37-39.

王新安，馬愛軍，黃智慧，劉慶明，郭黎，楊志，曲江波. 2011. 大菱鲆（Scophthalmus maximus）大規模家系選育優良家系的生長發育規律[J]. 海洋與湖沼，42（4）：266-273.

Wang X A，Ma A J，Huang Z H，Liu Q M，Guo L，Yang Z，Qu J B. 2011. Growth patterns of selectively bred turobot Scophthalmus maximus[J]. Oceanology et Limnology Sinica，42（4）：266-273.

項云，王偉. 2000. 不同育雛方式下雛鴿的生長曲線[J]. 浙江農業學報，12（4）：213-216.

Xiang Y，Wang W. 2000. Growth curves of pigeons squab fed with different feeding ways[J]. Acta Agriculturae Zhejiange-

nsis，12（4）：213-216.

肖俊，凌正寶，唐瞻楊，羅永巨，郭忠寶，郭恩彥，嚴欣，張明，甘西. 2012. 尼羅羅非魚（Oreochromis niloticusi）生長相關分析及生長模型構建[J]. 海洋與湖沼，43（6）：1272-1278.

Xiao J，Ling Z B，Tang Z Y，Luo Y J，Guo Z B，Guo E Y，Yan X，Zhang M，Gan X. 2012. Study on correlation analysis and growth model in nile tilapia Oreochromis niloticusi[J]. Oceanology et Limnology Sinica，43（6）：1272-1278.

徐如衛，楊福生，吳翩，李倩，俞奇力，劉士力，王雨辰，胡廷尖. 2014. 泥鰍循環水集約化養殖的技術要點[J]. 安徽農學通報，20（1-2）：128-129.

Xu R W，Yang F S，Wu P，Li Q，Yu Q L，Liu S L，Wang Y C，Hu T J. 2014. Circulating water intensive techniques of Misgurnus anguillicaudatus[J]. Anhui Agricultural Science Bulletin，20（1-2）：128-129.

楊運清，繆堯源. 1992. 動物生長曲線擬合方法的探討[J]. 畜牧獸醫學報，23（3）：219-224.

Yang Y Q，Miao R Y. 1992. Preliminary discussion on the fitting methods of animal growth curves[J]. Acta Veterinaria et Zootechnica Sinica，23（3）：219-224.

殷名稱. 1995. 魚類仔魚期的攝食和生長[J]. 水產學報，19（4）：335-342.

Yin M C. 1995. Feeding and growth of the larva stage of fish[J]. Journal of Fisheries of China，19（4）：335-342.

朱光來，趙子明. 2014. 泥鰍的人工繁殖與苗種培育初步研究[J]. 中國農學通報，30（5）：61-67.

Zhu G L，Zhao Z M. 2014. The preliminary study on artificial propagation and fry rearing of Misgurnus anguillicaudatus[J]. Chinese Agricultural Science Bulletin，30（5）：61-67.

Knízetová H，Hyánek J，Kníze B，Roubícek J. 1991. Analysis of growth curves of fowl. I. Chickens[J]. British Poultry Science，32（5）：1027-1038.

Kwakkel R P，Ducro B J，Koops W J. 1993. Multiphasic analysis of growth the body and its chemical components in white Leghorn pullets[J]. Poultry Science，72（8）：1421-1432.

Sales J，Du Preez J J，van Niekerk S A. 1997. Chemical composition of the pearl grey guinea fowl[J]. British Poultry Science，38（2）：223.

（責任編輯 蘭宗寶）