丰度分层解析等面积岛屿边界几何对时间驱动细菌群落演替的影响

《Microorganisms》：Abundance-Stratified Responses to Island Geometry During Time-Driven Microbial Succession

　　

理解群落演替轨迹的确定性及其塑造因素不仅是生态学的核心问题，也直接关系到食品加工、工业材料和生态系统管理中的微生物管理。在此框架下，由栖息地空间结构引起的扩散差异，以及这些差异如何随不同丰度类群而变化，为检验群落演替的空间异质性提供了重要视角。研究人员构建了6个面积相同但周长和边界复杂性不同的人工琼脂岛屿；使用共享微生物来源和相同培养条件，在1至96 h的16个时间点进行破坏性采样，每个岛屿×时间组合设置3个独立重复（n = 288），并通过16S rRNA基因扩增子测序表征细菌群落动态。群落在96 h内持续重组，ASV丰富度总体下降，Simpson指数（Simpson index）和Pielou均匀度（Pielou evenness）在后期总体升高。全群落PERMANOVA（permutational multivariate analysis of variance，置换多元方差分析）显示，Time解释51.27%的组成变异（p < 0.001），而Island仅解释1.45%且不显著（p = 0.0734）；Time × Island交互作用亦不显著。Island周长与3个演替速率指标或6个轨迹指标均无显著单调关系。群落持续保持长尾丰度结构。在747个具有可检测时间信号的ASV中，岛屿轨迹一致性与全局ASV丰度秩显著负相关（ρ = ?0.534，p < 0.001），表明随丰度降低，跨岛屿时间一致性总体减弱。按丰度分层后，Island对优势类群（平均相对丰度≥1%，p = 0.1532）不显著，但对中间类群（≥0.1%且<1%，p = 0.0106）和稀有类群（<0.1%，p < 0.001）显著；Time × Island交互作用仅在稀有层显著（p = 0.0014）。总体而言，时间相关变化构成全群落变异的主轴，而边界长度未产生稳定的全群落演替梯度。高丰度成员在岛屿间表现更强的时间同步性，低丰度类群显示更大的跨岛屿轨迹异质性。这些结果支持共享时间相关过程与高丰度核心成员的同步变化之间更密切对应；低丰度成员中岛屿相关差异是否由边界介导的扩散或资源交换引起，尚需直接的过程测定。

研究背景与问题。群落演替是否在共享环境下沿相对一致方向推进，还是因类群性状和局域过程而分异，是群落组装的核心问题。Clementsian观点强调定向且可重复的演替，Gleasonian观点强调物种特异性响应以及环境过滤、扩散、漂变和种间互作共同塑造的轨迹；对快速生长和频繁更替的微生物群落，可在受控条件下以较短时间尺度直接检验。栖息地斑块与外部环境之间的界面可影响生物移动、资源输入和代谢物交换；在面积相同的斑块中，较长周长或更复杂边界具有更大边界暴露，可能增强微生物和外部物质输入以及内部代谢物扩散与移除，从而影响群落成员的生长、替换和共存。然而，边界几何是否足以改变全群落演替方向或速率，还是仅在特定成员或特定阶段留下可检测信号，仍缺乏严格控面积、高时间分辨率的实验检验。碎片化研究还提示斑块尺度和空间配置共同影响生物多样性响应，边缘效应可能受斑块面积、隔离和周围基质混淆，并表现明显时空异质性。微生物群落通常呈显著丰度不均：少数高丰度类群贡献大部分相对丰度并可能主导演替的宽时间结构；低丰度类群虽贡献小，却构成重要稀有生物圈和潜在物种储库，可能随环境变化快速增减。因此，共享时间轨迹与局域空间响应应按丰度分层解析。该研究利用等面积人工琼脂岛屿模型，在共同微生物源和相同培养条件下，以高时间分辨率破坏性采样和16S rRNA基因测序比较不同空间条件下的群落演替，检验Time与Island的相对重要性、边界周长梯度、长尾丰度结构以及单个ASV跨岛屿时间一致性是否随丰度下降而减弱。其结论是时间相关过程构成全群落变化主框架，而边界几何相关差异若存在，更可能是丰度依赖、局部和非单调的；这对理解微生物群落空间异质性及食品加工表面、工业材料等微生物管理具有意义。

关键技术方法（≤250字）。接种物来自自然发酵paocai汤（含白菜、白萝卜、姜、花椒、辣椒、冰糖、盐和纯净水，混合后陶瓷罐自然发酵7 d，均质过滤）。研究人员构建6个面积相同（6361.725 mm²）而周长不同（282.743–491.018 mm）的人工琼脂岛屿A–F，并计算形状复杂度C = P²/(4πA)。在1–96 h的16个时间点破坏性采样，每个岛屿×时间组合设3个独立微宇宙重复，共288个样本。进行16S rRNA基因V3–V4区扩增子测序，DADA2获得ASV（amplicon sequence variant，扩增子序列变体），样本抽平至42,248 reads。分析采用Bray–Curtis距离、PCoA（principal coordinates analysis，主坐标分析）、PERMANOVA/PERMDISP（permutational analysis of multivariate dispersions，置换多元离散度分析）、精确置换检验、丰度分层、SAD（species abundance distribution，物种丰度分布）拟合、Lin’s CCC（concordance correlation coefficient，一致性相关系数）和Kendall’s W（Kendall’s coefficient of concordance，肯德尔和谐系数）等。

研究结果。3.1 群落组成和α多样性呈现一致的时间演替。通过16S rRNA基因扩增子测序和α多样性分析发现，6个岛屿细菌群落从1至96 h均发生显著组成重组，ASV丰富度总体下降；Simpson指数和Pielou均匀度在后期总体升高，Shannon指数呈非线性变化，早期升高后多数岛屿稳定或略降。结论是最一致的变化为群落成员数减少且剩余成员相对丰度分布更均匀。3.2 时间进程是全群落变化主轴，边界长度未产生稳定演替梯度。通过共享PCoA、PERMANOVA/PERMDISP、距离初始参考和周转率分析发现，Time解释51.27%的全群落组成变异（p < 0.001），Island仅解释1.45%且不显著（p = 0.0734），Time × Island交互作用不显著；时间组内离散度显著而岛屿间离散度不显著，故Time项不能仅解释为时间质心的纯移动，但岛屿质心共同移动及距离初始参考和周转率结果提供补充证据。群落单位时间周转率早期最高，后期快速下降；相对1 h参考质心的Bray–Curtis距离随时间增加。周长与3个演替速率指标和6个轨迹指标的精确置换检验经FDR（false discovery rate，错误发现率）校正后均不显著。结论是在固定面积和96 h观察窗内未检测到随周长一致变化的全群落演替速率或轨迹梯度。3.3 群落维持长尾丰度结构且岛屿相关差异较弱。通过固定相对丰度分箱的组成向量分析发现，多数ASV位于最低丰度箱，≥1%的高丰度箱仅含少数成员。ASV丰富度组成受Time显著影响（R² = 23.45%，p < 0.001），Island效应小但显著（R² = 2.71%，p = 0.0141），Time × Island显著（R² = 25.00%，p = 0.0106）；丰度贡献组成中Time显著（R² = 21.33%，p < 0.001），Island主效应不显著（p = 0.1931），交互作用名义显著（p = 0.0402）。SAD模型比较显示Poisson log-normal模型AIC（Akaike information criterion，赤池信息量准则）低于log-series模型，中位AIC差116.94，中位Pearson相关0.923，且该偏好在288/288个单样本中成立。结论是长尾结构在6个岛屿间高度相似，但成员数在丰度箱间的分配存在弱且时间依赖的岛屿差异，而非随周长简单增强的分异。3.4 不同丰度层表现不同演替动态和跨岛屿时间一致性。共检测15,719个ASV；优势层17个ASV（≥1%）仅占ASV数0.108%却贡献82.89%相对丰度，中间层37个ASV贡献11.53%，稀有层15,665个ASV占99.66%却仅贡献5.58%。丰度分层PERMANOVA显示Time在三层均显著（p < 0.001），R²由优势层58.40%降至中间层42.18%和稀有层17.77%；Island在优势层不显著（p = 0.1532），在中间层（p = 0.0106）和稀有层（p < 0.001）显著但效应量小；Time × Island交互仅稀有层显著（p = 0.0014），且稀有层Island-PERMDISP显著（p = 0.0357）。以岛屿为配对区组的Friedman检验显示，距离初始参考斜率和周转率斜率在三层间均不同（FDR分别为0.00496和0.00570）；稀有层早期瞬时周转更高、后期下降更快，但相对1 h参考的共享方向累积变化更弱。Kendall’s W分析显示全部ASV平均W为0.175、中位0.167，128个ASV经BH校正后显著；高一致性示例ASV7和ASV6的W分别为0.942和0.929，低一致性示例ASV355和ASV430的W约0.067和0.068。对747个具有可检测时间信号的ASV，ITC（island trajectory concordance index，岛屿轨迹一致性指数）与全局ASV丰度秩显著负相关（Spearman ρ = ?0.5338，p = 2.95 × 10^?56），且在48种过滤组合中均保持负相关。结论是高丰度成员跨岛时间同步更强，低丰度成员跨岛轨迹异质性更大；但低丰度成员一致性降低本身不能证明其差异由边界几何直接造成。

讨论总结。共享时间相关过程构成群落演替的主要解释轴：共同微生物源、培养基和培养条件提供相似初始物种库和资源背景，后续底物耗竭、代谢物积累及局部理化条件变化可能与跨岛同步重组相伴；Time应理解为 incubation 过程中共同变化的一组生态过程的综合表征，而非时间本身的直接生态效应。早期快速周转、后期减缓提示阶段依赖的组装动态，支持早期快速重组后缓慢调整，而非96 h内达到稳定终点。边界几何未产生稳定单调全群落梯度：边缘效应依赖实际迁移和物质通量，而微尺度扩散受连通性、浓度梯度和扩散距离限制；本实验中圆形孔位于相同塑料板外缘，孔数和分布不随岛屿周长变化，所有岛屿持续暴露于同一外部汤液并通过孔连通，较高基线交换可能降低周长增加的边际效应；周长与曲率、凹凸性和内部至边界距离共变，不能隔离单一几何属性；中间层和稀有层的弱Island效应提示几何影响更可能是局部、成员特异或阶段依赖，而非全群落对周长的单调响应。丰度与跨岛时间一致性相关：群落持续保留少数高丰度成员和大量稀有成员的长尾结构，Poisson log-normal较log-series获更强支持，但SAD本身不能识别单一组装机制；高相对丰度成员跨岛轨迹更一致，低丰度成员更易受漂变、随机定殖、局部微环境变异、零膨胀和检测限制影响。局限包括仅6个几何水平、周长与形状复杂度共变、单一面积、共同微生物源和96 h观察窗；16S rRNA基因扩增子测序仅提供相对丰度，组成性数据不能直接解释绝对丰度变化；迁移和物质通量未直接测定。未来需正交操纵周长、曲率和连通性，结合示踪测量、局部理化变量、qPCR或流式细胞术绝对定量及更长时间序列直接检验。应用上，表面几何和边界结构可能影响定殖与演替动态，进而影响食品加工表面生物膜形成以及钢和其他材料的生物附着和微生物影响腐蚀。

研究结论翻译：本研究表明，在面积相同但边界几何不同的人工微生境中，细菌群落在96 h内沿共享时间方向发生持续重组，Time解释的全群落组成变异显著多于Island。群落相对于1 h参考状态的累积位移增加，而单位时间周转逐步下降；全群落以及优势、中间和稀有层均未显示随边界周长变化的稳定演替速率或轨迹梯度。同时，丰度结构和按丰度分层的群落表现出微弱但可检测的岛屿相关异质性，高丰度ASV在岛屿间表现更一致的时间变化，低丰度成员显示更大的跨岛屿轨迹异质性。总体而言，共享时间相关过程构成全群落变化的主要解释框架，而边界几何相关差异（若存在）更可能是丰度依赖、局部和非单调的。这些差异是否直接由边界介导的扩散或资源交换引起，需通过绝对定量和直接过程测定进一步检验。该机制框架也为管理与微生物相关的工业系统提供了科学依据，包括食品加工表面的微生物定殖、生物膜形成，以及钢和其他材料的生物污损或微生物影响腐蚀（microbiologically influenced corrosion）。