《Cell Reports Methods》：mChIP-seq for multiplex epigenomic profiling uncovers cancer-specific histone features in cellular and circulating nucleosomes
编辑推荐：
表观基因组分析阐明了组蛋白修饰的调控作用，并为发现非侵入性无细胞核小体表观遗传模型提供了关键途径。研究人员提出多重染色质免疫沉淀测序（mChIP-seq），该技术兼容细胞样本与无细胞样本，可同时解析多个样本的多因素表观遗传景观。通过将一步法样本索引与混合-拆分
表观基因组分析阐明了组蛋白修饰的调控作用，并为发现非侵入性无细胞核小体表观遗传模型提供了关键途径。研究人员提出多重染色质免疫沉淀测序（mChIP-seq），该技术兼容细胞样本与无细胞样本，可同时解析多个样本的多因素表观遗传景观。通过将一步法样本索引与混合-拆分策略结合，mChIP-seq提高了效率并降低了成本。利用mChIP-seq，研究人员对代表九种癌症类型的细胞系中H2A.Z及10种组蛋白修饰进行了分析，发现与非癌对照相比，癌细胞启动子区H2A.Z与H3K4me3的相关性降低。研究人员进一步扩展mChIP-seq，建立无细胞mChIP-seq（cf-mChIP-seq），用于绘制循环无细胞核小体。对血浆样本中H2A.Z、H3K4me3、H3K27ac和H3K9me3进行分析，发现乳腺癌与健康对照之间存在与组蛋白修饰相关的独特cfDNA片段模式，突显cf-mChIP-seq扩展液体活检方法学的潜力。总体而言，mChIP-seq为细胞与无细胞情境下的表观基因组分析提供了稳健、可扩展的平台。
研究背景与问题：核小体（nucleosome）是染色质基本重复单位，由DNA缠绕组蛋白八聚体构成；组蛋白翻译后修饰（PTMs）和组蛋白变体构成动态表观遗传景观，调控基因表达。细胞死亡后染色质释放至血液等细胞外液，保留一部分无细胞核小体（cf-nucleosomes），其组蛋白标记可编码细胞类型相关特征，为液体活检（liquid biopsy）提供非侵入性表观遗传标志物。但标准ChIP-seq（染色质免疫沉淀测序）起始材料需求高，现有样本多重化ChIP-seq步骤繁琐，且对无细胞样本适应性不足；cfChIP-seq仍受血浆起始量大和逐个免疫沉淀通量低限制。H2A.Z（组蛋白H2A变体）在多种癌症中过表达并与预后不良相关，但缺乏全面表观遗传景观数据。为此，研究人员建立mChIP-seq（多重染色质免疫沉淀测序）并扩展cf-mChIP-seq，用于细胞和无细胞样本的高通量、低输入、低成本表观基因组分析，并在泛癌及乳腺癌血浆中研究H2A.Z及相关组蛋白修饰。
关键技术方法：研究人员利用末端脱氧核苷酸转移酶（TdT）加尾和单链DNA（ssDNA）连接实现一步样本索引，将索引样本混合后拆分进行免疫沉淀和测序，保留常规ChIP-seq主要步骤；验证用人HEK293T、HeLa、小鼠3T3、酵母BY4741及猪肝胰组织。细胞分析使用24株覆盖9种癌症类型的NCI-60细胞系，包括ENCODE Tier 1 K562；血浆队列38例女性（21例乳腺癌患者、17例健康对照；广州医科大学妇女儿童医学中心伦理批准），cf-mChIP-seq每谱用25 μL血浆；并整合公开ENCODE、GEO和RNA-seq数据。
研究结果：
（1）mChIP-seq设计与验证：通过一步索引和混合-拆分，mChIP-seq将样本条形码步骤减少为一步，低至约3,000细胞/靶标，24样本1–12靶标时每文库成本较常规ChIP-seq降低3–16倍；多物种多修饰验证中超过95% reads分配至样本，交叉污染低，与常规ChIP-seq信号高度相关，并兼容组织Pol II（RNA聚合酶II）分析。
（2）泛癌组蛋白标记图谱：用mChIP-seq在24个癌细胞系中分析H2A.Z和10种H3修饰，获得高分配率、唯一比对率数据；H2A.Z、H3K4me1/2/3、H3K27ac、H3K18ac在TSS周围呈双峰，H3K36me3在基因体富集；K562数据与ENCODE重叠和相关较高，FRiP信噪比相当。
（3）癌细胞H2A.Z表观特征：鉴定188,464个H2A.Z富集位点，平均19,655/样本；启动子51.73%、内含子19.72%、远端基因间24.31%；H2A.Z与启动子标记、活性标记及ATAC-seq可及性正相关；启动子和远端H2A.Z与表达正相关，内含子H2A.Z与宿主基因表达降低相关，伴活性标记低、H3K27me3高、H3K36me3低，公共数据重现。
（4）癌细胞H2A.Z表观失调：癌细胞H2AZ1/H2AZ2表达高于非癌；启动子H2A.Z与H3K4me3相关性在癌细胞中变异更大且显著低于非癌，而H3K4me3与H3K27ac或H3K27me3关联无显著差异；HCT116中不同H3K4me3强度组H2A.Z信号相近，IMR-90中则相应变化，其他癌细胞可见类似模式。
（5）启动子H2A.Z积累对表达影响：按H3K4me3强度分组，最高H3K4me3对应最高表达；启动子有H2A.Z但无H3K4me3的基因在癌细胞表达更高，并富集于代谢、生物合成和细胞周期调控；SRCAP、EP400、ANP32E在癌细胞表达较高，其与H2AZ1/H2AZ2的相关模式提示反馈调控可能，但因果关系原文称需验证。
（6）cf-mChIP-seq：将流程适配血浆，25 μL血浆/谱即可同时分析多个修饰；H3K4me3基准中，25 μL cf-mChIP-seq回收cfChIP-seq约一半peaks，未回收者多为低富集区，片段大小呈单/二核小体模式，信号高度相关。
（7）乳腺癌血浆循环核小体分析：38例血浆、4种标记获得156个图谱；乳腺癌患者启动子H2A.Z与H3K4me3相关性低于健康对照；H2A.Z相关cfDNA在患者富集启动子更短；不同组蛋白标记相关cfDNA中位长度不同（H3K4me3最长，其次H3K27ac、H2A.Z、H3K9me3，均长于Input）；癌症样本H2A.Z、H3K27ac、H3K9me3的二核小体比例显著降低；基于全基因组片段模式的弹性网分类器对H2A.Z和H3K9me3的AUC分别达0.988和0.963。
讨论与结论：讨论指出mChIP-seq相较现有多重方法有ssDNA连接效率高、一步索引减少步骤和样本损失、兼容cfDNA单链/缺口/超短片段等优势；H2A.Z调控作用依赖基因组位置，癌症中H2A.Z过表达可能驱动启动子异常沉积并与H3K4me3相关性降低，低H3K4me3基因表达升高涉及代谢、生物合成和细胞周期过程，但需过表达/敲低实验验证。cf-mChIP-seq将高通量表观分析扩展至循环核小体；不同组蛋白标记相关cfDNA片段模式可区分乳腺癌与健康对照，H3K9me3回收受染色质间距和Tn5空间位阻影响；癌症相关片段更短，需更大队列和多癌种验证。局限性包括抗体特异性依赖、低输入导致较高PCR重复、细胞系对交联/超声敏感性差异，以及仅乳腺癌小样本限制机器学习模型性能。研究结论翻译：总之，研究人员建立了mChIP-seq及其衍生cf-mChIP-seq，作为高效、低成本、低输入需求的组蛋白修饰表观基因组分析工具；研究工作揭示癌细胞中H2A.Z调控景观改变，表现为启动子区与H3K4me3相关性降低，以及乳腺癌患者血浆中H2A.Z富集cfDNA的独特片段模式。研究人员预期mChIP-seq有潜力支持ENCODE等大规模项目，cf-mChIP-seq可扩展液体活检方法学并用于多种临床场景。
