Page 35 - 《广西植物》2024年第7期
P. 35
7 期 席守鸿等: 南亚热带乡土树种与桉树人工林土壤真菌群落多样性和功能类群的比较 1 2 3 5
表 1 5 个林分的立地条件和林分特征
Table 1 Site conditions and stand characteristics of five stands
海拔 坡度 胸径 树高 林分密度
林分类型 坡向
Altitude Slope gradient DBH Tree height Stand density
Stand type Slope aspect
 ̄2
(m) (°) (cm) (m) (planthm )
PMP 540 SW 31 27.5 19.7 435
MMP 550 NE 35 17.4 18.0 1 225
MLP 550 NE 30 19.2 19.7 1 208
24.7(I) 18.8(I) 340(I)
CHP 530 S 30
10.2(Ⅱ) 8.5(Ⅱ) 612(Ⅱ)
EUGP 475 S 15 5.8 7.5 2 300
注: PMP. 马尾松林ꎻ MMP. 火力楠林ꎻ MLP. 米老排林ꎻ CHP. 红锥林ꎻ EUGP. 尾巨桉林ꎮ I. 主林层ꎻ Ⅱ. 次林层ꎮ 下同ꎮ
Note: PMP. Pinus massoniana plantationꎻ MMP. Michelia macclurei plantationꎻ MLP. Mytilaria laosensis plantationꎻ CHP. Castanopsis
hystrix plantationꎻ EUGP. Eucalyptus urophylla × E. grandis plantation. I. Main forest layerꎻ Ⅱ. Secondary forest layer. The same below.
20 μL 反应体系为 5 × FastPfu Buffer 4 μL、 dNTPs 据库来解析各样品真菌群落的营养型和功能类
 ̄1  ̄1 群ꎮ 为了保证解读真菌功能类群的可靠性ꎬ只保
(2.5 mmolL )2 μL、Forward Primer(5 μmolL )
 ̄1 留置信 度 为 “ 极 可 能” ( highly probable) 和 “ 很 可
0.8 μL、Reverse Primer(5 μmolL )0.8 μL、FastPfu
Polymerase 0. 4 μL、BSA 0. 2 μL、Template DNA 10 能”(probable)2 个等级(Nguyen et al.ꎬ 2016)ꎮ
ngꎬ补ddH O至 20 μLꎮ 95 ℃ 预变性 3 minꎬ之后 95 1.5 数据处理
2
℃ 30 sꎬ55 ℃ 30 sꎬ72 ℃ 45 sꎬ35 个循环ꎬ然后 72 利用 Kruskal ̄Wallis 秩和检验法对不同林分土
℃ 延伸 10 minꎬ10 ℃ 至停止ꎬ达到 PCR 扩增条件ꎮ 壤真菌门和目进行差异显著性检验ꎮ 土壤样品真
每个样本 3 个重复ꎬ将同一样本的 PCR 产物混合后 菌群落的 α 多样性采用基于 OTU 的 Chao1 指数、
用 2%琼脂糖凝胶电泳检测ꎬ使用 AxyPrep DNA 凝 Shannon 指数和 Simpson 指数来表征ꎬ其中 Simpson
胶回收试剂盒( AXYGEN 公司) 切胶回收 PCR 产 指数越小代表多样性越大ꎬ由 Mothur 软件完成计
物ꎬ并将 PCR 产物用 QuantiFluor  ̄ST 蓝色荧光定 算(Patrick et al.ꎬ 2017)ꎮ 不同林分之间土壤理化
TM
量系统(Promega 公司)进行检测定量ꎮ 委托美吉生 性质、真菌群落 α 多样性的差异显著性采用单因
物科 技 ( 上 海) 有 限 公 司 ( https:/ / www. majorbio. 素方差分析( one ̄way ANOVA) 检测ꎬ并用 Duncan
com)用 Illumina MiSeq PE300 平台进行文库构建以 法进行多重比较ꎻ使用 Spearman 秩相关系数检验
及高通量测序ꎮ 全部样品的高通量测序结果已提 α 多样性与土壤理化性质因子之间的相关性ꎮ 以
交至 NCBI SRA(https:/ / www.ncbi.nlm.nih.gov / )ꎬ编 上 计 算 均 使 用 SPSS 26. 0 软 件 ( SPSSꎬ Incꎬ
号为 PRJNA936188ꎮ Chicagoꎬ IL)完成ꎮ
1.4 生物信息学分析 基于 OTU 的 Bray ̄Curtis 距 离ꎬ 利 用 R 软 件
使用 Trimmomatic 软件对测序读数( reader) 的 vegan 程序包中的“metaMDS ( )” 函数对不同土壤
接头(adapter)序列和低质量序列进行过滤ꎮ 利用 样品的真菌群落 β 多样性进行非度量多维标度
USEARCH 软 件 ( vsesion 7. 1ꎬ http: / / drive5. com / (non ̄metric multidimensional scalingꎬNMDS)分析ꎻ进
uparse / )按照 97%相似性对非重复序列(不含单序 而用基于 999 次置换检验的置换多元方差分析
列)进行 OTU 聚类ꎬ在聚类过程中去除嵌合体ꎬ得 ( permutational multivariate analysis of varianceꎬ
到 OTU 的 代 表 序 列ꎮ 在 Qiime1 平 台 采 用 RDP PERMANOVA)法检测各土壤样品真菌群落 β 多样
classifier 贝叶斯算法对 97%相似水平的 OTU 代表 性的差异显著性ꎬ计算由 R 软件 vegan 程序包中的
序 列 比 对 Silva ( Release123ꎬ http: / / www. arb ̄ “adonis ( )”函数完成ꎮ 以加权 Bray ̄Curtis 非相似
silva.de)数据库(置信度阈值为 0.7)获得各土壤样 性矩阵作为响应变量ꎬ以土壤理化性质为预测变
品真菌的分类学信息ꎬ将 OTU 表和物种信息表合 量ꎬ采用广义非相似模型 ( generalized dissimilarity
并后ꎬ利用 Python 3.8 软件比对 FUNGuild v1.1 数 modelꎬ GDM)(Ferrier et al.ꎬ 2007)分析土壤理化性