Page 164 - 《广西植物》2025年第11期
P. 164
2 1 0 8 广 西 植 物 45 卷
土壤细菌基因组 DNAꎬ用 1%琼脂糖凝胶电泳检测 视化均 采 用 R ( Version 3. 3. 1) 实 现ꎬ 并 在 Adobe
DNA 提取质量ꎮ 将提取的总 DNA 作为模板ꎬ用引 Illustrator ( Version 2021 ) 与 Adobe Photoshop
物 338F ( 5′ - ACTCCTACGGGAGGCAGCAG - 3′) 与 (Version 2022)中进行图片处理与组合ꎮ
806R(5′-GGACTACHVGGGTWTCTAAT-3′) 对 16S
rDNA V3-V4 区进行 PCR 扩增(Sun et al.ꎬ 2013)ꎮ 2 结果与分析
PCR 扩 增 体 系 为 ( 20 μL) 模 板 DNA ꎬ 其 中 10 ×
Bufferv 2 μLꎬ2.5 mmolL dNTPs 2 μLꎬ双向引物 2.1 土壤理化性质
 ̄1
 ̄1 单因素方差分析结果(表 1) 显示ꎬ3 个研究区
各(5 μmolL ) 0.8 μLꎬrTaq 聚合酶 0.2 μLꎬBSA
0.2 μLꎬ模板 DNA 10 ngꎬ补 ddH O 至 20 μLꎮ PCR 的根际土壤理化性质存在显著性差异( P<0.05)ꎬ
2
反应程序:95 ℃ 3 minꎻ95 ℃ 30 sꎬ56 ℃ 30 sꎬ72 ℃ ZR 处根际土壤的含水量、电导率、有机 物、速 效
45 sꎬ25 个循环ꎻ72 ℃ 10 minꎬ10 ℃停止试验ꎮ 反应 钾、速效磷、铵态氮均显著高于 PR 和 SHꎮ
产物用 1%琼脂糖凝胶电泳检测(李敬科等ꎬ2020)ꎮ 2.2 根际土壤细菌多样性指数分析
样本送至上海美吉测序公司通过 Illumina HiSeq 平 测序共得到 9 个样本的原始序列 20 280 063
台进行高通量测序ꎮ 测序数据上传至国家微生物 条ꎬ质控后得到优化序列 506 328 条ꎬ平均序列长
科学数据中心 ( National Microbiology Data Centerꎬ 度为 446 bpꎬ对 97%相似水平下的 OTU 进行生物
NMDC)ꎬ编号为 NMDC10019540ꎮ 信息统计分析得到 490 个 OTU 集ꎬ22 门 42 纲 85
1.5 数据分析 目 156 科 238 属 335 种ꎮ 随着测序深度的加大ꎬ观
用 FLASH( Version 1.2.7) 软件进行双端测序 测到的细菌数不再增加ꎬShannon 稀释曲线区域平
序列 拼 接ꎬ Trimmomatic ( Version 0. 33) 软 件 过 滤 缓ꎬ平均覆盖度为 99.87%ꎬ说明本次测序深度已
Raw TagsꎬUCHIME(Version 4.2)软件去除嵌合体ꎬ 经足够ꎬ样品测序结果可用于后续分析ꎮ
Usearch 软件对相似度在 97%水平下的 Tags 进行 通过 α 多样性指数评估 3 个研究区根际土壤
聚类、获得 OTU( operational taxonomic unit)ꎬ并以 细菌 多 样 性 与 丰 富 度 特 征ꎬ 包 括 Shannon 指 数、
SILVA 为参考数据 库 对 特 征 序 列 进 行 分 类 学 注 Simpson 指数、Ace 指数和 Chao1 指数ꎮ 结果(图 1)
释ꎮ 通过单因素方差分析( one ̄way ANOVA) 检验 显示ꎬ3 个研究区的根际土壤细菌 Shannon 指数、
土壤 理 化 性 质 的 组 间 差 异ꎮ 用 Mothur ( Version Simpson 指数、Ace 指数无显著性差异(P>0.05)ꎬSH
1. 30 ) 进 行 微 生 物 群 落 多 样 性 ( Shannon 和 与 ZR 的根际土壤细菌 Chao1 指数存在显著性差异
Simpson)、群落丰富度(Chao1 和 Ace) 的 α 多样性 (P<0.05)且 ZR>SHꎬ说明 ZR 处根际土壤细菌丰富
分析ꎬQiime(Version 2020.2.0) 计算 β 多样性距离 度指数显著多于 SHꎻShannon 指数显示ꎬPR 的根际
矩阵ꎬ并采用 vegan( Version 2.4.3) 包进行非度量 土壤细菌多样性高于 ZR 和 SHꎬ表明 PR 处拥有更
多 维 尺 度 ( non ̄metric multidimentional scalingꎬ 高的物种多样性和均匀度ꎻSimpson 指数 ( 基于 λ
NMDS)分析ꎮ 使用 Kruskal ̄Wallis H test 进行组间 值ꎬ数值越高代表多样性越低)排序结果与 Shannon
差异显著性检验ꎬFdr 多重检验校正ꎮ 利用 QMME 指数相反ꎬ进一步验证了 Shannon 指数结果ꎻAce 指
软件生成不同分类水平的物种丰度表ꎬR( Version 数与 Chao 1 指数结果均显示 ZR>PR>SHꎬ表明 ZR
3.3.1)绘制样品各分类学水平下的群落结构ꎮ 通 处根际土壤细菌物种丰富度最高ꎮ
过 典 范 对 应 分 析 ( canonical correspondence 进一步通过 NMDS 分析对比组间差异ꎬ结果
analysisꎬCCA) 检验菌群与环境因子之间的关系ꎬ (图 2)显示ꎬStress<0.2ꎬ表明结果具有解释意义ꎬ3
采用 Spearman 相关性系数计算环境因子与细菌属 个研究区根际土壤细菌样本存在明显分离ꎬ说明 3
水平 之 间 的 关 系ꎬ 并 用 热 图 进 行 可 视 化ꎮ 使 用 个研究区的根际土壤细菌群落组成存在显著性差
Networkx(Version 1.11)进行共线性网络分析ꎮ 使 异(P<0.05)ꎬR = 0.777 8ꎬ表明组间差异大于组内
用 PICRUSt1 和 PICRUSt2 进行功能预测ꎬ比对数 差异ꎮ
据库为 COG 和 KEGG 数据库ꎬ并采用单因素方差 2.3 根际土壤细菌群落组成分析
分析检验优势功能基因的组间差异( 冼康华等ꎬ 2.3.1 门水平群落组成 对藓状雪灵芝根际土壤
2022ꎻ李继琼等ꎬ2025)ꎮ 以上数据分析和结果可 细菌 优 势 门 分 析 发 现 ( 图 3: A )ꎬ 变 形 菌 门

