IF50.5｜毛颚动物独特体型的基因组起源

毛颚动物通过大规模基因丢失、串联重复和转拼，重塑了细胞类型与调控网络。

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

研究背景

动物门类的出现往往伴随着新体型结构的建立，但驱动这些演化的基因组基础仍不清楚。毛颚动物（箭虫）是一类自由生活的掠食性海洋无脊椎动物，其系统发育位置长期存在争议，直到最近才被归入有颚动物支系（Gnathifera），作为螺旋卵裂动物（spiralians）中其他冠轮动物（lophotrochozoans）的姐妹群。毛颚动物具有许多独特的形态特征，如几丁质捕食棘、纤毛感觉器官和复杂的神经系统，但其体型自寒武纪以来几乎没有变化。相比之下，同属有颚动物的轮虫体型微小且基因组演化迅速。此前对轮虫的研究揭示了快速的基因组演化与无性生殖的关联，但毛颚动物的基因组演化模式及其对细胞类型和体型的影响尚不清楚。

科学问题

卡在哪里： 尽管已知毛颚动物属于有颚动物，但其独特的体型如何从基因组层面产生仍不清楚。此前缺乏染色体级别的基因组和单细胞数据，无法解析基因丢失、重复和调控变化对细胞类型演化的贡献。

本文要回答： 本文旨在通过染色体级别基因组、单细胞图谱和表观基因组数据，揭示毛颚动物独特体型的基因组起源。

技术路线

作者主要用了：染色体级别基因组组装、Hi-C、ATAC-seq、单细胞RNA-seq、甲基化组、原位杂交。

作者首先组装了毛颚动物 Paraspadella gotoi 的染色体级别基因组，并利用 Hi-C 数据推断染色体三维结构和着丝粒位置。随后，通过比较基因组学分析基因家族得失和染色体融合事件，揭示有颚动物支系的加速演化。为了将基因组变化与细胞类型联系起来，作者生成了涵盖幼体和成体的单细胞转录组图谱，鉴定出约30种细胞类型，并通过跨物种比较区分祖先和毛颚动物特有的细胞类型。最后，结合甲基化组和 ATAC-seq 数据，分析 DNA 甲基化重编程和转拼（trans-splicing）对基因调控的影响。

图 1　Fig. 1。图1展示了毛颚动物基因组特征和演化。a为成体照片；b为基因家族重建，显示有颚动物支系大量基因丢失；c为染色体图，显示DNA转座子密度、甲基化水平和基因密度，着丝粒区域富集转座子且高甲基化；d为丢失基因家族的GO富集，着丝粒组装相关通路显著；e为祖先连锁群保留情况，毛颚动物染色体由多个祖先连锁群融合而成；f为Hi-C互作热图，显示染色体臂间的区室化。

核心亮点

亮点 1｜有颚动物支系基因丢失与染色体融合

毛颚动物基因组仅257 Mb，注释了22,072个蛋白编码基因。与其它螺旋卵裂动物比较发现，有颚动物共同祖先经历了2,542个基因家族的丢失，其中与神经和感觉功能相关的家族显著收缩。尤为突出的是，着丝粒组装的关键基因 CenH3 和 CENP-T 在有颚动物中完全丢失，但毛颚动物仍具有局部着丝粒，这些区域富集DNA转座子并高度甲基化，提示新着丝粒采用了替代的分子机制。此外，毛颚动物染色体由多个祖先连锁群融合而成，保留了部分螺旋卵裂动物特有的融合事件，但轮虫中祖先连锁群完全被打乱，表明有颚动物支系染色体演化速率加快。

图 2　Fig. 2。图2为单细胞图谱。a为UMAP降维图，显示约30个细胞簇；b为标记基因点图，突出用于原位杂交验证的基因；c为SAMap跨物种比较，显示与其它原口动物共享的细胞类型；d为原位杂交结果，验证了神经元、肌肉、表皮等标记基因的表达模式。

亮点 2｜单细胞图谱揭示细胞类型创新

单细胞转录组分析鉴定了约30种细胞类型，包括神经元、肌肉、表皮、纤毛细胞、生殖细胞等。跨物种比较显示，神经元、肌肉、肠道等细胞类型与其它原口动物共享保守的转录因子表达，但毛颚动物特有的捕食棘、感觉乳突和纤毛感觉器官等细胞类型在其它动物中无对应。这些特有细胞类型表达大量毛颚动物特异性基因和串联重复基因，如机械感觉离子通道 Piezo 的8个拷贝中有一个在感觉乳突中特异表达。功能富集分析表明，重复基因富集于离子跨膜转运和发育相关通路，提示基因重复与细胞类型创新密切相关。

图 3　Fig. 3。图3聚焦基因重复与细胞类型创新。a为重复基因的GO富集，离子转运和发育通路显著；b为祖先连锁群在染色体上的分布，重复基因不遵循祖先连锁模式；c为同义替换率分布，显示单一峰，支持单次重复事件；d为染色体间重复基因对数量，无全基因组复制特征；e为重复基因对在染色体上的位置，部分相距很近，符合串联重复。

亮点 3｜串联重复爆发而非全基因组复制

毛颚动物经历了3,379个基因家族的扩张，规模与脊椎动物2R事件相当，但染色体分布和同义替换率分析不支持全基因组复制。重复基因对在染色体上随机分布，约25%位于同一染色体，且没有染色体对显示相互富集。重复基因间距离呈双峰分布，部分相距约10 kb，符合串联重复特征。此外，重复基因的内含子数量未减少，排除了逆转录转座爆发。Hox基因簇的扩张是串联重复的典型例子：毛颚动物拥有14个Hox基因，其中中段和后段基因显著扩增，且MedPost基因位于中段和后段基因之间，支持其作为有颚动物共衍征的地位。

图 4　Fig. 4。图4展示Hox基因簇扩张和表达。a为螺旋卵裂动物Hox基因组成比较，毛颚动物中段和后段基因显著扩增；b为六个Hox基因的整体原位杂交，显示前段基因在神经系统广泛表达，中后段基因在腹神经索呈阶梯式表达，MedPost在感觉乳突中表达；c为表达模式总结图。

亮点 4｜DNA甲基化重编程至转座子

毛颚动物的DNA甲基化模式与大多数动物不同，基因体几乎不甲基化，而转座子（尤其是LTR逆转录转座子）被高度甲基化。这种模式类似于真菌和少数动物（如早期线虫和栉水母）。DNA甲基化工具酶发生结构域简化：DNMT1和UHRF1的旁系同源物丢失了CXXC锌指或泛素结构域，DNMT3的6个拷贝均丢失ADD和PWWP结构域，后者通常负责将甲基化引导至基因体。这些变化可能导致DNA甲基化功能从转录调控转向转座子沉默。

图 5　Fig. 5。图5展示DNA甲基化和转拼。a为甲基化水平在基因和转座子上的分布，显示转座子特异性甲基化；b为DNA甲基化工具酶的结构域简化；c为转拼和操纵子的基因组分布；d为不同剪接前导序列的使用情况；e为操纵子基因的功能富集；f为转拼基因的表达特异性；g为生殖细胞中标记基因的转拼比例；h为转拼基因的ATAC-seq峰数量；i为转拼基因上游的富集基序。

亮点 5｜转拼与操纵子广泛存在

毛颚动物中近一半基因（10,197/21,072）发生转拼，18%的基因位于由转拼解析的操纵子中。转拼基因表达水平更高，但细胞类型特异性较低，且启动子区域ATAC-seq峰更少，提示其调控复杂度降低。操纵子基因富集于氮和蛋白质代谢等管家功能，而发育和神经系统基因被排除在外。转拼基因在生殖细胞中富集，可能有助于在基因体甲基化丢失的情况下稳定管家基因表达。

生信可带走

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

• 方法栈： 染色体级别基因组组装、Hi-C、ATAC-seq、单细胞RNA-seq、甲基化组、原位杂交
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 研究动物演化、基因组结构、单细胞图谱和表观遗传调控的生信与演化发育生物学研究者。

带走一句

基因丢失和串联重复可以重塑细胞类型，而DNA甲基化重编程和转拼可能补偿调控复杂度的降低。