刺参又名仿刺参(Apostichopus japonicus)，广泛分布于太平洋西部的中国、俄罗斯、韩国与日本沿海，含有丰富的蛋白质，粘性多糖及多种微量元素，具有很高的药用价值和营养保健价值，是我国海水养殖中重要的经济品种(廖玉麟, 1997)。自20世纪80年代刺参人工繁育技术突破以后，刺参养殖产业快速发展，截止2017年，全国刺参养殖面积达21.9万hm²，养殖产量21.9万t，直接产值逾300亿元(中国渔业统计年鉴, 2018)。然而，随着产业的迅猛发展，刺参养殖产业出现生长速度慢、病害频发、养殖成活率低等种质退化现象。而在韩国和俄罗斯等刺参自然分布海域，由于过度捕捞和环境污染的加剧，野生刺参自然资源趋于枯竭，根据世界自然保护联盟的数据，近30~50年以来，野生刺参的生物资源量，日本下降了30%，韩国下降了40%，俄罗斯下降了80%，中国下降了95%以上。刺参已被世界自然保护联盟收录到濒危物种红色名录的濒危Endangered(EN)等级(Purcell et al, 2013)。因此，开展刺参种质资源评估和遗传结构分析，不仅可以为刺参良种选育提供科学数据，还可为客观评估刺参的种质资源现状提供技术支撑。

目前，用于种质资源评估的分子标记主要包括微卫星DNA和线粒体DNA等。其中，线粒体DNA(mtDNA)是细胞核外具有自主复制、转录和翻译能力的遗传因子，与核DNA相比，具有分子小、结构简单、演化速度快和母系遗传等优点，是用于进行种群遗传结构分析的理想标记。目前，在水产动物中，用于遗传结构解析的mtDNA片段主要包括16S rDNA、D-loop和COⅠ基因。本研究采集中国、韩国和俄罗斯沿海6个海域8个不同地理群体的刺参，分别采用16S rDNA、D-loop和COⅠ这3个片段对所采集的地理群体种质资源和遗传结构进行评估，以期为科学评估目前刺参种质资源现状提供数据支撑，并为种质资源保护及开发利用提供理论依据。

实验用野生刺参分别采集于中国青岛、烟台、俄罗斯符拉迪沃斯托克、韩国浦项、韩国木浦和韩国群山共6个海域，根据所采集刺参体色的差异，韩国浦项采集的刺参又分为红参、黄参和黑参3个群体，每个群体随机选取10头刺参作为实验刺参，合计采集6个海域8个不同地理群体刺参80头。采集的刺参采用冷冻方式运回实验室，剖取纵肌用于DNA的提取。不同群体的采样地点及信息见图 1和表 1。

图 1(Fig. 1)

图 1 刺参采样地点 Fig.1 Locations of sea cucumber A. japonicus sampling

表 1(Tab.1)

表 1 不同地理种群刺参样本采样信息 Tab.1 Information of A. japonicus sampling 群体名称 Population name 采样时间 Collection date 位置 Location 样本量 Number 青岛太平角海参QD 2017-05-26 36°02'N, 120°21'E 10 韩国浦项黑参SK-PX-B 2017-06-09 36°05'N, 129°27'E 10 韩国浦项红参SK-PX-R 2017-06-09 36°05'N, 129°27'E 10 韩国浦项黄参SK-PX-Y 2017-06-10 36°05'N, 129°27'E 10 韩国群山黑参SK-QS-B 2017-06-09 36°34'N, 126°23'E 10 韩国木浦黑参SK-MP-B 2017-06-09 34°07'N, 126°17'E 10 烟台长岛海参YT 2017-07-03 38°17'N, 120°49'E 10 俄罗斯海参RU 2018-04-19 35°39'N, 119°49'E 10

表 1 不同地理种群刺参样本采样信息 Tab.1 Information of A. japonicus sampling

利用Omega Mollusc DNA kit提取海参基因组DNA，利用NanoDrop 1000及1%琼脂糖凝胶电泳检测DNA纯度及完整性后，置于–20℃冰箱保存。利用刺参16S rDNA、COⅠ和D-loop序列特异引物对3个mtDNA片段进行PCR扩增，对PCR产物测序，3对引物序列见表 2(Kerr et al, 2005; Arndt et al, 1996;姬南京等, 2014)，引物由青岛擎科生物公司合成。PCR总体积为50 μl：2×Tsingke Master mix 25 μl，10 μmol/L正反向引物各1 μl，模板DNA 1 μl (50 ng/μl)，加ddH₂O补足至50 μl。PCR反应扩增程序：95℃预变性60 s、95℃ 30 s、52℃ 30 s，72℃ 60 s，经35个循环后，72℃延伸10 min，4℃保存。PCR产物经1.5%的琼脂糖凝胶电泳检测后，送青岛擎科生物公司进行双向测序。

表 2(Tab.2)

表 2 刺参线粒体DNA扩增用引物序列 Tab.2 Primer sequences of mtDNA of A. japonicus 引物 Primer 序列 Sequences(5'–3') 引物出处 Reference 16S ar CGCCTGTTTATCAAAAACAT Kerr(2005) 16S br CTCCGGTTTGAACTCAGATCA COⅠ ef ATAATGATAGGAGGRTTTGG Arndt等(1996) COⅠ er GCTCGTGTRTCTACRTCCAT D-loop-F CAAGCCCTAAACCCAAGT 姬南京等(2014) D-loop-R ATTCCCTAGATGTAAGTTCG

表 2 刺参线粒体DNA扩增用引物序列 Tab.2 Primer sequences of mtDNA of A. japonicus

测序所得序列采用ClustalX1.83软件编辑、校对和排序(Thompson et al, 1997)。采用DnaSP 5.0软件计算多态位点数目(Number of polymorphicsites)、单倍型数目(Number of haplotypes)、单倍型多样性(Haplotype diversity)、核苷酸多样性(Nucleotide diversity)、平均核苷酸差异数(Average number of nucleotide differences)、两两群体间的遗传分化系数(F_st)(Librado et al, 2009)。利用MEGA 6.0软件统计碱基组成、变异位点数，计算群体内及群体间遗传距离，采用Kimura 2- Parameter方法构建N-J(Neighbor-joining)系统发育树(Saitou et al, 1987)。利用Network 5.0软件构建基于Median-Joining的单倍型网络图。

经PCR扩增、测序，获得刺参线粒体16S rDNA、COⅠ和D-loop序列，利用ClustalX1.83软件进行序列比对、校正去除序列两端非保守碱基后获得80个个体的相应序列信息。16S rDNA序列长度为543 bp，碱基组成比例相对稳定(表 3)，A的平均含量为27.5%，T的平均含量为29.1%，G的平均含量为21.7%，C的平均含量为21.7%，G与C的平均含量及各群体含量均相等，A+T的含量为56.6%，略高于G+C含量。COⅠ的序列长度为656 bp，A、T、G、C的平均含量分别为32.4%、26.5%、21.7%、19.4%，A+T的含量为58.9%，高于G+C含量。D-loop的序列长度为509~527 bp，A、T、G、C的平均含量分别为34.5%、25.3%、27.0%、13.3%，A+T的含量为59.7%，显著高于G+C含量。A+T含量在不同基因片段中差异显著，3个片段中D-loop序列最高，按照不同地理群体进行相应数据统计发现，在同一序列不同群体中差异不显著。

表 3(Tab.3)

表 3 刺参不同地理群体16S rDNA、COⅠ和D-loop序列的碱基组成(%) Tab.3 Base compositions of 16S rDNA, COⅠ and D-loop sequences in different populations of A. japonicus(%) 群体 Population 16S rDNA COⅠ D-loop A T G C A+T A T G C A+T A T G C A+T 青岛太平角QD 27.5 29.1 21.7 21.7 56.6 32.3 26.5 21.8 19.4 58.8 34.5 25.4 26.5 13.5 60.0 烟台长岛YT 27.5 29.1 21.7 21.7 56.6 32.3 26.5 21.8 19.4 58.8 34.6 25.4 26.9 13.1 60.0 韩国木浦黑参SK-MP-B 27.5 29.1 21.7 21.7 56.6 32.4 26.6 21.6 19.3 59.0 34.6 25.2 26.8 13.4 59.8 韩国群山黑参SK-QS-B 27.4 29.1 21.7 21.7 56.5 32.4 26.5 21.7 19.4 58.9 34.6 25.2 26.9 13.3 59.8 韩国浦项黄参SK-PX-Y 27.5 29.1 21.7 21.7 56.5 32.4 26.6 21.6 19.3 59.0 34.5 25.3 27.0 13.2 59.8 韩国浦项红参SK-PX-R 27.5 29.1 21.7 21.7 56.6 32.2 26.5 21.9 19.4 58.7 34.1 25.1 27.3 13.5 59.2 韩国浦项黑参SK-PX-B 27.5 29.1 21.7 21.7 56.6 32.3 26.6 21.8 19.3 58.9 34.3 25.1 27.2 13.3 59.4 俄罗斯海参RU 27.4 29.1 21.7 21.7 56.5 32.5 26.5 21.7 19.4 59.0 34.5 25.5 27.0 13.0 59.9 平均值Average 27.5 29.1 21.7 21.7 56.6* 32.4 26.5 21.7 19.4 58.9* 34.5 25.3 27.0 13.3 59.7* *表示差异较显著(P < 0.01) * represents significant difference (P < 0.01)

表 3 刺参不同地理群体16S rDNA、COⅠ和D-loop序列的碱基组成(%) Tab.3 Base compositions of 16S rDNA, COⅠ and D-loop sequences in different populations of A. japonicus(%)

16S rDNA的543个位点中，变异位点有16个(图 2)，占总位点数的2.9%，其中没有发现插入/缺失现象，转换多于颠换。在COⅠ基因的656个位点中，变异位点62个(图 3)，占总位点数的9.5%，其中，没有发现插入/缺失位点，转换多于颠换。D-loop的总长度为551 bp (包含缺失位点)，变异位点200个(图 4)，占总位点数的36.3%，其中，68个位点存在碱基的插入/缺失位点。

图 1(Fig. 1)

图 2 刺参16S rDNA序列变异位点 Fig.2 Variable site of 16S rDNA in A. japonicus 上方数字表示变异位点在单倍型中对应的位置，左侧表示单倍型名称，右侧表示个体数目 Numbers on the top represent the locations of mutative nucleotide acids, on the left are the names of haplotypes, on the right side represent the number of individuals

图 1(Fig. 1)

图 3 刺参COⅠ基因序列变异位点 Fig.3 Variable site of COⅠ gene in A. japonicus 上方数字表示变异位点在单倍型中对应的位置，左侧表示单倍型名称，右侧表示个体数目，N代表序列简并 Numbers on the top represent the locations of mutative nucleotide acids, on the left are the names of haplotypes, on the right side represent the number of individuals, N represents sequence degeneracy

图 1(Fig. 1)

图 4 刺参D-loop序列变异位点 Fig.4 Variable site of D-loop sequence in A. japonicus 上方数字表示变异位点在单倍型中对应的位置，左侧表示单倍型名称，右侧表示个体数目，N代表序列简并 Numbers on the top represent the locations of mutative nucleotide acids, on the left are the names of haplotypes, on the right side represent the number of individuals, N represents sequence degeneracy

8个群体间的16S rDNA、COⅠ、D-loop的核苷酸遗传多样性参数见表 4。从8个群体80个个体的16S rDNA序列中共检测到16个多态位点、16种单倍型，单倍型多样性指数为0.629，核苷酸多样性指数为0.0016，平均核苷酸差异数0.880，整体遗传多样性较低。8个群体中80个个体的COⅠ序列共检测到62个多态位点、38种单倍型，单倍型多样性指数为0.958，核苷酸多样性指数为0.0073，平均核苷酸差异数为4.796，说明8个群体的COⅠ序列单倍型及多态位点比较丰富，核苷酸差异数较大，其中，韩国浦项地区3个不同体色的刺参单倍型多样性比其他群体更高。从8个群体中80个个体的D-loop序列共检测到200个多态位点、61种单倍型，单倍型多样性指数为0.933，说明8个群体的D-loop序列的多态位点数及单倍型数最为丰富，单倍型多样性最高，核苷酸多样性指数为0.0157，平均核苷酸差异数为6.834，核苷酸差异数最大。对比3个片段的遗传多样性结果可以看出，COⅠ基因和D-loop的单倍型多样性、多态位点数及核苷酸差异数都显著高于16S rDNA。鉴于NCBI数据库中收录的海参纲COⅠ基因数据量高，而D-loop序列较少，为了保证后期数据分析的准确性和可比性，本研究选取COⅠ基因用于后续不同地理群体刺参遗传结构分析。

表 4(Tab.4)

表 4 刺参不同群体16S rDNA、COⅠ和D-loop片段的遗传多样性指数 Tab.4 The genetic diversity indices of partial 16S rDNA, COⅠ and D-loop sequence in different population of A. japonicus 序列片段 Gene fragment 群体 Population 样本数 Sample number 单倍型数H Number of haplotype 单倍型 多样性Hd Haplotype diversity 多态位 点数S Number of polymorphic 平均核苷酸差异数k Average number of nucleotide difference 核苷酸 多样性指数Pi Nucleotide diversity 16S rDNA 青岛太平角QD 10 3 0.511 2 0.556 0.0010 烟台长岛YT 10 5 0.756 5 1.156 0.0021 韩国木浦黑参SK-MP-B 10 3 0.600 2 0.667 0.0012 韩国群山黑参SK-QS-B 10 4 0.533 4 0.800 0.0015 韩国浦项黄参SK-PX-Y 10 4 0.644 4 1.067 0.0020 韩国浦项红参SK-PX-R 10 4 0.533 4 0.956 0.0018 韩国浦项黑参SK-PX-B 10 4 0.644 3 0.756 0.0014 俄罗斯海参RU 10 3 0.733 2 1.000 0.0018 总体Overall 80 16 0.629 16 0.880 0.0016 COⅠ 青岛太平角QD 10 8 0.933 16 4.600 0.0070 烟台长岛YT 10 5 0.867 10 3.733 0.0057 韩国木浦黑参SK-MP-B 10 9 0.978 16 4.156 0.0063 韩国群山黑参SK-QS-B 10 8 0.956 13 3.111 0.0047 韩国浦项黄参SK-PX-Y 10 8 0.956 22 5.378 0.0082 韩国浦项红参SK-PX-R 10 8 0.956 22 5.867 0.0089 韩国浦项黑参SK-PX-B 10 8 0.956 24 6.600 0.0101 俄罗斯海参RU 10 5 0.867 10 4.067 0.0062 总体Overall 80 38 0.958 62 4.796 0.0073 D-loop 青岛太平角QD 10 9 0.978 59 17.311 0.0345 烟台长岛YT 10 7 0.933 20 5.156 0.0109 韩国木浦黑参SK-MP-B 10 9 0.978 30 7.533 0.0159 韩国群山黑参SK-QS-B 10 7 0.933 47 14.756 0.0290 韩国浦项黄参SK-PX-Y 10 10 1.000 75 22.067 0.0441 韩国浦项红参SK-PX-R 10 10 1.000 58 16.200 0.0323 韩国浦项黑参SK-PX-B 10 9 0.978 49 13.778 0.0272 俄罗斯海参RU 10 7 0.933 37 16.511 0.0324 总体Overall 80 61 0.922 200 6.834 0.0157

表 4 刺参不同群体16S rDNA、COⅠ和D-loop片段的遗传多样性指数 Tab.4 The genetic diversity indices of partial 16S rDNA, COⅠ and D-loop sequence in different population of A. japonicus

根据所测定的8个群体的COⅠ序列信息，采用MEGA 6.0软件计算8个群体内和群体间的遗传距离的结果见表 5。从表 5可以看出，群体内遗传距离为0.0023~0.0084，群体内遗传距离最大的是韩国浦项红参群体(0.0084)，其次是韩国浦项黑参群体(0.0075)，俄罗斯符拉迪沃斯托克群体的群体内遗传距离最小(0.0023)。8个群体的群体间遗传距离为0.0027~ 0.0081，韩国浦项红参群体与韩国浦项黑参群体间的遗传距离最大(0.0081)，青岛海参群体、烟台群体与韩国浦项海域的3个海参群体遗传距离较大(0.0054~ 0.0071)，俄罗斯符拉迪沃斯托克群体与韩国群山群体的遗传距离最小(0.0027)。

表 5(Tab.5)

表 5 基于COⅠ序列的8个刺参群体内(对角线)和群体间的遗传距离(对角线下方)及遗传分化系数(对角线上方) Tab.5 Genetic distances within (diagonal) and between eight different populations (below diagonal), the fixation index (above diagonal) of A. japonicus 群体 Population COⅠ 青岛 QD 烟台 YT 韩国木浦黑参 SK-MP-B 韩国群山黑参 SK-QS-B 韩国浦项黄参 SK-PX-Y 韩国浦项红参 SK-PX-R 韩国浦项黑参 SK-PX-B 俄罗斯 RU 青岛QD 0.0054 -0.0532 0.0424 0.0454 0.0612 0.0362 0.0606 0.0931 烟台YT 0.0045 0.0040 0.1087 0.1171 0.1206 0.0826 0.1063 0.1598 韩国木浦黑参SK-MP-B 0.0047 0.0043 0.0036 -0.0168 -0.0192 0.0354 0.0225 0.0249 韩国群山黑参SK-QS-B 0.0044 0.0040 0.0033 0.0030 0.0016 0.0472 0.0404 0.0278 韩国浦项黄参SK-PX-Y 0.0058 0.0054 0.0045 0.0043 0.0055 0.0137 0.0053 -0.0125 韩国浦项红参SK-PX-R 0.0071 0.0068 0.0062 0.0060 0.0070 0.0084 0.0192 0.0246 韩国浦项黑参SK-PX-B 0.0069 0.0065 0.0057 0.0055 0.0065 0.0081 0.0075 0.0269 俄罗斯RU 0.0042 0.0037 0.0030 0.0027 0.0038 0.0054 0.0050 0.0023

表 5 基于COⅠ序列的8个刺参群体内(对角线)和群体间的遗传距离(对角线下方)及遗传分化系数(对角线上方) Tab.5 Genetic distances within (diagonal) and between eight different populations (below diagonal), the fixation index (above diagonal) of A. japonicus

通过DnaSP软件基于COⅠ序列信息分析刺参群体间的遗传分化的结果见表 5。遗传分化系数(F_st)反映了群体间的遗传分化程度，当F_st值在0~0.05时，群体间遗传差异很小，不存在分化；在0.05~0.15时，存在中等程度的遗传分化；0.15~0.25时，遗传分化程度较大；大于0.25时，则遗传分化很大(Wright, 1931)。俄罗斯群体与烟台群体间的遗传分化系数最大(0.1598)，遗传分化程度较高，烟台群体与韩国浦项、群山、木浦群体间均存在中等的遗传分化，青岛群体与韩国浦项黄参群体、韩国浦项黑参群体间也存在着中等遗传分化，而其他群体间差异较小，不存在分化。韩国浦项地区3个不同颜色的群体虽然遗传多样性较大，但遗传分化系数却较小，没有形成显著的遗传分化。

基于COⅠ基因构建的进化树(图 5)显示青岛海参群体与烟台海参群体聚为一支，再与韩国群山黑参群体聚为一支，与韩国木浦黑参群体聚在一起，再与俄罗斯群体聚为一支，向外依次为韩国浦项的3个不同体色的群体，这与其地理分布具有一定的相关性。韩国浦项由于其群体内遗传分化较大，在进化树中处于最外侧，可能与受其洋流影响有关。

图 1(Fig. 1)

图 5 8个群体刺参COⅠ基因的NJ系统发育树 Fig.5 Neighbor-joining tree of eight A. japonicus groups based on COⅠ gene

利用Network软件对所获得的COⅠ的38个单倍型构建单倍型网络图发现(图 6)，37个单倍型构成以Hap3为中心呈放射状分布的网络状，此外，Hap4、Hap6、Hap9、Hap12、Hap13、Hap18、Hap27、Hap32和Hap38成为外部辐射状网络节点，而其余单倍型则成为辐射图外部端点独立存在。Hap3存在于除俄罗斯群体和韩国浦项黄参群体外的群体中，为6个群体所共有，在37个单倍型中有10个单倍型为共享单倍型，占总单倍型的26.3%，28个单倍型为某个群体所特有。从网络图中难以分析出不同单倍型与地理位置的关系，因为其聚类无显著地理差异性，不同单倍型在不同群体间交叉分布，不同群体的单倍型具有共享性又具有特异性。

图 1(Fig. 1)

图 6 刺参COⅠ基因单倍型网络结构 Fig.6 The haplotype network of A. japonicus based on COⅠ gene

遗传多样性是指地球上所有生物所携带的遗传信息的总和，是评价群体资源状况的一个重要依据。mtDNA作为核外转录因子，具有演化速度快、母系遗传等优点，常被用于种群遗传结构分析(Liu et al, 2007)。其中，常用到的mtDNA标记包括16S rDNA、Cytb、COⅠ、D-loop。这些标记曾应用于不同种类鱼(郝君等, 2013;李大命等, 2017;孙超等, 2018)、甲鱼(Onychostoma rara)(彭珊等, 2009)、海胆属(Temnopleurus)(曾晓起等, 2012;刘晓慧等, 2007)、水母(Aurelia sp.)(张珰妮等, 2015;王建艳等, 2013)等的遗传结构分析。刺参的线粒体全长包含37个基因和3个非编码区(李云峰等, 2012)。前期多位研究者分别用16S rDNA、Cytb、COⅠ、D-loop等mtDNA片段对不同地理群体的遗传多样性进行了分析。如姬南京等(2014)利用D-loop序列的多态性对不同地理群体仿刺参的遗传结构进行了分析；李颖等(2006)利用16S rDNA、COⅠ、IrRNA-COⅠ基因对中国烟台、威海、莱州刺参进行遗传多样性分析及遗传结构分析。为评价不同片段分析刺参遗传多样性的效果，本研究应用16S rDNA、COⅠ、D-loop 3个mtDNA片段对6个地理种群8个群体进行了分析。3个片段对群体的遗传多样性检测结果表明，基于3个序列的各群体的核苷酸多样性指数为0.0007~0.0441，16S rDNA序列最为保守，表现出较低的遗传多样性，而COⅠ和D-loop遗传多样性较高。这与陈丽梅等(2008)研究的海参16S rDNA遗传多样性低于COⅠ基因的遗传多样性一致。相关结果均表明，16S rDNA不适于用于同一物种的种内不同群体的遗传结构分析，COⅠ和D-loop更适于进行种内群体间遗传多样性及遗传结构分析。在海参纲中，相比于D-loop序列，COⅠ基因信息收录于NCBI数据库中的种类及数目较多，更适合于进行比较分析，因而，本研究后续分析选用COⅠ基因。

本研究检测16S rDNA、D-loop和COⅠ序列的核苷酸多样性指数分别为0.0016、0.0157和0.0073，多样性指数较低。单倍型分别是16种、61种和38种，单倍型多样性分别为0.629、0.922和0.958，单倍型多样性较高，这种核苷酸多样性指数较低，单倍型多样性较高的现象，在其他物种海参线粒体基因片段中也广泛存在(Soliman et al, 2016; Uthicke et al, 2010)。对所获得的8个群体的遗传多样性分析结果表明，韩国浦项3个不同体色的野生刺参群体遗传多样性指数最高，相比之下，中国青岛和烟台野生群体刺参的遗传多样性较低，这与王鹤等(2016)的研究结果相一致，即韩国刺参群体的遗传多样性比中国群体的遗传多样性丰富。这种状况出现的原因可能是多方面的，首先地理位置差异导致其聚居区域洋流差异及海水温度差异较大，另一方面可能由于在中国青岛和烟台沿海刺参养殖较多，在排水过程中存在苗种逃逸现象，养殖刺参与野生刺参存在一定程度上的基因交流，近亲繁殖现象严重，从而导致青岛和烟台野生刺参出现遗传多样性较低的现象。

利用COⅠ基因对所采集的8个群体进行遗传距离分析并构建系统进化树，基于COⅠ基因序列的遗传距离为0.0023~0.0084。系统进化分析显示，不同群体的遗传结构与地理位置具有一定的相关性，即青岛、烟台、韩国木浦群体和韩国群山的刺参群体采样点为环渤海、黄海、东海相连通的海域地区，其中，青岛、烟台距离最近，在进化树中，这2个地理群体的刺参聚为一支，然后再与韩国群山、木浦群体聚为一支。韩国浦项群体位于韩国东海，俄罗斯群体采样点位于太平洋沿岸穆拉维约夫–阿穆尔斯基半岛的南端，临近日本海，地理距离与其他群体相对较远，这两地群体位于进化树外侧。本研究结果与潘传燕等(2012)利用微卫星标记进行的相关分析一致，即种群间相似程度与其地理分布具有一定的相关性。

本研究对采自于韩国浦项的3种不同颜色的刺参群体遗传多样性结果分析表明，该海域的刺参具有较高的遗传多样性，这可能与其所处的位置受日本暖流及千岛寒流影响，海水在次表层和中层2个层面中不断流动，运输着来自不同海域的营养物质，使海域中的浮游生物丰富度增加(王燕, 2015)，致使此地刺参群体变异比较大，具体的作用机制有待于进一步研究。战爱斌(2007)研究认为，洋流会促使幼虫甚至成体发生交换，会促使不同地域的海洋生物的生殖细胞、幼虫甚至成体发生交换，从而影响不同地理群体的遗传分化(Pante et al, 1992)。韩国浦项3个体色群体的遗传分化较小，还未达到分化成不同种的水平。Zhang等(2016)利用线粒体全基因组分析得出不同颜色的刺参属于同一物种，体色差异是由于海洋气候及洋流变化造成的。Sun等(2010)通过对红色、绿色、黑色刺参的线粒体基因组比较发现，不同体色刺参分化亦未达到不同物种的分化水平，这与本研究的结果是一致的。

本研究基于刺参3种线粒体DNA序列对中、韩、俄沿海野生刺参进行遗传多样性和遗传结构分析结果表明，COⅠ基因更适于进行刺参不同地理群体遗传结构分析；不同地理群体的进化关系跟其地理位置和所处海域的洋流关系密切；处于同一海域不同体色的刺参虽然其体色差异显著，但遗传分化不显著。