水庫放養模式下日本鰻鱺養成品形態性狀對體重和凈體重的影響效應

斯烈鋼，袁向陽，楊　磊，朱衛東，，申屠琰，申屠基康，王志錚

（1.浙江海洋學院水產學院，浙江舟山　316022；2.寧波市海洋與漁業研究院，浙江寧波　315012；3.余姚市水產技術推廣中心，浙江余姚　315400）

水庫放養模式下日本鰻鱺養成品形態性狀對體重和凈體重的影響效應

斯烈鋼1，2，袁向陽1，楊磊1，朱衛東1，3，申屠琰3，申屠基康1，2，王志錚1

（1.浙江海洋學院水產學院，浙江舟山316022；2.寧波市海洋與漁業研究院，浙江寧波315012；3.余姚市水產技術推廣中心，浙江余姚315400）

隨機選取25尾水庫放養模式養成的同生群日本鰻鱺為研究對象，逐一測量其全長（X1）、肛長（X2）、頭長（X3）、吻長（X4）、下頜長（X5）、頭高（X6）、肛部體高（X7）、頭寬（X8）、眼間距（X9）、肛部體寬（X10）、肛圍（X11）、體重（BW）、凈體重（NW）等13項表型性狀后，采用相關分析、通徑分析和偏回歸分析方法研究了日本鰻鱺形態性狀對體重和凈體重的影響效應。結果表明：（1）所測11項形態性狀與體重、凈體重的相關系數均達到顯著水平（P＜0.05）；（2）經通徑分析，被保留的肛部體寬（X10）、肛長（X2）、頭寬（X8）、全長（X1）和下頜長（X5）等5項形態性狀對體重、凈體重相關指數分別為0.985和0.986。其中，對體重和凈體重的直接作用均以肛部體寬（X10）為最大，分別為0.545和0.527，對體重和凈體重產生的間接作用均以頭寬（X8）為最大，分別為0.715和0.699；（3）經多元回歸分析，建立被保留的形態性狀與體重和凈體重間的回歸方程分別為BW=17.669X10+1.785X2+8.292X8+ 0.153X1-13.900X5-658.897和NW=15.505X10+1.643X2+8.209X8+0.139X1-12.030X5-606.031。

形態性狀；體重；凈體重；日本鰻鱺；水庫放養模式

體重系水產養殖動物最直接的目標選育性狀之一，受限于選育過程中的測量誤差大、受測對象易受損及凈體重無法直接獲取等原因，借助多元分析方法構建形態性狀與體重、凈體重間的回歸方程，以切實提高體重選育精度，已受到業內學者的廣泛關注和重視，目前已涉及的養殖魚類主要有牙鲆Paralichthys olivaceus［1］、草魚Ctenopharyngodon idellus［2］、黑斑狗魚Esox reicherti［3］、大菱鲆Scophthalmusmaximus［4］、泥鰍Misgurnus anguillicaudatus［5］、大鱗副泥鰍Paramisgurnus dabryanus［6］等。其中，鰻鱺屬種類僅安麗等［7］報道了1齡澳洲長鰭鰻Anguilla reinhardtii形態性狀對體重的影響效果。

日本鰻鱺Anguilla japonica隸屬于硬骨魚綱Osteichthyes、鰻鱺目Anguilliformes、鰻鱺科Anguillidae、鰻鱺屬Anguilla，主要分布于日本、臺灣、中國大陸沿海、朝鮮半島等海域，是我國重要的養殖經濟魚類之一。目前，國內有關日本鰻鱺的研究多集中于捕食行為學［8］、繁殖生物學［9-10］、人工養殖技術［11］及不同養殖模式下養成品質的比較［12-14］等方面，而有關日本鰻鱺形態性狀對其體重、凈體重影響效應的研究則迄今尚未見報道?；诖?，本文作者以水庫放養模式養成的日本鰻鱺同生群為研究對象，采用相關分析、通徑分析和偏回歸分析，建立形態性狀對體重和凈體重的多元回歸方程，以期為提高日本鰻鱺親本的選育精度提供科學依據。

1　材料與方法

1.1樣品來源及其形態性狀的測定

本研究實驗用樣品為購自余姚市東超漁業有限公司的4齡日本鰻鱺，養殖模式同王志錚等［12］所述的水庫放養模式。運回實驗室停食馴養1 d后，任選其中無病無傷、體形完整、肛長相近的25尾作為實驗對象，逐一測量其全長（X1）、肛長（X2）、頭長（X3）、吻長（X4）、下頜長（X5）、頭高（X6）、肛部體高（X7）、頭寬（X8）、眼間距（X9）、肛部體寬（X10）、肛圍（X11）、體重（BW）、凈體重（NW）等13項表型性狀，測量方法亦同王志錚等［12］。

1.2分析方法

采用SPSS 19.0軟件計算實驗樣品各表型性狀的均值、標準差和變異系數，并就各形態性狀與體重和凈體重進行相關分析的基礎上，經通徑分析保留那些通徑系數達到顯著水平的形態性狀，進而計算它們對體重、凈體重的決定系數和復相關指數，通過偏回歸分析，建立上述被保留的形態性狀對體重和凈體重的回歸方程，并進行統計學檢驗。

2　結果

2.1表型性狀的參數估計值

由表1可知，在本研究所測各項表型性狀指標中變異系數最大的為體重，達52.41%，凈體重次之，為51.03%，兩者變異系數值為其它形態表型的2.1～9.6倍，表明水庫放養模式下日本鰻鱺養成品的體重和凈體重均具較大的選擇潛力。

表1　日本鰻鱺各形態性狀的統計量描述Tab.1　The apparent statistics of various phenotypic and morphologic traits for Anguilla japonica

2.2相關分析

由表2可知，本研究所涉形態性狀中，除全長（X1）、下頜長（X5）與體重和凈體重的相關系數均僅達到顯著水平（P＜0.05）外，其它形態性狀與體重和凈體重的相關系數均達到極顯著水平（P＜0.01）。其中，與體重、凈體重的相關系數均以肛部體寬為最大，分別為0.974和0.973，與體重、凈體重的相關系數均以下頜長為最小，分別為-0.431和-0.428。

表2　日本鰻鱺各表型性狀間的相關系數（n=25）Tab.2　The phenotype correlation coefficient between the traits of A.japonica（n=25）

2.3通徑分析

由表3可知，全長（X1）、肛長（X2）、下頜長（X5）、頭寬（X8）和肛部體寬（X10）均為與體重和凈體重通徑系數達到顯著水平的形態性狀，它們與體重和凈體重間的相關系數、直接作用和間接作用均依次呈X10＞X8＞X2＞X1＞X5、X10＞X8＞X2＞X1＞X5和X8＞X2＞X10＞X1＞X5。其中，X10對體重和凈體重的直接作用最大，分別為0.545和0.527，X8對體重和凈體重的間接作用最大，分別為0.715和0.699；除肛部體寬（X10）直接作用大于間接作用外，其余性狀直接作用均小于間接作用。綜上可知，肛部體寬（X10）為影響體重和凈體重的核心變量，頭寬（X8）為影響體重和凈體重的重要變量。

表3　日本鰻鱺形態性狀對體重、凈體重影響的通徑分析Tab.3　The path analysis of morphometric traits on body weight and net weight of A.japonica

2.4決定程度分析

由表4可知，單個被保留形態性狀對體重、凈體重的決定系數分別呈X10（0.297 4）＞X8（0.043 0）＞X2（0.041 1）＞X1（0.010 5）＞X5（0.006 4），X10（0.277 4）＞X8（0.051 0）＞X2（0.042 1）＞X1（0.010 5）＞X5（0.005 8），性狀兩兩交互對體重、凈體重的決定系數分別呈X10-X8＞X10-X2＞X8-X2＞X10-X1＞X10-X5＞X1-X2＞X1-X8＞X5-X8＞X5-X2＞X1-X5和X10-X8＞X10-X2＞X8-X2＞X10-X1＞X10-X5＞X1-X8＞X1-X2＞X5-X8＞X5-X2＞X5-X1，即上述被保留形態性狀在對體重、凈體重的決定作用次序上具高度一致性。肛部體寬（X10）不僅在單一性狀對體重、凈體重的決定程度上均遠高于其它被保留的形態性狀，而且X10-X8組合對體重、凈體重的決定程度也均遠高于其它被保留形態性狀組合，再次表明肛部體寬（X10）為決定體重、凈體重的核心變量，頭寬（X8）為決定體重、凈體重的重要變量。

表4　日本鰻鱺形態性狀對體重、凈體重的決定系數Tab.4　Determinant coefficients of morphometric traits on body weight and net weight of A.japonica

2.5建立多元回歸方程

由表5可知，本研究所涉實驗測定群體中經通徑分析后被保留的形態性狀對體重、凈體重的復相關指數分別為0.985和0.986，表明它們是影響體重和凈體重的關鍵形態性狀組合。經偏回歸分析，建立上述形態性狀與體重、凈體重間的多元回歸方程分別為BW=17.669X10+1.785X2+8.292X8+0.153X1-13.900X5-658.897和NW=15.505X10+1.643X2+8.209X8+0.139X1-12.030X5-606.031（表6）；方差分析表明，所建多元回歸方程具有統計學意義（表7）。

表5　日本鰻鱺形態性狀與體重、凈體重的復相關分析Tab.5　The multiple-correlation coefficients of morphometric traits to the weight and net weight of A.japonica

表6　偏回歸系數和回歸常數顯著性檢驗Tab.6　Significance test of partial regression coefficient and regression constants

表7　多元回歸方程的方差分析Tab.7　Analysis of variance of the multiple regression equations for A.japonica

3　討論

3.1影響水庫放養模式日本鰻鱺養成品體重和凈體重的關鍵性狀組合的確定

由表2可見，本研究所涉11項形態性狀與體重、凈體重的相關系數均達到顯著水平（P＜0.05），但經通徑分析將其剖分為直接作用和間接作用后，達到顯著水平的形態性狀僅為肛部體寬（X10）、肛長（X2）、頭寬（X8）、全長（X1）和下頜長（X5）等5項形態性狀（表3），表明簡單相互關系并未能真實反映形態性狀與體重、凈體重間的內在關聯。劉小林等［15］指出，在表型性狀相關分析的基礎上，僅當相關指數R2或決定系數加和值大于或等于0.85時，表明影響因變量的主要自變量已經找到。由表4可見，本研究經通徑分析后被保留的形態性狀對體重、凈體重的決定系數加和值均為0.986，表明本研究被保留的形態性狀為影響日本鰻鱺體重、凈體重的關鍵形態性狀組合。

3.2關于決定水庫放養模式日本鰻鱺養成品體重和凈體重關鍵形態性狀組合的成因分析

一般而言，魚類隨著養殖時間的推移，其個體生長往往需經歷由以體型增長為主的幼魚階段逐漸向以體型增粗為主的養成階段轉變的過程，即與體重和凈重增長最為密切的變量往往首先是以體長為代表的橫向維度性狀，故本研究中全長（X1）和肛長（X2）被選入影響日本鰻鱺養成品體重和凈重增長的形態性狀應屬情理之中。泥鰍與日本鰻鱺體型均屬棍棒型。已有研究表明，決定4月齡泥鰍體重的關鍵形態性狀為全長［5］，而決定1年齡大鱗副泥鰍體重的關鍵形態性狀則為體厚［6］。這一情形與安麗等［7］報道的決定1年齡澳洲長鰭鰻體重的關鍵形態性狀為體長，而決定本研究所涉4年齡日本鰻鱺養成品體重的關鍵形態性狀則為肛部體寬（X10）的結果相近，表明本研究所涉日本鰻鱺養成品已處于以增寬為主的個體生長階段，故頭寬（X8）順勢成為在決定實驗對象體重上僅次于肛部體寬（X10）的重要變量，這也與體型較寬有利于營養物質積貯而體重相應較重的實際生產經驗相一致［3-4］。另，口裂寬［7］和下頜長分別被選入影響1年齡澳洲長鰭鰻和本研究所涉日本鰻鱺養成品體重增長的形態性狀的情形，無疑表征了攝食器官幾何空間大小在影響鰻鱺屬種類體重增長上的重要性。

同期研究還發現，與本研究肛長相仿的池塘專養模式和青蝦池塘套養模式養成的日本鰻鱺，其體重和凈體重與本研究所涉形態性狀間的相關指數均遠低于0.85，即無法從本研究所涉形態性狀中找尋到決定其體重和凈體重的關鍵形態性狀組合。究其原因，可能與池塘專養模式和青蝦池塘套養模式養成的日本鰻鱺因長期處于食物保障程度相對較高狀態，致使其體型均質化水平遠高于食物保障程度長期處于匱乏狀態的水庫放養模式，這一情形也側面印證了水庫放養模式日本鰻鱺養成品在形質［12］，血清生化指標和臟器消化酶、抗氧化酶活力［13］，以及體色和肌肉品質［14］上與池塘專養模式和青蝦池塘套養模式養成的日本鰻鱺具顯著差異的事實。

［1］嚴福升，王志剛，劉旭東，等.3月齡牙鲆形態性狀對體質量的通徑分析［J］.漁業科學進展，2010，31（2）:45-50.

［2］李璽洋，白俊杰，樊佳佳，等.二齡草魚形態性狀對體質量影響效果的分析［J］.上海海洋大學學報，2012，21（4）:535-541.

［3］霍堂斌，袁美云，張麗榮，等.黑斑狗魚形態性狀對體重的影響效果分析［J］.水產學雜志，2011，24（2）:26-30.

［4］王新安，馬愛軍，許可，等.大菱鲆幼魚表型形態性狀與體重之間的關系［J］.動物學報，2008，54（3）:540-545.

［5］王磊，凌去非，郝小鳳，等.4月齡泥鰍形態性狀對體重的影響分析［J］.水生態學雜志，2011，32（4）:97-102.

［6］張樊苗，呂鳳瓊，王磊，等.1齡大鱗副泥鰍形態性狀對體重的影響［J］.安徽農業科學，2012，40（12）:7 182-7 185.

［7］安麗，孟慶磊，董學颯，等.澳洲長鰭鰻各形態性狀對體重的影響效果分析［J］.中國農學通報，2012，28（2）:60-64.

［8］王志錚，楊磊，朱衛東，等.日本鰻鱺（Anguilla japonica）對日本沼蝦（Macrobrachium nipponense）的捕食效應［J］.海洋與湖沼，2011，42（1）:107-113.

［9］王廣軍，謝駿，余德光，等.鰻鱺繁殖生物學研究進展［J］.南方水產，2005，1（1）:71-75.

［10］張潔明，柳凌，郭峰，等.人工誘導日本鰻鱺性腺發育組織學研究［J］.中國水產科學，2007，14（4）:593-601.

［11］黃大明，陳世群.鰻鱺的生活史和人工育苗技術探討［J］.動物學雜志，1997，32（3）:39-48.

［12］王志錚，楊磊，朱衛東，等.三種養殖模式下日本鰻鱺養成品的形質差異［J］.應用生態學報，2012，23（5）:1 385-1 392.

［13］王志錚，趙晶，楊磊，等.三種養殖模式下日本鰻鱺（Anguilla japonica）養成品血清生化指標和臟器消化酶、抗氧化酶活力的差異［J］.海洋與湖沼，2013，44（2）:403-408.

［14］王志錚，付英杰，楊磊，等.三種養殖模式下日本鰻鱺（Anguilla japonica）養成品體色和肌肉品質的差異［J］.海洋與湖沼，2013，44（4）:1 042-1 049.

［15］劉小林，吳長功，張志懷，等.凡納對蝦形態性狀對體重的影響效果分析［J］.生態學報，2004，24（4）:857-863.

Influence Effect of Morphometric Traits on Body Weight and Net Weight of Adult Anguilla japonica from Reservoir in Stocking Model

SI Lie-gang1，2，YUAN Xiang-yang1，YANG Lei1，et al
（1.Fishery School of Zhejiang Ocean University，Zhoushan316022；2.Marine and Fishery Research Institute of Ningbo，Ningbo315012，China）

Influence effect of morphometric traits on body weight and net weight of adult Anguilla japonica was studied by using correlation analysis，path analysis and regression analysis.25 A.japonica were randomly sampled for experiment.Total length（TL），snout-vent length，head length，snout length，mandible length，head height，body depth，head width，interorbital space，rectal width at withers，periproct，body weight and netweight were measured for all individuals.The results showed that:（1）All correlation coefficients between phenotyoic traits and body weight，net weight were at significant level（P＜0.05）；（2）The correlation indexes of the five selected morphometric traits to body weight and net weight were 0.985 and 0.986，respectively.The utmost direct actions of rectal width at withers to body weight and net weight were 0.545 and 0.527，respectively.The utmost indirect actions of head width to body weight and net weight were 0.715 and 0.699，respectively；（3）Multiple regression equation used to estimate body weight，and net weight of adult A.japonica were established，respectively.

morphometric traits；body weight；net weight；Anguilla japonica；stocking model in reservior

S965.199

A

1008-830X（2015）02-0119-06

2014-12-20

國家科技部星火計劃項目（2013GA701022）；浙江省重大科技專項農業重點項目（2008C12083）；寧波市農業創新創業項目（2013C920017）

斯烈鋼（1976-），男，浙江諸暨人，高級工程師，研究方向：水域生態評價與水產品安全養殖技術.

王志錚，教授，研究員.E-mail:wzz_1225@163.com