《Frontiers in Microbiology》：Beyond the billion: dose-response immunophenotyping and machine learning classification of live vs. heat-treated gram-positive probiotic strains in human peripheral blood mononuclear cells—a single-donor proof-of-concept study
编辑推荐：
益生菌研究受到三个普遍存在但未受到充分质疑的假设的制约：发挥功效所需的最低菌落形成单位（CFU）为10亿、肠道定植的必要性，以及活菌制剂相对于热处理制剂具有固有优越性。本研究通过全因子实验设计解决了这些缺口，该设计在ISO 19344:2015标准下，以流式细
益生菌研究受到三个普遍存在但未受到充分质疑的假设的制约：发挥功效所需的最低菌落形成单位（CFU）为10亿、肠道定植的必要性，以及活菌制剂相对于热处理制剂具有固有优越性。本研究通过全因子实验设计解决了这些缺口，该设计在ISO 19344:2015标准下，以流式细胞术验证的三个浓度（105、106和107细胞/孔）评估了10株革兰氏阳性益生菌株的活菌（活性荧光单位，AFU）和热处理（总荧光单位，TFU）两种形式，实验对象为来自单一健康男性高加索供体（58岁）的原代人外周血单个核细胞（PBMC），并通过BioPlex悬浮阵列同时定量17种细胞因子。活菌制剂在17种分析物中的16种上诱导了显著大于热处理制剂的绝对细胞因子反应，而热处理制剂在17种细胞因子中的13种上表现出更强的剂量-反应相关性（平均株内Spearman ρ最高达1.00），研究人员将该发现归因于活菌在24小时共培养期间不受控制的增殖压缩了有效浓度范围。10株活菌制剂中有6株呈现单调递增曲线；两株表现出非单调的钟形动力学，在106 AFU/孔时活性达到峰值，并在107时显著衰减，这直接否证了剂量递增均匀增加免疫学活性的假设。MCP-1是配对分析中唯一在活力状态之间无显著差异的细胞因子（q = 0.98，几何均数比0.99，95% CI 0.89–1.12，排除了任一方向上约11%以上的任何差异），提供了内部方法学对照。在该单供体模型中，无监督层次聚类识别出三个免疫表型簇，这些分组在被视为可推广的生物学表型之前需要多供体验证，随机森林分类实现了83.3%的内部划分恢复一致性（95% CI 66.4–92.7%；簇源自相同数据；并非对新型菌株泛化能力的估计），而随机水平为33.3%；在该单供体实验中，IL-12p70、IL-13和IFN-γ（而非IL-6或IL-1β）是菌株身份的主要区分因子；该排序的泛化性需要多供体验证。热处理制剂在107 TFU/孔时相对于活菌制剂在140个可评估的菌株×分析物组合中的80个中实现了≥70%的功能等效性，占可评估组合的57.1%（功能等效剂量，FED70），而一株菌株在热处理形式下在所有浓度均表现出强烈减弱的活性，该发现受限于以下注意事项：实验设计中未包含阳性对照刺激物来正式验证实验当日PBMC的功能活性。这些发现建立了一个方法学框架，整合了流式细胞术标准化、多重免疫表型分析和机器学习，用于益生菌研究中的循证剂量表征、后生元功能等效性评估和数据驱动的菌株分类（所有p值均为单供体实验背景下的描述性指标）。
论文解读
一、研究背景与问题
益生菌领域长期受制于三个未被充分质疑的核心假设：其一，益生菌发挥功效需要至少10亿菌落形成单位（CFU）的最低剂量；其二，肠道定植是益生菌起效的必要前提；其三，活菌制剂在免疫调节效果上天然优于热处理（灭活）制剂。然而，这些假设缺乏系统性的实验验证。与此同时，该领域还存在一个深层的计量学危机——传统的CFU计数法仅能反映可培养的活菌数量，无法准确量化热处理后仍保留结构完整性的细菌细胞，导致活菌与灭活制剂之间缺乏可比性。此外，益生菌免疫调节研究普遍存在单剂量设计、缺乏剂量-反应关系证据等问题，Christensen等人（2006年）在71名受试者的多剂量人体研究中未发现任何剂量水平的显著免疫变化，这一阴性结果可能正反映了CFU计数法带来的剂量不确定性。基于上述背景，研究人员设计了这项单供体概念验证研究，旨在回答四个核心问题：同一菌株的活菌与热处理制剂在匹配浓度下是否存在可测量的免疫调节差异；在明确的浓度范围内能否证明剂量-反应关系的存在及其形状是否因菌株而异；如何测量细菌剂量以使活菌与热处理制剂具有可比性；机器学习能否从多重细胞因子数据中提取单变量分析无法揭示的菌株水平结构。
二、研究设计与技术方法
该研究采用全因子实验设计，使用来自Lonza公司的单一健康男性高加索供体（58岁）的冻存外周血单个核细胞（PBMC）作为细胞模型。研究人员将10株革兰氏阳性益生菌株（编码为A至L，排除J和K）分别制备为活菌和热处理两种形式，每种形式设置105、106和107细胞/孔三个浓度，每个条件设置3个技术重复，共180个处理观察值，并以未刺激PBMC培养物作为阴性对照。关键技术方法包括：依据ISO 19344:2015标准，采用硫唑橙（TO）/碘化丙啶（PI）双荧光染色流式细胞术进行细菌计数，活菌计为活性荧光单位（AFU），热处理菌计为总荧光单位（TFU），从而实现两种制剂在相同物理单位下的浓度匹配比较；热处理采用85°C、60分钟标准化水浴方案；共培养24小时后，采用Bio-Plex Pro Human Cytokine 17-plex Assay（Bio-Rad）同时定量17种细胞因子；数据分析采用Wilcoxon符号秩检验（配对比较）、Spearman秩相关（剂量-反应分析）、主成分分析（PCA）、K-means和层次聚类、随机森林（RF）等机器学习方法，以及基于Jeffreys无信息先验的贝叶斯估计框架。
三、主要研究结果
3.1 活菌诱导的细胞因子反应显著大于热处理制剂
配对分析显示，17种分析物中有16种在活菌与热处理条件之间存在显著差异，其中TNF-α的差异最大（几何均数比310.8倍，95% CI 36.4–2,653），其次为MIP-1β（69.4倍）、IL-8（25.8倍）、G-CSF（15.5倍）和IL-1β（13.3倍）。所有16种显著差异的细胞因子均同时超过统计学阈值（FDR < 0.05）和生物学意义阈值（1.5倍）。值得注意的是，MCP-1是唯一在两种活力状态间无显著差异的分析物（q = 0.98，几何均数比0.99，95% CI 0.89–1.12），这一有界等价性排除了任一方向约11%以上的差异，构成内部方法学对照。
3.2 剂量-反应关系在活力状态间存在定性差异
热处理制剂在17种细胞因子中的13种上表现出更强的剂量-反应相关性（平均株内Spearman ρ最高达1.00），而活菌制剂尽管在每个浓度均诱导了更高的绝对细胞因子水平，但大多数分析物的剂量-反应相关性较弱。研究人员将该不对称性归因于活菌在24小时共培养期间不受控制的增殖——即使最低接种量（105 AFU/孔）也可能在实验期间大幅扩增，压缩了三个名义剂量水平之间的有效差异。Kruskal-Wallis检验显示热处理条件下17种细胞因子中有14种存在整体浓度效应。
3.3 菌株特异性剂量-反应谱存在显著异质性
菌株A和B在两种活力状态下对大多数细胞因子均表现出单调递增反应（ρ = +1.0）；菌株E在活菌条件下几乎缺乏剂量-反应关系（平均ρ = 0.12），提示其免疫调节效应在较低浓度即已饱和。跨整个面板，热处理制剂中17种分析物有15种达到浓度-反应相关显著性，而活菌制剂中仅为9种。
3.4 非单调剂量-反应曲线：钟形动力学
菌株D和E表现出非单调的钟形剂量-反应动力学，在106 AFU/孔时活性达到峰值，在107时显著衰减（菌株E的G-CSF减少83%，IL-6减少25%）。这一发现直接否证了剂量递增均匀增加免疫学活性的假设，与Hua等人（2010年）提出的毒物兴奋效应概念一致。
3.5 菌株C：强活力依赖性免疫调节谱
菌株C在活菌状态下具有活性，但热处理后活性强烈减弱，在所有浓度下均表现出热处理形式的显著衰减，是唯一一株具有该特征的菌株。
3.6 活菌与热处理的功能等效性：FED70矩阵
研究人员引入功能等效剂量（FED70）概念，定义为热处理制剂达到对应活菌制剂在106 AFU/孔参考剂量下免疫学活性≥70%所需的最低TFU/孔浓度。在140个可评估的菌株×分析物组合中，80个（57.1%）在107 TFU/孔时达到FED70。部分菌株在特定分析物上表现出显著的功能趋同：菌株A在107 TFU/孔时产生相当于活菌106 AFU/孔参考IL-6活性的266%；菌株E热处理制剂在107时超过活菌106参考水平（105%）。
3.7 细胞因子间相关结构由全局反应幅度主导
全数据集Spearman相关分析显示，IL-2与IL-5（ρ = 0.97）、IL-4与IL-5（ρ = 0.95）、IL-2与IL-4（ρ = 0.94）之间存在近单位相关。这种全数据集水平的生态学相关主要由活菌制剂高剂量条件下的全局反应幅度驱动，而非细胞因子对之间的直接转录共调控。
3.8 数据驱动的菌株特异性免疫信号特征
对60个条件均值矩阵的PCA分析显示，前两个主成分解释了85.2%的总方差（PC1：75.5%，PC2：9.7%）。无监督层次聚类（Ward连接，欧氏距离）识别出三个功能簇：簇1（菌株A、B）——高幅度、广谱促炎细胞因子诱导；簇2（菌株E、F、G、H、I、L）——中等反应且动力学谱多样；簇3（菌株C、D）——低幅度活菌反应。随机森林10菌株分类（留一法交叉验证，LOO-CV）达到40.0%的准确率（95% CI 24.6–57.7%），约为10%随机水平的四倍。特征重要性分析确定IL-12p70（重要性0.110）、IL-13（0.103）、MCP-1（0.090）、IFN-γ（0.090）和IL-7（0.079）为菌株分类的主要区分因子，这些因子与绝对倍数变化最高的细胞因子（IL-6、IL-1β、G-CSF）不同，表明菌株身份信息存在于免疫极化方向而非激活幅度。3簇分类在留一条件交叉验证下达到83.3%的划分恢复一致性（95% CI 66.4–92.7%），在留一菌株交叉验证下达到66.7%的泛化准确率（随机水平33.3%）。活菌与热处理的二元制备类型分类仅达到63.3%的LOO准确率（95% CI 50.7–74.4%），略高于50%随机水平。
四、讨论与结论
研究人员指出，活菌与热处理之间TNF-α的310.8倍差异不仅是效力大小的体现，更要求机制性解释。活菌制剂的更强免疫刺激可能源于四种互补机制：代谢活跃细菌通过分裂、分泌和自溶持续释放微生物相关分子模式（MAMP）；活菌产生短链脂肪酸、乳酸、过氧化氢和细菌素等免疫活性代谢产物；活力依赖性结构（如被TLR8识别的细菌mRNA）优先激活特定模式识别受体（PRR）通路；活菌代谢导致的培养基酸化独立于受体介导信号调节免疫细胞功能。然而，这些使活菌更具免疫刺激性的机制同时也使其作为实验变量不可控——24小时共培养期间活菌的增殖使实验结束时的细菌与PBMC比例未知且不可知，保守估计初始接种量可能扩增16至256倍。这一解释与观察到的热处理制剂剂量-反应相关性更强的现象一致，但需要终点细菌计数实验直接验证。
关于剂量问题，研究人员在6/10菌株中证明了大多数细胞因子面板的剂量依赖性诱导，但单调剂量-反应缩放假设不能推广至所有菌株。菌株D和E的钟形动力学可能与IL-10介导的自分泌反馈（通过STAT3/SOCS3抑制NF-κB）、TLR脱敏或负性先天信号调节因子（IRAK-M、A20/TNFAIP3、TOLLIP）诱导有关，但单时间点设计无法区分这些机制。
研究的主要局限性包括：单供体设计无法估计个体间变异，需要至少5–8个生物学供体进行多供体验证；活菌增殖的不受控制性；缺乏全基因组序列；部分分析物存在大量超出定量范围的值需要插补；三点浓度范围无法进行正式药理学曲线拟合（EC50、Hill系数）；24小时单时间点无法捕捉反应的时间动态；未包含阳性对照刺激物。
结论方面，该研究建立了整合流式细胞术标准化、多重免疫表型分析和机器学习的循证方法学框架，为益生菌研究中的剂量表征、后生元功能等效性评估和数据驱动的菌株分类提供了新范式。研究直接否证了剂量递增均匀增加免疫活性的假设，引入FED70作为后生元剂量选择的操作性筛选指标，并证明菌株身份信息主要存在于免疫极化方向而非激活幅度。所有结论均限于单供体实验背景，需多供体验证后方可推广至群体水平。
