方 偉 ，周勝杰 ，趙 旺，楊 蕊，胡 靜，于 剛，馬振華

（1.中國水產(chǎn)科學(xué)研究院南海水產(chǎn)研究所熱帶水產(chǎn)研究開發(fā)中心，海南 三亞 572018；2.農(nóng)業(yè)農(nóng)村部南海漁業(yè)資源開發(fā)利用重點試驗室，廣東 廣州 510300；3.三亞熱帶水產(chǎn)研究院，海南 三亞 572018）

黃鰭金槍魚 (Thunnus albacores) 隸屬鱸形目、鯖亞目、鯖科、金槍魚屬，俗稱鮪魚、吞拿魚[1]，是一種生活在海洋中上層、距離海洋表層水溫8 ℃以內(nèi)水域的高度跨洋性洄游魚類[2-4]，主要分布于太平洋、印度洋、大西洋的熱帶、亞熱帶及溫帶廣闊海域，以赤道附近最多，國內(nèi)的黃鰭金槍魚主要分布在南海和臺灣附近海域[5-9]。其肉質(zhì)鮮嫩、營養(yǎng)豐富、綠色無污染[10-11]，已成為當(dāng)今世界海洋漁業(yè)發(fā)展中最具開發(fā)價值的魚種之一。目前針對黃鰭金槍魚的研究已涉及種群多樣性[12-14]、生理生化[15-17]、營養(yǎng)[18-20]、基因組學(xué)[21-23]等方面，但其人工繁育卻鮮有報道。缺乏系統(tǒng)的優(yōu)良性狀選育研究導(dǎo)致其人工繁育效果不理想且難以推廣應(yīng)用。體質(zhì)量是動物遺傳親本選育的常見指標(biāo)，結(jié)合形態(tài)性狀對體質(zhì)量進行通徑分析可有效提高育種效率[24-27]。目前，國內(nèi)通過研究形態(tài)性狀與體質(zhì)量的相互關(guān)系來實現(xiàn)優(yōu)勢性狀的篩選已廣泛應(yīng)用于動植物的苗種繁育，尤其是水產(chǎn)類的親本選育。魏海軍等[28]在仙女蛤 (Callista erycina) 選育過程中發(fā)現(xiàn)，各表型形態(tài)性狀間均呈極顯著相關(guān)，通徑分析結(jié)果確定殼寬為主要選擇性狀；劉文廣等[29]發(fā)現(xiàn)華貴櫛孔扇貝 (Mimachlamys nobilis) 不同貝齡的形態(tài)性狀對體質(zhì)量相關(guān)系數(shù)最高的是尾殼長，而通徑系數(shù)最小的是尾殼寬；吳水清等[30]發(fā)現(xiàn)雜交石斑魚 [斜帶石斑魚 (Epinephelus coioides) ×赤點石斑魚 (E.akaara)]與親本的生長存在差異，并確定雜交子代形態(tài)性狀更偏向于母本。在魚類選育的研究中往往以成魚為對象進行多元分析，但不同生長階段的選育指標(biāo)可能存在明顯差異。因此，研究幼魚形態(tài)性狀對體質(zhì)量的影響有利于盡早開展幼魚階段的選育工作，減少育種工作量并提高育種效率。本研究對5 月齡黃鰭金槍魚進行隨機采樣，獲取其體質(zhì)量和形態(tài)性狀的相關(guān)數(shù)據(jù)，采用多元回歸分析確定了影響體質(zhì)量的主要性狀，并建立形態(tài)性狀對體質(zhì)量的回歸方程，揭示其形態(tài)性狀與體質(zhì)量的相互關(guān)系，以期為黃鰭金槍魚優(yōu)良形態(tài)性狀的選育提供理論參考。

1 材料與方法

1.1 實驗材料

本實驗隨機選取生長狀態(tài)良好的65 尾黃鰭金槍魚 (由中國水產(chǎn)科學(xué)研究院南海水產(chǎn)研究所熱帶水產(chǎn)研究開發(fā)中心提供)，循環(huán)水養(yǎng)殖于室內(nèi)水泥池 (直徑5 m、高3 m，圓筒形)，每天換水量100%，餌料主要以冰鮮雜魚為主，投喂量為魚體質(zhì)量的2%～5%，養(yǎng)殖密度約5 kg·m?3，水溫28～32 ℃，鹽度27～34，pH 7.5～8.5，溶解氧質(zhì)量濃度≥7.0 mg·L?1，氨氮質(zhì)量濃度≤0.2 mg·L?1，亞硝酸鹽質(zhì)量濃度≤0.02 mg·L?1。

1.2 測量方法

測定體質(zhì)量 (y) 和14 個形態(tài)性狀 (圖1)：體長(x1)、全長 (x2)、體高 (x3)、頭長 (x4)、眼徑 (x5)、上顎長 (x6)、下顎長 (x7)、胸鰭長 (x8)、臀鰭長(x9)、尾鰭長 (x10)、尾鰭寬 (x11)、第Ⅰ背鰭長(x12)、第Ⅱ背鰭長 (x13)、尾柄高 (x14)。體質(zhì)量稱量使用電子天平 (精確度為0.01 g)，測量長度指標(biāo)使用游標(biāo)卡尺 (精確度為0.02 mm)。

圖1 黃鰭金槍魚各形態(tài)特征指標(biāo)測量示意圖x1.體長；x2.全長；x3.體高；x4.頭長；x5.眼徑；x6.上顎長；x7.下顎長；x8.胸鰭長；x9.臀鰭長；x10.尾鰭長；x11.尾鰭寬；x12.第Ⅰ背鰭長；x13.第Ⅱ背鰭長；x14.尾柄高；下同F(xiàn)igure 1 Measurement for each morphological index of T.albacoresx1.Body length;x2.Total length;x3.Body height;x4.Head length;x5.Eye diameter;x6.Palate length;x7.Jaw length;x8.Pectoral fin length;x9.Anal fin length;x10.Caudal fin length;x11.Caudal fin width;x12.1st dorsal fin length;x13.2nd dorsal fin length;x14.Caudal height;the same below

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

采用軟件SPSS 19.0 對體質(zhì)量和各形態(tài)性狀等相關(guān)數(shù)據(jù)進行統(tǒng)計分析。首先獲得各性狀的描述性統(tǒng)計結(jié)果，然后通過相關(guān)分析和通徑分析確定各形態(tài)性狀對體質(zhì)量的直接與間接作用大小，最后經(jīng)回歸分析構(gòu)建黃鰭金槍魚形態(tài)性狀對體質(zhì)量的多元回歸方程及偏回歸系數(shù)。

2 結(jié)果

2.1 黃鰭金槍魚各性狀的描述性結(jié)果

黃鰭金槍魚y的變異系數(shù)最大 (38.96%)，其他性狀的變異系數(shù)介于10.04%～31.45% (表1)，說明y在黃鰭金槍魚的各性狀中具有較大的選擇潛力，可作為優(yōu)質(zhì)親本選育的目標(biāo)性狀。其他形態(tài)性狀變異系數(shù)從大到小依次為x13、x7、x12、x4、x6、x9、x8、x5、x14、x11、x10、x3、x2和x1。

表1 黃鰭金槍魚體質(zhì)量與形態(tài)性狀的描述性結(jié)果Table 1 Descriptive results of body mass and morphological traits of T.albacores

2.2 黃鰭金槍魚各性狀間的相關(guān)分析

黃鰭金槍魚各形態(tài)性狀間的相關(guān)性均達極顯著水平 (P<0.01，表2)，其中x4與x6的相關(guān)系數(shù)最大 (0.989)，x7與x10的最小 (0.653)。此外，各形態(tài)性狀與體質(zhì)量的相關(guān)系數(shù)各不相同，其中最大的為x1(0.981)，其次為x14(0.961)，最小為x7(0.827)。

2.3 黃鰭金槍魚形態(tài)性狀對體質(zhì)量的通徑分析

黃鰭金槍魚的通徑分析結(jié)果見表3，其中體質(zhì)量為因變量，形態(tài)性狀為自變量。經(jīng)通徑系數(shù)的顯著性差異檢驗，剔除通徑系數(shù)不顯著的變量，保留4 個形態(tài)變量 (x2、x4、x6、x8)。根據(jù)相關(guān)系數(shù)的組成效應(yīng)，將黃鰭金槍魚形態(tài)性狀與體質(zhì)量的相關(guān)系數(shù)分解為各性狀的直接作用和間接作用兩部分。從直接作用來看，影響體質(zhì)量的直接作用系數(shù)最大的為x8(0.507,P<0.01)，其次依次為x2(0.307)、x4(0.181)，其中x6(?0.015) 對體質(zhì)量的直接作用系數(shù)為負值，說明其對體質(zhì)量的影響為負向作用，但x6通過與其他性狀產(chǎn)生較大的間接作用，抵消了負向作用；從間接作用來看，除x8外，黃鰭金槍魚各形態(tài)性狀對體質(zhì)量的直接作用均小于間接作用。其他性狀對體質(zhì)量的間接作用大小為0.650～0.968，主要通過x8間接影響體質(zhì)量。

2.4 黃鰭金槍魚形態(tài)性狀對體質(zhì)量的決定程度

黃鰭金槍魚主要形態(tài)性狀對體質(zhì)量的直接決定系數(shù)與間接決定系數(shù)的總和為0.942 (表4)，表明選取的x2、x4、x6和x8是影響黃鰭金槍魚體質(zhì)量的重要性狀，其他性狀對體質(zhì)量的影響較小。此外，以上4 個性狀對體質(zhì)量的直接決定程度存在差異，其中x8的直接決定系數(shù)最大 (0.257)，x6最小(0.001)。x2和x8對體質(zhì)量決定系數(shù)的共同作用最大 (0.299)。

表2 黃鰭金槍魚各形態(tài)性狀間的相關(guān)系數(shù)Table 2 Correlation coefficients among morphological traits of T.albacores

表3 黃鰭金槍魚各形態(tài)性狀對體質(zhì)量的通徑分析Table 3 Path analysis of morphological traits to body mass of T.albacores

表4 黃鰭金槍魚各形態(tài)性狀對體質(zhì)量的決定系數(shù)Table 4 Coefficient of determination of morphological traits to body mass of T.albacores

2.5 多元回歸方程的構(gòu)建

通過多元回歸分析剔除偏回歸系數(shù)不顯著的形態(tài)性狀，利用偏回歸系數(shù)顯著的形態(tài)性狀與體質(zhì)量建立黃鰭金槍魚形態(tài)性狀與體質(zhì)量的多元回歸方程[31]：y=?462.621+1.157x2+1.253x4+0.424x6+4.707x8。方差分析結(jié)果顯示，多元回歸方程的回歸關(guān)系達極顯著水平 (P<0.01,R2=0.921)。經(jīng)顯著性檢驗該回歸方程的偏回歸系數(shù)，所選的x2、x4、x6和x8對體質(zhì)量的偏回歸系數(shù)達顯著或極顯著水平 (x2:t=3.578,P<0.05;x4:t=1.198,P<0.01;x6:t=8.146,P<0.01;x8:t=1.974,P<0.01)。

2.6 曲線模型擬合結(jié)果

分別以x2、x4、x6和x8為自變量、y為因變量進行曲線模型擬合 (表5)。4 個形態(tài)性狀與體質(zhì)量的最優(yōu)擬合模型分別為二次函數(shù)、指數(shù)函數(shù)、冪函數(shù)和二次函數(shù)，其公式分別為y=?6.598x2+0.016x+795.394、y=40.412e0.026x、y=1.229x1.706、y=?4.937x2+0.104x+174.317，對應(yīng)的R2分別為0.972、0.952、0.937 和0.956。

表5 黃鰭金槍魚形態(tài)性狀與體質(zhì)量的曲線模型擬合Table 5 Curve model fitting of morphological traits and body mass of T.albacores

3 討論

目前，在魚類品種選育和經(jīng)濟性狀選擇中，體質(zhì)量作為重要的評價指標(biāo)已被廣泛應(yīng)用。在已報道的有關(guān)魚類體質(zhì)量的研究中發(fā)現(xiàn)，體質(zhì)量一般受形態(tài)性狀影響，并通過形態(tài)性狀反映出來[32-34]。黃小林等[35]采用通徑分析研究了影響青龍斑 (Epinephelusspp.) 和淡水鯊魚 (Pangasianodon hypophthalmus) 體質(zhì)量的主要性狀，揭示其軀干長、頭高、眼間距等是影響魚類體質(zhì)量的重要性狀；李培倫等[36]對大麻哈魚 (Oncorhynchus keta) 的研究表明，叉長、頭長和尾柄長是影響體質(zhì)量的關(guān)鍵性狀；馮冰冰等[37]揭示了體長是影響大鱗鲃 (Luciobarbus capito) 體質(zhì)量的主要形態(tài)性狀。本研究測量了5 月齡黃鰭金槍魚的體質(zhì)量和14 個形態(tài)性狀指標(biāo)，相關(guān)分析顯示各形態(tài)性狀與體質(zhì)量極顯著相關(guān)，但不能確定影響體質(zhì)量的主要形態(tài)性狀。經(jīng)通徑分析量化形態(tài)性狀與體質(zhì)量間的相互關(guān)系后發(fā)現(xiàn)，胸鰭長對體質(zhì)量的直接作用最大，其次為全長、頭長、上顎長，而上顎長對體質(zhì)量的直接影響為負向作用，其影響主要通過全長、頭長、胸鰭長來呈現(xiàn)。這4 個性狀對體質(zhì)量的共同決定系數(shù)之和為0.942，大于0.85，表明這些性狀是影響體質(zhì)量的主要性狀[38]。前人研究發(fā)現(xiàn)體長、全長和頭長等是最常見的關(guān)鍵影響性狀，而胸鰭長是黃鰭金槍魚的主要形態(tài)性狀，說明不同品種魚類的形態(tài)性狀對體質(zhì)量的影響存在差異。趙旺等[39]在斜帶石斑魚的通徑分析中指出不同品種的魚類，因為地理群體、生長需求、形態(tài)結(jié)構(gòu)等差異，其對體質(zhì)量產(chǎn)生主要影響的形態(tài)性狀也存在差異。對于遠洋魚類，尤其高度跨洋性洄游的黃鰭金槍魚，需要充足的動力來擴大生存空間，胸鰭展現(xiàn)出優(yōu)秀的力學(xué)特征滿足了這一需求。張照煌和李魏魏[40]在座頭鯨(Megaptera novaeangliae) 胸鰭動力學(xué)特性研究中發(fā)現(xiàn)，胸鰭仿生在流動控制及葉片增效方面有著極大的研究價值，進一步證實了胸鰭的重要性，且偏回歸系數(shù)的顯著性檢驗結(jié)果顯示，全長、頭長、上顎長和胸鰭長與體質(zhì)量的偏回歸系數(shù)達顯著或極顯著水平。郭華陽等[41]研究發(fā)現(xiàn)擬合曲線模型在一定程度上能更準確地反映單一形態(tài)性狀與體質(zhì)量的相互關(guān)系。本研究的曲線模擬結(jié)果顯示全長、頭長、上顎長和胸鰭長這4 個形態(tài)性狀與黃鰭金槍魚體質(zhì)量的6 種曲線模型擬合結(jié)果均達極顯著水平 (P<0.01)，說明擬合所獲得的曲線模型均有意義，從R2和F值來看，不同性狀的曲線模型存在差異，以上4 個性狀與體質(zhì)量的最優(yōu)擬合模型分別為二次函數(shù)、指數(shù)函數(shù)、冪函數(shù)和二次函數(shù)。因此，本研究認為在品種選育時應(yīng)以胸鰭長為最主要的選擇性狀，全長、頭長及上顎長為輔助選擇性狀。這與實際生產(chǎn)相符，黃鰭金槍魚為肉食性遠洋魚類，位于魚體前中部的胸鰭越長，不僅可提高身體的靈活性，還可擴大活動區(qū)域，從而捕獲更多的食物來滿足機體對能量的需求。

據(jù)一些研究報道，影響魚類體質(zhì)量的形態(tài)性狀在不同生長時期或不同養(yǎng)殖環(huán)境可能存在差異。佟雪紅等[42]發(fā)現(xiàn)牙鲆 (Paralichthys olivaceus) 20 日齡時影響其體質(zhì)量的主要性狀是體高，而在其他生長階段則是體長；劉賢德等[43]在大黃魚 (Larimichthys crocea) 的選育工作中發(fā)現(xiàn)，早期 (13 月齡) 應(yīng)以體高和體長為主，而中期階段則要考慮體高、體長和全長。此外，不同生長 (或養(yǎng)殖) 環(huán)境對形態(tài)性狀之間的相關(guān)性也存在差異。李俊偉等[44]在四指馬鲅 (Eleutheronema tetradactylum) 的研究中發(fā)現(xiàn)，室內(nèi)循環(huán)水養(yǎng)殖群體與體質(zhì)量相關(guān)程度較大的形態(tài)性狀為叉長、體長及頭長，而在池塘養(yǎng)殖群體中則為全長、尾柄高。本研究確定了全長、頭長、上顎長和胸鰭長這4 個密切影響黃鰭金槍魚幼魚體質(zhì)量的形態(tài)性狀，但在其他生長階段中形態(tài)性狀與體質(zhì)量的關(guān)系如何、是否仍是影響體質(zhì)量的主要性狀，還有待進一步研究。

4 結(jié)論

黃鰭金槍魚體質(zhì)量的變異系數(shù)是所有性狀中最大的，這說明在親本選育過程中，體質(zhì)量作為常見的參考數(shù)據(jù)具有一定的指導(dǎo)意義，但僅以體質(zhì)量作為參照指標(biāo)進行親本選育，往往會受到環(huán)境及其他因素的影響而產(chǎn)生較大的系統(tǒng)誤差。本研究認為在對黃鰭金槍魚5 月齡幼魚進行以體質(zhì)量為主要目標(biāo)的親本選育時，應(yīng)以胸鰭長為主要選擇性狀，全長、頭長、上顎長為輔助選擇性狀作為參照，這不僅能夠提高優(yōu)良親本選育的精準度，也為進一步的實際生產(chǎn)作業(yè)提供理論依據(jù)。