IF31.7｜真核生物6mA甲基化广泛存在且与转录激活相关

6mA在真核生物中并非罕见，而是与转录激活和H3K4me3共定位的古老表观标记。

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

研究背景

DNA甲基化是调控真核基因组的关键机制，最常见的是5-甲基胞嘧啶（5mC）。相比之下，N6-甲基腺嘌呤（6mA）在真核生物中的存在和功能一直存在争议，部分原因是方法学伪影导致结果相互矛盾。然而，一些单细胞谱系，包括纤毛虫、早期分化真菌和藻类Chlamydomonas，显示出稳健的6mA信号，这引发了关于其起源和进化作用的问题。在纤毛虫中，6mA由MT-A70家族的AMT1甲基转移酶催化，与转录起始位点下游的核小体间连接DNA相关，提示其与转录激活有关。但缺乏来自许多真核谱系的数据，限制了对6mA进化和功能的重建。

科学问题

卡在哪里： 此前研究受限于分类群采样不足和方法学伪影，导致对6mA在真核生物中的分布、进化起源和功能存在争议。特别是缺乏对主要真核超群的代表性物种进行碱基分辨率水平的系统分析，难以确定6mA是否广泛存在、其与转录的关系是否保守，以及其与5mC的进化关联。

本文要回答： 利用Oxford Nanopore测序在18种单细胞真核生物中绘制碱基分辨率的6mA图谱，重建6mA的进化历史，并阐明其与转录和染色质状态的关系。

技术路线

作者主要用了：Oxford Nanopore测序、系统发育分析、ChIP–seq、RNA-seq、抑制剂处理。

作者首先通过系统发育分析追溯MT-A70甲基转移酶的进化历史，发现AMT1通路可追溯到真核生物最后共同祖先（LECA）。然后选择18种代表主要真核超群的单细胞生物，包括15种编码AMT1的物种和3种缺乏AMT1的阴性对照，利用Oxford Nanopore测序获得全基因组6mA和5mC图谱。进一步结合RNA-seq和ChIP–seq数据，分析6mA与转录水平、组蛋白修饰（H3K4me3、H3K27ac）及核小体定位的关系。最后通过抑制剂实验检验H3K4me3对6mA沉积的影响。

图 1　Fig. 1。图1展示了MT-A70甲基转移酶的异源二聚体结构和系统发育树。AMT1与AMT6/7形成类似METTL3/METTL14的复合物，p1/p2参与靶向。系统发育树显示真核生物MT-A70分为六个家族，且与细菌序列聚类，支持单次细菌起源。

核心亮点

亮点 1｜AMT1通路起源于LECA并反复丢失

系统发育分析显示，MT-A70甲基转移酶家族在真核生物中分为六个家族（AMT1、AMT6/7、AMT2/5、METTL3、METTL14、METTL4），且大多数超群中都存在，表明它们在LECA之前已经复制。AMT1与AMT6/7及p1/p2蛋白共出现，形成保守的异源二聚体复合物。AMT1仅在具有稳健6mA信号的物种中发现，而缺乏AMT1的物种6mA水平接近背景。这一发现将6mA的起源追溯到LECA，并揭示其在多细胞和单细胞谱系中反复丢失。

图 2　Fig. 2。图2显示18种真核生物中6mA在不同序列背景下的全局水平。编码AMT1的物种在AT背景（尤其SATS）中6mA显著富集，而缺乏AMT1的物种6mA水平低且无序列偏好。在AMT1物种中，6mA在基因TSS下游富集，且与转录水平正相关。

亮点 2｜6mA富集于转录起始位点下游并与转录水平相关

在所有编码AMT1的物种中，6mA主要出现在ApT二核苷酸上，且偏好SATS序列背景。6mA在基因转录起始位点（TSS）下游富集，与转录活性正相关：高表达基因具有更高的6mA水平和更多甲基化ApT位点。相反，缺乏AMT1的物种6mA均匀分布，无TSS富集。这一模式在多个超群中保守，表明6mA与转录激活的关联是古老的。

图 3　Fig. 3。图3比较了五个物种在不同生活史阶段或温度条件下的6mA变化。大多数基因的6mA水平在阶段间差异极小，尽管有大量差异表达基因。例如S. punctatus中一些基因在某一阶段沉默但仍保留高6mA，表明6mA动态有限。

亮点 3｜6mA在转录动态变化中保持稳定

对五种具有不同生活史阶段或温度处理的物种（Acanthamoeba castellanii、Spizellomyces punctatus、Creolimax fragrantissima、Chromosphaera perkinsii、Naegleria gruberi）进行比较，发现尽管数千个基因在阶段间差异表达，但6mA水平变化极小。例如，在S. punctatus中，一些基因在一个阶段几乎完全沉默，但仍保留与活跃阶段相似的高6mA水平。这表明6mA与稳态转录相关，但对转录动态变化并非必需。

图 4　Fig. 4。图4展示6mA与5mC的分布关系。在同时具有两种修饰的物种中，5mC可能限于转座子或基因体，而6mA主要位于TSS下游。在具有基因体5mC的物种中，6mA和5mC在基因体上呈互补分布，6mA在5'端富集，5mC在3'端富集。

亮点 4｜6mA与5mC在基因组中占据不同功能域

在同时具有6mA和5mC的物种中，观察到两种模式：一种中5mC限于转座子而6mA标记表达基因（如Acanthamoeba、Naegleria、Spizellomyces、Chlamydomonas）；另一种中5mC也存在于基因体，但6mA主要占据基因体5'端，5mC在6mA减弱处达到峰值（如chlorophytes和Amoebidium）。这表明6mA和5mC在功能上分离，6mA与转录激活相关，而5mC与抑制或基因体甲基化相关。

图 5　Fig. 5。图5显示6mA与H3K4me3和H3K27ac的共定位。在多个物种中，6mA与H3K4me3在TSS下游重叠，且6mA的周期性模式与核小体间区域对应。在Naegleria中，6mA仅存在于少数基因，但这些基因具有最高的H3K4me3信号。

亮点 5｜6mA与H3K4me3共定位但不受其直接调控

ChIP–seq分析显示，6mA与H3K4me3在TSS下游紧密重叠，且6mA的沉积宽度与H3K4me3峰宽一致。在大多数物种中，6mA呈现周期性模式，与核小体间连接DNA对应。然而，在A. castellanii中使用WDR5抑制剂（OICR-9429和Piribedil）降低H3K4me3水平后，6mA水平无明显变化，表明H3K4me3并非6mA沉积的直接上游调控因子。

生信可带走

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

• 方法栈： Oxford Nanopore测序、系统发育分析、ChIP–seq、RNA-seq、抑制剂处理
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事表观基因组学、进化生物学或单细胞真核生物研究的科研人员。

带走一句

6mA是古老且保守的转录相关表观标记，与H3K4me3共定位，但并非转录动态所必需。分析6mA时应优先考虑AMT1通路的存在，并注意其与5mC的功能分离。