全球广布红藻 Asparagopsis taxiformis Lineage-2（Florideophyceae, Rhodophyta）低遗传多样性的系统地理学研究

《Algal Research》：Low genetic diversity in the cosmopolitan Asparagopsis taxiformis Lineage-2 (Florideophyceae, Rhodophyta)

　　

世界广布种 Asparagopsis taxiformis Lineage-2（Florideophyceae, Rhodophyta）表现出低的遗传多样性。研究人员分析了来自韩国的样本，发现所有韩国样本均属于 Lineage-2。全球单倍型分析显示，地中海种群，尤其是意大利，拥有最多的私有单倍型。基于线粒体细胞色素c氧化酶亚基I（cox1）的分析解析出东亚地区三个优势单倍型并呈区域结构。尽管 Lineage-2 在不同区域间表现出低至中等水平的细胞器基因组分化，但其整体遗传多样性较低。这些结果提示，意大利东部和朝鲜半岛南部可能曾在冰期充当避难所。

本文解读《世界广布种 Asparagopsis taxiformis Lineage-2 低遗传多样性的系统地理学研究》。

研究背景与意义：红藻属 Asparagopsis 因具有海洋生态入侵潜力和作为家畜甲烷减排饲料添加剂的功效而备受关注。该属历史上包含8个种，后经分类修订归并为3个种，其中 A. taxiformis 和 A. armata 具有异型二倍单倍体生活史。以往基于线粒体 cox2–3 间隔区（cox2–3 spacer）和部分 cox1 基因的系统发育分析显示，A. taxiformis 包含6个单系谱系（L1–L6），但形态差异不明显，被视为隐存种复合体。隐存种复合体会掩盖真实的分布格局、原产/引入状态及谱系特异的生态影响，因此需要整合基因组尺度数据开展系统地理学评估。本研究旨在利用韩国群体样本，结合全球已有数据，解析韩国 A. taxiformis 的谱系组成和遗传多样性，并深入刻画广布谱系 L2 的全球系统地理格局。

研究人员从韩国16个地点采集39个 A. taxiformis 个体，结合全球404个个体（含365个已鉴定谱系的样本）进行分析。通过 cox2–3 spacer 和 cox1 扩增，确定所有韩国样本均属于 L2。在此基础上，组装了韩国、意大利、美国三个 L2 株系的完整线粒体和质体基因组。结果显示，L2 虽广布多个洋盆，但细胞器遗传多样性极低；地中海区域特别是意大利拥有最高私有单倍型；cox1 在东亚解析出 H1、H2、H3 三个优势单倍型并呈区域结构；比较基因组分析显示 L2 的核苷酸多样性在跨洋盆广布种中最低。中性检验和贝叶斯天际线分析提示，意大利东部和朝鲜半岛南部可能为冰期避难所。该研究为 A. taxiformis 隐存种复合体的分类修订和谱系特异的入侵及甲烷减排能力评估提供了重要基线。

关键方法：研究使用了韩国16个地点39个野外样本、意大利那不勒斯参考培养株（CCAP 1341/1）以及美国加州公开数据（SRR12672538）。方法包括：Sanger测序获取线粒体 cox2–3 spacer 和 cox1 序列；Illumina Novaseq6000 平台测序，用 GetOrganelle 组装线粒体和质体基因组；MAFFT 比对与 IQ-TREE 构建最大似然树；TCS/PopART 构建单倍型网络；DnaSP 计算多样性指数；Tajima's D 和 Fu's FS 检验种群扩张；BEAST 进行贝叶斯天际线分析；利用 R 比较红藻纲物种的核苷酸多样性并做 Wilcoxon 秩和检验。

研究结果：

3.1 线粒体谱系的全球分布与系统地理结构：基于404个个体，共恢复 L1–L6 六个谱系。L2、L3、L4 跨多个洋盆广布，L1、L5、L6 局限于夏威夷、法属波利尼西亚和澳大利亚沿岸。所有韩国样本均归入 L2，显示 L2 是最广布的谱系。

3.2 系统发育分析：从37条 cox2–3 spacer 和32条 cox1 序列构建的最大似然树中，韩国序列均靠近 L2，但 L2 不呈单系；基于连接的 cox2–3 spacer 和 cox1 数据集（24个唯一单倍型）恢复 L1、L2、L3、L4、L6 五个谱系，其中 L1/L2/L6 形成强支持分支，L2 恢复为单系（UFBoot=61%）并与 L1 构成姐妹支。韩国序列在 cox1 分析中分为 H1（18个个体）、H2（11个个体）、H3（1个个体）三个单倍型。

3.3 基于 cox2–3 spacer 的遗传多样性与系统地理学：TCS 单倍型网络（n=402）将 L1–L6 分隔2–9个突变步。L2 呈星状拓扑，中心单倍型 H1 广布，私有单倍型在意大利最多（10个），韩国仅1个私有单倍型。中性检验显示整体（D=-2.350，FS=-26.032）和地中海群体（D=-1.927，FS=-10.549）显著负值；贝叶斯天际线分析显示地中海群体约0.175 Ma以来增长，东亚相对稳定。

3.4 基于 cox1 的遗传多样性与系统地理学：TCS 网络（n=121）显示 L2 内呈较网状结构，三个主要单倍型 H1/H2/H3 构成框架，H1 在韩国20个、H2 11个、H3 1个。东亚三个单倍型具有地理分隔：H1 集中于东海和济州东部，H2 在韩国南部和西部海岸及济州，H3 主要分布于琉球群岛。地中海仅含 H1/H2，夏威夷含 H1/H3。中性检验显示整体 FS=-5.176、东亚 FS=-2.501、地中海 FS=-2.979 均显著负值；贝叶斯天际线分析显示整体和东亚约0.03–0.02 Ma以来持续增长。

3.5 线粒体和质体基因组分析：三个 L2 个体（韩国、意大利、美国）均组装出完整环状细胞器基因组。线粒体基因组大小为26,073–26,174 bp，含52个基因（25个蛋白编码基因、3个 rRNA、24个 tRNA）；质体基因组大小为177,083–177,095 bp，含238个编码基因（202个蛋白编码基因、3个 rRNA、1个 tmRNA、1个 ncRNA、31个 tRNA）。三者基因内容和基因顺序一致，仅线粒体短重复拷贝数存在差异。

3.6 细胞器基因组比较变异：L2 内配对 SNP 极小，韩国-意大利线粒体16个 SNP（0.06%）、质体49个 SNP（0.03%）；韩国-美国线粒体119个 SNP（0.46%）、质体396个 SNP（0.22%）。在红藻纲（Rhodymeniophycidae）比较中，L2 的线粒体和质体核苷酸多样性处于低至中等水平，在跨三个洋盆的广布种中最低；与多数物种差异显著，但与 Gracilaria vermiculophylla 和 Gracilaria salicornia 的线粒体核苷酸多样性无显著差异。

讨论与结论：讨论部分指出，L2 虽广布但细胞器遗传多样性低，结合显著负的 Fu's FS 值，支持近期种群扩张事件。地中海区域高私有单倍型富集和贝叶斯天际线分析支持意大利作为长期避难所，可能与墨西拿盐度危机有关，但不能排除经苏伊士运河引入的观点，需要更多印度-太平洋和红海样本检验。东亚 cox1 结构显示三个优势单倍型的地域分化，可能与末次盛冰期低海平面和洋流格局相关，南部朝鲜半岛附近可能为冰期避难所，冰后期随洋流扩张。研究人员假设 L2 经历长期进化历史后，发生反复的区域瓶颈和冰后期扩张。最终结论可翻译为：本研究确认韩国所有 A. taxiformis 属于 L2，并揭示 L2 在全球尺度上具有异常低的细胞器遗传多样性；意大利东部和朝鲜半岛南部可能为 L2 的冰期避难所；未来需要更广泛采样和基因组数据来验证历史动态假说。