小鼠组织的三元编码DNA甲基化图谱

A ternary-code DNA methylome atlas of mouse tissues

作者信息Sol Moe Lee, David C Goldberg, Cameron Cloud, Jared B Parker, Christopher Krapp, Christian E Loo, Elliot Kim, Ivan Zhao, Chengcheng Jin, Rishi Porecha, Marisa S Bartolomei, Rahul M Kohli, Wanding Zhou
PMID41057935
发布时间2025-10-07
DOI10.1186/s13059-025-03808-y

实验完整度

研究包含多组织、多年龄、两种性别的群体样本,并通过bACE阵列、ACE-seq验证、随机森林模型及多组学关联分析,形成独立证据层级。

主要模型

C57BL/6J小鼠 29种小鼠组织 小鼠组织样本 转录组数据

重点核对

样本数量:265个组织样本,来自32只C57BL/6J小鼠 组织类型:29种,涵盖大脑皮层、小脑、皮质下脑、肝、血、皮肤、眼、股骨、心、肺、睾丸、视神经等 年龄范围:8-76周 性别:两种性别 技术验证:与TAB-seq和ACE-seq比较,Spearman相关系数达0.78

摘要

背景:DNA胞嘧啶修饰,包括5-甲基胞嘧啶(5mC)和5-羟甲基胞嘧啶(5hmC),是具有不同功能的关键表观遗传调控因子。由于传统基于亚硫酸氢盐转化的分析方法存在局限性,解析跨组织和细胞类型的三元编码(C、5mC、5hmC)仍然是一个关键优先事项。结果:在此,我们利用结合亚硫酸氢盐和酶(bACE)转化与小鼠甲基化BeadChip,生成了覆盖29种小鼠组织类型、年龄8–76周、两种性别的265个5mC和5hmC碱基分辨率三元编码修饰图谱。我们的图谱揭示了5hmC分布的复杂语法,受细胞有丝分裂活性、染色质状态以及与同一和邻近CpG位点的5mC相互作用共同塑造。值得注意的是,我们证明5hmC显著补充了基于5mC的生物标志物在描绘大脑和非大脑组织细胞身份方面的作用。每种修饰状态,包括单独的5hmC,都能准确区分组织类型,实现高精度的机器学习分类表观遗传身份。此外,三元甲基化组的变异广泛涉及基因转录变异,与年龄相关的变化与基因表达以组织依赖性方式相关。结论:我们的工作揭示了组织、性别和年龄如何共同调控两种胞嘧啶修饰的动态,扩展了DNA修饰生物标志物发现的范围,并为探索发育和疾病中的表观遗传动态提供了参考图谱。

实验结论

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

研究问题
如何精确解析组织中5mC和5hmC的分布,并探索其与组织身份、基因表达和衰老的关系?
核心机制
5hmC由5mC氧化生成,其分布受细胞增殖率、染色质状态和5mC底物可用性影响,在长寿命组织中积累,并与基因表达正相关。
主要证据
基于265个样本的bACE阵列图谱,结合转录组数据,通过相关性分析、机器学习分类和染色质状态富集分析,证实5hmC独立或辅助5mC定义细胞身份。
研究意义
为发育和疾病中的表观遗传机制研究提供参考图谱,扩展DNA修饰生物标志物发现范围。

研究路径

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

1

构建bACE阵列图谱

生成覆盖多组织、多年龄、性别的三元编码甲基化图谱

从32只C57BL/6J小鼠中收集29种组织共265个样本,提取DNA后进行亚硫酸氢盐转化和bACE转化,使用Infinium Mouse Methylation BeadChip进行甲基化芯片分析,并通过插值校正获得5modC和5hmC值。

2

验证bACE阵列准确性

评估bACE阵列测量5hmC的准确性

将bACE阵列数据与已发表的TAB-seq和ACE-seq数据进行比较,并采用低通量ACE-seq对四种组织进行正交验证。

3

分析全局修饰水平及局部邻域特征

揭示5mC和5hmC的全局分布及其局部空间相关性

通过三元图、相关分析、局部邻域分析等方法,解析不同组织中C、5mC、5hmC的比例及5mC和5hmC在CpG位点的空间分布。

4

分析染色质与基因上下文

探索5mC和5hmC在不同染色质状态和基因区的分布

将CpG位点与18状态染色质模型相交,分析不同组织中5hmC在基因体、增强子、启动子等的分布,并关联基因表达。

5

组织身份预测

评估5hmC单独能否定义细胞身份

基于5hmC数据,使用随机森林模型对8种组织类型进行训练和预测,计算SHAP值分析特征重要性。

6

识别组织特异性修饰位点

找出各组织特异性的高甲基化和低甲基化位点

使用一对多(OvR)方法计算每个CpG的AUC,挑选Δβ值超过阈值的位点作为组织特异性标记。

7

关联转录组和修饰

探究5mC和5hmC与基因表达的关系

将甲基化数据与公开转录组数据整合,计算CpG甲基化与基因表达的相关性,并分析年龄相关变化。

8

研究衰老相关修饰变化

揭示年龄相关的5mC和5hmC动态及其与基因表达的关系

分析不同年龄样本的甲基化变化,计算线性回归斜率,并关联基因表达变化。

研究方法

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

产品清单

实验环节名称品牌货号
Qubit 4.0荧光计Invitrogen--
dsDNA高灵敏度检测试剂盒InvitrogenQ33231
AMPURE XP磁珠Beckman CoulterA63881
EZ DNA甲基化试剂盒Zymo ResearchD5001
EZ-96 DNA甲基化MagPrepZymo ResearchD5040
A3A酶文中未明确说明--
二甲基亚砜Sigma-AldrichD2650-100ml
寡核苷酸纯化浓缩试剂盒ZymoD4060
Platinum II PCR预混液(2X)Thermo Fisher14000012
ROXInvitrogen12-223-012
EvaGreen染料EMSCONC0521178
384孔PCR板Applied BiosystemsA36931
QuantStudio 5Thermo Fisher--
Infinium小鼠甲基化BeadChip试剂盒Illumina--
Covaris超声破碎仪Covaris--
NEBNext UltraII DNA文库制备试剂盒NEBE7645
bGTNEBM0357
UDP-葡萄糖NEBS2200
CutSmart缓冲液NEB--
甲酰胺Thermo Fisher17899
APOBEC 10X缓冲液NEBE7134
牛血清白蛋白NEBE7135
APOBEC酶NEBE7133
KAPA HiFi U+聚合酶Roche50-196-5287
用于酶促甲基化测序的NEBNext多重寡核苷酸NEBE1720
高灵敏度D1000 ScreenTapeAgilentNC1786959
Illumina MiSeq测序仪Illumina--
异丙醇MilliporeSigmaEM1.09634.1011
GlycoBlue共沉淀剂InvitrogenAM9515
乙酸铵溶液Sigma-AldrichA2706
乙醇MilliporeSigmaEM1.00983.1011
Tris缓冲液pH 8.0VWR97062-674
CpG甲基化酶ZymoE2010
未甲基化Lambda DNAPromegaD1521
基因组DNA清洁浓缩试剂盒-10ZymoD4010
HpaII限制性内切酶NEBR0171S
琼脂糖VWR0710-500G

关键环节

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

环节核对要点
样本收集与DNA提取
样本数量、小鼠品系、年龄、性别、组织类型、DNA提取方法
阅读提示:Methods: Sample preparation and DNA extraction
亚硫酸氢盐和bACE转化
转化试剂盒、A3A浓度、反应条件、纯化方法
阅读提示:Methods: Bisulfite conversion of genomic DNA; A3A deamination of BCD for bACE array
芯片分析
芯片类型、数据处理流程、插值标准
阅读提示:Methods: BS array and bACE array by Illumina MM285 array and data pre processing; Methylome of Illumina MM285 array interpolation and evaluation
组织身份预测模型
训练/测试集划分、特征选择标准、模型参数、交叉验证方法
阅读提示:Methods: Tissue prediction
组织特异性签名分析
AUC计算、Δβ阈值、缺失值处理
阅读提示:Methods: Tissue-specific hyper- and hypomethylation signatures
甲基化与基因表达关联分析
转录组数据集、基因筛选标准、CpG筛选标准、统计方法
阅读提示:Methods: DNA methylation-gene expression correlations
衰老分析
年龄分组、线性回归模型、p值阈值
阅读提示:Methods: DNA methylation-gene expression correlations