《Microorganisms》：Sampling Design and Temporal Variation Constrain Inference of Human Gut Microbiome Recovery
编辑推荐：
人类肠道微生物组（gut microbiome）恢复估计取决于样本采集时间、所采用的参考状态以及群落差异的度量方式。研究人员重新分析了三个扰动队列和两个纵向比较队列的公开图谱，并在标量工作模型下评估采样设计。在12名抗生素暴露参与者中，第180天晚期相对于观察
人类肠道微生物组（gut microbiome）恢复估计取决于样本采集时间、所采用的参考状态以及群落差异的度量方式。研究人员重新分析了三个扰动队列和两个纵向比较队列的公开图谱，并在标量工作模型下评估采样设计。在12名抗生素暴露参与者中，第180天晚期相对于观察到的急性期位移比的中位数为0.874（95%参与者自举区间（bootstrap interval）0.810–0.992）。环境暴露发现取决于距离几何；23次腹泻事件后14天内随访未能确定恢复。较长采样间隔与九个Broad Institute-OpenBiome Microbiome Library（BIO-ML）供者和19名Vandeputte参与者中更大的相对组成差异相关。稀疏的个体采样方案无法识别所有自由恢复参数。在一个说明性模拟中，恢复尺度为10天、残余位移为0.20、高斯噪声（Gaussian noise）标准差为0.08，将六次随访从第30–180天重新分配到第2–180天，使估计值落在真值25%以内的比例从25.4%提高到69.0%。置信集（confidence sets）可包含真值，但仍过宽而不适于有用推断。弱信号和模型误设（model misspecification）可能使额外观测失效。这些发现支持队列特异性解释和情景依赖规划，但未建立普适恢复律或充分样本量。
研究背景与问题：人类肠道微生物组（gut microbiome）在无记录干预时也随时间变化，纵向研究已报道饮食、环境转变、旅行、肠道感染和抗生素后的个体化反应；部分维度快速回到基线，其他维度可持续改变。抗生素研究报道数周至数月内组成大幅回归，同时存在物种丢失、耐药变化和较长期暴露关联。核心问题在于：这些反应在按扰动幅度和时间尺度归一化后是否遵循共同恢复轨迹仍未解决；恢复估计又依赖采样时间、参考状态和群落差异度量。现有研究关注共享动态、丰度波动和生物动机恢复模型，但队列特异位移、滞后对比及个体恢复参数可识别性未同时评估，因此研究人员开展多队列再分析和模拟。
研究开展与意义：研究人员分别分析抗生素暴露、环境转变和旅行相关腹泻队列，并匹配分类学—抗性组（resistome）关联，进行度量和预处理敏感性分析；用两个独立观察队列BIO-ML九名供者和Vandeputte 19名参与者进行个体内滞后比较；四个重叠BIO-ML宏基因组供者和匹配定量微生物组图谱（QMP）检验处理与丰度表示；再用可识别性分析和已知参数模拟比较采样方案。结论是观察到的恢复对比取决于队列、参考状态和距离度量，相对组成差异随间隔变化；观察时间和急性参考处理影响恢复尺度估计准确性和精度；现有数据不能建立跨扰动普适恢复律或充分样本量。意义在于支持队列特异性解释和情景依赖的采样设计评估。
关键技术方法：二次分析样本来自SPIRE v1 mOTUs3重处理Palleja抗生素队列（12人）、Sun环境队列（14人）、Boolchandani腹泻队列（23人物种、18人抗微生物耐药基因（ARG））、BIO-ML（9名供者；4名重叠宏基因组供者）和Vandeputte（19人；相对微生物组图谱（RMP）与QMP），并以David Subject A作敏感性。方法包括组成性预处理与中心对数比（CLR）变换、Aitchison位移（Aitchison displacement）、Jensen–Shannon距离（Jensen–Shannon distance）、Bray–Curtis距离（Bray–Curtis distance）、参与者自举（bootstrap）置信区间、恢复比、Jacobian秩可识别性（identifiability）分析及ADEMP框架下已知参数模拟。
研究结果：
3.1 采样方案决定哪些恢复量可估计：对三种扰动队列和两个纵向队列比较基线/急性参考可用性并用Jacobian秩分析，代表性两点抗生素方案固定β时秩为2、单点方案秩为1；允许β变化时两约束对应三个个体参数，故仅部分恢复量可估计。
3.2 抗生素暴露后伴随部分分类学回归，但合并时间尺度不确定：SPIRE重处理显示12人中位Aitchison位移从急性最大2.383降至D42 2.283和D180 2.089；D180晚期/急性比0.874（95% CI 0.810–0.992；p=0.00684，BH q=0.0246），Jensen–Shannon和Bray–Curtis比0.768和0.707；条件合并τ随名义/日历时间和β变化，不识别个体恢复常数。
3.3 环境反弹取决于距离度量；腹泻早期随访不能确定恢复：14名学生T6/T3 Aitchison比0.846（CI 0.608–1.047；p=0.153），Jensen–Shannon和Bray–Curtis比0.815和0.812显著，ARG Aitchison比0.862；23例腹泻14天内物种Aitchison比1.175（CI 0.954–1.298；p=0.142），校正后无显著对比。
3.4 环境暴露后支持分类学—抗性组关联，腹泻后不支持：环境队列恢复比在Aitchison、Jensen–Shannon和Bray–Curtis下均正相关（ρ=0.908、0.912、0.873；q低至4.38×10
?5）；腹泻18人Aitchison ρ=0.350（p=0.155；q=0.206），六项关联未通过校正。
3.5 参与者水平组成差异随采样间隔增加：BIO-ML九名供者>90天/1–7天中位Aitchison距离比几何均值1.261（CI 1.196–1.332；p=0.00391），Vandeputte 19人31–90天比1.134（CI 1.088–1.192；p=1.91×10
?5）；交替14天块验证中BIO-ML短/长滞后超限7.9%/29.4%，Vandeputte为8.1%/25.8%，差值分别21.5和17.8百分点。
3.6 稀疏峰后采样使个体恢复参数不可识别：对归一化恢复族作局部可识别性分析，一个峰后观测仅给出一标量约束；固定β时兼容的τ与d
∞形成连续体，放开β增加未约束方向；r(90)=0.82示例中τ可至约453.5天，D365预测约0.447至0.82。
3.7 预处理敏感性保留队列特异性对比：九种流行率/零替换设置下抗生素D180中位比0.816–0.904，环境0.800–0.876，腹泻1.119–1.175；BIO-ML长/短比1.221–1.307，Vandeputte 1.104–1.155；留一参与者仍>1，方向稳定但不构成独立重复。
3.8 已知参数模拟量化时间安排收益与不确定性校准限度：在τ=10天、残余0.20、噪声0.08的模拟中，spread six较uniform six降低均方根对数误差（log-RMSE）90.6%，25%内估计由25.4%升至69.0%；τ=90天时spread six较差；覆盖、区间可用性和精度分离，模型误设时增加观测不能修复。
讨论部分总结：作者强调经验发现条件于队列和观察窗；抗生素支持相对采样急性参考的位移降低，环境证据随度量而变，腹泻早期随访不确定；分类学—抗性组比值仅在环境队列相关；纵向比较涉及区间依赖的相对组成；这些结果不能建立普适耦合、完全恢复或共同恢复律，样本对数量不能补偿有限生物学重复。模拟说明观察时间在指定模型下影响参数估计，快恢复时早期测量有用，慢变化需后期随访；增益取决于时间尺度、形状、信号和误差结构，覆盖率与精度可分离，程序可频繁包含真值但区间过宽。合并可估计共享参数但须假设共享曲线、时间尺度或残余状态，不能为个体增加观测时间点；两个峰后观测加三个自由参数使个体问题欠定。急性参考误差可方向相反：漏峰使表观恢复偏弱，含噪急性最大值膨胀分母使表观恢复偏强；需围绕急性窗采样、保留所有急性观测和日历时间。参考校准中端点不交叠块验证防止样本复用但相邻块仍可相关；更长训练隔离后可用贡献者过少。相对与绝对丰度解释不同，ARG家族组成、通路潜力和活性不同。局限包括队列小、非人群代表、筛查不排后续暴露、重叠子集、同日样本非技术重复、自举假设独立可交换并条件于日历和预处理、小样本校准不确定、Vandeputte无>61天、腹泻终点14天内、处理管线不统一和符号秩对称假设。模拟假设已知急性起点和单调回到残差状态，不能代表所有过冲、复发或晚期位移增加；固定β=1是估计假设，0.7/1.3为压力场景；高斯噪声不生成分类单元或测序计数；数值检查不验证机制或充分样本量。规划应始于生物学问题和可观察终点，预先定义主要度量、基线、急性窗规则、分析单位和多重性族，获取重复事件前观测，保留原始急性测量并记录实际采集时间和暴露；所需基线/急性样本数未确定。
结论部分翻译：观察到的恢复对比取决于队列、参考状态和距离度量，而相对组成差异随样本间隔变化。在所测试模型下，观察时间安排和急性参考处理影响恢复尺度估计的准确性和精度。因此，应根据明确的准确性、覆盖率和区间宽度目标来评估采样设计。现有数据不能建立跨扰动的普适恢复律。
