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

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

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

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

文章編號(hào)：1002-1302（2015）03-0218-03

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

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

1 材料與方法

1.1 供試材料

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

1.2 測(cè)定方法

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

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

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

2 結(jié)果與分析

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

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

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

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

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

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

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

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

3 結(jié)論與討論

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

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

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

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

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

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

參考文獻(xiàn)：

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

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

[3]唐逸民，吳常文. 曼氏無(wú)針烏賊生物學(xué)特性及漁場(chǎng)分布變化[J]. 浙江水產(chǎn)學(xué)院學(xué)報(bào)，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)性狀對(duì)體重影響的分析[J]. 海洋水產(chǎn)研究，2008，29（1）：44-50.

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

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

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

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

[13]于 飛，張慶文，孔 杰，等. 大菱鲆測(cè)量性狀對(duì)體重的影響效果分析[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)性狀對(duì)體重的影響[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