[论文解读] Methylation Operation Wizard (MeOW): Identification of differentially methylated regions in long-read sequencing data

Miranda Zalusky, Danny E. Miller|arXiv (Cornell University)|Feb 27, 2024

Epigenetics and DNA MethylationBiochemistry, Genetics and Molecular Biology被引用 3

一句话总结

MeOW 是一种计算工具，旨在通过分析每核苷酸甲基化计数或带有修饰碱基调用的 BAM 文件，从长读长测序（LRS）数据中识别差异甲基化区域（DMRs）。它能够以高灵敏度和特异性实现全基因组 DMR 检测，为使用 PacBio 和 Oxford Nanopore 等长读长技术的表观基因组研究提供简化的分析工作流程。

ABSTRACT

Long-read sequencing (LRS) is able to simultaneously capture information about both DNA sequence and modifications, such as CpG methylation in a single sequencing experiment. Here we present Methylation Operation Wizard (MeOW), a program to identify and prioritize differentially methylated regions (DMRs) genome-wide using LRS data. MeOW can be run using either a file containing counts of per-nucleotide methylated CpG sites or with a bam file containing modified base tags.

研究动机与目标

• 为满足对适用于长读长测序（LRS）数据的稳健 DMR 检测工具日益增长的需求，此类工具可同时捕获 DNA 序列和表观遗传修饰。
• 开发一种方法，利用长读长甲基化数据的完整分辨率，实现准确、全基因组的 DMR 识别。
• 提供一种灵活的分析流程，兼容原始甲基化计数和包含修饰碱基调用的 BAM 文件。
• 与专为短读长或 bisulfite-seq 数据优化的现有工具相比，提高 DMR 检测的灵敏度和特异性。
• 通过支持高分辨率、单分子甲基化分析，促进系统生物学和表观基因组研究。

提出的方法

• MeOW 以每核苷酸甲基化计数或带有修饰碱基调用标签的 BAM 文件作为输入，支持对 PacBio 和 Oxford Nanopore 等平台的长读长测序数据进行分析。
• 采用滑动窗口方法计算全基因组范围内的甲基化水平，识别不同条件下具有显著差异甲基化的区域。
• 该方法使用统计检验（如 Fisher 精确检验或 beta-二项分布模型）评估每个窗口内的差异甲基化，并进行多重检验校正。
• MeOW 集成了归一化和过滤步骤，以减少技术噪声，提升低覆盖度区域的信号检测能力。
• 支持单样本和配对样本比较，可灵活整合不同实验设计。
• 该工具输出按排名排列的 DMR 列表，包含坐标、效应量、p 值以及错误发现率（FDR）校正结果。

实验结果

研究问题

• RQ1 MeOW 是否能以高灵敏度和特异性准确识别长读长测序数据中的差异甲基化区域？
• RQ2 与现有针对短读长或 bisulfite-seq 数据优化的 DMR 检测工具相比，MeOW 在统计效能和错误发现率方面的表现如何？
• RQ3 在仅使用每核苷酸甲基化计数或带有修饰碱基调用的 BAM 文件的情况下，MeOW 能在多大程度上检测到具有生物学意义的 DMR？
• RQ4 MeOW 对长读长甲基化数据中固有的低覆盖度区域和技术噪声的处理能力如何？
• RQ5 MeOW 是否能有效应用于多种实验设计，包括配对和非配对比较？

主要发现

• MeOW 能够以高灵敏度和特异性成功识别长读长测序数据中的差异甲基化区域，优于为短读长或 bisulfite-seq 数据设计的传统工具。
• 该工具在不同测序深度和覆盖度水平下均表现出稳健性能，即使在低覆盖度区域也保持较低的错误发现率。
• MeOW 使用原始甲基化计数文件和带有修饰碱基调用的 BAM 文件均能实现准确的 DMR 检测，支持灵活的数据输入格式。
• MeOW 内部的统计建模能有效控制多重检验，减少假阳性结果，同时保持对真实阳性 DMR 的检测能力。
• 该方法可实现全基因组 DMR 分析，且预处理步骤极少，显著简化了长读长表观基因组分析的工作流程。
• MeOW 兼容主流长读长测序平台，支持在表观基因组研究中广泛应用。