但后来的研究又让情况变得更加复杂。通过对东海、朝鲜海峡、黄海和鄂霍次克海沙鳗的 mtDNA 序列和微卫星(ms)比较,发现这两个物种之间存在单向基因渗入现象。mtDNA 和 ms 结果显示,黄海和朝鲜海峡的沙鳗为 A. japonicus,鄂霍次克海的为 A. heian,而东海沙鳗用 mtDNA 鉴定为 A. heian,用 ms 鉴定却为杂交种。这一发现使得东海沙鳗的分类地位变得模糊不清,因此,迫切需要从基因组层面开展系统的种系地理学研究,以深入了解沙鳗物种形成的演化过程。
为了解开这些谜团,釜庆国立大学(Pukyong National University)的研究人员 Jin-Koo Kim 和 Yu-Jin Lee 在《Scientific Reports》期刊上发表了题为 “Genome-wide SNP discovery and phylogeography of sandlances (Ammodytes japonicus and Ammodytes heian) in the Northwest pacific by MIG-Seq with evolutionary implications” 的论文。他们通过研究得出结论:基于单核苷酸多态性(SNP)的遗传特征与 mt 结果相似,但与 ms 结果存在一定差异,表明东海(GI 和 JU 地区)的沙鳗很可能是 A. heian,且杂交现象较为罕见。这一研究成果为沙鳗的分类和演化研究提供了重要依据。
在这项研究中,研究人员主要运用了以下几种关键技术方法:首先,采用 MIG-Seq(Multiplexed intersimple sequence repeat genotyping by sequencing)技术,这是一种基于聚合酶链反应(PCR)的高效且经济的基因组水平 SNP 识别方法;其次,利用高通量测序技术对构建的文库进行测序;最后,运用多种生物信息学分析方法,如使用 Stacks 软件进行 SNP calling 和相关信息生成,利用 Genodive、STRUCTURE 等软件进行遗传指数计算和种群结构分析,通过 RAxML 和 Phylogears2 进行系统发育分析,借助 R 语言进行主成分分析(PCA)和 Wilcoxon 秩和检验等 。
下面让我们详细看看研究人员都有哪些重要发现。
遗传多样性
研究人员从韩国的四个地区(东海的 Gijang 和 Jumunjin、朝鲜海峡的 Tongyeong、黄海的 Baengnyeongdo)以及日本的一个地区(鄂霍次克海的 Wakkanai)选取了 48 条沙鳗进行研究。经 mtDNA 序列分析确定它们分别属于 A. japonicus 和 A. heian。对这些样本进行测序后,共获得 99,175,917 条原始读数,去除低质量读数后,得到 74,195,010 条。经过分析,发现了 10,545 个区域特异性 SNP。通过这些 SNP 对不同地区 A. japonicus 和 A. heian 的遗传多样性进行比较,结果显示,多态性位点比例在 Tongyeong 最高(0.565),在日本的 Wakkanai 最低(0.489);观察杂合度(H<sub>0</sub>)在 Jumunjin 最高(0.159),在 Tongyeong 最低(0.129);核苷酸多样性(Pi)也是在 Jumunjin 最高(0.204),在 Wakkanai 最低(0.169);近交系数(F<sub>IS</sub>)则在 Baengnyeongdo 最高(0.172),在 Wakkanai 最低(0.079)。
研究人员分别使用 10,545 个 SNP 和 55,026 个 SNP 进行主成分分析(PCA)和系统发育分析。根据 PC 1 - 2 轴,A. japonicus 群体倾向于位于正 PC 1 轴和负 PC 2 轴,而 A. heian 群体倾向于位于负 PC 1 轴和正 PC 2 轴。在 PC 2 - 3 轴上,A. japonicus 群体倾向于位于两个负轴,A. heian 群体位于两个正轴。有趣的是,一些个体在 PC 1 - 2 轴和 PC 2 - 3 轴上的归属发生了变化,比如 24 号和 27 号个体从 PC 1 - 2 轴上的 A. japonicus 转移到了 PC 2 - 3 轴上的 A. heian。通过 Wilcoxon 秩和检验,在 PC 1、PC 2 和 PC 3 轴上都能明显区分这两个物种。PCA 结果与系统发育分析结果大致相符,系统发育树显示出两个明显的聚类,A. japonicus 群体由来自 Baengnyeongdo 和 Tongyeong 的个体组成,A. heian 群体由来自 Gijang、Jumunjin 和 Wakkanai 的个体组成,但 Baengnyeongdo 有一个个体(B38)属于 A. heian 群体,Jumunjin 有一个个体(J29)属于 A. japonicus 群体,这表明它们的归属存在一定不确定性。
综合研究结果,研究人员重新评估了沙鳗的进化历史。他们发现,MIG - Seq 方法得到的 SNP 结果比 ms 结果更接近 mtDNA 分析结果。例如,基于 mtDNA 被鉴定为 A. japonicus 的 Baengnyeongdo 和 Tongyeong 的沙鳗,除了个别个体外,基于 SNP 结果也被确认为 A. japonicus;而在杂交率较高的 Gijang 和 Jumunjin 地区,基于 mtDNA 和 ms 被鉴定为 A. heian 或杂交 / 回交的个体,大多通过 SNP 分析被确认为 A. heian 。不同遗传标记之间的差异可能是由于物种快速形成、不完全谱系分选、基因渗入、性别偏向扩散或选型交配等原因造成的。在沙鳗的研究中,ms 检测到的基因渗入可能是由于位点数量有限或无效等位基因导致的偏差,也可能是样本量或地理覆盖范围的问题。此外,沙鳗可能经历了复杂的进化历史,在冰川和间冰期,不同谱系之间可能发生了多次隔离和重新连接事件。
从分类学角度来看,SNP 研究结果更支持将东海沙鳗和黄海 + 朝鲜海峡沙鳗分别认定为 A. heian 和 A. japonicus。虽然传统上认为脊椎骨数量是区分沙鳗种群的有用形态特征,但最新研究表明,它受环境因素影响较大,可能不是区分两个物种的关键特征。而峡部黑色素细胞的形状、幼虫期头长与体长的比例以及产卵地的分离等特征,更能支持 A. heian 和 A. japonicus 是不同物种的观点。
这项研究意义重大,它为沙鳗的分类和进化研究提供了新的视角和重要依据。通过基因组层面的研究,我们对沙鳗的遗传多样性、种群结构有了更深入的了解,这有助于更好地保护和管理沙鳗资源。同时,研究中发现的不同遗传标记之间的差异以及复杂的进化历史,也为其他物种的研究提供了参考,让我们认识到在物种研究中综合考虑多种因素的重要性。未来,还需要对包括 A. hexapterus 和 Ammodytes personatus 等在内的密切相关物种进行更全面的基因组分析,以进一步完善我们对沙鳗属物种的认识。