【研究意义】选育具有优良生产性能的新品种或者开发新物种是振兴珍珠产业的一种有效途径。【前人研究进展】合浦珠母贝（Pinctada fucata）、黑珠母贝（Pinctada nigra）和射肋珠母贝（Pinctada radiata）3种珍珠贝同属于珍珠贝目（Pterioida）珍珠贝科（Pteriidae），主要分布于热带和亚热带海区的岩石和珊瑚礁上［1］。3种珍珠贝同为小型珍珠贝，成贝个体间大小相近，形态相似，有着共同的地理分布和生理特征。然而，合浦珠母贝作为四大主养的海洋珍珠贝类之一，与大珠母贝（Pinctada maxima）［2］、珠母贝（Pinctada margaritifera）［3］及企鹅珍珠贝（Pteria penguin）［4］已在世界范围内被广泛进行人工增养殖。自20世纪70年代我国成功突破合浦珠母贝人工繁育以来，该贝已成为重要的水产养殖品种，由合浦珠母贝培育的珍珠占海水珍珠养殖量的95%以上，因此合浦珠母贝的珍珠养殖也曾一度成为部分沿海地区支柱性产业之一［5］。然而，同属于小型珍珠贝的黑珠母贝和射肋珠母贝，与合浦珠母贝的境遇却大不相同，除自然海区有一定分布外，并未得到其他开发和利用。有关黑珠母贝和射肋珠母贝的研究报道也相应较少，仅表现为初步尝试插核育珠实验［6-7］。而对合浦珠母贝的利用和研究则较为透彻，涉及到遗传育种［8-13］、苗种繁育［14］、摄食生理［15］、免疫学［16］及分子生物学［17-19］等各方面。但近些年来，由于合浦珠母贝长期不当的人工近亲繁殖，导致养殖性状退化，加上自然海区污染，使养殖过程中出现的生长缓慢、死亡率高和珍珠优珠率低等现象已严重制约着其产业的发展［20］。生物体的表型性状受基因型、环境以及基因型与环境互作的影响，而在实际的养殖生产中，又多以重量作为生物体的指标来衡量其生产性能［21］。对于合浦珠母贝、黑珠母贝和射肋珠母贝，其湿重、软体重和壳重是衡量其生产性能的重要经济性状。然而，相比经济性状，其形态性状则较为直观，更容易测量。因此发挥形态性状直观可辨的优势，研究其对经济性状的决定和影响程度，并以形态性状作为选育参数，筛选具有优良生产性能的群体，对指导选育工作具有重要意义［22-24］。【本研究切入点】本试验以捕捞的野生合浦珠母贝、黑珠母贝和射肋珠母贝为研究对象，比较3种野生珍珠贝的表型性状，并将其形态性状对经济性状进行多元回归分析。通过相关分析和通径分析，获取影响经济性状的主要形态性状，从而建立有效的回归方程。【拟解决的关键问题】分析其对经济性状的直接作用和间接作用，及其单独或共同对经济性状的决定系数，为珍珠贝新品种的开发及亲贝的挑选工作提供数据参考。
1 材料与方法
1.1 试验材料
3种野生珍珠贝（合浦珠母贝、黑珠母贝和射肋珠母贝）均为2018年5月采集于海南省三亚市西岛附近海域，从总数量超过5 000只的3种野生珍珠贝中挑选体型端正、规格较大的100只个体。
1.2 试验方法
对上述采集的3种野生珍珠贝清除其表面附着物并刷洗干净，用精确度为0.01 mm的数显卡尺测量其壳长（x1，mm）、壳高（x2，mm）、壳宽（x3，mm）和绞合线长（x4，mm），用精确度为0.01 g的电子称测量其湿重（y1，g）、软体重（y2，g）及壳重（y3，g）。3种野生珍珠贝如图1所示，各表型性状的测量方法参考邓正华等［25］。
图1 3种野生珍珠贝
Fig.1 Three species of wild pearl oysters
1.3 数据整理与分析
测量数据用Excel 2016进行处理，获得各性状表型参数统计量。数据采用SPSS 21统计软件进行单因素方差分析（One-Way ANOVA），LSD法多重比较相同性状不同种类间的差异显著性（P＜ 0.05）。分析各性状间的表型相关性，形态性状对经济性状的通径分析及决定系数，剖析各形态性状对经济性状的直接作用和间接作用。剔除偏回归系数不显著的形态性状，利用多元回归分析建立偏回归系数显著的形态性状与经济性状之间的多元回归方程，并进行显著性检验。相关计算公式如下：
含肉率（%）=软体重（g）/湿重（g）×100
决定系数di：di=pi2，式中，pi为形态性状对经济性状的通径系数。
多参数共同决定系数dij：dij=2rij pi pj，式中，rij为两个形态性状间的相关系数，pi、pj分别为单一参数对经济性状的通径系数。
经济性状与各形态性状间的线性回归方程：y = b0 + b1x1 + b2x2 + ... bixi，式中，b0 为常数，b1、b2、...bi为对应的偏回归系数。
2 结果与分析
2.1 3种野生珍珠贝表型性状统计及比较
从3种野生珍珠贝的表型性状统计及比较结果（表1）可见，合浦珠母贝的壳长、壳高、壳宽、绞合线长、湿重、软体重、壳重及含肉率的均值（
± SD）均显著大于其他两种珍珠贝；而射肋珠母贝的壳高、壳宽、绞合线长以及含肉率又显著大于黑珠母贝。3种珍珠贝间各表型性状的极大值（Max）与均值的比较结果相一致，合浦珠母贝的各表型性状均大于其他两种珍珠贝，而射肋珠母贝除绞合线长和壳重的极大值略小于黑珠母贝，其他的表型性状都略大于黑珠母贝。对于变异系数（CV），合浦珠母贝与射肋珠母贝的经济性状均大于其他形态性状；而黑珠母贝的各表型性状变异系数中经济性状大于除壳宽以外的形态性状。在3种野生珍珠贝中，其经济性状的变异系数大小顺序均为软体重＞壳重＞湿重。
表1 3种野生珍珠贝的表型性状统计及比较
Table 1 Statistics and comparison of phenotypic traits among three species of wild pearl oysters
注：同列数据后小写英文字母不同者表示差异显著（P ＜ 0.05）。
Note: Different lowercase letters in the same column represented significant difference（P ＜ 0.05）.
含肉率Tissue rate（%）合浦珠母贝P. fucata种类Species 壳长Shell length（mm）壳高Shell height（mm）壳宽Shell width（mm）绞合线长Hinge length（mm）湿重Wet weight（g）软体重Tissue weight（g）壳重Shell weight（g）± SD 69.32±6.21a 70.11±5.30a 22.25±1.78a 59.97±5.14a 30.23±6.11a 8.59±2.15a 21.60±4.48a 28.38±3.59a Max 90.14 85.41 28.07 76.36 52.24 17.44 37.61 40.00 CV 8.96 7.56 7.99 8.57 20.22 25.03 20.75 12.66黑珠母贝P. nigra± SD 60.16±5.82b 61.30±5.41c 15.01±2.96c 45.47±5.58c 20.43±3.26b 5.14±1.02b 15.27±2.51b 25.13±2.72c Max 75.69 74.23 19.95 61.13 29.60 7.93 23.15 32.00 CV 9.67 8.83 19.72 12.28 15.98 19.83 16.45 10.83射肋珠母贝P. radiata± SD 60.93±5.72b 62.78±4.77b 17.19±1.73b 47.07±4.52b 20.13±3.63b 5.24±1.05b 14.82±2.81b 26.08±2.88b Max 77.05 81.72 21.99 57.08 29.85 8.85 22.76 35.00 CV 9.39 7.59 10.07 9.60 18.03 19.97 18.95 11.03
2.2 3种野生珍珠贝各表型性状间的相关性分析
3种野生珍珠贝各表型性状间的相关系数及相关显著性结果见表2，其中除射肋珠母贝的壳高与壳宽间的相关性不显著外，其他性状间均表现为极显著水平。在相关系数中，3种野生珍珠贝的经济性状间的相关系数较大，其中湿重与壳重间的相关系数最大，均大于0.96；其次是湿重与软体重间的相关系数，均大于0.83。在3种珍珠贝中与经济性状相关系数最大的形态性状有所不同。其中，对于合浦珠母贝，与湿重和软体重相关系数最大的形态性状是壳宽，与壳重最大的是壳高。对于黑珠母贝，与湿重和壳重相关系数最大的形态性状是壳高，与软体重最大的是壳长。对于射肋珠母贝，与湿重和壳重相关系数最大的形态性状是壳长，与软体重最大的是绞合线长。单一形态性状与经济性状间的相关系数表示该性状对经济性状的总影响，其中除该性状的直接影响外，还包括其他性状间的影响，因此仅通过相关系数不能确定变量间的重要性，还需对相关系数作进一步分析。
表2 3种野生珍珠贝表型性状的相关系数
Table 2 Correlation coefficients of the phenotypic traits of three species of wild pearl oysters
注：“**”表示极显著相关（P ＜ 0.01）；“*”表示显著相关（P ＜ 0.05）。
Notes: “**” means extremely significant correlation（P ＜ 0.01）; “*” means significant correlation（P ＜ 0.05）.
壳重Shell weitht（g）合浦珠母贝P. fucata种类Species参数Parameters壳高Shell height（mm）壳宽Shell width（mm）绞合线长Hinge length（mm）湿重Wet weight（g）软体重Tissue weight（g）壳长 0.686** 0.650** 0.608** 0.601** 0.561** 0.575**壳高 0.569** 0.412** 0.636** 0.546** 0.622**壳宽 0.601** 0.426** 0.444** 0.393**绞合线长 0.503** 0.541** 0.436**湿重 0.841** 0.972**软体重 0.705**射肋珠母贝P.radiata壳长 0.703** 0.519** 0.573** 0.566** 0.438** 0.564**壳高 0.490** 0.564** 0.621** 0.427** 0.650**壳宽 0.465** 0.625** 0.519** 0.607**绞合线长 0.586** 0.479** 0.569**湿重 0.843** 0.965**软体重 0.676**黑珠母贝P. nigra壳长 0.682** 0.246** 0.455** 0.719** 0.518** 0.739**壳高 0.091 0.378** 0.636** 0.387** 0.686**壳宽 0.454** 0.536** 0.541** 0.500**绞合线长 0.639** 0.569** 0.627**湿重 0.837** 0.976**软体重 0.717**
2.3 3种野生珍珠贝形态性状对经济性状的通径分析
根据形态性状与经济性状间的相关系数，使用通径分析理论分解性状间的相关系数，并采用逐步回归的方法得出壳长、壳高、壳宽和绞合线长对经济性状的直接通径系数和其他性状间的间接作用之和。在相关系数中扣除其他形态性状与经济性状的直接作用，得到其他性状中每一个形态性状对经济性状的间接作用。
2.3.1 形态性状对湿重的通径分析 3种珍珠贝形态性状对湿重的通径分析结果（表3）表明：（1）对合浦珠母贝，被保留下来的形态性状有壳高、壳宽和绞合线长，其中对湿重直接影响最大的是壳宽（0.362），壳高通过壳宽对湿重的间接影响最大（0.177）；（2）对于黑珠母贝，被保留下来的形态性状为壳高和绞合线长，其中对湿重直接影响最大的形态性状为壳高（0.517），绞合线长通过壳高对湿重的间接影响最大（0.213）；（3）对于射肋珠母贝，所有形态性状均被保留下来，其中直接影响湿重最大的形态性状为壳长（0.337），壳高通过壳长对湿重的间接影响最大（0.230）。
表3 3种野生珍珠贝各形态性状对湿重的通径分析
Table 3 The path analysis on the morphometric traits to the wet weight of three species of wild pearl oysters
绞合线长Hinge length合浦珠母贝P. fucata种类Species参数Parameters相关系数rij Relative coefficient直接作用Pi Direct effect∑ 壳长Shell length壳高Shell height壳宽Shell width壳长 0.566**壳高 0.621** 0.305** 0.316 0.177 0.139壳宽 0.625** 0.362** 0.263 0.149 0.114绞合线长 0.586** 0.246** 0.340 0.172 0.168黑珠母贝P. nigra壳长 0.601**壳高 0.636** 0.517** 0.119 0.119壳长 0.719** 0.337** 0.381 0.198 0.079 0.104壳高 0.636** 0.291** 0.346 0.230 0.029 0.087壳宽 0.536** 0.323** 0.213 0.083 0.026 0.104绞合线长 0.639** 0.229** 0.410 0.153 0.110 0.147壳宽 0.426**绞合线长 0.503** 0.290** 0.213 0.213射肋珠母贝P. radiata
注：“**”表示极显著相关（P ＜ 0.01）；“*”表示显著相关（P ＜ 0.05）。以第一行为例，∑代表壳长通过壳高、壳宽和绞合线长对湿重的间接作用之和，如此类推。
Notes: “**” means extremely significant correlation（P ＜ 0.01）; “*” means significant correlation（P ＜ 0.05）. Take the first row as an example, ∑ represents the indirect effects of the shell length on the wet weight through the shell height, shell width and hinge length constant, which are the same below.
2.3.2 形态性状对软体重的通径分析 3种野生珍珠贝的形态性状对软体重的通径分析结果（表4）表明：（1）对于合浦珠母贝，被保留下来的形态性状有壳宽和绞合线长，其中直接影响软体重最大的形态性状是壳宽（0.379），绞合线长通过壳宽对软体重的间接影响最大（0.176）；（2）对于黑珠母贝，被保留下来的形态性状为壳高和绞合线长，其中直接影响软体重最大的形态性状为壳高（0.388），绞合线长通过壳高对软体重的间接影响最大（0.160）；（3）对于射肋珠母贝，被保留下来的形态性状为壳长、壳宽和绞合线长，其中直接影响软体重最大的形态性状为壳宽（0.340），绞合线长通过壳宽对软体重的间接影响最大（0.154）。
表4 3种野生珍珠贝各形态性状对软体重的通径分析
Table 4 The path analysis on the morphometric traits to the tissue weight of three species of wild pearl oysters
注：“**”表示极显著相关（P ＜ 0.01）；“*”表示显著相关（P ＜ 0.05）。以第一行为例，∑代表壳长通过壳高、壳宽和绞合线长对软体重的间接作用之和，如此类推。
Notes: “**” means extremely significant correlation（P ＜ 0.01）; “*” means significant correlation（P ＜ 0.05）. Take the first row as an example, ∑ represents the indirect effects of the shell length on the tissue weight through the shell height, shell width and hinge length constant, which are the same below.
绞合线长Hinge length合浦珠母贝P. fucata种类Species参数Parameters相关系数rij Relative coefficient直接作用Pi Direct effect∑ 壳长Shell length壳高Shell heigh壳宽Shell width壳长 0.438**壳高 0.427**壳宽 0.519** 0.379** 0.141 0.141绞合线长 0.479** 0.303** 0.176 0.176黑珠母贝P. nigra壳长 0.561**壳高 0.546** 0.388** 0.157 0.157绞合线长 0.541** 0.381** 0.160 0.160壳长 0.518** 0.309** 0.209 0.084 0.125壳高 0.387**壳宽 0.541** 0.340** 0.200 0.076 0.124绞合线长 0.569** 0.274** 0.295 0.141 0.154壳宽 0.444**射肋珠母贝P. radiata
2.3.3 形态性状对壳重的通径分析 3种野生珍珠贝的形态性状对壳重的通径分析结果（表5）表明：（1）对于合浦珠母贝，被保留下来的形态性状有壳高、壳宽和绞合线长，其中对壳重直接影响最大的形态性状为壳高（0.373），绞合线长通过壳高对壳重的间接影响最大（0.211）；（2）对于黑珠母贝，被保留下来的形态性状为壳长和壳高，其中对壳重直接影响最大的形态形状为壳高（0.430），壳长通过壳高（0.295）对壳重的间接影响最大；（3）对于射肋珠母贝，所有形态性状均被保留下来，其中对壳重直接影响最大的形态性状为壳高（0.355），壳长通过壳高对壳重的间接影响最大（0.242）。
表5 3种野生珍珠贝各形态性状对壳重的通径分析
Table 5 The path analysis on the morphometric traits to the shell weight of three species of wild pearl oysters
注：“**”表示极显著相关（P ＜ 0.01）；“*”表示显著相关（P ＜ 0.05）。以第一行为例，∑代表壳长通过壳高、壳宽和绞合线长对壳重的间接作用之和，如此类推。
Notes: “**” means extremely significant correlation（P ＜ 0.01）; “*” means significant correlation（P ＜ 0.05）. Take the first row as an example, ∑ represents the indirect effects of the shell length on the shell weight through the shell height, shell width and hinge length constant, which are the same below.
绞合线长Hinge length合浦珠母贝P. fucata种类Species参数Parameters相关系数rij Relative coefficient直接作用Pi Direct dffect∑ 壳长Shell length壳高Shell height壳宽Shell width壳长 0.739** 0.329** 0.410 0.242 0.072 0.096壳高 0.686** 0.355** 0.330 0.224 0.026 0.080壳宽 0.500** 0.291** 0.209 0.081 0.032 0.096绞合线长 0.627** 0.211** 0.416 0.150 0.134 0.132壳长 0.575** 0.279** 0.295 0.295壳高 0.622** 0.430** 0.191 0.191壳宽 0.393**绞合线长 0.436**射肋珠母贝P. radiata壳长 0.564**壳高 0.650** 0.373** 0.277 0.161 0.116壳宽 0.607** 0.329** 0.278 0.183 0.095绞合线长 0.569** 0.205* 0.364 0.211 0.153黑珠母贝P. nigra
2.4 3种野生珍珠贝形态性状对经济性状的决定程度
根据公式计算形态性状对经济性状的决定系数。表6～表8对角线上的数据是单一形态性状的决定系数，对角线以上的数据是一个性状通过其他性状间接决定经济性状的系数。
2.4.1 形态性状对湿重的决定程度 3种野生珍珠贝各形态性状对湿重的决定系数（表6）表明：（1）对于合浦珠母贝，单一性状决定系数最大的是壳宽（0.131），壳高通过壳宽的间接决定系数最大（0.108）；（2）对于黑珠母贝，单一性状决定系数最大的是壳高（0.267），壳高通过绞合线长的间接决定系数最大（0.123）；（3）对于射肋珠母贝，单一性状决定系数最大的是壳长（0.114），壳长通过壳高的间接决定系数最大（0.133）。
表6 3种野生珍珠贝各形态性状对湿重的决定系数
Table 6 The determinant coefficients of morphometric traits to the wet weight of three species of wild pearl oysters
种类Species 参数Parameters绞合线长Hinge length合浦珠母贝P. fucata壳长Shell length壳高Shell height壳宽Shell width壳长壳高 0.267 0.123壳宽绞合线长 0.084射肋珠母贝P. radiata壳长壳高 0.093 0.108 0.085壳宽 0.131 0.083绞合线长 0.061黑珠母贝P. nigra壳长 0.114 0.133 0.053 0.070壳高 0.085 0.017 0.051壳宽 0.104 0.067绞合线长 0.052
2.4.2 形态性状对软体重的决定程度 3种野生珍珠贝各形态性状对软体重的决定系数（表7）表明：（1）对于合浦珠珠母贝，单一性状对软体重决定系数最大的形态性状为壳宽（0.144），壳宽通过绞合线长的间接决定系数最大（0.107）；（2）对于黑珠母贝，单一性状对软体重决定系数最大的形态性状为壳高（0.151），壳高通过绞合线长的间接决定系数最大（0.122）；（3）对于射肋珠母贝，单一性状对软体重决定系数最大的形态性状为壳宽（0.116），壳宽通过绞合线长的间接决定系数最大（0.084）。
表7 3种野生珍珠贝各形态性状对软体重的决定系数
Table 7 The determinant coefficients of morphometric traits to the tissue weight of three species of wild pearl oysters
绞合线长Hinge length合浦珠母贝P. fucata种类Species参数Parameters壳长Shell length壳高Shell height壳宽Shell width壳长壳高 0.151 0.122壳宽绞合线长 0.145射肋珠母贝P. radiata壳长壳高壳宽 0.144 0.107绞合线长 0.092黑珠母贝P. nigra壳长 0.095 0.052 0.077壳高壳宽 0.116 0.084绞合线长 0.075
表8 3种野生珍珠贝各形态性状对壳重的决定系数
Table 8 The determinant coefficients of morphometric traits to the shell weight of three species of wild pearl oysters
种类Species 参数Parameters绞合线长Hinge length合浦珠母贝P. fucata壳长Shell length壳高Shell height壳宽Shell width壳长 0.078 0.165壳高 0.185壳宽绞合线长射肋珠母贝P. radiata壳长壳高 0.139 0.120 0.087壳宽 0.108 0.063绞合线长 0.042黑珠母贝P. nigra壳长 0.108 0.159 0.047 0.063壳高 0.126 0.018 0.057壳宽 0.085 0.056绞合线长 0.045
2.4.3 形态性状对壳重的决定程度 3种野生珍珠贝各形态性状对软体重的决定系数（表8）表明：（1）对于合浦珠珠母贝，单一形态性状对壳重决定系数最大的形态性状为壳高（0.139），壳高通过壳的间接决定系数最大（0.120）；（2）对于黑珠母贝，单一形态性状对壳重决定系数最大的是壳高（0.185），壳长通过壳高的间接决定系数最大（0.165）；（3）对于射肋珠母贝，单一形态性状对壳重决定系数最大的是壳高（0.126），壳长通过壳高的间接决定系数最大（0.159）。
2.5 3种野生珍珠贝形态性状对经济性状的回归方程
2.5.1 形态性状对湿重的回归分析 3种野生珍珠贝被保留下来的形态性状对湿重的偏回归系数统计结果见表9。结果表明，各形态性状对湿重均产生极显著影响，说明保留的形态性状与湿重之间有极显著的线性关系。对3种野生珍珠贝建立以形态性状为自变量、湿重为因变量的回归方程：
表9 3种野生珍珠贝形态性状对湿重的偏回归系数和回归常数的显著性检验
Table 9 Significance test of the partial regression coefficients and constants of morphometric traits to the wet weight of three species of wild pearl oysters
参数 
偏回归系数Partial regression coefficient t P合浦珠母贝P. fucata常数 -6.389 2.881 -2.218 0.029壳高 0.312 0.049 6.402 0.000绞合线长 0.169 0.047 3.588 0.001射肋珠母贝P. radiata常数 -39.626 6.352 -6.239 0.000壳高 0.351 0.100 3.528 0.001壳宽 1.244 0.277 4.496 0.000绞合线长 0.293 0.101 2.895 0.005黑珠母贝P. nigra常数 -27.092 3.001 -9.029 0.000壳长 0.213 0.047 4.512 0.000壳高 0.222 0.055 4.044 0.000壳宽 0.678 0.123 5.496 0.000绞合线长 0.184 0.051 3.570 0.001
2.5.2 形态性状对软体重的回归分析 3种野生珍珠贝被保留下来的形态性状对软体重的偏回归系数统计结果见表10。结果表明，各形态性状对软体重均产生极显著的影响，说明保留的形态性状与软体重之间有极显著的线性关系。对3种野生珍珠贝建立以形态性状为自变量、软体重为因变量的回归方程：
表10 3种野生珍珠贝形态性状对软体重的偏回归系数和回归常数的显著性检验
Table 10 Significance test of the partial regression coefficients and constants of morphometric traits to the tissue weight of three species of wild pearl oysters
偏回归系数Partial regression coefficien t P合浦珠母贝P. fucata参数Parameter
常数 -2.512 0.947 -2.653 0.009壳高 0.073 0.016 4.571 0.000绞合线长 0.070 0.016 4.488 0.000射肋珠母贝P. radiata常数 -9.201 2.513 -3.661 0.000壳宽 0.458 0.112 4.069 0.000绞合线长 0.127 0.039 3.259 0.002黑珠母贝P. nigra常数 -4.739 1.024 -4.627 0.000壳长 0.057 0.015 3.814 0.000壳宽 0.206 0.049 4.206 0.000绞合线长 0.063 0.020 3.105 0.002
2.5.3 形态性状对壳重的回归分析 3种野生珍珠贝被保留下来的形态性状对壳重的偏回归系数统计结果见表11。结果表明，除合浦珠母贝的绞合线长对壳重产生显著影响外，其他各形态性状对壳重均产生极显著的影响，说明保留的形态性状与壳重之间有显著的线性关系。对3种野生珍珠贝建立以形态性状为自变量、壳重为因变量的回归方程：
表11 3种野生珍珠贝形态性状对壳重的偏回归系数和回归常数的显著性检验
Table 11 Significance test of the partial regression coefficientsand constants of morphometric traits to the shell weight of three species of wild pearl oysters
参数 
偏回归系数Partial regression coefficient t P合浦珠母贝P. fucata常数 -29.684 4.661 -6.368 0.000壳高 0.315 0.073 4.313 0.000壳宽 0.829 0.203 4.084 0.000绞合线长 0.179 0.074 2.412 0.018黑珠母贝P. nigra常数 -4.233 2.296 -1.843 0.068壳长 0.121 0.046 2.647 0.009壳高 0.200 0.049 4.076 0.000射肋珠母贝P. radiata常数 -22.457 2.235 -10.049 0.000壳长 0.162 0.035 4.588 0.000壳高 0.209 0.041 5.130 0.000壳宽 0.472 0.092 5.135 0.000绞合线长 0.131 0.038 3.423 0.001
多元线性回归的方差分析结果表明，3种野生珍珠贝的形态性状与经济性状之间的回归均达到极显著水平。其中3种野生珍珠贝的形态性状与湿重间的方差分析结果为，合浦珠母贝：F=40.670 ＞ F0.01（3，96）；黑珠母贝：F= 43.837 ＞F0.01（2，97）；射肋珠母贝：F= 69.920 ＞ F0.01（4，95）。3种野生珍珠贝的形态性状与软体重间的方差分析结果为，合浦珠母贝：F= 25.186 ＞ F0.01（2，97）；黑珠母贝：F= 34.904 ＞ F0.01（2，97）；射肋珠母贝：F= 32.086 ＞ F0.01（3，96）。3 种野生珍珠贝的形态性状与壳重间的方差分析结果为，合浦珠母贝：F= 40.520 ＞ F0.01（3，96）；黑珠母贝F= 36.312 ＞ F0.01（2，97）；射肋珠母贝：F=77.354 ＞ F0.01（4，95）。以上方差分析结果均符合P ＜ 0.01，方差分析结果存在极显著差异。
3 讨论
3种野生珍珠贝中，合浦珠母贝各表型性状的均值和极大值都大于黑珠母贝和射肋珠母贝；而射肋珠母贝的壳高、壳宽、绞合线长及含肉率显著大于黑珠母贝。正是由于合浦珠母贝相较于黑珠母贝和射肋珠母贝所具有的优良生产性状，从而能够在世界范围内被广泛接受并进行大规模人工养殖的重要原因之一。生物体的表型性状是受基因型、环境以及基因型与环境共同作用的影响［21］，本研究中，3种野生珍珠贝同处于共同的环境之下，因此可以排除环境因素影响对其表型性状分析所造成的差异，三者的表型性状差异可以认为是由自身基因型所决定的。喻达辉等［26］曾对珠母贝属的合浦珠母贝、黑珠母贝、射肋珠母贝、大珠母贝、珠母贝、企鹅珍珠贝等其他珍珠贝类进行遗传距离的分析，发现相比于其他珍珠贝，合浦珠母贝、黑珠母贝以及射肋珠母贝三者被分为一支，其遗传距离最小，证明三者的亲缘关系较近。 从表型性状分析，3种野生珍珠贝经济性状的变异系数基本都大于其形态性状，结果与栉孔扇贝（Chlamys farreri）［21］、文蛤（Meretrix meretrix）［23］、长肋日月贝（Amusium pleuronects）［25］、香港巨牡蛎（Crassostrea hongkongensis）［27］等相同，而在湿重、软体重及壳重间，三者又以软体重的变异系数最大。王雨等［28］认为，成熟的贝类一般四季都可以繁殖，排卵（精）受潮水周期、气候变化以及个体生理状况差异的影响，因此软体重变异系数较大。壳宽是珍珠贝插核育珠过程中一个极其重要的形态性状，3种珍珠贝中以黑珠母贝的壳宽均值最小，然而其变异系数是最大的，因此黑珠母贝在以壳宽为选育目标时还有着巨大的提高潜力。
相关系数是变量间相互关系的综合体现，3种野生珍珠贝中，除射肋珍珠贝的壳高与壳宽没有显著相关外，其他种类的表型性状间均表现出极显著的相关水平，这也表明射肋珠母贝同种类个体间的体型差异较大。本研究中，3种野生珍珠贝的形态性状与经济性状间的关系均达到极显著水平，但并不能因此而说明所有的形态性状是影响其经济性状的主要因素，因此有必要通过通径分析，从而找到实际影响其湿重、软体重及壳重的主要形态性状。在该研究中，经逐步回归剔除法，只有相关系数水平达到显著水平才被保留下来，对于3种珍珠贝中与经济性状间被保留的形态性状在种间有所差异。形态性状对湿重的通径分析中，对于合浦珠母贝，被保留的形态性状为壳高、壳宽和绞合线长，直接影响最大的是壳宽。对于黑珠母贝，被保留的形态性状为壳高和绞合线长，直接影响作用最大的是壳高。而对于射肋珠母贝，所有形态性状均被保留下来，直接影响作用最大的是壳长。可见，不同种类的珍珠贝影响其湿重作用最大的形态性状有所不同，3种野生珍珠贝中，不同种类的珍珠贝影响其软体重和壳重的作用最大的形态性状也不尽相同。不仅在不同种类生物间表现出有此差异，在同种类生物的野生群体与养殖群体间［20］、不同生长阶段［21-22］、不同壳色间［23］、不同地理群体间［24-25］及雌雄个体间［27］都存在着一定差异。在相关分析的基础上进行通径分析时，只有当复相关系数或各自变量对因变量的单独决定系数及两两共同决定系数的R2≥0.850时，才表明影响因变量的主要变量已经找到。本研究中，3种野生珍珠贝被保留下来的形态性状与湿重、软体重和壳重的复相关系数均小于0.85。因此，说明影响3种野生珍珠贝湿重、软体重和壳重的形态性状除了被保留的性状外还存在其他性状。
4 结论
综上，对于3种野生珍珠贝，合浦珠母贝各表型性状均大于其他两种珍珠贝。3种珍珠贝间影响其湿重最大的形态性状有所不同，影响其软体重最大的形态性状在合浦珠母贝与射肋珠母贝间相同但与黑珠母贝不同，影响其壳重最大的形态性状在三者间均相同。因此，在对特定物种、特定种群进行人工选育时，应根据具体被选择的对象来进行特定的分析，以此来确定最佳的测量指标。
［1］杨文，蔡英亚，邝雪梅. 中国南海经济贝类原色图谱［M］. 北京:中国农业出版社，2013：168-171.YANG W, CAI Y Y, KUANG X M. Color atlas of economic mollusca from the South China Sea［M］. Beijing: China Agriculture Press,2013:168-171.
［2］MAMANGKEY N G F, SOUTHGATE P C. Regemeration of excised mantle tissue by the silver-lip pearl oyster, Pinctada maxima（Jameson）［J］. Fish & Shellfish Immunology, 2009, 27: 164- 174. doi: 10.1016/j.fsi.2009.03.009
［3］邓正华，李海梅，张博，苏家齐，刘宝锁，范嗣刚，周凤，吴开畅，喻达辉.珠母贝幼虫运动、自然沉降及摄食行为［J］. 中国水产科学，2018，25（1）：116-123. doi: 10.3724/SP.J.1118.2018.17044.DENG Z H, LI H M, ZHANG B, SU J Q, LIU B S, FAN S G,ZHOU F, WU K C, YU D H. Swimming, natural sedimentation and ingestion of pearl oyster Pinctada margaritifera larvae［J］. Journal of Fishery Sciences of China, 2018, 25（1）: 116-123.doi: 10.3724/SP.J.1118.2018.17044.
［4］SOUTHGATE P C, BEER A C, NGALUAFE P. Hatchery culture of the winged pearl oyster, Pteria penguin, without living microalgae［J］. Aquaculture, 2016, 45: 121-124. doi: 10.1016/j.aquaculture.2015.09.007.
［5］牛志凯，刘宝锁，张东玲，谭才钢，张博，陈明强，范嗣刚，姜松，黄桂菊，李有宁，喻达辉. 合浦珠母贝3个地理群体杂交后代生长性状和闭壳肌拉力的比较分析［J］. 南方水产科学，2015，11（1）：26-32. doi:10.3969/j.issn.2095-0780.2015.01.004.NIU Z K, LIU B S, ZHANG D L, TAN C G, ZHANG B, CHEN M Q, FAN S G, JIANG S, HUANG G J, LI Y N, YU D H. Comparative analysis of growth traits and shell-closing strength among hybrid populations from three geographical grouos of pearl oyster（Pinctada fucata）［J］. South China Fisheries Science, 2015, 11（1）: 26-32. doi:10.3969/j.issn.2095-0780.2015.01.004.
［6］蒙钊美，张勤奋，李有宁，邢孔武. 黑珠母贝的核位解剖及其插核育珠初步研究［J］. 中国水产科学，1994，（1）：25-30.MENG Z M, ZHANG Q F, LI Y N, XING K W. A study on anatomy of sites of nucleus-insertion and pearl culture in Pinctada nigra（Gould）［J］. Journal of Fishery Sciences of China, 1994（1）: 25-30.
［7］李有宁，邢孔武，陈明强. 射肋珠母贝插核试验［J］. 上海水产大学学报，1998，7（2）：162-166.LI Y N, XING K W, CHEN M Q. A test for nucleus-insertion of Pinctada radiata［J］. Journal of Shanghai University, 1998, 7（2）: 162-166.
［8］王学颖，高远镇，杜晓东，邓岳文，卢婧. 马氏珠母贝金黄壳色F3和基础群体遗传结构比较 ［J］. 海洋通报, 2012, 31（3）：324-328.WANG X Y, GAO Y Z, DU X D, DENG Y W, LU J. Genetic structure of the third generation yellow colored line and base stock of pearl oyster Pinctada martensii as revealed by SSR marker system［J］. Marine Science Bulletion, 2012,31（3）: 324-328.
［9］陈静，刘志刚，孙小真，王辉，杜晓东. 马氏珠母贝4种壳色选育系F3的生长特性和遗传多样性比较［J］. 上海海洋大学学报，2010, 19（5）：588-595.CHEN J, LIU Z G, SUN X Z, WANG H, DU X D. Analysis on growth characteristic and genetic diversity in the selected lines（F3） of four shell colors of Pinctada martensii［J］. Journal of Shanghai Ocean University, 2010, 19（5）: 588-595.
［10］王庆恒，邓岳文，杜晓东，符韶. 马氏珠母贝4个壳色选育系F1幼虫的生长比较［J］. 中国水产科学，2008，15（3）：488-492.WANG Q H, DENG Y W, DU X D, FU S. Growth comparision of F1 larvae between four shell color lines of pearl oyster Pinctada martensii［J］. Journal of Fishey Sciences of China, 2008, 15（3）: 488-492.
［11］陈明强，刘宝锁，严俊贤，郭华阳，吴开畅，李有宁. 合浦珠母贝4种壳色选育系主要性状的比较分析［J］. 南方水产科学，2016, 12（5）：118-122. doi: 10. 3969/ j. issn. 2095 - 0780. 2016. 05. 015.CHEN M Q, LIU B S, YAN J X, GUO H Y, WU K C, LI Y N.Comparsion analysis of major traits among four shell color selective lines of pearl oyster（Pinctada fucata）［J］. South China Fisheries Science, 2016, 12（5）: 118-122. doi: 10. 3969/ j. issn. 2095 - 0780.2016. 05. 015.
［12］吴曼，刘宝锁，黄桂菊，陈明强，郭奕惠，陈素文，喻达辉. 合浦珠母贝全同胞家系贝壳珍珠质颜色分析［J］. 南方水产科学，2014，10（6）：44-50. doi: 10. 3969 /j. issn. 2095-0780. 2014. 06. 006.WU M, LIU B S, HUANG G J, CHEN M Q, GUO Y H, CHEN S W, YU D H. Analysis of shell nacre color in full-sib families of pearl oyster（Pinctada fucata）［J］. South China Fishery Science, 2014, 10（6）: 44-50. doi: 10. 3969 /j. issn. 2095-0780. 2014. 06. 006.
［13］MENG Z H, ZHANG B, LIU B S, LI H M, FAN S G, YU D H. High carotenoids content can enhance resistance of selected Pinctada fucata families to high temperature stress［J］. Fish & Shellfish Immunology,2017,61: 211-218. doi: org/10.1016/j.fsi.2016.12.032.
［14］WANG H, ZHU X W, WANG Y A, LUO M M, LIU Z G. Determination of optimum temperation and salinity for fertilization and hatching in the Chinese pearl oyster Pinctada martensii（Dunker）. Aquaculture, 2012,358-359: 292-297. doi: 10.1016/j.aquaculture.2012.04.050.
［15］邓正华，姜松，张博，刘宝锁，黄桂菊，喻达辉. 合浦珠母贝对不同种类及浓度的单胞藻摄食与消化效果研究［J］. 南方水产科学, 2016,12（3）: 112-118. doi: 10. 3969/ j. issn. 2095 - 0780. 2016. 03. 015.DENG Z H, JIANG S, ZHANG B, LIU B S, HUANG G J, YU D H.Ingestion and digestion of pearl oyster（Pinctada fucata） on microalgae of different types and concentrations［J］. South China Fishery Science,2016, 12（3）: 112-118.doi: 10. 3969/ j. issn. 2095 - 0780. 2016. 03. 015.
［16］WEI J F, FAN S G, LIU B S, ZHANG B, SU J Q, YU D H. Transeriptome analysis of the immune reaction of the pearl oyster Pinctada fucata to xenograft from Pinctada maxima［J］. Fish & Shellfish Immunology,2017, 67: 331-345.doi: org/10.1016/j.fsi.2017.06.030.
［17］LI H M, WANG D Q, DENG Z H, HUANG G J, FAN S G, ZHOU D Z, LIU B S, ZHANG B. YU D H. Molecular characterization and expression analysis of chitinase from the pearl oyster Pinctada fucata［J］. Comparative Biochemistry and physiology Part B, 2017, 203:141-148. doi: org/10.1016/j.cbpb.2016.10.007.
［18］GUAN Y Y, HE M X, WU H B. Differential mantle transcriptomics and characterization of growth-related genes in the diploid and triploid pearl oyster Pinctada fucata［J］. Marine Genomics, 2017.33: 31-38.doi: org/10.1016/j.margen.2017.01.001.
［19］范嗣刚，周代志，刘宝锁，邓正华，郭奕惠，喻达辉. 合浦珠母贝BMP7b基因的克隆与表达分析［J］. 南方水产科学，2018，14（1）：121-126. doi: 10.3969/j.issn.2095 - 0780.2018.01.016.FAN S G, ZHOU D Z, LIU B S, DENG Z H, GUO Y H, YU D H.Molecular cloning and expression analysis of BMP7b from Pinctada fucata［J］. South China Fishery Science, 2018, 14（1）: 121-126.doi:10.3969/j.issn.2095 -0780.2018.01.016.
［20］刘志刚，王辉，孙小真，杜晓东，李雷斌. 马氏珠母贝经济性状对体重决定效应分析［J］. 广东海洋大学学报，2007，27（4）：15-20.LIU Z G, WANG H, SUN X Z, DU X D, LI L B. Analysis of determination effects of economic traits on body weight of Pinctada martensii （Dunker）［J］. Journal of Guangdong Ocean University, 2007,27（4）: 15-20.
［21］杜美荣，方建光，高亚平，房景辉，蒋增杰. 不同贝龄栉孔扇贝数量性状的相关性和通径分析［J］. 水产学报，2017，41（4）：580-587.doi: 10.11964/jfc.20160610463.DU M R, FANG J G, GAO Y P, FANG J H, JIANG Z J. Correlation and path analysis of quantitative traits of different-age Chlamys farreri ［J］.Journal of Fisheries of China, 2017, 41（4）: 580-587.doi: 10.11964/jfc.20160610463.
［22］边力，刘长琳，陈四清，张乐乐，葛建龙，刘克奉，谭杰. 不同生长期圆斑星鲽形态性状对体重影响的通径分析［J］. 中国水产科学，2017，24（6）：1168-1175. doi: 10.3724/sp.j.1118.2017.17144.BIAN L, LIU C L, CHEN S Q, ZHANG L L, GE J L. LIU K F, TAN J.Path analysis of effects of morphometric traits on body weight in spotted halibut Verasper variegatus at different growth stages［J］. Journal of Fishery Sciences of China, 2017, 24（6）: 1168-1175. doi: 10.3724/sp.j.1118.2017.17144.
［23］方军，肖国强，张炯明，柴雪良，滕爽爽. 两种壳色文蛤壳形态性状对活体质量的影响［J］. 大连海洋大学学报，2017，32（3）：310-315. doi:10.16535 /j.cnki.dlhyxb.2017.03.009.FANG J, XIAO G Q, ZHANG J M, CAI X L, TENG S S. Effect of morphological traits on live body weight in clam Meretrix meretrix with two kind color shells［J］. Journal of Dalian Ocean University,2017,32（3）: 310-315. doi:10.16535 /j.cnki.dlhyxb.2017.03.009.
［24］黄有辉，范斌，李一鸣，刘志权，杨洋，吴东蕾，范武江，赵云龙. 日本沼虾五群体形态性状对体质量的通径分析［J］. 水产学报，2016,40（8）: 1173-1185. doi: 10.11964/jfc.20150709986.HUANG Y H, FAN B, LI Y M, LIU Z Q, YANG Y, WU D L, FAN W J,ZHAO Y L. Path analysis for the correlations between morphorogical traits and body weight in five Macrobrachium nipponense populations［J］. Journal of Fisheries of China, 2016, 40（8）: 1173-1185. doi:10.11964/jfc.20150709986.
［25］邓正华, 陈明强, 李有宁, 吴开畅, 王雨. 5个地区长肋日月贝形态性状对主要经济性状的相关性和通径分析［J］. 中国渔业质量与标准，2018，8（2）：17-28. doi:10.3969/j.issn.2095-1833.2018.02.003.DENG Z H, CHEN M Q, LI Y N, WU K C, WANG Y. Correlation and path analysis between the morphological traits and the main economic traits of Amusium pleuronects from 5 regions［J］. Chinese Fishery Quanlity and Standards, 2018,8（2）: 17-28. doi:10.3969/j.issn.2095-1833.2018.02.003.
［26］喻达辉，朱嘉濠. 珠母贝属的系统发育：核rDNA ITS序列证据［J］.生物多样性，2005，13（4）：315-323. doi: 10.1360/biodiv.050050.YU D H, ZHU J H. Phylogenetics of the common pearl oysters in the genus Pinctada: evidence from nrDNA ITS sequence［J］. Biodiversity Science, 2005, 13（4）: 315-323. doi: 10.1360/biodiv.050050.
［27］肖述，符政君，喻子牛. 香港巨牡蛎雌雄群体的数量性状通径分析［J］. 南方水产科学，2011，7（4）：1-9. doi: 10. 3969 /j. issn. 2095-0780. 2011. 04. 001.XIAO S, FU Z J, YU Z N. Path analysis of quantitative traits of male and female Hong Kong oyster Crassostrea hongkongensis［J］. South China Fisheries Science, 2011, 7（4）: 1-9.doi: 10. 3969 /j. issn. 2095-0780. 2011. 04. 001.
［28］王雨，叶乐，陈旭，杨其彬，温为庚，吴开畅. 海南野生长肋日月贝形态性状与重量性状的通径分析［J］. 安徽农业科学，2009，37（8）：3570-3572. doi :10.13989/j .cnki .0517-6611.2009.08.155.WANG Y, YE L, CHEN X, YANG Q B, WEN W G, WU K C. Path analysis on the morphological and weight characters of wild Amusium pleuronectes in Hainan［J］. Journal of Anhui Agriculture, 2009,37（8）:3570-3572.doi :10.13989/j .cnki .0517-6611.2009.08.155.