Genome Res：用于多样本ChIP-seq数据分组定量比较的MAnorm2计算模型

　　

染色质免疫共沉淀测序（ChIP-seq）实验被广泛用于刻画转录因子结合和组蛋白修饰的全基因组分布。比较来自不同细胞类型的ChIP-seq样本是刻画细胞分化和病变过程中动态转录和表观调控的关键基础。2012年，邵振与张一婧等合作在Genome Biology发表了用于两个ChIP-seq样本之间进行一对一定量比较的MAnorm模型。近年来，随着实验技术的发展和测序成本的不断降低，在ChIP-seq样本组（而非单个样本）之间进行比较分析已成为越来越常见的研究需求。一方面，研究者会产生同一实验的多个生物学重复来提高实验结果的可信度。另一方面，通过将来自不同个体的样本根据特定标签（如年龄、性别、患病与否、疾病亚型等）分组进行比较，研究者能够控制个体差异造成的影响，更可靠地识别与该标签关联的差异结合位点。然而，由于ChIP-seq实验固有的高复杂度和高噪声水平，以及不同比较场景所特有的技术困难，现阶段对多样本ChIP-seq数据进行分组定量比较仍然是一个巨大的计算方法学挑战。

11月19日，Genome Research在线发表了中国科学院上海营养与健康研究所中科院计算生物学重点实验室（马普伙伴计算生物学研究所）邵振课题组的方法学论文“MAnorm2 for quantitatively comparing groups of ChIP-seq samples”，报道了其开发的新一代MAnorm2计算模型。该模型能够对多样本ChIP/ATAC-seq数据按照特定标签分组进行统计建模和组间定量比较，可靠地在样本组层面鉴定组间显著差异的ChIP/ATAC-seq信号。

在ChIP-seq数据标准化这一步，MAnorm2沿用了MAnorm的核心假设，通过重构其信号强度变换体系，新发展了以参照样本为基准的多样本并行ChIP-seq信号标准化流程。进一步，针对多样本分组比较的需求，MAnorm2搭建了一个理论上适应任意树状分组结构的层级化多样本标准化策略。在完成标准化后，MAnorm2接下来针对每个基因组区域上观察到的ChIP-seq信号组间差异进行统计检验。在通常组内样本数较少的局限下（2-3个重复本），为了能更准确地衡量每一个基因组区域上的组内样本间ChIP-seq信号变化水平（within-group variability），MAnorm2设计了一个经验贝叶斯框架，利用拟合均值-方差曲线来给单个区域的组内变化水平赋予一个先验分布，并进一步通过平衡先验和后验观测来更准确地估计ChIP-seq信号的组内变化水平，从而提高对组间差异ChIP-seq信号的灵敏度。

与已有的其他经验贝叶斯方法相比，MAnorm2最大的优势在于考虑了不同样本组的组内ChIP-seq信号变化水平可能存在系统性差别。这一情形在正常人和癌症患者之间的比较中经常出现：由于肿瘤组织或血液样本本身的异质性以及癌症亚类型和不同患病阶段的多样性，癌症样本组的组内信号变化水平往往显著高于正常样本组。为了解决这个问题，MAnorm2通过在建模过程中引入一个方差比率因子，首先把不同样本组的全局组内信号变化水平修正至一致，然后使用修正后的方差进行均值-方差曲线的拟合和参数估计。研究人员将MAnorm2与现有的其他ChIP-seq差异分析工具进行了系统地比较，发现MAnorm2展现了明显更优越的使用性能，尤其是当进行比较的样本组拥有明显不同的组内变化水平时，例如癌症和正常样本相比较。

此外，该模型的应用场景和统计模型具有良好的可扩展性。在正文中，研究人员不仅展示了MAnorm2在ATAC-seq数据差异分析上同样适用，还将其统计模型扩展到可以同时比较任意多个样本组，并发现其使用效果优于传统的ANOVA方法。

该研究由中国科学院上海营养与健康研究所等多家机构合作完成。中科院营养健康所博士后涂世奇为该论文第一作者，邵振研究员为通讯作者。中科院植物生理生态研究所张一婧研究员、美国西南医学中心徐剑教授和波士顿大学的David J. Waxman教授对该研究工作的提出和完善也做出了重要贡献。该研究获得了国家自然科学基金委（31871280 和31701140）、科技部、中科院等机构的资助。

原文标题：

https://genome.cshlp.org/content/early/2020/11/18/gr.262675.120