王成东,聂鸿涛、2,鹿瑶,赵力强,宋亚亚,闫喜武
(1.大连海洋大学辽宁省贝类良种繁育工程技术研究中心,辽宁大连116023;2.大连海洋大学农业部北方海水增养殖重点实验室,辽宁大连116023)
薄片镜蛤Dosinia corrugate隶属于软体动物门、瓣鳃纲、帘蛤目、帘蛤科、镜蛤属,俗称灰蚶子、黑蛤、蛤叉,主要生活在潮间带中、低潮区的泥沙底质中。在中国北起辽宁,南至海南琼山北港,在国外从日本陆奥湾到九州、朝鲜半岛、菲律宾、印度尼西亚均有分布[1]。薄片镜蛤是一种经济价值很高的海产贝类,具有出肉率高、肉味鲜美、营养丰富等优点,深受大众喜爱,发展潜力很大。目前,国内对薄片镜蛤生物学及人工养殖方面的研究较少,仅见闫喜武等[2]对薄片镜蛤的人工育苗试验进行过报道,薄片镜蛤的生物学基础理论亟需完善。
海洋贝类的壳形态和质量性状作为其重要的经济性状,是养殖生产选育苗种过程中的重要指标[3-4]。在选育过程中,主要侧重于质量性状的筛选,但有时质量性状不易测量,如软体部质量需要对贝类进行解剖,这样就对亲贝造成了伤害,在养殖和育种研究中增加了难度,而贝类形态性状的测量相对较为容易。目前,主要采用相关分析和多元回归分析方法研究贝类形态性状与质量性状的关系[3-5],了解壳形态性状对活体质量和软体部质量的影响,便可通过对壳形态性状的选择达到选择育种的目的,这在实际生产中具有重要的应用价值[5]。
利用相关分析和多元回归分析的方法研究贝类的壳形态性状和质量性状的关系,国内外已有较多报道。Ahmed等[5]运用相关分析和通径分析方法研究了贝类幼龄期壳形态性状和质量性状相关的生长参数。在海洋贝类中,李朝霞等[6]通过通径分析对不同日龄海湾扇贝的自交和杂交后代表型性状对质量性状的影响进行了深入的研究。此外,也有关于魁蚶[7]、珍珠贝[8-9]、紫石房蛤[10]、 青蛤[11-12]、硬壳蛤[13]、大竹蛏[14]、鲍[15]、广大扁玉螺[16]、虾夷扇贝[17-18]等形态性状和质量性状相关性和通径分析的报道。但对镜蛤属的相关研究仅见张伟杰等[19]对日本镜蛤有过报道。本研究中,以大连庄河薄片镜蛤野生群体为对象,通过对其数量性状间相关性及通径分析研究,查明各数量性状间的相关性及影响产量的因素组成,旨在为薄片镜蛤的遗传育种提供理论依据。
试验用薄片镜蛤为2013年11月在大连市庄河近海采集的野生薄片镜蛤群体,随机抽取125个个体为样本,取回后在辽宁省贝类良种繁育工程技术研究中心实验室暂养1 d后进行体尺和体质量性状的测量。
1.2.1 薄片镜蛤性状的测量 用电子游标卡尺(精度为0.01 mm)测量薄片镜蛤的壳长 (SL)、壳高 (SH)、壳宽 (SW)、壳顶至壳前 (a)、壳顶至壳后 (b)、韧带长 (c),测量位点见图1。用电子天平 (精度为0.01 g)测量薄片镜蛤的活体质量 (W),将薄片镜蛤解剖后用吸水纸吸去壳内外表面、软体部表面的水分,再测量软体质量(MW)、壳质量 (Sw)。
图1 薄片镜蛤壳型性状测量位点Fig.1 Landmark points of morphological measurement in clam Dosinia corrugate
1.2.2 分析方法 用Excel软件对测得的薄片镜蛤各性状数据进行整理和简单计算,得到各性状的平均值、标准差和变异系数。数据分析参照Huo等[20]的方法,使用 SPSS 17.0软件进行分析,对所测各性状进行相关分析、通径分析,并计算决定系数,采用逐步回归法剔除偏回归系数不显著的性状,保留偏回归系数显著的形态性状对体质量建立回归方程。设置0.05为显著性水平,0.01为极显著性水平。
相关系数 (r)、通径系数 (Pi)、相关指数(R2)和决定系数 (d)的计算公式如下:
其中:式 (2)中,bxi为自变量的偏回归系数,σxi为自变量的标准差,σy为因变量的标准差;式(3)中,rxiy为自变量与因变量之间的相关系数;式 (4)中,di为单一自变量对因变量的决定系数,dij为两个自变量对因变量的共同决定系数,rij为两个自变量之间的相关系数,Pi和Pj分别为两个自变量对因变量的通径系数。
从表1可见:质量性状的变异系数较大,其中以软体部质量的变异系数最大,为30.96%;质量性状的变异系数都高于形态性状,说明质量性状的选择潜力大于表型性状。对各性状数据进行正态性检验,均符合正态分布。
表1 薄片镜蛤各性状的描述统计量 (n=125)Tab.1 Statistics of various traits of clam Dosinia corrugate(n=125)
对各性状进行相关性分析,结果如表2所示。
从表2可见:各性状间均呈极显著相关 (P<0.01),说明对本试验中所选的各性状进行相关性分析是具有实际意义的;从相关度强度[21]来看,大部分性状为高度相关 (r≥0.7),少数性状为中度相关 (0.4<r<0.7);各表型性状中,壳长和壳高与质量性状的相关系数较大,韧带长与质量性状的相关系数最小。
表2 薄片镜蛤各性状之间的相关系数Tab.2 Correlation coefficients among various traits of clam Dosinia corrugate
选择壳形态性状作为自变量,质量性状为因变量,运用SPSS软件进行回归分析,得到对因变量影响显著的各自变量的通径系数 (表3)。根据显著性检验结果,保留通径系数达到显著水平的性状,剔除不显著的性状。在软体质量和壳质量的分析中,剔除了壳宽、韧带长、壳顶至壳前、壳顶至壳后;对总质量进行回归分析时,剔除了韧带长、壳宽、壳顶、壳前,得到薄片镜蛤形态性状对质量性状的最优回归方程:
表3 薄片镜蛤表型性状的偏回归系数检验Tab.3 Multiple regression equations between body weight traits and morphological traits of clam Dosinia corrugate
各方程回归关系均达到了极显著水平 (P<0.01)(表4),说明所得方程成立,可以进行通径分析。经回归预测,估计值与实际观测值差异不显著 (P>0.05),说明上述方程可以真实客观地反映薄片镜蛤形态性状和质量性状的关系。
表4 薄片镜蛤表型性状与质量性状间多元回归方程的方差分析Tab.4 ANOVA analysis of multiple regression equations between body weight traits and morphological traits of clam Dosinia corrugate
通过回归分析,差异不显著的表型性状被剔除,计算得到形态性状对质量性状的通径系数和相关指数 (表5)。从表5可见,壳高对软体质量的直接作用大于壳长,壳长通过壳高对软体质量的间接作用也较大,达到了0.429,可见壳高是影响软体部质量的主要因素。在育种过程中,通过对壳高的选择可以有效地增加软体质量 (MW);壳长和壳高对壳质量的直接影响均较大,且二者通过另一性状对壳质量的间接影响也较大,可见壳长和壳高都是影响壳质量的主要因素;对总质量的直接影响中,壳长最大,为0.648,显著高于其他表型性状(P<0.05),因此,壳长是影响总质量的主要因素,壳高通过壳长的间接作用最大,为0.614。
表5 表型性状对质量性状影响的通径分析Tab.5 Direct and indirect effects of morphological traits on body weight traits of clam Dosinia corrugate
根据式 (4)计算得到表型性状间对质量性状的决定系数 (表6)。从表6可见:壳高对软体质量的相对决定程度大于壳长,但壳高和壳长的共同决定系数最大,为0.335,远大于自身对软体质量的决定系数;形态性状对壳质量的决定程度与对软体质量的决定程度相似,但壳长的相对决定程度大于壳高;形态性状对总质量的决定程度中,壳长的决定系数最大 (0.420),壳顶至壳后的决定系数最小 (0.035),共同决定系数中,壳高和壳长的共同决定系数最大 (0.412),大于其他共同决定系数。
表6 表型性状对质量性状的决定系数Tab.6 Determinant coefficients of morphological traits on body weight traits of clam Dosinia corrugate
各表型性状对软体质量决定系数总和为0.693,对壳质量决定系数总和为0.739,对总质量决定系数总和为0.690,分别等于对应的相关指数R2值。说明找到的形态性状是影响薄片镜蛤软体质量、壳质量和总质量的主要表型性状,其他被剔除的性状对质量性状的影响较小。
在贝类养殖过程中,经济贝类的软体质量和总质量作为主要的经济性状对提高养殖效益具有重要作用。壳质量并非直观的经济性状,但在以往的研究中发现,除年龄、性腺质量、活体质量等因素外,性别、壳质量也是影响滩涂贝类软体质量的重要因素[22]。贝壳在生态学中也有重要意义,刘慧等[23]、张继红等[24]研究发现,在海洋碳汇中贝壳有着重要的作用。所以本试验中选择了薄片镜蛤的软体质量、壳质量和总质量3个质量性状进行分析。袁志发等[25-26]指出,性状间的相关系数是变量间相互关系的综合作用,它包含了两者的直接关系和通过其他变量的间接关系。目前,在选择育种中,绝大多数以增加生长速度、个体质量等数量性状作为育种的主要目标,但由于质量难以准确获取,就要考虑通过其他相关性较高的目标性状加以间接选择,以达到更好的选育效果。本试验中,随机抽取了125个薄片镜蛤,分析各个性状间的相关性及形态性状对质量性状的影响,为养殖选种过程中经济性状的选择提供了理论依据;本试验中确定的质量性状回归方程,还可以解决不易直接测量的质量性状的选择问题。
本试验中,薄片镜蛤的各性状间均存在极显著的表型相关 (P<0.01),这与已有的对其他多种贝类的研究结果一致[27-29]。通过进一步进行的逐步回归分析,剔除了与质量性状相关性不显著的形态性状,求得形态性状对各质量性状的多元回归方程,经检验所得方程可靠,可为薄片镜蛤选育提供理论依据和合理的测量指标。
普遍认为,进行通径分析时,只有当相关指数R2或各自变量对因变量的决定系数的总和大于或等于0.85时,才表明影响因变量的主要自变量已经找到[30]。本试验中,薄片镜蛤的质量性状和形态性状的相关指数R2均未达到0.85,说明本试验中所选的形态性状还不够全面,除试验所选形状外还有其他表型性状对质量性状影响较大。由本试验结果可知,壳形态性状对质量性状的影响主要是壳长和壳高,二者的共同决定程度分别达到了33.5%、35.9%、41.2%。而且本试验中决定系数分析的结果与通径分析结果相同。由此可见,壳长和壳高是决定薄片镜蛤软体部质量、壳质量和总质量的最主要性状,这一结果与其他贝类的研究结果相近。吴彪等[7]研究了魁蚶的7个表型性状对质量性状的影响效果,结果表明,壳长、壳高和壳宽是影响体质量的主要表型性状。吴杨平等[14]的研究结果表明,大竹蛏壳长等9个表型性状中壳长对质量性状的决定程度最大。
本试验中确定了影响薄片镜蛤质量性状的主要表型性状,能够为今后开展薄片镜蛤的选择育种和选育过程提供理论参考。但由于所选表型性状数量有限,其他性状也可能会对质量性状有较大影响,后续研究中应测量更多的表型性状,从而进行更全面地分析,以期为薄片镜蛤的遗传育种研究提供更详尽的参考资料。
[1]庄启谦.中国动物志 软体动物门[M].北京:科学出版社,2001:153-157.
[2]闫喜武,左江鹏,张跃环,等.薄片镜蛤人工育苗技术的初步研究[J].大连水产学院学报,2008,23(4):268-272.
[3]Kvingedal R,Evans B S,Lind C E,et al.Population and family growth response to different rearing location,heritability estimates and genotype×environment interaction in the silver-lip pearl oyster(Pinctada maxima)[J].Aquaculture,2010,304(1):1-6.
[4]Wang H,Chai X,Liu B.Estimation of genetic parameters for growth traits in cultured clam Meretrix meretrix(Bivalvia:Veneridae)using the Bayesian method based on Gibbs sampling[J].Aquaculture Research,2011,42(2):240-247.
[5]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.
[6]李朝霞,王春德.海湾扇贝自交与杂交子代的生长比较和通径分析[J].中国农学报,2009,25(6):282-285.
[7]吴彪,杨爱国,刘志鸿,等.魁蚶两个不同群体形态性状对体重的影响效果分析[J].渔业科学进展,2010,31(6):54-59.
[8]Deng Yuewen,Du Xiaodong,Wang Qingheng,et al.Correlation and path analysis for growth traits in F1population of pearl oyster Pinctada martensii[J].Marine Science Bulletin,2008,10(2):68-73.
[9]栗志民,刘志刚,王辉,等.企鹅珍珠贝主要经济性状对体质量的影响效果分析[J].海洋与湖沼,2011,42(6):798-803.
[10]黎筠,王昭萍,于瑞海,等.紫石房蛤壳性状对活体重影响的定量分析[J].海洋水产研究,2008,29(6):71-77.
[11]高玮玮,袁媛,潘宝平,等.青蛤贝壳形态性状对软体重的影响分析[J].海洋与湖沼,2009,40(2):166-169.
[12]孙同秋,曾海祥,郑小东,等.丹东青蛤野生群体数量性状的相关分析[J].大连海洋大学学报,2013,28(2):171-173.
[13]宋坚,张伟杰,常亚青,等.硬壳蛤形态性状对活体重的影响效果分析[J].安徽农业大学学报,2010,37(2):273-277.
[14]吴杨平,陈爱华,姚国兴,等.大竹蛏表型性状通径和回归分析[J].南京师大学报:自然科学版,2012,35(2):97-102.
[15]Luo X,Ke C H,You W W.Estimates of correlations for shell morphological traits on body weight of interspecific hybrid abalone(Haliotis discus hannai and Haliotis gigantea)[J].Journal of Shellfish Research,2013,32(1):115-118.
[16]Zhao Liqiang,He Yangyang,Yang Feng,et al.Correlation and path analysis of morphological and weight traits in marine gastropod Glossaulax reiniana[J].Chinese Journal of Oceanology and Limnology,2014,324:821-827.
[17]于德良,丁君,郝振林,等.不同养殖群体虾夷扇贝数量性状的相关性与通径分析[J].大连海洋大学学报,2013,28(4):350-354.
[18]赵鹏,丁君,常亚青.两种壳色虾夷扇贝壳体尺性状对活体重影响效果的分析[J].大连海洋大学学报,2011,26(1):1-5.
[19]张伟杰,常亚青,丁君,等.日本镜蛤(Dosinia japonica Reeve)壳尺寸与重量性状的相关与回归分析[J].海洋与湖沼,2013,44(3):796-800.
[20]Huo Z M,Yan X W,Zhao L Q,et al.Effects of shell morphological traits on the weight traits of Manila clam(Ruditapes philippinarum)[J].Acta Ecologica Sinica,2010,30(2):251-256.
[21]袁卫,庞皓,曾五一.统计学[M].北京:高等教育出版社,2000:142-184.
[22]闫喜武,王琰,郭文学,等.四角蛤蜊形态性状对重量性状的影响效果分析[J].水产学报,2011,35(10):1513-1517.
[23]刘慧,唐启升.国际海洋生物碳汇研究进展[J].中国水产科学,2011,18(3):695-702.
[24]张继红,方建光,唐启升.中国浅海贝藻养殖对海洋碳循环的贡献[J].地球科学进展,2005,20(3):359-365.
[25]袁志发,周敬芋.多元统计分析[M].北京:科学出版社,2002:130-l31.
[26]袁志发,周敬芋,郭满才,等.决定系数——通径系数中的决策指标[J].西北农林科技大学学报:自然科学版,2001,29(5):131-133.
[27]刘志刚,王辉,孙小真,等.马氏珠母贝经济性状对体重决定效应分析[J].广东海洋大学学报,2007,27(4):15-20.
[28]张存善,常亚青,曹学彬,等.虾夷扇贝体形性状对软体重和闭壳肌重的影响效果分析[J].水产学报,2009,33(1):87-94.
[29]刘小林,常亚青,相建海,等.栉孔扇贝壳尺寸性状对活体重的影响效果分析[J].海洋与湖沼,2002,33(6):673-678.
[30]刘小林,吴长功,张志怀,等.凡纳对虾形态性状对体重的影响效果分析[J].生态学报,2004,24(4):857-862.