IF50.5｜陆生动物基因组趋同演化模式与时间框架

11次独立登陆事件揭示基因功能趋同，但半陆生与全陆生路径迥异。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

动物从水生到陆生的转变是生命史上最重大的事件之一，在节肢动物、脊椎动物、轮虫、软体动物、环节动物、线虫、缓步动物和有爪动物等多个门类中独立发生。这些平行演化实验为研究基因组适应中的偶然性与趋同性提供了独特机会。已有研究表明，表型适应如保水皮肤或角质层、免疫系统改造、骨骼与运动方式变化、代谢率提升、发育适应以及视觉适应等，在陆生谱系中广泛趋同。在基因型层面，基因创新、复制和丢失等基因组变化对重大演化转变至关重要，特定基因（如水通道蛋白）已被关联到多个类群的登陆过程。然而，与陆生植物相比，动物登陆的基因组基础仍缺乏系统表征。

科学问题

卡在哪里： 此前研究多聚焦于单个谱系或少数基因，缺乏跨动物界的系统比较。不同陆生事件中基因家族的功能趋同程度、基因丢失的普遍性及其在适应中的角色，以及半陆生与全陆生谱系在基因组适应模式上的差异，均未得到充分解析。此外，动物登陆的时间框架仍不明确，缺乏基于分子钟的统一标定。

本文要回答： 本文旨在通过比较154个动物基因组，系统揭示11次独立登陆事件中基因获得与丢失的趋同模式，区分半陆生与全陆生谱系的适应差异，并重建动物登陆的时间框架。

技术路线

作者主要用了：比较基因组学、同源基因家族推断、CAFE5、GO/Pfam功能注释、分子钟。

作者收集了154个基因组（151个动物、3个非动物全息动物），覆盖21个动物门，重点采样围绕登陆节点。通过OrthoFinder聚类得到483,458个同源基因家族（HGs），并利用Paps和Holland的解析脚本推断祖先节点的新颖、核心新颖和丢失HGs，用CAFE5推断扩张和收缩HGs。随后对每个登陆事件的新颖基因进行GO和Pfam功能注释，并通过InterEvo框架计算不同节点间共享的生物学功能交集，以识别跨谱系的趋同功能。同时，利用943个保守直系同源基因构建超矩阵，用IQ-TREE建树，再用MCMCTree进行分子钟分析，估计登陆事件的分歧时间。

图 1　Fig. 1。图1展示了154个物种的系统发育树，绿色文字标注11个陆生事件节点。每个节点旁用箱线图显示新颖、核心新颖、扩张、收缩和丢失HGs的数量。分支长度按时间缩放，并标注了陆生节点的分歧时间和95%最高后验密度区间。该图直观呈现了陆生事件在动物演化树上的分布及其基因周转规模。

核心亮点

亮点 1｜登陆节点基因周转显著

所有陆生谱系相对于其最近水生祖先均表现出大量基因获得（新颖基因和扩张）和基因减少（丢失和收缩），但不同类群幅度差异大。新颖基因在蛭态轮虫、线虫、四足类和陆生腹足类中特别高，而蛛形纲和六足类的新颖基因较少。基因丢失在线虫、缓步动物和有爪动物中最为显著。通过置换检验，陆生节点的新颖基因速率显著高于水生节点（P=0.0015），表明登陆事件伴随高水平的基因组可塑性。蛛形纲和六足类的低周转提示其演化可能更多依赖基因共择而非新基因产生。

图 2　Fig. 2。图2以柱状图对比每个陆生事件节点与其三个最近祖先的五类HGs数量。可见大多数陆生节点的新颖和扩张HGs高于祖先，而丢失和收缩HGs也普遍增加，但蛛形纲和六足类例外。该图强调登陆事件伴随高水平的基因获得与丢失，但不同谱系幅度不同。

亮点 2｜基因获得驱动功能趋同

在11个陆生节点中，新颖基因家族独立演化出相似功能。至少10个节点共享118个GO术语和26个核心新颖GO术语，涉及渗透调节、脂肪酸代谢、生殖、解毒、感觉接收和刺激反应。Pfam结构域分析同样发现神经递质门控离子通道、跨膜受体和细胞色素P450等结构域在多个节点独立出现。扩张基因家族中也发现趋同，如细胞色素P450、含黄素单加氧酶、谷胱甘肽S-转移酶和G蛋白偶联受体在至少4个节点扩张，这些功能与解毒、氧化应激和感觉刺激响应相关。

图 3　Fig. 3。图3a的UpSet图显示至少10个陆生节点共享的GO术语数量，新颖HGs有118个共享GO术语，核心新颖HGs有26个。图3b用树图将55个最特异的GO术语分为生物过程、细胞组分和分子功能三大类，突出离子转运、刺激反应、代谢和生殖等功能。图3c展示共享Pfam结构域，如离子通道和细胞色素P450。

亮点 3｜基因丢失标记陆生适应

基因丢失在陆生节点中普遍且数量超过基因获得。Dbl同源结构域基因家族在11个节点中的8个丢失，pleckstrin同源结构域基因家族在7个节点丢失，这两个结构域是Rho GTP酶鸟苷酸交换因子的组分，与再生和伤口愈合相关。叶绿素酶基因家族和Shugoshin C端结构域蛋白也在多个节点丢失，可能反映食性转变和生殖适应。收缩基因家族中，氯离子通道、碳水化合物磺基转移酶和褪黑素相关受体在至少4个节点收缩，涉及渗透调节、细胞外通讯和昼夜节律。

图 4　Fig. 4。图4a和4b的UpSet图分别展示跨陆生节点共享的扩张和收缩HGs。扩张HGs中，细胞色素P450、含黄素单加氧酶、谷胱甘肽S-转移酶和G蛋白偶联受体在至少4个节点扩张；收缩HGs中，氯离子通道、碳水化合物磺基转移酶和褪黑素相关受体在至少4个节点收缩。图中用绿色、蓝色和橙色区分全陆生、半陆生和混合组合。

亮点 4｜半陆生与全陆生路径分化

基于Jaccard距离的PCoA显示，半陆生与全陆生类群在GO术语和Pfam组成上部分分离，PERMANOVA检验显著（GO: R2=0.0995, P<0.01; Pfam: R2=0.0992, P<0.01）。半陆生类群的新颖基因功能富集于角质层重塑、视觉发育和应激反应，呈现广泛而多功能的工具包；全陆生类群则集中于神经元发育和离子膜稳态，功能集更精简。半陆生谱系间共享多种功能，而全陆生谱系间共享极少，表明半陆生适应具有趋同性，全陆生则更多依赖谱系特异性路径。

图 5　Fig. 5。图5的PCoA基于Jaccard距离展示61个陆生物种在GO术语（a）和Pfam（b）组成上的分布。半陆生类群（橙色椭圆）与全陆生类群（绿色椭圆）部分分离，PERMANOVA检验显著，表明两组在功能组成上存在差异，但组内离散度无显著差异。

亮点 5｜三个登陆时间窗口

分子钟分析揭示动物登陆发生在三个主要时间窗口：第一个窗口在中寒武世至中奥陶世（约515-436 Ma），线虫、多足类、六足类和蛛形纲率先登陆，与早期陆生植物出现时间重叠；第二个窗口在晚泥盆世至早石炭世（约464-338 Ma），环带纲环节动物和四足类独立适应陆地；第三个窗口在白垩纪（约181-39 Ma），蛭态轮虫和陆生腹足类登陆。每个窗口与特定生态背景相关，如早期隐花植物覆盖、泥盆纪洪泛平原和白垩纪温室气候。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： 比较基因组学、同源基因家族推断、CAFE5、GO/Pfam功能注释、分子钟
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事比较基因组学、进化发育生物学或动物适应性演化研究的研究者。

带走一句

比较基因组学中，功能层面的趋同比基因家族层面的趋同更普遍；分析陆生适应时，应同时关注基因获得与丢失，并区分半陆生与全陆生谱系的不同演化模式。