河川沙塘鱧表型性狀對(duì)體質(zhì)量影響的通徑分析

徐宇+史楊白+鄭嘯宇+霍春林

摘要： 測(cè)定了500尾河川沙塘鱧（Odontobutis potamophila）體質(zhì)量、體長(zhǎng)、體高、頭長(zhǎng)、尾柄長(zhǎng)、尾柄高、頭寬、吻長(zhǎng)等8個(gè)性狀，計(jì)算各性狀間的相關(guān)系數(shù)，采用通徑分析方法，以形態(tài)性狀為自變量，體質(zhì)量為依變量的通徑系數(shù)決定系數(shù)，并對(duì)各形態(tài)性狀與體質(zhì)量的相關(guān)性進(jìn)行分析。結(jié)果顯示，各形態(tài)性狀與體質(zhì)量之間的相關(guān)性均達(dá)到極顯著水平（P<0.01）；體長(zhǎng)對(duì)體質(zhì)量的直接作用最大（0.511），剔除通徑系數(shù)檢驗(yàn)中不顯著的自變量，利用逐步回歸的方法，建立以體長(zhǎng)、體高、頭長(zhǎng)、頭寬、尾柄高為自變量估計(jì)體質(zhì)量的最優(yōu)多元線性回歸方程為：y=3.773x1+3.997x2+4.118x3+2.206x4+4.431x5-48.676。結(jié)果表明，體長(zhǎng)、體高、頭長(zhǎng)、頭寬、尾柄高等5個(gè)形態(tài)性狀對(duì)體質(zhì)量直接作用達(dá)極顯著水平。

關(guān)鍵詞： 河川沙塘鱧；形態(tài)性狀；相關(guān)分析；通徑分析；多元回歸

中圖分類號(hào)： S961.2 文獻(xiàn)標(biāo)志碼： A 文章編號(hào)：1002-1302（2016）03-0262-03

河川沙塘鱧（Odontobutis potamophila）隸屬于鱸形目，是我國(guó)特有的品種，主要分布于長(zhǎng)江中下游及沿江各支流、錢(qián)塘江水系、閩江水系，偶見(jiàn)于黃河水系[1]。該魚(yú)個(gè)體小，自然資源量少，但含肉量高、肉味鮮美，深受人們的喜愛(ài)，其價(jià)格更是一路攀升，是極具發(fā)展?jié)摿Φ乃a(chǎn)養(yǎng)殖新品種。形態(tài)性狀、體質(zhì)量是遺傳育種與種苗繁育研究的重要指標(biāo)。在實(shí)際生產(chǎn)中，活體體質(zhì)量稱量不如形態(tài)指標(biāo)簡(jiǎn)單準(zhǔn)確。因此，分析了解各形態(tài)性狀之間的相關(guān)性以及與產(chǎn)量之間的關(guān)系，可以為遺傳育種提供理論依據(jù)。近年來(lái)，利用相關(guān)分析、通徑分析、多元回歸分析等評(píng)價(jià)生物形態(tài)性狀與體質(zhì)量的關(guān)系已被廣泛應(yīng)用于水產(chǎn)生物如魚(yú)[2-7]、蝦[8-9]、蟹[10]、貝類[11-12] 中。關(guān)于鴨綠沙塘鱧選育研究已有報(bào)道[13]。本研究選取河川沙塘鱧為試驗(yàn)材料，運(yùn)用回歸分析、相關(guān)分析、通徑分析方法，評(píng)估影響河川沙塘鱧體質(zhì)量的主要形態(tài)性狀，建立估計(jì)體質(zhì)量的最優(yōu)回歸方程，以期為河川沙塘鱧選育工作提供理論依據(jù)。

1 材料與方法

1.1 材料

2014年7月20日從江蘇省淡水水產(chǎn)研究所祿口基地隨機(jī)選取河川沙塘鱧500尾，體質(zhì)量12.80～69.60 g，體長(zhǎng) 839～18.18 cm。

1.2 測(cè)量方法

8個(gè)形態(tài)性狀（體長(zhǎng)、體高、尾柄高、尾柄長(zhǎng)、頭長(zhǎng)、頭寬、吻長(zhǎng)、體質(zhì)量）定義參考文獻(xiàn)[14]。使用游標(biāo)卡尺測(cè)量各長(zhǎng)度性狀，精確到0.01 cm；用電子天平測(cè)量體質(zhì)量，精確到 0.01 g。

1.3 數(shù)據(jù)分析

利用數(shù)據(jù)統(tǒng)計(jì)分析軟件SPSS19.0中的K-S單樣本檢驗(yàn)（Kolomogorov-Simirnov one sample test） 檢驗(yàn)各性狀是否符合正態(tài)分布。計(jì)算各形態(tài)性狀的均值、標(biāo)準(zhǔn)差、變異系數(shù)，并在表型性狀相關(guān)分析以及形態(tài)性狀對(duì)總體質(zhì)量通徑分析的基礎(chǔ)上，分析其通徑系數(shù)達(dá)到顯著水平的形態(tài)性狀對(duì)體質(zhì)量的直接作用與間接作用，計(jì)算決定系數(shù)，確定影響體質(zhì)量的關(guān)鍵形態(tài)性狀；采用逐步線性回歸法剔除不顯著的表型性狀，最終建立表型性狀對(duì)質(zhì)量性狀的多元回歸方程。參照杜家菊等的方法[15]，用SPSS 19.0軟件分析數(shù)據(jù)。

2 結(jié)果與分析

2.1 各性狀的參數(shù)

河川沙塘鱧8個(gè)性狀的表型數(shù)據(jù)見(jiàn)表1。單樣本K-S檢驗(yàn)結(jié)果顯示，P值均在0.05以上，說(shuō)明本研究所分析性狀的分布均沒(méi)有顯著偏離正態(tài)分布，可以對(duì)這些性狀進(jìn)行進(jìn)一步的通徑分析。性狀變異系數(shù)越大，表明該性狀的選擇潛力越大[3]。由表1可知，體質(zhì)量的變異系數(shù)最大，說(shuō)明其具有較大的選擇潛力，其次是頭寬、尾柄長(zhǎng)，體長(zhǎng)、頭長(zhǎng)、體高的變異系數(shù)較小，選擇潛力較弱，其余性狀變異系數(shù)中等。

2.2 性狀間表型相關(guān)系數(shù)

對(duì)河川沙塘鱧8個(gè)性狀進(jìn)行表型相關(guān)分析，獲得各性狀間的表型相關(guān)系數(shù)（表2）。結(jié)果顯示，各性狀間均呈極顯著相關(guān)（P<0.01），表明將所選的指標(biāo)用來(lái)進(jìn)行相關(guān)分析具有實(shí)際意義。其中體長(zhǎng)與體質(zhì)量的相關(guān)系數(shù)最大，為0.851，頭長(zhǎng)與頭寬的相關(guān)系數(shù)最小，為0.321，各形態(tài)性狀與體質(zhì)量相關(guān)系數(shù)大小依次為：體長(zhǎng)>頭長(zhǎng)>尾柄高>體高>頭寬>吻長(zhǎng)>尾柄長(zhǎng)。

2.3 各性狀對(duì)體質(zhì)量的影響

在表型相關(guān)基礎(chǔ)上，利用SPSS 19.0軟件計(jì)算河川沙塘鱧各形態(tài)性狀對(duì)體質(zhì)量的通徑系數(shù)，經(jīng)顯著性檢驗(yàn)，保留顯著性較高的體長(zhǎng)、體高、頭長(zhǎng)、頭寬、尾柄高5個(gè)性狀（表3），結(jié)果顯示，體長(zhǎng)對(duì)體質(zhì)量的直接作用（0.511）最大，且大于其間接作用（0.342），是影響體質(zhì)量的主要因素，其他4個(gè)表型性狀對(duì)體質(zhì)量的直接作用均小于間接作用，主要通過(guò)體長(zhǎng)對(duì)體質(zhì)量間接產(chǎn)生影響。

2.4 各性狀對(duì)體質(zhì)量的決定程度

表4表明，各性狀對(duì)體質(zhì)量單獨(dú)的決定系數(shù)中體長(zhǎng)的決定程度最大，為0.261，尾柄高的決定程度最小，為0.008；兩兩共同決定系數(shù)中，體長(zhǎng)和頭長(zhǎng)的共同決定程度最大，為0171，頭寬和尾柄高的共同決定程度最小，為0.013。單獨(dú)的決定系數(shù)和兩兩共同決定系數(shù)的總和等于總的決定系數(shù)，5個(gè)性狀對(duì)體質(zhì)量的總決定系數(shù)為0.817。

2.5 多元回歸方程

根據(jù)SPSS線性多元回歸中的逐步回歸分析的結(jié)果（表5），剔除通徑系數(shù)檢驗(yàn)不顯著的變量，對(duì)剩余變量建立以體質(zhì)量為依變量的最優(yōu)線性回歸方程：

y=3.773x1+3.997x2+4.118x3+2.206x4+4.431x5-48.676。

式中：y為體質(zhì)量（g），x1、 x2、 x3、x4、 x5分別為體長(zhǎng)（cm）、體高（cm）、頭長(zhǎng)（cm）、頭寬（cm）、尾柄高（cm）。多元回歸關(guān)系和各個(gè)偏回歸系數(shù)的顯著性檢驗(yàn)結(jié)果表明，回歸關(guān)系達(dá)到極顯著水平（P<0.01）（表6）。

3 結(jié)論與討論

3.1 影響河川沙塘鱧體質(zhì)量主要性狀的確定

本研究結(jié)果表明，河川沙塘鱧形態(tài)性狀與體質(zhì)量的表型相關(guān)系數(shù)均達(dá)到極顯著水平，但由于表型相關(guān)系數(shù)中并未剔除其他性狀的間接影響，不能準(zhǔn)確反映各形態(tài)性狀與體質(zhì)量之間的關(guān)系，因此，利用通徑分析和多元線性回歸分析進(jìn)一步探討各形態(tài)性狀與體質(zhì)量的真實(shí)關(guān)系。結(jié)果顯示，河川沙塘鱧的體長(zhǎng)、體高、頭長(zhǎng)、頭寬、尾柄高5個(gè)形態(tài)性狀對(duì)體質(zhì)量的直接作用達(dá)極顯著水平，其中體長(zhǎng)對(duì)體質(zhì)量的決定系數(shù)（0.261）最大。在通徑分析和決定程度分析時(shí)，只有當(dāng)各形態(tài)性狀對(duì)體質(zhì)量的單獨(dú)決定系數(shù)及兩兩共同決定系數(shù)的總和大于或等于0.85時(shí)才能確定影響體質(zhì)量的主要形態(tài)性狀均已找到[3]。本研究中，系數(shù)總和為0.817，略低于0.85，表明本研究所選擇的性狀還不夠全面，原因可能是不同體型的魚(yú)類對(duì)體質(zhì)量造成影響的性狀也會(huì)有所差異，如大口黑鱸的體寬、體長(zhǎng)、眼間距是影響體質(zhì)量的主要形態(tài)性狀[16]。沙塘鱧同為鱸形目，是小型魚(yú)類，體型敦實(shí)，后續(xù)研究應(yīng)增加體寬、眼間距等表型數(shù)據(jù)進(jìn)行通徑分析。此外，試驗(yàn)動(dòng)物月齡的差異也是可能原因之一，目前已有研究證實(shí)，不同月齡水產(chǎn)動(dòng)物的性狀主成分有所不同，且影響其體質(zhì)量的主要形態(tài)性狀也有所不同[17-19]。對(duì)與河川沙塘鱧體型類似的鴨綠沙塘鱧的形態(tài)性狀與體質(zhì)量進(jìn)行通徑分析表明，全長(zhǎng)、體長(zhǎng)、體高、尾柄高是影響體質(zhì)量的主要性狀[13]。因此本研究獲得的這5個(gè)形態(tài)性狀也可作為該階段的理想測(cè)度選育指標(biāo)。

3.2 對(duì)選育的指導(dǎo)意義

在水產(chǎn)養(yǎng)殖方案優(yōu)化和生產(chǎn)量估計(jì)中常用到相關(guān)分析、通徑分析、多元回歸分析[20]。魚(yú)類生長(zhǎng)性狀的選育中，一般以體長(zhǎng)、體質(zhì)量作為選育指標(biāo)。研究表明，體長(zhǎng)、體高等形態(tài)性狀對(duì)體質(zhì)量的增長(zhǎng)具有重要的影響，因此在河川沙塘鱧選育中應(yīng)對(duì)體質(zhì)量以及對(duì)與體質(zhì)量具有決定作用的形態(tài)性狀進(jìn)行間接選擇[21]。此外，在后續(xù)研究中，也可考慮將體寬、眼間距等性狀作為重要測(cè)量指標(biāo)。

參考文獻(xiàn)：

[1]伍漢霖. 中國(guó)沙塘鱧屬魚(yú)類的整理和新種的敘述[J]. 上海水產(chǎn)大學(xué)學(xué)報(bào)，1993，2（1）：52-58.

[2]楊慧贊，林 勇，唐章生，等. 吉富羅非魚(yú)生長(zhǎng)性狀的相關(guān)與通徑分析[J]. 華北農(nóng)學(xué)報(bào)，2011（增刊）：264-268.

[3]曾 蘭，林 勇，張永德，等. 奧利亞羅非魚(yú)形態(tài)性狀與體重的通徑分析[J]. 西南農(nóng)業(yè)學(xué)報(bào)，2012，25（1）：295-301.

[4]嚴(yán)福升，王志剛，劉旭東，等. 3月齡牙鲆形態(tài)性狀對(duì)體質(zhì)量的通徑分析[J]. 漁業(yè)科學(xué)進(jìn)展，2010，31（2）：45-50.

[5]佟雪紅，袁新華，董在杰，等. 建鯉自交及與黃河鯉正反雜交子代的生長(zhǎng)比較和通徑分析[J]. 水產(chǎn)學(xué)報(bào)，2008，32（2）：182-189.

[6]佟廣香，匡友誼，許凌雪，等. 哲羅鮭形態(tài)性狀與體重的相關(guān)性分析[J]. 水產(chǎn)學(xué)雜志，2011，24（2）：31-36.

[7]馮冰冰，梅肖樂(lè)，余 寧，等. 龍睛金魚(yú)表型性狀與體質(zhì)量的相關(guān)與通徑分析[J]. 江蘇農(nóng)業(yè)科學(xué)，2013，41（11）：255-257.

[8]張成松，李富花，相建海.脊尾白蝦形態(tài)性狀對(duì)體質(zhì)量影響的通徑分析[J]. 水產(chǎn)學(xué)報(bào)，2013，37（6）：809-815.

[9]楊長(zhǎng)明，何 銅，劉小林，等. 凡納對(duì)蝦形態(tài)性狀對(duì)體質(zhì)量的逐步回歸分析[J]. 西北農(nóng)業(yè)學(xué)報(bào)，2011，20（2）：15-20.

[10]耿緒云，王雪惠，孫金生，等. 中華絨螯蟹（Eriocheir sinensis）一齡幼蟹外部形態(tài)性狀對(duì)體重的影響效果分析[J]. 海洋與湖沼，2007，38（1）：49-54.

[11]蔣壽佳，張建設(shè)，陳永久，等. 岱山等邊淺蛤形態(tài)與重量性狀參數(shù)的相關(guān)性及通徑分析[J]. 安徽農(nóng)業(yè)科學(xué)，2015，43（5）：122-124.

[12]陸雅鳳，趙 晟，徐梅英，等. 東極厚殼貽貝養(yǎng)殖群體表型性狀的相關(guān)與通徑分析[J]. 安徽農(nóng)業(yè)科學(xué)，2015，43（3）：147-150.

[13]杜 華，李 軍，宋文華. 鴨綠沙塘鱧形態(tài)性狀對(duì)體質(zhì)量影響的通徑分析[J]. 遼東學(xué)院學(xué)報(bào)：自然科學(xué)版，2012，19（4）：237-241.

[14]伍漢霖，鐘俊生. 中國(guó)動(dòng)物志：硬骨魚(yú)綱鱸形目（五）蝦虎魚(yú)亞目[M]. 北京：中國(guó)科學(xué)出版社，2008：150-152.

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

[16]何小燕，劉小林，白俊杰，等. 大口黑鱸形態(tài)性狀對(duì)體重的影響效果分析[J]. 水產(chǎn)學(xué)報(bào)，2009，33（4）：597-603.

[17]何 銅，劉小林，楊長(zhǎng)明，等. 凡納濱對(duì)蝦各月齡性狀的主成分與判別分析[J]. 生態(tài)學(xué)報(bào)，2009，29（4）：2134-2142.

[18]唐瞻楊，肖 俊，李莉萍，等. 尼羅羅非魚(yú)（Oreochromis niloticus）不同月齡性狀的主成分與判別分析[J]. 海洋與湖沼，2012，43（2）：288-293.

[19]安 麗，劉 萍，李 健，等. “黃海1號(hào)”中國(guó)明對(duì)蝦形態(tài)性狀對(duì)體質(zhì)量的影響效果分析[J]. 中國(guó)水產(chǎn)科學(xué)，2008，15（5）：779-786.

[20]Turker H，Eversole A G. Evaluation of nondestructive method for determining body composition of crayfish[J]. Journal of Shellfish Research，1998，17（1）：339.

[21]董世瑞，孔 杰，萬(wàn)初坤，等. 中國(guó)對(duì)蝦形態(tài)性狀對(duì)體重影響的通徑分析[J]. 海洋水產(chǎn)研究，2007，28（3）：15-22.