土壤中细菌多样性的高通量测序检测实验报告

土壤中细菌多样性的高通量测序检测实验报告一、实验材料与方法1.1土壤样品采集本实验选取3种不同利用类型的土壤作为研究对象，分别为东北黑土区玉米连作农田土壤（采样点：黑龙江省哈尔滨市呼兰区，东经126°11′、北纬45°58′，采样深度0~20cm）、南方红壤区林地土壤（采样点：福建省三明市沙县，东经117°47′、北纬26°24′，采样深度0~20cm）、西北荒漠区草原土壤（采样点：内蒙古自治区阿拉善左旗，东经105°43′、北纬38°38′，采样深度0~20cm）。每个采样点设置5个重复样方，样方间距10m，去除表面凋落物和石块后，采用五点混合法采集土壤样品，装入无菌离心管后置于4℃冰盒中运输至实验室，过2mm筛后一部分用于土壤理化性质测定，另一部分保存于-80℃冰箱用于DNA提取。1.2土壤理化性质测定土壤pH采用电位法（土水比1:2.5）测定；有机碳采用重铬酸钾容量法-外加热法测定；全氮采用凯氏定氮法测定；速效磷采用碳酸氢钠浸提-钼锑抗比色法测定；速效钾采用醋酸铵浸提-火焰光度法测定；土壤含水量采用烘干法（105℃，12h）测定。3种土壤的基本理化性质如下：东北黑土pH6.82、有机碳28.73g/kg、全氮1.79g/kg、速效磷22.46mg/kg、速效钾158.32mg/kg、含水量22.15%；南方红壤pH5.17、有机碳16.25g/kg、全氮0.98g/kg、速效磷8.72mg/kg、速效钾89.64mg/kg、含水量18.42%；西北荒漠草原土pH8.36、有机碳6.41g/kg、全氮0.32g/kg、速效磷3.15mg/kg、速效钾215.78mg/kg、含水量7.69%。1.3土壤微生物总DNA提取采用MoBioPowerSoil®DNAIsolationKit（MOBIOLaboratories,Carlsbad,CA,USA）提取0.5g土壤样品中的总DNA，提取步骤严格按照试剂盒说明书进行。提取后的DNA用1%琼脂糖凝胶电泳检测完整性，并用NanoDrop2000分光光度计（ThermoFisherScientific,Wilmington,DE,USA）测定浓度和纯度，要求DNA浓度≥50ng/μL，OD260/OD280值在1.8~2.0之间，OD260/OD230值在1.7~2.2之间，合格的DNA样品保存于-20℃备用。1.416SrRNA基因PCR扩增与文库构建选择细菌16SrRNA基因V3-V4可变区作为扩增靶点，引物序列为341F（5'-CCTAYGGGRBGCASCAG-3'）和806R（5'-GGACTACNNGGGTATCTAAT-3'），引物5'端添加Barcode序列用于区分不同样品。PCR反应体系为50μL：2×TaqPCRMasterMix25μL，上下游引物（10μmol/L）各2μL，DNA模板（20ng/μL）3μL，无酶水18μL。PCR反应程序：95℃预变性5min；95℃变性30s，55℃退火30s，72℃延伸45s，共28个循环；72℃终延伸10min。PCR产物用2%琼脂糖凝胶电泳检测，切胶回收目的片段（约460bp），采用AxyPrepDNAGelExtractionKit（AxygenBiosciences,UnionCity,CA,USA）纯化。纯化后的PCR产物用QuantiFluor™-ST蓝色荧光定量系统（Promega,Madison,WI,USA）进行定量，按照每个样品的测序量要求，将不同样品按相应比例混合，构建测序文库，文库合格后采用IlluminaNovaSeq6000平台进行双端测序（PE250）。1.5测序数据处理与生物信息学分析首先对下机的RawReads进行质量控制，使用Trimmomatic软件（v0.39）去除低质量碱基（质量值<20）、引物序列和Barcode序列，利用FLASH软件（v1.2.11）对双端reads进行拼接，拼接条件为最小重叠长度10bp，错配率≤0.2。拼接后的CleanTags采用UCHIME软件（v8.1）去除嵌合体序列，得到有效序列（EffectiveTags）。使用QIIME2软件（v2022.2）对有效序列进行操作分类单元（OperationalTaxonomicUnits,OTUs）聚类，以97%的序列相似性为阈值划分OTUs，同时选取每个OTU的代表性序列。采用SILVA数据库（v138.1）对OTU代表性序列进行物种注释，注释阈值为70%，得到每个OTU在界、门、纲、目、科、属、种水平的分类信息。稀释曲线（RarefactionCurve）、Alpha多样性指数（包括Chao1指数、ACE指数、Shannon指数、Simpson指数、Coverage指数）使用Mothur软件（v1.46.1）计算。Beta多样性分析基于Bray-Curtis距离矩阵进行主坐标分析（PrincipalCoordinatesAnalysis,PCoA）和非度量多维尺度分析（Non-metricMultidimensionalScaling,NMDS），采用ANOSIM检验不同组间群落结构差异的显著性。物种差异分析采用LEfSe（LinearDiscriminantAnalysisEffectSize）分析，LDA阈值设置为4.0，筛选不同土壤类型中的差异显著物种。功能预测采用PICRUSt2软件（v2.4.1），基于KEGG数据库（KyotoEncyclopediaofGenesandGenomes）预测细菌群落的代谢功能通路。二、实验结果2.1测序数据质量评估本次实验共测定15个土壤样品（3种土壤类型，每种5个重复），共得到1247623条RawReads，经过质量控制和拼接后，共获得1198547条EffectiveTags，每个样品的有效序列数在76234~84519之间，平均长度为428bp。所有样品的Coverage指数均在99.2%以上，稀释曲线趋于平缓，表明本次测序深度足够，能够覆盖土壤中绝大多数细菌物种，测序结果可以用于后续分析。2.2土壤细菌群落Alpha多样性分析Alpha多样性指数结果显示，不同利用类型土壤的细菌多样性存在显著差异（P<0.05）。东北黑土的Chao1指数为4862.37±215.46，ACE指数为4917.52±231.84，Shannon指数为9.72±0.31，Simpson指数为0.994±0.002；南方红壤的Chao1指数为3521.89±187.23，ACE指数为3586.74±196.57，Shannon指数为8.45±0.27，Simpson指数为0.987±0.003；西北荒漠草原土的Chao1指数为2147.63±124.35，ACE指数为2198.41±132.76，Shannon指数为6.89±0.22，Simpson指数为0.962±0.005。整体来看，土壤细菌丰富度和多样性表现为东北黑土>南方红壤>西北荒漠草原土，这与土壤有机碳、全氮等养分含量的变化趋势一致，表明土壤养分水平是影响细菌多样性的关键因素。2.3土壤细菌群落组成分析在门水平上，共检测到32个细菌门，其中相对丰度>1%的优势门共8个，分别为变形菌门（Proteobacteria）、酸杆菌门（Acidobacteria）、放线菌门（Actinobacteria）、绿弯菌门（Chloroflexi）、拟杆菌门（Bacteroidetes）、厚壁菌门（Firmicutes）、疣微菌门（Verrucomicrobia）、芽单胞菌门（Gemmatimonadetes），这8个门的总相对丰度占所有细菌的92.7%~96.4%。不同土壤中优势门的相对丰度存在明显差异：东北黑土中变形菌门相对丰度最高（31.26%±2.15%），其次为酸杆菌门（18.73%±1.42%）、放线菌门（15.24%±1.08%）；南方红壤中酸杆菌门相对丰度最高（27.42%±1.87%），其次为变形菌门（22.58%±1.64%）、绿弯菌门（14.36%±1.12%），这与红壤酸性环境适合酸杆菌门生长有关；西北荒漠草原土中放线菌门相对丰度最高（28.67%±2.03%），其次为变形菌门（24.15%±1.79%）、厚壁菌门（12.43%±0.97%），这与荒漠土壤干旱、高盐的环境下，放线菌和厚壁菌的抗逆性更强有关。在属水平上，共检测到628个细菌属，相对丰度>1%的属共21个。东北黑土中的优势属为鞘氨醇单胞菌属（Sphingomonas，4.82%±0.37%）、假单胞菌属（Pseudomonas，4.27%±0.32%）、根瘤菌属（Rhizobium，3.94%±0.29%），这些属大多具有固氮、解磷功能，与农田土壤较高的养分水平和作物根系分泌物输入有关；南方红壤中的优势属为酸杆菌属（Acidobacterium，6.35%±0.42%）、伯克霍尔德菌属（Burkholderia，3.72%±0.28%）、链霉菌属（Streptomyces，3.18%±0.24%），酸杆菌属能够适应酸性低养分环境，伯克霍尔德菌属具有降解复杂有机质的能力；西北荒漠草原土中的优势属为节杆菌属（Arthrobacter，7.24%±0.46%）、芽孢杆菌属（Bacillus，6.87%±0.41%）、诺卡氏菌属（Nocardia，4.35%±0.33%），这些属多数为革兰氏阳性菌，细胞壁较厚，能够耐受干旱和高渗透压环境。2.4土壤细菌群落Beta多样性分析基于Bray-Curtis距离的PCoA分析结果显示，第一主坐标（PC1）解释了48.27%的群落变异，第二主坐标（PC2）解释了18.63%的总变异，3种土壤类型的细菌群落在PCoA图上明显分开，同一类型土壤的重复样品聚集在一起。ANOSIM检验结果显示，不同土壤类型之间的细菌群落结构差异极显著（R=0.927，P=0.001），表明土壤利用类型是驱动细菌群落结构分异的主要因素。NMDS分析结果同样显示3组样品之间存在明显的分离趋势，应力值（Stress）为0.087，小于0.1，表明排序结果可靠。2.5差异显著物种分析LEfSe分析结果显示，共有47个物种在不同土壤中存在显著差异（LDA>4.0）。东北黑土中差异显著的生物标志物共有16个，主要包括变形菌门的α-变形菌纲、根瘤菌目、鞘氨醇单胞菌科、假单胞菌科，拟杆菌门的黄杆菌纲、鞘脂杆菌目，这些物种大多参与土壤碳氮循环、作物养分供应；南方红壤中差异显著的生物标志物共有14个，主要包括酸杆菌门的酸杆菌纲、酸杆菌目、酸杆菌科，绿弯菌门的厌氧绳菌纲、暖绳菌目，这些物种具有耐酸、降解难降解有机质的特性；西北荒漠草原土中差异显著的生物标志物共有17个，主要包括放线菌门的放线菌纲、棒杆菌目、诺卡氏菌科，厚壁菌门的芽孢杆菌纲、芽孢杆菌目、芽孢杆菌科，这些物种普遍具有较强的抗逆性，能够在干旱贫瘠的环境中生存。2.6细菌群落功能预测基于PICRUSt2的功能预测结果显示，土壤细菌群落的功能主要集中在代谢（Metabolism）、遗传信息处理（GeneticInformationProcessing）、环境信息处理（EnvironmentalInformationProcessing）、细胞过程（CellularProcesses）四大类，其中代谢类功能的相对丰度最高，占总功能的78.3%~81.2%。具体代谢通路中，碳水化合物代谢（CarbohydrateMetabolism）相对丰度为15.27%~16.84%，氨基酸代谢（AminoAcidMetabolism）为12.63%~13.92%，能量代谢（EnergyMetabolism）为7.84%~8.72%，辅助因子和维生素代谢（MetabolismofCofactorsandVitamins）为5.36%~6.15%。不同土壤类型的细菌功能存在显著差异：东北黑土中参与氮代谢（NitrogenMetabolism）、磷代谢（PhosphorusMetabolism）、次生代谢物生物合成（BiosynthesisofSecondaryMetabolites）的功能基因相对丰度显著高于其他两种土壤（P<0.05），这与农田土壤较高的养分循环需求有关；南方红壤中参与氙生物降解（XenobioticsBiodegradation）、萜类和聚酮类化合物代谢（MetabolismofTerpenoidsandPolyketides）的功能基因相对丰度显著较高（P<0.05），这与红壤中存在较多难降解的木质素、单宁等物质有关；西北荒漠草原土中参与渗透调节（OsmoticStressRegulation）、孢子萌发（SporeGermination）、DNA修复（DNARepair）的功能基因相对丰度显著较高（P<0.05），这些功能有助于细菌应对干旱、强紫外线等胁迫环境。三、讨论本研究通过高通量测序技术系统解析了3种不同利用类型土壤的细菌多样性特征，结果显示土壤养分含量、pH和水分条件是驱