海洋科学  2020, Vol. 44 Issue (12): 23-31   PDF    
http://dx.doi.org/10.11759/hykx20200328001

文章信息

钟凯乐, 宋小含, 段德麟, 胡自民. 2020.
ZHONG Kai-le, SONG Xiao-han, DUAN De-lin, HU Zi-min. 2020.
中国黄渤海真江蓠(Agarophyton vermiculophyllum)群体遗传多样性研究
Population genetic diversity of the red alga Agarophyton vermiculophyllum (Gigartinales, Rhodophyta) in the Yellow-Bohai Sea
海洋科学, 44(12): 23-31
Marine Sciences, 44(12): 23-31.
http://dx.doi.org/10.11759/hykx20200328001

文章历史

收稿日期：2020-03-28
修回日期：2020-04-20
引用本文
钟凯乐, 宋小含, 段德麟, 胡自民. 2020. 中国黄渤海真江蓠(Agarophyton vermiculophyllum)群体遗传多样性研究[J]. 海洋科学, 44(12): 23-31.
ZHONG Kai-le, SONG Xiao-han, DUAN De-lin, HU Zi-min. 2020. Population genetic diversity of the red alga Agarophyton vermiculophyllum (Gigartinales, Rhodophyta) in the Yellow-Bohai Sea[J]. Marine Sciences, 44(12): 23-31.
中国黄渤海真江蓠(Agarophyton vermiculophyllum)群体遗传多样性研究
钟凯乐1,2,3, 宋小含1,2,3, 段德麟1,2, 胡自民1,2     
1. 中国科学院海洋研究所海洋大科学研究中心 实验海洋生物学重点实验室, 山东 青岛 266071;
2. 青岛海洋科学与技术试点国家实验室 海洋生物学与生物技术功能实验室, 山东 青岛 266237;
3. 中国科学院大学, 北京 100049
收稿日期：2020-03-28；修回日期：2020-04-20
基金项目：国家自然科学基金项目（41761144057，31971395）
作者简介：钟凯乐(1995-), 女, 汉族, 山东德州人, 硕士, 研究方向为海藻系统进化与多样性, E-mail:18865550027@163.com.
通信作者：胡自民, 通信作者, 男, 博士, 研究员, 研究方向为海藻系统进化与多样性, E-mail:huzm@qdio.ac.cn.
摘要：真江蓠是原产于西北太平洋的重要经济红藻。我们利用10对微卫星引物检测中国黄渤海地区真江蓠的群体遗传多样性和结构。10个微卫星位点在12个群体中共检测到65个等位基因，每个位点的等位基因（Na）为1~28，有效等位基因（Ne）为1.0~9.6。每个群体的平均等位基因（Na）、平均有效等位基因（Ne）、平均香浓指数（I）、平均观察杂合度（Ho）和平均预期杂合度（He）分别为2.4、1.6、0.419、0.133和0.227，显示较低的群体遗传多样性。中国黄渤海12个真江蓠群体间遗传分化较大（Fst=0.398 7），基因流有限（Nm=0.377 1），近交系数为正（Fis=0.391 3，Fit=0.634 0），表明可能存在近交和杂合子缺失现象。Structure和UPGMA系统进化分析一致将12个群体分为两个遗传组，并在黑石礁群体（HS）和石岛群体（SD）中发现明显的遗传混杂现象。AMOVA分析显示遗传变异主要来自于群体内（73.27%）。该研究可为黄渤海地区真江蓠自然资源保护和管理提供科学依据。
关键词：真江蓠(Agarophyton vermiculophyllum)    微卫星    遗传多样性    基因流    黄渤海    
Population genetic diversity of the red alga Agarophyton vermiculophyllum (Gigartinales, Rhodophyta) in the Yellow-Bohai Sea
ZHONG Kai-le1,2,3, SONG Xiao-han1,2,3, DUAN De-lin1,2, HU Zi-min1,2     
1. Key Laboratory of Experimental Marine Biology, Center for Ocean Mega-Science, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China;
2. Laboratory for Marine Biology and Biotechnology, Pilot National Laboratory for Marine Science and Technology(Qingdao), Qingdao 266237, China;
3. University of Chinese Academy of Sciences, Beijing 100049, China
Received: Mar. 28, 2020
Foundation: the National Natural Science Foundation of China, No. 41761144057, No. 31971395
Abstract: Agarophyton vermiculophyllum is one of the most economically important red algae native to the northwest Pacific. We surveyed the population genetic diversity and the structure of A. vermiculophyllum in the Yellow-Bohai Sea using 10 microsatellites. These microsatellites detected 65 alleles among 12 populations. The number of alleles (Na) varied from 1 to 28, and the number of effective alleles (Ne) varied from 1.0 to 9.6. The average number of alleles (Na), the average number of effective alleles (Ne), the average Shannon's diversity index (I), the average observed heterozygosity (Ho), and the average expected heterozygosity (He) were 2.4, 1.6, 0.419, 0.133, and 0.227, respectively, indicating a low population genetic diversity. There was a high genetic differentiation among the 12 A. vermiculophyllum populations (Fst=0.398 7) and a limited gene flow (Nm=0.377 1). The average inbreeding coefficient was positive (Fis=0.391 3, Fit=0.634 0), which indicates inbreeding and heterozygotic deletion. Structure analysis and UPGMA tree analysis consistently divided the 12 populations into two genetic groups, and obvious genetic mixing was found in the Heishijiao, Dalian and Shi Island, Weihai populations. Analysis of molecular variation showed that the genetic variation in the Yellow-Bohai Sea mainly occurred at the intra-population level (73.27%). This study may provide a scientific basis for the conservation of natural resources and management of A. vermiculophyllum in the Yellow-Bohai Sea.
Key words: Agarophyton vermiculophyllum    microsatellite    genetic diversity    gene flow    Yellow-Bohai Sea    

气候变化和近岸环境变迁是影响潮间带海洋生物群体遗传多样性和分化的重要因素。黄渤海作为太平洋西部的边缘海, 形成于第四纪更新世晚期[1], 在末次盛冰期(Last Glacial Maximum, LGM)经历了剧烈的气候和海岸环境变迁。更新世冰期导致全球海平面下降120~140 m[2], 使得西北太平洋边缘海面积缩小, 黄渤海陆架几乎完全暴露成为陆地, 东海向东缩减变为狭长的冲绳海槽[2-3]。边缘海面积的缩减导致海洋物种大量减少, 残遗物种只能在避难所生存, 因此避难所的不连续性成为不同群体之间隔离和分化的重要因素[4], 其中冲绳海槽被证实是多种海洋生物的冰期避难所[3, 5]。在更新世间冰期, 气温回暖海平面上升, 残遗物种从避难所扩展到新的栖息地。这种群体扩张可能导致来自不同避难所的群体之间的遗传交换, 消除了隔离造成的群体分化[6]。除此之外, 黄渤海的洋流系统对海洋生物的遗传多样性和分化也具有重要影响[7]。主导黄渤海海域的洋流主要是北向的黑潮分支黄海暖流和南向的中国沿岸流, 黄海暖流和中国沿岸流驱动的群体遗传均质性在鼠尾藻(Sargassum thunbergii)[7-8]和角叉菜(Chondrus ocellatus)[5]等多种海洋生物中已有报道。

真江蓠(Agarophyton vermiculophyllum)隶属于红藻门(Rhodophyta), 真红藻纲(Florideophyceae), 江蓠目(Gracilariales), 江蓠科(Gracilariaceae), Agarophyton属, 是原产于亚洲-西北太平洋的一种食用红藻, 也是一种在世界范围内广泛分布的入侵种[9]。真江蓠不仅是提取琼胶的主要原料, 而且是重要的生态奠基物种, 因此具有重要的经济和生态价值[10-11]。除此之外, 真江蓠还可作为生物过滤器净化水质, 在海水富营养化修复和多营养层次综合水产养殖系(Integrated Multi-Trophic Aquaculture, IMTA)中都具有广泛应用[12-13]。

目前, 相对于真江蓠的入侵, 对该物种的群体遗传学研究仍然较少, 且大多基于单一的线粒体分子标记。Krueger-Hadfield等[9]基于线粒体cox1 (~1 200 bp)检测到黄渤海地区真江蓠群体为一个单一遗传谱系, 但该研究主要涉及朝鲜半岛西海岸群体, 中国沿海采样点极少。刘若愚等[14]基于cox1(641 bp)在中国黄渤海沿岸真江蓠群体中检测到10个单倍型, 为一个共同的遗传谱系, 群体间遗传分化较弱[14]。Hu等[15]基于AFLP(Amplified Fragment Length Polymorphism)发现黄渤海地区真江蓠群体为共同的遗传谱系, 并且在青岛地区群体中发现小尺度的遗传分化。

微卫星具有遗传信息丰富、多态性高、中性选择、共显性遗传等特点, 重复性好、易于检测、稳定性高, 在海藻群体遗传多样性研究中具有广泛应用[16-17]。本研究共收集中国黄渤海12个群体的196株真江蓠样品, 利用微卫星数据旨在深入分析该地区真江蓠的群体多样性特征, 同时检测是否在小地理尺度存在遗传结构和遗传变异的时空分布模式。该研究有助于理解气候变化与当代环境的相互作用对海藻遗传变异的影响, 从而为真江蓠自然资源保护和管理提供科学依据。

1 材料与方法 1.1 真江蓠样本采集

我们采集到来自辽东半岛和山东半岛潮间带12个地点的真江蓠样本(表 1)。在每个采样点, 随机采集6~20个真江蓠个体, 确保个体之间距离大于10 m以最大程度避免采集到遗传一致的个体。样品清洗干净后置于硅胶中干燥保存。

表 1 黄渤海真江蓠样本采集信息 Tab. 1 Sampling information regarding A. vermiculophyllum in the Yellow-Bohai sea
简称 地理位置 经纬度 样本量
ZZ 辽宁大连獐子岛 39.04°N, 122.74°E 20
HS 辽宁大连黑石礁 38.87°N, 121.56°E 8
HN 辽宁大连黄泥川 38.82°N, 121.56°E 20
CD 山东烟台长岛 37.94°N, 120.72°E 16
JM 山东威海鸡鸣岛 37.45°N, 122.48°E 6
DC 山东威海东楮岛 37.04°N, 122.56°E 20
SD 山东威海石岛 36.91°N, 122.41°E 13
LR 山东青岛石老人 36.09°N, 120.49°E 20
YH 山东青岛银海国际 36.06°N, 120.42°E 20
YY 山东青岛一浴 36.05°N, 120.34°E 16
EY 山东青岛二浴 36.05°N, 120.34°E 17
SY 山东青岛三浴 36.05°N, 120.36°E 20
表选项
1.2 DNA提取与检测

取30 mg真江蓠样本, 在液氮中冷冻并研磨成粉末, 使用诺贝莱一步法多糖多酚植物基因组DNA抽提试剂盒DNE379(Nobelab)提取DNA, 具体操作见说明书。1%琼脂糖凝胶电泳(140 V, 30 min)检测DNA质量。

1.3 PCR扩增和电泳检测

我们使用Kollars等[18]和Krueger-Hadfield等[19]开发的10对真江蓠微卫星引物进行片段扩增(表 2), 引物由北京擎科生物科技有限公司青岛分公司合成。本实验我们采用双重荧光检测, F正向引物的5′端连接荧光基团(FAM或HEX)。PCR扩增反应体系15 μL:模板DNA 1.0 μL(约15 ng), 2*Tsingke Master mix 7.5 μL, 正向引物和反向引物各1.0 μL(10 μmol/L), 灭菌水4.5 μL。PCR反应程序: 98℃预变性2 min; 98℃变性10 s, 56℃退火10 s, 72℃延伸10 s, 30个循环; 最终72℃延伸5 min。PCR扩增产物在3730型遗传分析仪(ABI, USA)上进行毛细管电泳检测(北京擎科生物科技有限公司青岛分公司)。

表 2 真江蓠10个微卫星位点的信息 Tab. 2 Information about A. vermiculophyllum obtained using 10 microsatellites
位点 序列号 核心序列 引物序列(5′ to 3′) TM 片段大小 荧光 Na Ne D
Gverm 5276 KT232089 (AC)10 F: GGAGAGCAGCACGTTTTAGG
R: CTGCTTAGTTCCACGATCGAC		 57.4℃
57.6℃		 282~316 FAM 28 9.614 –0.783
Gverm 6311 KT232090 (AG)9 F: GCGTCATTCCACTGAATGTG
R: GATGAACCTCAATGCCTCGT		 55.4℃
55.4℃		 203~223 FAM 10 4.563 –0.652
Gverm 8036 KT232091 (AC)12 F: GCCCTTTTAAGGATGCAACA
R: GGGGTAAACGACCACAGAGA		 53.4℃
57.4℃		 213~251 FAM 4 1.612 –0.772
Gverm 3003 KT232092 (AG)11 F: CATCTTGCTTCCTTGCTTCC
R: TTGAAAGCCGGAATTTATCG		 55.4℃
51.3℃		 198~230 HEX 3 1.058 –0.721
Gverm 1203 KT232093 (AAG)8 F: CTCCTGGTGCACAAGCAATA
R: ACATTCTGCGCACCTTTCTT		 55.4℃
53.4℃		 284~308 HEX 4 2.106 –0.699
Gverm 1803 KT232094 (AC)11 F: GCGTGCACGATGTCTACACT
R: GACAGCAACAAGTGGGGTTT		 57.4℃
55.4℃		 352~356 HEX 3 1.142 –0.017
Gverm 804 KT232095 (AAG)8 F: TGTAGGATTGCTCTCCTGGTG
R: CAGGCTGGCCAAAATAACAT		 57.6℃
53.4℃		 182~188 FAM 4 2.189 –0.634
Gverm 10367 KT232096 (AG)8 F: GCTGAGAAATGAAGCGAAGG
R: GCAAACCTGCCTTGTTTGTT		 55.4℃
53.4℃		 198~200 FAM 4 1.132 –0.824
Gverm 2790 KT232097 (AATGC)5 F: GAACAATGCGGGAAAACATT
R: GGAAGAGGCTCAAAAGCAGA		 51.3℃
55.4℃		 262~267 HEX 4 1.542 0.190
Gverm 7969 KT232099 (ATC)6 F: CGATCCTTCCTCCTGTGGTA
R: AATTGGGATACGCAATACGG		 57.4℃
53.4℃		 约260 HEX 1 1.000 0.000
表选项
1.4 数据分析

利用软件Gene mapper v4.1读取数据并汇总。Popgene v1.32[20]和GenAlex v6.501[21]计算各微卫星位点和群体的遗传参数。包括等位基因数(Na), 有效等位基因数(Ne), 香浓指数(I), 观察杂合度(Ho), 预期杂合度(He), 哈迪-温伯格平衡偏离指数(D=Ho– He/He), F-统计量(Fis, Fit和Fst), 基因流(Nm), Nei’s遗传距离和相似性。

用Poptree v2[22]构建基于遗传距离Da的UPGMA聚类树, 设置自展值为(bootstrap)1 000。在Structure v2.3.4[23]中基于贝叶斯聚类法对群体遗传结构进行分析。参数设置为: Length of Burnin period为105, Number of MCMC Reps after Burnin为106, K值设置为1~5, 每个K值重复运算10次。基于Delta K=mean(|L″(K)|)/ sd(L(K)), 在Structure Harvester中确定最佳K值, 即最佳分类群数。

在Arlequin v3.5.1.3[24]中进行分子方差分析(AMOVA), 基于104次重复抽样检测群体间和群体内的遗传结构和遗传变异的分布。

2 结果 2.1 遗传参数

10个微卫星位点共检测出65个等位基因, 每个微卫星位点等位基因(Na)为1~28, 有效等位基因(Ne)为1.000~9.614。其中Gverm 5276和Gverm 7969分别具有最高和最低的等位基因和有效等位基因数目(表 2)。

每个群体在10个微卫星位点的遗传参数见表 3。各群体的等位基因数(Na)为1.4~3.5, 平均等位基因数为2.4;有效等位基因数(Ne)为1.123~2.346, 平均有效等位基因数为1.644;香浓指数(I)为0.120~ 0.649, 平均香浓指数为0.419;观测杂合度(Ho)和预期杂合度(He)分别为0.012~0.276和0.073~0.349, 平均观测杂合度和预期杂合度分别为0.133和0.227, 各群体的观测杂合度均低于预期杂合度。12个群体中哈迪-温伯格偏离指数(D)均为负值, 表现为杂合子缺失。

表 3 10个微卫星位点在12个真江蓠群体的遗传多样性参数 Tab. 3 Genetic diversity parameters calculated using 10 microsatellites in 12 A. vermiculophyllum populations
群体 参数
Na Ne I Ho He D
ZZ 1.9 1.239 0.258 0.072 0.145 –0.503
HS 1.8 1.184 0.242 0.038 0.130 –0.708
HN 2.1 1.557 0.478 0.276 0.300 –0.024
CD 1.4 1.123 0.120 0.013 0.073 –0.822
JM 1.4 1.296 0.216 0.000 0.139 –1.000
DC 3.2 2.346 0.649 0.155 0.330 –0.530
SD 2.0 1.493 0.382 0.108 0.223 –0.516
LR 3.1 1.755 0.538 0.160 0.276 –0.420
YH 2.5 1.776 0.399 0.140 0.205 –0.317
YY 2.9 2.168 0.547 0.181 0.282 –0.358
EY 3.0 1.852 0.625 0.276 0.349 –0.209
SY 3.5 1.936 0.580 0.180 0.274 –0.343
平均 2.4 1.644 0.419 0.133 0.227 —
表选项
2.2 群体遗传结构

12个真江蓠群体的F-统计量及基因流分析(表 4)显示, 近交系数Fis为–0.506 4~0.855 8, 平均近交系数Fis为0.391 3, 其中Gverm 1803和Gverm 5276的近交系数为负值; 各位点的总群体近交系数Fit为–0.144 8~0.770 4, 平均值为0.634 0, Gverm 5276为负值; 度量群体间遗传差异程度的Fst为0.000 0~ 0.730 6, 平均值为0.398 7, 遗传分化程度较高; 基因流Nm为0.092 2~2.214 4, 平均值为0.377 1, 不同位点的基因流值变化较大。

表 4 F-统计量和基因流 Tab. 4 F statistics and gene flow obtained from 10 microsatellites
位点 Fis Fit Fst Nm
Gverm 5276 0.730 9 0.802 6 0.266 4 0.688 6
Gverm 6311 0.415 6 0.655 9 0.411 2 0.357 9
Gverm 8036 0.325 9 0.818 4 0.730 6 0.092 2
Gverm 3003 0.766 9 0.854 6 0.376 4 0.414 2
Gverm 1203 0.441 4 0.738 6 0.532 1 0.219 8
Gverm 1803 –0.187 8 0.034 1 0.186 8 1.088 1
Gverm 804 0.380 8 0.662 7 0.455 3 0.299 1
Gverm 10367 0.855 8 0.870 4 0.101 4 2.214 4
Gverm 5276 –0.506 4 –0.144 8 0.240 1 0.791 4
Gverm 6311 — — 0.000 0 —
平均 0.391 3 0.634 0 0.398 7 0.377 1
表选项

12个真江蓠群体Nei's遗传相似性及遗传距离见表 5。各群体间遗传相似性指数以银海国际(YH)和三浴(SY)之间最大(0.987), 鸡鸣岛(JM)和银海国际(YH)之间最小(0.497);相应的, Nei's遗传距离以鸡鸣岛(JM)和银海国际(YH)之间最大(0.699), 银海国际(YH)和三浴(SY)之间最小(0.013)。

表 5 12个真江蓠群体间的Nei's 遗传相似性(上三角)和遗传距离(下三角) Tab. 5 Nei's genetic identity (above diagonal) and genetic distance (below diagonal) among 12 A. vermiculophyllum populations
ZZ HS HN CD JM DC SD LR YH YY EY SY
ZZ — 0.773 0.720 0.852 0.851 0.734 0.751 0.641 0.574 0.753 0.663 0.607
HS 0.257 — 0.812 0.735 0.692 0.889 0.777 0.793 0.741 0.887 0.833 0.774
HN 0.329 0.208 — 0.677 0.639 0.865 0.766 0.849 0.777 0.885 0.876 0.826
CD 0.160 0.307 0.390 — 0.839 0.727 0.648 0.624 0.549 0.734 0.653 0.581
JM 0.162 0.368 0.448 0.175 — 0.649 0.632 0.583 0.497 0.668 0.586 0.524
DC 0.309 0.117 0.145 0.319 0.433 — 0.898 0.956 0.932 0.980 0.960 0.948
SD 0.286 0.253 0.267 0.433 0.458 0.107 — 0.912 0.888 0.905 0.864 0.896
LR 0.445 0.232 0.163 0.472 0.539 0.045 0.092 — 0.964 0.949 0.960 0.979
YH 0.555 0.300 0.252 0.600 0.699 0.071 0.119 0.037 — 0.908 0.949 0.987
YY 0.283 0.120 0.122 0.310 0.403 0.020 0.100 0.052 0.096 — 0.948 0.930
EY 0.411 0.183 0.133 0.427 0.534 0.040 0.147 0.041 0.052 0.054 — 0.962
SY 0.498 0.256 0.191 0.542 0.647 0.053 0.110 0.021 0.013 0.073 0.038 —
表选项

12个真江蓠群体的遗传分化系数(Fst)和基因流(Nm)见表 6。群体间最大的遗传分化在长岛(CD)和银海国际(YH)之间(0.536), 最小在石老人(LR)和三浴(SY)之间(0.022)。最大基因流发生在三浴(SY)和石老人(LR)之间(11.166), 最小基因流发生在银海国际(YH)和长岛(CD)之间(0.216)。山东半岛南侧的群体之间遗传分化较小, 基因交流频繁。

表 6 12个真江蓠群体间的遗传分化系数Fst (下三角)和基因流Nm(上三角) Tab. 6 Fst (below diagonal) and Nm (above diagonal) among 12 A. vermiculophyllum populations
ZZ HS HN CD JM DC SD LR YH YY EY SY
ZZ — 0.579 0.658 0.905 0.752 0.758 0.640 0.580 0.275 0.683 0.554 0.375
HS 0.302 — 1.008 0.430 0.417 1.751 0.858 1.094 0.530 1.702 1.142 0.744
HN 0.275 0.199 — 0.503 0.498 2.144 1.107 1.745 0.918 2.069 2.447 1.360
CD 0.216 0.368 0.332 — 0.577 0.634 0.404 0.484 0.216 0.544 0.460 0.294
JM 0.249 0.375 0.334 0.302 — 0.563 0.348 0.430 0.229 0.492 0.442 0.300
DC 0.248 0.125 0.104 0.283 0.308 — 2.645 5.799 2.346 10.063 7.650 4.359
SD 0.281 0.226 0.184 0.382 0.418 0.086 — 3.064 1.643 2.737 2.000 2.471
LR 0.301 0.186 0.125 0.340 0.368 0.041 0.075 — 5.081 4.774 5.600 11.166
YH 0.476 0.320 0.214 0.536 0.522 0.096 0.132 0.047 — 1.517 3.042 10.275
YY 0.268 0.128 0.108 0.315 0.337 0.024 0.084 0.050 0.142 — 4.064 2.870
EY 0.311 0.180 0.093 0.352 0.361 0.032 0.111 0.043 0.076 0.058 — 5.493
SY 0.400 0.251 0.155 0.460 0.455 0.054 0.092 0.022 0.024 0.080 0.044 —
表选项

Structure分析显示黄渤海12个真江蓠群体最佳遗传聚类数K=2(图 1), 与UPGMA系统发育树结构一致(图 2)。獐子岛(ZZ)、长岛(CD)和鸡鸣岛(JM)为一个遗传组, 其他群体聚为第二个遗传组, 其中黑石礁(HS)和石岛(SD)具有明显的遗传混杂现象。AMOVA分析显示群体间遗传变异占36.73%, 群体内变异占73.27%, 表明群体内遗传变异是引起黄渤海真江蓠群体变异的主要因素(表 7)。

图 1 基于10个微卫星的12个真江蓠群体的遗传结构图(底图审图号: GS(2019)1711号) Fig. 1 Genetic structures of 12 A. vermiculophyllum populations based on 10 microsatellites
图选项

图 2 基于Nei's遗传距离构建的12真江蓠群体的UPGMA聚类树 Fig. 2 UPGMA tree of 12 A. vermiculophyllum populations based on Nei's genetic distance 注:分支上的数字表示检验重复1 000次所得的大于50%的支持率
图选项

表 7 12个真江蓠群体的分子方差分析 Tab. 7 AMOVA of 12 A. vermiculophyllum populations
变异来源 自由度 平方和 遗传变异组成 百分比
群体间 11 268.282 0.714 38 36.73%
群体内 380 467.700 1.230 79 63.27%
总体 391 735.982 1.945 17
表选项
3 讨论

等位基因、有效等位基因、杂合度等是判断群体遗传多样性的重要指标, 种群的遗传多样性越高其环境适应能力越强, 相应的进化潜力也越大。本研究中12个真江蓠群体的Na(2.4), Ne(1.644), I(0.419), Ho(0.133)和He(0.227)等遗传参数较低, 与刘若愚等[14]利用cox1(Hd=0.300, π=0.050×10–2)和Hu等[15]利用AFLP(Na, 0.336~1.224; Ne, 1.057~1.142; I, 0.059~0.161)检测的结果基本吻合, 一致表明中国黄渤海地区真江蓠群体遗传多样性和变异程度较低, 环境适应能力可能较差。Fis(0.391 3)、Fit(0.634 0)和哈迪-温伯格平衡偏离指数D(D < 0)表明该地区真江蓠群体普遍存在杂合子缺失现象。这些遗传特征可能与真江蓠的繁殖方式和繁殖特点有关。在自然界, 相对于四分孢子体植株, 江蓠科海藻雌雄配子体植株数量较少[25-26]。真江蓠等红藻进行有性繁殖时, 释放的四分孢子和精子存活时间较短, 没有鞭毛结构且不具自主游动能力, 只能随着水流在近距离的雌配子果胞上完成受精[27], 这种有限距离的扩散降低了群体间遗传交换, 导致群体遗传多样性偏低。此外, 真江蓠也可通过断枝方式进行营养繁殖[28], 而营养繁殖的方式也大大降低了真江蓠的杂合度和群体遗传多样性。

通过群体间遗传分化系数Fst、Nei's遗传距离和基因流Nm可知黄渤海地区真江蓠群体之间遗传分化较大, 基因流有限。Krueger-Hadfield等[9]利用微卫星在日本地区真江蓠群体间检测到类似的较高的群体间遗传分化模式, 并证实遗传分化系数Fst与地理距离呈正相关。Hu等[15]利用AFLP也发现黄渤海真江蓠群体间较高的遗传分化, 与本研究结果一致。除了真江蓠的繁殖特点外, 群体间栖息地的不连续性造成的地理隔离可能是另一个驱动因素[29]。真江蓠可适应多种栖息环境, 在河口和海岸潮间带环境均有分布, 生境的不连续性和异质性(例如:泥滩、沙滩、礁石和石砾)可能导致了群体水平的遗传分化。此外, 这种小地理尺度的遗传分化可能也与微地理环境差异有关[30], 例如盐度与潮汐。黄渤海不同季节和地区盐度变化较大[31], 而真江蓠可以适应5‰~ 60‰的盐度变化[32]。已有报道显示盐度梯度可对等位基因频率产生强选择(例如, Furcellaria lumbricalis)[33], 对无融合生殖和表型变异产生影响并最终导致基因型分化(例如, Ulva prolifera)[34]。山东半岛和辽东半岛毗邻黄海, 属于温带季风气候, 全年处于不规则半日潮中。高滨岸和低滨岸之间的潮汐强度也可能对真江蓠产生选择压力, 从而产生微地理尺度的遗传分化。

通过Structure遗传结构可知, 中国黄渤海12个真江蓠群体分化为2个遗传组, 但是并没有发现如Hu等[15]在青岛地区检测到的特有遗传分化现象。Hu等[15]检测到黄渤海真江蓠地理障碍发生在JM和ZZ群体间, 而本研究发现这种遗传隔离发生在山东半岛的东南部(JM和DC群体之间), 类似的现象在该地区并不常见。众多研究表明, 由于更新世冰期海平面下降, 现今黄渤海地区的海洋生物大多起源于冲绳海槽避难所, 并在间冰期随着黑潮、黄海暖流和中国沿岸流向黄渤海和东海地区扩张迁移, 故该地区的海洋生物大多属于同一遗传谱系, 具有较强的遗传均质性[5, 7-8]。因此, 我们推测这种遗传隔离可能与山东半岛形成的地理屏障和真江蓠有限传播能力共同造成的基因交流障碍有关。

刘若愚[35]检测到真江蓠的基因流方向为辽东半岛至山东半岛至青岛, 与该地区洋流运动相一致, 这可能是山东半岛南部真江蓠群体(SD)出现北部基因型(遗传组1)的重要原因。我们在辽东半岛的群体中(HS和HN)发现山东半岛南部群体的主要基因型(遗传组2), 推测这种逆洋流方向的遗传混杂可能与辽东半岛和山东半岛之间频繁的海上航运等人为因素有关。虽然真江蓠的孢子和精子运动能力有限, 但作为一个世界范围内的入侵种它可以附着在船体、缆绳、鱼竿、虾和蟹笼等硬基质上, 通过海上运输、休闲渔业和水产养殖活动等传播扩散[36]。

4 结论

本研究基于10对微卫星引物对12个黄渤海真江蓠群体进行群体遗传多样性分析, 发现该地区真江蓠群体遗传多样性较低, 杂合子缺失, 并在山东半岛东南部检测到罕见的遗传隔离。该研究对黄渤海海藻多样性形成和自然资源保护提供了科学基础。

参考文献
[1]
耿秀山. 黄渤海地貌特征及形成因素探讨[J]. 地理学报, 1981, 36(4): 423-434.
Geng Xiushan. The geomorphological features and forming factors of submarine relief in the Bohai Sea and the Yellow Sea[J]. Acta geographica sinica, 1981, 36(4): 423-434. DOI:10.3321/j.issn:0375-5444.1981.04.008
[2]
Lambeck K, Esat T M, Potter E K. Links between climate and sea levels for the past three million years[J]. Nature, 2002, 419(6903): 199-206. DOI:10.1038/nature01089
[3]
Wang P X. Response of Western Pacific marginal seas to glacial cycles:Paleoceanographic and sedimentological features[J]. Marine Geology, 1999, 156(1-4): 5-39. DOI:10.1016/S0025-3227(98)00172-8
[4]
Ni G, Li Q, Kong L F, et al. Comparative phylogeography in marginal seas of the northwestern Pacific[J]. Molecular Ecology, 2014, 23(3): 534-548. DOI:10.1111/mec.12620
[5]
Hu Z M, Li J J, Sun Z M, et al. Phylogeographic structure and deep lineage diversication of the red alga Chondrus ocellatus Holmes in the Northwest Pacific[J]. Molecular Ecology, 2015, 24(19): 5020-5033. DOI:10.1111/mec.13367
[6]
Xue D X, Wang H Y, Zhang T, et al. Population genetic structure and demographic history of Atrina pectinata based on mitochondrial DNA and microsatellite markers[J]. PLoS One, 2014, 9: e95436. DOI:10.1371/journal.pone.0095436
[7]
Li J J, Hu Z M, Gao X, et al. Oceanic currents drove population genetic connectivity of the brown alga Sargassum thunbergii in the north-west Pacific[J]. Journal of Biogeography, 2017, 44: 230-242. DOI:10.1111/jbi.12856
[8]
Liu J, Zhan D, Zhang R, et al. AFLP analysis revealed a north to south genetic break in the brown alga Sargassum thunbergii along the coast of China[J]. Journal of Applied Phycology, 2018, 30(4): 1-9.
[9]
Krueger-Hadfield S A, Kollars N M, Strand A E, et al. Genetic identification of source and likely vector of a widespread marine invader[J]. Ecology and Evolution, 2017, 7(12): 4432-4447. DOI:10.1002/ece3.3001
[10]
史升耀, 张燕霞, 范晓, 等. 碱处理对中国江蓠属海藻所含琼胶的作用[J]. 水产学报, 1988, 12(2): 145-155.
Shi Yaosheng, Zhang Yanxia, Fan Xiao, et al. The effects of alkali treatment on agars from chinese species of Gracilaria[J]. Journal of fisheries of China, 1988, 12(2): 145-155.
[11]
Davoult D, Surget G, Stiger-Pouvreau V, et al. Multiple effects of a Gracilaria vermiculophylla invasion on estuarine mudflat functioning and diversity[J]. Marine Environmental Research, 2017, 131: 227-235. DOI:10.1016/j.marenvres.2017.09.020
[12]
徐姗楠, 温珊珊, 吴望星, 等. 真江蓠(Gracilaria verrucosa)对网箱养殖海区的生态修复及生态养殖匹配模式[J]. 生态学报, 2008, 28(4): 1466-1475.
Xu Shannan, Wen Shanshan, Wu Wangxing, et al. Biremediation of caged fish aquaculture by the red alga Gracilaria verrucosa in an integrated multi-trophic aquaculture system[J]. Acta Ecologica Sinica, 2008, 28(4): 1466-1475. DOI:10.3321/j.issn:1000-0933.2008.04.013
[13]
Abreu M H, Pereira R, Yarish C, et al. IMTA with Gracilaria vermiculophylla:Productivity and nutrient removal performance of the seaweed in a land-based pilot scale system[J]. Aquaculture, 2011, 312(1-4): 77-87. DOI:10.1016/j.aquaculture.2010.12.036
[14]
刘若愚, 孙忠民, 姚建亭, 等. 中国近海重要生态建群红藻真江蓠的群体遗传多样性[J]. 生物多样性, 2016, 24(7): 781-790.
Liu Ruoyu, Sun Zhongmin, Yao Jianting, et al. Genetic diversity of the habitat-forming red alga Gracilaria vermiculophylla along Chinese coasts[J]. Biodiversity Science, 2016, 24(7): 781-790.
[15]
Hu Z M, Liu R Y, Zhang J, et al. A unique genetic lineage at the southern coast of China in the agar-producing Gracilaria vermiculophylla (Gracilariales, Florideophyceae)[J]. Algae, 2018, 33: 269-278. DOI:10.4490/algae.2018.33.8.30
[16]
Tatarenkov A, Jönsson R B, Kautsky L, et al. Genetic structure in populations of Fucus vesiculosus (Phaeophyceae) over spatial scales from 10 m to 800 km[J]. Journal of Phycology, 2007, 43(4): 675-685. DOI:10.1111/j.1529-8817.2007.00369.x
[17]
Muhlin J F, Engel C R, Stessel R, et al. The influence of coastal topography, circulation patterns, and rafting in structuring populations of an intertidal alga[J]. Molecular ecology, 2008, 17(5): 1198-1210. DOI:10.1111/j.1365-294X.2007.03624.x
[18]
Kollars N M, Krueger-Hadfield S A, Byers J E, et al. Development and characterization of microsatellite loci for the haploid-diploid red seaweed Gracilaria vermiculophylla[J]. PeerJ, 2015, 3: e1159. DOI:10.7717/peerj.1159
[19]
Krueger-Hadfield S A, Kollars N M, Byers J E, et al. Invasion of novel habitats uncouples haplo-diplontic life cycles[J]. Molecular Ecology, 2016, 25(16): 3801-3816. DOI:10.1111/mec.13718
[20]
Yeh F C, Yang R C, Boyle T. Popgene Version 1.31 Quick User Guide[M]. Canada: University of Alberta and Centre for International Forestry Research, 1999.
[21]
Peakall R. Smouse P E. GenAlEx 6.5:genetic analysis in Excel. Population genetic software for teaching and research-an update[J]. Bioinformatics, 2012, 28: 2537-2539. DOI:10.1093/bioinformatics/bts460
[22]
Takezaki N, Nei M, Tamura K. POPTREE2:Software for constructing population trees from allele frequency data and computing other population statistics with Windows interface[J]. Molecular Biology and Evolution, 2010, 27: 747-752. DOI:10.1093/molbev/msp312
[23]
Pritchard J K, Stephens M, Donnely P. Inference of population structure using multilocus genotype data[J]. Genetics, 2000, 155: 945-959.
[24]
Excoffier L, Laval G, Schneider S. Arlequin (ver. 3.0):An integrated software package for population genetics data analysis[J]. Evolutionary Bioinformatics online, 2005, 1: 47-50.
[25]
Thornber C S, Gaines S D. Population demographics in species with biphasic life cycles[J]. Ecology, 2004, 85(6): 1661-1674. DOI:10.1890/02-4101
[26]
Dyck J L, Dewreede R E. Reproduction and survival in Mazzaella splendens (Gigartinales, Rhodophyta)[J]. Phycologia, 2006, 45(3): 302-310. DOI:10.2216/04-80.1
[27]
Destombe C, Godin J, Remy J M. Viability and dissemination of spermatia of Gracilaria verrucosa (Gracilariales, Rhodophyta)[J]. Hydrobiologia, 1990, 204/205: 219-233. DOI:10.1007/BF00040237
[28]
任国忠, 陈美琴. 江蓠主枝不同部位的生长观察[J]. 海洋科学, 1987, 11(3): 39-42.
Ren Guozhong, Chen Meiqin. Observations on the growth of different part of plant segment of Gracilaria asiasica sp. nov.[J]. Marine Sciences, 1987, 11(3): 39-42.
[29]
Zardi G I, Nicastro K R, Serräo E A, et al. Closer to the rear edge:eco-logical and genetic diversity down the core-edge gradient of a marine macroalga[J]. Ecosphere, 2015, 6: 1-25.
[30]
Krueger-Hadfield S A, Roze D, Mauger S, et al. Intergametophytic selfing and microgeographic genetic structure shape populations of the intertidal red seaweed Chondrus crispus[J]. Molecular Ecology, 2013, 22: 3242-3260. DOI:10.1111/mec.12191
[31]
黄玥.末次冰期以来南海及日本海硅藻及其古环境变化[D].上海: 华东师范大学, 2009.
Huang Yue. Diatom response to changes in palaeoenvironments of the South China Sea and the East Sea (Sea of Japan) since the Last Glacial Maximum[D]. Shanghai: East China Normal University, 2009.
[32]
Yokoya N S, Kakita H, Obika H, et al. Effects of environmental factors and plant growth regulators on growth of the red alga Gracilaria vermiculophylla from Shikoku island, Japan[J]. Hydrobiologia, 1999, 398: 339-347. DOI:10.1023/A%3A1017072508583
[33]
Kostamo K, Korpelainen H, Olsson S. Comparative study on the population genetics of the red algae Furcellaria lumbricalis occupying different salinity conditions[J]. Marine Biology, 2012, 159(3): 561-571. DOI:10.1007/s00227-011-1835-z
[34]
Ogawa T, Ohki K, Kamiya M. High heterozygosity and phenotypic variation of zoids in apomictic Ulva prolifera (Ulvophyceae) from brackish environments[J]. Aquatic Botany, 2014, 120: 185-192.
[35]
刘若愚.中国近海重要生态建群红藻真江蓠的群体遗传多样性[D].青岛: 中国科学院研究院海洋研究所, 2016.
Liu Ruoyu. Genetic diversity of the habitat-forming red alga Gracilaria vermiculophylla along Chinese coasts[D]. Qingdao: Institute of Oceanology, Chinese Academy of Sciences, 2016.
[36]
Freshwater D W, Montgomery F, Greene J K, et al. Distribution and identification of an invasive Gracilaria species that is hampering commercial fishing operations in southeastern North Carolina, USA[J]. Biological Invasions, 2006, 8(4): 631-637. DOI:10.1007/s10530-005-1809-5