芽殖物种形成、线粒体捕获和基因渗入在Asellus aquaticus物种复合体地表与洞穴谱系之间

《Cladistics》：Budding speciation, mitochondrial capture and introgression between surface and cave lineages in the Asellus aquaticus species complex

　　

祖先与演化衍生种群的重复对为验证关于演化确定性定律的假说提供了机会。Asellus aquaticus物种复合体是一个无脊椎动物模型系统，具有多次独立的地表到洞穴转换以及复杂且未解决的演化历史。本研究（研究人员）对所有已知的主要谱系和地表-洞穴种群对进行了采样，并利用全基因组双酶切限制性位点关联DNA（ddRAD）测序来提高系统发育分辨率。研究人员提出了一个稳健的系统发育骨架，这对于检验生态转换、性状演化和比较分析的假说至关重要。研究结果支持从地表祖先到洞穴无脊椎动物演化的一个新视角。首先，不同地表谱系之间的频繁混合可能妨碍为每次洞穴定殖确定一个单一的地表姐妹谱系。其次，在更新世（Pleistocene）零星杂交事件中，一些洞穴谱系的线粒体基因组（mitogenomes）已被邻近地表种群的线粒体基因组所取代。由此产生的A. aquaticus新系统发育视图是，连续混合的地表谱系与从普遍存在的地表祖先中芽殖出来的几个隔离洞穴谱系共存。一些洞穴谱系已被描述为独立物种。因此，物种复合体的分类结构是由生态驱动的芽殖物种形成塑造的，而不是通常的二分法谱系分裂序列。

研究背景：Asellus aquaticus物种复合体是一个无脊椎动物模型系统，包含多次独立的地表到洞穴（surface-to-cave）生态转换，但此前主要基于线粒体DNA（mtDNA）的系统发育研究导致谱系间关系不明确，且未能区分独立起源与杂交事件。这种不一致性阻碍了对洞穴演化机制（如选择、漂变和基因流）的深入理解。因此，研究人员旨在利用全基因组数据构建稳健的系统发育骨架，以澄清谱系间关系、检测基因渗入（introgression）并量化独立洞穴起源的次数。

研究内容及结论：研究人员对所有已知主要谱系及地表-洞穴种群对进行了大规模采样，通过双酶切限制性位点关联DNA（ddRAD）测序获得全基因组单核苷酸多态性（SNP）数据，结合线粒体细胞色素c氧化酶亚基I（COI）序列，采用种群结构分析、系统发育重建、分子定年和祖先性状重建等方法，揭示了芽殖物种形成（budding speciation）模式、多次线粒体捕获（mitochondrial capture）事件，并识别出四次独立的地表到洞穴转换。该研究为Asellus aquaticus模型系统提供了可靠的演化框架，纠正了以往基于mtDNA的误解，强调了在洞穴演化研究中结合核基因组与线粒体数据的重要性。论文发表在《Cladistics》。

关键技术方法：研究使用了双酶切限制性位点关联DNA（ddRAD）测序技术，从17个种群（11个地表，6个洞穴）的88个个体中获取全基因组SNP数据；同时扩增和测序线粒体COI基因（共259个个体，其中166个来自GenBank）。分析工具包括Admixture（推断遗传簇）、fineRADstructure（共祖矩阵）、RAxML和MrBayes（最大似然和贝叶斯系统发育树）、TNT（简约法）、MCMCtree和BPP（分子定年，含地质校准点）、Dsuite（ABBA-BABA检验及Fbranch分析用于基因渗入检测）以及corHMM（祖先性状重建）。样本涵盖欧洲主要分布区，包括斯洛文尼亚、意大利、罗马尼亚、阿尔巴尼亚、希腊等地的地表与洞穴种群。

研究结果：
**核和线粒体系统地理结构**：通过Admixture分析（K=16最为支持）和fineRADstructure共祖矩阵，识别出16个遗传簇，主要聚为五个谱系：跨阿尔卑斯谱系（Trans-Alpine Lineage）、巴尔干谱系（Balkan Lineage）、中欧谱系（Central European Lineage）、北迪纳里克谱系（Northern Dinaric Lineage）和经典地下谱系（Classical Subterranean Lineage）。其中经典地下谱系为核基因组树中新发现的谱系，包含三个洞穴种群（Old Subterranean Pivka、Subterranean Rak和Asellus kosswigi），而线粒体树未能恢复该谱系的单系性。Pivka河混合谱系（Pivka River hybrid lineage）因存在显著混合（核基因组主要来自跨阿尔卑斯谱系，另有16–35%来自北迪纳里克谱系）而被排除于后续核系统发育分析。

**核-线粒体系统发育不一致及基因渗入重建**：Shimodaira-Hasegawa检验显著支持核树与线粒体树不一致（ΔL=8722.5，P<0.0001）。基因一致性因子分析表明，不一致主要源于线粒体捕获而非随机误差。具体检测到三次基因渗入事件：1）Pivka河混合谱系的所有个体携带北迪纳里克谱系的线粒体单倍型，表明完全的线粒体捕获；2）经典地下谱系中的Old Subterranean Pivka和Subterranean Rak种群也携带北迪纳里克谱系的线粒体DNA，推测为过去杂交事件中的捕获；3）罗马尼亚洞穴物种A. infernus的线粒体DNA来自地理邻近的地表种群，但未完全取代。Dsuite的Fbranch分析还揭示了北迪纳里克谱系内两个地表种群（Cerknica湖与Ljubljansko Barje）间的近期基因流。

**分歧时间和地表-洞穴转换的时间**：利用地质校准点（晚中新世托尔顿期约11–7百万年作为根年龄；Sne?nik–Risnjak地块抬升5–3百万年；中欧谱系年龄1.65±0.55百万年）进行分子定年。结果表明，跨阿尔卑斯谱系与其他谱系的分歧约在7百万年前。祖先性状重建（等速率模型，单速率类别AICc最低）结合75%似然差异阈值，识别出四次独立的地表到洞穴转换：1）经典地下谱系（约4–1.5百万年前）；2）A. infernus（最大~2.9百万年前）；3）Young Subterranean Pivka（最大~1.1百万年前）；4）Subterranean Krka（最大~0.8百万年前）。Dollo模型（不允许洞穴到地表转换）给出相同结果。

**洞穴和地表谱系的分子多样性**：基于ddRAD基因型位点（249,871个过滤位点）的多样性估算显示，洞穴种群的观测杂合度（HO，平均0.0010）和核苷酸多样性（π，平均0.0010）均显著低于地表种群（HO平均0.0022，π平均0.0025），约降低一半（F_13,2=27.95，R²=0.6581，P<0.0001）。线粒体捕获和基因渗入的方向一致地从高多样性的地表谱系流向低多样性的洞穴谱系。

讨论总结：讨论部分强调了芽殖物种形成在分类结构中的主导作用：地表种群虽在核基因上高度分化，但因广泛杂交而保持为单一物种，而每次洞穴定殖均导致局部地下谱系芽殖，使地表祖先成为并系群。这一过程与植物的芽殖物种形成类似，不同于传统二分法分裂。新的系统发育视图要求重新评估以往基于线粒体数据的分类和演化结论，例如不应将谱系分裂时间直接等同于洞穴定殖时间（如Saclier等人估计A. kosswigi的洞穴转换年龄实际是物种复合体的冠群年龄）。此外，线粒体捕获和杂交事件的检测提供了对地下演化的新认识：尽管有过去的地表到洞穴基因渗入，但现有洞穴物种与地表祖先之间的杂交似乎很罕见（如同地分布的A. kosswigi与地表种群间无基因流）。研究最终澄清了真正独立的洞穴-地表种群对数量，纠正了以往因线粒体系统发育误导的估计。

研究结论翻译：因此，物种复合体的分类结构是由生态驱动的芽殖物种形成塑造的，而不是通常的二分法谱系分裂序列。在目前最可能的演化场景下，三个经典地下谱系（Classical Subterranean lineages）均未保留现存的祖先地表种群。剩余无争议的洞穴-地表种群对为：（i）Subterranean Krka与相邻的北迪纳里克谱系地表种群；（ii）A. infernus与罗马尼亚相邻地表种群；（iii）匈牙利Molnár János洞穴与相邻地表种群。