大概两年前的这个时候我手里攥着一批好不容易分选出来的小鼠T细胞亚群总数不到3万却要做转录因子结合谱。按ChIP-seq的常规要求这个数量连最乐观的实验方案都撑不起来。后来查文献看到了CUTRUNCleavage Under Targets and Release Using Nuclease靶向切割与释放不到两周就把这批样本跑完了出来数据的信噪比让我想把之前ChIP-seq的结果全部重做一遍。这不是夸张CUTRUN这几年在表观遗传领域的扩散速度已经说明了一切——它把蛋白质-DNA互作研究的细胞量门槛从百万级拉低到了万级背景还低一个数量级整个流程从两天缩短到一天。这篇文章我不打算复述说明书想把原理、实操、生信和踩坑经验串起来讲一遍给想上手或者正在纠结技术选型的人一个完整的参考。1. ChIP-seq的老问题CUTRUN到底在解决什么1.1 三个绕不开的痛点先回忆一下ChIP-seq的标准流程甲醛交联固定细胞超声或者酶切把染色质随机打断成片段然后用目标蛋白的抗体做免疫沉淀洗掉非特异结合的片段解交联、纯化DNA建库测序。这套流程发展了十几年非常成熟但它的底层逻辑有先天缺陷。第一个痛点是细胞量。常规ChIP-seq需要100万到1000万个起始细胞因为超声打断和免疫沉淀都是损耗型工艺每一步都在丢东西。哪怕用优化过的微量ChIP方案低于10万个细胞做出来的文库也经常一片空白。可真实课题里流式分选出来的原代细胞、类器官里的稀有亚群、临床穿刺样本动辄只有几万个甚至几千个细胞这时候ChIP-seq基本等于宣判实验做不了。第二个痛点是背景。超声打断染色体是不分青红皂白的所有区域的DNA都被随机打碎。免疫沉淀时抗体抓到的是目标蛋白-DNA复合物但同时也有大量不带目标蛋白的裸露DNA被磁珠或抗体非特异吸附。所以ChIP-seq数据里很大一部分reads是背景噪音只能靠测序深度去堆信号这也是为什么ChIP-seq通常要求50M甚至更多reads的原因。第三个痛点比较隐形是超声偏好性。超声波打断开放染色质区域时会更容易导致某些区域被过度打断某些异染色质区域几乎打不动这种系统偏差会直接污染后续的峰识别和差异比较。1.2 从先打断再富集到先标记再切割CUTRUN的思路完全反过来了。它不再先打断染色质而是让抗体在完整的细胞核里找到目标蛋白然后把一个能切割DNA的酶拴在抗体上最后加钙离子激活这个酶让它只切抗体所在位置的DNA片段。这里的核心工具是一个融合蛋白pA-MNaseprotein A-MNase。protein A来自金黄色葡萄球菌能特异性结合大多数抗体的Fc段MNase是微球菌核酸酶能在钙离子存在时切割DNA。你只需要把抗体孵育到细胞核里它结合到目标蛋白上再加入pA-MNase它会像磁铁一样贴到已经结合在染色质上的抗体屁股后面这时候清洗掉多余未结合的pA-MNase加入钙离子触发切割。因为游离的MNase已经被洗掉了真正留在体系里能切DNA的酶几乎只存在于抗体定位的那个局部区域。切下来的DNA片段由于比较短亚核小体级别或单核小体级别会从细胞核里扩散出来进入上清。离心收集上清纯化DNA建库测序。整个过程没有超声没有免疫沉淀没有交联大多数情况下这就是它背景低、样本量要求低的结构性原因。1.3 为什么说背景能低一个数量级我见过很多第一次接触CUTRUN的人以为它只是把ChIP-seq里的超声换成了酶切其实差别比这深刻得多。关键在于切割的靶向性。超声是全局性的、物理性的随机剪切免疫沉淀的富集效率决定了信噪比而CUTRUN的切割是抗体引导的、化学计量学层面的局部事件。没有抗体结合的区域即使MNase偶然扩散过去也因为没有稳定结合而很快被洗掉真正有效的切割只围绕抗体-靶蛋白-染色质复合物发生。结果就是非靶区域的DNA几乎原封不动留在核内最后进入上清的一般只占总DNA的1%-2%而这1%-2%里目标区域的占比又极高。打个比方ChIP-seq是把整栋楼炸碎然后趴在地上筛想找的砖头CUTRUN是用定点爆破拆除目标房间直接用车拉走。前者需要很多楼细胞量、很多筛子清洗和IP步骤、很高的测序深度后者只需要少量样本就能拿到干净得多的产物。2. CUTRUN实操流程全拆解2.1 样品准备与ConA磁珠固定CUTRUN的第一步是制备细胞悬液和细胞核。贴壁细胞用胰酶或者细胞刮刀收下来离心清洗悬浮细胞相对简单直接计数取用组织样本则需要先消化成单细胞悬液或者提取细胞核。关键要求是细胞状态好、碎片少、死细胞比例低——死细胞会把基因组DNA降解成大分子粘稠物严重影响后续的磁珠结合和切割。接下来是ConA刀豆蛋白A磁珠结合。ConA识别细胞核膜上的甘露糖残基把完整的细胞核粘到磁珠上。这里有个技术细节ConA结合需要钙离子和锰离子的辅助所以Binding Buffer里必须加1 mM的CaCl₂和1 mM的MnCl₂否则结合效率会大幅下降。结合时间一般10-15分钟室温即可。结合完成后的清洗操作要格外轻柔。细胞核靠着糖蛋白-凝集素的非共价相互作用挂在磁珠上剧烈吹打会把核弄破一旦核膜破裂染色质外泄整个实验就废了。我的习惯是磁分离后用枪头缓慢吸走上清不要吹悬磁珠而是加溶液后轻弹管底让沉淀松动。2.2 抗体孵育与pA-MNase结合决定成败的两步CUTRUN的抗体孵育和ChIP-seq有个重要差异它基本不做交联识别的是天然构象的表位。这意味着ChIP级别的抗体不一定适合CUTRUN而一些在免疫荧光、流式或native ChIP中表现好的抗体反而可能更合适。同一支抗体在这两种实验里表现可能差异巨大所以刚上手时务必做抗体预实验。抗体用量一般参考说明书常规体系5×10⁵细胞常用0.5-2 µg/反应。孵育条件有两种4°C过夜或室温2小时。室温孵育信号更高但非特异背景也会略上升4°C过夜则更温和适合对温度敏感的转录因子。我建议第一次跑可以从室温1-2小时起步用IgG对照同时做。孵育结束洗掉未结合抗体后加入pA-MNase1:100到1:200稀释具体浓度需按批次优化。这里特别提醒很多商业化的pA-MNase推荐在4°C结合1小时千万别在室温放太久不然酶还没被激活就可能降解核内的DNA。另外如果你的一抗是鼠源单抗你会发现protein A对小鼠IgG的亲和力不太够这时候要用pAG-MNaseprotein A/G融合它同时结合兔源和人源、鼠源的各类IgG亚型。注意pA-MNase的浓度不是越大越好。酶多了即使在0°C也会产生非特异切割酶少了目标区域切不充分。首次上手最好做一个1:50、1:100、1:200、1:400的梯度用qPCR检已知靶位点的富集倍数来定。2.3 钙离子激活切割与片段释放pA-MNase结合并清洗干净后加入CaCl₂到终浓度2 mM在冰上进行30-60分钟的切割。MNase是一个钙依赖核酸酶没有钙离子它完全不工作加了钙之后立刻激活。冰上处理的目的是让酶和染色质都处于低温低扩散状态最大程度限制切割发生的位置减少背景。切割结束后加入含EGTA的Stop Buffer终止反应。EGTA是钙离子的螯合剂能瞬间把游离钙离子络合掉让MNase失活。在Stop Buffer这一步骤要加入E. coli染色质或者E. coli细胞作为spike-in内参——这是CUTRUN体系里一个非常重要的实验设计后面数据分析部分我会专门展开。随后将反应管离心600 g5分钟或用磁力架分离磁珠上清就是释放出来的DNA片段。这里有个容易被忽略的点部分片段可能还停留在核内没完全释放可以选择加一次释放缓冲液再离心把产量提高20%左右。上清转移到新的低吸附离心管进入DNA纯化。2.4 从DNA到文库小体系操作的细节CUTRUN产物极其微量一个转录因子实验通常只能拿到几纳克到几十纳克的DNA。这个数量下任何管壁吸附都可能是致命的。我第一次做的时候用的普通EP管最后纯化洗脱下来跑Qubit显示只有0.3 ng/µL当时一度怀疑整个流程出了问题。后来换成DNA LoBind管同批次样品直接翻了三倍。DNA纯化可以用酚氯仿抽提加乙醇沉淀也可以用专门的微量DNA纯化柱。如果用柱子务必选用对短片段150 bp回收效率高的产品很多通用胶回收柱对小片段的结合效率很低。沉淀时加GlycoBlue糖原类载体能显著提高微量DNA的回收率。文库构建方面CUTRUN产物不需要机械打断直接做末端修复、加A尾、连接接头、PCR扩增。由于起始量极低PCR循环数通常在12-14个循环这里要平衡两个问题循环数太少文库量不够上机循环数太多PCR重复增多、复杂度下降。遇到产量特别低的样本宁可接受稍低文库产量也不要盲目加到16个循环以上否则后续去重后有效reads会少得可怜。片段大小筛选也要注意做组蛋白修饰时目标产物是单核小体约150 bp及少量二核小体做转录因子时目标产物主要是亚核小体片段30-120 bp可以优先回收小片段。切胶或双轮Ampure XP磁珠都可以实现但我更推荐后者操作稳定且损耗小。3. CUTRUN、CUTTag与ChIP-seq三种路线到底怎么选3.1 核心差异对照CUTRUN火了之后紧接着又出现了CUTTagCleavage Under Targets and Tagmentation。很多人把它们混为一谈其实分工并不一样。三个技术可以放在一张表里看维度ChIP-seqCUTRUNCUTTag切割工具超声物理打断pA-MNase定点酶切pA-Tn5转座酶起始细胞量10⁶-10⁷10⁴-10⁶可低至千级10²-10⁵可单细胞交联通常需要甲醛交联多数不需要通常不需要建库方式片段化后连接接头纯化DNA后连接接头转座酶直接插入接头背景信号较高低较低实验周期2-3天1-2天1-2天适用抗体已广泛验证需重新验证需重新验证成本构成测序深度高、耗材常规酶/试剂盒较贵、测序便宜试剂盒贵、测序便宜CUTTag的核心是把Tn5转座酶融合到protein A上抗体定位之后转座酶直接在目标区域把测序接头标进染色质片段化和建库一步完成。它的细胞量门槛比CUTRUN更低做到单细胞也是它在拓展而且不需要单独纯化DNA操作更简单。3.2 我实际使用中的选型建议先给结论如果你的细胞量在5万以上目标是转录因子或组蛋白修饰且希望拿到比较纯粹的蛋白-DNA互作信号CUTRUN是首选。它的切割机制决定了信号更接近真实的结合位点分辨率高峰型尖锐背景极低。如果细胞量在1万以下特别是几百到几千个细胞CUTTag更合适。Tn5转座酶的效率远高于MNase在极低起始量下建库成功率更高。如果你的目标最终是单细胞水平那基本只有CUTTag一个选项。如果样本是临床FFPE组织或者你的抗体库全部是已发表的ChIP-seq抗体且没有时间去重新验证那ChIP-seq仍然是相对稳妥的路径。CUTRUN虽然是升级版但升级意味着规则变了抗体、内参、生信都要重新适配这个成本要提前算进去。还有一个容易被忽视的维度如果你后续要做大量的组蛋白修饰全基因组图谱并且要和自己之前采集的ChIP-seq公共数据做大规模比较跨技术比较的可信度是个问题。CUTRUN的信号分布特征和ChIP-seq不一样差异分析时可能会引入批次效应这种情况除非所有新数据都用同一种技术补齐否则谨慎更换方法。4. 数据分析的关键环节比对、归一化、找峰4.1 比对管线的参数优化很多人觉得CUTRUN的生信分析可以照搬ChIP-seq流程实操后就会发现不是那么回事。CUTRUN的插入片段集中在30-150 bp双端测序时两条reads的比对特征和普通ChIP-seq很不一样。推荐管线是质量评估FastQC→ 去接头cutadapt或trim_galore→ 比对bowtie2→ 过滤低质量比对和重复 → 生成bigWig可视化 → 寻找peaks。bowtie2的参数建议这样设置bowtie2 -p 8 --end-to-end --very-sensitive --no-unal \ --no-mixed --no-discordant --phred33 -I 10 -X 700 \ -x /ref/bowtie2/hg38 -1 clean_R1.fq -2 clean_R2.fq -S sample.sam--no-mixed和--no-discordant是为了过滤掉非一致性比对因为CUTRUN片段短理论上双端reads应该都能唯一定位到同一区域。比对完成后用-q 10过滤掉低质量比对再去除PCR重复。加性别决定染色体区域去除可选但至少要把线粒体reads去掉它们经常占比非常高。4.2 Spike-in归一化跨样本比较的基准CUTRUN实验里不同样本之间的磁珠结合效率、清洗损失、酶切效率都会带来系统差异。如果不做归一化两个样本的绝对信号高度无法直接比较。E. coli spike-in就是用来解决这个问题的。操作层面你在Stop Buffer里加入恒定数量的E. coli细胞或E. coli染色质它们会经历和你的样本相同的后续所有步骤。E. coli基因组上没有你目标蛋白的结合位点理论上不会产生特异性信号但它们的DNA也会被环境中残余的MNase非特异切割、释放、纯化、建库。所以样本间spike-in reads的比例就反映了整个实验流程的效率差异。分析时先把reads分别比对到目标基因组如hg38和E. coli基因组如K12 MG1655统计两个比对文件里的有效reads数。归一化系数可以用# 计算样本spike-in比例以最小spike-in reads样本为基准 scale_factor min_ecoli_reads / sample_ecoli_reads然后用这个系数对bam文件缩放生成可视化trackbamCoverage -b sample.sorted.bam -o sample.scaled.bw \ --scaleFactor 1.32 --binSize 10 --normalizeUsing nonespike-in reads一般占总reads的1%-5%比较理想。如果某个样本的E. coli比例特别高说明酶切过度或者清洗不充分如果特别低说明整个流程效率差这个样本的信号可信度就要打问号。4.3 SEACR与MACS2的取舍找峰的软件选择是CUTRUN数据分析里争论最多的点。Henikoff实验室专门开发了SEACRSparse Enrichment Analysis for CUTRUN核心思路是用全基因组信号分布来替代背景模型不依赖IgG对照。MACS2虽然是ChIP-seq时代的标准工具但在CUTRUN上直接使用会产生很多假峰因为两者的负二项分布假设并不匹配。我的经验是转录因子窄峰建议优先试SEACR的stringent模式组蛋白修饰宽峰可以试relaxed模式。SEACR需要输入bedgraph文件命令如下# 生成bedgraph bedtools genomecov -bg -scale 0.25 -ibam sample.sorted.bam sample.bedgraph # SEACR找峰stringent模式 SEACR_1.3.sh sample.bedgraph IgG.bedgraph stringent sample_stringent # SEACR找峰relaxed模式适合宽峰 SEACR_1.3.sh sample.bedgraph IgG.bedgraph relaxed sample_relaxed如果要用MACS2需要调整参数适配CUTRUN的短片段特征macs2 callpeak -t sample.bam -c IgG.bam -f BAMPE -g hs \ --keep-dup all -q 0.05 --nomodel --shift -100 --extsize 200 \ -n sample实际使用中我会把SEACR和MACS2都跑一遍然后在IGV里目视检查一组已知位点。两个工具结果重合的部分基本可信只有其中一个检出的峰需要重点怀疑。质量评估还要看reads in peaks比例FRiPCUTRUN做得好时这个比例能到20%-60%远高于ChIP-seq的5%-10%。如果FRiP低于5%先别急着找峰回头排查抗体和酶切环节。5. 我踩过的坑以及一份排查清单5.1 高背景问题往往出在酶和温度实验室接CUTRUN的第一个月某个组蛋白修饰样本跑出来的reads在基因组上均匀分布完全看不到峰。当时第一反应是抗体不行结果排查一圈发现是pA-MNase浓度过高加室温孵育时间太长导致非特异切割全面开花。这件事给我的教训是CUTRUN的低背景不是自动获得的它的先决条件是酶被严格控制在抗体标记的局部。压背景有四个关键旋钮pA-MNase浓度宁低勿高、CaCl₂激活温度必须冰上、激活时间不超过60分钟、清洗次数抗体和酶结合后的清洗不能省。出现高背景时先检查这几个因素甚至可以把清洗次数从两次加到三次背景会显著下降。5.2 文库产量极低可能死于管壁和洗脱有一次给一个只有1万个细胞的稀有样本做CUTRUN最后Qubit检不出DNA浓度。当时浪费了一整个辛苦分选的样本。后来复盘发现两个致命问题一是用了普通离心管微量DNA被管壁吸附二是DNA纯化时用乙醇沉淀但温度没达到-20°C沉淀时间也不够。从那以后我立了几条规矩全程用LoBind管乙醇沉淀至少-20°C过夜并加GlycoBlue纯化完上机前不反复冻融文库PCR循环数严格控制在14以内如果产量不够宁可选更灵敏的建库试剂盒比如NEBNext Ultra II的低起始量方案也不盲目增加扩增循环。5.3 抗体不兼容CUTRUN有自己的抗体脾气CUTRUN对抗体的要求是天然构象下能识别目标蛋白。一支说明书标注ChIP Grade的抗体可能在甲醛交联后暴露的表位结合良好但在未交联的天然染色质上完全失效反过来有些IF级别抗体反而在CUTRUN里出奇地好。我建议正式实验前做一次抗体小筛用你的目标细胞和一支已知阳性抗体比如H3K4me3或H3K27ac这种在CUTRUN中很稳定的抗体做平行测试同时用IgG做对照跑一次低深度测序或至少qPCR验证靶位点富集。如果阳性抗体的FRiP明显高、IgG背景低说明整个体系是通的剩下的就是换目标抗体了。这一步看起来很费时间实际上能帮你少走两个月的弯路。5.4 质控指标速查表现象可能原因处理建议FRiP 5%抗体无效/浓度低、酶量不足、细胞状态差换抗体、做抗体梯度、检查细胞活力Spike-in比例过高15%清洗不充分、酶切过度增加清洗次数、降低pA-MNase浓度Spike-in比例过低0.5%流程损失大、细胞量太少用LoBind管、延长沉淀时间文库无PCR产物微量DNA丢失、接头连接失败加GlycoBlue、检查接头连接体系Peak集中在少数位点PCR重复过多、测序深度不足增加起始量或循环数提高reads数组蛋白修饰信号杂乱酶切时间过长、钙离子过量缩短激活时间、检查Stop Buffer浓度6. 上手CUTRUN之前值得想清楚的几件事6.1 你的样本类型是否适配不是所有样本都适合直接套用标准CUTRUN流程。植物细胞有细胞壁需要先制备原生质体或纯化细胞核临床组织样本需要先做单细胞化或核提取组织消化过程中的细胞死亡会直接影响ConA结合和后续切割冻存细胞原则上能做但冻存后细胞膜通透性和核膜状态都会变化建议先用活细胞做一版完整流程熟悉体系再拿冻存样本优化。另外CUTRUN的清洗步骤依赖磁珠高脂含量或黑色素含量高的样本会干扰磁分离和后续酶活这类样本光做细胞悬液准备可能就要额外花一周调试要做好心理准备。6.2 抗体验证的优先级如果让我用一句话总结CUTRUN最值得投入的地方我会说把70%的优化精力放在抗体上。实验流程是标准化的酶是商业化的数据分析管线也有成熟的工具唯独抗体是个变量。一支好抗体的CUTRUN数据和一支差抗体的数据差距比任何参数调整带来的差异都大。我验证抗体的标准流程是先用Western blot确认特异性再看是否有已发表的CUTRUN/CUTTag数据使用记录然后做一个起点为2 µg/反应、逐步降低浓度的测试最后在IGV里目视检查已知阳性位点和阴性区域。三个环节都过了这支抗体才真正入库。6.3 预算与外包的选择CUTRUN的试剂盒比如CST的CUTRUN Kit或EpiCypher的CUTANA系列价格不低单样本成本通常几百到上千元但它省掉的测序费是惊人的——CUTRUN单样本一般测10-20M reads就够了只有ChIP-seq的一半甚至三分之一。如果你的实验室没有现成的实验体系又急着出结果外包给有成熟CUTRUN平台的测序服务商也是一个务实选择。前提是确认服务商会提供spike-in归一化、SEACR分析这类CUTRUN专属的生信流程而不是拿ChIP-seq的流程糊弄你。最后再分享一个小技巧CUTRUN上机前可以用qPCR快速预检几个已知阳性位点和阴性位点的富集倍数。阳性位点比阴性位点高几十倍以上再送测序基本稳了。这个小步骤只花两个小时却能避免烧掉一次几千块的测序试错成本。实验做得越多越明白表观遗传学领域真正稀缺的不是花哨的技术而是能稳定复现的结果。
大概两年前的这个时候我手里攥着一批好不容易分选出来的小鼠T细胞亚群总数不到3万却要做转录因子结合谱。按ChIP-seq的常规要求这个数量连最乐观的实验方案都撑不起来。后来查文献看到了CUTRUNCleavage Under Targets and Release Using Nuclease靶向切割与释放不到两周就把这批样本跑完了出来数据的信噪比让我想把之前ChIP-seq的结果全部重做一遍。这不是夸张CUTRUN这几年在表观遗传领域的扩散速度已经说明了一切——它把蛋白质-DNA互作研究的细胞量门槛从百万级拉低到了万级背景还低一个数量级整个流程从两天缩短到一天。这篇文章我不打算复述说明书想把原理、实操、生信和踩坑经验串起来讲一遍给想上手或者正在纠结技术选型的人一个完整的参考。1. ChIP-seq的老问题CUTRUN到底在解决什么1.1 三个绕不开的痛点先回忆一下ChIP-seq的标准流程甲醛交联固定细胞超声或者酶切把染色质随机打断成片段然后用目标蛋白的抗体做免疫沉淀洗掉非特异结合的片段解交联、纯化DNA建库测序。这套流程发展了十几年非常成熟但它的底层逻辑有先天缺陷。第一个痛点是细胞量。常规ChIP-seq需要100万到1000万个起始细胞因为超声打断和免疫沉淀都是损耗型工艺每一步都在丢东西。哪怕用优化过的微量ChIP方案低于10万个细胞做出来的文库也经常一片空白。可真实课题里流式分选出来的原代细胞、类器官里的稀有亚群、临床穿刺样本动辄只有几万个甚至几千个细胞这时候ChIP-seq基本等于宣判实验做不了。第二个痛点是背景。超声打断染色体是不分青红皂白的所有区域的DNA都被随机打碎。免疫沉淀时抗体抓到的是目标蛋白-DNA复合物但同时也有大量不带目标蛋白的裸露DNA被磁珠或抗体非特异吸附。所以ChIP-seq数据里很大一部分reads是背景噪音只能靠测序深度去堆信号这也是为什么ChIP-seq通常要求50M甚至更多reads的原因。第三个痛点比较隐形是超声偏好性。超声波打断开放染色质区域时会更容易导致某些区域被过度打断某些异染色质区域几乎打不动这种系统偏差会直接污染后续的峰识别和差异比较。1.2 从先打断再富集到先标记再切割CUTRUN的思路完全反过来了。它不再先打断染色质而是让抗体在完整的细胞核里找到目标蛋白然后把一个能切割DNA的酶拴在抗体上最后加钙离子激活这个酶让它只切抗体所在位置的DNA片段。这里的核心工具是一个融合蛋白pA-MNaseprotein A-MNase。protein A来自金黄色葡萄球菌能特异性结合大多数抗体的Fc段MNase是微球菌核酸酶能在钙离子存在时切割DNA。你只需要把抗体孵育到细胞核里它结合到目标蛋白上再加入pA-MNase它会像磁铁一样贴到已经结合在染色质上的抗体屁股后面这时候清洗掉多余未结合的pA-MNase加入钙离子触发切割。因为游离的MNase已经被洗掉了真正留在体系里能切DNA的酶几乎只存在于抗体定位的那个局部区域。切下来的DNA片段由于比较短亚核小体级别或单核小体级别会从细胞核里扩散出来进入上清。离心收集上清纯化DNA建库测序。整个过程没有超声没有免疫沉淀没有交联大多数情况下这就是它背景低、样本量要求低的结构性原因。1.3 为什么说背景能低一个数量级我见过很多第一次接触CUTRUN的人以为它只是把ChIP-seq里的超声换成了酶切其实差别比这深刻得多。关键在于切割的靶向性。超声是全局性的、物理性的随机剪切免疫沉淀的富集效率决定了信噪比而CUTRUN的切割是抗体引导的、化学计量学层面的局部事件。没有抗体结合的区域即使MNase偶然扩散过去也因为没有稳定结合而很快被洗掉真正有效的切割只围绕抗体-靶蛋白-染色质复合物发生。结果就是非靶区域的DNA几乎原封不动留在核内最后进入上清的一般只占总DNA的1%-2%而这1%-2%里目标区域的占比又极高。打个比方ChIP-seq是把整栋楼炸碎然后趴在地上筛想找的砖头CUTRUN是用定点爆破拆除目标房间直接用车拉走。前者需要很多楼细胞量、很多筛子清洗和IP步骤、很高的测序深度后者只需要少量样本就能拿到干净得多的产物。2. CUTRUN实操流程全拆解2.1 样品准备与ConA磁珠固定CUTRUN的第一步是制备细胞悬液和细胞核。贴壁细胞用胰酶或者细胞刮刀收下来离心清洗悬浮细胞相对简单直接计数取用组织样本则需要先消化成单细胞悬液或者提取细胞核。关键要求是细胞状态好、碎片少、死细胞比例低——死细胞会把基因组DNA降解成大分子粘稠物严重影响后续的磁珠结合和切割。接下来是ConA刀豆蛋白A磁珠结合。ConA识别细胞核膜上的甘露糖残基把完整的细胞核粘到磁珠上。这里有个技术细节ConA结合需要钙离子和锰离子的辅助所以Binding Buffer里必须加1 mM的CaCl₂和1 mM的MnCl₂否则结合效率会大幅下降。结合时间一般10-15分钟室温即可。结合完成后的清洗操作要格外轻柔。细胞核靠着糖蛋白-凝集素的非共价相互作用挂在磁珠上剧烈吹打会把核弄破一旦核膜破裂染色质外泄整个实验就废了。我的习惯是磁分离后用枪头缓慢吸走上清不要吹悬磁珠而是加溶液后轻弹管底让沉淀松动。2.2 抗体孵育与pA-MNase结合决定成败的两步CUTRUN的抗体孵育和ChIP-seq有个重要差异它基本不做交联识别的是天然构象的表位。这意味着ChIP级别的抗体不一定适合CUTRUN而一些在免疫荧光、流式或native ChIP中表现好的抗体反而可能更合适。同一支抗体在这两种实验里表现可能差异巨大所以刚上手时务必做抗体预实验。抗体用量一般参考说明书常规体系5×10⁵细胞常用0.5-2 µg/反应。孵育条件有两种4°C过夜或室温2小时。室温孵育信号更高但非特异背景也会略上升4°C过夜则更温和适合对温度敏感的转录因子。我建议第一次跑可以从室温1-2小时起步用IgG对照同时做。孵育结束洗掉未结合抗体后加入pA-MNase1:100到1:200稀释具体浓度需按批次优化。这里特别提醒很多商业化的pA-MNase推荐在4°C结合1小时千万别在室温放太久不然酶还没被激活就可能降解核内的DNA。另外如果你的一抗是鼠源单抗你会发现protein A对小鼠IgG的亲和力不太够这时候要用pAG-MNaseprotein A/G融合它同时结合兔源和人源、鼠源的各类IgG亚型。注意pA-MNase的浓度不是越大越好。酶多了即使在0°C也会产生非特异切割酶少了目标区域切不充分。首次上手最好做一个1:50、1:100、1:200、1:400的梯度用qPCR检已知靶位点的富集倍数来定。2.3 钙离子激活切割与片段释放pA-MNase结合并清洗干净后加入CaCl₂到终浓度2 mM在冰上进行30-60分钟的切割。MNase是一个钙依赖核酸酶没有钙离子它完全不工作加了钙之后立刻激活。冰上处理的目的是让酶和染色质都处于低温低扩散状态最大程度限制切割发生的位置减少背景。切割结束后加入含EGTA的Stop Buffer终止反应。EGTA是钙离子的螯合剂能瞬间把游离钙离子络合掉让MNase失活。在Stop Buffer这一步骤要加入E. coli染色质或者E. coli细胞作为spike-in内参——这是CUTRUN体系里一个非常重要的实验设计后面数据分析部分我会专门展开。随后将反应管离心600 g5分钟或用磁力架分离磁珠上清就是释放出来的DNA片段。这里有个容易被忽略的点部分片段可能还停留在核内没完全释放可以选择加一次释放缓冲液再离心把产量提高20%左右。上清转移到新的低吸附离心管进入DNA纯化。2.4 从DNA到文库小体系操作的细节CUTRUN产物极其微量一个转录因子实验通常只能拿到几纳克到几十纳克的DNA。这个数量下任何管壁吸附都可能是致命的。我第一次做的时候用的普通EP管最后纯化洗脱下来跑Qubit显示只有0.3 ng/µL当时一度怀疑整个流程出了问题。后来换成DNA LoBind管同批次样品直接翻了三倍。DNA纯化可以用酚氯仿抽提加乙醇沉淀也可以用专门的微量DNA纯化柱。如果用柱子务必选用对短片段150 bp回收效率高的产品很多通用胶回收柱对小片段的结合效率很低。沉淀时加GlycoBlue糖原类载体能显著提高微量DNA的回收率。文库构建方面CUTRUN产物不需要机械打断直接做末端修复、加A尾、连接接头、PCR扩增。由于起始量极低PCR循环数通常在12-14个循环这里要平衡两个问题循环数太少文库量不够上机循环数太多PCR重复增多、复杂度下降。遇到产量特别低的样本宁可接受稍低文库产量也不要盲目加到16个循环以上否则后续去重后有效reads会少得可怜。片段大小筛选也要注意做组蛋白修饰时目标产物是单核小体约150 bp及少量二核小体做转录因子时目标产物主要是亚核小体片段30-120 bp可以优先回收小片段。切胶或双轮Ampure XP磁珠都可以实现但我更推荐后者操作稳定且损耗小。3. CUTRUN、CUTTag与ChIP-seq三种路线到底怎么选3.1 核心差异对照CUTRUN火了之后紧接着又出现了CUTTagCleavage Under Targets and Tagmentation。很多人把它们混为一谈其实分工并不一样。三个技术可以放在一张表里看维度ChIP-seqCUTRUNCUTTag切割工具超声物理打断pA-MNase定点酶切pA-Tn5转座酶起始细胞量10⁶-10⁷10⁴-10⁶可低至千级10²-10⁵可单细胞交联通常需要甲醛交联多数不需要通常不需要建库方式片段化后连接接头纯化DNA后连接接头转座酶直接插入接头背景信号较高低较低实验周期2-3天1-2天1-2天适用抗体已广泛验证需重新验证需重新验证成本构成测序深度高、耗材常规酶/试剂盒较贵、测序便宜试剂盒贵、测序便宜CUTTag的核心是把Tn5转座酶融合到protein A上抗体定位之后转座酶直接在目标区域把测序接头标进染色质片段化和建库一步完成。它的细胞量门槛比CUTRUN更低做到单细胞也是它在拓展而且不需要单独纯化DNA操作更简单。3.2 我实际使用中的选型建议先给结论如果你的细胞量在5万以上目标是转录因子或组蛋白修饰且希望拿到比较纯粹的蛋白-DNA互作信号CUTRUN是首选。它的切割机制决定了信号更接近真实的结合位点分辨率高峰型尖锐背景极低。如果细胞量在1万以下特别是几百到几千个细胞CUTTag更合适。Tn5转座酶的效率远高于MNase在极低起始量下建库成功率更高。如果你的目标最终是单细胞水平那基本只有CUTTag一个选项。如果样本是临床FFPE组织或者你的抗体库全部是已发表的ChIP-seq抗体且没有时间去重新验证那ChIP-seq仍然是相对稳妥的路径。CUTRUN虽然是升级版但升级意味着规则变了抗体、内参、生信都要重新适配这个成本要提前算进去。还有一个容易被忽视的维度如果你后续要做大量的组蛋白修饰全基因组图谱并且要和自己之前采集的ChIP-seq公共数据做大规模比较跨技术比较的可信度是个问题。CUTRUN的信号分布特征和ChIP-seq不一样差异分析时可能会引入批次效应这种情况除非所有新数据都用同一种技术补齐否则谨慎更换方法。4. 数据分析的关键环节比对、归一化、找峰4.1 比对管线的参数优化很多人觉得CUTRUN的生信分析可以照搬ChIP-seq流程实操后就会发现不是那么回事。CUTRUN的插入片段集中在30-150 bp双端测序时两条reads的比对特征和普通ChIP-seq很不一样。推荐管线是质量评估FastQC→ 去接头cutadapt或trim_galore→ 比对bowtie2→ 过滤低质量比对和重复 → 生成bigWig可视化 → 寻找peaks。bowtie2的参数建议这样设置bowtie2 -p 8 --end-to-end --very-sensitive --no-unal \ --no-mixed --no-discordant --phred33 -I 10 -X 700 \ -x /ref/bowtie2/hg38 -1 clean_R1.fq -2 clean_R2.fq -S sample.sam--no-mixed和--no-discordant是为了过滤掉非一致性比对因为CUTRUN片段短理论上双端reads应该都能唯一定位到同一区域。比对完成后用-q 10过滤掉低质量比对再去除PCR重复。加性别决定染色体区域去除可选但至少要把线粒体reads去掉它们经常占比非常高。4.2 Spike-in归一化跨样本比较的基准CUTRUN实验里不同样本之间的磁珠结合效率、清洗损失、酶切效率都会带来系统差异。如果不做归一化两个样本的绝对信号高度无法直接比较。E. coli spike-in就是用来解决这个问题的。操作层面你在Stop Buffer里加入恒定数量的E. coli细胞或E. coli染色质它们会经历和你的样本相同的后续所有步骤。E. coli基因组上没有你目标蛋白的结合位点理论上不会产生特异性信号但它们的DNA也会被环境中残余的MNase非特异切割、释放、纯化、建库。所以样本间spike-in reads的比例就反映了整个实验流程的效率差异。分析时先把reads分别比对到目标基因组如hg38和E. coli基因组如K12 MG1655统计两个比对文件里的有效reads数。归一化系数可以用# 计算样本spike-in比例以最小spike-in reads样本为基准 scale_factor min_ecoli_reads / sample_ecoli_reads然后用这个系数对bam文件缩放生成可视化trackbamCoverage -b sample.sorted.bam -o sample.scaled.bw \ --scaleFactor 1.32 --binSize 10 --normalizeUsing nonespike-in reads一般占总reads的1%-5%比较理想。如果某个样本的E. coli比例特别高说明酶切过度或者清洗不充分如果特别低说明整个流程效率差这个样本的信号可信度就要打问号。4.3 SEACR与MACS2的取舍找峰的软件选择是CUTRUN数据分析里争论最多的点。Henikoff实验室专门开发了SEACRSparse Enrichment Analysis for CUTRUN核心思路是用全基因组信号分布来替代背景模型不依赖IgG对照。MACS2虽然是ChIP-seq时代的标准工具但在CUTRUN上直接使用会产生很多假峰因为两者的负二项分布假设并不匹配。我的经验是转录因子窄峰建议优先试SEACR的stringent模式组蛋白修饰宽峰可以试relaxed模式。SEACR需要输入bedgraph文件命令如下# 生成bedgraph bedtools genomecov -bg -scale 0.25 -ibam sample.sorted.bam sample.bedgraph # SEACR找峰stringent模式 SEACR_1.3.sh sample.bedgraph IgG.bedgraph stringent sample_stringent # SEACR找峰relaxed模式适合宽峰 SEACR_1.3.sh sample.bedgraph IgG.bedgraph relaxed sample_relaxed如果要用MACS2需要调整参数适配CUTRUN的短片段特征macs2 callpeak -t sample.bam -c IgG.bam -f BAMPE -g hs \ --keep-dup all -q 0.05 --nomodel --shift -100 --extsize 200 \ -n sample实际使用中我会把SEACR和MACS2都跑一遍然后在IGV里目视检查一组已知位点。两个工具结果重合的部分基本可信只有其中一个检出的峰需要重点怀疑。质量评估还要看reads in peaks比例FRiPCUTRUN做得好时这个比例能到20%-60%远高于ChIP-seq的5%-10%。如果FRiP低于5%先别急着找峰回头排查抗体和酶切环节。5. 我踩过的坑以及一份排查清单5.1 高背景问题往往出在酶和温度实验室接CUTRUN的第一个月某个组蛋白修饰样本跑出来的reads在基因组上均匀分布完全看不到峰。当时第一反应是抗体不行结果排查一圈发现是pA-MNase浓度过高加室温孵育时间太长导致非特异切割全面开花。这件事给我的教训是CUTRUN的低背景不是自动获得的它的先决条件是酶被严格控制在抗体标记的局部。压背景有四个关键旋钮pA-MNase浓度宁低勿高、CaCl₂激活温度必须冰上、激活时间不超过60分钟、清洗次数抗体和酶结合后的清洗不能省。出现高背景时先检查这几个因素甚至可以把清洗次数从两次加到三次背景会显著下降。5.2 文库产量极低可能死于管壁和洗脱有一次给一个只有1万个细胞的稀有样本做CUTRUN最后Qubit检不出DNA浓度。当时浪费了一整个辛苦分选的样本。后来复盘发现两个致命问题一是用了普通离心管微量DNA被管壁吸附二是DNA纯化时用乙醇沉淀但温度没达到-20°C沉淀时间也不够。从那以后我立了几条规矩全程用LoBind管乙醇沉淀至少-20°C过夜并加GlycoBlue纯化完上机前不反复冻融文库PCR循环数严格控制在14以内如果产量不够宁可选更灵敏的建库试剂盒比如NEBNext Ultra II的低起始量方案也不盲目增加扩增循环。5.3 抗体不兼容CUTRUN有自己的抗体脾气CUTRUN对抗体的要求是天然构象下能识别目标蛋白。一支说明书标注ChIP Grade的抗体可能在甲醛交联后暴露的表位结合良好但在未交联的天然染色质上完全失效反过来有些IF级别抗体反而在CUTRUN里出奇地好。我建议正式实验前做一次抗体小筛用你的目标细胞和一支已知阳性抗体比如H3K4me3或H3K27ac这种在CUTRUN中很稳定的抗体做平行测试同时用IgG做对照跑一次低深度测序或至少qPCR验证靶位点富集。如果阳性抗体的FRiP明显高、IgG背景低说明整个体系是通的剩下的就是换目标抗体了。这一步看起来很费时间实际上能帮你少走两个月的弯路。5.4 质控指标速查表现象可能原因处理建议FRiP 5%抗体无效/浓度低、酶量不足、细胞状态差换抗体、做抗体梯度、检查细胞活力Spike-in比例过高15%清洗不充分、酶切过度增加清洗次数、降低pA-MNase浓度Spike-in比例过低0.5%流程损失大、细胞量太少用LoBind管、延长沉淀时间文库无PCR产物微量DNA丢失、接头连接失败加GlycoBlue、检查接头连接体系Peak集中在少数位点PCR重复过多、测序深度不足增加起始量或循环数提高reads数组蛋白修饰信号杂乱酶切时间过长、钙离子过量缩短激活时间、检查Stop Buffer浓度6. 上手CUTRUN之前值得想清楚的几件事6.1 你的样本类型是否适配不是所有样本都适合直接套用标准CUTRUN流程。植物细胞有细胞壁需要先制备原生质体或纯化细胞核临床组织样本需要先做单细胞化或核提取组织消化过程中的细胞死亡会直接影响ConA结合和后续切割冻存细胞原则上能做但冻存后细胞膜通透性和核膜状态都会变化建议先用活细胞做一版完整流程熟悉体系再拿冻存样本优化。另外CUTRUN的清洗步骤依赖磁珠高脂含量或黑色素含量高的样本会干扰磁分离和后续酶活这类样本光做细胞悬液准备可能就要额外花一周调试要做好心理准备。6.2 抗体验证的优先级如果让我用一句话总结CUTRUN最值得投入的地方我会说把70%的优化精力放在抗体上。实验流程是标准化的酶是商业化的数据分析管线也有成熟的工具唯独抗体是个变量。一支好抗体的CUTRUN数据和一支差抗体的数据差距比任何参数调整带来的差异都大。我验证抗体的标准流程是先用Western blot确认特异性再看是否有已发表的CUTRUN/CUTTag数据使用记录然后做一个起点为2 µg/反应、逐步降低浓度的测试最后在IGV里目视检查已知阳性位点和阴性区域。三个环节都过了这支抗体才真正入库。6.3 预算与外包的选择CUTRUN的试剂盒比如CST的CUTRUN Kit或EpiCypher的CUTANA系列价格不低单样本成本通常几百到上千元但它省掉的测序费是惊人的——CUTRUN单样本一般测10-20M reads就够了只有ChIP-seq的一半甚至三分之一。如果你的实验室没有现成的实验体系又急着出结果外包给有成熟CUTRUN平台的测序服务商也是一个务实选择。前提是确认服务商会提供spike-in归一化、SEACR分析这类CUTRUN专属的生信流程而不是拿ChIP-seq的流程糊弄你。最后再分享一个小技巧CUTRUN上机前可以用qPCR快速预检几个已知阳性位点和阴性位点的富集倍数。阳性位点比阴性位点高几十倍以上再送测序基本稳了。这个小步骤只花两个小时却能避免烧掉一次几千块的测序试错成本。实验做得越多越明白表观遗传学领域真正稀缺的不是花哨的技术而是能稳定复现的结果。
 的 API 用法、生命周期回调与双通道上传原理)
