导语： N6-甲基腺苷（m6A）是大脑中含量最丰富的mRNA修饰之一，在神经发育和脑功能中具有重要作用。然而，由于技术限制，一直以来无法在构成大脑的单个细胞类型中全局分析m6A位点。2024年9月24日，《自然·神经科学》在线发表了一项来自杜克大学医学院的重要方法学资源论文（2024年第27卷2512–2520页）。研究者开发了一种小鼠模型（DART-seq小鼠），能够在任何感兴趣的组织中以单细胞分辨率进行转录组范围的m6A检测。利用该模型，他们绘制了小鼠不同脑区以及皮层单个细胞内的m6A图谱，发现不同脑区和细胞类型之间m6A甲基化具有高度共享性，但也识别出少量在神经元中差异甲基化的mRNA（编码神经元信号传导的重要调控因子），并发现小胶质细胞的m6A水平低于其他细胞类型。最后，他们还在老年小鼠中进行了单细胞m6A图谱分析，鉴定出许多随年龄发生变化的转录本。该方法为研究RNA修饰在神经系统发育、功能及衰老中的细胞类型特异性作用提供了强大工具。

研究背景：m6A修饰在脑中的重要性与技术瓶颈

m6A的生理功能

m6A是真核生物mRNA中最常见的内部修饰，由甲基转移酶复合物（METTL3/METTL14） 催化写入，可被阅读蛋白（YTHDF家族等） 识别，并被去甲基酶（FTO、ALKBH5） 动态擦除。在大脑中，m6A参与：

• 神经干细胞自我更新与分化

• 少突胶质细胞成熟与髓鞘形成

• 突触可塑性、学习与记忆

• 应激反应

• 其异常与阿尔茨海默病、帕金森病、胶质母细胞瘤等相关

技术瓶颈

传统m6A检测方法（如MeRIP-seq、miCLIP）存在以下局限：

• 需要大量RNA（无法应用于单个细胞）

• 依赖抗体，特异性差，背景噪音高

• 无法区分细胞类型（大脑由多种细胞类型组成）

近年来开发的DART-seq（利用APOBEC1-YTH融合蛋白，在m6A位点相邻胞嘧啶诱导C-to-U编辑）可在少量细胞中检测m6A，但此前仍无法达到真正的单细胞分辨率。

核心方法学突破：可诱导DART-seq小鼠模型

模型设计（图1a）

• 转基因小鼠：在CAG启动子控制下表达APOBEC1-YTH融合蛋白（含YTH结构域，特异性结合m6A）或APOBEC1-YTHmut（突变体，无法结合m6A，作为阴性对照）

• 诱导系统：Cre/loxP-STOP盒 + 他莫昔芬诱导Cre重组酶 → 实现时空可控的APOBEC1-YTH表达

• 原理：APOBEC1-YTH结合m6A位点后，在相邻胞嘧啶（通常为C+2位置）催化C-to-U脱氨。通过RNA测序检测C-to-U编辑位点，即可反推m6A位置

验证（Extended Data Fig. 1-3）

• 诱导效率：他莫昔芬注射5天后，C-to-U编辑率在脑组织中达15-40%，且编辑依赖于m6A结合（YTHmut几乎无编辑）

• 组织特异性：在皮层、海马、小脑、肝脏、肺、肾、心脏中均可诱导（Extended Data Fig. 2a-c）

• 无脱靶效应：非m6A位点编辑率极低（Extended Data Fig. 2d-e）

• 不影响全局m6A水平：METTL3表达不变（Extended Data Fig. 1e-f）

核心发现之一：脑区间的m6A图谱与保守性

数据集

• 三个脑区：皮层、海马、小脑（各4只小鼠）

• 鉴定m6A位点数：皮层17,878个、海马13,461个、小脑11,391个

• 编辑率分布：中位数约15-20%，小脑编辑率略低（Extended Data Fig. 4e）

脑区间共享与差异（Extended Data Fig. 4j, m）

• 大多数m6A位点（约60%）在所有三个脑区中均存在，提示m6A在脑区之间具有高度保守性

• 少数差异甲基化位点：主要富集于突触信号传导、离子转运、学习相关通路

• 独立验证：通过RT-qPCR验证了差异甲基化位点（如Paqr8, Ube2v1, Nova2）（Extended Data Fig. 4i）

核心发现之二：皮层单细胞m6A图谱的绘制

单细胞DART-seq（scDART-seq）

• 方法：结合split-pool条形码技术（而非微流控），对27,400个皮层细胞进行单细胞RNA-seq，同时检测m6A编辑

• 细胞类型识别：通过基因表达特征鉴定出11种细胞类型（兴奋性神经元、GABA能神经元、星形胶质细胞、少突胶质细胞、小胶质细胞、Cajal-Retzius细胞、室管膜细胞、内皮细胞、周细胞、外周巨噬细胞、平滑肌细胞）（Extended Data Fig. 5d-e）

细胞类型间m6A的共享性（图2, Extended Data Fig. 5g）

• 大多数m6A位点在所有细胞类型中均存在（约70-80%）

• 少数细胞类型特异性位点：主要出现在兴奋性神经元和小胶质细胞中

核心发现之三：小胶质细胞具有低m6A水平

关键数据（图3, Extended Data Fig. 6）

• 小胶质细胞的m6A位点数量显著低于其他细胞类型（即使校正测序深度后仍如此）（图3c-d, Extended Data Fig. 6d-g）

• 并非由于APOBEC1-YTH表达水平低（蛋白水平与其他细胞相似）（Extended Data Fig. 6j-l）

• 可能原因：小胶质细胞中m6A甲基转移酶复合物（METTL3、METTL14、WTAP）mRNA表达水平较低（图3h, Extended Data Fig. 6m）

• 功能提示：小胶质细胞是大脑的免疫细胞，低m6A可能与其独特的免疫应答和吞噬功能需求相关

核心发现之四：细胞类型间差异甲基化mRNA的鉴定

神经元 vs. 星形胶质细胞（图4, Extended Data Fig. 7）

• 鉴定出214个差异甲基化mRNA（主要在兴奋性神经元中高甲基化）

• 差异甲基化位点主要位于3‘ UTR和非编码RNA中（Extended Data Fig. 6q-r）

• 基因本体富集：突触传递、离子转运、神经发育、学习记忆（图4b）

• 代表性例子：

  • Gria1（AMPA受体亚基GluA1）：在神经元中高甲基化（图4d, Extended Data Fig. 7g-i）

  • Gria2（GluA2）：同样在神经元中高甲基化（Extended Data Fig. 7j-l）

  • Grin2b（NMDA受体亚基NR2B）：在神经元中高甲基化

  • Atp2b2（质膜钙ATP酶，与自闭症关联）：在神经元中高甲基化

兴奋性神经元亚型之间的差异甲基化（Extended Data Fig. 8）

• 不同皮层层次（L2/3、L4/5、L6）的端脑内（IT）神经元之间，以及锥体束（PT）神经元与IT神经元之间，存在差异甲基化位点

• 例子：Fos（即刻早期基因）和Nr4a2（核受体转录因子）在L2/3 IT神经元中甲基化水平低于L4/5和L6 IT神经元（Extended Data Fig. 8b-e）

核心发现之五：基于m6A谱的神经元异质性

m6A聚类揭示功能相关的亚群（图5, Extended Data Fig. 9）

• 根据m6A甲基化谱，兴奋性神经元可划分为4个亚群（图5a）

• 这些亚群并非简单地对应于已知的转录组亚型（L2/3 IT、L4/5 IT、L6 IT、PT、CT）（Extended Data Fig. 9a-b）

• 差异甲基化位点所在的基因富集于：钙离子信号通路、突触囊泡循环、轴突导向（Extended Data Fig. 9f）

• 例子：

  • Grin2b A284位点：在cluster 4中甲基化水平最高（图5f-g）

  • Gria2 A1397位点：在cluster 3中甲基化水平最低（图5h-i）

• 功能提示：m6A修饰的异质性可能精细调控同一类型神经元中的不同功能状态（如突触可塑性水平）

核心发现之六：年龄依赖性m6A变化（图6）

实验设计

• 老年小鼠：20月龄（n=4），与年轻小鼠（2月龄）比较

• 分析细胞类型：兴奋性神经元（占大多数）

关键结果

• 鉴定出2,055个随年龄差异甲基化的mRNA（图6a）

• 变化模式：

  • 约60%的位点甲基化水平随年龄增加

  • 约40%的位点甲基化水平随年龄降低

• 功能富集（图6b）：

  • 高甲基化mRNA富集于：线粒体功能、代谢、应激反应

  • 低甲基化mRNA富集于：突触传递、离子通道、学习记忆、神经发育

• 代表性例子：

  • App（淀粉样前体蛋白，AD关键基因）：随年龄甲基化水平降低（图6c-e, Extended Data Fig. 10i-j）

  • Lrp1（LDL受体相关蛋白，参与Aβ清除）：甲基化水平降低

  • Grin2b、Gria1、Gria2等突触受体：甲基化水平均降低

独立验证（Extended Data Fig. 10k）：通过体外DART-seq验证了App甲基化随年龄降低的趋势（虽未达统计学显著性，可能与样本量小有关）。

综合模型：m6A修饰的细胞类型特异性和年龄依赖性调控

保守性与特异性的平衡

• 普遍存在的m6A：参与基础细胞功能（转录、剪接、翻译）的mRNA在所有细胞类型中均甲基化

• 细胞类型特异性m6A：参与特化功能的mRNA（如神经元中的突触蛋白、小胶质细胞中的免疫应答蛋白）差异甲基化

小胶质细胞的低m6A

• 可能反映了其独特的转录后调控需求：

  • 需要快速响应免疫刺激（m6A可能增加mRNA稳定性，但过高的稳定性可能阻碍快速降解）

  • 可能依赖其他RNA修饰（如m5C、m1A）或非编码RNA进行调控

年龄相关m6A变化的可能后果

• 突触蛋白mRNA甲基化降低 → mRNA稳定性下降 → 突触蛋白合成减少 → 突触功能衰退

• 线粒体蛋白mRNA甲基化升高 → 可能影响线粒体动力学和能量代谢（与衰老相关）

• App低甲基化 → 可能影响APP mRNA翻译效率或稳定性，参与Aβ产生（需进一步验证）

对神经科学的启示

资源可及性

• DART-seq小鼠模型：可通过合作或商业渠道获取（具体见原文）

• 测序数据：GEO（PRJNA996152，多个子数据集）

• 蛋白质组数据：MassIVE（MSV000094375）

• m6A位点搜索数据库：http://dartsource.duke.edu

• 分析代码：https://github.com/mflamand/Bullseye

局限性与未来方向

1. C-to-U编辑的非绝对定量性：编辑率受APOBEC1-YTH表达水平、RNA丰度、测序深度等影响，需谨慎解释

2. split-pool条形码的低mRNA捕获效率：与10x Genomics相比，每个细胞的转录本数量较少，可能遗漏低丰度m6A位点

3. 老年小鼠样本量较小（n=4），且未进行认知行为测试，无法直接关联m6A变化与功能衰退

4. 缺乏功能验证：未进行m6A位点突变（如通过CRISPR将A突变为C，消除甲基化）来验证其对mRNA稳定性/翻译的影响

5. 小胶质细胞数量少：仅分析54个小胶质细胞，需扩大样本量确认低m6A现象的普遍性

未来方向：

• 功能丧失/获得实验：在特定细胞类型中敲除METTL3，结合scDART-seq和行为测试，验证m6A在突触可塑性和记忆中的因果作用

• 疾病模型应用：在AD小鼠模型（5xFAD、APP/PS1）中，绘制疾病进程中不同细胞类型的m6A动态变化

• 空间m6A图谱：结合空间转录组学，将m6A信息映射到组织解剖位置

• 转译至人类：尝试在人脑类器官或死后人脑组织中应用DART-seq，验证细胞类型特异性m6A模式的跨物种保守性

结语

这项研究开发了一种可诱导、单细胞分辨率的m6A检测小鼠模型（DART-seq小鼠），并首次绘制了小鼠大脑不同脑区和细胞类型（特别是11种皮层细胞类型）的m6A图谱。研究发现：m6A在脑区和细胞类型之间高度共享，但也存在少量具有细胞类型特异性的甲基化位点，主要富集于突触信号传导相关基因；小胶质细胞具有显著降低的m6A水平，可能与其免疫特性相关；基于m6A谱可将兴奋性神经元进一步分为功能相关的亚群；在衰老过程中，突触相关mRNA的m6A水平普遍降低，而线粒体/代谢相关mRNA的m6A水平升高。该模型为研究RNA修饰在神经系统发育、可塑性、衰老和疾病中的细胞类型特异性功能提供了强大的、可广泛应用的平台。

原始论文：Tegowski, M., Prater, A.K., Holley, C.L. et al. Single-cell m6A profiling in the mouse brain uncovers cell type-specific RNA methylomes and age-dependent differential methylation. Nat Neurosci 27, 2512–2520 (2024). https://doi.org/10.1038/s41593-024-01768-3

资源：https://github.com/mflamand/Bullseye ； http://dartsource.duke.edu