平洪領(lǐng)　王天明　呂振明等

摘要：選取40日齡曼氏無針烏賊100只，測定了胴長（X1）、全長（X2）、眼間距（X3）、體寬（X4）、觸腕長（X5）、體質(zhì)量（Y）等6項(xiàng)性狀指標(biāo)。采用相關(guān)分析和通徑分析方法，計算了以測量性狀為自變量對因變量體質(zhì)量的相關(guān)系數(shù)、通徑系數(shù)、決定系數(shù)、相關(guān)指數(shù)，運(yùn)用逐步剔除方法后，建立測量性狀對體質(zhì)量的回歸方程。結(jié)果表明：曼氏無針烏賊的5個體型性狀均與體質(zhì)量呈極顯著相關(guān)；胴長、體寬對體質(zhì)量的通徑系數(shù)達(dá)到顯著水平，是直接影響體質(zhì)量的重要性狀；決定系數(shù)分析結(jié)果與通徑分析結(jié)果相同；采用逐步線性回歸的方法建立了以形態(tài)性狀指標(biāo)值估計體質(zhì)量的回歸方程：Y=-78.44+13.98X1+6.21X4；選取的性狀與體質(zhì)量的復(fù)相關(guān)指數(shù)為0.895，定量分析了形態(tài)性狀對體質(zhì)量的影響效果。

關(guān)鍵詞：曼氏無針烏賊；形態(tài)性狀；通徑分析；多元回歸方程

中圖分類號： S917.4 文獻(xiàn)標(biāo)志碼： A

文章編號：1002-1302（2015）03-0218-03

曼氏無針烏賊（Sepiella japonica）屬軟體動物門（Mollusca）頭足綱（Cephalopoda）烏賊目（Sepioidea）烏賊科（Sepiidae）無針烏賊屬（Sepiella），為中型烏賊。曼氏無針烏賊是我國重要的水產(chǎn)經(jīng)濟(jì)動物，廣泛分布于中國沿海、俄羅斯遠(yuǎn)東海、日本瀨戶內(nèi)海、朝鮮西南沿海等海區(qū)，是我國東海漁場傳統(tǒng)的“四大海產(chǎn)”之一。歷史上浙江省曼氏無針烏賊最高年產(chǎn)量達(dá)6萬t，占該省海洋捕撈總產(chǎn)量的9.3%[1]。

國內(nèi)外學(xué)者已對曼氏無針烏賊的生物學(xué)、繁殖習(xí)性、洄游分布、資源狀況以及增養(yǎng)殖開展過研究[2-5]。形態(tài)性狀和體質(zhì)量指標(biāo)是遺傳育種與種苗繁殖研究的重要依據(jù)。自20世紀(jì)20年代Wright首次提出通徑系數(shù)的分析方法[6]以來，多元分析法已被廣泛應(yīng)用于水產(chǎn)養(yǎng)殖育種方案的優(yōu)化和生產(chǎn)量的估計[7-8]，例如劉小林等對櫛孔扇貝的體質(zhì)量與多個測量性狀進(jìn)行了多元分析[9]；平洪領(lǐng)等利用逐步線性回歸法實(shí)現(xiàn)天津厚蟹表型性狀與體質(zhì)量的通徑分析[10]。隨著曼氏無針烏賊野生資源遭到極大的破壞[11]，人工養(yǎng)殖規(guī)模的擴(kuò)大以及國家對曼氏無針烏賊放流投入的加大，對曼氏無針烏賊的人工選育工作提出了迫切需求。本研究用曼氏無針烏賊的體質(zhì)量與測量性狀進(jìn)行多元分析，利用相關(guān)分析、通徑分析、回歸分析方法，給出了影響曼氏無針烏賊體質(zhì)量的主要測量性狀及其直接作用和間接作用的大小，建立了估計曼氏無針烏賊體質(zhì)量的多元回歸方程，以期為曼氏無針烏賊的選育工作提供理論依據(jù)。

1 材料與方法

1.1 供試材料

供試曼氏無針烏賊來源于福建省漳州市自然海區(qū)的野生受精卵，均由浙江大海洋科技有限公司培育。隨機(jī)選取40日齡曼氏無針烏賊100只，分別測定其胴長（X1）、全長（X2）、眼間距（X3）、體寬（X4）、觸腕長（X5）、體質(zhì)量（Y）等6項(xiàng)指標(biāo)。

1.2 測定方法

采用游標(biāo)卡尺測量曼氏無針烏賊胴長、全長、眼間距、體寬、觸腕長等5項(xiàng)體型性狀的指標(biāo)（圖1），精確到0.01 cm。用電子天平稱量曼氏無針烏賊體質(zhì)量，精確到0.01 g。

1.3 數(shù)據(jù)統(tǒng)計與分析

應(yīng)用SPSS 17.0、Excel 2010軟件對曼氏無針烏賊形態(tài)性狀、體質(zhì)量數(shù)據(jù)進(jìn)行初步統(tǒng)計處理，利用各表型參數(shù)統(tǒng)計量，分別進(jìn)行表型相關(guān)分析、測量性狀各指標(biāo)對體質(zhì)量的通徑分析和決定系數(shù)的計算。依據(jù)各性狀對體質(zhì)量的直接作用和間接影響建立回歸方程。運(yùn)用SPSS 17.0軟件只需“Analyze→Regression→Linear”程序便可獲得通徑系數(shù)、相關(guān)系數(shù)以及顯著性檢驗(yàn)等信息[12]。利用SPSS 17.0軟件的獨(dú)立樣t檢驗(yàn)對體質(zhì)量進(jìn)行差異顯著性檢驗(yàn)，置信區(qū)間設(shè)為95%。

2 結(jié)果與分析

2.1 曼氏無針烏賊各性狀參數(shù)統(tǒng)計

由表1可知，所測的5個形態(tài)指標(biāo)中，觸腕長最大，為13.533 cm；眼間距最小，為3.198 cm。體質(zhì)量變異程度最高，為21.624%，其次依次為體寬>觸腕長>全長>胴長>眼間距。形態(tài)指標(biāo)間的變異范圍在7.974%～12.551%，表明所選擇的樣本均一性較高。

2.2 曼氏無針烏賊性狀間的相關(guān)分析

由表2可知，曼氏無針烏賊胴長、全長、體寬、眼間距、觸腕長、體質(zhì)量之間的相關(guān)性均達(dá)到了極顯著水平（P<0.01），表明將所選指標(biāo)用來進(jìn)行相關(guān)分析具有可靠性。其中體質(zhì)量與胴長的相關(guān)系數(shù)最大，為0.883；眼間距與體寬的相關(guān)系數(shù)最小，為0.460；體質(zhì)量與其他各性狀間均呈極顯著正相關(guān)，與其他性狀的相關(guān)程度大小依次為胴長>眼間距>全長>觸腕長>體寬。

2.4 測量性狀對體質(zhì)量的通徑分析

為了揭示各性狀間的本質(zhì)聯(lián)系，應(yīng)對各性狀間的相關(guān)系數(shù)進(jìn)行剖析，即找出從原因到結(jié)果的直接影響和通過其他原因到該結(jié)果的間接影響，正是在間接影響的作用下，使得相關(guān)系數(shù)不能直接體現(xiàn)變量間的相關(guān)關(guān)系。由表7可見，曼氏無針烏賊胴長、體寬對體質(zhì)量的通徑系數(shù)達(dá)到了極顯著水平，它們的通徑系數(shù)分別是0.671、0.276，可以剔除其他通徑系數(shù)不顯著的性狀。經(jīng)顯著性檢驗(yàn)，由表5可知，保留的2個形態(tài)性狀胴長、體寬都達(dá)到了顯著水平（P<0.05）。其中對體質(zhì)量直接作用的大小關(guān)系為胴長>體寬；間接作用的大小關(guān)系為體寬>胴長（表7）。因此，胴長對體質(zhì)量的直接作用較大，而間接作用較小， 說明胴長主要是通過直接作用對體質(zhì)量產(chǎn)生影響，而體寬則主要是通過間接作用對體質(zhì)量產(chǎn)生影響。

2.5 各測量性狀對體質(zhì)量的決定程度分析

單個性狀對體質(zhì)量的決定系數(shù)di=P2iy，2個性狀對體質(zhì)量的共同決定系數(shù)dij=2rijPiyPjy，由此計算出形態(tài)性狀間協(xié)同對體質(zhì)量的決定系數(shù)。由表8可見，胴長、體寬單獨(dú)對體質(zhì)量的決定系數(shù)分別為0.450、0.076，0.370為胴長與體寬共同對體質(zhì)量的決定系數(shù)。決定系數(shù)總和∑d=0.896，這與表4中決定系數(shù)（r2）相同。表明本研究所列形態(tài)性狀是影響體質(zhì)量的重點(diǎn)性狀，其他性狀影響相對較小。endprint

3 結(jié)論與討論

3.1 相關(guān)分析與回歸分析的聯(lián)系

從相關(guān)性分析結(jié)果可以看出，盡管所測的5個形態(tài)性狀均與體質(zhì)量呈極顯著相關(guān)，但經(jīng)過逐步引入剔除法（Stepwise）進(jìn)行回歸分析時，大部分性狀指標(biāo)被剔除掉，只保留了胴長和體寬。這說明即便相關(guān)系數(shù)很大，回歸關(guān)系也不一定會顯著。這可能是因?yàn)檫@些被剔除掉的形態(tài)參數(shù)與體質(zhì)量的回歸關(guān)系不夠顯著，而且它們分別與入選的自變量存在很強(qiáng)的相關(guān)性，不能在回歸方程中共存[13]。

3.2 影響體質(zhì)量重點(diǎn)性狀的確定

進(jìn)行通徑分析和決定系數(shù)分析時，只有當(dāng)決定指數(shù)（r2）或各自變量對因變量的單獨(dú)決定系數(shù)及兩兩共同決定系數(shù)的總和在數(shù)值上≥0.85時，才表明依變量的主要自變量已找到。本研究得到r2=0.896，說明所保留的曼氏無針烏賊的形態(tài)性狀是影響體質(zhì)量的重點(diǎn)性狀，進(jìn)一步說明通徑系數(shù)分析結(jié)果能夠反映形態(tài)性狀與體質(zhì)量之間的真實(shí)關(guān)系。在選擇育種中，胴長、體寬是理想的測度指標(biāo)。因此，在曼氏無針烏賊生長快速品種的選育中，在以體質(zhì)量為主要選育性狀的前提下，可以考慮以胴長、體寬等性狀作為輔助選擇性狀，綜合考慮以加快生長快的新品種選育。

3.3 曼氏無針烏賊數(shù)量性狀選擇育種的意義

運(yùn)用通徑分析和多元分析對魚類[14]、貝類[15]、蝦蟹類[16-17]進(jìn)行數(shù)量性狀選育的報道較多，但是對曼氏無針烏賊進(jìn)行性狀選育的報道較少。本研究對所測性狀數(shù)據(jù)進(jìn)行生物學(xué)分析統(tǒng)計后發(fā)現(xiàn)，體質(zhì)量的變異系數(shù)最大，一方面說明對曼氏無針烏賊的體質(zhì)量性狀進(jìn)行選育潛力很大，另一方面又說明依據(jù)體質(zhì)量對曼氏無針烏賊直接選育，可能會因?yàn)榄h(huán)境因素的改變而產(chǎn)生較大的系統(tǒng)誤差。而胴長與體寬的變異系數(shù)相對較小，因此可以根據(jù)以上形態(tài)性狀對體質(zhì)量進(jìn)行間接選育，這樣可以減少環(huán)境所產(chǎn)生的影響，確保選育的真實(shí)性與準(zhǔn)確性。本研究通過通徑分析和多元分析確定了影響體質(zhì)量的最優(yōu)回歸方程：Y=-78.44+13.98X1+6.21X3。因此，在曼氏無針烏賊的選育工作中，可以直接測量胴長與體寬，不需要測量其他形態(tài)性狀，這樣既保證了選育的效果，又可以減少工作量。

參考文獻(xiàn)：

[1]劉雅丹. 2001年世界頭足類市場的特點(diǎn)[J]. 中國漁業(yè)經(jīng)濟(jì)，2002（2）：50.

[2]張 炯，盧偉成. 曼氏無針烏賊繁殖習(xí)性的初步觀察[J]. 水產(chǎn)學(xué)報，1965，2（2）：35-43.

[3]唐逸民，吳常文. 曼氏無針烏賊生物學(xué)特性及漁場分布變化[J]. 浙江水產(chǎn)學(xué)院學(xué)報，1986，5（2）：165-170.

[4]Moltschaniwskyj N A，Jackson G D. Growth and tissue composition as a function of feeding history in juvenile cephalopods[J]. Journal of Experimental Marine Biology and Ecology，2000，253（2）：229-241.

[5]Correia M，Domingues P M，Sykes A，et al. Effects of culture density on growth and broodstock management of the cuttlefish，Sepia officinalis（Linnaeus，1758）[J]. Aquaculture，2005，245（1/4）：163-173.

[6]Wright S. Correlation and causation[J]. Journal of Agriculture Research，1921，20（7）：557-585.

[7]Deboski P，Dobosz S，Robak S，et al. Fat level in body of juvenile Atlantic salmon（Salmo salar L.），and sea trout（Salmo trutta m. trutta L.），and method of estimation from morphometric data[J]. Archives of Polish Fisheries，1999，7（2）：237-243.

[8]高保全，劉 萍，李 健，等. 三疣梭子蟹形態(tài)性狀對體重影響的分析[J]. 海洋水產(chǎn)研究，2008，29（1）：44-50.

[9]劉小林，常亞青，相建海，等. 櫛孔扇貝殼尺寸性狀對活體重的影響效果分析[J]. 海洋與湖沼，2002，33（6）：673-678.

[10]平洪領(lǐng)，李玉全. 逐步線性回歸法實(shí)現(xiàn)天津厚蟹（Helice tientsinensis）表型性狀與體重的通徑分析[J]. 海洋與湖沼，2013，44（5）：1353-1357.

[11]吳常文，趙淑江，徐蝶娜.舟山漁場針烏賊（Sepia andreana）的生物學(xué)特性及其漁場分布變遷[J]. 海洋與湖沼，2006，37（3）：231-237.

[12]杜家菊，陳志偉. 使用SPSS線性回歸實(shí)現(xiàn)通徑分析的方法[J]. 生物學(xué)通報，2010，45（2）：4-6.

[13]于 飛，張慶文，孔 杰，等. 大菱鲆測量性狀對體重的影響效果分析[J]. 海洋水產(chǎn)研究，2008，06（6）：33-39.

[14]Kora H，Tsuchimoto M，Miyata K，et al. Estimation of body fat content from standard body length and body weight on cultured Red Sea bream[J]. Fisheries Science，2000，66（2）：365-371.

[15]Yang H S，Zhang T，Wang J，et al. Growth characteristics of Chlamys farreri and its relation with environmental factors in intensive raft-culture areas of Sishiliwan Bay，Yantai[J]. Journal Journal of Shellfish Research，1999，18（1）：71-76.

[16]鄧平平，戴習(xí)林，臧維玲，等. 羅氏沼蝦形態(tài)性狀對體重的影響[J]. 江蘇農(nóng)業(yè)科學(xué)，2012，40（7）：205-211.

[17]Thomas M M. Age and growth，length weight relation and relative condition factor of Penaeus semiculcatus[J]. Fish，1975，22：133-142.endprint