编辑推荐：
尽管益生菌和活体生物治疗产品具有强有力的机制原理并被广泛使用，但其临床获益常常不一致。研究人员提出，这种异质性反映了给药时肠道的生态状态，此时引入的菌株必须入侵一个已被占据的生态系统，该生态系统由宿主理化过滤器和常驻群落相互作用塑造而成。研究人员引入微生物生态
尽管益生菌和活体生物治疗产品具有强有力的机制原理并被广泛使用，但其临床获益常常不一致。研究人员提出，这种异质性反映了给药时肠道的生态状态，此时引入的菌株必须入侵一个已被占据的生态系统，该生态系统由宿主理化过滤器和常驻群落相互作用塑造而成。研究人员引入微生物生态位张力(microbial niche tension)作为一个可试验操作的状态变量，用以捕捉由排斥力（资源竞争、拮抗作用、空间拥挤以及宿主免疫和化学屏障）和促进力（生态位空缺、交叉喂养以及短暂的可允许窗口）共同决定的定植净阻抗。研究人员进一步定义了生态位张力指数(niche tension index, NTI)，作为一个定量框架，用于跨空间、时间和宿主背景评估定植可允许性。研究人员综合了证据，表明生态位张力受肠道转运时间(gut transit)、pH值、氧梯度、饮食、抗生素以及包括性别和激素状态在内的宿主因素所结构化。最后，研究人员概述了精准益生菌的设计规则，包括生态位预处理、合作性菌群构建、性别意识分层以及数据驱动的预测工具，以提高可重复性和临床疗效。
引言
人类肠道微生物群(gut microbiota, GM)通过与脑、肝、肺等多个器官形成轴，在维持健康方面发挥关键作用。保持其健康和正常功能依赖于GM发酵营养物质、产生必需维生素以及协助训练免疫系统的能力。对GM作用的认识催生了操纵肠道生态系统以改善多种疾病治疗结局的可能性，例如癌症、糖尿病、肥胖、抑郁症、多囊卵巢综合征、心力衰竭以及昼夜节律下的应激反应。益生菌(probiotics)是微生物组操纵最突出的领域之一，定义为“当给予足够数量时，能给宿主带来健康益处的活微生物”，被广泛营销和消费，声称有益于消化健康到降低早产儿死亡率、预防急性腹泻以及改善癌症治疗等。然而，益生菌在对照试验中的总体记录却充满矛盾，多项研究显示其在原本有希望的适应症中无完全疗效、效果甚微或无效果。新的临床指南反映了这一差异——美国胃肠病学会近期建议不要在临床试验之外常规使用益生菌治疗溃疡性结肠炎、克罗恩病或肠易激综合征，理由是缺乏高质量证据且结果不一致。因此，益生菌的前景虽引人注目，但目前仅部分实现，期望与现实之间的差距提出了一个根本性问题：为什么益生菌在某些情况下有效，而在其他情况下却失败？这一难题的关键可能在于肠道的动态生态学。当服用益生菌时，实际上是将新物种引入一个已被占据的生态系统，这是一个能够抵抗入侵者的活跃群落，这种现象称为定植抵抗(colonization resistance)。一个健康、平衡的微生物组擅长排除新来者——这对宿主是一种进化优势，因为它能阻止病原体定植，但这同时也使益生菌干预复杂化。抗生素扰动生动地说明了这一点：抗生素破坏原生菌群并降低多样性后，定植抵抗减弱，益生菌（或机会主义者）更容易立足。相反，在未受扰动、高多样性的肠道中，即使被认为有益的菌株也往往无法找到开放的生态位长期栖息。此外，益生菌定植中一个经常被忽视的因素是性别的作用、家庭成员共同生活、不同药物（包括益生菌）的服用以及饲养宠物，这些也会改变个人GM。男性和女性在众多生理和免疫参数上存在差异，新兴证据表明这些差异延伸到了肠道微生物组的组成和动态。肠道细菌可以代谢性激素并调节其可用性，表明内分泌系统和微生物组之间存在双向串扰。此外，关于性别转换过程中GM变化的新见解正在涌现。这些差异可能意味着女性肠道的生态景观与男性肠道并不相同，影响着益生菌在每个性别内的竞争或合作方式。然而，迄今为止大多数益生菌试验和研究尚未按性别分层结果或检查性别特异性结局。因此，需要一种对性别敏感的益生菌定植生态学，认识到激素、免疫反应乃至典型饮食或药物暴露等因素可能会调节特定益生菌的成功率和效果。此外，肠道内的高微生物多样性和功能冗余意味着大多数营养和空间生态位已被填满；进入的细菌通常找不到空置的生态位，并在营养和附着位点竞争中被淘汰。在益生菌背景下，同样的效应可能导致定植不良——益生菌存活通过转运并到达肠道，但随后在很大程度上未能建立立足点，因为固有微生物和宿主因素创造了一个敌对的竞争环境。而且，给定益生菌菌株稳定定植的能力似乎取决于宿主及其微生物组的个体化特征。例如，一项针对代谢综合征成人的益生菌混合物对照试验发现，只有一部分个体——那些具有特定预先存在的微生物组和饮食特征的个体——获得了代谢益处，而其他人（具有不同肠道生态）则未见改善甚至出现生物标志物的不利变化。这些结果强调，益生菌的使用需要对患者概况和肠道概况进行更详细的分析才能使其存活或发挥作用。此外，肠道环境本身施加了严格的过滤器——从胃酸和胆汁盐的严峻考验到肠道中的防御素(defensins)和免疫球蛋白A(IgA)——大多数吞服的益生菌无法存活。那些存活下来的还必须应对原生微生物群的竞争策略——争夺营养和附着位点、产生抗菌化合物以及调节局部pH值和氧气水平。所有这些因素共同创造了肠道中所谓的“微生物生态位张力”。这个术语并非要确立定植抵抗的概念，而是对其加以补充。定植抵抗通常被视为常驻微生物群和宿主屏障的一个涌现结果——即观察到的肠道生态系统阻止外来微生物建立的能力。相比之下，生态入侵理论描述了从到达扩散的过程，受繁殖体压力(propagule pressure)、非生物过滤(abiotic filtering)和生物相互作用(biotic interactions)支配，但它很少以一种对活体生物治疗剂立即“可试验操作”的方式书写。这一框架与基础微生物入侵生态学一致。土壤微观世界的实验工作表明，常驻微生物多样性可以决定对引入的细菌病原体的入侵抵抗，支持了已填充生态位和功能冗余限制建立的观念。概念性和综述性的微生物入侵框架进一步将入侵定义为一个涉及引入、建立、生长、扩散和影响的分阶段过程，受扩散或繁殖体压力、环境过滤选择、生态漂变、多样化以及常驻群落性质支配。人类益生菌研究进一步表明成功高度依赖情境，存在“个性化”的黏膜定植抵抗以及个体间不同的抗生素后恢复轨迹。基于这些基础，研究人员引入微生物生态位张力作为一个状态变量——即在特定时间和肠道位置，引入菌株遇到的净生态阻抗——整合(i)排斥力（竞争、拮抗、宿主理化/免疫过滤器）和(ii)促进力（生态位空缺、交叉喂养、支持性饮食/扰动窗口），它们共同决定定植是可允许的还是抵抗的，这与生态定植框架以及微生物组干预后人菌株级定植决定因素一致。因此，生态位张力指的是益生菌菌株试图在人体肠道中占据生态位时所遇到的多方面生态阻力。它概括了引入微生物的定植驱动力与反对它的各种力量（常驻微生物加宿主介导因素）之间持续拉锯战的概念。鉴于这些观点，本综述旨在通过生态位张力的视角，详细审视益生菌何时、何地以及对谁可能有效。通过采取全面的、个性化的、性别特异性的生态学方法，旨在为未来几年使益生菌干预更加可预测、个性化和有效提供路线图。
绘制人体肠道中的生态位张力
微生物生态位张力作为定植抵抗和定植的决定因素
微生物生态位张力是人类相关微生物组中微生物间竞争性和合作性相互作用的动态生态平衡，其中有限的空间、营养和宿主衍生因素驱动冲突（如竞争性排斥、代谢拮抗）和共存（如互惠共生、交叉喂养），从而塑造群落结构。这种张力构成了肠道定植抵抗的基础——常驻微生物群通过消耗共享资源、产生抑制性化合物和启动宿主防御来抵御入侵者的能力。一个多样化、填充良好的群落倾向于以功能冗余占据可用生态位，为新来者留下很少空间。事实上，数据驱动模型表明，微生物组的初始组成可以预测引入物种是否会建立以及以何种丰度存在。一致地，实验发现单一益生菌菌株通常提供最小保护，而由许多共栖菌组成的联合体可以将病原体抑制数个数量级。这些发现指向一个可测量的竞争能力“张力指数”，由(i)物种多样性/丰富度、(ii)功能冗余和(iii)与潜在入侵者的资源重叠等因素驱动——这些参数反映了生态位被填充的紧密程度以及固有者抵抗生态新来者的力度。这些参数反映了早期入侵生态学研究中描述的微生物可入侵性的核心决定因素：常驻群落多样性和均匀度、资源可用性、入侵者供应和环境过滤。从这个意义上说，NTI并未取代定植抵抗或入侵生态学，而是将其核心决定因素转化为可在给药前估计的可试验测量变量。为了使这个概念可测量，研究人员将NTI定义为作用于给定宿主中候选菌株的排斥力E和促进力F之间的归一化平衡：NTI = (E - F) / (E + F + ε) ∈ [-1, 1]，其中ε是一个小的正常数，添加到分母以防止除以零（或不稳定值）。NTI ≈ +1表示高张力、抗入侵的生态位；NTI ≈ 0表示有争议/中性的生态位；NTI ≈ -1表示可允许的低张力生态位。在实践中，E可以通过可试验测量的代理参数化，如群落占有率/多样性、功能冗余、资源重叠和宿主过滤代谢物（如胆汁酸谱），这些都与弹性/定植屏障和定植抵抗有经验联系。相反，F捕捉生态位空缺和促进作用（包括交叉喂养和“优先效应”，即早期定植者可以抢占或修改生态位以抑制——有时促进——后来的到达者），支持从生态组装理论到活体生物治疗产品/益生菌试验中的应答者分层和生态位预处理的定量桥梁。NTI可以通过逻辑回归模型与定植结果联系起来：P(engraft) = 1 / (1 + exp(α · NTI + β))，其中e是欧拉数约2.718，α是拟合斜率/缩放参数（α > 0时，更高的NTI降低P），β是拟合截距。一个最小的、可试验测量的代理集用于估计排斥性E和促进性F组分，应捕捉群落占有率/竞争压力、代谢可替代性、资源景观兼容性和宿主过滤。具体而言，E可以从基线α-多样性和均匀度作为占有率指标、引入菌株最近分类和/或功能邻居的丰度作为生态位抢占的代理，以及从宏基因组或组成特征推断的群落规模生长潜力指标来近似。额外的排斥压力可通过功能冗余（操作为与s预测生态位重叠的代谢途径冗余度）和资源重叠（定义为s预测底物利用谱与群落代谢能力的重叠）来量化。宿主衍生的约束贡献于E应包括胆汁酸谱、粪便pH值、肠道转运时间的代理以及反映免疫和理化过滤的炎症标志物。相反，F可以使用生态位空缺指标（低丰度或缺乏竞争类群）、交叉喂养潜力（匹配s需求的代谢物/底物可用性）、支持s的配对益生元或膳食底物的可用性以及扰动时机（如暂时增加定植可允许性的抗生素后窗口）来估计。还应指出，空间组织进一步调节生态位张力。肠道微生物在肠道各段分布不均；相反，不同的区域和微栖息地拥有独特的群落。例如，体内肠道图谱显示，中结肠的一段独特受免疫信号调控，形成了一个与其上方或下方截然不同的封闭微生物群落。此外，对人类粪便的高分辨率空间采样识别出了微米级的“空间枢纽”——反复出现在个体间的共现细菌稳定簇（通常以拟杆菌科、瘤胃球菌科/毛螺菌科为主）。当环境发生变化（如在饮食中添加可发酵纤维菊粉(inulin)）时，这些微尺度群落会重组并形成新的局部伙伴关系，一旦扰动结束便恢复原状。时间因素也塑造着生态位张力。生命早期，肠道微生物组经历可预测的演替：如双歧杆菌属（由乳汁滋养）在婴儿期占主导，然后随着饮食多样化让位于厌氧菌（如粪杆菌属、毛螺菌科）。破坏这一序列（通过剖腹产、抗生素等）可能扰乱正常的组装，而向婴儿接种关键种（如双歧杆菌或拟杆菌）则将群落引导至平衡、有弹性的状态。这些优先效应——早期到达者影响后来者——意味着时机至关重要。例如，一个次优的定植者（长双歧杆菌）可能最初占据一个生态位，但如果存在正确的营养（母乳低聚糖），后来引入的一个更适应的物种（婴儿双歧杆菌）可以超越它。同样，在成人中，引入一个新的有益菌株通常只有在生态窗口开放时才成功（如抗生素清除后）；否则，根深蒂固的微生物群会排斥入侵者。因此，人体肠道中的生态位张力源于微生物间的竞争性排斥和合作性互动，随空间和时间变化，从根本上决定了群落稳定性和宿主对扰动的抵抗力。重要的是，最近的微生物组干预研究表明，定植不应被解释为临床结局的直接替代指标。例如，尽管在溃疡性结肠炎治疗期间，较高的供体和受体α-多样性与丁酸盐相关类群的富集与更好的反应相关，但供体微生物群定植本身不能清晰地与临床缓解联系起来，这表明转移群落的生态和功能质量可能比供体获取率更重要。这一警示得到了其他研究的加强，该研究表明，除非包含安慰剂测序对照，否则粪菌移植(FMT)后的表观定植可能被大幅高估，表明一些“定植”信号可能反映的是分析噪声而非生物学上有意义的定植。一致地，在转移性肾细胞癌的随机2期TACITO试验中，特定菌株的获得或丢失（而非总定植）与主要临床终点相关，这表明结局取决于哪些功能和分类群被转移，而非单纯的批量定植。同样，PERFORM试验将抗炎类群和代谢功能的持久定植与降低的免疫相关毒性和改善的反应联系起来，而供体来源促炎酶的富集则与毒性相关，再次表明定植根据其功能载荷可能是有益的、中性的或有害的。此外，最近一项更仔细研究粪菌移植联合免疫疗法效果的癌症研究表明，通过消除致病微生物而不是引入新菌株可以获得更好的结果。总之，这些数据支持这样的观点：有临床意义的定植不仅仅是在粪便中检测到供体，而是持久的、空间相关的、菌株和功能特异性的生态位占用。
生态位张力与定植抵抗及入侵生态学术语的区别
在经典术语中，定植抵抗是已建立的肠道微生物群（和宿主屏障）阻止外来微生物立足的涌现现象。它是事后观察到的（入侵者成功还是失败？），并且不具体说明原因——可能是由于生态位被占据、抗菌代谢物或免疫因素。相比之下，生态位张力在这里被定义为一种定植前的生态状态：在特定肠道部位和时间，阻碍入侵者的因素（如资源重叠、拮抗代谢物、胆汁/氧气/免疫过滤）与促成入侵者的因素（如空置生态位、交叉喂养机会或近期抗生素诱导的“可允许窗口”）之间的平衡。一个实际困难是，许多提出的生态位张力代理通常在接触产品后才测量，此时它们可能已经反映了活体生物治疗活性、失败的定植或生态位重塑，而非宿主生态系统的真正可允许状态。这个问题尤其重要，因为益生菌在人类中的定植高度个体化，取决于宿主-微生物组基线特征，而粪便排出不一定能捕捉到黏膜定植抵抗。因此，生态位张力应在操作上定义为首次接触活体生物治疗产品前测量的暴露前状态变量。人体试验应纳入一个基线导入期，进行重复的给药前采样，理想情况下包括给药前几天收集的粪便，因为肠道群落结构和代谢物谱表现出显著的时间变异性，并受粪便稠度、粪便含水量和肠道转运时间的影响。药物和饮食元数据应前瞻性地捕获，因为质子泵抑制剂、泻药、抗生素和二甲双胍是肠道微生物组组成和功能最强的药物相关修饰因子之一。基线分析应结合菌株分辨宏基因组学与靶向代谢组学、胆汁酸谱、短链脂肪酸(SCFA)定量、粪便pH值、含水量和炎症标志物如钙卫蛋白(calprotectin)。在可行的情况下，直肠拭子、黏膜活检或管腔抽吸物应补充粪便采样，因为黏膜和粪便群落可能不同，且黏膜采样可能更好地反映目标生态位的定植抵抗。对于大型益生菌或活体生物治疗试验，这并不意味着每个参与者都应进行侵入性采样。一个可行的替代方案是分层采样设计，每个采样层回答不同的问题。核心采样层应包括所有参与者的重复粪便收集，因为粪便仍然是基线生态筛查、NTI估计、纵向菌株追踪和代谢组学分析最具可扩展性的基质。第二层空间验证层应应用于预先指定的代表性参与者子集，使用非侵入性可摄入采样胶囊，因为肠道微生物、蛋白质组和胆汁酸谱可能与粪便有显著差异，并可能揭示粪便样本中不可见的区域特异性生态状态。第三层机制性子研究层应为应答者/非应答者比较、高优先级验证队列或已接受临床指征内镜检查的参与者保留直肠拭子、黏膜活检或管腔抽吸物，这与粪便排出不一定能捕捉到经验性益生菌的个体化黏膜定植抵抗的证据一致。这种设计保留了多中心试验的可行性，同时测试了粪便来源的NTI估计是否反映了决定益生菌定植的黏膜和节段特异性生态位。这些给药前变量应前瞻性地建模为后期菌株持久性、产物衍生代谢物和临床反应的预测因子，而给药后相同变量的变化应作为定植的中介或后果单独分析。这种时间分离得到了粪菌移植研究的支持，这些研究显示供体菌株定植因疾病而异，受受体基线生态学和抗生素暴露影响，并与临床成功相关。它也得到代谢组学研究的支持，这些研究显示胆汁酸和SCFA谱在微生物转移后发生变化。因此，这一区分对于避免循环推理至关重要，即将定植后的生态变化错误分类为定植前的决定因素。换句话说，定植抵抗是一个结果；生态位张力是一个可量化的属性，可以在给药前测量或估计。由于定植抵抗和微生物生态位张力部分重叠，它们的区别在表1中阐明。入侵生态学提供了更广阔的视角，将微生物入侵定义为“外来微生物类型在常驻群落中的建立”。基础的微生物入侵框架将此过程描述为引入、建立、生长、扩散和影响的序列，受扩散或繁殖体压力、非生物过滤、生物相互作用、生态漂变和多样化塑造。实验性微生物入侵研究进一步表明，常驻多样性和群落结构可以决定入侵抵抗，而繁殖体压力在某些条件下可以改变入侵成功率。然而，这些框架本身并非患者可操作的：它们识别一般的生态过程，但未指定在益生菌或活体生物治疗试验中给药前应测量哪些基线宿主-微生物组变量。生态位张力明确将这些思想应用于临床试验：它强调了应指导设计的患者特异性生态变量。例如，与其将“抵抗”视为统一的，不如根据目标生态位的基线占有率（通过现有相似菌株的丰度衡量）对参与者进行分层或选择，并检查底物可用性（膳食纤维、聚糖）或近期扰动（抗生素、疾病）。这些区别改变了试验策略。只有当生态位未被完全占据或其他屏障较低时，增加剂量或添加菌株才有帮助。三个最近的试验说明了生态位张力如何比传统术语更好地解释结果。在一项为期12周的巴氏灭菌阿克曼菌(Akkermansia muciniphila, A. muciniphila AKK-WST01)在肥胖2型糖尿病患者中的随机对照试验(RCT)中，整个队列与安慰剂相比无显著差异。关键在于，在基线阿克曼菌低的患者中，益生菌确实定植了，这些患者的体重、脂肪量和HbA1c显著降低；而在基线阿克曼菌高的患者中，益生菌定植差且无代谢益处。因此，基线阿克曼菌占有率——而非剂量——决定了定植成功和临床效果。在一项为期2个月的对120名服用非布司他的痛风患者的试验中，联合使用多菌株益生菌(Probio-X)平均显著降低了血清尿酸和痛风发作。然而，事后分析显示，只有“益生菌应答”亚组实现了这些改善；“非应答者”亚组与安慰剂无区别。研究人员甚至构建了一个宏基因组分类器，表明基线肠道微生物组组成可以预测谁会应答。换句话说，患者间的生态差异——而非益生菌依从性——决定了定植和效果。此外，在3期新生儿重症监护病房(NICU)试验中，早产儿被随机分配接受为期28天的多物种益生菌（包括婴儿双歧杆菌BB-02）或安慰剂。益生菌将微生物群转向更“益生”的轮廓，但未能减少多重耐药菌的定植（主要终点）。令人震惊的是，宏基因组分析发现相同的婴儿双歧杆菌菌株存在于49%的安慰剂婴儿中。其存在因地点高度可变，并且与拥有接受治疗的兄弟姐妹密切相关。这意味着环境暴露和兄弟姐妹传播（NICU中的生态位可用性）——而非试验分组——驱动了该菌株的定植。传统报告可能称之为试验方案的“违反”，但生态位张力将其归因于某些医院环境中生态位的开放和高繁殖体压力。在另一项为期18周的多菌株益生菌在代谢综合征成人中的RCT中，未见整体益处。相反，益生菌组中一个“应答者”亚组显示出甘油三酯和血压的显著下降，而非应答者未见益处甚至血糖恶化。重要的是，这些应答者差异与饮食一致：应答者在基线时有独特的饮食模式和肠道微生物组。这表明益生菌需要特定的营养底物或肠道条件（生态位张力的一个方面）才能定植并发挥作用，而这只有部分参与者具备。这些案例有一个共同主题：益生菌的成功或失败不能仅用剂量、菌株或依从性来解释，还取决于患者特异性生态学。认识到这一点，生态位张力框架将引导试验筛选或调节那些生态因素：选择基线占有率低的受试者、共同施用填充所需底物的益生元或饮食、多点采样（如使用肠道胶囊）以验证定植，并使用定植/读数终点。
NTI标准化、加权和验证的分析流程
为了使NTI从一个概念框架转变为可试验操作的工具，其计算必须遵循预定义且可重现的分析流程。该流程应规定如何将原始生态学、功能和宿主来源的测量值转换为标准化的、方向一致的排斥力代理E和促进力代理F，如何对这些代理进行重新缩放、加权和组合，以及如何将所得指数与观察到的定植或临床反应结局进行验证。如果没有这样的标准化，NTI估计可能会因测量尺度、测序深度、缺失元数据、批次效应、组成约束、异常值、低生物量污染和队列特异性技术噪声的差异而失真。第一步是在分析前定义候选代理集。排斥性代理应捕捉使定植不太可能的生态屏障，包括高常驻群落占有率、高α-多样性和均匀度、引入菌株的近缘分类或功能竞争对手的丰度、代谢途径冗余度、资源重叠、抗菌基因负担、胆汁酸介导的过滤、不利的粪便pH值、快速或延迟的肠道转运、炎症和宿主免疫激活。促进性代理应捕捉增加生态位可允许性的因素，包括竞争类群的低丰度、降低的功能冗余、抗生素暴露后的生态位空缺、兼容膳食底物或益生元的可用性、交叉喂养潜力、低拮抗压力以及由扰动或疾病状态创造的短暂生态窗口。每个代理应在建模前分配预期的方向性：例如，与常驻分类群的更高资源重叠预计会增加E，而引入菌株特异性所需底物的可用性预计会增加F。在代理整合之前，原始数据需要进行质量控制。对于基于测序的变量，应使用预定义的阈值过滤样本，包括读长深度、测序质量、宿主DNA污染、阴性对照信号、批次完整性和样本处理元数据。这一步对于低生物量样本尤其重要，因为试剂污染物和交叉污染可能主导微生物信号并产生虚假的分类群-结局关联。因此，应尽可能包含阴性对照、提取空白、模拟群落和重复样本。流行率极低的分类群或功能应被移除或单独分析，因为稀疏特征可能引入不稳定性、增加多重检验负担并增加假阳性关联的风险。对于代谢组学、胆汁酸、细胞因子、pH值或转运相关变量，应在转换前应用测定特异性质量控制。低于检测限的值应根据预定义规则处理，如用检测限的一半替换、删失数据建模或多重插补，具体取决于测定方法和缺失结构。在不同分析平台上测量的变量不应直接合并，直到批次效应和平台特异性缩放被评估。缺失数据应被明确表征而非忽略。缺失应尽可能分类为完全随机缺失、随机缺失或非随机缺失。当测量的缺失本身可能反映生物学或临床状态时，应考虑非随机缺失，例如因腹泻、抗生素暴露、疾病严重程度或低生物量导致未能收集样本。缺失率非常高的变量应从主要NTI模型中排除或留作探索性分析。中等缺失率可通过多重插补、k近邻插补或基于模型的插补来处理，方法根据变量类型、分布和生物学合理性选择。重要的是，插补必须在模型开发期间的训练集内进行，以避免信息泄露到验证集。应报告每个代理的缺失比例、插补策略和完整病例敏感性分析。异常值应以结构化的方式评估。极端值可能反映真实的生物学状态，如严重炎症、非常高的胆汁酸浓度或抗生素后高度紊乱的微生物群，但也可能反映技术伪影。因此，应首先使用样本元数据、测定质量指标、重复一致性和批次信息评估异常值。技术异常值应使用预定义标准排除。生物学上合理的极端值可以保留、缩尾(winsorized)或稳健缩放，但所选方法应在敏感性分析中测试。当分布经变换后仍非正态时，使用中位数和四分位距的稳健缩放可能更可取。连续变量应在适当时进行变换以减少偏态。高度偏斜的代理，包括代谢物浓度、细胞因子和分类群丰度，可在缩放前进行对数变换。然后可使用Z评分进行标准化：zi = (xi - μ) / σ，其中xi是观测的代理值，μ是队列均值，σ是标准差。此变换将不同单位的变量置于可比尺度上，并防止数值大的变量主导指数。二元变量可编码为0/1，有序变量可进行秩变换或作为有序因子处理，有界变量可在其解释依赖于固定生物学范围时进行最小-最大缩放。如果NTI代理集包含组成性微生物组变量，如来自16S rRNA测序、鸟枪法宏基因组学或宏转录组学的相对丰度，则不应直接在相对丰度尺度上分析这些变量。由于测序数据受固定文库大小或总和约束的限制，一个分类群的明显增加必然导致其他分类群的明显减少，即使它们的绝对丰度没有变化。这种闭合效应可能产生虚假的相关性并扭曲下游的生态位占有率、功能冗余和资源重叠估计。因此，组成性特征应在纳入NTI流程前进行变换。一种常用的方法是中心对数比(centered log-ratio, CLR)变换：CLR(xj) = ln(xj / g(x))，其中g(x) = (∏j=1k xj)1/k是所有(k)个成分在样本(x)中的几何平均值。这里，xj表示成分j的丰度，如分类群、基因家族或代谢途径；k是测量成分的总数。CLR变换将每个成分表示为相对于完整组成几何平均值的比率，从而减少了单位总和约束引起的失真。由于零值无法计算对数，应在变换前使用预定义的伪计数或零替换方法处理零值。所有样本应一致应用相同的零处理规则，敏感性分析应测试不同的伪计数或零替换策略是否实质性改变NTI估计。这对于稀有分类群、低生物量样本和稀疏宏基因组特征尤为重要。预处理后，所有代理应在聚合前进行方向对齐并置于非负尺度上。每个排斥力代理的方向应设置为较高值表示更强的定植抵抗，而每个促进力代理的方向应设置为较高值表示更大的生态位可允许性。例如，密切相关竞争者的高丰度、高功能冗余和高资源重叠应增加E，而生态位空缺、兼容底物可用性或交叉喂养潜力应增加F。方向性应事先指定，以避免事后加权决策人为改善模型性能。由于Z评分变量可能取负值，标准化的代理应在组合成E和F之前转换为非负的有向分数。这可以通过最小-最大缩放、百分位秩变换或其他预定义的单调变换来实现。由此产生的有向分数可记为smE和snF，其中较高的值始终表示更强的排斥力或促进力。然后可将复合排斥力得分E和促进力得分F计算为这些非负有向分数的加权和：E = Σm=1M wmE smE，F = Σn=1N wnF snF。权重应为非负且在每个代理类别内归一化，使得Σm=1M wmE = 1且Σn=1N wnF = 1。等权可作为透明基线模型使用，但不应假定其在生物学上最优，因为不同的生态屏障可能在不同的菌株-宿主背景下占主导。例如，资源重叠可能是营养特化菌的主要屏障，而胆汁酸过滤、氧气暴露或上皮炎症可能对严格厌氧菌占主导。加权可通过三种互补策略进行。首先，当存在强机制先验知识时，可应用专家知情加权，例如当引入菌株的代谢需求已知时。其次，当目标是降维而非预测时，可使用无监督数据驱动加权。主成分分析或因子分析可将相关代理总结为正交成分，尽管这种方法可能降低生物学可解释性。第三，当定植结局可用时，可使用监督加权。逻辑回归、随机森林、梯度提升模型或其他监督学习器可识别与菌株持久性或功能性定植相关的代理。然而，非线性模型的特征重要性应解释为预测性的而非因果性的，不应被视为直接的生物学权重。当候选代理数量大、代理相关或定植事件数量有限时，惩罚回归提供了一个有用的框架。对于二元定植结局，惩罚逻辑回归最小化负对数似然加上惩罚项。LASSO应用L1惩罚，岭回归应用L2惩罚，弹性网络结合两种惩罚。弹性网络可能特别适合NTI建模，因为生态代理如α-多样性、常驻类群占有率、功能冗余和资源重叠很可能是相关的。一旦计算出E和F作为非负复合得分，NTI可定义为排斥力和促进力之间的归一化对比：NTI = (E - F) / (E + F + ε)，其中ε是一个小的正常数，防止E+F接近零时的数值不稳定性。在此约定下，正的NTI值表示排斥力占主导的抗性生态位，接近零的值表示平衡或有争议的生态位，负值表示促进力占主导的可允许生态位。ε的值应在分析前固定并在敏感性分析中评估，因为当E和F都很小时，不同的选择可能影响NTI估计。在特定生物学背景下，替代的NTI公式可能合适，但指数的方向性必须事先固定并针对观察到的定植结局进行验证。当NTI代理在不同测序批次、中心、队列、提取试剂盒、代谢组学平台或分析批次中测量时，需要技术协调。应使用阴性对照、模拟群落、重复样本、批次着色排序图和批次相关方差估计来评估批次效应。经验贝叶斯方法如ComBat可以校正批次的均值和方差特异性偏移，同时保留生物学协变量，但仅在确认生物学组别不与批次完全混淆后应用。在更简单的情况下，批次特异性中心和缩放可能就足够了。当全局分布预期在各批次间可比时，可考虑分位数归一化，但应谨慎使用，因为强制所有样本共享相同的边际分布可能去除真正的生物学差异。当数据集无法在原始数据级别协调时，可使用随机效应荟萃分析来跨队列组合代理-定植关联。NTI代理之间的共线性应在模型拟合前评估，因为许多候选变量可能捕捉重叠的生态信息。例如，高α-多样性可能与高功能冗余、高常驻类群占有率和减少的生态位空缺相关。方差膨胀因子(VIF)可计算每个代理：VIFi = 1 / (1 - Ri2)，其中Ri2是通过将代理i对所有其他候选代理进行回归获得的。VIF值显著高于5-10通常用作警告阈值，尽管此类截断值应在语境中解释而非机械应用。在高维设置中，当候选代理数量接近或超过参与者数量时，VIF估计本身可能不稳定；在这种情况下，应首选相关性聚类、降维或惩罚建模。如果检测到有问题的共线性，有几种补救措施可用：主成分回归可用正交成分替换相关预测因子；岭或弹性网络回归可稳定相关系数；变量聚类可将相关代理分组到潜在生态域中；层次建模可在域和域内级别分配权重。当假设因果关系结构时，结构方程模型可能有助于区分直接效应（如资源重叠）和由饮食、炎症或群落多样性介导的间接效应。一旦这些来自预处理、批次相关和共线性驱动的失真源被解决，NTI可以作为预测变量而非仅作为描述性生态得分进行评估。NTI与定植概率之间的关系随后可使用逻辑回归建模：P(engraftment) = 1 / (1 + exp[-(α + β · NTI)])，其中α是截距，β估计每单位NTI增加时定植概率的变化。如果NTI经过缩放使得较高值表示较强的生态抵抗，则β预期为负，因为较高的NTI应降低定植概率。然而，此关系应作为可检验假设而非指数的假定属性。在早期应用中，NTI应作为连续预测变量和分层变量进行评估，例如通过比较低、中、高张力层。这将确定NTI是否提供定植概率的单调预测因子，或者是否需要临床上有用的阈值。初步验证应使用现有的微生物组干预数据集，其中基线微生物组数据和干预后菌株级定植均可获得。合适的数据集包括具有个体化粪便或黏膜定植结局的对照益生菌定植研究、抗生素后益生菌恢复研究、具有供体-受体菌株追踪的FMT数据集，以及具有基线宏基因组学和纵向菌株检测的活体生物治疗试验。人类益生菌定植研究尤其相关，因为它们显示益生菌持久性可以高度个体化并与基线的宿主和微生物组特征相关。抗生素后恢复研究可以测试NTI是否捕捉由生态扰动创造的短暂可允许窗口。FMT菌株追踪数据集可用于评估预处理受体生态学、供体菌株丰度和供体-受体系统发育相关性是否预测定植。验证应在几个层面进行。首先，模型应测试NTI是否预测二元定植，定义为引入菌株在预定义丰度或流行率阈值以上的可检测持久性。其次，NTI应针对定量定植进行评估，如给药后菌株丰度、持久性持续时间或菌株-丰度曲线下面积。第三，当可用时，NTI应针对功能性定植进行测试，如菌株特异性代谢途径的恢复、目标代谢物的产生或宿主生物标志物的调节。最后，模型应确定NTI是否在更简单的基线特征（如α-多样性、相关分类群的基线丰度、近期抗生素暴露、粪便pH值或膳食纤维摄入量）之上增加了预测价值。理想情况下，NTI模型应在一个数据集上训练，并在另一个由不同中心、平台或干预类型生成的独立数据集上进行外部测试。模型性能应使用区分度和校准度进行评估。区分度可使用受试者工作特征曲线和曲线下面积(AUROC)进行评估。当定植不常见或定植者比例高度不平衡时，精确率-召回率曲线和曲线下面积(AUPRC)可能更具信息性。校准度应使用观测值与预测值定植图、校准截距、校准斜率和布里尔分数(Brier score)进行评估，因为模型可以正确地对参与者进行排名，但仍高估或低估绝对的定植概率。决策曲线分析可用于确定基于NTI的分层是否比治疗所有参与者、不治疗任何参与者或使用更简单的资格标准提供临床上有用的净收益。如果充分协调的外部数据集不可用，则应使用基于模拟的验证作为临时策略。可以生成具有预定义多样性、常驻类群占有率、功能冗余、资源重叠、交叉喂养潜力、胆汁酸过滤、pH值约束、抗生素扰动和批次噪声水平的合成微生物群落。然后可以在已知参数值下模拟定植结局，使研究者能够测试NTI流程是否恢复了生态抵抗与菌株建立之间的预期关系。模拟场景应变化样本量、定植流行率、效应量、测序深度、缺失率、零膨胀、代理噪声、批次效应、预测因子之间的相关性以及未测量混杂的强度。由于真实参数值在模拟中是已知的，该方法可以估计NTI模型在受控条件下的偏差、方差、校准误差、过拟合、假阳性率和稳健性。基于模拟的验证也应测试失效模式。例如，模拟应检查当一个代理类系统性缺失时、当一个批次与治疗组部分混淆时、当稀有分类群主导CLR变换的特征集时，或当NTI与定植之间的真实关系是非线性时，NTI是否保持稳定。非线性效应可使用受限三次样条、广义加性模型或基于树的模型进行评估，但这些方法在使用于试验分层前应进行外部验证。如果模拟结果表明NTI在现实的缺失、稀疏性或批次结构下不稳定，则应在前瞻性临床应用前简化代理集。未来的基于NTI的试验应包括基于模拟的功效分析，而非依赖通用的样本量规则。所需的样本量将取决于定植流行率、预期NTI效应量、候选代理数量、代理相关性、缺失率、批次效应、预期模型复杂度以及终点是二元定植、定量菌株丰度、功能性定植还是临床反应。对于预测模型开发，样本量规划应考虑候选预测参数的数量、预期结局分数和预期模型性能，而非仅依赖简单的事件每变量规则。如果样本量有限，应事先减少NTI代理集，或使用降维和惩罚建模以减少过拟合。功效模拟应估计检测NTI与定植概率之间预设关联所需的最小参与者数量，同时控制假发现率并维持可接受的校准度。对于二元定植终点，模拟应变化每单位NTI增加的预期比值比、基线定植概率、缺失数据率和候选代理数量。对于定量菌株丰度，模拟应变化预期效应量、丰度离散度、零膨胀和测量误差。对于试验分层，模拟应测试NTI定义的低张力组和高张力组在定植概率上是否存在足够差异，以证明富集、自适应随机化、生态位预处理或靶向给药窗口的合理性。最后，NTI模型在用于前瞻性试验决策之前，应经过严格的内部和外部验证。内部验证可使用自助法、重复交叉验证或训练-留存拆分来估计乐观度并调整模型参数，如惩罚回归中的惩罚参数λ。外部验证应在独立队列中测试最终锁定的模型，而不重新选择代理或重新拟合权重。在此阶段，NTI应被视为一个假设生成和试验设计框架，而非临床可操作的生物标志物。在NTI用于指导患者选择、给药窗口、生态位预处理或临床决策之前，其代理集、加权方案、终点定义和分析计划应预先注册，并在多个独立队列中进行验证，理想情况下跨越不同的中心、测序平台、干预类型和宿主人群。临床使用应要求证明锁定的NTI模型在这些队列中保持区分度、校准度和决策曲线效用，并在更简单的基线变量（如α-多样性、近期抗生素暴露、相关分类群的基线丰度、粪便pH值或膳食纤维摄入量）之上增加预测价值。报告应遵循既定的预测模型报告原则，包括清晰描述候选预测因子、缺失数据处理、模型规格、内部验证、外部验证和校准度。
