深渊端足类Paralicella的分布悖论——P. tenuipes与P. caperesca物种复合体中广泛分布谱系缺乏清晰地理分区

《Zoologica Scripta》：The Paralicella Paradox: Broadly Distributed Lineages Without Clear Geographic Partitioning in Paralicella tenuipes and the Paralicella caperesca Species Complex

　　

深渊环境（3000–6000 m海洋深度）常被认为因其广阔空间范围、表观环境均一性和基因流障碍少而促进跨洋盆的广泛遗传连通性（genetic connectivity）。虽然部分深海分类群符合这一预期，但越来越多证据显示，即使在广泛分布物种中也可存在显著遗传结构。本研究利用扩展的全球深渊取样多基因座数据（multilocus data；COI，细胞色素c氧化酶亚基I，cytochrome c oxidase subunit I；16S rRNA，16S核糖体RNA；28S rRNA，28S核糖体RNA），重新评估两种广泛分布深渊端足类（Amphipoda）Paralicella tenuipes Chevreux, 1908和Paralicella caperesca Schulenberger & Barnard, 1976内的谱系结构。物种界定分析（species delimitation analyses）在P. tenuipes中恢复反复出现的线粒体谱系结构（mitochondrial lineage structure），在P. caperesca中恢复更强的多基因座分歧，与后者为物种复合体（species complex）的解释一致。TCS单倍型网络（haplotype networks）未显示一致的地理或深度相关分隔，恢复的谱系和主要单倍型分布于多个洋盆。人口统计和中性检验（demographic and neutrality tests）表明P. tenuipes具有中等至高的线粒体多样性，但无一致的全种近期人口变化信号；P. caperesca显示更高的线粒体多样性和更强的、依赖标记的人口统计信号。总之，这些结果表明遗传结构不与地理或深度强烈关联。取样的地理范围表明Paralicella是迄今分布最广的深海端足类属之一。这种广泛深渊分布与非地理关联谱系结构的组合表明，在这些端足类中产生遗传多样性的过程并不与常用于推断深渊生态系统多样化的地理单元整齐对应。

研究背景：深渊区（abyssal zone，约3000–6000 m）是地球表面最广阔、地形上较连续的底栖环境，常被假设因环境均一和物理障碍少而促进跨洋盆基因流。然而，越来越多分子证据显示，一些广泛分布深渊分类群仍存在明显遗传结构。Paralicella为深渊诱饵陷阱中常见且丰度高的端足类属，其中P. tenuipes Chevreux, 1908和P. caperesca Shulenberger and Barnard, 1976形态保守、全球调查常被鉴定，但分子研究在两个名义种内均发现分化谱系：P. tenuipes多被解释为广布种内地理广布线粒体谱系，P. caperesca则反复发现深分歧谱系，可能为物种复合体。此前研究受地理取样不足、标记覆盖差异、物种界定方法和谱系解释不同限制，谱系多样性、连通性和物种界限仍不确定。因此，研究人员扩展全球深渊取样和多基因座数据，重新评估两名义种的谱系结构、地理分区和连通性。

研究开展、结论与意义：研究人员组装P. tenuipes和P. caperesca扩展深渊数据集，纳入四个洋盆此前未取样区域的新材料，整合COI、16S rRNA和28S rRNA数据，应用三种物种界定方法于单基因和串联数据集，并结合系统地理与人口统计分析。结果表明，P. tenuipes存在反复但较浅的线粒体谱系结构，宜解释为广布种内四个地理广布线粒体谱系；P. caperesca具更深多基因座分歧，支持其为至少两个主要深分歧谱系组成的物种复合体。两名义种的谱系和主要单倍型跨多个洋盆分布，无一致地理或深度分隔。该研究确认Paralicella为迄今分布最广的深海端足类属之一，并提示深渊遗传多样性形成过程不简单对应常用地理单元，为后续基因组和整合分类研究提供基线。

关键技术方法（≤250字）：样品由RV Sonne、DSSV Pressure Drop、MV Pangaea和RV Dagon等航次的诱饵着陆器采集，覆盖太平洋、大西洋、印度洋和南大洋深渊区（北极未获Paralicella），新选168个个体（P. tenuipes n=119；P. caperesca n=49），并整合GenBank公开序列。提取DNA后扩增COI、16S rRNA和28S rRNA；用Geneious/MAFFT比对并去冗余代表序列；BEAST和IQ-TREE建树；GMYC（Generalized Mixed Yule-Coalescent，广义混合Yule-合并模型）、PTP（Poisson tree process，泊松树过程；含bPTP、mPTP）和ABGD（Automated Barcode Gap Discovery，自动条形码间隙发现）进行物种界定；PopART构建TCS单倍型网络，DNAsp计算多样性、中性检验和错配统计；MEGA计算未校正p距离（uncorrected p-distance）和K2P距离（Kimura 2-parameter）；GEBCO（General Bathymetric Chart of the Oceans，全球海洋通用测深图）深度数据估算潜在栖息地。

研究结果：
3.1 物种界定（Species Delimitation）：物种界定输出以分子操作分类单元（MOTU，Molecular Operational Taxonomic Unit）表示，反复恢复的MOTU再解释为分子谱系。P. tenuipes的实体估计随标记、方法和建树方式变化，单实体和四实体解释反复出现；仅COI+16S数据集的GMYC显著支持四实体，mPTP在该数据集BEAST树上推断四实体，ABGD对COI恢复四组、对16S恢复一组。总体不支持正式种级分离，但支持广布种内反复线粒体谱系结构。P. caperesca结构更强：单基因估计较高（GMYC COI为21、16S为10），含28S的串联数据集估计减少；COI+28S和COI+16S+28S的GMYC显著，mPTP在两种拓扑上均恢复两实体，ABGD对COI恢复三组、16S四组。多基因座分析支持P. caperesca物种复合体内两个主要深分歧分子谱系，单基因高估计提示额外未解决线粒体结构。

3.2 成对距离（Pairwise Distances）：P. tenuipes最大未校正p距离为COI 4.03%、16S rRNA 2.34%、28S rRNA 0.99%，K2P距离相近；谱系内COI p距离较低（0.32%–1.11%），谱系间较高（2.10%–3.91%），最大分歧在lineage 2与lineage 4之间，支持浅至中等谱系结构。P. caperesca最大p距离为COI 13.85%、16S rRNA 11.86%、28S rRNA 1.98%，K2P最高达15.55%、12.91%、2.01%；两个主要谱系的COI+16S+28S谱系内p距离为0.65%–0.79%，谱系间为6.37%–7.08%，支持物种复合体内深分歧谱系。

3.3 系统地理与人口统计分析（Phylogeographic and Demographic Analysis）：TCS网络显示两物种均有广布单倍型和分歧单倍型簇，但未按洋盆或深度清晰结构化。P. tenuipes的COI网络多个簇由推断突变步连接，大簇跨洋盆共享；COI+16S网络将四个谱系分开，但各谱系不限于单一洋盆或深度区间，支持地理广布线粒体谱系而非区域隔离。P. caperesca的COI网络更复杂，低频单倍型多且地理混合；COI+16S+28S网络显示两个主要组分由长突变路径分隔，亦不限于单洋盆或深度，支持两个主要深分歧谱系。人口统计显示，P. tenuipes线粒体单倍型多样性中至高（COI Hd=0.848；16S rRNA Hd=0.751），多数中性检验不显著，28S rRNA多样性低但Tajima's D显著负，无一致全种近期人口变化信号；lineage 3具显著Fu's F和R²。P. caperesca线粒体多样性很高（COI h=134、Hd=0.944；16S rRNA h=28、Hd=0.821），COI的Fu's F显著负，16S rRNA的Fu and Li's D*和F*显著负，28S rRNA无显著偏离，提示线粒体稀有变异过量及标记依赖的人口信号。

3.4 潜在栖息地估算（Potential Habitat Estimates）：利用810个个体的深度记录，P. tenuipes记录深度为2627–7502 m，P. caperesca为2090–6180 m，覆盖lower bathyal、abyssal和upper hadal。以GEBCO深度数据、ArcGIS和Mollweide等面积投影计算太平洋、大西洋、印度洋和南大洋中深度定义的潜在栖息地，两分类单元合计对应评估洋盆面积的82.3%（P. tenuipes潜在面积约284.0 million km²，P. caperesca约283.7 million km²）。该估算为深度代理的潜在栖息地面积，并非完整生态位模型或实际分布；北极和地中海无记录。

讨论总结：讨论指出，尽管扩展了地理范围、个体数和位点，仍未将遗传分组对应洋盆、水深或地理，说明当前线粒体和核糖体标记可识别宽泛遗传结构，但不足以可靠解决物种边界或地理分区。P. tenuipes为广布种内四个反复、弱至中等分歧且全球混合的线粒体谱系；P. caperesca为至少两个广泛分布、深分歧且缺乏清晰地理分隔的谱系组成的物种复合体。与Eurythenes gryllus sensu lato、Bathycallisoma schellenbergi等显示地理或深度结构不同，Paralicella谱系跨洋盆混合，与Abyssorchomene distinctus的广泛单倍型共享部分相似。其潜在栖息地超过评估四大洋面积的82%，并超过对Alicella gigantea的推断，但栖息地广度不预测遗传结构，形成“Paralicella paradox”。北极缺失可能涉及地形过滤、洋流、生态位占用或取样不足，尚不能定论。未来需ddRAD-seq、WGS-SNPs、UCEs等基因组尺度数据和整合形态学重新检验隐藏结构。

研究结论翻译：P. tenuipes与P. caperesca在全球深渊展示复杂且动态的谱系结构。跨标记和分析方法，物种界定、单倍型网络和人口统计结果一致恢复两个名义种内地理广布的遗传谱系，但这些谱系与洋盆、地理区域、深度或明显外部形态对应很少。两分类单元进化解释不同：P. tenuipes最符合广布种伴反复线粒体谱系结构；P. caperesca更适合视为由深分歧遗传谱系组成的物种复合体。结合广泛深渊深度范围及太平洋、大西洋、印度洋和南大洋的广泛潜在栖息地，这些持续模式确立P. tenuipes和P. caperesca物种复合体为目前深海已知分布最广的端足类分类单元之一。结果不支持该范围内简单遗传均一，而显示广布深渊清道夫可维持广泛单倍型共享并蕴藏分歧遗传谱系。解决产生和维持这种隐藏结构的过程，需要未来将广泛地理取样与基因组尺度数据和整合分类学方法结合，以更直接检验清道夫端足类的物种边界、连通性和隐藏谱系多样性。