水库放养模式下日本鳗鲡养成品形态性状对体重和净体重的影响效应

斯烈钢，袁向阳，杨　磊，朱卫东，，申屠琰，申屠基康，王志铮

（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