A Reproducible Meta-Analytic Pipeline for Quantifying Ecological Bottlenecks, Alpha Diversity Decay, and Core Biomarker Fluctuations in Ex Vivo Human Gut Microbiota Fermentations

技术摘要：用于量化离体肠道发酵中生态瓶颈的可重复元分析流水线

问题陈述
离体连续流生物反应器系统是机制微生物组研究的标准工具，允许在受控参数下评估化合物的发酵情况。然而，将复杂的、多门类的类人肠道生态系统从“体内”环境过渡到人工容器的过程，会触发即时且严重的生态瓶颈。这种转变涉及宿主粘膜保护载体的突然缺失、流体剪切应力的引入、转运时间的改变以及化学生长培养基的选择性。这些因素创造了一种选择压力，使具有韧性的广谱种（generalist species）获得优势，同时迅速排挤掉特化型的稀有分类单元（specialized, rare taxa）。本文认为，准确量化这种初始的“生物反应器漂移”（bioreactor drift）和 $\alpha$ 多样性崩溃的程度，是验证任何下游实验处理（如纤维或益生元干预）的关键前提。

方法论
本研究展示了一个基于 R 语言开发的独立、可重复的统计验证工作流，旨在分离并测量基准生物反应器漂移。该方法依赖于以下步骤：

• 数据摄取与清洗： 分析利用了来自 Bressure 等人 (2025) 补充文件中公开的宏基因组物种泛基因组 (MSP) 相对丰度矩阵。数据通过 readxl 包从原始 Excel 工作簿中摄取。流水线包含了算法清洗，例如对元数据进行空白修剪 (trimws)、动态解析样本 ID 以区分“新鲜粪便”（基准）和“48小时生物反应器”样本，以及严格执行缺失分类数据的零填充，以确保行总数计算的准确性。
• 统计框架： 使用 vegan R 包 (v2.6-4) 计算生态指数。评估了两个主要的 $\alpha$ 多样性参数：
  1. 观测分类丰富度 (Observed Taxonomic Richness)： 每个样本中独特的、非零 MSP 分类单元的绝对计数。
  2. 香农多样性指数 (Shannon Diversity Index, H')： 一种结合了丰富度和物种分布均匀性的群落熵度量。
• 定向追踪： 构建了一个特定矩阵，用于剖析核心共生菌群及广谱指标菌（包括 Escherichia coli、Bacteroides ovatus、Faecalibacterium prausnitzii [菌株 1 & 4]、Akkermansia muciniphila 和 Bifidobacterium longum）的存活情况。
• 假设检验： 鉴于微生物相对丰度的偏态分布，研究拒绝了参数检验假设。研究采用两侧非参数 Wilcoxon 秩和检验（Mann-Whitney U），通过自动化循环计算所有指标的精确 p 值和 W 统计量。

关键结果
分析证实，从宿主到容器的转变诱发了一个显著的筛选事件，其特征是群落复杂性的剧烈下降：

• $\alpha$ 多样性崩溃： 观测分类丰富度从新鲜供体粪便的基准中位数 281.00 个独特 MSP，骤降至生物反应器中的 34.00 个 MSP 简化核心。这代表了 87.9% 的丰富度崩溃 ($W = 216, p = 0.0036$)。
• 群落重构： 香农多样性指数显示出显著变化 ($W = 0, p = 0.0037$)，中位数从新鲜粪便中的 $2.21 \times 10^{-5}$ 增加到容器中的 $1.39 \times 10^{-4}$。作者将其解释为群落熵的基本重构，即稀有分类单元的丧失改变了剩余群落的均匀度。
• 目标生物标志物波动：
  • 广谱种 (Generalists)： Escherichia coli 和 Bacteroides ovatus 在生物反应器中表现出相对丰度的增加，尽管统计显著性各异 ($p=0.05$ 和 $p=0.10$，分别为后者)。
  • 特化种 (Specialists)： 高度特化的分类单元，包括 Faecalibacterium prausnitzii（两种菌株）、Akkermansia muciniphila 和 Bifidobacterium longum，在生物反应器环境中的中值相对丰度为 0.00（IQR 为 0.00）。然而，由于容器组中完全缺乏信号，Wilcoxon 检验得出的 p 值不具显著性 ($p \geq 0.10$)，这表明虽然中值暗示了消除现象，但统计检验并不支持存在完全消除所有生物学重复样本的确定性结论。

意义与主张
本文将自身定位为并非发现新的生物学机制，而是建立了稳健的开放科学框架，用于下游干预的基准测试。其主要贡献在于该可重复的流水线，它能将“基准生物反应器漂移”从实验效应中分离出来。

作者声称，生物反应器提供的是一个“功能性浓缩版的原生生态系统”，而非直接的生理镜像。研究结论指出，由于缺乏宿主载体，稀有本土物种无法生存，而具有韧性的广谱种则趋于稳定。通过量化这 88% 的分类单元损失，该流水线允许研究人员将此类固有局限性与后续实验处理的具体效应区分开来。代码和数据均已公开，以确保可重复性，并促进未来微生物组研究中生态基准评估的标准化。