结肠癌微生物组的生态驱动因素：确定性生态位选择优于随机定植的证据

《Frontiers in Oncology》：Ecological drivers of the colon cancer microbiome: evidence for niche-based selection over stochastic colonization

　　

结肠癌是一种受遗传、环境和微生物因素影响的全球主要恶性肿瘤。本研究整合了生态学理论，包括Hubbell的中性理论，并由Ning等人（2019，PNAS）的随机性分析框架以及经典生态位理论加以补充，以评估随机定植和确定性生态位选择在结肠癌组织中塑造微生物组组装和动态变化的贡献。研究人员利用Levins生态位宽度分析鉴定了专性物种（specialist species）。通过重新分析Nejman等人的数据集（包括22例正常邻近组织（NT）和22例原发肿瘤（PT）样本），研究发现：（i）确定性生态位选择似乎在NT和PT微生物组的群落组装和多样性维持中占主导地位，归一化随机性比值（NSR）显示NT为0.267，PT为0.287，均低于0.5（Wilcoxon检验PT>NT，p=0.027）。仅有18.2%的NT样本和27.3%的PT样本符合中性模型，证实随机漂变作用有限，生态位选择在组织微生物组结构中占主导地位。（ii）通过生态位宽度分析，研究人员鉴定出106种专性物种，其中14种是已知的结直肠癌（CRC）相关类群，包括具核梭杆菌（Fusobacterium nucleatum）、大肠杆菌（Escherichia coli）、粪肠球菌（Enterococcus faecalis）和厌氧消化链球菌（Peptostreptococcus anaerobius）。这些物种在炎症、基因毒性、免疫调节和化疗耐药中具有明确的作用。（iii）这些发现表明，确定性生态位选择而非随机定植似乎主导着结肠癌组织中的微生物群落组装。（iv）研究人员推测，肿瘤微环境施加了生化过滤条件——包括缺氧、乳酸积累和pH值改变——从而筛选出能够适应这些条件的专性物种。（v）尽管样本量有限且需要在独立队列中进行验证，但本研究展示了一个用于理解结肠癌微生物组动态变化的整合框架。研究人员为未来研究提出了可检验的假设：（i）专性物种丰度可能与肿瘤分期或患者生存率相关；（2）在临床前模型中，专性物种的耗竭可能改变肿瘤进展；（3）混合中性-生态位模型可能更适合较大的纵向数据集。通过将生态学理论与微生物组研究相结合，本研究为理解驱动肿瘤组织中微生物群落组装的力提供了见解，并突出了基于微生物组的诊断和治疗方法的潜在靶点。

研究背景与科学问题

结直肠癌（CRC）是全球最常见和最致命的恶性肿瘤之一，其发生发展受遗传、环境和微生物因素的复杂相互作用驱动。肠道微生物组在肿瘤起始、进展和治疗反应中扮演活跃角色，尤其是在微生物密集的结肠环境中。尽管已有大量研究揭示了CRC相关的微生物失调特征，例如促炎细菌扩增和有益丁酸产生菌减少，以及特定病原菌如Fusobacterium nucleatum和产毒素Bacteroides fragilis（ETBF）的致癌机制，但关于微生物群落在癌症背景下如何组装的基本生态学规则仍不清楚。核心争议在于：微生物群落的构建是由确定性（生态位驱动）过程主导，还是由随机（中性）过程主导。经典生态位理论认为每个物种占据独特的生态位，而中性理论则假设生态等价性，预测群落结构由随机的出生、死亡和迁移驱动。此前在其他癌症类型中的研究（如肺癌中确定性肿瘤驱动选择、乳腺癌中恶性转化过程中中性物种比例下降）提示肿瘤微环境可能施加选择性压力，但在结肠癌组织驻留微生物组中直接比较这两种过程的证据仍然有限。因此，本研究旨在应用生态学理论框架，量化结肠癌微生物组中随机性与确定性的平衡，揭示肿瘤微环境对微生物群落组装的生态学影响。

研究方法概述

研究人员使用了Nejman等人公开发表的数据集，该数据集包含来自结肠癌患者的22例正常邻近组织（NT）和22例原发肿瘤（PT）的微生物组样本。主要分析方法包括：采用Ning等人（2019）提出的归一化随机性比值（NSR）框架量化生态随机性程度；应用Hubbell中性理论的Etienne精确中性检验评估样本是否符合中性模型的预期；运用Levins生态位宽度（B_N）分析结合Finn等人的定量限（LOQ）标准鉴定专性物种（specialist species）和泛化物种（generalist species）。

研究结果

NSR评估正常与肿瘤微生物组组装的随机性

研究人员首先采用NSR框架量化NT和PT微生物组中确定性与随机性过程的相对贡献。结果显示，NT微生物组的NSR值为0.267，PT微生物组为0.287，两者均低于0.5的判定阈值，表明确定性过程在两种组织的微生物群落组装中均占主导地位。Wilcoxon检验显示PT的NSR显著高于NT（p=0.027），提示肿瘤组织中的随机性略有增加。

Hubbell中性理论检验

应用Etienne精确中性检验后，研究人员发现22例NT样本中仅有4例（18.2%）符合中性模型的预测，22例PT样本中有6例（27.3%）符合中性模型。绝大多数样本的观测对数似然值落在中性模型95%置信区间之外，表明随机过程单独无法完全解释结肠癌微生物组的群落组装。图1展示了代表性中性拟合样本的等级-丰度曲线，其中观测数据（红色实线）与中性模拟数据（黑色虚线）的吻合情况得到直观呈现。

生态位宽度与专性物种鉴定

鉴于中性模型和NSR分析均指向确定性过程的主导作用，研究人员进一步通过Levins生态位宽度分析鉴定专性物种。在R=2（代表NT和PT两种资源状态）的框架下，共鉴定出106种专性物种，其特征为B_N接近0.5、校正后p<0.05且丰度高于LOQ。其中14种为已知的CRC相关类群，包括Fusobacterium nucleatum（B_N=0.604，参与炎症、NF-κB激活、化疗耐药和不良预后）、Escherichia coli（B_N=0.583，通过colibactin产生基因毒性DNA损伤）、Enterococcus faecalis（B_N=0.567，引起炎症和氧化性DNA损伤）等。值得注意的是，尽管这106种专性物种仅占NT和PT队列中超过10000个总物种的一小部分，但其累积相对丰度在某些样本中超过50%（如图2所示），突显了少数专性物种在微生物群落结构中的生态主导地位。有趣的是，大多数符合Hubbell中性模型的样本显示出较低的专性物种相对丰度，暗示随机过程主导的样本中确定性选择压力较弱。

讨论与结论

在讨论部分，研究人员指出本研究表明确定性生态位选择在结肠癌微生物组群落组装中发挥重要作用，这与肺癌等其他癌症类型的研究发现一致。肿瘤微环境（TME）通过特定的生化过滤条件——包括缺氧（选择兼性厌氧菌如E. coli和F. nucleatum）、乳酸积累（偏好能代谢替代碳源的物种）和pH值改变（限制酸敏感类群）——驱动确定性选择。极低的中性模型符合率（44个样本中仅10个）和低于0.3的NSR值共同支持这一结论。此外，Sloan近中性模型在所有比较中拟合失败（R2<0），进一步证实了生态位选择的优势地位。研究人员承认研究的局限性，包括样本量较小（每组n=22）、缺乏独立验证队列、肿瘤组织低生物量带来的挑战、未完全考虑所有随机过程（如物种形成和扩散），以及未直接测量TME中的生化因子。未来工作应验证三个假设：（i）专性物种丰度与肿瘤分期或患者生存率相关；（ii）在动物模型中耗竭专性物种可改变肿瘤进展；（iii）混合中性-生态位模型更适合较大规模的纵向队列。

结论翻译：本研究表明，确定性生态位选择似乎在结肠癌微生物群落组装中占主导地位，随机过程作用有限。一小群专性物种——包括已知的CRC相关类群如F. nucleatum、E. coli和E. faecalis——驱动正常组织和肿瘤组织的群落结构。这些发现为解读肿瘤-微生物组相互作用提供了一个生态学框架，并确定了值得进一步机制研究和转化研究的候选专性物种。