3337 个 OTU × 999 次置换：mBio 高分沉积物群落构建全 R 管线拆解（附核心图）
请添加图片描述
核心图一览（Fig 1--4）
下图从上到下依次对应分析管线的四个关键环节：环境因子与 VFA 响应（Fig 1）→ 群落 α/β 多样性，含 NMDS（Fig 2）→ 共现网络与保留模块（Fig 3）→ 群落构建过程定量（Fig 4）。
导读：Fig 2 的面板 D 就是 NMDS 排序图（binary Jaccard 距离），而 adonis / ANOSIM / MRPP 三种组间差异检验的 P 值收在补充表 Table S1/S2，不与主图同框------这正是很多教程拼接不齐的地方。
分析管线拆解
一、数据来源与样本设计
- 研究系统：缺氧海洋沉积物富集微宇宙，用钼酸盐（molybdate）抑制硫酸盐还原菌（SRB）
- 分组：对照组（Control, C）× SRB 抑制组（Molybdate, M） ，跨越 5 个孵育时间点 D00 / D05 / D12 / D21 / D30
- 规模：每个处理 × 时间点 20 个微宇宙 ，合计 207 个 16S 样本
- 测序：16S rRNA V3--V4（338F/806R），Illumina MiSeq PE300
- 样本标签约定：脚本用 4 字符前缀定位组别与时间（如
T00C、T05M...T30M），后续所有分析都靠它切子集
二、数据处理与特征计算
- 上游聚类：UPARSE v7.0 ，97% 相似度 → 3,337 个 OTU（每样本 341--1,781）
- 物种注释：Silva v138，置信度阈值 0.7
- β 多样性距离：binary Jaccard（OTU 表先二值化，做"存在/缺失"而非丰度）
- 网络构建：仅保留在 ≥ 11 个样本中出现的 OTU 进入 MENs，保证共现稳健性
三、统计模型总览
- 核心分组变量：处理（C/M）× 孵育时间点
- 对比切分方式：脚本固定 12 组对比------"同一天 C vs M"（跨处理）与"同一组内跨天"（跨时间）两类
- 没有混合效应模型：差异检验以置换检验 为主，事后比较用 Tukey HSD
- 不确定性来源：置换 999 次 （adonis / anosim / mrpp）+ iCAMP 100 次 bootstrap 重采样
四、主要分析目标（Aim 1--3）
Aim 1：抑制 SRB 是否改变群落组成与结构
- α 多样性（Observed / Simpson / Shannon）：ANOVA + Tukey HSD
- β 多样性：NMDS 非度量多维尺度排序（binary Jaccard，二维，报告 stress）
- 组间差异三连：adonis（PERMANOVA）+ ANOSIM + MRPP，各 999 次置换
- β 多样性成因拆解：betapart 拆成 nestedness（物种丢失）与 turnover（物种替换）
Aim 2：群落结构与网络稳定性如何变化
- MENs 共现网络 → netbiov / igraph 提取模块与节点角色
- 模块保留性：Fisher 精确检验
- 稳定性模拟：随机移除 / 靶向移除模块枢纽，计算 robustness 与 vulnerability
Aim 3：群落构建过程如何被改写
- iCAMP 框架拆分 HeS（异质选择）/ HoS（同质选择）/ HD / DL（扩散限制）/ DR（漂变）
- 效应量：Cohen's d （分 L/M/S 等级）；组间比较用 Mann-Whitney U
五、结果呈现逻辑
- 抑制 SRB 后，VFA 持续累积、理化因子在 D12--D30 与对照组显著分化
- α 与 β 多样性在对照 / 抑制组之间显著分离，NMDS 上两组沿时间轴各成一簇
- β 多样性差异主要来自 nestedness（物种丢失），而非 turnover（替换）
- 抑制 SRB 使共现网络稳健性下降、脆弱性上升，网络结构从"稳定-复杂度正相关"被改写
- 构建过程层面，同质选择（HoS）与扩散限制（DL）占比被显著改变，随机性增强
- 钼酸盐耐受表型与 Marinilabiliales 类群下降相互印证，形成从环境到群落到表型的闭环
六、一句话总结
16S 原始 reads → UPARSE 3337 OTU → NMDS + adonis/anosim/mrpp 分辨群落分异 → MENs 网络稳定性 → iCAMP 定量构建过程。
写在最后
框架已经全部摊开，但每一步的参数细节、避坑点、结果核对清单没法在一篇短文里穷尽。比如"距离到底用 binary Jaccard 还是 Bray-Curtis""999 次置换怎么设""多个对比要不要补多重校正"，这些脚本里的取舍才是真正卡住复现的地方。
如果你正在做：
- 16S / 宏基因组数据，手头有 OTU 或 ASV 表，却找不到一条从多样性到群落构建的完整管线
- 毕业论文里要复现 NMDS + PERMANOVA/ANOSIM/MRPP 这套"β 多样性三件套"并迁移到自己的样本
- 想系统学 iCAMP 群落构建过程定量，但被中间结果和重采样逻辑绕晕
👉 微生物学分析 ，都在淘宝店「万汇泛科技」（拍前可先私信你的数据类型，帮你判断能不能套这套管线）。