MicroMap：人类微生物组代谢网络可视化资源的构建与应用

《npj Biofilms and Microbiomes》：MicroMap: a network visualisation resource for human microbiome metabolism

本刊推荐：为填补大规模微生物组代谢网络可视化空白，研究团队开发了MicroMap——首个手工绘制的人类微生物组代谢网络图谱。该资源整合AGORA2和APOLLO数据库中257,429个微生物基因组尺度代谢模型（GEM），涵盖5,064个独特反应和3,499个代谢物（含98种药物代谢）。通过约束性重构与分析（COBRA）框架，支持代谢通量可视化、菌群代谢能力比较及动态动画展示，为营养-宿主-微生物组-疾病轴研究提供机制性假说生成平台。

　　

人类微生物组如同一个精密的代谢工厂，其化学反应网络深刻影响着宿主的健康与疾病。尽管基因组尺度代谢模型（Genome-scale metabolic models, GEMs）已成为研究微生物组功能的重要数学工具，但科学界长期缺乏直观展示微生物代谢网络的可视化平台。现有数据库如KEGG或BRENDA的通路图覆盖范围过广，缺乏对人类微生物组的特异性聚焦，且无法与主流的约束性重构与分析（Constraint-Based Reconstruction and Analysis, COBRA）建模框架直接整合。这种可视化工具的缺失，使得研究人员难以直观探索不同微生物的代谢差异，也无法有效展示计算模型得出的动态通量变化。

为解决这一瓶颈，来自爱尔兰高威大学等机构的研究团队在《npj Biofilms and Microbiomes》发表了名为“MicroMap”的人类微生物组代谢网络可视化资源。该研究通过手工绘制方式，构建了覆盖25万多个微生物GEM的代谢图谱，包含5,064个独特反应和3,499个代谢物，其中特别整合了98种常用药物的代谢途径。MicroMap采用城市地图风格的设计理念，将相关反应聚类展示，并设置337个区位标识指引生物化学子系统，使研究人员能够像查阅城市地图般直观导航复杂的微生物代谢网络。

研究团队开发的关键技术方法主要包括：基于CellDesigner软件的手工代谢网络绘制、利用COBRA工具箱的代谢地图学功能进行质量校验、通过高性能服务器并行生成257,429个微生物重构可视化图谱，以及开发新型MATLAB函数实现反应存在热图和通量动画可视化。研究使用的微生物基因组数据来源于AGORA2数据库的7,302个菌株水平GEM和APOLLO数据库的247,092个宏基因组组装基因组（MAG）衍生重构。

微生物代谢能力的可视化比较

通过将不同微生物的代谢重构可视化在MicroMap上，研究人员能够直观比较各类微生物的代谢特征。例如，对比芽孢杆菌纲（Bacilli）和疣微菌门（Verrucomicrobia）的代谢图谱时，可清晰发现前者具有更广泛的药物代谢能力。这种可视化方法使微生物代谢知识对非生物信息学背景的用户也变得易于理解，特别是在教育场景中具有重要价值。

热图揭示微生物间相对反应存在差异

研究团队进一步开发了反应存在热图功能，能够可视化比较不同微生物重构间的代谢能力差异。以AGORA2中包含的14个假单胞菌物种为例，热图分析显示这些物种在氟尿嘧啶代谢能力上存在明显差异，有些物种能够介导已知的微生物组相关氟尿嘧啶代谢反应，而其他物种则缺乏这种能力。热图通过线条粗细和颜色深浅（从细线海军蓝到粗线棉花糖色）直观展示反应的相对存在程度。

代谢建模结果的可视化揭示动态通量变化

MicroMap支持COBRA建模结果的通量可视化，能够展示代谢物在网络中的流动情况。研究人员通过纵向时间序列分析，生成了一系列通量矢量可视化图谱，并进一步创建帧-by-帧动画以突显随时间变化的通量动态。这种动画能够清晰显示通量方向和大小的变化，例如在谷氨酸代谢反应P5CD中观察到的动态波动，为识别潜在候选通路提供了直观手段。

机制性生化见解的生成

通过整合上述可视化方法，MicroMap能够帮助生成微生物组代谢的新机制性见解。以硫化氢（H₂S）代谢为例，案例研究显示MicroMap能够将特定代谢物或反应置于更广泛的网络环境中进行背景分析，包括上游前体、下游产物及相关生化子系统。研究发现H₂S代谢主要涉及硫、蛋氨酸/半胱氨酸和丝氨酸代谢子系统，存在多种直接和间接的H₂S形成途径。特别观察到反复出现的H₂S与亚硫酸根（SO₃^2-）相互转化模式，这种模式将有机和无机硫代谢途径相互连接。

研究还通过可视化发现，不同微生物在特定代谢过程中的参与能力存在显著差异。例如，脱硫弧菌（Desulfovibrio desulfuricans）表现出广泛的H₂S代谢能力，而巴氏乳杆菌（Lactobacillus pasteurii）则完全没有。通量平衡分析结果的可视化进一步揭示了不同条件下通路使用的差异，当没有直接半胱氨酸来源时，通量会通过合成半胱氨酸的反应进行路由，这表明脱硫弧菌在给定条件下内源性产生半胱氨酸。

该研究的创新性在于首次提供了大规模微生物组代谢网络的可视化平台，填补了该领域长期存在的技术空白。MicroMap与虚拟代谢人类数据库（Virtual Metabolic Human, VMH）完全集成，可与人类代谢重构图（ReconMap）形成互补，为研究营养-肠-脑-健康轴等人类-微生物组相互作用提供了重要工具。除了研究应用外，MicroMap还具有教育价值，已在Virtuome数字健康暑期学校等教学场景中成功应用，帮助非专业受众理解微生物组代谢的系统性特征。

尽管现有自动绘图工具如Fluxer能够快速生成GEM衍生网络，但其自动布局在生化背景下难以解释，且对自定义COBRA建模适应的灵活性有限。研究团队认为，手工创建MicroMap是实现COBRA背景下大规模微生物组代谢网络可视化的必要第一步，这将催化与现有平台的进一步整合和完善。随着微生物组研究的深入，MicroMap有望成为探索微生物组与人类健康复杂关系的重要桥梁，推动个性化医疗和营养干预策略的发展。