《Nature Communications》：mist: a hierarchical Bayesian framework for detecting differential DNA methylation dynamics in single-cell data
编辑推荐：
单细胞DNA甲基化（scDNAm）测序技术揭示了细胞异质性，但缺乏沿假时间轴分析甲基化动态变化的方法。本研究开发了层级贝叶斯框架MIST，用于在scDNAm数据中识别沿伪时间轨迹的差异甲基化特征。应用于小鼠胚胎发育和人类大脑发育多组学数据集，MIST成功鉴定了与谱系转换一致的关键发育调控因子。该工具填补了方法学空白，为理解发育和疾病中的表观遗传动力学提供了新途径。
在生命科学的微观世界里，每个细胞都携带着独特的身份密码，这不仅由基因序列决定，更受到一层精密的“化学修饰”调控，这就是表观遗传学。DNA甲基化（DNA methylation）作为其中一种关键修饰，如同基因上的“开关标签”，深刻影响着基因的活跃与否，从而主导着细胞的命运，如分化、发育乃至癌变。近年来，单细胞DNA甲基化（scDNAm）测序技术的突破，让我们得以以前所未有的分辨率窥见细胞群体中复杂的表观遗传“风景画”，理解细胞间的异质性。科学家们还发展出“伪时间轨迹分析”技术，能够像给细胞拍一部动态电影一样，将看似静态的细胞群体按照其发育或分化过程排序，从而追踪基因组变化的连续过程。
然而，在这幅日益清晰的画卷中，仍有一块关键的拼图缺失：尽管我们能将细胞沿着伪时间轴排列，但目前却没有任何专门的方法，能够系统地从scDNAm数据中，分析和建模DNA甲基化沿这条假想时间轴的动态变化规律。具体来说，我们缺乏有效的工具来精准识别在连续生物过程（如胚胎发育）中，哪些基因组位点或区域发生了显著且具有生物学意义的甲基化变化。这一方法学的空白，限制了我们深入挖掘表观遗传动态在细胞状态转变中的驱动机制。为了解决这一紧迫问题，一个名为“mist”的新方法应运而生，相关研究发表在《Nature Communications》上。
为开展研究，研究人员主要运用了单细胞DNA甲基化测序技术来获取基础数据。研究核心是提出了一个创新的层级贝叶斯统计框架。该研究使用了模拟数据验证方法的准确性，并应用于两个真实世界的多组学数据集：一是小鼠胚胎发育数据集，二是发育中的人类大脑数据集，以验证其发现生物学意义的能力。
研究结果
1. MIST：一个用于单细胞DNA甲基化轨迹分析的层级贝叶斯框架
研究人员开发了名为“mist”的计算方法。其核心创新在于建立了一个概率统计模型，该模型能够同时考虑单细胞甲基化数据的高噪音特性、细胞间的异质性，并建模甲基化水平沿用户定义的伪时间轨迹的连续变化。通过引入层级先验分布，MIST能够更稳健地估计参数，并识别出在轨迹上表现出显著动态甲基化模式的基因组区域。
2. 在模拟数据中验证性能
通过构建复杂的模拟数据集，研究人员将MIST与现有的、非专门为轨迹分析设计的差异甲基化分析方法进行比较。结果显示，MIST在识别沿伪时间变化的差异甲基化基因方面，表现出更高的准确性和统计效力，证实了其模型在处理此类问题上的优越性。
3. 揭示小鼠胚胎发育中的表观遗传动力学
将MIST应用于小鼠早期胚胎发育的scDNAm数据后，研究成功绘制了DNA甲基化重编程过程沿发育轨迹的动态图谱。分析鉴定出了一系列在特定发育阶段（例如，从受精卵到囊胚的谱系转换过程中）发生剧烈甲基化变化的关键基因组位点。这些位点富集于发育调控基因和印记控制区域，表明MIST能够捕捉到与细胞命运决定密切相关的表观遗传事件。
4. 解析人脑发育中的甲基化变化轨迹
在发育中的人类大脑前额叶皮层scDNAm数据中，MIST被用于分析神经发生和皮质层形成的伪时间过程。研究发现，许多在轨迹上显示差异甲基化的基因，是已知的神经发育关键调控因子。这些基因的甲基化动态变化模式与神经元的分化成熟阶段高度同步，为理解人类大脑皮层构建过程中精确的表观遗传时序调控提供了新线索。
结论与讨论
本研究成功开发并验证了首个专门为分析单细胞DNA甲基化数据沿伪时间轨迹变化而设计的计算框架——MIST。该框架通过层级贝叶斯模型，有效解决了从高噪声、稀疏的单细胞数据中推断连续甲基化动态的统计学挑战。研究不仅证明了MIST在模拟环境中优于现有方法，更重要的是，通过将其应用于小鼠胚胎和人类大脑发育这两个关键生物学体系，揭示了一系列与细胞谱系转换和发育阶段紧密相关的差异甲基化基因与调控元件。这些发现证实了MIST在发现具有生物学功能的表观遗传动态事件方面的强大能力。
该研究的重要意义在于填补了当前单细胞表观基因组学分析方法中的一个重要缺口。MIST为研究人员提供了一个强大的工具，能够从动态的、非线性的生物学过程视角，系统探索DNA甲基化在发育、分化、疾病演进（如癌症进展）等连续状态转变中的作用。它推动单细胞表观遗传学研究从静态的“细胞图谱”绘制，迈向动态的“调控过程”解析。作为开源R/Bioconductor软件包，MIST的发布将极大促进该领域的研究，有望在未来帮助科学家们更深入地解码复杂生物学过程中表观遗传调控的精确时空调控密码。
