小鼠大脑中的单细胞m6A谱分析揭示细胞类型特异性RNA甲基化组和年龄依赖性差异甲基化

Single-cell m6A profiling in the mouse brain uncovers cell type-specific RNA methylomes and age-dependent differential methylation

作者信息Matthew Tegowski, Anna K Prater, Christopher L Holley, Kate D Meyer
PMID39317796
发布时间2024-12
DOI10.1038/s41593-024-01768-3

实验完整度

研究包含DART小鼠模型建立、跨脑区bulk DART-seq、单细胞scDART-seq、GLORI、质谱验证、功能验证(RNA甲基化聚类)和体内衰老分析等多个独立证据层级。

主要模型

转基因DART小鼠(APOBEC1-YTH) 小鼠大脑皮层单细胞(scDART-seq) 小鼠大脑皮层、海马、小脑组织 纯化的小鼠皮层神经元、小胶质细胞和星形胶质细胞

重点核对

DART小鼠APOBEC1-YTH转基因的诱导(他莫昔芬注射剂量20mg/kg,连续5天) 单细胞scDART-seq的细胞分群(标记基因鉴定细胞类型) m6A位点鉴定标准(Bullseye算法,RAC基序,至少3个细胞支持) 差异甲基化分析标准(Wilcoxon秩和检验,FDR<0.05,log2FC>0.58) 年龄依赖的差异甲基化比较(年轻8-9周龄 vs 老年14-15月龄)

摘要

N6-甲基腺苷(m6A)是大脑中丰富的mRNA修饰,在神经发育和脑功能中发挥重要作用。然而,由于技术限制,无法对构成大脑的各个细胞类型中的m6A位点进行全局分析。在这里,我们开发了一种小鼠模型,能够在任何感兴趣的组织中以单细胞分辨率进行转录组范围的m6A检测。我们利用这些小鼠绘制了不同脑区和单个小鼠皮层细胞中的m6A图谱,发现脑区和细胞类型之间存在高度共享的甲基化。然而,我们也发现神经元中少数差异甲基化的mRNA编码神经元信号传导的重要调节因子,并且发现小胶质细胞的m6A水平低于其他细胞类型。最后,我们在老年小鼠中进行了单细胞m6A作图,并鉴定出许多m6A随年龄变化而变化的转录本。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
大脑中不同脑区和细胞类型的m6A分布是否具有特异性,以及随着年龄增长m6A是否发生细胞类型特异的改变?
核心机制
m6A在不同脑区和大多数细胞类型中高度共享,但部分mRNA存在细胞类型特异性甲基化,小胶质细胞呈现低m6A水平,且谷氨酸能神经元中m6A模式的异质性独立于基因表达。
主要证据
通过DART转基因小鼠进行bulk和单细胞DART-seq,结合GLORI和质谱验证;发现小胶质细胞m6A减少,神经元与星形胶质细胞之间存在大量差异甲基化RNA,并鉴定出衰老过程中神经元特异性App mRNA甲基化降低。
研究意义
DART小鼠为体内全局m6A分析提供了首个遗传编码模型,可用于研究多种条件下m6A的变化,并揭示了单细胞水平上大脑m6A分布的新特征。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

模型构建与验证

生成他莫昔芬诱导表达APOBEC1-YTH的DART转基因小鼠,并验证m6A依赖的C-to-U编辑

将APOBEC1-YTH-HA序列插入Rosa26位点,与UBC-CreER杂交,通过他莫昔芬诱导表达;通过Western blot和Sanger测序验证表达和编辑。

2

跨脑区m6A谱分析

绘制不同脑区(皮层、海马、小脑)的转录组m6A位点,并比较区域间甲基化异同

对DART小鼠进行bulk DART-seq,使用Bullseye鉴定m6A位点,并进行差异甲基化分析。

3

单细胞m6A谱分析

绘制小鼠皮层单细胞水平的m6A图谱,比较不同细胞类型的甲基化模式

对年轻成年小鼠皮层进行scDART-seq(7,702个细胞),使用Parse Biosciences split-pool barcoding建库,并通过Seurat聚类鉴定细胞类型,使用Bullseye鉴定单细胞m6A位点。

4

验证小胶质细胞低甲基化

验证小胶质细胞m6A水平低于其他细胞类型的发现

通过流式分选纯化小胶质细胞和星形胶质细胞,使用GLORI和UPLC-MS定量m6A水平。

5

细胞类型间差异甲基化分析

比较不同皮质细胞类型之间的m6A差异,特别是谷氨酸能神经元和星形胶质细胞

使用Wilcoxon秩和检验比较每对细胞类型的平均单细胞%C2U值,鉴定差异甲基化RNA。

6

神经元亚型特异的m6A分析

探究不同皮质层谷氨酸能神经元亚型之间是否存在m6A差异

对谷氨酸能神经元亚型进行鉴定(基于标志基因),并比较各亚型间的甲基化谱。

7

基于m6A的细胞聚类

检验m6A模式是否构成独立于基因表达的细胞异质性

使用Signac对谷氨酸能神经元的m6A数据进行聚类,并比较m6A聚类与基因表达聚类的关系。

8

年龄相关的m6A变化分析

比较年轻和老年小鼠皮层不同细胞类型中的m6A差异

对老年小鼠(14-15月龄)进行scDART-seq,与年轻(8-9周龄)数据集成,使用Wilcoxon秩和检验鉴定年龄相关的差异甲基化RNA。

9

验证App mRNA的年龄相关去甲基化

验证神经元中App mRNA m6A水平随年龄降低的发现

通过RT-qPCR-based m6A定量和纯化神经元样品的验证,确认App甲基化下降的特异性。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
pCMV-APOBEC1-YTHAddgene131636
Ai9Addgene22799
他莫昔芬Sigma-Millipore--
TRIzolInvitrogen--
iScript逆转录超混合液Bio-rad1708840
CloneAmp 2x HiFi预混液Takara--
Qubit蛋白分析试剂盒Thermo Fisher--
抗GAPDH抗体Proteintech10494-1-AP
抗HA抗体Cell Signaling Technologies3724
山羊抗兔IgG-HRPAbcamab6721
山羊抗小鼠IgG-HRPInvitrogen62-6520
Amersham ECL Prime试剂GE Amersham--
Chemidoc MP成像系统BioRad--
Thermo Vanquish NeoThermo--
Thermo Orbitrap AstralThermo--
Spectronaut----
NEBNext Ultra II 定向RNA文库制备试剂盒NEB--
NEBNext多重寡核苷酸(Illumina)NEB--
NovaSeq6000Illumina--
蔗糖----
CM3050S-3-1-1Leica--
兔抗HA抗体Cell SignalingC29F4
鸡抗Iba1抗体Synaptic Systems234 009
豚鼠抗GFAP抗体Synaptic Systems173 004
AlexaFluor488偶联山羊抗鸡抗体Thermo-FisherA-11039
AlexaFluor568偶联山羊抗豚鼠抗体InvitrogenA-11075
AlexaFluor647偶联山羊抗兔抗体Thermo-FisherA-21245
Leica DMi8倒置STED和共聚焦显微镜Leica--
Imaris x64Bitplane--
抗Cd11b-Alexa Fluor 488抗体Biolegend101217
抗Acsa2-APC抗体Miltenyi130-099-138
TruStainFcX PlusBiolegend--
FACS Aria A02BD--
抗NeuN-AlexaFluor 488抗体Abcamab190195
FACSymphony S6BD--
Dynabeads mRNA纯化试剂盒Invitrogen--
核酸酶P1Sigma--
南极磷酸酶NEB--
Xevo TQ-S系统Waters--
RNeasy MinElute纯化试剂盒Qiagen--
BstI聚合酶NEB--
SuperScript III逆转录酶Invitrogen--
全细胞固定试剂盒Parse Biosciences--
Evercode WT试剂盒Parse Biosciences--
Evercode WT Mini试剂盒Parse Biosciences--
SMARTer链特异性总RNA测序试剂盒v3 -Pico输入哺乳动物Takara--
NovaSeqXIllumina--

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
DART小鼠诱导
他莫昔芬剂量(20mg/kg)、注射频率(每日一次,连续5天)、诱导后组织收集时间(7天后)
阅读提示:Methods:Tissue collection for Sanger sequencing, immunoblotting, and bulk DART-seq
bulk DART-seq
每个脑区的生物重复数(4个)、RNA起始量(1μg)、文库制备试剂盒(NEBNext Ultra II)、测序平台(NovaSeq6000)
阅读提示:Methods:Bulk DART-seq library preparation;m6A位点鉴定标准(Bulk DART-seq m6A site identification)
scDART-seq
细胞固定和建库试剂盒(Evercode WT Kit)、测序平台(NovaSeq6000)、单细胞m6A位点标准(覆盖≥3 reads,编辑率0.5%-95%,至少3个细胞支持)
阅读提示:Methods:Preparation of scDART-seq libraries from cortex of young mice;scDART-seq m6A site identification
细胞类型鉴定
Seurat分析参数(如PC数、分辨率)、标记基因列表
阅读提示:Methods:Identification of cell types using scDART-seq data
差异甲基化分析
统计检验(Wilcoxon秩和检验)、FDR阈值(0.05)、log2FC阈值(0.58)、最小细胞数(10)
阅读提示:Methods:Differential methylation analysis of cell types and neuronal subtypes;Differential methylation between the young and aged cortex
小胶质细胞低压验证
细胞分选方案(anti-Cd11b和anti-Acsa2)、GLORI实验条件(5-10ng RNA)、UPLC-MS样本量(n=2)
阅读提示:Methods:Isolation of microglia and astrocytes;m6A profiling using GLORI;Mass spectrometry for m6A quantification
神经元亚型分析
神经元亚型的标志基因(如Slc30a3、Rorb、Bcl6等)
阅读提示:Methods:Identification of cell types using scDART-seq data(神经元亚型部分)