中国鲚属鱼类物种有效性及刀鲚种群遗传结构研究

摘要

分析了9个群体214尾刀鲚Coilianasus、湖鲚C.nasustaihuensis、短颌鲚C.brachygnathus、七丝鲚C.grayii及凤鲚C.mystus个体的COI基因DNA条形码序列变异。结果显示，密码子碱基组成存在明显的偏倚性，GC含量的高低顺序为1st>2nd>3rd，群体间密码子碱基组成差异不明显。所有个体共包含71种单倍型，单倍型共享现象在刀鲚、湖鲚、短颌鲚群体问频繁出现。单突变位点主要集中在100bp和600bp附近。刀鲚、短颌鲚和湖鲚群体间的遗传距离在0.253％-0.557％之间，这与群体内的遗传差异相近，却显著低于COⅠ基因DNA条形码鉴别不同物种2％的遗传距离阈值，也不符合种间差异大于种内差异10倍的物种鉴定规则，因此表明这3个群体应为同一物种。但是凤鲚两群体间的遗传距离为5.080％，远大于2％的鉴别阈值，显示凤鲚两群体出现明显的分化，其程度至少已达亚种级差异水平。以日本鲲Engraulisjaponicus为外群，用邻接法(NJ)、最大简约法(MP)和最大似然法(ML)构建了分子系统树显示，刀鲚、湖鲚和短颌鲚群体聚在一起，未能各自形成单系；凤鲚根据地理分布聚为两支；七丝鲚则聚成单系。表明COⅠ基因条形码技术可用于我国鲚属物种的鉴定。在确定刀鲚、湖鲚和短颌鲚为同一物种之后，进一步对刀鲚6个群体(舟山ZS、靖江JJ、太湖TH、九段沙JD、鄱阳湖PYH、有明海AR)的遗传结构和演化历史进行了研究。结果显示：种群分化指数Fst和基因交流值Nm分别在0.17％-48.1％和0.54-288范围内，鄱阳湖群体与其余群体问的Fst值显著大于，Nm值显著小于其余群体之间的数值。K-2-P遗传距离、AMOVA分析结果以及基于COⅠ基因DNA条形码序列构建的NJ树的结果均表明，长江中游的鄱阳湖群体与下游群体发生较为明显的分化，而长江中下游及日本有明海刀鲚群体间未发生明显的地理分化。中性检验和网络亲缘关系分析皆表明，刀鲚群体整体经历过种群扩张的历史，扩张时间约在更新世末期的0.028-0.038百万年前，这可能与末次盛冰期时海平面升降有关。
　　 采用PIMA法(PCRisolationofmicrosatellitearrays)分离多态性较高的刀鲚微卫星标记，结果得到102个含有5次以上重复的微卫星的序列，设计并合成了21对微卫星引物，经严格筛选得到11对多态性高且稳定的微卫星引物，利用这11对引物对6个刀鲚群体(舟山ZS、靖江JJ、太湖TH、九段沙JD、鄱阳湖PYH、有明海AR)的遗传结构和种群历史进行分析。11对微卫星引物的平均多态信息指数PIC为0.7485，每对引物扩增出的等位基因数5-14个，平均8.45个，Micro-checker检测未发现无意等位基因(NullAllele)。6个群体中，PYH、AR群体的平均等位基因数和PIC值(分别为5.273、0.5983和5.455、0.6082)显著低于其余4个群体的平均值(7.128、0.7177)，表现出了较低的遗传多样性。11个微卫星位点中，除了micro639位点外，其余位点均在不同群体中出现偏离Hardy-Weinberg平衡的现象，舟山ZS、靖江JJ、太湖TH、九段沙JD、鄱阳湖PYH、有明海AR6个群体偏离Hardy-Weinberg平衡的位点个数分别为4、4、4、3、4、1，AR群体偏离平衡的位点数显著小于其他群体。群体间分化指数Fst和基因交流值Nm分别在0.39％-13.40/0和1.5223-63.5255范围内，其中PYH、AR群体与其他群体之间的Fst平均值(12.7％)远大于ZS、JD、TH、JJ4个群体之间Fst平均值(0.97％)；而Nm均值(1.72)则远小于4个群体间的Nm均值(29.98)。另外，AMOVA分析和基于Bayesian模型的聚类结果均表明，PYH群体、AR群体与长江下游的其他群体发生了一定程度的分化，而且PYH群体的分化趋势更明显，但是长江下游的4个群体间交流频繁，未出现分化迹象。Bottleneck分析得出刀鲚群体整体曾经历过强烈的瓶颈效应，而其中的AR群体则几乎未受其影响。我们推测现今日本有明海附近的刀鲚可能是由原本分布在长江口及外海广大区域的刀鲚群体在更新世末期海平面下降时受困于日本九州岛西南侧深水区域，后经历间冰期海平面上升和对马暖流的阻隔从而逐渐与长江口群体失去交流而形成。