整合系统发育基因组学揭示Dinaric喀斯特拟蝎的隐藏多样性与重复洞穴定殖

《Cladistics》：Integrative phylogenomics uncovers hidden diversity and repeated cave colonizations in Dinaric pseudoscorpions

　　

地下生态系统拥有由极端环境适应塑造的非凡生物多样性，但仍探索不足。洞穴生物局限于小而隔离的种群，易灭绝，因此整合物种界定对生物多样性评估和保护至关重要。Dinaric Karst是全球地下多样性热点，也是首个被描述的洞穴拟蝎Neobisium spelaeum的产地，其演化历史仍未知。研究人员通过整合线粒体和基因组数据与群体遗传学、物种界定、种群统计学建模和形态学分析，调查了N. spelaeum及其近缘种N. stygium的多样化动态。DNA条形码（DNA barcode）分析提示洞穴栖居Neobisium存在重复地下定殖。在该多系（polyphyletic）组合中，N. spelaeum和N. stygium形成一个支系（clade）。全基因组数据揭示八个深度分化谱系，对应受渐新世–中新世地质屏障塑造的喀斯特区域，另有七个单例（singletons），而线粒体-核不一致（mito-nuclear discordance）使演化轨迹重建复杂化。谱系间形态分化细微但具统计学显著性。种群统计学分析揭示270 000–470 000年前瓶颈后种群扩张，与中更新世（Middle Pleistocene）冰期旋回期间区域洞穴定殖事件相容。研究人员的发现突出Dinaric Karst作为地下特有性主要中心，强调保护喀斯特景观及其特有动物群的必要性，并为分类修订和保护规划提供框架。

研究背景：地下生态系统具有高度特有性，但探索不足；洞穴生物种群小、隔离且易灭绝，整合物种界定对生物多样性评估和保护至关重要。Dinaric Karst是欧洲最大连续喀斯特区域和全球地下生物多样性热点，其复杂造山过程形成持久扩散屏障，促进隔离分化、高物种丰富度、局域特有性和重复演化辐射。拟蝎是常见洞穴捕食者但研究不足，Neobisium spelaeum是首个被描述的洞穴拟蝎，N. stygium与其形态相似，亚种和异名关系经修订后仍显示分类边界不确定，线粒体数据亦提示隐藏多样性。问题在于形态趋同可掩盖演化轨迹，单基因或短片段可能不足界定物种且高遗传结构可能夸大隐存种数量，连续成种又模糊谱系边界。因此，研究人员以该复合群为模型，检验多次独立洞穴定殖与单次定殖–隔离分化假说，厘定基因组谱系，评估pedipalpal chela形态的分类效用，并识别受保护关注的狭域特有谱系。

研究内容与结论：研究人员采用两尺度策略，先基于COI（cytochrome c oxidase I，细胞色素c氧化酶I）构建Dinaric拟蝎区域系统发育，检验洞穴定殖是否在不同谱系中独立发生并建立目标复合群的演化背景；再对N. spelaeum–N. stygium复合群利用基因组SNP（single-nucleotide polymorphism，单核苷酸多态性）进行精细解析，整合群体遗传学、物种界定、种群统计学建模和形态学分析。结论是洞穴Neobisium呈多系，支持重复地下定殖；目标复合群形成支持良好的支系；全基因组数据揭示八个深度分化谱系和七个单例，线粒体-核不一致使演化轨迹重建复杂化；谱系间形态差异细微但显著；部分谱系在270 000–470 000年前经历瓶颈后扩张，与中更新世冰期旋回中的区域洞穴定殖相容。该研究突出Dinaric Karst为地下特有性中心，强调保护喀斯特景观和特有动物群，并为分类修订与保护规划提供框架。

关键技术方法（≤250字）：样本覆盖N. spelaeum和N. stygium分布区100余个洞穴点，超500号标本用于形态观察，110个体用于DNA提取，96个体用于ddRAD-seq（double-digest restriction-site associated DNA sequencing，双酶切限制性位点相关DNA测序），最终87个体纳入；同时分析维也纳NHM与柏林MfN的19件模式标本，并整合Hlebec et al. (2023)已发表COI数据。主要方法包括COI条形码与区域系统发育、ddRAD-seq SNP分型、群体结构分析（STRUCTURE、ADMIXTURE、dPCA、VAE）、单基因物种界定（ASAP、bPTP、GMYC）、ASTRAL物种树与BPP多物种溯祖界定、IBD（isolation by distance，距离隔离）与F_ST、Stairway Plot 2种群历史重建，以及线性和几何形态测量。

研究结果：
Molecular dataset summary and phylogenetic inference：COI条形码数据集含484个单倍型，其中239个变异位点、200个简约信息位点；最大似然树显示troglomorphic Neobisium为多系，N. spelaeum–N. stygium复合群形成支持良好的支系（UFBoot=95），并与其它洞穴种一起靠近地表种。复合群185条序列折叠为153个单倍型；种内p-distance（未校正成对遗传距离）最高达0.1618，四种形态种间最大p-distance为0.1664。贝叶斯树高度结构化且节点支持强，但四个形态种均非单系；ASAP、GMYC、bPTP分别界定54、51、61个实体，显示单基因界定过度分裂。ddRAD数据平均保留3.25×10⁶ reads/样本，最终获得207 927个SNP，其中12 895个非连锁SNP，缺失率32.71%；最大似然树按地理区域形成深度分化群，L5在树中为并系。

Population genomics and species boundaries：ADMIXTURE和STRUCTURE均识别K=8为最优遗传簇，定义为L1–L8；多数个体近完全归入单一簇，少数呈混合 ancestry。L5在ML树中并系，ADMIXTURE与STRUCTURE对其处理不同，研究人员依ADMIXTURE赋值。另有七个单例占据孤立系统发育位置。LDA（linear discriminant analysis，线性判别分析）显示与地理分布和地质边界一致的强遗传结构；VAE（variational autoencoder，变分自编码器）将高维SNP降至二维并分离出八个遗传群，L2、3、4、5、8显示地理亚结构。ASTRAL恢复八谱系关系但部分节点支持弱，尤其L5；BPP在MSC（multispecies coalescent，多物种溯祖）模型下以PP=1.0支持八个谱系，无谱系对达到合并阈值（gdi≤0.3）；姐妹对7–6的gdi=0.78–0.84，2–1的gdi=0.44–0.53；迁移模型推断7–6和1–2间存在低但非可忽略的双向基因流；N_e异质，θ=0.000567–0.004629；7–6和2–1分化时间相近且为最近分裂事件。

Isolation by distance and population differentiation：全部个体的遗传距离与地理距离显著相关（Spearman ρ=0.530，P<0.0001；Mantel r=0.535，P<0.0001），表明强IBD；遗传距离132–304，地理距离0.01–216 km。谱系间F_ST（Weir & Cockerham）为0.096–0.264，均值0.187±0.045；Hudson F_ST均值0.302，范围0.168–0.458，与前者强相关，说明种群结构在不同估计量下稳健。

Demographic analyses：Stairway Plot显示八个谱系均发生种群扩张，扩张后总体稳定；中位N_e从L5的2.5×10⁶至L8的7.7×10⁶。L1、3、4、7、8显示轻度下降：L8约270 000年前，L4约350 000年前，L1和L3约400 000年前，L7约470 000年前。

Phenotypic data and phylogenetic signal：线性形态测量PCA中PC1解释60.35%变异，主要为体型，PC2解释5.73%为形状；CVA（canonical variate analysis，典型变量分析）显示N. stygium csikii与N. meuseli形成明显形态簇，其他类群重叠较大。21个性状中5个宽度性状具显著性别差异，雌性更大。几何形态测量显示谱系间pedipalpal chela形状显著分化（置换MANOVA：F=3.21，R²=0.254，P=0.0001），性别间不显著；centroid size（中心大小）在谱系间差异显著。但系统发育信号弱且不显著（K_mult=0.676，P=0.348），phylomorphospace中多数谱系重叠，L2和L8较独特，L1、5、6、7在中央形态空间重叠。

讨论总结：COI树中Dinaric洞穴Neobisium非单系提示地下适应具有演化重复性，但研究人员指出未显式检验性状独立演化，不能假定八个基因组谱系各自独立获得troglomorphic性状；文中“多次定殖”指由地理和系统发育结构推断的重复洞穴进入，而非每个谱系独立获得洞穴形态。地质证据与区域重复定殖相符：中生代连续碳酸盐台地经渐新世–早中新世（34–20 Ma，峰期30–20 Ma）破碎，始新世flysch、上新世现代喀斯特和更新世洞穴系统形成；谱系分布追踪主要地质单元，指示祖先地表种群在破碎化后分化，并于中更新世区域进入洞穴。五个谱系在270 000–470 000年前非同步瓶颈后扩张，与冰期旋回相容；强IBD也可由祖先结构化地表种群继承解释，且时间不足以完全解耦。跨类群分布一致支持地质气候事件塑造区域地下多样性。形态挑战在于基因组揭示八谱系而形态分化细微；线粒体界定约55个实体，较基因组数据约七倍高估；形态相似可能反映获得地下表型后有限分化，而非目标复合群内重复趋同，需地表近缘种检验。线性形态主轴为体型，常用诊断性状无谱系结构；非体型变异多为性二型；几何形态虽显著但系统发育信号弱，chela ratios同塑，需分子数据指导物种界定。线粒体-核不一致使界定复杂：洞穴缺氧对线粒体氧化代谢的正选择、母系线粒体较小N_e加速排序、拟蝎线粒体基因重排和碱基偏好均可降低COI可靠性；COI仍可作条形码，建议将COI MOTUs视作基因组谱系内的操作单元而非物种，否则全球洞穴物种清单可能被高估。分类上八个谱系可能各代表独立种，但正式修订前不作分类动作；clade 8含N. spelaeum istriacum拓扑标本，单独提升会使N. spelaeum并系。多个谱系为微特有，限于单洞或洞系，应独立保护评估；隐藏多样性需纳入地下生物多样性评估并修订西北Dinaric Karst保护优先级。研究结论翻译：通过整合系统发育基因组学、种群统计学、形态测量学和生物地理学证据，研究人员揭示了地质破碎化、气候振荡和地下适应如何共同在世界上最丰富的地下特有性中心之一产生了非凡的隐藏多样性。