引言
慢生根瘤菌（Bradyrhizobium）是一类能与大豆等植物建立共生固氮关系的根瘤菌，可减少化学氮肥使用，因此在农业微生物接种研究中被视为模式生物。除了固氮，部分慢生根瘤菌携带编码N2O还原酶的nosZ基因，有助于降低温室气体N2O排放。然而，微生物接种对土著土壤群落的生态影响尚不明确，且田间试验受降雨、土壤异质性等因素干扰，难以直接评估接种菌的存活与竞争。为此，本研究采用封闭土壤微宇宙系统（无植物）以控制环境条件，并将分析尺度从整体群落缩窄至功能或分类学定义的亚群落。
材料与方法
研究选用三种农田土壤：两种火山灰土（ANT和AT）和一种潜育土（G）。接种菌为两株高N2O还原能力的慢生根瘤菌：Bradyrhizobium ottawaenseSG09与B. diazoefficiensUSDA110的nasS突变体，均携带链霉素/壮观霉素抗性及DsRed2标记以便追踪。每克土壤接种约105–106CFU，在25°C黑暗条件下培养249天，于第0、14、28、57、249天取样进行CFU计数与DNA分析。DNA提取后，分别对16S rRNA基因和nosZ支系I基因进行扩增子测序，采用带有独特分子标识符（UMI）的MAUI-seq技术以降低误差。序列经质量过滤、去重和聚类（1个错配阈值）得到OTU，并进行分类学注释。
为分析群落结构，计算了Bray–Curtis和加权UniFrac距离矩阵，并运用四种降维方法生成排序坐标：主坐标分析（PCoA）、Sammon映射、非度量多维尺度分析（metaMDS）以及均匀流形近似与投影（UMAP）。此外，还采用稳健主成分分析（RPCA）和系统发育RPCA基于rCLR转换的计数数据生成嵌入坐标。所有方法均在多个维度（d = 2, 3, 4, 7, 10, 15, 20, 30, 70）下进行计算，以便在高维空间中进行几何分析。具体分析包括：计算群落轨迹的方向相似性（余弦相似度cosθ）、平滑度与线性度；评估信噪比（SNR）、k最近邻距离及其变异系数（CV）、霍普金斯统计量等聚类相关指标；并通过PERMANOVA分析各因子（土壤、时间、接种处理）对群落变异的解释程度。
结果
慢生根瘤菌在三种土壤中的长期存活
CFU计数显示，接种菌株SG09_1在所有土壤中均能长期存活。接种后初期（至第28天）CFU下降，随后稳定在105–106CFU/g水平。AT和G土壤中的CFU约为ANT的三倍。在ANT土壤中比较四株菌（SG09与USDA110的各两个衍生株）的存活性，发现初期USDA110的CFU更高，但到第249天时各菌株间无显著差异，表明长期存活性相似。
整体群落分析
Alpha多样性（Shannon指数、Chao1丰富度）不受接种影响。PERMANOVA显示，全局模型中土壤类型解释的变异最大（R2≈ 0.72），时间次之，接种处理的影响虽小但显著（R2≈ 0.009–0.020）。土壤特异性分析中，接种效应在G土壤最强，AT中等，ANT中未检测到。
通过二维排序可视化发现，Bray–Curtis距离下PCoA和metaMDS主要按土壤类型分组，但难以分辨更细的变化；Sammon映射改善了土壤内部模式的分辨率；UMAP则进一步增强了处理组间的分离。加权UniFrac距离由于包含系统发育信息，减少了土壤间差异，从而提高了土壤内部的分辨率。其中，Sammon映射最佳地保留了原始UniFrac矩阵中观察到的土壤内离散度相对关系。
高维几何分析揭示的群落动态
在加权UniFrac–Sammon空间中，计算了从第0天到第249天的控制组与接种组轨迹间的余弦相似度（cosθ）。二维空间中，ANT和AT的cosθ值很高（0.98和0.96），表明轨迹方向近乎平行；而G土壤的cosθ较低（0.31），暗示响应发散。但随着维度升高，cosθ值收敛至约0.6–0.8，表明低维空间存在压缩效应。分段cosθ分析揭示了土壤间不同的时间模式：ANT和AT土壤在早期区间（0–14天）cosθ低，晚期区间（57–249天）cosθ高，显示后期方向一致性增强；G土壤则随进程cosθ逐渐降低，表明轨迹方向越来越发散。
轨迹平滑度分析显示，ANT对照样本的平滑度值始终较低（约-0.5），反映轨迹方向频繁变化；而G土壤在低至中维度下表现出相对较高的平滑度。轨迹线性度在中等维度（约d = 3–7）出现峰值，随后在更高维度下降。聚类相关指标（SNR、kNN距离的CV、霍普金斯统计量）同样在4–10维达到峰值，更高维度则导致均质化。这些结果表明，中等维度空间在保留群落空间结构保真度方面达到最佳平衡。
基于RPCA坐标重建距离矩阵的PERMANOVA分析显示，低维度下土壤效应主导，时间和接种处理的贡献很小且常不显著。随着维度增加，时间效应的解释力在约d = 4时急剧上升，接种处理效应在约d = 7时上升并变得显著。土壤特异性分析中，接种效应在不同土壤中显现的维度尺度不同：G土壤最早（d = 7），其次是AT（d = 10）和ANT（d = 15）。
nosZ携带群落分析
与16S结果相反，nosZ群落的PERMANOVA显示接种处理效应在所有土壤中均很强（全局R2≈ 0.18–0.25）。土壤特异性分析中，AT土壤的响应最强。值得注意的是，即使接种菌株对应的OTU在接种249天后仍占nosZ群落的20%–50%，但若将该接种源OTU从数据集中排除，则PERMANOVA中接种处理的解释力（R2）下降约90%，且剩余土著慢生根瘤菌群落结构与对照相比变化有限。这表明检测到的大部分变异源于接种菌株本身的加入，而非对土著群落结构的重大扰动。
亚群落分析：基于16S和nosZ的慢生根瘤菌群落动态
在16S数据集中，所有土壤的土著慢生根瘤菌群落均以与B. elkanii和B. arachidis相关的两个OTU为主。接种导致接种源OTU相对丰度在初期急剧上升，随后随时间逐渐下降，但在第249天在所有接种土壤中仍可检测到。PERMANOVA表明接种显著改变了慢生根瘤菌群落结构，即使在ANT土壤中（整体群落水平未检测到效应）也是如此。
在nosZ数据集中，第0天对照土壤中注释为慢生根瘤菌的nosZOTU相对丰度为5%–29%，远高于16S数据集估计的慢生根瘤菌相对丰度。各土壤的OTU组成差异大，优势OTU不同。接种显著改变了携带nosZ的慢生根瘤菌群落结构。同样，排除接种源OTU后，剩余土著OTU的组成与对照相比基本不变。nosZ分析成功区分了SG09和USDA110菌株：在ANT土壤中，SG09菌株存活率高（约21%），而USDA110较低（约1.9%），且USDA110接种样本的坐标随时间逐渐向未接种对照靠拢。
讨论
本研究采用的微宇宙系统确保了接种菌的密闭环境和灵敏检测，通过抗生素抗性标记实现了在富含土著慢生根瘤菌的土壤中对活菌的特异性追踪。这有助于将回收的细胞与16S和nosZ扩增子测序鉴定的特定OTU直接关联。
在分析方法上，加权UniFrac通过弱化浅层系统发育周转，更有效地凸显了接种菌株的影响。在降维方法中，Sammon映射对本地结构差异更敏感，能更好解析土壤内部变异和接种效应。研究强调了将降维不仅作为可视化工具，更作为高维嵌入框架用于几何分析的价值。中等维度（4–10维）能最佳平衡信号分辨与均质化，在此维度下，轨迹几何形态和聚类结构得到最清晰解析。RPCA分析进一步将排序维度与统计推断联系起来，显示土壤、时间和接种效应在不同维度尺度上被检测到，其中接种效应可检测的维度因土壤而异。
从生态背景看，土壤理化性质（总氮、无机氮、pH、铁含量）无一与接种结果呈单调关系。相反，系统发育相似性和土著种群初始丰度可能调节接种效应的表现。例如，ANT土壤的土著群落与接种菌株系统发育距离最近，且土著慢生根瘤菌相对丰度最高，而接种效应却最弱；G土壤则相反。这些因素的影响仅在特定分析空间和维度尺度下变得明显，凸显了分析维度在定量揭示土壤依赖性接种效应中的重要性。
结论与意义
本研究证实，在受控微宇宙条件下，慢生根瘤菌接种能在三种不同土壤中长期存活（249天），且在不严重扰乱土著群落结构的前提下，增强N2O还原功能潜力。接种菌株在nosZ携带群落中能保持20%–50%的优势度。研究创新性地将降维技术重新诠释为高维嵌入框架，用于后续几何分析，从而揭示了被低维空间扭曲或极高维空间均质化所掩盖的生态信号。该方法通过增强微小但有意义效应的分辨率，为微生物接种研究中的菌株选择、应用策略和风险评估提供了稳健的分析框架。这些发现不仅推进了对微生物接种生态影响的理解，也为在背景变异复杂的田间条件下解析微生物动态提供了新的方法论视角。
