遗传多样性是物种生存与进化的物质基础，开展鱼类遗传多样性研究可以为遗传育种和种质资源保护提供重要参考(Ward, 2000)。动物线粒体属于母系遗传，无遗传物质的重组，且其进化速度比核基因快，较少的样本就可较好地反映该种群的遗传结构，是动物遗传多样性研究中理想的分子标记(Xu et al, 2014)。而且，动物线粒体的COⅠ序列具有一定的进化变异，可以较好地反映同种生物的不同群体间的遗传差异，除了被广泛应用于DNA条形码的研究中，也适用于遗传多样性的分析研究(柳淑芳等, 2016; Chang et al, 2016; Xu et al, 2012; 李大命等, 2017)。

前鳞 (Liza affinis)隶属鲻形目(Mugiliformes)、鲻科(Mugilidae)，是分布于西北太平洋的日本南部及中国台湾海域等的洄游性鱼类。前鳞为杂食性鱼类，主要以藻类、有机碎屑为食，大多栖息在河口处，具有广盐性且对淡水的耐受性较强(孙典荣等, 2013; Lin et al, 1999)。同时，前鳞的生态位低、生长快、肉味鲜甜，具有可观的渔业价值(Richardson et al, 2011; 孟玮等, 2007)。目前，关于中国前鳞的研究主要为集中在物种鉴定方面，其遗传多样性研究仅涉及浙江沿岸(许则滩, 2015)。为了解中国前鳞的遗传背景，为其种质资源的可持续开发利用提供科学依据，该研究测定了采自中国4省5个地点共计82尾样本的COⅠ序列，并展开相关的遗传分析。

为了研究中国前鳞群体在不同海域(东海与南海)中的遗传多样性，及海峡(琼州海峡与台湾海峡)隔离对其遗传结构的影响，结合采样实际情况，集82尾前鳞样本，分别为广西东兴(26尾)、广东闸坡(6尾)、广东吴川(4尾)、福建平潭(23尾)、浙江岱山(10尾)和浙江洞头(13尾)。参照乐小亮等(2010)方法提取82尾样本的DNA，使用实验室自行设计的引物COⅠF: 5'-TGTAAAACGACGGCCAGTCCTGTGGC AATYACDCGCTGAT，COⅠR: 5'-CAGGAAACAGC TATGACNACYTCNGGRTGNCCRAAGAA与其匹配并进行PCR扩增，将条带明亮的PCR产物送华大基因公司测序。

运用MEGA 7.0(Kumar et al, 2016)对测序峰图进行人工校对，并计算碱基组成、颠换率。通过DnaSP 5.10计算单倍型数、单倍型多样性(H_d)、核苷酸多样性(π) (Librado et al, 2009)。使用Network (Bandelt et al, 1999)分析单倍型的分布情况，通过Arlequin 3.5进行AMOVA分析、核苷酸错配分析、Tajima's D和Fu's Fs中性检验，计算群组间的F_st值、SSD值、R_g值、τ值等参数(Excoffier et al, 2010)。运用SAMOVA 2.0进行空间分子变异分析。根据公式[T=(τ/2μk)×代时]计算种群扩张时间，其中，τ、μ、k、T依次代表种群扩张时间参数、序列的变异速率[线粒体COⅠ基因进化速率为(1%~3%)/百万年(Wang et al, 2013)]、序列长度以及种群扩张时间(Rogers et al, 1992)。是属的代表性鱼类，个体较大，与前鳞的生态习性相近，其性成熟需要2~3年(石琼等, 2015)。大多数鱼类因性成熟而迎来个体生长的停滞(钱国英, 2000)，同属的前鳞个体较小，故其性成熟所需时间可能短于鱼，但目前尚没有确切时间的报道，估算其大概为1~3年，即公式中的代时取值为1~3。

1.3 中国前鳞群体划分

为了方便比较分布在不同地理位置的前鳞群体间的遗传差异，该研究将按以下方式对中国前鳞群体进行划分、分析。按采样地点把中国前鳞划分为东兴群体、平潭群体、闸坡群体、吴川群体、岱山群体和洞头群体，其中，因闸坡和吴川的地理位置相邻且两地的样本数量较少，故将在两地所采集的样本合并为1个采样群体，在下文简称为广东群体。根据群体所处海域将中国前鳞划分为东海群体(岱山、洞头和平潭群体)及南海群体(广东和东兴群体)。根据海峡的隔离将中国前鳞划分为东兴群体、广东群体和东海群体。

2.1 前鳞线粒体COⅠ序列特征和遗传多样性

在全长为712 bp的序列中，没有发现碱基的插入和缺失，碱基A、T、C和G的含量分别为23.5%、29.4%、28.8%和18.3%，A+T含量(52.9%)略高于C+G含量(47.1%)，与其他硬骨鱼类COⅠ序列碱基组成特征基本一致。其中，共有18个多态位点、6个简约信息位点，颠换与转换的比值为4.7，表明序列的突变尚未达到饱和，适合进行系统发育分析(Kumar, 2005)。此外，还发现21个单倍型，总体遗传多样性偏低(H_d=0.4840±0.0700，π=0.0010±0.0002)，其中，洞头群体最高，岱山群体最低，东海群体的遗传多样性高于南海群体，具体的遗传多样性情况见表 1。

表 1(Tab.1)

表 1 中国沿海前鳞群体的遗传多样性 Tab.1 Genetic diversity statistics of Liza affinis in coastal waters of China 群体所处的地理位置 Region of population 数量 Size 变异位点数 Variable site 单倍型数量 Nh 单倍型多样性 Hd 核苷酸多样性 π 岱山Daishan 10 2 5 0.2000±0.1540 0.0006±0.0004 洞头Dongtou 13 9 12 0.8080±0.1130 0.0021±0.0006 平潭Pingtan 23 8 10 0.5850±0.1220 0.0013±0.0004 吴川, 闸坡Wuchua, Zhapo 10 3 6 0.3780±0.1810 0.0008±0.0005 东兴Dongxing 26 4 10 0.3510±0.1170 0.0005±0.0002 东海East China Sea 46 15 23 0.5790±0.0890 0.0014±0.0003 南海South China Sea 36 6 13 0.3540±0.1020 0.0006±0.0002 总体Total 82 18 32 0.4840±0.0700 0.0010±0.0002

表 1 中国沿海前鳞群体的遗传多样性 Tab.1 Genetic diversity statistics of Liza affinis in coastal waters of China

单倍型最大简约性网络图(图 1)呈现星状结构，其中，Hap1的数量最多(59)，由全部群体共享。Hap3的数量为3，由东兴群体与平潭群体共享。Hap11的数量为2，由平潭群体与洞头群体共享，其余均为各个群体特有的单倍型。也就是说，不同地理来源的个体相互聚集分布，不同的地区中存在相同的单倍型，表明前鳞在中国沿海范围内没有明显的地理聚群。

图 1(Fig. 1)

图 1 前鳞单倍型网络图 Fig.1 Parsimony network of Liza affinis COⅠhaplotypes 图中每个圆圈分别代表一种单倍型，枝上短线数量代表突变步骤的数量 Each circle represents a unique haplotype, the number of short lines on the branch represents the number of mutation steps

在AMOVA分析(表 2)中，按照上文的群体划分进行分析：东兴、平潭、闸坡、广东、洞头等5个群体间的F_st=-0.0221~0.0141(P > 0.05)，组间的遗传变异比例达-0.26%。东海群体与南海群体间的F_st= 0.0020 (P > 0.05)，组间的遗传变异比例达0.46%。东兴、广东、东海3个群体间的F_st= -0.0159~0.0141 (P > 0.05)，组间的遗传变异比例为-0.6%。上述各类划分中，群体间的F_st值皆小于0.05，均表现出不显著水平(P > 0.05)，且组间的遗传变异比例均属于极低水平，表明中国前鳞群体在各个地区、海域、海峡间没有明显的遗传分化，遗传变异主要集中在群体中的个体间。而SAMOVA分析(表 2)中，82条序列被软件自动划分为2组，即洞头群体独立归为1组，其余群体为1组，测得组间的变异比例只有2.87%，同样表明组间分化程度很低，与F_st及AMOVA分析的结果相吻合。

表 2(Tab.2)

表 2 前鳞种群结构的分子方差分析及SAMOVA分析 Tab.2 AMOVA and SAMOVA analysis of Liza affinis 分析方法 Method of analysis 分组方式 Method of grouping 分析参数 Parameter of analysis 群组间 Among groups 群组中群体间 Among populations within groups 群体内部 Within populations 总体 Total AMOVA Group (5个群体 5 populations) 自由度df 4.0000 77.0000 81.0000 平方和Sum of squares 1.4280 28.6450 30.0730 变异组成Variance components 0.0010 0.3720 0.3710 变异比例Percentage of variation –0.2600 100.2600 Group1(东海 East China Sea) 自由度df 1.0000 3.0000 77.0000 81.0000 平方和Sum of squares 0.4000 1.0270 28.6450 30.0730 Group2 (南海 South China Sea) 变异组成Variance components 0.0017 0.0021 0.3720 0.3717 变异比例Percentage of variation 0.4600 –0.5500 100.0900 Group1(东兴Dongxing) Group2(吴川Wuchua闸坡Zhapo) Group3 (东海East China Sea) 自由度df 2.0000 2.0000 77.0000 81.0000 平方和Sum of squares 0.6640 0.7640 28.6450 30.0730 变异组成Variance components –0.0022 0.0007 0.3720 0.3704 变异比例Percentage of variation –0.6000 0.1800 100.4100 SAMOVA Group1(洞头Dongtou) Group2(其余4个群体Other 4 populations) 自由度df 1.0000 3.0000 77.0000 81.0000 平方和Sum of squares 0.5430 0.8850 28.6450 30.0730 变异组成Variance components 0.0109 –0.0048 0.3720 0.3781 变异比例Percentage of variation 2.8700 –1.2600 98.3800

表 2 前鳞种群结构的分子方差分析及SAMOVA分析 Tab.2 AMOVA and SAMOVA analysis of Liza affinis

2.2 前鳞的种群历史动态

中性检验的分析结果(表 3)如下：Tajima's D方面，洞头、平潭、广东、东兴4个群体均为显著性负值，其余为非显著性负值。Fu's Fs方面，洞头、广东、东兴3个群体均呈显著性负值，总体也是显著性负值(F_S=-20.3900，P=0)。Ramos-Onsins & Rozas's R²检验的分析显示(表 3)：R₂的范围是0.1610~0.1620 (P < 0.05)，SSD值(0.0010~0.0220)和R_g值(0.0500~ 0.7200)均较小，且多数呈不显著性(P > 0.05)，表明中国前鳞群体没有显著偏离扩张模型。此外，核苷酸错配图(图 2)呈现明显单峰、单倍型网络图(图 1)的星状结构与Fu's Fs的显著性负值均表明中国前鳞在近期经历过快速扩张，其扩张时间参数τ为0.6370，计算出中国前鳞种群扩张大约发生在13.4199~1.4911万年前，处于晚更新世时期(刘嘉麒等, 1997)。

表 3(Tab.3)

表 3 前鳞的中性检验、核苷酸不配对分析 Tab.3 Neutrality tests and mismatch analysis of Liza affinis 分析范围 Range analysis 群体所处的地理位置 Region of population Tajima's D值(P值) Fu's Fs值(P值) R2值(P值) SSD值(P值) Rg值(P值) 东海 East China Sea 岱山Daishan –1.4000(0.0880) 0.5860(0.4240) 0.1620(0.0000) 0.0220(0.0180) 0.7200(0.7200) 洞头Dongtou –1.8670(0.0090) –4.7410(0.0010) 0.1620(0.0000) 0.0070(0.6000) 0.0830(0.5900) 平潭Pingtan –1.8480(0.0110) –6.5790(0.0000) 0.1620(0.0000) 0.0020(0.7500) 0.0500(0.9000) 南海 South China Sea 吴川, 闸坡 Wuchua, Zhapo –1.5600(0.0400) –0.4590(0.1580) 0.1620(0.0000) 0.0040(0.7000) 0.2220(0.7300) 东兴Dongxing –1.7070(0.0220) –3.4830(0.0000) 0.1610(0.0000) 0.0030(0.3200) 0.1940(0.3800) 总体Total –2.3250(0.0000) –20.3900(0.0000) 0.1620(0.0000) 0.0010(0.8600) 0.0840(0.8800)

表 3 前鳞的中性检验、核苷酸不配对分析 Tab.3 Neutrality tests and mismatch analysis of Liza affinis

图 1(Fig. 1)

图 2 前鳞核苷酸错配分析图 Fig.2 Nucleotide mismatch distribution of Liza affinis

3.1 中国前鳞的遗传多样性

某物种的遗传多样性丰富与否和该物种所反映的单倍型多样性(H_d)及核苷酸多样性(π)密切相关(Bonin et al, 2010)。该研究中，中国前鳞的遗传多样性(H_d=0.4840±0.0700，π=0.0010±0.0002)呈现出以洞头(H_d=0.8080±0.1130，π=0.0021±0.0006)为中心向南北逐渐减少的趋势，与许则滩(2015)报道的浙江前鳞遗传多样性由南向北递减的趋势相符。但该研究中洞头群体和岱山群体的遗传多样性有别于许则滩(2015)的研究结果(温州：H_d=0.6170±0.1350，π= 0.0014±0.0004；岱山：H_d=0.3950±0.1280，π=0.0007± 0.0003)，可能是由于二者分析的序列长度及群体中个体数量的不同而致[许则滩(2015)研究所采用的序列长度为623 bp，群体数量为温州16，岱山23]。与其他同样分布在日本南部及中国台湾海域等的鱼类基于COⅠ序列的遗传多样性相比，中国前鳞的遗传多样性处于较低水平：鱼(Miichthys miiuy) (H_d= 0.7310±0.0470, π=0.0028±0.0019) (Xu et al, 2014)，黄姑鱼(Nibea albiflora) (H_d=0.6970±0.0023, π=0.0017± 0.0002)(Chang et al, 2016)，银鲳(Pampus argenteus) (H_d=0.8010, π=0.0040) (Xu et al, 2012)。此外，东海前鳞群体的遗传多样性(H_d=0.5790±0.0890, π=0.0014±0.0003)明显高于南海群体(H_d=0.3540±0.1020, π= 0.0006±0.0002)，其原因可能有：1)东海为前鳞的主要产区，有更为丰富的前鳞资源。曹艳等(2015)的研究发现，蓝点马鲛(Scomberomorus niphonius)在东海中的遗传多样性较高，推测其与东海为蓝点马鲛的主要产区有关。该研究结果与曹艳(2015)的研究情况相似，故推测东海前鳞群体的较高遗传多样性与东海为前鳞的主要产区有关。2)南海种群属于边缘种群。把中国沿海由南向北分布的各个前鳞种群作为一个整体，东海种群处于中心位置，而南海种群属于边缘区。边缘种群易受瓶颈效应、近交、遗传漂变的影响而导致遗传多样性呈下降趋势(Lesica et al, 1995)，该研究情况与杨喜书等(2017)对华南攀鲈(Anabas testudineus)群体的研究情况相似。3)南海的鱼类资源衰退。受过度捕捞及环境污染的影响，南海的鱼类资源处于被过度开发的状态(鞠海龙, 2012)，该情况可能是导致南海前鳞群体遗传多样性较低的原因之一。

3.2 前鳞鱼类的遗传结构和历史动态

该研究中AMOVA分析所得的组间变异比例(–0.6%~0.46%)与许则滩(2015)对浙江沿岸前鳞的研究结果(0.25%)相似，均属极低水平。结合单倍型网络图、群体间F_st值及SAMOVA的分析结果，可知中国前鳞群体在各个地区、海域、海峡间无明显的遗传分化，推测该现象的形成可能与以下因素有关：1)频繁的基因交流。前鳞属洄游性鱼类，可在湖泊和江河索饵场洄游到海中产卵场进行繁殖(陆忠康, 1985)，所产鱼卵主要漂浮在中国沿岸浅海处，随洋流漂流。已知冬季黑潮会把海洋有机物从菲律宾经东海送至中国大陆沿海，再结合中国南海的夏季暖流流向浙江近海，且目前暂无关于东海及南海的前鳞群体间在繁殖习性上存在差异的报道。故中国前鳞可能受到洄游产卵和海洋洋流的共同影响(Nakajima et al, 2014)，导致分布在不同地理位置的群体间有着频繁的基因交流，没有足够时间来积累遗传变异。2)近期种群扩张事件。中国前鳞的单倍型网络图、Tajima's D和Fu's Fs中性检验、核苷酸错配分析等结果都表明中国前鳞整体在近期经历过种群的快速扩张，估算出其扩张大约发生在13.4199~1.4911万年前，处于晚更新世时期(刘嘉麒等, 1997)。这一时期，冰期中海平面的升降造成大部分近海鱼类都发生过种群的快速扩张(Ni et al, 2014)，推测中国前鳞的快速扩张可能也与海平面的变化有关。在距今13.4199~1.4911万年的时间跨度中，北半球经历了2次冰期及1次间冰期(Stirling et al, 1995; Winograd et al, 1997)，其交替过程如下：距今约13.4199~13.2000万年前，北半球处于冰期，而在13.2000~11.5000万年前，北半球处于末次间冰期，至11.5000~1.5000万年前，北半球处于末次冰期。冰期中的气温很低，中国大陆出现了大量冰原，中国沿海的海平面出现明显的下降(下降幅度最大可达150 m)，中国沿海大陆骨架大面积露出(Wang, 1999)，大大压缩了前鳞的生活空间。而末次间冰期是近15万年来全球最暖的时期，大陆的冰川开始消退，海平面逐渐上升，使前鳞的生活空间得以扩大，故前鳞较强的鱼卵扩散能力及活动能力伴随海侵导致种群发生快速扩张。资料显示，冰期中海平面下降时，孑遗的物种会集中于避难所中形成随机交配的祖先群体，故在冰期后经种群扩张至不同地理位置的群体间的遗传分化不明显(李双, 2015)。该研究中中国前鳞的扩张时间与之吻合，满足共同祖先假说。

3.3 中国前鳞的种质资源保护

鱼类作为水生生态系统的高级消费者，处于食物链的较高位置，是生态系统中物质循环、能量流动的重要参与者，与水生环境有着密切的联系，故各种水体污染和不适当的人类活动都可能对鱼类的生存产生不利的影响。而正常的生态系统功能的维持，需要制定相关的种质资源保护方案。仅就该研究中，中国前鳞没有形成明显的地理聚群，且各地理群体间的分化程度较低，可以将其视作一个进化显著单位。但各地区群体的遗传多样性差异较大，其中，洞头群体的遗传多样性(H_d=0.8080±0.1130, π=0.0021±0.0006)明显高于其他地理群体，建议优先对其进行保护。由于采样限制，该研究只涉及中国前鳞部分群体，且线粒体COⅠ基因的进化速度不及线粒体控制区序列，不一定能如后者那样反映中国前鳞群体间可能存在的细微分化。另外，线粒体仅为母系遗传，也可能没有体现双亲遗传的核基因所能反映的遗传信息。因而在日后的工作中，需要结合测序控制区序列和基于核基因的多样性分析技术，补充更多地理群体，以全面了解中国前鳞的遗传背景及其基因流通道的分布，为种质资源的保护和开发利用提供科学依据。