山東省淡水漁業(yè)研究院山東省淡水水產(chǎn)遺傳育種重點(diǎn)實(shí)驗(yàn)室,山東 濟(jì)南 250117
大鱗鲃(Barbuscapito)屬鯉形目鯉科鲃亞科鲃屬,主要分布于黑海、咸海和里海流域,由于其具有食性廣、生長(zhǎng)快、抗逆強(qiáng)、耐鹽堿、肉鮮味美等優(yōu)點(diǎn),深受當(dāng)?shù)厝说臍g迎和喜愛。自2003年黑龍江水產(chǎn)研究所首次將大鱗鲃從國(guó)外引入我國(guó)以來,國(guó)內(nèi)科研工作者對(duì)大鱗鲃的研究日益增多。徐偉等[1]、丁辰龍等[2]、蔡靈等[3]、張棋等[4]和于振海等[5]對(duì)大鱗鲃的親魚培育、人工繁殖、胚胎發(fā)育等方面進(jìn)行了研究和報(bào)道;徐敏等[6]、楊建[7]、耿龍武等[8]對(duì)大鱗鲃耐鹽堿性方面進(jìn)行了研究和報(bào)道;藺玉華等[9]、許紅等[10]和王常安等[11]對(duì)大鱗鲃肌肉營(yíng)養(yǎng)成分及幼魚氨基酸、蛋白質(zhì)需要量等方面進(jìn)行了研究和報(bào)道;藺玉華等[12]、葛彥龍等[13]和耿龍武等[14]對(duì)大鱗鲃種質(zhì)特性進(jìn)行了研究和報(bào)道。為了進(jìn)一步提高育苗及人工養(yǎng)殖的效率,尋找一種有效的評(píng)價(jià)、選擇親魚及各個(gè)養(yǎng)殖階段魚體的方法顯得十分必要。體重指標(biāo)是苗種繁育和遺傳育種的最直接的目標(biāo)性狀,而形態(tài)性狀是與體重密切相關(guān)的重要指標(biāo),因此可以通過尋找與體重相關(guān)的形態(tài)性狀間接地實(shí)現(xiàn)對(duì)大鱗鲃的選擇。
現(xiàn)階段已對(duì)很多養(yǎng)殖魚類的形態(tài)性狀與體重關(guān)系進(jìn)行了研究探討[15~17],確定重要的形態(tài)指標(biāo),為其生長(zhǎng)評(píng)估、良種選育提供重要依據(jù),但對(duì)大鱗鲃形態(tài)性狀與體重關(guān)系的研究還未見報(bào)道。為此,筆者通過對(duì)大鱗鲃的形態(tài)性狀和體重關(guān)系的相關(guān)分析、通徑分析和逐步線性回歸分析,確定與體重密切相關(guān)的主要形態(tài)性狀,建立最優(yōu)回歸方程,以期為大鱗鲃的評(píng)價(jià)和選擇提供依據(jù)。
試驗(yàn)用魚選用山東省淡水漁業(yè)研究院養(yǎng)殖基地2.5齡人工養(yǎng)殖大鱗鲃,隨機(jī)取出30尾健康無損傷的個(gè)體,使用60μg/L的丁香油將其麻醉,進(jìn)行稱重和測(cè)量。
使用直尺測(cè)量全長(zhǎng)(X1)和體長(zhǎng)(X2);使用游標(biāo)卡尺測(cè)量體高(X3)、體寬(X4)、眼間距(X5)、頭長(zhǎng)(X6)、頭高(X7)、頭寬(X8)、吻長(zhǎng)(X9)、尾柄長(zhǎng)(X10)和尾柄高(X11)等11個(gè)形態(tài)指標(biāo),其中X1和X2精確到1cm,X3~X11精確到0.01cm。用電子天平稱體重(Y),精確到0.1g。
使用SPSS 16.0對(duì)形態(tài)性狀和體重進(jìn)行統(tǒng)計(jì)分析,通過相關(guān)分析和通徑分析確定各形態(tài)性狀對(duì)體重的直接與間接作用,再通過逐步回歸分析,構(gòu)建形態(tài)性狀對(duì)體重的多元回歸方程。
所測(cè)大鱗鲃各形態(tài)性狀數(shù)據(jù)結(jié)果見表1。由表1結(jié)果來看,體重的變異系數(shù)最大,達(dá)36.57%,表明大鱗鲃體重的變化范圍最大,具有較大的選擇空間,實(shí)際選育過程中應(yīng)作為最主要的目標(biāo)性狀。
表1 各性狀表型統(tǒng)計(jì)結(jié)果
大鱗鲃各形態(tài)性狀間的相關(guān)系數(shù)見表2。由檢驗(yàn)結(jié)果可知,各形態(tài)性狀之間的相關(guān)程度均達(dá)到了極顯著水平(P<0.01),各性狀與體重之間的相關(guān)程度也均達(dá)到了極顯著水平(P<0.01)。其中,體高與體重的相關(guān)程度最大,頭高與體重的相關(guān)程度最小,但均達(dá)到極顯著水平,說明選擇與體重相關(guān)的指標(biāo)進(jìn)行選擇具有可行性。
運(yùn)用逐步回歸分析的方法建立多元回歸方程,最終回歸方程中只包括了2個(gè)形態(tài)性狀X1和X3,結(jié)果見表3。從回歸系數(shù)看,這2個(gè)自變量與體重的回歸關(guān)系都達(dá)到極顯著水平(P<0.01),說明這2個(gè)性狀對(duì)體重的影響是顯著的,而其他形態(tài)性狀對(duì)體重的影響不顯著,所以被全部剔除。大鱗鲃形態(tài)性狀對(duì)體重的多元回歸方程式為Y=-204.669+0.641X1+4.187X3。
回歸預(yù)測(cè)結(jié)果表明,回歸關(guān)系達(dá)到極顯著水平(P<0.01),其校正復(fù)相關(guān)指數(shù)為0.956,因此該方程可用于大鱗鲃的評(píng)價(jià)和選擇。
表3 偏回歸系數(shù)和回歸常數(shù)檢驗(yàn)
通徑系數(shù)定義為變量標(biāo)準(zhǔn)化的偏回歸系數(shù),用來表示相關(guān)變量因果關(guān)系的統(tǒng)計(jì)量。經(jīng)顯著性差異檢驗(yàn),剔除通徑系數(shù)不顯著的變量,保留X3和X1這2個(gè)形態(tài)性狀變量,這2個(gè)形態(tài)性狀對(duì)體重的通徑系數(shù)分析結(jié)果見表4。保留下的2個(gè)變量X3和X1對(duì)Y的直接作用P值分別為0.563和0.439。相關(guān)系數(shù)由各性狀的直接作用和通過其他性狀的間接作用組成,體高對(duì)體重的直接作用大于全長(zhǎng)對(duì)體重的直接作用,則體高通過其他性狀對(duì)體重的間接作用小于全長(zhǎng)通過其他性狀對(duì)體重的間接作用。
表4 保留性狀對(duì)體重的通徑系數(shù)
表5 保留性狀對(duì)體重的決定系數(shù)
2個(gè)保留性狀對(duì)體重的決定系數(shù)見表5。結(jié)果顯示,體高對(duì)體重的決定系數(shù)(0.317)大于全長(zhǎng)對(duì)體重的決定系數(shù)(0.193),直接決定系數(shù)和間接決定系數(shù)總和為0.961(0.317+0.451+0.193=0.961),表明選取的這2個(gè)性狀是影響2.5齡大鱗鲃體重的重要性狀。
一般來說,商品個(gè)體體重越大,經(jīng)濟(jì)價(jià)值越高,體重是決定商品品質(zhì)和市場(chǎng)價(jià)格的重要性狀,因此在選擇育種過程中通常被認(rèn)為是最重要的目標(biāo)性狀。但實(shí)際測(cè)量中體重的變異系數(shù)最大,準(zhǔn)確性不高,而體長(zhǎng)、體高、體寬等形態(tài)性狀變異系數(shù)相對(duì)較小,且均與體重存在顯著或極顯著相關(guān)關(guān)系,因此為了確保準(zhǔn)確性,減少誤差,可通過其他性狀間接選擇以獲得理想的選擇效果。本研究相關(guān)分析結(jié)果顯示各性狀與體重之間的相關(guān)程度均達(dá)到了極顯著水平(P<0.01),但是它反映的是變量間的綜合關(guān)系,不能反映變量間的真實(shí)關(guān)系以及自變量對(duì)因變量產(chǎn)生影響的大小。因此在相關(guān)分析的基礎(chǔ)上進(jìn)行了通徑分析,結(jié)果表明,體高對(duì)2.5齡大鱗鲃體重的直接影響最大,全長(zhǎng)通過體高間接影響體重。保留的2個(gè)性狀對(duì)2.5齡大鱗鲃體重的直接和間接決定系數(shù)之和為0.961(≥0.85),說明體高和全長(zhǎng)是影響2.5齡大鱗鲃體重的主要形態(tài)性狀。
隨著魚類育種研究的進(jìn)行,形態(tài)性狀多元回歸分析已廣泛應(yīng)用其中。李炎璐等[18]通過多元回歸分析建立了以體高、體長(zhǎng)、尾柄高為自變量,以體重為因變量的黃鰭金槍魚的最優(yōu)多元回歸方程。趙旺等[19]通過多元回歸分析建立了4月齡斜帶石斑魚形態(tài)性狀與體質(zhì)量的回歸方程。劉峰等[20]應(yīng)用逐步多元回歸分析建立了以體質(zhì)量為因變量,以體長(zhǎng)、軀干長(zhǎng)、尾柄高和體高為自變量的小黃魚的回歸方程。研究通過逐步引入-剔除的方法對(duì)2.5齡大鱗鲃形態(tài)性狀和體重進(jìn)行多元回歸分析,以體高和全長(zhǎng)為自變量,以體重為因變量構(gòu)建了多元回歸方程Y=-204.669+0.641X1+4.187X3,構(gòu)建的多元回歸方程經(jīng)檢驗(yàn)達(dá)到極顯著水平,表明回歸方程有意義,可用于對(duì)2.5齡大鱗鲃的預(yù)測(cè)和估計(jì)。