海洋贝类中普遍存在壳色多态性。丰富多彩的壳色品种在受到消费者青睐的同时也提高了商品的经济价值，颜色已成为一种重要的经济性状。大量研究表明，壳色主要由遗传因素决定^[1-2]，并且与生长、存活等性状有密切关联^[3-6]。因此，壳色作为一种可稳定遗传的性状，已经成为贝类遗传育种的一个研究热点。目前，以壳色为目标性状进行纯色系的定向选育，在多种经济贝类中已取得了重要进展，例如：在马氏珠母贝(Pinctada martensii)培育出红色、红褐色、黄色和白色4种壳色选育系^[7]，在文蛤(Meretrix meretrix.)培育出青灰色、淡黄色、灰白色等壳色群体^[8]，在海湾扇贝(Argopecten irradians)有灰色、棕色、黄色、橙色、紫色等壳色群体^[9]。

长牡蛎(Crassostrea gigas)又称太平洋牡蛎，具有适应性强、生长快、肉质细腻且营养丰富等优点，是世界上养殖范围最广、产量最高的经济贝类。2015年中国牡蛎年产量高达457万t，占全国贝类产量的33.7%^[10]，牡蛎养殖在我国海产贝类养殖业中占有重要地位。目前长牡蛎养殖中存在性状退化、品质下降等问题，而培育高产抗逆长牡蛎新品种是解决这一问题的有效途径。因此，利用壳色作为遗传标记开展种质改良，为长牡蛎的遗传育种提供了新方向。经过继代选育，已培育出壳色稳定遗传的金色、黑色、白色和紫色长牡蛎群体(见图 1)，并开展了壳色性状遗传参数评估及与生长性状的关联分析，比较了不同壳色家系子代的表型性^[11-12]，而对壳色性状量化评价的研究报道较少。

本研究以长牡蛎金色、黑色、白色和紫色4种壳色群体为材料，利用国际照明委员会(CIE)L^*、a^*、b^*颜色模型对4种壳色性状进行量化分析，并与普通养殖群体进行比较，旨在为长牡蛎壳色性状判定提供一种有效的方法，为长牡蛎壳色群体选育以及壳色性状分子基础研究提供技术支撑。

1 材料与方法 1.1 实验材料

实验所用的金色、白色、紫色、黑色4种壳色长牡蛎样品的亲本是自2010年开始，经3代连续选育的壳色性状纯化个体。2013年7月，建立4种壳色群体，同时以普通养殖长牡蛎为亲本建立对照组群体。2014年7月，随机选取4种壳色以及普通壳色各38个，将左壳用清水清洗干净后干燥备测，由于右壳面积小，波状层纹密集，颜色不明显，所以右壳颜色未测。

1.2 数据处理及分析

壳色的测量参照王雪磊等^[11]的方法，计算机视觉系统(Computer Vision System, CVS)，该系统由CCD摄像头或照相机、扫描仪、图像采集卡、实验箱、稳压电源、光源以及相应的图像处理软件等组成。相机拍摄时，佳能相机(Canon EOS 60D)左右两侧各设置2个20 W的白炽灯，以黑色塑胶为背景采集图像。光圈值f/5.6，曝光时间1/160 s，焦距92 mm。图像采集完成后，使用Photoshop CS6(Adobe System Incorporated)图像处理软件对采集的牡蛎整个左壳图像进行处理。将每个贝壳颜色参数L、a及b值进行记录。L^*、a^*、b^*依据Chakraborty^[12]提出的方法换算，公式如下：

$ \begin{array}{l} {L^*} = \frac{L}{{255}} \times 100;\\ {a^*} = \frac{{240a}}{{255}} - 120;\\ {b^*} = \frac{{240b}}{{255}} - 120。\end{array} $

式中：L^*表示明度，为正值，L^*越大越接近于白色，达到100视为全白，完全黑则视为0；a^*为红绿轴色品指数，大于0，a^*值越大颜色越偏向红，小于0时，a^*值越小越偏向绿；b^*为黄蓝轴色品指数，大于0时，b^*值越大越偏向黄色，小于0时，b^*值越小越偏向蓝色^[13-14]。

采用16.0SPSS统计软件对数据进行统计分析，将黑、紫、金、白4种壳色群体的L^*、a^*、b^*值分别与普通壳色的颜色参数进行对比及显著性检验，并进行主成分分析。运用逐步判别法对4种壳色进行判别分析，建立4个判别方程，将L^*、a^*、b^*值带入4个公式，可分别计算出Y，观测个体颜色覆盖最多的颜色为Y值最大的判别函数所对应的群体颜色。对4种壳色群体所有个体进行逐个验证，判别准确率的计算公式如下^[15]:

判别准确率P₁(%)=判别正确的个数/实测数×100%；

判别准确率P₂(%)=判别正确的个数/判别数×100%；

判别准确率$ P = \sum\limits_{i = 1}^k {{A_i}/} \sum\limits_{i = 1}^k {{B_i}} $。

2 结果 2.1 长牡蛎4种壳色群体及对照群体壳色参数统计量

长牡蛎各壳色参数：其中，L^*值最高的为白壳群体，平均值为89.04，最低的为黑壳群体，平均值为11.43；a^*值最高的为紫壳群体，平均值为2.89，最低的为白壳群体，平均值为-1.11，b^*值最高的为金壳群体，最低的为黑壳群体，平均值分别为18.4和-3.42。L^*、a^*、b^*在4种壳色群体中的值与普通对照群体均有显著性差异(P < 0.05)。在黑壳、紫壳、白壳、金壳和对照群体中，方差最大的均为L^*，分别为1.58、4.56、3.91、2.02和15.01。

2.2 长牡蛎4种壳色群体及对照群体壳色参数主成分分析

利用长牡蛎4种壳色群体和普通对照群体的L^*、a^*、b^*值进行主成分分析，提取了2个主成分，2个主成分累计贡献率达到91.9%(见表 2)，表明这2个主成分可以解释5个群体间的主要颜色变异。第一主成分解释总变异的53.09%，其中a^*在PC1中负荷值最大，为0.94。第二主成分解释了总变异的38.89%，L^*值在PC2中负荷最大(见表 2)。

基于PC1和PC2两个主成分的散点图显示了4种壳色群体的散点分布与对照群体明显不同。金壳、黑壳、白壳与紫壳群体散点分布均较为集中，对照群体散点分布分散；紫壳与金壳群体散点较黑壳与白壳群体分散(见图 2)。

2.3 长牡蛎4种壳色群体判别分析

运用黑壳、紫壳、金壳和白壳群体的颜色参数，采用逐步判别分析法对4个壳色群体进行判别分析，基于Fisher’s法分别建立了4个壳色选育群体的判别公式，判别参数如表 3。

4个选育群体的判别公式如下：

$\begin{array}{*{20}{l}} {{Y_b} = 1.37{x_1} + 1.49{x_2} - 1.78{x_3} - 12.35;}\\ {{Y_p} = 1.68{x_1} + 4.70{x_2} - 2.14{x_3} - 21.98;}\\ {{Y_g} = 5.11{x_1} - 3.77{x_2} + 2.33{x_3} - 165.28;}\\ {{Y_w} = 9.52{x_1} - 0.79{x_2} - 5.18{x_3} - 434.31}。\end{array} $

式中：x₁，x₂，x₃分别代表L^*、a^*、b^*，对4个壳色群体按上述判别函数进行预测分类，验证判别公式的实用性(见表 4)。黑壳群体的判别准确率P1为100%，P2为86.36%；紫壳群体的判别准确率P1为84.21%，P2为100%；白壳群体与金壳群体判别准确率为100%，综合判别率为96.05%。

3 讨论

长牡蛎壳色一般由色系深浅程度不一的黄、白、紫、黑等多种壳色构成，色泽分布一致的个体较少。传统的颜色判定主要是借助于肉眼的识别进行颜色等级标准的制定，其结果因人而异，容易造成对壳色的判定结果的不准确。本研究采用CIE的L^*、a^*、b^*色度系统是分析颜色的国际标准，利用颜色分析系统研究分析壳色，能使对壳色模糊的描述信息认识变成数据化的客观指标，具有科学性和可比性。目前色度系统主要应用于合浦珠母贝(Pinctada martensi)、三角帆蚌(Hyriopisis cumingii)、马氏珠母贝等产珠贝类的珍珠质颜色的测量^[16-18]。

本研究结果表明，白色壳色L^*值最大，均值达89.04金色；黑色壳色的L^*值最小，平均为11.42。珍珠颜色通过L^*a^*b^*评估表明，白色珍珠L^*值一般大于90，黑色珍珠L^*值一般小于55^[18]。这与我们的实验测量结果相似。本研究中b^*值最高的为金壳群体，这是因为b^*值越大表示越偏向金黄色，L^*值在黑壳群体最小，紫壳群体次之，与L^*值越低颜色越深的判断标准一致。与普通对照群体相比，4种壳色群体的L^*、a^*、b^*值的方差均小于对照组，这表明4种壳色群体的L^*、a^*、b^*值较对照群体波动小，且L^*、a^*、b^*值较对照组均有显著性变化(P < 0.05)，表明壳色性状经过连续多代选育，已经得到了比较稳定的遗传。对马氏珠母贝及三角帆蚌的壳色研究也表明通过选育可以获得壳色的稳定遗传^{[17, 19]}。在对长牡蛎金壳性状颜色参数进行遗传评估中发现，参数L^*遗传力低，h²为0.13，而a^*、b^*均有中高水平遗传力^[11]。本研究中发现，L^*值在每种壳色中的方差比a^*、b^*的方差高，表明每种壳色中L^*较a^*和b^*波动大，这可能与L^*的低遗传力有关。

主成分分析散点分布显示，4种壳色选育系的壳色均较为一致，这也进一步证实长牡蛎壳色性状经过连续多代选育，已经得到了高程度的纯化。而紫壳与金壳群体的散点分布范围较白壳与黑壳群体分散，说明金壳与紫壳群体的壳色性状有进一步纯化的空间。此外，主成分分析显示壳色间差异主要是由参数a^*造成的。通过L^*、a^*、b^*值对3种不同壳色三角帆蚌的颜色进行量化分析，结果也显示参数a^*是引起颜色变化的主要原因^[20]。

判别分析是在分类确定的条件下，根据某一研究对象的各种特征值判别其类型归属的统计分析方法。该方法在鱼类的形态差异分析中通常有较多的应用，例如梁前进等^[21]等对野生鲫(Carassius auratus)和5个金鱼(Carassius auratus)品种进行判别分析，各品种正确率达85%以上；唐瞻杨等^[22]采用判别分析的方法对不同月龄性状的尼罗罗非鱼(Oreochromis niloticus)进行研究，综合判别正确率达到99.25%；窦亚琪等^[23]对翘嘴鳜(Siniperca chuatsi)不同月龄性状进行判别分析，判别正确率达到98.87%。在贝类壳色分析中，判别分析已经应用于三角帆蚌壳色的分析中，研究运用判别分析及L^*、a^*、b^*值测量建立了2个判别方程，可通过a^*, b^*值代入方程，成功将“白贵妃”与“紫皇后”2个品种分开，综合判别正确率达100%^[20]。本研究建立的4种壳色的判别方程中，将四种纯化壳色进行颜色判别分析结果显示，综合验证率为96.05%，判别效果较为理想，说明该判别函数在实际应用时参考价值较大。通过判别分析，可判断未经壳色选育的长牡蛎何种壳色覆盖度最高，有利于长牡蛎的壳色选育。

4 结语

长牡蛎4种壳色群体的L^*、a^*、b^*值与普通对照群体相比均有显著性差异，表明壳色性状经过连续多代选育，已经得到了比较稳定的遗传；主成分分析进一步证明了该结论，同时也显示壳色间差异主要由参数a^*造成的。研究建立了4个壳色的判别方程，壳色判别方程的建立可为长牡蛎壳色的判断提供一种行之有效的方法，为长牡蛎的壳色选育提供基础。