李 莉 , 鄭永允 , 徐科鳳, 邱兆星 , 郭萍萍

(1.山東省海洋生物研究院, 山東 青島 266002; 2.青島市淺海底棲漁業(yè)增殖重點實驗室, 山東 青島266002; 3.青島國家海洋科學研究中心, 山東 青島 266071)

毛蚶(Scapharca subcrenataLischke)隸屬于軟體動物門(Mollusca)、瓣鰓綱(Lamellibranchia)、蚶目(Arcoida)、蚶科(Arcidae)、毛蚶屬(Scapharca), 廣泛分布于日本、朝鮮、中國沿海, 在我國北起鴨綠江、南至廣西均有分布, 以萊州灣、渤海灣、遼東灣、海州灣等淺海區(qū)資源最為豐富[1]。20世紀90年代以來,由于過度捕撈、生態(tài)環(huán)境破壞等因素, 毛蚶資源逐年衰退, 已遠遠不能滿足市場需求, 價格不斷攀升。因此, 近幾年毛蚶增養(yǎng)殖工作越來越受到人們的重視,并已逐漸成為我國南北方沿海重要的增養(yǎng)殖經(jīng)濟貝類之一。

生物學特征研究對于開展一個物種的增養(yǎng)殖是十分必要的, 尤其是有關(guān)水產(chǎn)經(jīng)濟動物經(jīng)濟性狀相關(guān)性研究的報道較多, 所涉及的物種有魚類[2-3]、貝類[4-5]以及甲殼類[6-8]等數(shù)十種。如劉賢德等[9]對不同生長時期大黃魚的研究表明, 早期選種應(yīng)注重對體高和體長的挑選, 而后期選種時全長也是保證選擇效果的考慮因素; 宋堅等[10]對長蛸的研究建立了體質(zhì)量預(yù)測凈肉質(zhì)量的最優(yōu)模型。這些研究通過分析形態(tài)學性狀與體質(zhì)量的關(guān)系, 為水產(chǎn)經(jīng)濟動物的遺傳育種、增養(yǎng)殖方案優(yōu)化和生產(chǎn)量估計等提供理論依據(jù)。

關(guān)于毛蚶各形態(tài)性狀與質(zhì)量性狀間的相關(guān)分析已有研究報道[11-13], 但上述研究對象均為毛蚶成貝,不同貝齡時期毛蚶各性狀的相互關(guān)系, 尤其是殼性狀對體質(zhì)量的影響效果尚未見報道。本研究分別分析了 3個不同貝齡期毛蚶殼形態(tài)性狀與體質(zhì)量之間的相互關(guān)系, 以確定不同貝齡毛蚶不同殼形態(tài)性狀對體質(zhì)量的影響效果, 完善毛蚶的生物學資料, 同時為毛蚶的遺傳育種及其資源的合理利用等提供基礎(chǔ)資料。

1 材料與方法

所用毛蚶來源于2011年10月底播增殖個體, 底播于山東省廣饒縣海區(qū)。底播苗種為當年人工繁育苗種(1齡), 2012年10月、2013年10月回捕個體分別為 2齡、3齡。隨機抽取 3個貝齡期毛蚶個體各100粒, 測量其殼長(SL)、殼高(SH)、殼寬(SW)和體質(zhì)量(BW)。殼形態(tài)指標用數(shù)顯測量尺測量, 精確到0.01 mm。活體質(zhì)量用電子天平稱量, 精確到0.01 g。

使用SPSS16.0統(tǒng)計軟件進行表型性狀間的相關(guān)分析和通徑分析。通徑分析參照杜家菊和陳志偉[14]的方法進行, 設(shè)BW為因變量, SL、SH、SW為自變量, 因變量實施正態(tài)性檢驗后, 采用逐步多元線性回歸法建立殼性狀對質(zhì)量性狀的線性回歸方程。

依據(jù)相關(guān)分析結(jié)果和通經(jīng)系數(shù)計算決定系數(shù),參照劉春雷等[15]的計算方法。

2 結(jié)果

2.1 各性狀的表型參數(shù)統(tǒng)計量

各性狀的表型參數(shù)統(tǒng)計信息見表1, 不同貝齡期毛蚶殼形態(tài)性狀和體質(zhì)量等表型性狀的變異趨勢是一致的, 均表現(xiàn)為體質(zhì)量>殼寬>殼高>殼長, 其中1齡貝的變異系數(shù)最大, 2齡貝的變異系數(shù)最小。

表1 不同貝齡期各性狀表型的參數(shù)統(tǒng)計Tab.1 The descriptive of morphometric traits of S.subcrenata at different years of age

2.2 各性狀間的相關(guān)系數(shù)

各性狀的相關(guān)分析結(jié)果見表2, 不同貝齡毛蚶各性狀間呈正相關(guān), 且均達到極顯著水平(P<0.001),表明所選指標進行相關(guān)分析具有重要的現(xiàn)實意義。從毛蚶殼形態(tài)性狀與體質(zhì)量的相關(guān)性來看, 1齡殼形態(tài)性狀中殼長與體質(zhì)量之間的相關(guān)系數(shù)最大, 殼寬與體質(zhì)量之間的相關(guān)系數(shù)最小; 2齡與3齡均表現(xiàn)為殼寬與體質(zhì)量之間的相關(guān)系數(shù)最大。

表2 不同貝齡期各性狀間的相關(guān)系數(shù)Tab.2 Correlation coefficients among quantitative of S.subcrenata at different years of age

2.3 殼形態(tài)性狀對體質(zhì)量的通徑分析

以體質(zhì)量為因變量, 利用單個樣本Kolmogorov-Smirnov Test進行正態(tài)性檢驗。1齡、2齡、3齡體質(zhì)量K-S分析的Z值依次為0.847、0.793和0.706, 顯著性水平P值(0.469、0.555、0.701)大于0.05, 表明因變量服從正太分布, 可以進行通徑分析。

不同殼形態(tài)性狀對體質(zhì)量的影響作用見表3。不同貝齡期毛蚶殼性狀對體質(zhì)量的影響作用不同。1齡時, 殼長是影響體質(zhì)量的主要因素, 殼高、殼寬的間接影響作用大, 且均是通過殼長間接影響體質(zhì)量; 2齡時, 殼長、殼寬是影響體質(zhì)量的主要因素, 殼高的間接作用最大; 3齡時, 殼寬是影響體質(zhì)量的主要因素, 殼長的間接作用最大。

表3 不同貝齡期各形態(tài)性狀對體質(zhì)量的通徑分析Tab.3 Path analysis of morphological traits to total weight of S.subcrenata at different years of age

2.4 殼形態(tài)性狀對體質(zhì)量的決定程度分析

不同殼形態(tài)性狀對體質(zhì)量的決定程度分析見表4。不同貝齡期毛蚶殼形態(tài)性狀對體質(zhì)量的決定作用不同。1齡時, 殼長是決定體質(zhì)量的主要因素, 其直接決定系數(shù)最大(0.202), 與殼高對體質(zhì)量的共同決定系數(shù)最大(0.219); 2齡時, 殼長、殼寬是決定體質(zhì)量的主要因素, 直接決定系數(shù)依次為0.143和0.140,且兩者的共同決定系數(shù)最大(0.244); 3齡時, 殼寬是決定體質(zhì)量的主要因素(0.158), 其與殼高的共同決定系數(shù)最大(0.230)。不同貝齡期毛蚶殼形態(tài)性狀對體質(zhì)量的決定系數(shù)總和依次為0.890、0.891和0.921,與相關(guān)性指數(shù)基本一致, 均大于0.85, 表明所選3個殼形態(tài)性狀是影響不同貝齡期毛蚶體質(zhì)量的主要性狀, 其他性狀的影響相對較小。

表4 不同貝齡期各形態(tài)性狀對體質(zhì)量的決定系數(shù)Tab.4 Coefficient of determination of morphological traits to total weight of S.subcrenata at different years of age

2.5 多元回歸方程的建立

以體質(zhì)量為依變量, 殼形態(tài)性狀為自變量, 采用逐步回歸的方法, 進行多元線性回歸分析。多元共線性診斷以方差膨脹因子(VIF)為指標。通常情況下,當VIF≥10時, 認為自變量間存在嚴重共線性。從表5看, 剔除存在嚴重共線性關(guān)系的自變量, 模型1中所列性狀的方差膨脹因子均小于10。

表5 方差膨脹因子Tab.5 Variance inflation factors

建立 3個不同貝齡期毛蚶的線性回歸方程, 其線性回歸系數(shù)都達到極顯著水平(P< 0.001), 表明各自變量對于因變量存在顯著的線性關(guān)系, 擬合的線性統(tǒng)計模型及相應(yīng)參數(shù)的估計可靠。

式中, BW1、BW2、BW3依次為 1齡、2齡、3齡毛蚶的體質(zhì)量, SL為殼長, SH為殼高, SW為殼寬。

3 討論

體質(zhì)量的增長是水生經(jīng)濟動物養(yǎng)殖的關(guān)鍵, 它是養(yǎng)殖產(chǎn)量與周期的體現(xiàn), 而體質(zhì)量生長性能優(yōu)劣往往通過體型表現(xiàn)出來, 這是遺傳育種的理論依據(jù),也是選苗的依據(jù)。利用多元回歸分析研究形態(tài)性狀與體質(zhì)量之間的關(guān)系已廣泛應(yīng)用于水生動物研究中。Ahmed和 Abbas[16]對魚和貝類幼齡期體長、體質(zhì)量相關(guān)的生長參數(shù)進行了多元相關(guān)分析; 劉春雷等[15]運用相關(guān)分析研究了轉(zhuǎn)大麻哈魚生長激素基因鯉與對照鯉10個表型性狀對體質(zhì)量的影響程度; 孫振興等[17]對扁玉螺表型形態(tài)性狀與質(zhì)量性狀等 8個參數(shù)進行了相關(guān)分析。本研究中利用多元回歸分析研究了海區(qū)底播增殖毛蚶不同貝齡殼長、殼高、殼寬等 3個殼形態(tài)性狀與體質(zhì)量的相關(guān)性, 相關(guān)系數(shù)均呈正相關(guān), 并達到極顯著水平(P< 0.001), 研究表明, 進一步的統(tǒng)計分析具有現(xiàn)實意義。

變異系數(shù)是衡量各觀測值變異程度的一個統(tǒng)計量。本研究中2齡貝各性狀的變異系數(shù)小于3齡貝,表現(xiàn)為 3個不同齡中的最小值, 分析其主要原因可能有兩個方面: 一是在底播后苗種受自然環(huán)境變化與敵害生物捕食等的影響, 經(jīng)歷了自然淘汰; 二是在拖捕過程中網(wǎng)具的使用可能將毛蚶進行了分選,使得捕獲毛蚶規(guī)格趨于整齊。

性狀間相關(guān)性是變量間相互關(guān)系的綜合作用,包含了兩者的之間關(guān)系和通過其他變量的間接關(guān)系[18-19]。一般采用相關(guān)系數(shù)、通徑系數(shù)、決定程度及多元回歸等綜合分析, 從而提高結(jié)果的可信度[20]。本研究中, 在表型相關(guān)分析的基礎(chǔ)上進行通徑分析和決定系數(shù)分析, 不同貝齡毛蚶殼形態(tài)性狀對體質(zhì)量的相關(guān)指數(shù)均大于 0.85, 表明所列性狀是影響毛蚶體質(zhì)量的主要性狀[21], 并且影響不同貝齡期毛蚶體質(zhì)量的主要因素是不同的。殼長是影響 1齡毛蚶體質(zhì)量的主要因素, 殼高、殼寬均是通過殼長間接影響體質(zhì)量; 2齡時, 殼長、殼寬是影響體質(zhì)量的主要因素, 且兩者的共同決定系數(shù)最大; 3齡時, 殼寬是影響體質(zhì)量的主要因素, 其與殼高的共同決定系數(shù)最大。宋堅等[22]對 3種規(guī)格硬殼蛤的研究表明小規(guī)格硬殼蛤的殼寬、中規(guī)格和大規(guī)格硬殼蛤的殼長分別是影響體質(zhì)量的的主要因素; 劉文廣等[23]對 5種規(guī)格華貴櫛孔扇貝的研究表明殼高是影響體質(zhì)量的主要指標。說明貝殼的形狀隨著貝類的生長是不斷發(fā)生變化的, 而不同雙殼貝類由于其生物學特征不同, 其貝殼生長變化的趨勢不同, 影響其體質(zhì)量的貝殼性狀也不同。

同時, 本研究結(jié)果與毛蚶相關(guān)研究報道之間存在相同之處, 也存在不同之處。如王輝等[11]對2齡南海毛蚶雌雄兩個性別的分析表明, 殼高是影響雄性毛蚶體質(zhì)量的主要因素, 殼寬是影響雌性的主要因素。陳麗梅等[12]對渤海地區(qū)毛蚶的分析表明, 殼高對體質(zhì)量的直接作用和綜合決定系數(shù)最大, 是影響體質(zhì)量的主要因素。郝偉等[13]對不同性別 3齡池塘底播養(yǎng)殖毛蚶的分析表明, 殼長、殼高是影響雌性組體質(zhì)量的主要因素, 殼長、殼寬是影響雄性組的主要因素, 殼長對混合組的影響最大?？梢? 不同的遺傳背景、年齡、生長環(huán)境與狀況等要素均會影響毛蚶貝殼形態(tài)性狀對體質(zhì)量的影響效果。該現(xiàn)象也常見于其它經(jīng)濟貝類的研究中, 如鄭懷平等[24]對 1齡華貴櫛孔扇貝分析認為, 選育指標應(yīng)以殼高為主, 同時加強對殼長的協(xié)同選擇; 而劉志剛等[25]則認為應(yīng)以殼高作為第一選育目標性狀。

本研究中通過對1齡、2齡以及3齡毛蚶殼形態(tài)性狀與體質(zhì)量的相關(guān)性分析, 進一步摸清了毛蚶殼長、殼高、殼寬等主要殼形態(tài)性狀隨著生長的變化及其對體質(zhì)量的影響, 利用多元線性回歸分析得到了 3個回歸方程, 可分別應(yīng)用于毛蚶不同養(yǎng)殖時期增養(yǎng)殖方案的制定與生產(chǎn)量的估計, 同時為下一步毛蚶選育工作的開展提供了理論依據(jù)。

[1] 陳建華, 閆斌倫, 高煥.毛蚶生物學特性及其研究進展[J].河北漁業(yè), 2006, 9: 24-25, 49.

[2] Deboski P, Dobosz S, Robak S, et al.Fat level in body of juvenile Atlantic salmon (Salmo salarL.), and sea trout (Salmo truttaM.truttaL.), and method of estimation from morphometric data [J].Archives of Polish Fisheries, 1999, 7(2): 237-243.

[3] Tomomi H, Tetsuya T, Mitsuhiro N, et al.Maternal effects and larval survival of marbled solePseudopleuronectes yokohamae[J].Journal of Sea Research, 2007, 58(1): 78-89.

[4] 吳彪, 楊愛國, 劉志鴻, 等.魁蚶兩個不同群體形態(tài)性狀對體質(zhì)量的影響效果分析[J].漁業(yè)科學進展,2010, 31(6): 54-59.

[5] 白臨建, 楊愛國, 周麗青, 等.櫛江珧形態(tài)性狀對重量性狀的影響[J].漁業(yè)科學進展, 2012, 33(6): 87-92.

[6] 李洋, 劉萍, 李健, 等.脊尾白蝦形態(tài)性狀對體重的相關(guān)性及通徑分析[J].漁業(yè)科學進展, 2012, 33(6):60-65.

[7] 邊力, 鐘聲平, 劉洪濤, 等.兩月齡日本囊對蝦形態(tài)性狀對體質(zhì)量的通徑分析[J].廈門大學學報, 2013,52(3): 427-432.

[8] 孫苗苗, 陳百堯, 楊其彬, 等.斑節(jié)對蝦形態(tài)性狀對體質(zhì)量影響的分析[J].海洋科學, 2013, 37(5): 50-54.

[9] 劉賢德, 蔡明夷, 王志勇, 等.不同生長時期大黃魚形態(tài)性狀與體重的相關(guān)性分析[J].熱帶海洋學報,2010, 29(5): 159-163.

[10] 宋堅, 肖登兵, 張偉杰, 等.長蛸體尺與重量性狀間相關(guān)關(guān)系的研究[J].中國農(nóng)學通報, 2012, 28(23):118-122.

[11] 王輝, 劉志剛, 符世偉.南海毛蚶形態(tài)特征對體重的相關(guān)分析[J].熱帶海洋學報, 2007, 26(6): 58-61.

[12] 陳麗梅, 劉利華, 秦傳新, 等.渤海地區(qū)毛蚶形態(tài)性狀對活體重的影響效果[J].安徽農(nóng)業(yè)科學, 2012,40(28): 13813-13814.

[13] 郝偉, 徐國成, 龐作寬, 等.毛蚶殼性狀對體質(zhì)量和軟體部質(zhì)量的影響效果分析[J].淮海工學院學報,2012, 21(3): 81-87.

[14] 杜家菊, 陳志偉.使用 SPSS線性回歸實現(xiàn)通徑分析的方法[J].生物學通報, 2010, 45(2): 4-6.

[15] 劉春雷, 常玉梅, 梁利群, 等.轉(zhuǎn)大麻哈魚生長激素基因鯉表型性狀與體質(zhì)量的相關(guān)性及通徑分析[J].應(yīng)用生態(tài)學報, 2011, 22(7): 1893-1899.

[16] Ahmed M, Abbas G.Growth parameters of finfish and shellfish juveniles in the tidal waters of Bhanbhore′Korangi Creek and Miani Hor Lagoon [J].Pakistan Journal of Zoology, 2000, 32(1): 21-26.

[17] 孫振興, 常林瑞, 徐建鵬.扁玉螺(Neverita didyma)表型性狀對體重和軟體部重的影響效應(yīng)分析[J].海洋與湖沼, 2010, 41(4): 513-518.

[18] Jonah P M, Adeniji O T, Wammanda D T.Genetic correlations and path analysis in bambara groundnut(Vigna subterranea) [J].Journal of Agriculture & Sicial Sciences, 2010, 6: 1-5.

[19] Reza A.Correlations betwwen some horticultural traits in walnut [J].Hort Science, 2010, 45(11): 1690-1694.

[20] 高瑋瑋, 袁媛, 潘寶平, 等.青蛤貝殼形態(tài)性狀對軟體部重的影響分析[J].海洋與湖沼, 2009, 40(2): 166-169.

[21] 朱波, 青長樂, 牟樹森.紫色土外源鋅、鎘形態(tài)的生物有效性[J].應(yīng)用生態(tài)學報, 2012, 13(5): 555-558.

[22] 宋堅, 張偉杰, 常亞青, 等.硬殼蛤形態(tài)性狀對活體重的影響效果分析[J].安徽農(nóng)業(yè)大學學報, 2010,37(2): 273-277.

[23] 劉文廣, 林堅士, 何毛賢.不同貝齡華貴櫛孔扇貝數(shù)量性狀的通徑分析[J].南方水產(chǎn)科學, 2012, 8(1): 43-48.

[24] 鄭懷平, 孫澤偉, 張濤, 等.華貴櫛孔扇貝 1齡貝數(shù)量性狀的相關(guān)性及通徑分析[J].中國農(nóng)學通報, 2009,25(20): 322-326.

[25] 劉志剛, 章啟忠, 王輝.華貴櫛孔扇貝主要經(jīng)濟性狀對閉殼肌重的影響效果分析[J].熱帶海洋學報, 2009,28(1): 61-66.