《Microbiology Spectrum》：SparCC co-occurrence networking reveals intracommunity dynamics of the microbiome following colorectal surgery
编辑推荐：
本研究运用新兴的SparCC共现网络分析方法，深入探讨了结肠镜检查与结直肠手术对患者肠道菌群的生态扰动与恢复过程。通过比较三种不同干预（单独MBP、SBP非切除术、SBP切除术）下患者术后长达6个月的菌群网络拓扑结构变化，文章系统评估了网络拓扑参数（如模块性、边密度）的临床解释价值，并识别了在手术恢复关键时间点起枢纽作用的物种（如Enterocloster、Ruthenibacter），为理解肠道菌群在术后伤口愈合中的作用提供了新的生态学视角。结果表明，网络分析可作为传统多样性分析（如β多样性）的补充工具，为优化围手术期管理、预防手术部位感染（SSI）提供了新思路。
肠道菌群在术后伤口愈合中扮演着关键角色，但围手术期预防性抗生素的使用可能对其产生不利影响。明尼苏达结直肠手术观察性研究（MN MiCROS）旨在通过16S rRNA扩增子测序，评估不同扰动程度下患者肠道菌群的围手术期及长期变化。本研究运用SparCC共现网络分析，深入探究了微生物群落内部的相互作用动态，以揭示更具生物学相关性的模式，并评估该方法的转化应用价值。
研究方法与设计
该研究纳入了81名成年患者，分为三组：仅接受机械性肠道准备（MBP）的结肠镜检查组、接受手术肠道准备（SBP）的非结直肠切除术组，以及接受SBP的结直肠切除术组。研究在基线、手术/操作当天（POD0）、术后10天（POD10）、术后3-6周（POD30）和术后6个月（POD180）收集粪便样本。利用16S rRNA测序数据，在属水平上构建了SparCC共现网络，并使用NetCoMi软件包计算网络拓扑参数和枢纽物种，通过置换检验进行统计比较。
组内网络比较反映β多样性变化
在结肠镜检查组中，网络相似性与β多样性的变化基本一致。基线网络与POD0网络在统计学上无法区分（ARI = 0.007, P = 0.814），表明拓扑结构存在差异。然而，POD10、POD30和POD180的网络与基线网络在统计学上相似，且相似性随时间增加，这与之前观察到的α多样性恢复情况相符。虽然网络拓扑参数在所有网络中无显著差异，但基线网络与POD0或POD10之间，中心节点的度中心性、接近中心性和特征向量中心性的相似性显著低于随机预期，这可能反映了群落内部连接的重组。
在非切除术组中，β多样性在基线与POD10和POD30之间存在显著差异。与结肠镜组类似，基线网络与POD0网络在统计上随机相似。不同于结肠镜组，POD10和POD30的网络与基线相比仍为随机相似，这支持了β多样性的分析结果。基线网络与POD180网络在统计学上相似。在拓扑参数方面，POD10时的网络模块性相对于基线显著增加。在中心节点方面，POD10和POD30时的中介中心性，以及POD10时的接近中心性和特征向量中心性，与基线的差异小于随机预期。
在切除术组中，基线群落组成与所有后续时间点均存在显著差异。出乎意料的是，所有网络在拓扑结构上与基线均具有统计学相似性。在拓扑参数上，POD10时的正相关边百分比高于基线，POD30时的聚类系数下降而模块性增加，这个时间点也对应了β多样性变化的拐点。在中心节点方面，POD10时的中介中心性与基线的差异小于随机预期。
组间比较揭示干预差异
在基线时，三组间的β多样性无差异。结肠镜检查组的网络与非切除术组和切除术组的网络在统计学上相似，但两个手术组之间的网络没有可检测的统计相似性。
在POD0时，结肠镜检查组与两个手术组的网络在统计上没有可检测的相似性，两个手术组之间的网络也呈随机相似。切除术组的网络组件数量显著多于结肠镜组，正相关边百分比更低。
在POD10时，结肠镜检查组网络与两个手术组网络在统计学上相似，两个手术组网络之间也相似。与POD0相反，切除术组的正相关边百分比显著高于结肠镜组。多组中心节点间的中心性度量也显示出多种小于或大于随机预期的差异。
在POD30时，结肠镜检查组网络与两个手术组网络相似，两个手术组网络之间也相似。切除术组网络相较于结肠镜组网络具有更高的模块性，但聚类系数和边密度显著更低。
在POD180时，结肠镜检查组与非切除术组的网络在拓扑上相似，但切除术组网络与这两者均无统计学相似性。
枢纽物种的识别与意义
网络分析识别出的枢纽物种并不总是与相对丰度最高的属相对应。在结肠镜检查组，拟杆菌目成员在术后恢复期成为枢纽。在非切除术组，双歧杆菌、乳酸杆菌等共生菌在POD10和POD30被识别为枢纽物种。在切除术组，低丰度成员如Ruthenibacterium和Eggerthella在POD0成为枢纽，而Faecalibacterium和毛螺菌科物种在恢复期成为枢纽。尤其值得注意的是，较少被报道的属Enterocloster和Ruthenibacter在手术恢复相关的时间点被识别为枢纽物种，暗示它们在伤口愈合中可能扮演着新的角色。
研究者特别关注了与手术部位感染相关的链球菌属和肠球菌属。尽管链球菌属是整个研究中多个模块的主要成员，但从未被识别为枢纽物种。有趣的是，研究发现链球菌与Blautia始终存在正相关。Blautia的丰度此前被报道可预测更高的瓜氨酸浓度，并显示出对耐万古霉素肠球菌的定植抵抗力。这一发现提示，预防手术部位感染的优化策略或许应侧重于增加能抵抗感染的共生菌，而非使用抗生素清除病原体及共生菌群。
讨论与结论
本研究通过结合网络拓扑参数的统计评估，扩展了关于微生物组稳定性的理论框架。结果显示，在结肠镜检查组和非切除术组中，网络相似性通常与β多样性的差异相匹配。中心性度量的不相似可能反映了网络中不同类群相对重要性的变化。在两组手术队列中，模块性在术后显著增加，这可能反映了与伤口愈合相关的功能分化。尽管切除术组存在显著的β多样性差异，但其所有网络在拓扑结构上均与基线相似，这可能反映了维持基本生态系统过程所必需的组织结构。
研究认为，共现网络分析是探究微生物群落内部动态的补充工具。边密度可能是群落稳定性的拓扑替代指标，而模块性在术后伤口愈合期间增加，暗示了微生物组的功能多样化。研究还识别出Ruthenibacterium和Enterocloster作为枢纽属，可能在伤口愈合和术后恢复中扮演着此前未被预料和探索的角色。链球菌、肠球菌与Blautia的正相关关联也值得进一步探究。
需要谨慎的是，网络特性基于相关性，可能无法准确反映潜在的生态动力学。因此，阐明这些动态的机制需要专门的体外和体内研究。总体而言，这项工作为开始定量利用基于网络的方法研究微生物生态学迈出了基础性的一步。
