《Foods》：Rare Biosphere Reveals a Decoupling Between Microbial Abundance and Intrinsic Physiological Potential in Shanxi Aged Vinegar Fermentation
编辑推荐：
原位微生物丰度与实际代谢性能之间的差异是解释宏组学数据中微生物功能的一个关键挑战。在此，研究人员整合了宏基因组（metagenomic）和宏转录组（metatranscriptomic）测序，以研究山西老陈醋（SAV）固态发酵中微生物丰度与基于培养的生理潜力之
原位微生物丰度与实际代谢性能之间的差异是解释宏组学数据中微生物功能的一个关键挑战。在此，研究人员整合了宏基因组（metagenomic）和宏转录组（metatranscriptomic）测序，以研究山西老陈醋（SAV）固态发酵中微生物丰度与基于培养的生理潜力之间的这种解耦现象。Lactobacillus acetotolerans 在基因组（40.89%）和转录组（55.36%）水平上均主导群落，而 Pediococcus acidilactici 仅占 0.11%——这是一个典型的稀有生物圈（rare biosphere）成员。通过桑基（Sankey）分析的源追踪表明，涉及乙酸产生的基因主要归因于 Acetobacter pasteurianus，而涉及乳酸产生的基因主要与 Lactobacillus spp. 相关。然而，L. acetotolerans 表现出有限的酸耐受性和乳酸产生能力，而低丰度的 P. acidilactici AAF1-5 在发酵相关条件下表现出强大的胁迫耐受性和优越的乳酸产生能力——这是微生物丰度与生理性能之间的显著对比。代谢互作网络（metabolic interaction network）分析预测，P. acidilactici 可能受到 L. acetotolerans (Ixy = ?2.737, 资源竞争) 和 A. pasteurianus (Ixy = ?1.887, 酸胁迫) 的共同抑制。为了检验这种低丰度是否源于生态限制而非内在代谢能力，研究人员将来自 A. pasteurianus 的热休克共伴侣蛋白基因 grpE 异源表达于 P. acidilactici AAF1-5 中作为实验工具。重组菌株 P. acidilactici-grpE 在乙酸胁迫下表现出显著增强的存活能力，并且在模拟 SAV 发酵中，与野生型对照相比，乳酸含量增加了 23.63%。这些结果表明，宏组学（meta-omic）丰度不一定能预测生理性能，低丰度可能反映生态限制而非内在功能缺陷。本研究提供了一个生态学框架，用于在复杂发酵微生物组中超越基于序列的丰度推断，将微生物丰度与生理功能联系起来。
**论文解读**
## 研究背景、问题与意义
山西老陈醋（SAV）采用固态发酵（SSF）工艺，其中有机酸是风味的核心决定因素，乳酸和乙酸占总有机酸的85%以上，二者比例直接影响醋液的感官平衡。在发酵过程中，乳酸菌（LAB）和醋酸菌（AAB）被认为是核心功能类群，但传统宏组学（meta-omics）方法（如宏基因组学metagenomics和宏转录组学metatranscriptomics）仅能测量群落总基因丰度和转录本活性，无法直接反映单细胞的代谢通量或实际代谢速率。因此，高丰度物种是否一定具有高生理性能，低丰度微生物是否只是“生态冗余”，成为困扰微生物功能解析的关键问题。研究人员前期在SAV微生物组研究中发现，从醋醅中分离出的80余株LAB中，丰度最高的Lactobacillus acetotolerans在实验室培养条件下生长能力和酸耐受性有限，而丰度极低的Pediococcus acidilactici AAF1-5却表现出优异的热耐受性和乳酸产量。这提示生态丰度与生理性能之间可能存在解耦（decoupling），但缺乏系统验证。为此，本研究整合宏组学测序、代谢互作网络分析及实验验证，旨在定量揭示这种解耦现象，并检验低丰度是否源于生态限制而非内在代谢缺陷。该研究发表在《Foods》期刊上，为复杂发酵微生物组中丰度-功能关系的解读提供了生态学框架。
## 主要关键技术方法
研究人员从山西太原的SAV醋酸发酵（AAF）醋醅中采集样本（第1、3、5、7、9天，每个时间点3个平行缸，五点采样法混合）。主要技术方法包括：①宏基因组和宏转录组测序（Illumina NovaSeq 6000平台），数据基于前期研究（NCBI PRJNA689964和PRJNA602043），此次重新分析以进行基因组-转录组功能源追踪；②基于KEGG数据库的代谢通路注释与桑基图（Sankey diagram）源追踪，比较10个核心物种在12个EC（Enzyme Commission）功能类别上的编码潜力与表达活性；③代谢互作网络分析，采用Sung等（2017）的框架计算代谢影响指数I
xy，量化物种间的净促进/抑制效应；④异源表达（heterologous expression）策略：将A. pasteurianus AC2005的酸耐受基因grpE（编码热休克共伴侣蛋白）克隆至穿梭载体pNZ8148，电转化至P. acidilactici AAF1-5，构建重组菌株P. acidilactici-grpE；⑤模拟固态发酵（SSSF）验证：将重组菌株接种至新鲜醋醅中，监测乳酸和乙酸动态。
## 研究结果
### 3.1 有机酸动态
通过高效液相色谱（HPLC）分析，发现乳酸含量在发酵第3天达峰值（2823.4 mg/100 g醋醅）后逐渐下降至约1800-2000 mg/100 g，乙酸则持续积累至第9天的3255.9 mg/100 g，总有机酸升至5348.2 mg/100 g。这反映了早期乳酸主导、后期乙酸主导的微生物代谢演替。
### 3.2 整合宏组学群落结构
宏基因组鉴定出186个属和644个种，宏转录组鉴定出818个属和1786个转录活性种。Lactobacillus在基因组（52.46%）和转录组（70.67%）中均占主导，Acetobacter次之（10.86%和11.79%）。种水平上，L. acetotolerans丰度最高（基因组40.89%，转录组55.36%），而P. acidilactici仅占0.11%，属于典型稀有生物圈成员。
### 3.3 有机酸合成基因源追踪
桑基图显示，L. acetotolerans的乳酸脱氢酶基因（EC:1.1.1.27和EC:1.1.1.28）编码潜力占群落70-87%，但其转录活性仅占自身总转录本的~0.4%，反映出高丰度低表达。A. pasteurianus的乙酸代谢相关EC类目（如EC:1.2.1.22）转录活性较基因组编码潜力提升26倍，表明其代谢活跃。四种物种（如Lactobacillus crustorum）转录活性为零，可能处于休眠状态。
### 3.4 代谢互作网络揭示P. acidilactici受共同抑制
网络分析包含50个物种，L. acetotolerans和A. pasteurianus构成两个中心节点。L. acetotolerans对P. acidilactici的净抑制作用为I
xy = ?2.737（强抑制，资源竞争），A. pasteurianus对P. acidilactici的抑制为I
xy = ?1.887（中等抑制，酸胁迫）。P. acidilactici受到平均I
xy = ?2.312的共同抑制，说明其低丰度可能源于生态压制。
### 3.5 酸耐受重组菌株构建与表征
从A. pasteurianus中成功扩增并克隆了adhA、adhB、aatA和grpE四个酸耐受基因，其中grpE重组菌株P. acidilactici-grpE生长曲线与野生型最接近，峰值生物量提高15%。在乙酸胁迫下，P. acidilactici-grpE在1.5 g/100 g醋醅乙酸浓度时活菌数达4.2×10
6 CFU/g，乳酸产量1.23 g/100 g；在4 g/100 g高浓度时仍保持66 CFU/g和0.38 g乳酸/100 g，分别比野生型提高约47倍和3.5倍。
### 3.6 模拟发酵验证
模拟固态发酵显示，P. acidilactici-grpE组的乳酸含量在整个发酵过程中持续最高，末期达2.25 g/100 g醋醅，较野生型组（1.82 g/100 g）提高23.63%，较空白对照（1.64 g/100 g）提高37.20%。乳酸/乙酸比值较野生型组提高47.37%，而乙酸含量在后期略有下降，可能由于碳限制下P. acidilactici利用乙酸作为碳源。
## 总结讨论与结论
讨论部分指出，本研究核心发现是生态丰度与生理性能的解耦：L. acetotolerans虽丰度占优，但生理性能反而弱于稀有生物圈成员P. acidilactici。解耦的原因可从三个层面理解：方法学层面，宏组学测量基因/转录本丰度而非代谢通量，L. acetotolerans将大部分转录资源投入糖酵解和ATP合成等维持功能，乳酸脱氢酶基因仅占其转录本的~0.4%；生态系统层面，SSF的多重胁迫（高温、高酸、乙醇）选择具有广谱耐受性的“通才”而非专性代谢者；生物学层面，L. acetotolerans本身代谢范围窄，其高丰度可能源于早期生态位抢占而非持续高代谢输出。本研究不仅描述性指出稀有生物圈的功能重要性，还通过代谢互作网络定量识别出P. acidilactici受L. acetotolerans（资源竞争）和A. pasteurianus（酸胁迫）的共同抑制，并通过异源表达grpE缓解酸胁迫，在模拟发酵中证实了其功能潜力的释放。这证明低丰度不应被解释为低功能潜力，而是生态限制的结果。研究局限包括：模拟发酵系统与工业实际存在差异；grpE介导的酸保护机制需蛋白质组学解析；实际应用前需进行食品安全评估（质粒稳定性、标记物去除、食品级表达系统）。结论部分翻译如下：本研究证明，在SAV微生物群落中，宏组学丰度不应直接解释为功能重要性。生态丰度与生理性能的解耦——最丰富的物种（L. acetotolerans，40.89%）在生理上劣于稀有生物圈成员（P. acidilactici，0.11%）——与生态互作（特别是竞争和酸胁迫）可掩盖低丰度微生物内在生理潜力的观点一致。通过整合群落分析、互作网络和生理表征，本研究证明P. acidilactici具有潜在功能能力，该能力在其原生境中至少部分被酸胁迫所掩盖。作为实验工具异源表达grpE，为缓解特定环境约束（酸胁迫）可使该稀有生物圈成员表达其功能潜力提供了支持。需强调，该实验验证仅针对酸胁迫约束，资源竞争仍待未来实验检验。本研究为复杂发酵微生物组中丰度-功能关系的解释提供了生态学框架。从实践角度，识别受生态约束但功能强大的稀有类群，并结合缓解特定约束（如酸胁迫）的策略，可能为SAV生产中的发酵剂设计或生物强化提供参考，但需进行食品安全评估和工业规模验证。
