GENETIC DIVERSITY OF ATRINA PECTINATA FROM THE NORTHERN SOUTH CHINA SEA BASED ON MITOCHONDRIAL DNA COⅠ SEQUENCE
-
摘要: 为了解南海栉江珧(Atrina pectinata)的群体遗传变异特征, 研究利用线粒体DNA COⅠ(细胞色素氧化酶亚单位Ⅰ)基因部分序列对7个群体共191个栉江珧个体进行群体遗传多样性和遗传结构分析。结果显示: 在600 bp长的序列中, 共检测到113个核苷酸变异位点, 定义了73个单倍型。南海北部栉江珧总体呈现较高的单倍型多样性(0.8996)和较高的核苷酸多样性(0.0257), 但L1组内6个群体(汕头、阳江、湛江、海口、琼海、北海)呈现较高的单倍型多样性(0.8133—0.9286)和较低的核苷酸多样性(0.0033—0.0045)。单倍型系统发育树和网络支系图将7个群体划分为L1和L2(防城港群体)两大类群, 但L1组单倍型并未表现出与地理位置相对应的谱系结构。Fst分析显示, L1组内6个群体间不存在明显的遗传分化(Fst= –0.0200— –0.0055, P>0.05), 但L1组与L2组间存在极显著的遗传分化(Fst=0.8729—0.8821, P<0.01)。L1组的Tajima’s D检验(D= –2.3190, P=0)和Fu’s Fs检验(Fs=–26.8316, P=0)均为显著负值, 核苷酸不配对分布呈明显的单峰分布; L2组的Tajima’s D检验(D=–1.4320, P=0.0565)为不显著负值, Fu’s Fs检验(Fs=4.9540, P=0.9620)为不显著正值。以上数据说明, L1组和L2组栉江珧可能已经分化为两个群体, L1组内6个群体具有频繁的基因交流, 导致了较高的遗传同质性。Abstract: Atrina pectinata is a commercially important species in the Northern South China Sea. This study analzyed genetic diversity and genetic structure of A. pectinata in Northern South China Sea based on 191 individuals of 7 populations using Cytochrome Oxidase Ⅰ(COⅠ) sequence. A 600bp segment of COⅠ genes were sequenced, from which 113 polymorphic sites were tested and 73 haplotypes were defined. The total population in Northern South China Sea had high haplotype diversity (0.8996) and high nucleotide diversity (0.0257). However, 6 populations in L1 group (ST, YJ, ZJ, HK, QH, BH) had high haplotype diversity (0.8133—0.9286) and low nucleotide diversity (0.0033—0.0045). Based on neighbor-joining tree and haplotype network, 7 populations were classified into two groups, named L1 group and L2 group (FCG), and no clustering corresponding to sampling localities was found in L1 group haplotypes. Fst analysis showed no obvious genetic differentiation among 6 populations in L1 group (Fst= –0.0200— –0.0055, P>0.05) and significant genetic differentiation between L1 group and L2 group (Fst=0.8729—0.8821, P<0.01). Both Tajima’s D (D= –2.3190, P=0) and Fu’s Fs (Fs= –26.7990, P=0) of L1 group yielded significant negative values. Unimodal distribution was observed in L1 group. However, in L2 group, Tajima’s D (D= –1.4320, P=0.0565) was insignificant negative value, and Fu’s Fs (Fs=4.9540, P=0.9620) was insignificant positive value. These data implied L1 group and L2 group might have differentiated into two populations. The frequent gene exchange resulted in high genetic homogeneity among 6 populations in L1 group.
-
栉江珧(Atrina pectinata), 隶属于软体动物门(Mollusca), 瓣鳃纲(Lamellibranchia), 贻贝目(Mytiloida), 江珧科(Pinnidae), 江珧属(Atrina)[1], 有带子、大海红、海锨、大海荞麦、土杯、马蹄等俗称。栉江珧为一种广温、广盐的稳居性贝类, 广泛分布于温、热带泥沙质近海地区, 我国的渤海、黄海、东海和南海四大海域均有分布, 其中南海是我国栉江珧分布种类最多、产量最高的海区[1]。栉江珧闭壳肌肥大、肉质鲜嫩, 具有较高的经济价值和药用价值[2], 在我国南海沿岸地区备受欢迎。然而, 栉江珧的人工养殖技术尚不完善, 一般采用半人工养殖方法, 即采捕野生栉江珧苗种、然后选择合适海区进行集中养殖的方式[3]。多年的过度捕捞和人工养殖的相对滞后使栉江珧的野生资源遭到大量破坏, 种质资源明显退化[4]。这与南海沿岸日益上升的需求量相违背, 对南海栉江珧资源进行调查和保护成为亟待解决的科学问题。
细胞色素氧化酶亚单位Ⅰ(Cytochrome Oxidase Ⅰ, COⅠ)是一种线粒体DNA标记, 因其结构简单、进化速度快、母系遗传、易于扩增等特点, 被广泛应用于不同生物群体的遗传多样性分析[5]。近年来有研究者基于COⅠ序列对不同地理位置栉江珧的群体遗传性进行了分析, 如严加坤等[4]调查了长岛、日照、文登、湛江、海口五个城市的野生栉江珧资源, 发现南方群体和北方群体之间存在极大遗传分化, 而南北方群体内部几乎没有遗传分化。Xue等[6]利用COⅠ基因分析了辽宁獐子岛、山东蓬莱、山东刘公岛、山东荣成、山东红岛、山东日照、江苏连云港、浙江舟山、福建福州和韩国海州(Haeju)8个地区的栉江珧资源, 发现来自中国北方海域的7个群体虽然遗传多样性较高, 但没有遗传分化; 而中国和韩国样品间存在明显的遗传分化。以上利用COⅠ基因对栉江珧的研究多聚焦在中国北部海域, 或者更关注南北方海域群体的差异, 而对于栉江珧分布种类最多、产量最高的南海北部海域群体内部的资源调查研究较少。
本研究基于COⅠ基因序列对我国南海北部7个沿海城市的栉江珧群体进行群体遗传学研究, 从遗传多样性、遗传分化和历史动态探讨南海北部栉江珧的系统地理格局, 以期阐明南海北部栉江珧的群体遗传特征, 评估其遗传资源状况, 为资源保护、人工养殖和良种选育提供理论储备。
1. 材料与方法
1.1 样品采集
采集了汕头(ST)、阳江(YJ)、湛江(ZJ)、海口(HK)、琼海(QH)、北海(BH)、防城港(FCG)近海等7个海区野生群体, 共191只栉江珧个体(图 1和表 1)。对所有样品进行形态学鉴定后取其闭壳肌, –80℃保存备用。
![]()
图 1 南海北部栉江珧7个群体的采样地点Figure 1. Sample locations for seven A. pectinata populations from the Northern South China Sea表 1 南海北部栉江珧7个群体的采样信息和遗传多样性参数Table 1. The sampling information and genetic diversity indices for 7 A. pectinata populations from the Northern South China Sea群体
Population样品数量
Sample size单倍型数量
Haplotype number单倍型多样性
Haplotype diversity核苷酸多样性
Nucleotide diversity多态位点数
Polymorphic site汕头ST 25 12 0.8133±0.0804 0.0041±0.002559 18 阳江YJ 28 17 0.9286±0.0359 0.0045±0.002755 20 湛江ZJ 30 20 0.9126±0.0457 0.0043±0.002657 24 海口HK 25 14 0.8700±0.0615 0.0033±0.002129 14 琼海QH 29 14 0.8645±0.0525 0.0041±0.002519 21 北海BH 30 18 0.8483±0.0670 0.0039±0.002413 25 防城港FCG 24 9 0.8148±0.0722 0.0219±0.011376 77 总群体 Total 191 73 0.8996±0.0182 0.0258±0.012763 113 1.2 DNA提取和目的基因的扩增
用E.Z.N.A.R Tissue DNA试剂盒(OMEGA)提取各样品基因组DNA, 1%琼脂糖凝胶电泳检测DNA的纯度和浓度。COⅠ基因片段扩增所用引物序列为: LCO1490 (GGTCAACAAATCATAAA GATATTGG)和HCO2198 (TAAACTTCAGGGT GACCAAAAAATCA)。PCR反应体系为25 μL 2×EasyTaq PCR SuperMix (TransGen)、正反引物(10 μmol/L)各1 μL, 基因组DNA (100 ng/μL) 2 μL, ddH2O补齐50 μL。PCR反应程序为: 94℃预变性3min; 94℃变性30s, 57℃退火30s, 72℃延伸1min, 35个循环; 最后72℃延伸10min。PCR产物经1%琼脂糖凝胶电泳检测后由生工测序。
1.3 数据处理
利用DNAStar 7.1和Clustal X 1.81[7]进行序列校对和编辑; 利用DnaSP 5.10.01[8]计算单倍型数目、单核苷酸变异位点数和多态简约信息位点数; 利用Arlequin 3.5.1.2[9]统计序列的平均碱基组成, 计算单倍型多样性(h)和核苷酸多样性(π)等群体遗传学参数; 利用MEGA 6.06[10]计算转换百分比、颠换百分比和转换/颠换比率(Ts/Tv), 选择DNA序列的最佳替换模型。利用MrBayes 3.2.6构建单倍型系统进化树, 利用Network5.0.0.1(http://www.fluxus-engineering.com/share-net.htm)中的中介连接网络法构建单倍型网络支系图[11]
群体遗传结构采用3种方法进行分析: (1)采用分化固定指数Fst评价两两群体间的遗传差异, 检验结果的显著性; (2)根据样品来源、遗传学参数和Fst, 将7个栉江珧群体划分为2个组: 汕头(ST)、阳江(YJ)、湛江(ZJ)、海口(HK)、琼海(QH)和北海(BH)6个栉江珧群体为一组(L1组), 防城港群体(FCG)单独为一组(L2组), 通过AMOVA分析估算遗传变异在组群间、组群内群体间和群体内的分布, 进行显著性分析; (3)通过Exact检验检测单倍型在群体间分布频率的差异。
群体历史动态采用2种方法进行检测: (1)采用Tajima’s D检验和Fu’s Fs检验检测中性假说是否成立, 并进行显著性分析; (2)基于核苷酸不配对分布检验南海北部栉江珧群体是否存在数量和空间扩张, 并通过拟合优度检验(Goodness of fit test)检测核苷酸不配对的观测分布和预期分布之间的一致性, 并估算群体历史扩张时间。在本研究中, COⅠ基因的分化速率采用2.4%/MY (Million year, MY)[6], 世代时间设置为1年[12]。
2. 结果
2.1 序列分析
本研究共获得7个群体191条长度为600 bp的COⅠ基因部分序列, 其G、C、A、T的平均含量分别为24.05%、14.81%、19.18%和41.96%。共检测到113个核苷酸变异位点(18.8%), 包括14个单核苷酸变异位点, 99个多态简约信息位点, 无插入和缺失位点。转换数共有107次, 明显高于颠换数(19次), Ts/Tv=15.33。转换和颠换在G、C、A、T之间均有发生, 其中C-T (53.68%)转换多于G-A (40.88%)转换, T-G (1.79%)和T-A (1.66%)颠换稍多于C-G (1.05%)和C-A (0.92%)颠换。突变位点均匀地分布在COⅠ基因上, 但在群体间的分布差异较大, 77个多态位点来源于FCG群体。
2.2 单倍型和遗传多样性分析
COⅠ单倍型在7个群体中的分布见表 1和表 2。113个核苷酸变异位点共定义73个单倍型(GenBank登录号为MH536017—MH536089), 包括58个独有单倍型和15个共享单倍型。在共享单倍型中, 3个单倍型(H23、H36和H44)为6个群体所共有, 1个单倍型(H2)为5个群体所共享, 1个单倍型(H29)为3个群体所共享, 10个单倍型(H1、H18、H20、H34、H47、H51、H52、H54、H61和H62)为2个群体所共享。H44出现的频率最高(占总样品的31.41%), 这表明这个单倍型是南海北部栉江珧在长期进化过程中形成的较为稳定的优势基因型。值得注意的是, 防城港群体具有9个单倍型(H65-H73), 均为独有单倍型。
表 2 南海北部栉江珧7个群体的单倍型分布情况Table 2. Haplotype distribution of 7 A. pectinata populations from the Northern South China Sea单倍型Haplotype 汕头ST 阳江YJ 湛江ZJ 海口HK 琼海QH 北海BH 防城港FCG 总群体Total H1 2 1 3 H2 1 2 1 2 2 8 H3 1 1 H4 1 1 H5 1 1 H6 1 1 H7 1 1 H8 1 1 H9 1 1 H10 1 1 H11 1 1 H12 1 1 H13 1 1 H14 1 1 H15 1 1 H16 1 1 H17 1 1 H18 1 1 2 H19 1 1 H20 1 1 2 遗传多样性参数如表 1所示, 南海北部栉江珧总群体的单倍型多样性和核苷酸多样性分别为0.8996和0.0257, 表明南海北部栉江珧群体的遗传多样性目前仍处于较高的水平。在7个群体中, L1组6个群体的单倍型多样性为0.8133—0.9286, 核苷酸多样性为0.0033—0.0045, 表现出高h低π遗传多样性模式; 而防城港群体的单倍型多样性为0.8148, 核苷酸多样性较高(0.0219), 表现出高h高π遗传多样性模式。
2.3 群体遗传学关系分析
以旗江珧(Atrina vexillum) COⅠ基因序列为外群, 构建栉江珧73个单倍型的贝叶斯推理树(图 2)。结果显示, L1组6个群体的各单倍型交叉出现在系统发育树中, 并未表现出与地理位置相对应的谱系结构; L2组的独有单倍型(H65-H73)聚为一个独立的进化分支, 显示了较远的遗传距离。单倍型网络支系图(图 3)具有1个主体单倍型(H44), 该主体单倍型在6个群体中均有分布, 且出现频率较高, 说明它可能是起源于母系祖先的主体单倍型; 其他单倍型呈辐射状分布在主体单倍型周围, 其中L2组的9个独有单倍型作为一个独立分支, 与其他单倍型具有较远的遗传距离。这暗示着南海北部栉江珧北部湾的防城港群体与其他群体间存在显著的遗传结构。
![]()
图 2 南海北部栉江珧COⅠ单倍型贝叶斯推理树外类群为旗江珧(KU987243), 各分支上的数字为后验概率值, 分支上无数字表示1Figure 2. Bayesian inference tree for A. pectinata from the Northern South China Sea based on COⅠAtrina vexillum was used as outgroup. Posterior probability values were shown adjacent to each internode. 1 was not shown![]()
图 3 南海北部栉江珧COⅠ单倍型支系图Figure 3. Haplotype network of A. pectinata from the Northern South China Sea based on COⅠ2.4 群体遗传结构
由表 3可知, L1组6个群体间的Fst值均较小(–0.0200— –0.0055), 且统计学差异不显著(P>0.05); 但L2组与L1组6个群体间的Fst值均较大(0.8729—0.8821), 且统计学差异极显著(P<0.01)。分子方差分析(Analysis of molecular variance, AMOVA) 结果显示(表 4), 组群间、组群内群体间和群体内个体间的遗传变异比例分别为97.01%、–0.1%和3.09%, 组群间、组群内群体间和群体内个体间的固定系数Fct、Fsc和Fst分别为0.9701 (P>0.05)、–0.0324 (P>0.05)和0.9691 (P<0.01)。Fst值和AMOVA结果表明, L1组6个群体间在分子水平不存在明显的遗传分化, 但L2组与L1组6个群体间存在显著遗传分化。Exact检验结果显示, 单倍型在L1组6个群体间的分布频率并无显著差异, 但在L2组与L1组间差异显著。
表 3 南海北部栉江珧7个群体间的遗传分化指数FstTable 3. The pairwise Fst among seven populations of A. pectinata from the Northern South China Sea群体Populations 汕头ST 阳江YJ 湛江ZJ 海口HK 琼海QH 北海BH 防城港FCG 汕头ST 阳江YJ –0.0200 湛江ZJ –0.0165 –0.0099 海口HK –0.0156 –0.0123 –0.0065 琼海QH –0.0190 –0.0067 –0.0072 –0.0187 北海BH –0.0055 –0.0110 –0.0104 –0.0094 –0.0091 防城港FCG 0.8729 0.8755 0.8795 0.8772 0.8792 0.8821 表 4 南海北部栉江珧7个群体间的AMOVA分析Table 4. Analysis of molecular variance (AMOVA) among seven populations of A. pectinata from the Northern South China Sea变异来源
Source of variation自由度
df方差总和
Sum of square变异组分
Variance component变异比例
Percentage of variation (%)固定系数
Fixation index组群间Among groups 1 9862.6640 234.9600 Va 97.01 Fct=0.9701 组群内群体间Among populations within groups 5 4.7620 –0.2348 Vb –0.1 Fsc=–0.0324 群体内Within populations 184 1376.4930 7.4809 Vc 3.09 Fst=0.9691 2.5 群体历史动态
由于FCG群体与其他6个群体间具有显著的遗传分化, 故将FCG群体单独进行历史动态分析。中性检验结果显示, FCG群体所在L2组的Tajima’s D检验(D= –1.4320, P=0.0565)为不显著负值, Fu’s Fs检验(Fs=4.9540, P=0.9620)为不显著正值, 说明FCG群体进化符合中性假说。L1组6个群体间无显著的遗传分化, 故将其合并为一个大群体进行历史动态分析。结果显示, Tajima’s D检验(D= –2.3190, P=0)和Fu’s Fs检验(Fs= –26.8316, P=0)均为显著负值。基于数量扩张模型和空间扩张模型的核苷酸不配对分布均呈单峰分布(图 4), 且观测分布和期望分布拟合优度较好(SSD), Hri值小且不显著(表 5)。中性检验和核苷酸不配对分布结果均提示南海北部栉江珧L1组6个群体在历史上可能发生过数量和栖息地扩张事件。数量扩张模型和空间扩张模型下的τ值分别为1.2830和1.0880 (表 5), 据此推算南海北部栉江珧群体扩张时间约在75500—89000年前。
![]()
图 4 南海北部栉江珧总群体COⅠ核苷酸不配对分布图A. 数量扩张模型; B. 空间扩张模型Figure 4. Mismatch distributions of A. pectinata from the Northern South China Sea based on COⅠA. sudden expansion model; B. spatial expansion model表 5 南海北部栉江珧总群体COⅠ核苷酸不配对分布参数Table 5. Mismatch distribution indices of A. pectinata from the Northern South China Sea based on COⅠ模型Model 数量扩张Sudden expansion 空间扩张Spatial expansion 拟合优度检验
Test of goodness-of-fitSSD PSSD Hri PHri SSD PSSD Hri PHri 0.0014 0.7130 0.0176 0.6810 0.0015 0.6370 0.0176 0.7620 群体动态参数
Demographic parametersτ θ0 θ1 t/年 τ θ M t/年 1.2830 1.4310 21.7580 89097 1.0808 1.5680 170.1250 75555 3. 讨论
3.1 群体遗传多样性
单倍型多样性(h)和核苷酸多样性(π)是衡量物种遗传多样性的2个重要指标[13]。在本研究中, 南海北部栉江珧L1组群呈现出较高的单倍型多样性(0.8133—0.9286)和较低的核苷酸多样性(0.0033—0.0045), 即高h低π模式。这种高h低π的遗传参数特征通常是由一个较小的有效群体经过近期快速扩张成一个大的群体所引起的[14]。在末次冰盛期, 南海海平面下降约100—120 m, 变成一个半封闭的内陆海(图 1), 成为众多海洋生物的冰期避难所[15, 16]。在冰期恶劣的环境下, 南海栉江珧的分布范围和数量可能经历了一定程度的收缩; 随着间冰期气温回升, 海平面上升, 产生大量近岸栖息地, 冰期避难所内残留的栉江珧可能发生快速扩张重新进入南海大陆架。此过程南海栉江珧群体积累了单倍型的多样性, 但尚缺乏足够的时间积累核苷酸序列的多样性[11]。而防城港群体显示了较高的单倍型多样性(0.8148)和核苷酸多样性(0.0219), 这暗示着北部湾丰富的生态资源和相对稳定的生存环境对栉江珧保持遗传多样性产生了重要影响。
单倍型多样性和核苷酸多样性分析结果表明, 南海北部栉江珧群体的遗传多样性目前仍处于较高的水平, 这可能与栉江珧广阔的生境范围有关。栉江珧在我国东海、渤海、黄海和南海海域均有分布, 广阔的生境范围使其面临较小的自然选择压力, 积累较多的遗传变异。虽然自第四纪以来, 全球温度的大幅度升降使许多物种的遗传多样性损失严重[17], 但栉江珧作为一种底栖动物可能受到的影响较小, 保留了较高的遗传多样性。
3.2 群体遗传结构
在本研究中, AMOVA分析、Fst值、单倍型系统发育树和网络支系图结果均表明, L1组6个栉江珧群体的遗传差异不显著, 不存在显著的群体遗传结构。造成L1组6个群体遗传分化不显著的原因可能有三点: 一是扩散生活史特征。对于有幼虫阶段的海洋生物, 浮游期(Pelagic larval duration, PLD)对于基因交流和遗传结构具有重要影响[18, 19]。栉江珧的浮游期约30d左右, 比栉孔扇贝(Chlamys farreri)(约15d)、中国蛤蜊(Mactra chinensis)(约10d)等双壳类具有更长的浮游期[18], 这暗示着栉江珧具有较强的扩散能力, 可以越过地理环境和水深条件, 产生更多的基因交流。二是洋流的输送作用。南海沿岸栉江珧5月份进入成熟期和排放期, 排放期可一直延续到12月份, 在此期间栉江珧的性腺分批多次排放[20]。西南季风漂流、南海暖流、粤东沿岸流和海南岛以东沿岸流等北向海流[21, 22]有利于QH、ZJ等地方群体卵和幼虫的的北上扩散; 而东北季风漂流、广东沿岸流和黑潮南海分支等南向海流[21, 22]则可促进ST群体向南迁移, 从而进行基因交流。三是目前一般认为海洋生物是从海洋到陆地的迁移方式, 陆地周围的淡水、浊流和人类活动对生物体的扩散形成一个地理和生态的障碍[23]。但栉江珧适应能力强, 从潮下带到水下100 m、泥质或砂质海区皆可分布[24, 25], 这意味着栉江珧比其他贝类具有更强的克服地理和生态障碍的能力, 能够被洋流输送更远的距离, 产生更多的基因交流。
值得一提的是防城港群体与其他6个群体间产生了明显的遗传分化, 具有显著的遗传结构, 有3个可能的解释: 一是北部湾是一个半封闭的海域, 位于大陆架边缘, 于古近系断陷期形成, 南海栉江珧很可能在此时期发生了分化[26], 其中L1组群随后进行了大规模扩散, 而防城港群体由于地理位置所限, 扩散程度则相对较小, 经过多年积累形成了遗传分化。二是基于生态竞争不相容原理。当两个相近的物种或个体共存于一个相对狭小的生境中时, 他们可能会划分生存空间和食物资源, 以最大程度地减少生态位重叠和竞争[27, 28]。长期空间和食物的差异导致两个物种最终产生遗传分化, 尤其体现在防城港群体和北海群体间。三是北部湾的自然资源对于栉江珧群体而言非常丰富, 尤其是红河的汇入带来了丰富的海底沉积物。北部湾的红树林资源亦很丰富, 空间异质性强, 可以作为双壳类幼虫的天然屏障, 不但给幼虫提供了多样性的食物和小生境, 也减少其他物种对幼虫的掠夺[26, 29], 从而使北部湾狭小空间内产生多样性的系统分化成为可能。
3.3 资源管理
目前, 南海北部栉江珧种质资源明显退化[4], 而物种的遗传多样性是生存和进化的前提, 丰富的遗传多样性意味着较大的进化潜力, 也是生物适应环境变化的必要条件[30]。本研究表明, 北部湾防城港栉江珧群体COⅠ的遗传多样性水平较高; 汕头、阳江、湛江、海口、琼海和北海6个群体的遗传多样性相对偏低, 其中海口群体的遗传多样性最低。为避免栉江珧种质资源严重衰退的现象, 应当加强栉江珧的管理和保护, 尤其是面临较大捕捞压力的南海北部沿岸群体。
-
![]()
图 1 南海北部栉江珧7个群体的采样地点
Figure 1. Sample locations for seven A. pectinata populations from the Northern South China Sea
![]()
图 2 南海北部栉江珧COⅠ单倍型贝叶斯推理树
外类群为旗江珧(KU987243), 各分支上的数字为后验概率值, 分支上无数字表示1
Figure 2. Bayesian inference tree for A. pectinata from the Northern South China Sea based on COⅠ
Atrina vexillum was used as outgroup. Posterior probability values were shown adjacent to each internode. 1 was not shown
![]()
图 3 南海北部栉江珧COⅠ单倍型支系图
Figure 3. Haplotype network of A. pectinata from the Northern South China Sea based on COⅠ
![]()
图 4 南海北部栉江珧总群体COⅠ核苷酸不配对分布图
A. 数量扩张模型; B. 空间扩张模型
Figure 4. Mismatch distributions of A. pectinata from the Northern South China Sea based on COⅠ
A. sudden expansion model; B. spatial expansion model
表 1 南海北部栉江珧7个群体的采样信息和遗传多样性参数
Table 1 The sampling information and genetic diversity indices for 7 A. pectinata populations from the Northern South China Sea
群体
Population样品数量
Sample size单倍型数量
Haplotype number单倍型多样性
Haplotype diversity核苷酸多样性
Nucleotide diversity多态位点数
Polymorphic site汕头ST 25 12 0.8133±0.0804 0.0041±0.002559 18 阳江YJ 28 17 0.9286±0.0359 0.0045±0.002755 20 湛江ZJ 30 20 0.9126±0.0457 0.0043±0.002657 24 海口HK 25 14 0.8700±0.0615 0.0033±0.002129 14 琼海QH 29 14 0.8645±0.0525 0.0041±0.002519 21 北海BH 30 18 0.8483±0.0670 0.0039±0.002413 25 防城港FCG 24 9 0.8148±0.0722 0.0219±0.011376 77 总群体 Total 191 73 0.8996±0.0182 0.0258±0.012763 113 
下载: 导出CSV
表 2 南海北部栉江珧7个群体的单倍型分布情况
Table 2 Haplotype distribution of 7 A. pectinata populations from the Northern South China Sea
单倍型Haplotype 汕头ST 阳江YJ 湛江ZJ 海口HK 琼海QH 北海BH 防城港FCG 总群体Total H1 2 1 3 H2 1 2 1 2 2 8 H3 1 1 H4 1 1 H5 1 1 H6 1 1 H7 1 1 H8 1 1 H9 1 1 H10 1 1 H11 1 1 H12 1 1 H13 1 1 H14 1 1 H15 1 1 H16 1 1 H17 1 1 H18 1 1 2 H19 1 1 H20 1 1 2
下载: 导出CSV
表 3 南海北部栉江珧7个群体间的遗传分化指数Fst
Table 3 The pairwise Fst among seven populations of A. pectinata from the Northern South China Sea
群体Populations 汕头ST 阳江YJ 湛江ZJ 海口HK 琼海QH 北海BH 防城港FCG 汕头ST 阳江YJ –0.0200 湛江ZJ –0.0165 –0.0099 海口HK –0.0156 –0.0123 –0.0065 琼海QH –0.0190 –0.0067 –0.0072 –0.0187 北海BH –0.0055 –0.0110 –0.0104 –0.0094 –0.0091 防城港FCG 0.8729 0.8755 0.8795 0.8772 0.8792 0.8821
下载: 导出CSV
表 4 南海北部栉江珧7个群体间的AMOVA分析
Table 4 Analysis of molecular variance (AMOVA) among seven populations of A. pectinata from the Northern South China Sea
变异来源
Source of variation自由度
df方差总和
Sum of square变异组分
Variance component变异比例
Percentage of variation (%)固定系数
Fixation index组群间Among groups 1 9862.6640 234.9600 Va 97.01 Fct=0.9701 组群内群体间Among populations within groups 5 4.7620 –0.2348 Vb –0.1 Fsc=–0.0324 群体内Within populations 184 1376.4930 7.4809 Vc 3.09 Fst=0.9691
下载: 导出CSV
表 5 南海北部栉江珧总群体COⅠ核苷酸不配对分布参数
Table 5 Mismatch distribution indices of A. pectinata from the Northern South China Sea based on COⅠ
模型Model 数量扩张Sudden expansion 空间扩张Spatial expansion 拟合优度检验
Test of goodness-of-fitSSD PSSD Hri PHri SSD PSSD Hri PHri 0.0014 0.7130 0.0176 0.6810 0.0015 0.6370 0.0176 0.7620 群体动态参数
Demographic parametersτ θ0 θ1 t/年 τ θ M t/年 1.2830 1.4310 21.7580 89097 1.0808 1.5680 170.1250 75555
下载: 导出CSV
-
[1] 李明爽. 栉江珧微卫星标记的筛选及其遗传多样性的研究. 硕士学位论文, 南京农业大学, 南京. 2008 Li M S. Identification of microsatellite markers and analysis of genetic diversity of Atrina pectinata [D]. Thesis for Master of Science. Nanjing Agricultural University, Nanjing. 2008
[2] 王如才, 王昭萍. 海水贝类养殖学. 青岛: 中国海洋大学出版社. 2008, 532—534 Wang R C, Wang Z P. Science of Marine Shellfish Culture [M]. Qingdao: China Ocean University Press. 2008, 532—534
[3] 郑言鑫, 杨爱国, 吴彪, 等. 栉江珧(Atrina pectinata)催产方法及幼虫培养条件. 渔业科学进展, 2015, 36(6): 127—133 Zhen Y X, Yang A G, Wu B, et al. Studies on the spawning inducement and the conditions for the larval culture of Atrina pectinata [J]. Progress in Fishery Sciences, 2015, 36(6): 127—133
[4] 严加坤, 杨爱国, 周丽青, 等. 我国不同地理群体栉江珧遗传多样性及系统发生分析. 渔业科学进展, 2013, 34(2): 29—35 Yan J K, Yang A G, Zhou L Q, et al. The genetic diversity and phylogenetic studies of five populations of Atrina pectinata [J]. Progress in Fishery Sciences, 2013, 34(2): 29—35
[5] Lü J, Wu S, Zhang Y, et al. Assessment of four DNA fragments (COI, 16S rDNA, ITS2, 12S rDNA) for species identification of the Ixodida (Acari: Ixodida) [J]. Parasit Vectors, 2014, 7: 93
[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(4): e95436
[7] Thompson J D, Gibson T J, Plewniak F, et al. The Clustal_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools [J]. Nucleic Acids Research, 1997, 25(24): 4876—4882
[8] Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data [J]. Bioinformatics, 2009, 25(11): 1451—1452
[9] Excoffier L, Lischer H E L. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows [J]. Molecular Ecology Resources, 2010, 10(3): 564—567
[10] Tamura K, Stecher G, Peterson D, et al. MEGA6: molecular evolutionary genetics analysis version 6.0 [J]. Molecular Biology and Evolution, 2013, 30(12): 2725—2729
[11] 牛素芳, 苏永全, 王军, 等. 福建近海竹荚鱼线粒体DNA控制区和细胞色素b遗传多态性. 中国水产科学, 2011, 18(1): 66—74 Niu S F, Su Y Q, Wang J, et al. Genetic polymorphism of mitochondrial control region and cyt b in Trachurus japonicus from Fujian coastal waters [J]. Journal of Fishery Sciences of China, 2011, 18(1): 66—74
[12] 邓景耀, 赵传絪, 唐启升, 等. 海洋渔业生物学. 北京: 农业出版社. 1991, 485—518 Deng J Y, Zhao C Y, Tang Q S, et al. Marine Fisheries Biology [M]. Beijing: Agriculture Press. 1991, 485—518
[13] 陈骁, 杨圣云, 潘聪. 我国南部海域条纹斑竹鲨线粒体DNA控制区遗传多态性研究. 海洋学报, 2008, 30(6): 115—121 Chen X, Yang S Y, Pan C. Genetic polymorphism of mitochondrial control region in Chiloscyllium plagiosum from coastal waters of southern China [J]. Acta Oceanologica Sinica, 2008, 30(6): 115—121
[14] Grant W S, Bowen B W. Shallow population histories in deep evolutionary lineages of marine fishes: insights from sardines and anchovies and lessons for conservation [J]. Journal of Heredity, 1998, 89(5): 415—426
[15] Zhu Y X, Cheng Q Q, Rogers S M. Genetic structure of Scomber japonicus (Perciformes: Scombridae) along the coast of China revealed by complete mitochondrial cytochrome b sequences [J]. Mitochondrial DNA Part A, 2016, 27(6): 3828—3836
[16] He L J, Mukai T, Chu K H, et al. Biogeographical role of the Kuroshio current in the amphibious mudskipper Periophthalmus modestus indicated by mitochondrial DNA data [J]. Scientific Reports, 2015, 5: 15645
[17] Hewitt G M. Some genetic consequences of ice ages, and their role in divergence and speciation [J]. Biological Journal of the Linnean Society, 1996, 58(3): 247—276
[18] Selkoe K, Toonen RJ. Marine connectivity: a new look at pelagic larval duration and genetic metrics of dispersal [J]. Marine Ecology Progress Seriesm, 2011, 436: 291—305
[19] 宋娜, 王航宁, 刘淑德, 等. 青岛近海路氏双髻鲨群体遗传学分析. 水生生物学报, 2017, 41(3): 603—608 Song N, Wang H N, Liu S D, et al. Genetic study on Sphyrna lewini from Qingdao coastal waters [J]. Acta Hydrobiologica Sinica, 2017, 41(3): 603—608
[20] 王梅芳, 余祥勇, 叶富良. 北部湾及附近海域栉江珧性腺发育研究. 广西科学, 2000, 7(2): 140—143, 157 Wang M F, Yu X Y, Ye F L. The gonad development of Pinna (Atrina) pectinata in Beibu Gulf and adjoining coast [J]. Guangxi Sciences, 2000, 7(2): 140—143, 157
[21] 葛倩. 晚末次冰期以来南海古环境和古气候记录. 硕士学位论文, 中国地质大学, 武汉. 2010 Ge Q. Paleoenvironmental and Paleoclimatic records from the South China Sea since the late last glacical [D]. Thesis for Master of Science. China University of Geosciences, Wuhan. 2010
[22] 乔方利主编. 中国区域海洋学-物理海洋学. 北京: 海洋出版社. 2012, 310—350 Qiao F L (Eds.), Regional Oceanography of China Seas-Physical Oceanography [M]. Beijing: Ocean Press. 2012, 310—350
[23] Claremont M, Williams S T, Barraclough T G, et al. The geographic scale of speciation in a marine snail with high dispersal potential [J]. Journal of Biogeography, 2011, 38(6): 1016—1032
[24] 王祯瑞. 中国动物志: 软体动物门双壳纲贻贝目. 北京: 科学出版社. 1997, 362 Wang Z R. Fauna Sinica, Phylum Mollusca, Order Mytiloida [M]. Beijing: Science Press. 1997, 362
[25] An H S, Kim B H, Lee J W, et al. Comparison between wild and hatchery populations of Korean Pen Shell (Atrina pectinata) using microsatellite DNA markers [J]. International Journal of Molecular Sciences, 2011, 12(9): 6024—6039
[26] Liu J, Li Q, Kong L, et al. Cryptic diversity in the pen shell Atrina pectinata (Bivalvia: Pinnidae): high divergence and hybridization revealed by molecular and morphological data [J]. Molecular Ecology, 2011, 20(20): 4332—4345
[27] Chen G, Hare M P. Cryptic ecological diversification of a planktonic estuarine copepod, Acartia tonsa [J]. Molecular Ecology, 2008, 17(6): 1451—1468
[28] Zhao K, Duan Z Y, Peng Z G, et al. The youngest split in sympatric schizothoracine fish (Cyprinidae) is shaped by ecological adaptations in a Tibetan Plateau glacier lake [J]. Molecular Ecology, 2009, 18(17): 3616—3628
[29] Laegdsgaard P, Johnson C. Why do juvenile fish utilize mangrove habitats [J]? Journal of Experimental Marine Biology and Ecology, 2001, 257(2): 229—253
[30] 蒙子宁, 庄志猛, 金显仕, 等. 黄海带鱼、小带鱼RAPD和线粒体16S rRNA基因序列变异分析. 自然科学进展, 2003, 13(11): 1170—1176 Meng Z N, Zhuang Z M, Jin X S, et al. Variation analysis of RAPD and mitochondrial DNA 16S rRNA sequence of Trichiurus japonicus and Eupleurogrammus muticus from the Yellow Sea [J]. Progress in Natural Science, 2003, 13(11): 1170—1176
计量
- 文章访问数:
- HTML全文浏览量:
- PDF下载量: