《Systematic Entomology》：Little devils of Neotropical forests: A first molecular phylogeny of stygnid harvestmen with UCE-museomics and a gigamatrix approach (Arachnida, Opiliones, Stygnidae)
编辑推荐：
Stygnidae科是新热带区盲蛛目中一个具有独特形态特征且多样性显著的类群，长期以来吸引着分类学界的关注，然而其内部关系以及基于传统形态学的分类体系在很大程度上尚未经过分子数据的验证。本研究中，研究人员利用来源于新鲜标本和历史标本的超保守元件（UCE）位点，
Stygnidae科是新热带区盲蛛目中一个具有独特形态特征且多样性显著的类群，长期以来吸引着分类学界的关注，然而其内部关系以及基于传统形态学的分类体系在很大程度上尚未经过分子数据的验证。本研究中，研究人员利用来源于新鲜标本和历史标本的超保守元件（UCE）位点，结合整合了传统标记数据的gigamatrix矩阵，首次重建了Stygnidae科全范围的系统基因组学框架。研究人员获得了一个得到强支持的Stygnidae科单系群，该单系群包含六个主要谱系。尽管由短深层内节点处严重的谱系分选不完全（ILS）所驱动的系统发育冲突使得基部关系的最终解析变得复杂，但在后续分析中，传统的两个亚科划分并未得到支持。相反，研究结果支持一个由Stygninae、Heterostygninae和Auraninae subfam. nov.组成的三分支系统，其中后两者形成姐妹群关系：（Stygninae, (Heterostygninae, Auraninae subfam. nov.)））。考虑到传统外部性状中观察到的潜在同塑性（homoplasy），研究人员将分子拓扑结构与雄性生殖器形态学的全面检验相结合，提出了三个亚科的更新分类，并正式重新定义了各亚科的鉴别特征。此外，gigamatrix分析增加了属级覆盖度，使研究人员得以解析额外属的系统发育位置，并为该科内部近期基于形态学的分类修订提供了稳健的分子支持。从进化角度来看，整体系统发育模式表明Stygnidae科可能起源于北部亚马逊地区，而短深层内节点和长末端分支的存在可能指示了新热带区一次古老的快速多样化事件。小型Stygnid谱系之间形成对比的系统发育位置可能反映了小型化的反复演化或小型体型祖先形态的持续保留，尽管这些替代假说尚需进一步检验。通过提供一个稳定的系统发育框架，本研究为这些独特的新热带区盲蛛的形态演化和生物地理学演化研究奠定了基础。
Stygnidae科隶属于盲蛛目（Opiliones）中的Gonyleptoidea总科，是新热带区最具物种多样性的谱系之一，已描述有效种超过2095种，约占该目总多样性的30%。该科成员以其独特的形态特征著称，包括彼此分离的眼、亚矩形背甲以及异常延长的捕食性触肢。然而，尽管Stygnidae科在分类学上备受关注，其内部系统关系却长期缺乏分子水平的检验。传统的形态分类将该科划分为Heterostygninae和Stygninae两个亚科，但这一划分是否成立亟待独立的分子数据验证。
为解决这一问题，Kim、Alegre-Barroso和Giribet等研究人员开展了本项研究，论文发表于《Systematic Entomology》。研究团队利用超保守元件（UCE）测序技术，结合gigamatrix策略，首次建立起Stygnidae科首个全范围分子系统发育框架，并基于分子结果与雄性生殖器形态学的整合，提出了全新的三亚科分类体系， including the newly established Auraninae subfam. nov.。这一研究不仅为该科的分类学提供了稳定的系统发育基础，也为理解新热带区生物多样性的形成机制提供了重要见解。
本研究使用的关键技术方法包括以下六个方面。第一，样本采集方面，标本来源于哈佛大学比较动物学博物馆（MCZ）的无脊椎动物馆藏，涵盖1929年至2024年间采集的92份标本，包括保存在95%乙醇中并于≤?20°C冷存的新鲜标本，以及保存在70%乙醇中于室温存放的历史标本。第二，DNA提取与文库制备方面，新鲜标本使用Qiagen DNEasy试剂盒提取基因组DNA，历史标本则采用专门的二氧化硅珠提取方案；随后使用Kapa HyperPlus或Hyper Prep试剂盒构建索引测序文库。第三，UCE靶向富集与测序方面，使用Opiliones专用探针组在60°C下进行16小时杂交捕获，最终文库在哈佛大学Bauer核心设施使用Illumina NovaSeq S4平台（150 bp双端测序）进行测序。第四，序列处理与矩阵构建方面，利用PHYLUCE v.1.7.1流水线进行质量控制、contig组装和UCE位点匹配，经Gblocks和CIAlign清洗后，分别构建了50%和70%物种覆盖度的M50矩阵（675个位点）和M70矩阵（208个位点）；同时通过挖掘非靶标读段中的legacy标记，构建了整合UCE与传统标记的gigamatrix。第五，系统发育重建方面，分别采用最大似然法（IQ-TREE，含1000次超快速自助法重复）、串联贝叶斯推断（ExaBayes）和多物种溯组树方法（ASTRAL）进行系统发育分析。第六，系统冲突评估方面，运用基因一致性因子（gCF）、位点一致性因子（sCF）、四重奏一致性因子（qCF）以及四聚类似然映射（FcLM）等方法，量化评估不完全谱系分选（ILS）和系统误差对系统发育重建的影响，并进行了近似无偏（AU）拓扑检验和置换FcLM检验以区分真实的系统发育信号与系统假象。
研究结果部分可按以下小标题归纳。
基部关系的拓扑冲突与解析。研究在所有核心UCE系统发育树中均获得了得到最大支持的Stygnidae科单系群（UFB 100%/BPP 1.0/LPP 1.0），并一致识别出六个主要谱系：Stygnoplus分支、Auranus分支、PlA分支（Planophareus + cf. Actinostygnoides）、Stygnus分支、PrV分支（Protimesius + Verrucastygnus）和PiO分支（Pickeliana + Obidosus）。然而，这三个主要分支之间的深层关系存在显著的拓扑冲突，形成了三种竞争性假说：假说1（H1，三分支系统）支持Heterostygninae与Auraninae互为姐妹群，再与Stygninae构成姐妹关系；假说2（H2，传统两亚科系统）支持Heterostygninae为其他所有Stygnid类群的姐妹群；假说3（H3，Auraninae早期分歧）呈现梯形树状，Auraninae谱系依次基部发出。研究通过综合分析判定H1为最可能的进化历史，其依据在于：H1在更严格的M70矩阵中表现出高度的分析稳定性；M50矩阵中H3的出现与Auraninae谱系特有的缺失数据模式相关，提示其为系统假象；而ASTRAL物种树在所有矩阵中一致支持H1拓扑。
早期辐射过程中的严重ILS基因组特征。基因一致性因子和位点一致性因子分析揭示了清晰的时空对比模式：在属级和主要分支的冠节点处一致性因子较高，但在四个深层分支间分歧节点处急剧下降，gCF值降低而sCF值趋近于约33%——即严重ILS下硬多歧的理论期望值。四重奏频率分析显示，这些深层节点的q1频率较低（≤41%），且替代拓扑呈近乎对称分布（q2 ≈ q3），卡方检验未显示q2与q3之间存在显著差异。引入迁移指数（IH-i）为零，而ILS指数（ILS-i）高达88.0%–96.0%，表明ILS而非基因渐渗是系统发育冲突的主要驱动因素。基于溯组近似计算，这些节点的共同祖先有效种群大小（Ne）约为分歧间隔时间（T）的五倍，提供了快速进化辐射的强有力人口统计学信号。
传统亚科分类的否定。为在高度ILS背景下评估竞争性假说，研究人员进行了四聚类似然映射（FcLM）分析，结果明确反对传统的两亚科分类系统（H2）：在四聚类分析中，支持H2拓扑的四重奏比例显著低下（M70中为9.0%）。置换FcLM检验进一步揭示，当保留缺失模式但随机化系统发育信号时，H2拓扑的支持度反而高于真实数据，表明其在gappy M50矩阵中的出现是由非系统发育偏差驱动的假象。由此，研究结论认为"tarsal claw-based"的传统分类标准所依据的性状实际上是该科的祖征（plesiomorphic），而非界定Heterostygninae的共有衍征。
冲突信号的解析：长枝吸引。研究通过分析个体标本水平的四重奏支持分布，发现Stygnoplus属于最敏感的类群，其系统发育位置极易因四重奏中其他类群的组成而波动。通过将Stygninae分支拆解为其组成谱系进行FcLM分析，研究揭示了PiO分支（Pickeliana + Obidosus）是长枝吸引（LBA）的具体来源：当包含PiO分支时，数据大幅度转向支持H3拓扑；而排除PiO分支后，Heterostygninae + Auraninae节点的自助支持率从65%显著提升至97%。这一后验分析验证了Stygnoplus与PiO分支之间存在特异的LBA假象。
Gigamatrix系统发育与Heterostygninae的恢复。通过gigamatrix策略将六个legacy标记整合入UCE骨架，研究的属级覆盖度从30%提升至47%，新增了28个分类单元。该方法成功恢复了在UCE-only分析中因数据质量不足而被排除的Heterostygninae代表性类群，包括Stygnidius（Heterostygnus的 senior synonym，即Heterostygninae的模式属）、Eutimesius等，首次为Heterostygninae的单系性提供了分子证据。此外，gigamatrix还使研究人员能够从1929年采集的历史标本中恢复部分核基因序列，将Fortia属成功归入Heterostygninae，纠正了其此前被错误归入Stygninae的分类地位。
Auraninae与小型Stygnid类群。研究将(Auranus + PlA)分支确立为一个独特的主要谱系，并基于其分子和形态学综合特征正式建立新亚科Auraninae subfam. nov.。该类群以微小体型（通常<2 mm）、局限于圭亚那地盾和北部亚马逊的分布，以及独特的雄性生殖器形态为特征。研究同时揭示小型Stygnid并非单系群：Verrucastygnus被稳固地归入Stygninae的PrV分支；而先前假设独立演化的北部大西洋雨林（NAF）微型Stygnid（包括Gaibulus和Iguarassua等属）在gigamatrix分析中被恢复为PiO分支内的单系群。这些差异显著的系统发育位置表明，小型化在Stygnidae科内至少独立演化了两次或三次，作为对落叶层微生境的趋同适应。
Stygninae与近期系统学研究的符合。排除Auraninae后的Stygninae形成得到稳健支持的单系群，但其内部三个主要分支（Stygnus、PrV和PiO）之间的确切分支顺序仍未解决。尽管如此，研究结果显示了与近期基于形态学的分类修订的高度一致性：包括Obidosus属的重新确立、Verrucastygnus hoeferscovitorum的分类转移、以及Gaibulus和Iguarassua归入PiO分支并与Pickeliana近缘等，均获得分子数据的强有力支持。
生物地理学分析。研究讨论了Stygnidae科的生物地理格局。系统发育树呈现"早期爆发"式的多样化模式，其特征为极短的深层内节点后跟随长的末端分支。三个主要Stygninae分支和Auraninae在亚马逊北部存在显著的地理重叠，结合Yapacana属（形态上与Stygnoplus最近的类群）在该区域的分布，研究提出Stygnidae科共同祖先可能栖息于北部亚马逊的假说。随后，Stygninae从该区域早期广泛扩张，而Heterostygninae向其北的加勒比小安的列斯群岛扩散，Auraninae则主要留守于其祖先分布区。研究还观察到体型与分布范围之间的 intriguing 相关性：微型化类群表现出高度的局域性特有分布，而较大体型的属则具有更广泛的分布范围，提示形态进化——特别是体型大小——可能与扩散能力或生态耐受性相关。在Heterostygninae内部，Stygnoplus属在小安的列斯群岛呈现出复杂的岛屿生物地理模式，其拓殖历史不符合简单的由南至北跳板模型，暗示了复杂的扩散事件或古地理事件相关的隔离分化的历史。
系统学修订。基于分子系统发育结果，研究对Stygnidae科的亚科界限进行了全面修订。鉴于传统外部形态性状（如跗爪的位置和形态、第三至四步足跗节的形状、以及感觉毛的分布等）易受功能或生态因素影响而同塑进化，研究将分类诊断的重点转向雄性生殖器形态。三个亚科各自具有独特的生殖器结构模式：Heterostygninae的lamina parva远端薄而清晰，常具顶端U形缺刻和基部缢缩；Auraninae的lamina parva薄而顶端呈杯状（calyx-shaped）；Stygninae的lamina parva远端增厚，顶端平截呈屋顶状。研究正式建立了Auraninae Kim, Alegre & Giribet subfam. nov.，以Auranus为模式属，包含Auranus、Actinostygnoides、Jime和Planophareus四个属；同时更新了Heterostygninae和Stygninae的属组成和鉴别特征，并将Fortia从Stygninae转移至Heterostygninae。对于Phareus、Jabbastygnus、Metaphareus等七个属，由于其形态特征介于不同亚科之间或缺乏雄性标本，暂列为Stygnidae科地位未定的属（incertae sedis）。
在讨论与结论部分，研究人员指出本研究首次为Stygnidae科提供了全面的分子系统发育框架。尽管严重的ILS使得基部关系的完全解析仍具挑战性，但综合严格的M70矩阵、多物种溯组物种树、详尽的一致性因子分析、四聚类似然映射以及置换检验等多重证据线，研究明确否定了传统的两亚科分类系统，支持建立包含Auraninae subfam. nov.的三亚科系统。为促进这一新系统与形态学的整合，研究从根本上重新评估了亚科的鉴别特征，指出传统的跗爪形态等外部性状易于趋同演化，而以雄性生殖器共有衍征为中心的分类标准更为可靠。Gigamatrix分析策略被证明能有效整合历史数据和公共数据库序列，显著扩展了分类单元覆盖度，并成功解析了多个关键属的系统发育位置。研究结果还为近期多项基于形态学的分类修订提供了独立的分子验证。关于生物地理演化，研究提出的北部亚马逊起源假说以及"早期爆发"式快速辐射模式，为理解新热带区生物多样性的形成提供了新的视角。对于体型进化的讨论，研究指出小型化在Stygnidae中的多次独立起源，但究竟是趋同适应还是祖先状态的保留，尚需更密集的取样和定量性状重建来进一步检验。最后，研究人员展望未来的研究方向，强调通过更密集的分类取样以"打断"长枝，有望进一步缓解LBA效应，提高系统发育分辨率，最终为Stygnidae科提供一个更加稳健的系统发育框架。
研究结论部分原文翻译如下：研究人员呈现了Stygnidae科盲蛛的首个全面分子系统发育。尽管分析强烈支持当前界定的科及其六个主要谱系的单系性，固有的"早期爆发"拓扑结构——由严重的ILS所驱动——使得完全解析基部关系变得困难。通过剖析潜在的系统发育信号，研究人员发现传统的两亚科分类被似然映射方法所否定。相反，基于更严格的M70矩阵和多物种溯组物种树，研究人员提出了三分支系统并建立了Auraninae subfam. nov.。这一独特的圭亚那-亚马逊落叶层谱系分支与Heterostygninae形成分支，两者共同代表Stygninae的姐妹群。至关重要的是，为使这一分子拓扑与形态学协调，研究人员从根本上重新评估了亚科的鉴别特征。研究人员得出结论，传统的跗爪形态等外部性状似乎易于趋同演化。通过将分类学重点转移到雄性生殖器的共有衍征上，研究人员验证了三亚科系统并为每个亚科建立了更新的鉴别特征。Gigamatrix分析被证明对于整合历史序列数据和增加分类单元取样是有效的。具体而言，这一扩展的属级覆盖度使研究人员能够将Fortia属归入Heterostygninae，同时为该科内部多项近期基于形态学的系统学推断提供了关键的分子证据。尽管深层节点仍然难以解析——可能是由于该群体古老的快速辐射——研究人员预期未来的研究通过更密集的分类取样来打断长枝，将实现更高的分辨率并为Stygnidae科提供一个更加稳健的系统发育框架。
