2. 上海海洋大学水产与生命学院 上海 201306
2. College of Fisheries and Life Science, Shanghai Ocean University, Shanghai 201306
在水产动物育种中,通过对不同性状间的相关性进行分析可为育种方案提供必要数据。体质量和形态性状具有直观性和可度量性,是鱼类遗传育种最直接且重要的定向选择培育的指标,但体质量性状的测量受鱼肠道残留饲料、体表残留海水等因素,致使测量结果出现偏差,而形态性状在测量上准确度高,便于统一规格,因此,形态性状与体质量的关系研究尤为必要。通径分析是多元回归分析的扩展,通径系数又具有相关系数与偏回归系数的特点(Akintunde, 2012),所以通径分析可以通过通径系数来证明自变量和因变量之间的直接作用、间接作用以及作用的大小,通径分析因具有准确、直观等特点,被广泛应用于大菱鲆(Scophthalmus maximus) (王新安等, 2008)、翘嘴鳜(Siniperca chuatsi) (董浚键等, 2018)、脊尾白虾(Exopalamon carincauda) (李洋等, 2012)、斧文蛤(Meretrix lamarchii) (张嘉丽等, 2014)等多种水产动物育种。20世纪80年代,邓聚龙(1982)提出灰色关联度分析,这种方法具有分析样本的需求量少、得到信息量大、方法简便等优点,可以将众多测量指标综合成能够反应整体的数据,即灰色关联度。通过比较关联度的大小,来确定自变量对自变量、自变量对因变量的重要性,已广泛应用于种植业、畜牧业的新品系筛选和性状相关性研究中(孙峰成等, 2012; 韦春波等, 2010)。
许氏平鲉(Sebastes schlegelii),俗称黑鲪,隶属鲉形目(Scorpaeniformes)、鲉科(Scorpaenidae)、平鲉属(Sebastes),为冷温性近海底层卵胎生鱼类,在中国黄、渤海以及韩国和日本沿海都有较广泛的分布(Feng et al, 2014)。因具有适温范围广、生长快、抗逆性强、营养价值高等特点,已成为我国北方沿海资源增殖和网箱养殖的理想品种。开展人工选育工作可保证许氏平鲉养殖业的健康发展。目前,有关许氏平鲉选择育种的研究主要集中在家系构建(姜海滨等, 2014)、分子标记开发(Ma et al, 2015)、群体遗传结构分析(刘阳等, 2018)等。目前,采用通径和灰色关联度分析对许氏平鲉性状相关性研究未见报道。本研究采用通径和灰色关联度分析2种方法对许氏平鲉三代选育群体的11个形态性状与体质量的相对关系进行评价,为下一步许氏平鲉选育提供合理科学的目标评价指标性状,以期提高选育效率。
1 材料与方法 1.1 实验材料本课题组从2005~2017年在烟台泰华海洋科技有限公司采用群体选育和家系选育相结合的方法陆续开展了许氏平鲉的三代选育工作,2005年利用荣成优良野生群体通过控温控光技术在室内工厂化自由交尾的方式,2006年建立群体内选育子1代(F1),2010年和2016年通过自由交尾将待产雌鱼分池单独培育构建半同胞家系,参照姜海滨等(2014),分别于次年建立群体内的家系选育子2代(F2)、子3代(F3),2017年集体产仔,在相同的培育环境下,将仔鱼培育至8月龄,分别从三代选育群体中随机取样80尾测量12个性状指标:体质量(Y)、全长(X1)、体长(X2)、体高(X3)、头长(X4)、尾长(X5)、尾柄长(X6)、尾柄高(X7)、吻长(X8)、躯干长(X9)、眼径(X10)和眼后头长(X11)。用电子天平称量体质量,精确至0.01 g;其他形态性状采用放置刻度尺拍照,并用Image J软件测量,精确至0.01 cm。
1.2 分析方法 1.2.1 通径分析使用Excel 2010整理各形态性状和体质量的测量数据,统计各性状的最大值、最小值、平均值、标准差、变异系数等表型参数。利用软件SPSS 17.0在各性状间进行相关性分析,在表型性状相关的基础上,利用逐步法进行形态性状为自变量对体质量的多元回归分析,建立多元回归方程;通过标准化回归系数法计算自变量对因变量的直接通径系数,进一步分解,获得间接通径系数。结合所得各形态性状之间的相关系数,进而求得单性状的直接决定系数、两两性状间共同决定系数,计算方法参照Wright(1921),具体公式如下:
$ {间接通径系数\text{ = }}{r_{{x_i}{x_j}}} \times {P_{{x_j}}} $ | (1) |
$ 决定系数{d_{{x_i}}} = P_{{x_i}}^2 $ | (2) |
$ 共同决定系数{d_{{x_i}{x_j}}} = 2{r_{{x_i}{x_j}}} \times {P_{{x_i}}} \times {P_{{x_j}}} $ | (3) |
式中,
由于许氏平鲉形态性状和体质量的量纲不同,性状间无法进行直接比较。因此,本研究以体质量为参考序列X0,11个形态性状为比较序列Xi,i=1, 2, 3, …11,采用标准差法对原始数据进行无量纲化处理,参照公式(4)。根据灰色理论系统,将11个形态性状和体质量视为一个灰色系统,参照公式(5)和(6)计算灰色关联系数和关联度,最后依据灰色关联度的大小来评价体质量与各形态性状的密切程度,从而判断各形态性状对体质量的重要性,具体计算方法参照邓聚龙(1982)。相关公式如下:
$ {X'_i}(k) = \frac{{{X_i}(k) - {{\bar X}_i}}}{{{S_i}}} $ | (4) |
$ {\delta _i}(k){\text{ }} = \frac{{min\left| {{{X'}_0}(k) - {{X'}_i}(k)} \right| + \rho max\left| {{{X'}_0}(k) - {{X'}_i}(k)} \right|}}{{\left| {{{X'}_0}(k) - {{X'}_i}(k)} \right| + \rho max\left| {{{X'}_0}(k) - {{X'}_i}(k)} \right|}} $ | (5) |
$ {\gamma _i} = \frac{1}{n}\mathop \sum \limits_{k = 1}^n {\delta _i}(k) $ | (6) |
式中,
许氏平鲉三代选育群体的各性状表型参数值如表 1所示,三代选育群体中,体质量的变异系数均最大,分别为22.71%、24.28%和18.84%;体长的变异系数均最小,分别为9.32%、9.70%和7.85%。各形态性状间及与体质量的相关系数见表 2,其形态性状间及与体质量间的相关系数达到极显著水平(P < 0.01)或显著水平(P < 0.05),三代选育群体中全长和体长的相关性均最大,相关系数分别为0.985、0.984和0.972,体高对体质量的相关性最高,眼径与体质量之间的相关性均最低。
根据通径分析原理计算各形态性状对体质量的通径系数,结果见表 3。经显著性检验分析,剔除通径系数不显著的表型性状,F1代和F2代均保留了达到极显著水平的体高和躯干长2个性状,所得通径系数分别为0.662、0.336和0.756、0.245;F3代保留了极显著水平的体高和尾柄长2个性状,通径系数为0.683、0.317。三代选育群体的体高对体质量的直接作用均大于间接作用,可见体高是影响体质量的主要因素。而躯干长(F1、F2)和尾柄长(F3)2个性状对体质量的影响均是间接作用大于直接作用。各形态性状对体质量的决定系数见表 4。三代选育群体性状单独的决定系数和两性状两两共同决定系数之和分别为0.936、0.943和0.893,均大于0.85,表明影响许氏平鲉体质量的主要形态性状均被纳入研究中。
采用逐步回归的方法,剔除对体质量影响不显著的形态性状,建立许氏平鲉三代选育群体多元回归方程,依次为:
F1代:Y=−47.383+17.097X3+7.308X9
F2代:Y=−34.615+15.779X3+3.901X9
F3代:Y=−54.496+21.966X3+7.283X6
式中,Y为体质量(g),X3、X6和X9分别为体高(cm)、尾柄长(cm)和躯干长(cm)。对多元回归关系和各偏回归系数进行显著性检验(表 5),均达到极显著水平(P=0.000 < 0.001),说明建立的多元回归方程模拟效果较为准确。对所构建的多元回归方程采用F-检验(表 6),回归方程检验结果均达到显著水平(F1=248.991, P < 0.001; F2=157.772, P < 0.001; F3=72.631, P < 0.001),说明本研究建立的许氏平鲉三代选育群体多元回归方程可以估算个体的体质量,并进行实际应用。
许氏平鲉三代选育群体的形态性状与体质量的灰色关联系数和灰色关联度见表 7、表 8,灰色关联系数均值依次为0.745~0.866、0.764~0.866、0.705~ 0.793。体质量与形态性状的关联度大小排序依次为:体高 > 全长 > 体长 > 躯干长 > 尾柄高 > 尾长 > 头长 > 眼后头长 > 吻长 > 尾柄长 > 眼径(F1);体高 > 体长 > 全长 > 躯干长 > 尾长 > 尾柄高 > 眼后头长 > 头长 > 尾柄长 > 吻长 > 眼径(F2);体高 > 体长 > 全长 > 尾柄高 > 尾长 > 躯干长 > 尾柄长 > 眼径 > 头长 > 吻长 > 眼后头长(F3)。由此可知,三代选育群体的形态性状与体质量的关联程度存在一定差异,但三代的体质量与形态性状的关联度大小排在前3位的性状均为体高、全长、体长,且体高与体质量的关联度均最高(0.8660、0.8663和0.7932),说明许氏平鲉三代选育群体的体高这一性状对体质量影响最大。
体质量是水产动物遗传改良的直接目标性状,当体质量的遗传力下降时,通过对体质量的筛选难以实现遗传选育的预期目标,此时需要通过寻找与体质量相关性较高的其他形态性状进行间接筛选(Okamoto et al, 2006)。相关分析结果显示,各形态性状与体质量存在极显著(P < 0.01)或显著(P < 0.05)的表型相关关系,但表型相关未剔除其他变量的影响,所以不能准确反映自变量与因变量之间的关系(Yücel, 2004)。通径分析可以将性状相关剖分为直接作用和间接作用,从而探明性状间真实的关系。本研究通过通径分析和多元回归分析剔除了不显著的表型性状,进一步建立最优的多元回归方程,三代选育群体均保留极显著的2个形态性状,F1和F2代为体高和躯干长,F3代为体高和尾柄长。韩慧宗等(2016)对12和15月龄许氏平鲉形态性状对体质量的影响时发现,体长和体高是影响体质量的主要形态性状,而本研究中影响体质量的主要形态性状有所不同,可能是不同月龄影响体质量的形态性状存在差异,但两研究均存在体高这一相同的形态性状。刘贤德等(2008)研究13月龄和20月龄的大黄鱼形态性状对体质量影响发现体高对体质量的直接影响最大。白晓慧等(2015)研究斑石鲷(Oplegnathus punctatus)发现,全长对体质量的直接影响最大,而本研究与前者研究结果相同,与后者不同,可能是鱼类体型、种类等差异导致。进行决定程度分析时,当自变量对依变量的单独决定系数及两两共同决定系数的总和∑d≥0.85时,说明影响依变量的自变量的主要变量已经找到(李洋等, 2012)。本研究中,三代选育群体入选的形态性状对体质量的总决定系数分别为0.936、0.943和0.893,说明入选的形态性状是影响体质量的主要性状,而其他剔除或未测定性状对体质量影响较小。
3.2 许氏平鲉三代选育群体各形态性状与体质量的灰色关联度分析本研究首次将灰色关联度分析应用到许氏平鲉11个形态性状与体质量的相关性分析中。水产动物各形态性状与体质量之间存在灰色性,运用灰色系统理论可以探明它们之间相互依存关系,关联度大小表示与性状间的密切程度。谭才钢等(2015)利用灰色关联度分析得到的结果与Deng等(2008)和陈飞飞等(2012)利用相关分析、多元回归分析和通径分析得到的结果相类似,表明灰色关联度分析方法与其他相关分析、通径分析等统计学方法相比,将其运用到水产动物中进行性状间相关性分析是可行的。对于水产动物形态性状与体质量的相关性研究已有较多报道,如鲤鱼杂交后代(苏胜彦等, 2011)、合浦珠母贝(Pinctada martensii) (谭才钢等, 2015)、牙鲆(Paralichthys olivaceus) (刘永新等, 2014)等。刘峰等(2017)对4~5月龄小黄鱼(Larimichthys polyactis)灰色关联度分析发现,雌性、雄性小黄鱼形态性状与体质量的影响次序有所不同,均是全长和体长这2个形态性状对体质量的关联度最高,本研究结果与之类似。由于灰色关联度分析法计算简单扼要,不需要符合某个理论分布,能较好分析小样本数据,能对研究对象多个性状进行评估,有较强的可比性和可靠性,能找到间接选育性状。
3.3 2种分析方法比较及选育性状的确定本研究选取许氏平鲉三代选育群体为研究对象,通过通径分析和灰色关联度分析对形态性状与体质量的相关性进行研究,一方面通过两种方法的分析结果去查究影响三代选育群体体质量的主要形态性状的差异性,另一方面比较两种分析方法在许氏平鲉不同世代选育群体的适用性,最终确定目标选育形态性状。结果显示,前一种方法均得到2个形态性状是影响体质量的主要性状,分别为体高、躯干长或体高、尾柄长,后一种方法得到与体质量关联度最大的2个性状为体高、全长或体高、体长,这2种方法得到的结果相差较大,但排在第一位的主控性状均是体高,这与王凯等(2008)利用通径分析和刘永新等(2014)利用灰色关联度分析2种分析方法对牙鲆体质量与形态性状的相关性的研究结果相类似。原因是通径分析是在考虑自变量间相互关系的基础上进行的,而灰色关联度分析是将自变量独立对待,未考虑自变量的显著与否,因而2种方法在确定自变量的相对重要性上内涵是不同的(乔有明, 1997)。在农业上,灰色关联度分析经常与通径分析比较使用,例如陈剑锋等(2011)对小麦主要性状与产量、蔡树美等(2014)对西瓜(Citrullus lanatus)产量与养分含量进行通径与灰关联分析。在水产选择育种上,利用这2种方法进行比较利用未见报道。本研究结果说明,某一种统计分析方法并不是确定最优选育性状的方法,每一种统计分析方法因其应用原理不同有相应的特点,得出的结果也并不完全一致。所以,在今后确定选育性状工作时,需要针对实际情况选择合适的统计学方法,或利用2种以上的分析方法进行比较使用,确定最理想的选育性状。
本研究利用通径分析和灰色关联度分析2种方法对许氏平鲉形态性状对体质量的影响程度进行了研究,并建立了许氏平鲉三代选育群体主要形态性状与体质量的最优回归方程,确定体高可作为许氏平鲉选育的重要选育目标性状。本研究结果可为下一步许氏平鲉的选育提供科学合理的选育评价指标,以提高选育效率。
Akintunde AN. Path analysis step by step using Excel. Journal of Technical Science and Technologies, 2012, 1(1): 9-15 |
Bai XH, You HZ, Zhang W, et al. Path analysis of effects of morphometric traits on body weight of Oplegnathus punctatus. Journal of Economic Animal, 2015, 19(4): 191-194 [白晓慧, 尤宏争, 张韦, 等. 斑石鲷形态性状对体质量影响的通径分析. 经济动物学报, 2015, 19(4): 191-194] |
Cai SM, Zhu HT, Yu TY, et al. Grey correlation analysis and path coefficient analysis between nutrient content and yield of watermelon. Soils, 2014, 46(4): 651-655 [蔡树美, 诸海焘, 余廷园, 等. 西瓜产量与养分含量的灰色关联度和通径分析. 土壤, 2014, 46(4): 651-655] |
Chen FF, Huang GJ, Chen MQ, et al. Correlation and path analysis of growth-related traits in Sanya cultured population of pearl oyster Pinctada fucata. Guangdong Agricultural Sciences, 2012, 39(9): 122-125 [陈飞飞, 黄桂菊, 陈明强, 等. 合浦珠母贝三亚养殖群体生长性状的相关与通径分析. 广东农业科学, 2012, 39(9): 122-125 DOI:10.3969/j.issn.1004-874X.2012.09.040] |
Chen JF, Chen DY, Zhang QY. Path coefficient analysis and grey correlation analysis between main characters and yield in wheat. Fujian Science and Technology of Rice and Wheat, 2011, 29(4): 5-8 [陈剑锋, 陈登耀, 张秋英. 小麦主要性状与产量的通径和灰色关联度分析. 福建稻麦科技, 2011, 29(4): 5-8 DOI:10.3969/j.issn.1008-9799.2011.04.002] |
Deng JL. The grey control system. Journal of Huazhong (Central China) University of Science and Technology, 1982, 10(3): 9-18 [邓聚龙. 灰色控制系统. 华中工学院学报, 1982, 10(3): 9-18] |
Deng YW, Du XD, Wang QH, et al. Correlation and path analysis for growth traits in F1 population of pearl oyster Pinctada martensii. Marine Science Bulletin, 2008, 10(2): 68-73 |
Dong JJ, Sun CF, Tian YY, et al. Correlation analysis of the main morphological traits and body weight of mandarin fish (Siniperca chuatsi) and morphological traits between males and females. Progress in Fishery Sciences, 2018, 39(2): 76-84 [董浚键, 孙成飞, 田园园, 等. 翘嘴鳜主要形态性状与体重的相关性及雌雄形态性状差异分析. 渔业科学进展, 2018, 39(2): 76-84] |
Feng JR, Liu LM, Jiang HB, et al. Histological observation of germ cell development and discovery of spermatophores in ovoviviparous black rockfish (Sebastes schlegeli, Hilgendorf) in reproductive season. Journal of Ocean University of China, 2014, 13(5): 830-836 DOI:10.1007/s11802-014-2265-6 |
Han HZ, Jiang HB, Wang F, et al. Principal component and path analysis of morphological traits of selective groups at different month ages of Sebastes schlegelii. Journal of Fisheries of China, 2016, 40(8): 1163-1172 [韩慧宗, 姜海滨, 王斐, 等. 许氏平鲉不同月龄选育群体形态性状的主成分与通径分析. 水产学报, 2016, 40(8): 1163-1172] |
Jiang HB, Ma HT, Liu LJ, et al. A preliminary study of superior half-sib family selection of Sebastes schlegeli. Marine Sciences, 2014, 38(10): 70-75 [姜海滨, 马海涛, 刘丽娟, 等. 黑鲪半同胞家系选育的初步研究. 海洋科学, 2014, 38(10): 70-75 DOI:10.11759/hykx20130814001] |
Li Y, Liu P, Li J, et al. Correlation and path analysis of morphometric traits on body weight for Exopalamon carincauda. Progress in Fishery Sciences, 2012, 33(6): 59-65 [李洋, 刘萍, 李健, 等. 脊尾白虾形态性状对体重的相关性及通径分析. 渔业科学进展, 2012, 33(6): 59-65 DOI:10.3969/j.issn.1000-7075.2012.06.009] |
Liu F, Lou B, Chen RY, et al. Analysis of grey relationship between morphological traits and body weight in the small yellow croaker (Pseudosciaena polyactis). Journal of Shanghai Ocean University, 2017, 26(1): 131-137 [刘峰, 楼宝, 陈睿毅, 等. 小黄鱼形态性状与体质量的灰色关联分析. 上海海洋大学学报, 2017, 26(1): 131-137] |
Liu XD, Cai MY, Wang ZY, et al. The correlation and path analysis for growth-related traits of large yellow croaker Pseudosciaena crocea from Min-Yuedong Tribe. Periodical of Ocean University of China, 2008, 38(6): 916-920 [刘贤德, 蔡明夷, 王志勇, 等. 闽-粤东族大黄鱼生长性状的相关与通径分析. 中国海洋大学学报:自然科学版, 2008, 38(6): 916-920] |
Liu Y, Han HZ, Wang TT, et al. Analysis of growth performances and genetic characteristics of the fast-growing breeding populations of Sebastes schlegelii. Oceanologia et Limnologia Sinica, 2018, 49(2): 432-443 [刘阳, 韩慧宗, 王腾腾, 等. 许氏平鲉(Sebastes schlegelii)速生选育群体的生长性能与遗传特征分析. 海洋与湖沼, 2018, 49(2): 432-443] |
Liu YX, Liu YJ, Zhou J, et al. Grey relational analysis between main growth traits and body weight in Japanese flounder (Paralichthys olivaceus). Journal of Fishery Sciences of China, 2014, 21(2): 205-213 [刘永新, 刘英杰, 周勤, 等. 牙鲆主要生长性状与体质量的灰色关联度分析. 中国水产科学, 2014, 21(2): 205-213] |
Ma HT, Jia CF, Yang JM, et al. Development of novel microsatellite markers in the Korean rockfish Sebastes schlegeli. Genetics and Molecular Research, 2015, 14(2): 5099-5102 DOI:10.4238/2015.May.12.12 |
Okamoto C, Komaru A, Hayashi M, et al. Variation of shell-closing strength among several families in pearl oyster, Pinctada fucata martensii. Aquaculture Science, 2006, 54: 525-529 |
Qiao YM. Comparison of path-coefficent analysis and grey system relevant analysis. Journal of Animal Husbandry and Veterinary in Qinghai, 1997, 27(4): 22-25 [乔有明. 通径分析与灰色关联分析的比较. 青海畜牧兽医杂志, 1997, 27(4): 22-25] |
Su SY, Dong ZJ, Qu JQ, et al. Grey correlated degree analysis of growth trait of hybrids from 3 common carp (Cyprinus carpio L.) varieties and prediction on body weight of offspring from multiple crosses. Journal of Fisheries of China, 2011, 35(1): 20-26 [苏胜彦, 董在杰, 曲疆奇, 等. 3个鲤群体杂交后代生长性状的灰色关联及复合杂交后代的体重预测分析. 水产学报, 2011, 35(1): 20-26] |
Sun FC, Feng Y, Yu Z, et al. Grey relativity analysis on main agronomic characters of 12 maize populations with their yields and traits. Acta Agriculturae Boreali-Sinica, 2012, 27(1): 102-105 [孙峰成, 冯勇, 于卓, 等. 12个玉米群体的主要农艺性状与产量、品质的灰色关联度分析. 华北农学报, 2012, 27(1): 102-105 DOI:10.3969/j.issn.1000-7091.2012.01.019] |
Tan CG, Liu BS, Zhang DL, et al. Analysis of grey relationship between morphological traits and body weight of pearl oyster (Pinctada fucata). South China Fisheries Science, 2015, 11(2): 35-40 [谭才钢, 刘宝锁, 张东玲, 等. 合浦珠母贝主要形态性状与体质量的灰色关联分析. 南方水产科学, 2015, 11(2): 35-40 DOI:10.3969/j.issn.2095-0780.2015.02.005] |
Wang K, Liu HJ, Liu YX, et al. Mathematical analysis of effects of morphometric attributions on body weight for Paralichthys olivaceus. Journal of Shanghai Fisheries University, 2008, 17(6): 655-660 [王凯, 刘海金, 刘永新, 等. 牙鲆形态性状对体重的影响效果分析. 上海水产大学学报, 2008, 17(6): 655-660] |
Wang XA, Ma AJ, Xu K, et al. Relationship between morphometric attributes and body weight of juvenile turbots Scophthalmus maximus. Acta Zoologica Sinica, 2008, 54(3): 540-545 [王新安, 马爱军, 许可, 等. 大菱鲆幼鱼表型形态性状与体重之间的关系. 动物学报, 2008, 54(3): 540-545] |
Wei CB, Liu SS, Jia YQ, et al. The grey correlational analysis of DHI index and production performance of dairy cow. Journal of Heilongjiang Bayi Agricultural University, 2010, 22(1): 106-109 [韦春波, 李甡屾, 贾永全, 等. 奶牛DHI测定指标与其生产性能的灰色关联分析. 黑龙江八一农垦大学学报, 2010, 22(1): 106-109 DOI:10.3969/j.issn.1002-2090.2010.01.028] |
Wright S. Correlation and causation. Journal of Agricultural Research, 1921, 20(7): 557-585 |
Yücel C. Correlation and path coefficient analyses of seed yield components in the Narbon Bean (Vicia narbonensis L.). ). Turkish Journal of Agriculture and Forestry, 2004, 28(5): 371-376 |
Zhang JL, Wang QH, Deng YW, et al. Effects of morphological traits on the weight of Meretrix lamarchii. Progress in Fishery Sciences, 2014, 35(6): 110-113 [张嘉丽, 王庆恒, 邓岳文, 等. 斧文蛤(Meretrix lamarchii)形态性状对体质量的影响效果分析. 渔业科学进展, 2014, 35(6): 110-113] |