DRACH观察/期望值指标揭示两侧对称动物表观转录组调控的进化景观

《Annals of the New York Academy of Sciences》：A DRACH Observed/Expected Metric Reveals the Evolutionary Landscape of Epitranscriptomic Regulation Across Bilateria

N6-甲基腺苷(N6-methyladenosine, m6A)RNA修饰是一种动态的转录后调控机制，通过影响RNA稳定性、输出、剪接、翻译和降解来快速调节基因表达。在两侧对称动物中，m6A主要沉积于DRACH共有基序(consensus motif)，该基序包含18种序列变体。尽管DRACH基序的存在并不表明某个转录本当前或曾经被甲基化，但它反映了该转录本接受m6A介导调控的潜力。在此，研究人员引入了一种新的指标——DRACH观察/期望值(DRACH observed/expected, DRACH o/e)——用于独立于序列组成地量化DRACH基序的相对富集程度。研究人员应用该方法表征了太平洋牡蛎(Pacific oyster, Magallana gigas)这一m6A调控尚不明确的代表性冠轮动物(lophotrochozoan)物种中DRACH基序的丰度、分布和基因组组织。这些分析进一步扩展到涵盖主要两侧对称动物谱系并表现出对比生态特征的另外八个物种。比较分析揭示了m6A调控机制的强保守性，以及DRACH基序分布和功能富集的保守性与谱系特异性模式。本研究建立了一个稳健的框架，用于评估转录本m6A甲基化的潜力，同时克服了表观转录组的固有变异性，从而为探索其潜在的系统发育和/或生态驱动因素的比较分析铺平了道路。

论文解读文章

一、研究背景与科学问题

N⁶-甲基腺苷(N⁶-methyladenosine, m⁶A)是真核生物RNA中最普遍的修饰之一，由METTL3-METTL14异二聚体为核心的写入复合体(writer complex)催化，作用于DRACH共有基序(D = A/G/U, R = A/G, H = A/C/U)。m⁶A修饰通过招募不同的读取蛋白(reader proteins)调控RNA代谢的多方面过程，包括前体mRNA剪接、稳定性、输出、翻译和降解，并且可由擦除蛋白(eraser proteins)FTO和ALKBH5逆转。因此，m⁶A水平在不同发育阶段、细胞类型和生理状态下呈现动态波动。

然而，正是由于m⁶A标记的瞬时性和高度上下文依赖性，实际的m⁶A图谱仅提供了不断变化的表观转录组的一瞥，难以进行跨物种的直接比较。相比之下，底层的DRACH基序序列含量保持不变。此前研究表明，DRACH位点在人类和小鼠间高度保守，而m⁶A峰则不然。但DRACH基序分布在两侧对称动物(bilaterian)进化中的保守程度及其所揭示的表观转录组调控的进化历史仍不清楚。进展受限于来自系统发育远缘谱系（如冠轮动物）数据的稀缺，特别是与蜕皮动物(ecdysozoans)的大量文献相比。

太平洋牡蛎(Magallana gigas)作为一种双壳类模型，展现出高度保守的表观转录组通路，包括m⁶A机制、分布模式和生物学作用，特别是在发育过程、应激反应和剪接调控中。研究人员利用这一独特实验框架，开发了DRACH观察/期望值(DRACH o/e)指标，并应用于九种代表主要两侧对称动物支系的物种，以探究表观转录组调控的进化。

二、主要技术方法

研究人员开发了一种基于Python脚本的计算方法，用于计数FASTA格式基因组和转录组中的DRACH基序及其18种子基序(submotif)。为了消除GC含量和序列长度的偏差，研究人员创建了DRACH o/e指标，即观察到的DRACH数量除以基于随机分布和GC含量的期望值，并将该值标准化为z分数，以±2为阈值定义DRACH富集基因(DRACHen)和DRACH缺失基因(DRACHdep)。通过计算边缘度(edginess)指数量化基因在染色体上的位置偏好。利用NCBI保守结构域数据库(CDD)鉴定m⁶A机制蛋白的结构域架构。采用KOGMWU R包进行基于delta排名的真核直系同源群(euKaryotic Orthologous Groups, KOG)功能富集分析。研究涉及九个物种（见表1）：Acropora millepora、Apis mellifera、Ciona intestinalis、Culex pipiens pallens、Drosophila melanogaster、Homo sapiens、Magallana gigas、Mytilus trossulus和Ostrea edulis，其基因组数据均来自公共数据库。

三、研究结果

3.1 DRACH基序在M. gigas基因组和转录组中的特征**

在牡蛎基因组中，DRACH基序相对于非甲基化对照基序DRCCH和DRGCH显著富集，但与DRTCH无显著差异。子基序频率不均匀，AAACA及其对应变体显著过表达。所有DRXCH五聚体在染色体上均匀分布。在转录组中也观察到类似模式。DRACH基序在基因区和基因间区均比预期更丰富，而DRCCH和DRGCH则低于预期。这种组成偏向部分归因于牡蛎基因组的低GC含量(33.5%)。

3.2 牡蛎基因组中的DRACH o/e指标

DRACH o/e指标显示，牡蛎基因组中基因区的DRACH富集度(峰值约1.16)高于基因间区(约1.09)，支持DRACH基序因其在m⁶A介导的RNA调控中的作用而被进化保守的假说。通过z分数将DRACH o/e标准化后，定义了DRACHen(z > 2, 1685个基因)和DRACHdep(z < -2, 986个基因)两个群体。DRACHen基因优先定位于染色体末端附近，而非甲基化DRCCH基序则无此分布偏好。在转录本内部，DRACHen mRNA中的基序主要位于编码区(CDS)和3'非翻译区(3'UTR)。基因本体(Gene Ontology, GO)富集分析表明，DRACHen基因与基因表达调控、神经功能和细胞周期等相关，而DRACHdep基因仅与运输活性相关的少数功能有关。

3.3 m⁶A机制蛋白在后生动物中的保守性

全蛋白序列分析显示，核心写入组分METTL3和METTL14最为保守(序列一致性55%-70%)，而MACOM辅因子VIRMA和ZC3H13分歧最大。擦除蛋白ALKBH5在所分析的三种昆虫中缺失。读取蛋白YTHDC1、YTHDF2和YTHDC2呈中等保守性。结构域水平的保守性显著高于全长序列，其中METTL3的MT-A70结构域最为保守(79%-89%)。

3.4 DRACH o/e、DRACHen和DRACHdep在系统发育中的分布

DRACH并非在所有物种中都是最丰富的基序：在软体动物、C. intestinalis和C. pipiens中占优，而DRGCH在人类、蜜蜂、果蝇和珊瑚中占优。DRACH子基序丰度的种间成对比较显示出强相关性(Pearson r = 0.660-0.994)，高于非甲基化对照基序，表明甲基化基序受到更强的进化约束。基于DRACH子基序组成的聚类与系统发育关系一致。染色体分布方面，软体动物的DRACHen基因靠近染色体末端，人类的DRACHen基因靠近染色体中心，而昆虫中未检测到显著的位置偏好，表明DRACH富集基因的基因组组织具有支系特异性。

3.5 DRACH富集的功能意义与进化保守性

真核直系同源群(KOG)富集分析显示，细胞骨架组织相关基因在所有物种中均与DRACH富集基因关联，而无机离子转运和代谢类别则与DRACH缺失基因关联。其他功能类别呈现支系特异性模式：RNA加工、翻译和细胞周期控制在软体动物中富集于DRACHen基因，但在蚊子中则相反。Delta排名相关性分析显示，软体动物内高度相关(r = 0.838-0.876)，而与昆虫的相关性急剧下降，C. pipiens甚至接近零或负相关，表明功能部署受系统发育之外的其他因素影响。

四、讨论与结论

本研究提出DRACH o/e指标作为衡量转录本内在甲基化潜能(methylability)的指标，克服了m⁶A标记动态变异性的限制。在牡蛎中，DRACH基序在基因区域过代表明其被积极保守。DRACH富集基因位于染色体末端并与需快速紧密调控的过程（如细胞周期、基因表达和神经功能）相关，而DRACH缺失基因与看家功能相关，这与诱导型转录本通常比看家转录本具有更高m⁶A甲基化的观察一致。

跨物种比较表明，m⁶A调控机制在后生动物进化中深度保守，但DRACH富集基因的染色体分布呈现支系特异性模式，提示m⁶A-RNA甲基化可能不仅参与RNA加工，还可能与高阶染色质组织有关。DRACH子基序的强种间相关性表明甲基化基序受到更强的进化约束。功能分析显示，某些功能类别（如细胞骨架组织）的DRACH富集模式在所有物种中一致，而其他类别则呈现支系特异性，表明功能部署受额外进化或生态驱动因素塑造。研究未发现生殖策略(r/K选择)与DRACH相关特征之间的关联。

研究结论：本研究建立了比较框架，用于跨越进化尺度研究表观转录组调控，同时克服了RNA m⁶A标记的可变性和动态性。研究结果表明，虽然系统发育亲缘关系部分解释了功能类别内DRACH富集的保守性，但不能完全解释观察到的模式，强调了额外生态和进化因素的贡献。通过结合DRACH基序分布的比较分析和九种代表性后生动物物种m⁶A机制蛋白的综合表征，本研究为探究与RNA甲基化相关且可能由其驱动的进化约束提供了方法论框架。