WGBS+ChIP-seq+RNA-seq多组学解密H3K36me2调控早期胚胎DNA甲基化重建

熟悉我的人都知道，我这两年一直在表观遗传组学方向泡着。前几天易基因发布了颉伟、卢绪坤、张宇团队发表在Nature Cell Biology上的那篇文章解读，IF 19.1，标题很长：WGBS+ChIP-seq+RNA-seq等揭示早期胚胎发育过程中H3K36me2调控DNA甲基化重建机制。群里几个做生殖发育和干细胞研究的朋友直接炸了，问这个技术方案到底怎么搭的、微量样本怎么扛下来的、数据整合有没有能抄的作业。

这篇文章值得拆的东西太多了。表面看是一篇高水平论文，实际上是一套完整的表观遗传多组学打法。今天不说官话套话，我就从一个常年帮人跑组学项目、天天跟微量样本和数据质控打交道的人的角度，把这篇文章背后的技术逻辑、实验设计思路、实操细节和经验教训，一层一层剥开。不管你是想做胚胎发育、类器官、生殖医学还是干细胞相关课题，这篇博文里讲的东西，都能让你少走几个月的弯路。

1. 别看标题长，核心就一句话：DNA甲基化重建是谁管的

1.1 为什么要盯住“早期胚胎发育”这个窗口期

先聊背景，不然采样、建库、数据解读都是空中楼阁。哺乳动物的早期胚胎发育，是指从受精卵到囊胚形成的这段关键期。这个阶段极其特殊：两个完全分化的配子（精子和卵母细胞）要擦掉自己的表观记忆，重新建立一套全能性的表观状态，然后进一步分化出内细胞团和滋养层。

这个过程里，DNA甲基化是变化最剧烈、也最难捉摸的表观修饰之一。刚受精时，父源基因组会经历快速的大规模去甲基化，母源基因组稍微稳一些，但也会缓慢地被动去甲基化。到了囊胚阶段，DNA甲基化水平会降到一个谷底，之后又要经历重新甲基化，为后续的谱系分化做准备。这一整套流程，学术界叫“表观遗传重编程”，其中重新甲基化这一步，是胚胎从“几乎裸奔”的状态，重新获得稳定基因表达调控网络的核心环节。

问题来了：那么大规模的DNA甲基化重建，是谁在背后指挥？靠什么机制把甲基化信号精准地放回该放的位置？如果这一步出错，胚胎基本废掉，发育停滞、流产、印记异常都跟它有关。所以这个窗口期，是研究表观遗传调控的黄金模型。

1.2 颉伟/卢绪坤/张宇团队做了什么

颉伟团队一直做早期胚胎的表观遗传研究，实验室的核心技术底子就是微量细胞级别的染色质和DNA修饰检测。卢绪坤和张宇在这个课题里分别承担了湿实验建库和数据分析的关键环节。他们这次把注意力放到了组蛋白H3第36位赖氨酸的二甲基化修饰上，也就是H3K36me2，研究的任务很明确：搞清楚H3K36me2在DNA甲基化重建过程中扮演什么角色。

说句实话，H3K36me2在圈内一直是个“干活多但名声小”的修饰。实验室里大家天天提H3K27me3、H3K4me3这些明星分子，H3K36me2的检测也一直没有特别高分辨率的全基因组图谱，因为它的分布太弥散了。但这次研究把H3K36me2从配角拉到了C位。三个组学技术同时上阵，目标就一个：证明H3K36me2是早期胚胎DNA甲基化重建的关键调控因子，并且把机制链条打通。看完这篇文章的技术路线，你会发现这不是一篇靠测序堆出来的数据论文，而是一套有严格假设、有正交验证、有功能实验的完整逻辑闭环。

2. 科学问题拆解：H3K36me2凭什么成为嫌疑对象

2.1 重新甲基化不是“一刀切”，而是“精准滴灌”

要理解这个研究的巧妙之处，必须先弄明白DNA甲基化重建的生物化学逻辑。细胞里的DNA甲基转移酶（主要是DNMT3A和DNMT3B）负责把甲基基团加到CpG二核苷酸上，但它们不是满天撒网，而是有选择性的——比如基因间区要高度甲基化，启动子区域特别是CpG岛要保持低甲基化，基因体内部要有适度的甲基化。

问题来了：DNMT3家族有催化活性，但它的靶向特异性并不强。它怎么知道该去哪些区域干活？靠的就是跟组蛋白修饰的“对话机制”。之前大量研究已经证明，H3K36me3（三甲基修饰）能够招募DNMT3B参与基因体区域的甲基化维持，H3K36me2则被认为是DNMT3A识别的重要信号——特别是DNMT3A的PWWP结构域能直接结合H3K36me2，从而被招募到基因组特定区域。

但在早期胚胎这个场景里，H3K36me2和DNA甲基化重建的关系一直缺乏系统性验证。难点在于：第一，胚胎细胞量极少，技术门槛高，很难做高分辨率的H3K36me2图谱；第二，H3K36me2在基因组上的分布太宽泛，数据分析难度大；第三，早期胚胎中DNMT3A和DNMT3L的表达水平是动态变化的，需要精确定量。

2.2 关键假说的形成：从“相关性”走向“因果性”

研究团队的高明之处在于，他们没有上来就盲目全基因组扫描找差异区域，而是先依托已有的生化证据和表达谱数据，形成一个明确可检验的假说：如果H3K36me2真的参与指导DNA甲基化重建，那么在重新甲基化即将启动的早期胚胎阶段，H3K36me2的分布应该先于DNA甲基化建立，两者在全基因组范围内应该有显著正相关；当H3K36me2被扰动时，DNA甲基化重建应该出现区域性失败。

这个假说里有明确的时空调度关系：先有时间上的先后，再有空间上的共定位，然后再有表型上的因果。为了验证这三层逻辑，WGBS、ChIP-seq、RNA-seq三个组学恰好各管一段：WGBS给全基因组甲基化状态做高清快照，ChIP-seq提供H3K36me2以及相关组蛋白修饰的位置信息，RNA-seq提供甲基转移酶和其他关键因子的表达动态。这三块拼图合在一起，才能回答一个完整的机制问题。

这里有个经验值得你学：**高水平组学课题，不是靠“能测就测”堆出来的，而是要先把生物学问题翻译成可拆解的技术需求。**拿到一个研究构想后，先想清楚每个组学数据是为了回答哪一层问题，再决定要不要测、测多深、覆盖什么样本。

3. 多组学技术方案设计与实施：WGBS、ChIP-seq、RNA-seq如何分工

3.1 WGBS：全基因组甲基化的“底图”

WGBS的核心原理很简单：用亚硫酸盐处理DNA，把未甲基化的胞嘧啶转变成尿嘧啶，甲基化的胞嘧啶不变。测序后对照参考基因组，就能逐个CpG位点判断甲基化状态。这是DNA甲基化检测的金标准，覆盖度能铺满整个基因组。

但在早期胚胎样本上跑WGBS，难处全在“微量”两个字上。取到的卵母细胞、受精卵、桑葚胚，里面是几十到几百个细胞，DNA总量通常在纳克级甚至皮克级。常规WGBS建库要几百纳克DNA起步，直接套用是行不通的。

这篇文章能跑下来，说明建库环节踩准了几个关键点：

• 首先是DNA抽提和亚硫酸盐转化的衔接，要尽量减少DNA在化学处理中的降解。亚硫酸盐反应本身就是剧烈的化学过程，会打断DNA链，微量DNA经不起折腾，必须优化反应时间、温度、纯化方式，把损失降到最低。
• 其次是建库策略。用末端修复加A、接头连接之后再亚硫酸处理，还是先处理再建库？不同策略对微量起始量的耐受差别很大。前者叫post-bisulfite，连接效率受甲基化状态影响小；后者叫pre-bisulfite，对起始量要求更高。针对几十个细胞级别，选择预扩增和优化连接酶体系是关键中的关键。

实际上，组学服务商和PI反复打磨微量WGBS流程，核心指标就是三个：转化率（一般要求≥99%）、比对率、有效CpG位点覆盖深度。这篇文章的数据能支撑全基因组级别的甲基化重建分析，说明WGBS覆盖做到了足够的深度（至少每个样本几百万到上千万个甲基化位点）。

3.2 ChIP-seq：在微观尺度上给H3K36me2“画地图”

ChIP-seq是用来定位组蛋白修饰的标准手段：用抗体把带有特定修饰的染色质片段“钓”出来，测序之后就知道修饰发生在基因组的哪些位置。

H3K36me2的ChIP-seq比常见的H3K27me3、H3K4me3要难做。因为它的信号在全基因组分布比较均匀弥散，信噪比天然就低，富集倍数没有“峰”那么鲜明。再加上早期胚胎染色质含量极低，每个细胞的组蛋白总量都少，抗体效率稍差就拉不出信号。

我见过不少实验室在跑H3K36me2 ChIP-seq一票否决的翻车案例，主要死因有三个：

1. 抗体批次不稳，有的批次特异性和亲和力不够；
2. 染色质断裂方式不合适，如果使用超声打断不充分，目标区域富集效率差；
3. 起始细胞量太低，qPCR质检就不过关。

颉伟团队在微量ChIP-seq上的积累是有体系的，包括对文库扩增循环数的精准控制、去背景噪声的测序策略、以及基于单管反应减少样本丢失的操作流程。另一个技术选择也值得注意：这种针对弥散型组蛋白修饰的检测，使用可能采用了CUT&Tag或类似方案来替代传统ChIP-seq，但文章标题里明确写了ChIP-seq，说明他们用常规工作流跑出了高质量信号，微量染色质免疫沉淀的效率优化确实有门道。对样本量在几千乃至几百细胞级别的组蛋白修饰图谱检测来说，抗体滴定、IgG对照、生物学重复这三样东西，一个都不能省。

3.3 RNA-seq：动态表达谱修正“因果链条”

WGBS建好了甲基化底图，ChIP-seq定位了H3K36me2位置，接下来就需要知道那些关键的甲基转移酶辅因子在哪个时间点表达。RNA-seq在这里起到的是“时间坐标”功能。

如果用伪时序的思路去描述胚胎早期发育，精子和卵母细胞是两个不同的起点，受精后基因表达经历剧烈的启动变化，合子基因组激活（ZGA）是其中最大的一个转折点。DNMT3A、DNMT3B、DNMT3L等因子的表达水平在不同阶段完全不同。RNA-seq做的是在外显子和转录本水平上定量这些基因的表达变化，帮助建立“什么时候有酶、哪里有信号、什么时候开始加甲基”的时序对应。

这里必须提醒一个常见误区：**RNA-seq在组学整合中不只是用来找差异表达基因的。更多时候，它是用来做因果链条中的“表达验证节点”。**你光看到H3K36me2和DNA甲基化空间相关还不够，还要看到关键酶的RNA表达和蛋白活性在时间上是否同步。如果酶还没表达，甲基化重建已经开始，那你的因果解释就有问题。这篇文章里RNA-seq与ChIP-seq时序数据互相对照，才敢下“H3K36me2调控DNA甲基化重建”这个强的结论。

4. 核心机制与关键发现：H3K36me2如何管住DNA甲基化重建

4.1 从全基因组相关性中抓出“空间共定位”

这篇文章最核心的技术性证据，来自WGBS和ChIP-seq的整合分析。简单来说，就是把基因组按H3K36me2的信号强弱分成不同区段，然后看这些区段对应的DNA甲基化水平。如果H3K36me2确实指导了DNA甲基化重建，那么H3K36me2信号强的区域，DNA甲基化水平应当更高；反过来说，那些缺乏H3K36me2信号的区域，甲基化重建会延迟或失败。

结果也确实是这个方向：早期胚胎的重新甲基化区域，与H3K36me2的分布高度重叠。特别是基因间区、内含子区域这些原本H3K36me2密度比较高的区域，DNA甲基化恢复得比较早、也比较充分。而启动子区、CpG岛等H3K36me2标记较弱的位置，即使在高甲基化的背景下也维持着相对低的甲基化水平。

这种“空间共定位”是机制研究的基础，但还不够有力。所以研究团队继续深入了一层——看H3K36me2与DNA甲基化重建在发育进程中的动态关系。结果发现，H3K36me2信号建立的时间点早于DNA甲基化的恢复，存在清晰的时间先后顺序。这个先后关系是从“相关”迈向“因果”的重要一步。

4.2 调控链条：组蛋白修饰与甲基转移酶之间的“接头”

光有时间和空间的相关性还不够，要做实机制，必须把中间的分子桥梁找出来。这就是跨组学数据联动的价值所在。已知DNMT3A有PWWP结构域，能够结合H3K36me2。研究进一步支持了这样一个模型：

在卵母细胞到早期胚胎的转变过程中，H3K36me2如同铺设在地上的“标记线”，甲基转移酶DNMT3A/3B结合这些标记，被引导到需要重建DNA甲基化的区域。在H3K36me2缺失或者信号很弱的区域，甲基转移酶失去锚点，重新甲基化无法有效进行。

还有一个细节值得讲：研究里可能还对比了H3K36me2和H3K36me3的不同角色。三甲基化修饰H3K36me3在基因体区域富集，与转录延伸紧密相关；而H3K36me2更多覆盖基因间区和散在区域。两种修饰虽然都在H3K36上，但对甲基转移酶招募的偏好可能不一样。把两者区分开，才能避免把“H3K36me3维持基因体甲基化”的老结论直接挪用到胚胎重编程场景里。

4.3 功能验证：扰动H3K36me2后，DNA甲基化重建是否“翻车”

理想中的机制研究，最后一步一定要做功能验证。不是只在野生型样本里说“相关性强”，而是要把H3K36me2这个变量扰动掉，观察DNA甲基化重建是否出现可预测的缺陷。

这种扰动在胚胎里比细胞系复杂得多。实验室常见的做法是在卵母细胞或胚胎中干扰H3K36me2的催化酶（例如SETD2负责三甲基化，而NSD1/NSD2/NSD3负责二甲基化），或者直接干扰DNMT3A的表达。如果H3K36me2调控DNA甲基化重建的假说成立，那么干扰H3K36me2生成后，应该看到特定区域的DNA甲基化恢复受阻，基因表达也随之出现异常。文章中做的功能实验与测序数据形成闭环，不是只停留在“我们观察到了很多数据”这个层面。

这一点恰恰是很多组学文章被审稿人挑刺的重灾区。如果你写了“X调控Y的发生”，却不做任何扰动实验，审稿人两个字就给你打发了：关联。关联性研究和高水平机制研究的差别，在于有没有在因果层级上给出失活/过表达/回复证据。

5. 微量样本组学的实操要点与避坑指南

5.1 从样本收集到文库构建，每一个环节都在“掉数据”

我见过太多课题组，实验设计写得漂亮，最后死在样本量不足和样本质量不佳上。早期胚胎样本的性质决定了，你不可能像拿细胞系那样反复补样。一颗卵母细胞用掉了就没有了，一组囊胚用完就得等下批小鼠。这逼着你必须把每一步的回收率都做到极致。

从实验操作的角度，我列几个最容易被低估的损失点：

• 数次管壁吸附：微量DNA在低吸附管里也会被吸附，更别说普通EP管。DNA溶解液最好用加了Tween-20的缓冲液，能明显降低吸附损失。
• 纯化柱的洗脱体积：别贪图“浓缩”，每次洗脱用小体积但分两次操作，可以显著提高回收量，上柱前也不必盲目追求低洗脱体积，因为吸附柱自身就有死体积。
• 热循环管的封口性：微量建库的PCR体系在高温下容易蒸发，用薄膜封板或者热盖到最高档，别省钱。
• 建库扩增循环数：微量ChIP-seq文库必须比常规样本多扩几个循环。但循环数越高，重复序列的PCR偏向和冗余也越高。要看清楚文库的复杂度曲线，达到足够产物后立刻停。

5.2 ChIP-seq抗体和IgG对照的那些事

H3K36me2的ChIP-seq想做成功，抗体是生死线。市售抗体品牌很多，同一个靶标的抗体在不同批次间表现差异极大。我建议在建库前先做小规模dot blot或者Western blot测试，最好再跑一个已知富集区域的qPCR验证，不要看了官网说“ChIP级别”就直接上。很多所谓ChIP级别的抗体，实际上在微量样本上效率并不稳定。

再就是IgG对照。普通细胞样本的ChIP-seq可能用input做背景扣除就可以了，但微量样本不行，非特异背景的比例会被放大，不做IgG对照会导致大量假阳性峰。我们做微量ChIP-seq时，IgG对照的起始量要和IP一样，就哪怕细胞数量极少，也要硬挤出对照用的量。别嫌浪费，数据后期清洗时你会感谢这份对照。

5.3 WGBS的转化率与覆盖深度的平衡

WGBS分析甲基化水平时，最重要的一项质控就是亚硫酸盐转化率。转化率低，未甲基化胞嘧啶残留，会让所有甲基化水平虚高。理论上要求达到99%以上。怎么估？最常见的方式是加Lambda DNA等未甲基化的外部对照，或者利用非CpG上的CHH/CHG甲基化信号估算，因为哺乳动物非CpG甲基化水平本身极低，信号干净。

再强调一下覆盖深度：全基因组甲基化分析不是越深越好，太深费用爆炸，太浅又看不出差异。做胚胎发育这种课题，关键是覆盖率的均匀性，而不是单个位点的极端深度。建议至少满足每个样本全基因组有效CpG位点覆盖在10x以上，这样在分析中才能有信心地对启动子、基因体、重复区域的甲基化做统计推断。

6. 多组学整合分析：不是堆数据，是拼拼图

6.1 数据怎么“对齐”才不会张冠李戴

拿到WGBS、ChIP-seq、RNA-seq三套数据后，第一件事是把它们放在同一个坐标系里比对。这里的麻烦有两个：

1. 基因组版本必须统一，都是mm10就都mm10，别一个用mm10一个用GRCm38的小差异版本，细节能坑死人；
2. 峰/区间定义的统一，H3K36me2 ChIP-seq的区间不能用最简单的peak calling，因为弥散型修饰没有清晰峰形，需要采用“宽峰”模式或基于信号分箱（bin）的思路。

我更推荐的做法是：定义一套统一的基因组分箱（比如把基因组切成1kb的窗口），然后分别统计每个窗口的甲基化水平（WGBS）、H3K36me2信号（ChIP-seq）、基因表达量（RNA-seq）。然后在窗口层级做相关性分析和共定位分析。这种“bin-based”的整合方式可以很好地规避不同数据类型分辨率不匹配的问题。

6.2 差异区域与分层聚类怎么看

一个比较实用的分析路径是：先用WGBS数据找出重新甲基化阶段“恢复不足”或者“恢复过快”的差异甲基化区域（DMR），然后叠加上H3K36me2的ChIP-seq信号，评估这些DMR区域是被H3K36me2覆盖还是缺失。再往下，把基因注释挂上去，看这些DMR关联的基因表达是否受影响。

这里有个小技巧：在分析前，把全基因组按H3K36me2信号高低分成几个等级（比如高、中、低、无），做一个sankey图或者热图，直观展示H3K36me2等级与DNA甲基化水平之间的对应关系。这种图既是文章视觉的高光，也是你和合作者讨论时要绝对能拿出手的“定盘星”。

6.3 分辨率降维：从全基因组粗粒度到位点级别刻画的策略

甲基化和组蛋白修饰的分析分辨率天然不同：WGBS可以做单碱基分辨率，ChIP-seq则是富集区间级别。做机制研究时，你可以先在全基因组层面提出“共定位”结论，然后把焦点缩小到一个DNA甲基化重建很典型的区域，比如某个发育调控基因的基因间增强子区域，做高分辨率展示。这种从全局到局部的分析策略，比一上来就细抠某个位点更容易说服人。

7. 常见问题与排查技巧实录

7.1 问题速查表

7.2 实操心得：组学项目成功的几根“定海神针”

第一，一定要在项目启动前做一次全流程“试跑”。哪怕是拿几万个细胞的细胞系先模拟一遍微量样本的建库流程，也比直接拿胚胎练手强。试跑能暴露所有环节的问题：抗体能不能富集、WGBS转化率是否达标、RNA建库会不会降解。胚胎样本一旦消耗就没法补救，预实验的投入绝对值得。

第二，生物学重复的数量不能妥协。早期胚胎样本个体差异很大，不同小鼠、不同批次处理的胚胎，表观状态可能差很远。如果你只用1个或2个重复，审稿人大概率会质疑统计效力。理想情况下每个时间点至少3个生物学重复，而且尽量保持性别、品系一致。

第三，数据分析和湿实验要一直保持对话。不要湿实验跑完就把数据扔给生信人，然后等结果。最理想的状态是每个时间点的数据一出来，马上做快速QC，把质检报告拿给湿实验负责人看，共同判断是否需要补样。多组学项目最大的坑，是等到所有数据测完才发现有一个时间点的数据不合格，这时候已经晚了。

第四，组蛋白修饰和DNA甲基化的动态窗口，用什么时间点多长间隔，要非常讲究。胚胎发育是个连续过程，但测序只能取离散时间点。如果时间点隔得太远，会错过甲基化重建启动的关键窗口，数据就算测了也会发现变化早就发生完了。这个时间点的选择逻辑，建议和文献数据对照着定，别拍脑门。

8. 写在最后的经验补充

这篇文章的解读，我不想做成论文复述，更愿意把它当成一套方法论的案例。你在复现和借鉴的时候，有几个细节是我个人特别想再叮嘱一次的。

一个是弥散型组蛋白修饰的ChIP-seq，永远不要只用默认的peak caller参数。H3K36me2不是H3K4me3那种峰型清晰、信号锐利的分型修饰，你拿MACS2默认参数跑，大概率得到一堆碎片化的伪峰。建议做信号窗口化分析，把全基因组切成固定大小的bin，统计每个bin的归一化信号强度，再去和甲基化数据做关联。这种方法牺牲了一点点单点精度，但保住了全局视角的可靠性。

另一个是多组学整合分析前，先把批次效应和文库复杂度差异处理好。微量样本一般扩增循环数较高，文库复杂度低，在比对后要留意重复去除的比例。不同时间点的样本之间，测序深度可能差几倍，一定要先做归一化（比如CPM、RPKM，或者更严格的quantile normalization），再做跨样本比较。我见过太多人拿着没深度匹配的数据直接画热图，画出来的“差异”其实是测序深度差异。

最后，好的表观遗传研究永远是“数据驱动+假设驱动”双轮走。WGBS和ChIP-seq联合给出的相关图谱只是起点，真正让审稿人买账的，是把这条通路上的关键因子的表达、定位、功能扰动，全部串联上。各位做同类课题的时候，别急着堆样本跑数据，先花两周时间把实验设计的逻辑图谱画清楚，多想想每个数据块在机制故事里承担什么角色。磨刀不误砍柴工，这话放在大文章里，永远是真理。