《Biological Invasions》：Origin of the invasive sea anemone Metridium senile in the southern cone of South America
编辑推荐：
摘要：橙红列指海葵(Metridium senile)(Linnaeus, 1761)于2005年首次报道于南美巴塔哥尼亚中部，当时已呈高密度分布，随后沿智利峡湾区域广泛扩散。其爆发性增殖产生了显著的生态影响，尤其表现为排挤本地底栖生物群落及覆裹具有重要商业价
摘要：橙红列指海葵(Metridium senile)(Linnaeus, 1761)于2005年首次报道于南美巴塔哥尼亚中部，当时已呈高密度分布，随后沿智利峡湾区域广泛扩散。其爆发性增殖产生了显著的生态影响，尤其表现为排挤本地底栖生物群落及覆裹具有重要商业价值的海胆床。尽管其扩散迅速，M. senile在南半球的引入途径与生物地理起源仍不明确。本研究通过分析智利巴塔哥尼亚与阿根廷海岸M. senile种群的核基因(ITS)与线粒体基因(COI, COIII)序列数据，并结合使用六放珊瑚-v2-海葵目标捕获探针集(hexacoral-v2-actiniaria baitset)获取的全基因组SNP位点，探讨其起源与扩散历史。系统发育与群体结构分析识别出三个得到强支持的分支，分别对应主要的太平洋与大西洋谱系，且其中一个分支内存在显著分化，东北大西洋种群与西北大西洋种群彼此分离。智利样本主要聚类于北太平洋类群，而阿根廷样本与北大西洋种群亲缘关系密切，表明存在独立的引入途径。近似贝叶斯计算(Approximate Bayesian Computation, ABC)支持阿根廷种群源于东北大西洋，智利水域的入侵源自北太平洋(可能为单次引入)，推测经船舶压舱水(ballast water)运输，并以逐步扩散(stepping stone)模式向南蔓延至整个区域。研究结果证明南美南部锥体地区M. senile存在多次区域性独立的引入事件，最可能由海上航运介导。这些发现凸显了该物种较高的入侵潜力，强调需针对性监测以遏制其在南半球进一步扩散。
《Biological Invasions》刊载论文解读：南美南部锥体地区入侵海葵Metridium senile的起源与引入路径研究
一、研究背景与立项依据
橙红列指海葵Metridium senile(Linnaeus, 1761)(刺胞动物门Cnidaria：海葵目Actiniaria)为环北方分布种，原生分布于北半球，已被认定为南半球多处海域的入侵种，南非、福克兰群岛(马尔维纳斯群岛Islas Malvinas)及南美沿岸均有定殖记录。该物种2005年首现于巴塔哥尼亚中部并已形成高丰度种群，此后沿智利峡湾(fjord region)持续向南扩展，通过覆裹和排挤本地底栖群落特别是经济重要种海胆(sea urchin)，对北巴塔哥尼亚手工渔业造成经济影响。然而M. senile在南美南锥体(Southern Cone of South America)的引入途径与生物地理起源一直悬而未决。传统线粒体标记(COI、COIII)与核ITS因海葵类(Anthozoa)线粒体进化速率极低且M. senile盛行无性足盘裂解(pedal laceration)致克隆集聚、遗传多样性低，难以区分不同地理来源种群。已有全球六放珊瑚目标捕获测序(sequence-capture)研究建立了M. senile北半球各源区(日本、北美西岸、北美东岸、欧洲)系统发育框架，为本研究比对南美样本提供参照。鉴于此，研究人员采集智利巴塔哥尼亚太平洋沿岸8个位点及阿根廷大西洋沿岸1个位点的M. senile标本，联用传统基因测序与全基因组目标捕获SNP分析，结合系统发育重建、群体结构解析及近似贝叶斯计算(Approximate Bayesian Computation, ABC)模拟，旨在确认南美南部锥体M. senile身份、溯源其定殖祖先区域并重构引入扩散路线。
二、主要关键技术方法
研究人员于智利巴塔哥尼亚太平洋岸(Carelmapu、Guaitecas、Messier Channel、Punta Arenas等)8个位点及阿根廷大西洋岸Puerto Madryn采集M. senile组织标本，部分样本提取gDNA后进行ITS、COI(Folmer引物LCO1490/HCO2198)、COIII常规PCR与Sanger测序比对确认物种；另选40个样本(智利32份、阿根廷8份)采用六放珊瑚-v2-海葵探针集(hexacoral-v2-actiniaria baitset)进行目标序列富集(target enrichment)，Illumina NovaSeq PE150测序；同时整合已发表全球M. senile与近缘种拼接叠连群(contigs, NCBI PRJNA666413)。生物信息学流程包括SPAdes de novo组装、PHYLUCE匹配探针提取位点、MAFFT比对与Gblocks修剪，构建50%与75%分类单元占有率(taxon occupancy)两位点矩阵；最大似然(Maximum Likelihood, ML)系统树用RAxML-NG与IQ-TREE 2(GTR+FO+G4m模型)重建；SNP调用采用BWA比对至参考叠连群、Picard去重复、GATK3分型、VCFtools过滤为双等位、最小位间距800 bp、75%完整性阈值的非连锁SNP集，群体结构用STRUCTURE(K=1–12)分析并由Structure Harvester按Evanno法确定最优K；ABC模型选择以msprime模拟各入侵场景(单/多次引入、逐步扩散stepping-stone或独立引入)，用R abc包进行多项式Logistic回归(1%容忍度)及交叉验证遴选最佳入侵历史模型。
三、研究结果
DNA sequences（DNA序列分析）
研究人员从南美样本获得ITS(596 bp)53条、COIII(634 bp)49条、COI(553 bp)53条。系统树显示南美太平洋与大西洋样本均与北半球M. senile聚在一起，证实南锥体太平洋与大西洋沿岸定殖种为M. senile。但传统标记在南美样本内及与北半球间几乎无遗传分化，无法溯源，仅完成物种鉴定。
Sequence-capture data——Assembly, locus recovery, and alignment statistics（目标捕获数据——组装、位点回收与比对统计）
40个新样本测序原始reads均值约473万条，去低质量及去重复后均留存>98%。与已发表contigs合并，用hexacoral-v2海葵探针子集匹配，50%占有率矩阵得1,237个位点(309,781 bp)、75%占有率矩阵得462个位点(124,898 bp)，平均各位点长度827±232 bp，满足后续系统发育与SNP分析需求。
Phylogenetic Analysis（系统发育分析）
基于50%矩阵ML树将全球M. senile及近缘种分为三个高支持分支：Clade I为北太平洋东北冷温带(阿拉斯加、加州、华盛顿、BC省)；Clade II含西北太平洋(日本)、部分北太平洋东北种群及南美麦哲伦类群(智利、福克兰群岛太平洋侧集群)——本研智利样本(Carelmapu CAR、Guaitecas GUA、Punta Arenas PTA)归入此Clade II；Clade III细分IIIa(日本+阿拉斯加部分)、IIIb(北欧海域丹麦、瑞典、挪威、冰岛、英伦三岛、荷兰)、IIIc(西北大西洋加拿大与新英格兰+部分北欧德/挪威+南美麦哲伦大西洋侧阿根廷ARG及福克兰群岛Stanley港集群)——阿根廷样本归入Clade IIIc，支持智利与阿根廷种群分属太平洋与大西洋不同源谱。
SNP（SNP群体结构分析）
筛选960个非连锁SNP做STRUCTURE分析，Evanno法ΔK峰值K=2区分M. farcimen(Cluster 4)与M. senile(Clade II+III)，次峰K=5展现 finer划分：Cluster 1(深蓝，西北大西洋及南大西洋福克兰Stanley港、阿根廷部分)对应Clade IIIc，Cluster 2(粉红，北太平洋东岸主导，智利样本及福克兰偏远太平洋侧集群)对应Clade II，Cluster 3(浅蓝，东北大西洋)对应Clade IIIb，Cluster 4(灰，M. farcimen)对应Clade I，Cluster 5(红，北太平洋日/阿留申部分)对应Clade IIIa。智利样本主属Cluster 2(北太平洋来源)，少数具Cluster 5成分；阿根廷样本呈Cluster 1(~40%)与Cluster 3(~60%)混合特征，指示北大西洋(东北与西北)混合或混合源区定殖；福克兰Stanley港归Cluster 1(西北大西洋来源)，偏远区归Cluster 2(北太平洋来源)，证实福克兰存在两独立入侵来源。
Approximate bayesian computation——Pacific connectivity（近似贝叶斯计算——太平洋连通性）
检验北太平洋入迁智利后逐步向南扩散(Stepping-Stone, SS)、智利链逐步扩散+福克兰太平洋侧独立自北太平洋(2I)、逐步扩散+PTA与福克兰各自独立(3I)、仅北部智利逐步扩散余三者独立(4I)四模型。ABC模型选择最高后验概率归2I模型(p=0.500)：北太平洋单次引入智利北部→逐步向南经Carelmapu→Guaitecas→Messier Channel→Punta Arenas，福克兰偏远太平洋侧集群为北太平洋另一独立引入事件。
Approximate bayesian computation——Atlantic connectivity（近似贝叶斯计算——大西洋连通性）
排除西北大西洋ARG×东北大西洋福克兰组合(零接受模拟)，剩余五模型ABC优选：阿根廷(ARG)源于东北大西洋(Northeastern Atlantic)、福克兰Stanley港(FIS)源于西北大西洋(Northwestern Atlantic)之双独立引入模型具最高支持(p≈1.000，原始拒绝比例亦最高25%)。FST值佐证：FIS-西北大西洋=0.028(极近缘)、FIS-东北大西洋=0.235(远缘)；ARG-东北大西洋=0.118、ARG-西北大西洋=0.159，略倾向东北大西洋为主要源区，ABC量化支持东北大西洋为ARG主源。
四、讨论与结论翻译总结
传统COI/COIII/ITS标记确认南美太平洋与大西洋沿岸物种为M. senile，但因海葵目(Anthozoa)线粒体低突变率及无性繁殖致克隆斑块化低遗传多态性，此类PCR靶标标记无法分辨种群来源；全基因组SNP方揭示清晰地理渊源。智利巴塔哥尼亚种群遗传组分主要对应北太平洋集群(Cluster 2及部分Cluster 5)，ABC支持北太平洋单次引入智利北部后沿峡湾逐步向南扩散，福克兰群岛偏远太平洋侧种群为北太平洋又一独立引入——推断引入媒介为往来北太平洋港口船舶压舱水(ballast water transport)。阿根廷种群呈现西北与东北大西洋混合信号(Cluster 1+Cluster 3)，ABC最支持东北大西洋为主要定殖源，暗示阿根廷可能接受多源或混合源港口引入；福克兰Stanley港种群独立源自西北大西洋而非自阿根廷逐步扩散。福克兰存在太平洋源(北太平洋独立引入)与大西洋源(西北大西洋独立引入)两波互不关联入侵事件。研究证实南美南部锥体M. senile经海运导致至少三次区域性独立引入(智利太平洋侧北太平洋源逐步扩散、阿根廷大西洋侧东北大西洋源、福克兰两独立源)，凸显其强入侵潜能。传统线粒体及核ITS标记在此类入侵历史推断中力有未逮，目标捕获SNP方能有效解析。尽管入侵种群遗传多样性低，M. senile成功定殖扩布符合入侵遗传学"遗传悖论(genetic paradox of invasions)"，可能受表型可塑性(phenotypic plasticity)、引入后适应性演化、多源混种(admixture after multiple introductions)及表观遗传变化共同促成。结果强调针对航运压舱水管控与南半球沿海定向监测之必要性。
