土壤微生物多样性分析实验报告

土壤微生物多样性分析实验报告一、实验材料与方法（一）土壤样品采集本次实验共采集5份不同生态环境的土壤样品，具体信息如下：农田土壤（S1）：采集于江苏省南京市江宁区某水稻田，深度为0-20cm，采样时间为2026年4月10日。该区域长期种植水稻，每年施用氮磷钾复合肥约500kg/hm²，土壤类型为潴育型水稻土。森林土壤（S2）：采集于安徽省黄山市黄山风景区常绿阔叶林，深度为0-15cm，采样时间为2026年4月12日。该区域植被覆盖率达95%以上，主要树种为黄山松、青冈栎等，土壤类型为山地黄壤。草原土壤（S3）：采集于内蒙古自治区呼伦贝尔市呼伦贝尔草原，深度为0-10cm，采样时间为2026年4月15日。该区域以羊草、针茅为优势物种，年降水量约350mm，土壤类型为黑钙土。荒漠土壤（S4）：采集于新疆维吾尔自治区吐鲁番市库木塔格沙漠边缘，深度为0-5cm，采样时间为2026年4月18日。该区域年降水量不足50mm，植被稀疏，主要为骆驼刺、梭梭等耐旱植物，土壤类型为灰漠土。湿地土壤（S5）：采集于浙江省杭州市西溪国家湿地公园，深度为0-25cm，采样时间为2026年4月20日。该区域常年积水，植被以芦苇、香蒲等水生植物为主，土壤类型为沼泽土。每个采样点采用五点采样法，将5个采样点的土壤混合均匀后，用无菌自封袋密封，置于冰盒中带回实验室，去除石块、植物残体等杂质后，过2mm筛，分为两份，一份置于4℃冰箱保存用于土壤理化性质分析，另一份置于-80℃冰箱保存用于微生物多样性分析。（二）土壤理化性质测定pH值：采用电位法测定，将土壤与去离子水按1:2.5的比例混合，搅拌均匀后静置30min，使用pH计（雷磁PHS-3C）测定上清液pH值。有机质含量：采用重铬酸钾容量法-外加热法测定，称取0.5g过筛土壤样品，加入5mL0.8mol/L重铬酸钾溶液和5mL浓硫酸，在170-180℃油浴锅中加热5min，冷却后用0.2mol/L硫酸亚铁溶液滴定，根据消耗的硫酸亚铁溶液体积计算有机质含量。全氮含量：采用凯氏定氮法测定，称取1g过筛土壤样品，加入浓硫酸和催化剂（硫酸铜、硫酸钾、硒粉），在420℃下消煮，待溶液澄清后，用凯氏定氮仪（FOSSKjeltec8400）测定全氮含量。有效磷含量：采用钼锑抗比色法测定，称取2.5g过筛土壤样品，加入50mL0.5mol/L碳酸氢钠溶液，振荡30min后过滤，取滤液加入钼锑抗显色剂，在700nm波长下测定吸光度，根据标准曲线计算有效磷含量。速效钾含量：采用乙酸铵浸提-火焰光度法测定，称取5g过筛土壤样品，加入50mL1mol/L乙酸铵溶液，振荡30min后过滤，用火焰光度计（上海仪电分析FP640）测定滤液中钾离子浓度，计算速效钾含量。（三）土壤微生物DNA提取与PCR扩增DNA提取：采用E.Z.N.A.®SoilDNAKit（OmegaBio-tek）试剂盒提取土壤微生物总DNA，具体操作步骤如下：（1）称取0.5g冷冻土壤样品，加入到含有玻璃珠和裂解液的离心管中，涡旋振荡10min，使细胞充分裂解。（2）加入200μL蛋白酶K溶液，在55℃水浴中孵育30min。（3）加入200μL裂解液，涡旋振荡1min，在65℃水浴中孵育10min。（4）加入200μL无水乙醇，涡旋振荡1min，将溶液转移到吸附柱中，12000rpm离心1min，弃去滤液。（5）加入500μL洗涤液1，12000rpm离心1min，弃去滤液；加入500μL洗涤液2，12000rpm离心1min，弃去滤液；重复洗涤一次。（6）将吸附柱置于新的离心管中，加入100μL洗脱液，在65℃水浴中孵育5min，12000rpm离心1min，收集洗脱液，即为提取的土壤微生物总DNA。（7）采用Nanodrop2000分光光度计测定DNA浓度和纯度，要求OD260/OD280在1.8-2.0之间，OD260/OD230在2.0以上，同时采用1%琼脂糖凝胶电泳检测DNA完整性。PCR扩增：选取细菌16SrRNA基因V3-V4区进行扩增，引物序列为338F（5'-ACTCCTACGGGAGGCAGCAG-3'）和806R（5'-GGACTACHVGGGTWTCTAAT-3'）；选取真菌ITS1区进行扩增，引物序列为ITS1F（5'-CTTGGTCATTTAGAGGAAGTAA-3'）和ITS2R（5'-GCTGCGTTCTTCATCGATGC-3'）。PCR反应体系（25μL）：12.5μL2×TaqPCRMasterMix（天根生化），1μL上游引物（10μmol/L），1μL下游引物（10μmol/L），2μLDNA模板（10ng/μL），8.5μL无菌水。PCR反应程序：细菌16SrRNA基因：95℃预变性3min；95℃变性30s，55℃退火30s，72℃延伸45s，共30个循环；72℃终延伸10min。真菌ITS1区：95℃预变性3min；95℃变性30s，52℃退火30s，72℃延伸45s，共30个循环；72℃终延伸10min。PCR产物采用2%琼脂糖凝胶电泳检测，目的条带大小约468bp（细菌）和200-300bp（真菌），采用AxyPrepDNAGelExtractionKit（Axygen）试剂盒回收纯化PCR产物。（四）高通量测序与数据分析高通量测序：将纯化后的PCR产物等量混合，构建测序文库，采用IlluminaMiSeqPE300平台进行高通量测序，测序工作由上海美吉生物医药科技有限公司完成。数据分析：采用QIIME2（2022.2版本）软件进行数据分析，具体步骤如下：（1）原始数据质控：使用DADA2插件对原始测序数据进行质控，去除低质量碱基、引物序列和嵌合体，得到高质量的特征序列（ASVs，AmpliconSequenceVariants）。（2）物种注释：使用feature-classifier插件，将ASVs与Silva138数据库（细菌）和UNITE8.3数据库（真菌）进行比对，进行物种注释，注释水平包括界、门、纲、目、科、属、种。（3）α多样性分析：计算Chao1指数、Ace指数、Shannon指数、Simpson指数和Coverage指数，评估土壤微生物群落的丰富度和多样性。（4）β多样性分析：基于Bray-Curtis距离矩阵，进行主坐标分析（PCoA）和非度量多维尺度分析（NMDS），比较不同土壤样品微生物群落结构的差异；采用Adonis分析（置换次数为999）检验不同生态环境对土壤微生物群落结构的影响。（5）物种组成分析：统计不同分类水平下的物种相对丰度，绘制柱状图、热图和进化树，分析不同土壤样品微生物群落的物种组成特征。（6）环境因子关联分析：采用冗余分析（RDA）或典范对应分析（CCA），分析土壤理化性质与微生物群落结构的相关性。二、实验结果与分析（一）土壤理化性质分析不同生态环境土壤样品的理化性质测定结果如表1所示。由表1可知，5份土壤样品的理化性质存在显著差异（P<0.05）。农田土壤（S1）的pH值为6.82，呈微酸性，有机质含量为28.65g/kg，全氮含量为1.52g/kg，有效磷含量为35.26mg/kg，速效钾含量为128.53mg/kg，养分含量较高，这与长期施肥和种植作物有关。森林土壤（S2）的pH值为5.47，呈酸性，有机质含量最高，达45.32g/kg，全氮含量为2.18g/kg，有效磷含量为18.65mg/kg，速效钾含量为95.36mg/kg，这主要是因为森林植被茂密，枯枝落叶多，有机质积累丰富。草原土壤（S3）的pH值为7.85，呈微碱性，有机质含量为22.58g/kg，全氮含量为1.26g/kg，有效磷含量为12.35mg/kg，速效钾含量为156.89mg/kg，速效钾含量较高，可能与草原植被对钾的吸收利用较少有关。荒漠土壤（S4）的pH值为8.92，呈碱性，有机质含量最低，仅为3.25g/kg，全氮含量为0.21g/kg，有效磷含量为2.56mg/kg，速效钾含量为65.32mg/kg，养分极度匮乏，这与荒漠地区降水稀少、植被稀疏有关。湿地土壤（S5）的pH值为6.23，呈酸性，有机质含量为32.56g/kg，全氮含量为1.85g/kg，有效磷含量为25.68mg/kg，速效钾含量为89.56mg/kg，有机质和全氮含量较高，这是因为湿地常年积水，有机质分解缓慢，积累较多。表1不同生态环境土壤样品的理化性质样品编号pH值有机质（g/kg）全氮（g/kg）有效磷（mg/kg）速效钾（mg/kg）S16.82±0.12a28.65±1.23b1.52±0.08b35.26±2.15a128.53±5.62bS25.47±0.08c45.32±2.15a2.18±0.12a18.65±1.02c95.36±4.25cS37.85±0.15b22.58±1.05c1.26±0.06c12.35±0.85d156.89±6.35aS48.92±0.21d3.25±0.25e0.21±0.02e2.56±0.32e65.32±3.15dS56.23±0.10c32.56±1.56b1.85±0.10b25.68±1.56b89.56±4.05c注：同列数据后不同小写字母表示差异显著（P<0.05）。（二）土壤微生物α多样性分析不同生态环境土壤样品的细菌和真菌α多样性指数测定结果如表2所示。由表2可知，5份土壤样品的细菌和真菌α多样性指数均存在显著差异（P<0.05）。在细菌α多样性方面，森林土壤（S2）的Chao1指数、Ace指数和Shannon指数最高，分别为4256.32、4321.56和9.25，Simpson指数最低，为0.008，表明森林土壤细菌群落的丰富度和多样性最高；荒漠土壤（S4）的Chao1指数、Ace指数和Shannon指数最低，分别为856.23、872.15和4.25，Simpson指数最高，为0.256，表明荒漠土壤细菌群落的丰富度和多样性最低；农田土壤（S1）、草原土壤（S3）和湿地土壤（S5）的细菌α多样性指数介于两者之间。Coverage指数均在0.99以上，表明测序深度足够，能够反映土壤微生物群落的真实情况。在真菌α多样性方面，森林土壤（S2）的Chao1指数、Ace指数和Shannon指数也最高，分别为1256.32、1321.56和6.25，Simpson指数最低，为0.025，表明森林土壤真菌群落的丰富度和多样性最高；荒漠土壤（S4）的Chao1指数、Ace指数和Shannon指数最低，分别为156.23、172.15和1.25，Simpson指数最高，为0.656，表明荒漠土壤真菌群落的丰富度和多样性最低；农田土壤（S1）、草原土壤（S3）和湿地土壤（S5）的真菌α多样性指数介于两者之间。表2不同生态环境土壤样品的微生物α多样性指数样品编号细菌Chao1指数细菌Ace指数细菌Shannon指数细菌Simpson指数细菌Coverage指数真菌Chao1指数真菌Ace指数真菌Shannon指数真菌Simpson指数真菌Coverage指数S13256.23±125.65b3321.56±132.15b8.25±0.32b0.012±0.002b0.992±0.001a856.23±56.32b872.15±62.15b4.25±0.25b0.125±0.012b0.985±0.002aS24256.32±156.32a4321.56±162.15a9.25±0.42a0.008±0.001a0.995±0.001a1256.32±86.32a1321.56±92.15a6.25±0.32a0.025±0.005a0.990±0.001aS32856.23±105.65c2921.56±112.15c7.85±0.28c0.015±0.002b0.991±0.001a656.23±46.32c672.15±52.15c3.85±0.22c0.156±0.015b0.982±0.002aS4856.23±56.32e872.15±62.15e4.25±0.25e0.256±0.025d0.990±0.001a156.23±26.32e172.15±32.15e1.25±0.12e0.656±0.056d0.975±0.003aS53556.23±135.65b3621.56±142.15b8.55±0.35b0.010±0.002b0.993±0.001a956.23±66.32b972.15±72.15b4.55±0.28b0.112±0.010b0.988±0.002a注：同列数据后不同小写字母表示差异显著（P<0.05）。（三）土壤微生物β多样性分析基于Bray-Curtis距离矩阵的细菌和真菌群落主坐标分析（PCoA）结果如图1和图2所示。由图1可知，PCoA1和PCoA2的解释率分别为35.26%和22.15%，累计解释率为57.41%。5份土壤样品的细菌群落明显分为5组，其中森林土壤（S2）和湿地土壤（S5）的细菌群落较为接近，农田土壤（S1）和草原土壤（S3）的细菌群落较为接近，荒漠土壤（S4）的细菌群落与其他样品差异显著。Adonis分析结果显示，不同生态环境对土壤细菌群落结构的影响显著（R²=0.856，P=0.001）。由图2可知，PCoA1和PCoA2的解释率分别为42.15%和25.32%，累计解释率为67.47%。5份土壤样品的真菌群落也明显分为5组，其中森林土壤（S2）的真菌群落与其他样品差异最大，荒漠土壤（S4）的真菌群落与其他样品差异也较大，农田土壤（S1）、草原土壤（S3）和湿地土壤（S5）的真菌群落较为接近。Adonis分析结果显示，不同生态环境对土壤真菌群落结构的影响显著（R²=0.892，P=0.001）。（四）土壤微生物物种组成分析细菌物种组成分析在门水平上，不同生态环境土壤样品的细菌物种组成如图3所示。由图3可知，5份土壤样品的细菌群落主要由变形菌门（Proteobacteria）、酸杆菌门（Acidobacteria）、放线菌门（Actinobacteria）、厚壁菌门（Firmicutes）、绿弯菌门（Chloroflexi）等组成，但相对丰度存在显著差异。农田土壤（S1）中变形菌门的相对丰度最高，达35.26%，其次是放线菌门（22.15%）和厚壁菌门（15.68%）；森林土壤（S2）中酸杆菌门的相对丰度最高，达32.56%，其次是变形菌门（25.32%）和放线菌门（18.65%）；草原土壤（S3）中变形菌门的相对丰度最高，达38.56%，其次是放线菌门（20.15%）和厚壁菌门（16.89%）；荒漠土壤（S4）中变形菌门的相对丰度最高，达65.26%，其次是放线菌门（15.32%）和厚壁菌门（8.65%）；湿地土壤（S5）中变形菌门的相对丰度最高，达32.15%，其次是酸杆菌门（28.65%）和绿弯菌门（12.35%）。在属水平上，不同生态环境土壤样品的细菌物种组成存在更大的差异。农田土壤（S1）中优势属为假单胞菌属（Pseudomonas）、芽孢杆菌属（Bacillus）和节杆菌属（Arthrobacter）；森林土壤（S2）中优势属为酸杆菌属（Acidobacterium）、鞘氨醇单胞菌属（Sphingomonas）和红游动菌属（Rhodoplanes）；草原土壤（S3）中优势属为假单胞菌属（Pseudomonas）、芽孢杆菌属（Bacillus）和根瘤菌属（Rhizobium）；荒漠土壤（S4）中优势属为假单胞菌属（Pseudomonas）、盐单胞菌属（Halomonas）和节杆菌属（Arthrobacter）；湿地土壤（S5）中优势属为酸杆菌属（Acidobacterium）、鞘氨醇单胞菌属（Sphingomonas）和浮霉菌属（Planctomyces）。真菌物种组成分析在门水平上，不同生态环境土壤样品的真菌物种组成如图4所示。由图4可知，5份土壤样品的真菌群落主要由子囊菌门（Ascomycota）、担子菌门（Basidiomycota）、接合菌门（Zygomycota）和壶菌门（Chytridiomycota）等组成，但相对丰度存在显著差异。农田土壤（S1）中子囊菌门的相对丰度最高，达65.26%，其次是担子菌门（22.15%）和接合菌门（8.65%）；森林土壤（S2）中子囊菌门的相对丰度最高，达58.56%，其次是担子菌门（28.32%）和接合菌门（6.89%）；草原土壤（S3）中子囊菌门的相对丰度最高，达72.15%，其次是担子菌门（18.65%）和接合菌门（5.32%）；荒漠土壤（S4）中子囊菌门的相对丰度最高，达85.26%，其次是接合菌门（8.32%）和担子菌门（4.65%）；湿地土壤（S5）中子囊菌门的相对丰度最高，达62.15%，其次是担子菌门（25.65%）和壶菌门（8.32%）。在属水平上，不同生态环境土壤样品的真菌物种组成也存在较大的差异。农田土壤（S1）中优势属为镰刀菌属（Fusarium）、青霉属（Penicillium）和曲霉属（Aspergillus）；森林土壤（S2）中优势属为丝膜菌属（Cortinarius）、鹅膏属（Amanita）和牛肝菌属（Boletus）；草原土壤（S3）中优势属为镰刀菌属（Fusarium）、青霉属（Penicillium）和曲霉属（Aspergillus）；荒漠土壤（S4）中优势属为曲霉属（Aspergillus）、青霉属（Penicillium）和镰刀菌属（Fusarium）；湿地土壤（S5）中优势属为腐霉属（Pythium）、镰刀菌属（Fusarium）和青霉属（Penicillium）。（四）土壤微生物与环境因子的关联分析基于Bray-Curtis距离矩阵的细菌和真菌群落冗余分析（RDA）结果如图5和图6所示。由图5可知，RDA1和RDA2的解释率分别为28.65%和18.32%，累计解释率为46.97%。土壤pH值、有机质含量、全氮含量、有效磷含量和速效钾含量与细菌群落结构均存在显著相关性（P<0.05）。其中，pH值与变形菌门和放线菌门呈正相关，与酸杆菌门呈负相关；有机质含量、全氮含量、有效磷含量和速效钾含量与酸杆菌门、绿弯菌门呈正相关，与变形菌门和放线菌门呈负相关。由图6可知，RDA1和RDA2的解释率分别为32.15%和20.65%，累计解释率为52.80%。土壤pH值、有机质含量、全氮含量、有效磷含量和速效钾含量与真菌群落结构也均存在显著相关性（P<0.05）。其中，pH值与子囊菌门呈正相关，与担子菌门呈负相关；有机质含量、全氮含量、有效磷含量和速效钾含量与担子菌门、壶菌门呈正相关，与子囊菌门呈负相关。三、讨论（一）生态环境对土壤微生物多样性的影响本实验结果表明，不同生态环境土壤样品的微生物多样性存在显著差异，森林土壤的微生物多样性最高，荒漠土壤的微生物多样性最低，这与以往的研究结果一致。森林植被茂密，枯枝落叶多，有机质积累丰富，为微生物提供了充足的碳源和氮源，同时森林环境湿润、温度适宜，有利于微生物的生长繁殖，因此微生物多样性较高；荒漠地区降水稀少，气候干旱，土壤养分匮乏，环境条件恶劣，只有少数耐旱、耐贫瘠的微生物能够生存，因此微生物多样性较低。农田土壤和草原土壤由于受到人类活动或放牧的影响，微生物多样性介于森林