松嫩平原芦苇湿地退化与修复过程中土壤细菌和甲烷代谢微生物的群落结构

您所在的位置：网站首页 › 芦苇的结构 › 松嫩平原芦苇湿地退化与修复过程中土壤细菌和甲烷代谢微生物的群落结构

松嫩平原芦苇湿地退化与修复过程中土壤细菌和甲烷代谢微生物的群落结构

2024-02-01 12:42| 来源: 网络整理| 查看: 265

甲烷(CH4)是大气中第二大重要的温室气体, 其全球变暖潜力是二氧化碳(CO2)的34倍, 自前工业时代以来, 全球变暖大约有20%是由CH4排放造成的[1].湿地是全球CH4最大的自然来源[2], 大气中CH4的浓度取决于土壤中CH4生成和氧化之间的平衡, 产甲烷菌和甲烷氧化菌是驱动CH4生成和氧化过程的关键微生物, 在湿地甲烷循环过程中发挥重要作用[3~5].产甲烷菌是严格厌氧的古细菌, 目前发现的产甲烷菌全部属于广古菌门; 甲烷氧化菌包括甲烷好氧氧化菌和甲烷厌氧氧化菌, 甲烷好氧氧化菌主要由变形菌门和疣微菌门组成, 甲烷厌氧氧化菌则主要由甲烷厌氧氧化古菌(anaerobic methanotrophic archaea, ANME)、NC10门细菌和未培养细菌(marine benthic group D, MBGD)组成[6].有研究表明, 微生物群落对外界环境因子的变化非常敏感, 湿地产甲烷菌和甲烷氧化菌群落结构易受温度、pH、有机质、水位和植被等多种环境因子的影响[5, 6].

松嫩平原地处半干旱半湿润气候区[7], 是中国内陆盐碱湿地的集中分布区之一, 其盐碱湿地面积约为160×104 hm2[8].芦苇湿地是松嫩平原湿地的主要类型, 然而, 长期气候干旱及人类活动导致芦苇湿地出现了不同程度的退化[9~12].湿地植被退化会影响土壤微生物多样性并使湿地生态系统功能发生改变[13], 同时, 湿地的退化将直接导致控制湿地甲烷产生和排放的环境因素和生物因素的改变[14], 将对湿地CH4排放产生重要影响.近几年相关部门在松嫩平原扎龙、向海、科尔沁和莫莫格等重要湿地实施了局部生态修复和重建项目[10, 15], 目前对于松嫩平原退化湿地及修复湿地的相关研究, 主要集中于植被多样性保护、湿地功能提升等方面, 而针对土壤微生物群落结构、多样性差异及甲烷代谢微生物组成变化等方面的研究则鲜见报道.

本研究基于细菌16S rRNA基因、产甲烷菌mcrA基因和甲烷氧化菌pmoA基因的高通量测序, 解析松嫩平原原始未退化芦苇湿地土壤、退化的和正在修复的芦苇湿地土壤中细菌和甲烷代谢微生物的群落结构特征, 并探究细菌和甲烷代谢微生物多样性与环境因子之间的关系, 旨在为了解我国东北地区典型湿地退化与修复过程中的土壤质量演变和微生物动态变化提供参考.

1 材料与方法 1.1 土壤采集与处理

本研究所采集土壤位于松嫩平原(43°13′~48°40′N, 121°30′~127°00′E), 采样地地上植被均为芦苇, 根据湿地不同状态, 选择3种湿地土壤作为研究对象: ①BL(47°11′55.58″N, 124°14′8.20″E), 原始未退化湿地土壤, 位于扎龙自然保护区, 芦苇为优势种, 以原生型沼泽土和泥炭土为主; ②HL(47°22′33.68″N, 123°21′2.47″E), 退化的芦苇湿地土壤, 位于哈拉海湿地保护区外缘, 因农田灌溉排水导致湿地土壤干旱, 中生植物蒿子入侵芦苇湿地; ③ML(45°54′0.80″N, 123°36′38.06″E), 正在实施修复的芦苇湿地土壤, 位于莫莫格自然保护区, 曾因农耕导致土壤退化和旱生植物入侵.采样时间为2018年8月, 每种湿地土壤选取3个采样点作为生物学重复, 每个采样点为5 m×5 m的样方.按照“L”字形采样路线多点采集各样点表层(0~15 cm)土壤并混合均匀作为一份样本, 共计9份样本.将新鲜土壤样本在4℃保存下运送至实验室, 称取大约5 g土样置于-80℃冰箱中待提取DNA及后续实验, 剩余的土壤储存于-20℃冰箱中用于土壤理化性质测定.

1.2 土壤基本理化性质测定

土壤pH根据电位法利用pH计测定, 总碳(total carbon, TC)和总氮(total nitrogen, TN)利用元素分析仪进行测定, 土壤含水率(soil water content, SWC)通过烘干法进行测定.

1.3 土壤DNA提取与高通量测序

土壤总DNA的提取使用FastDNA SPIN试剂盒(MP Biomedicals, Santa Ana, California, USA), 按照说明书进行操作.DNA的提取质量通过1%的琼脂糖凝胶电泳进行检测, DNA的浓度及纯度使用NanoDrop 2000(Thermo scientific, Wilmington, DE, USA)进行测定.

使用通用引物338F(5′-ACTCCTACGGGAGGC AGCA-3′)和806R(5′-GGACTACHVGGGTWTCTA AT-3′)扩增原核生物16S rRNA基因的V3-V4高可变区[16]; 使用引物MLfF(5′-GGTGGTGTMGGAT TCACACARTAYGCWACAGC-3′)和MLrR(5′-TTCATTGCRTAGTTWGGRTAGTT-3′)扩增产甲烷菌mcrA基因[17]; 使用引物A189F(5′-GGNGACTGGGACTTCTGG-3′)和mb661R(5′-CCGGMGCAACGTCYTTACC-3′)扩增甲烷氧化菌pmoA基因[18]. PCR反应体系为: 5×TransStart FastPfu缓冲液4 μL, 2.5mmol·L-1 dNTPs 2 μL, 上游引物(5 μmol·L-1)0.8 μL, 下游引物(5μmol·L-1)0.8 μL, TransStart FastPfu DNA聚合酶0.4 μL, 模板DNA 10 ng, ddH2O补足至20 μL.每个样本3个重复.PCR扩增程序依次为: 95℃预变性3 min, 95℃变性30 s, 52℃退火30 s(mcrA基因和pmoA基因退火温度为60℃), 72℃延伸30 s, 循环30次, 之后72℃稳定延伸10 min.将同一样本的PCR产物混合后使用2%琼脂糖凝胶回收PCR产物, 利用AxyPrep DNA Gel Extraction Kit(Axygen Biosciences, Union City, CA, USA)进行回收产物纯化, 2%琼脂糖凝胶电泳检测, 并用QuantusTM Fluorometer(Promega, USA)对回收产物进行检测定量.使用NEXTFLEX® Rapid DNA-Seq Kit进行建库, 利用Illumina公司的MiSeq PE300平台进行测序(上海美吉生物医药科技有限公司).

1.4 测序结果处理

使用Fastp软件(https://github.com/OpenGene/fastp, version 0.20.0)[19]对原始测序序列进行质控, 使用FLASH软件(http://www.cbcb.umd.edu/software/flash, version 1.2.7)[20]进行拼接, 使用UPARSE软件(http://drive5.com/uparse/, version 7.1)[21]根据97%的相似度对序列进行OTU聚类并剔除嵌合体[21, 22].利用RDP classifier(http://rdp.cme.msu.edu/, version 2.2)[23]对每条序列进行物种分类注释, 16S rRNA基因序列比对SILVA数据库(https://www.arb-silva.de/, version 132)[23], mcrA和pmoA基因序列比对FunGene数据库(http://www.fungene-db.fr/, version 9.6)[24], 设置比对阈值为70%.

1.5 统计分析

本研究使用Galaxy数据分析平台(http://mem.rcees.ac.cn:8080/)[25]和Origin软件(version 2018)对测序数据进行分析处理, 分析项目包括α多样性分析、环境因子与α多样性的关联分析和群落结构组间差异分析.

采用Shannon指数和Simpson指数表征微生物α多样性[26], Shannon指数与微生物α多样性呈正相关, Simpson指数与微生物α多样性呈负相关[27].采用单因素方差分析(ANOVA)检验样品各项理化因子以及α多样性指数的组间差异是否显著.对土壤理化因子和α多样性指数进行Pearson相关性分析, 以确定土壤理化性质与土壤细菌、产甲烷菌和甲烷氧化菌α多样性之间的相关性[28].基于OTU水平对微生物群落进行PCoA分析, 利用ANOSIM检验方法判断群落组间差异是否显著[29].在不同分类学水平上统计各样本的物种丰度, 采用ANOVA方法检验优势菌群组间差异是否显著, 并采用LEfSe(linear discriminant analysis effect size)分析方法寻找细菌群落组间具有显著性差异的物种[30].

2 结果与讨论 2.1 土壤理化性质

3组湿地土壤样品的基本理化性质如表 1所示.各样本土壤pH均为碱性, 其中, 修复中的芦苇湿地土壤pH最高, 其次是原始湿地土壤, 退化湿地的土壤pH最低.原始湿地土壤总碳、总氮指标显著高于退化湿地, 这可能与退化湿地地表植物生物量降低导致土壤有机质含量降低有关.原始湿地芦苇生长于淹水沼泽, 故含水率最高, 平均含水率为67.41%; 退化的湿地由于植被退化导致湿地蓄水能力弱, 因此含水率下降, 平均含水率为44.54%; 同样, 修复湿地的芦苇土壤样品由于长期植被覆盖退化造成表土含水率显著低于原始芦苇湿地土壤样品, 平均值为24.3%.

表 1 (Table 1) 表 1 土壤样品基本理化性质1) Table 1 Basic physical and chemical properties of the soil samples 样品编号 pH TC/g·kg-1 TN/g·kg-1 SWC/% BL 8.77±0.11a 27.57±2.66a 1.52±0.17a 67.41±4.89a HL 8.38±0.03b 15.82±2.06b 0.81±0.25b 44.54±0.95b ML 8.85±0.08a 19.7±0.78b 1.2±0.07ab 24.3±3.02c 1)数据为平均值±标准偏差, n=3, 不同小写字母表示差异性显著水平表 1 土壤样品基本理化性质1) Table 1 Basic physical and chemical properties of the soil samples 2.2 土壤细菌和甲烷代谢微生物α多样性

各土壤样品细菌、产甲烷菌和甲烷氧化菌的Shannon指数和Simpson指数如表 2所示.原始湿地土壤的细菌和产甲烷菌群落Shannon指数高于退化湿地和修复湿地土壤样品; 对于甲烷氧化菌群落, 退化土壤的Shannon指数最高, 原始土壤次之, 修复土壤最低.而细菌群落的Simpson指数变化趋势与Shannon指数相反.结果表明湿地退化导致土壤细菌和产甲烷菌的多样性降低, 甲烷氧化菌的多样性增加.Pearson相关性分析表明, 细菌和产甲烷菌的α多样性与土壤含水率呈显著正相关关系(表 3), 这可能是由于土壤含水率越高, 厌氧程度越大, 而产甲烷菌属于严格厌氧菌, 因此含水率越高的土壤产甲烷菌的多样性也越高[2, 31].

表 2 (Table 2) 表 2 土壤样品微生物α多样性1) Table 2 Values of α-diversity in microbial communities of the soil samples 样品编号细菌产甲烷菌甲烷氧化菌 Shannon Simpson×10-3 Shannon Simpson×10-3 Shannon Simpson×10-3 BL 6.68±0.06a 3.37±0.81a 3.8±0.17a 51.29±13.83a 3.06±0.11a 113.54±7.62a HL 6.46±0.04b 3.64±0.4a 3.48±0.17ab 73.18±25.04a 3.2±0.17a 90.24±17.41a ML 5.25±0.39b 20.25±9.52a 3.05±0.22b 85.32±24.32a 2.19±0.84a 298.45±207.35a 1)数据为平均值±标准偏差, n=3, 不同小写字母表示差异性显著水平表 2 土壤样品微生物α多样性1) Table 2 Values of α-diversity in microbial communities of the soil samples 表 3 (Table 3) 表 3 土壤理化性质与微生物α多样性的Pearson相关性1) Table 3 Pearson correlation coefficients for environmental factors and microbial α-diversity 项目细菌产甲烷菌甲烷氧化菌 Shannon Simpson Shannon Simpson Shannon Simpson pH -0.462 0.496 -0.234 -0.00975 -0.457 0.446 TC 0.31 -0.193 0.508 -0.542 0.023 2 -0.016 4 TN 0.083 4 0.001 67 0.313 -0.477 -0.169 0.166 SWC 0.862** -0.727* 0.858** -0.568 0.601 -0.569 1)*和**表示差异性显著水平, *为P < 0.05, **为P < 0.01 表 3 土壤理化性质与微生物α多样性的Pearson相关性1) Table 3 Pearson correlation coefficients for environmental factors and microbial α-diversity 2.3 土壤细菌和甲烷代谢微生物的群落组成

PCoA分析结果(图 1)显示3组样本的微生物群落结构存在显著差异(P < 0.05).细菌群落的PCoA分析结果如图 1(a)所示, PC1和PC2分别解释了66.29%和24.02%的群落结构差异.产甲烷菌群落的PCoA分析结果如图 1(b)所示, PC1和PC2分别解释了51.62%和32.98%的群落结构差异.甲烷好氧氧化菌群落的PCoA分析结果如图 1(c)所示, PC1和PC2分别解释了55.28%和30.33%的群落结构差异.

图 1 Fig. 1 图 1 基于OTU水平的微生物群落的PCoA分析 Fig. 1 Principal coordinates analysis(PCoA)for microbial communities based at the OTU level

3组土壤样本门水平细菌群落组成如图 2所示.各样本相对丰度超过1%的菌门为Proteobacteria(变形菌门)、Actinobacteria(放线菌门)、Chloroflexi(绿弯菌门)、Acidobacteria(酸杆菌门)、Firmicutes(厚壁菌门)、Bacteroidetes(拟杆菌门)、Gemmatimonadetes(芽单胞菌门)、Rokubacteria和Nitrospirae(硝化螺菌门), 其中相对丰度最高的是Proteobacteria, 其次是Actinobacteria, 这两种菌门总相对丰度超过50%, 与已有的我国东北湿地细菌群落结构研究结果相符[32, 33].但退化湿地土壤中Proteobacteria和Actinobacteria的相对丰度低于原始土壤和修复土壤, Acidobacteria的相对丰度高于原始土壤和修复土壤.Proteobacteria在营养元素含量高的土壤中相对丰度更高[34, 35], 因此由于湿地退化导致土壤总碳、总氮水平降低, Proteobacteria的相对丰度也降低.而在修复中的芦苇湿地样品中, 地上植物凋落物的水解及转化会使土壤总碳、总氮水平升高, Proteobacteria的相对丰度也相应升高.而Acidobacteria相对丰度的变化可能与土壤pH有关, 有研究表明土壤中Acidobacteria的相对丰度与pH呈显著负相关[36], 因此在pH相对较低的退化土壤中其相对丰度较高.

图 2 Fig. 2 图 2 细菌门水平群落组成 Fig. 2 Community composition of bacteria at the phylum level

在纲水平至属水平上, 3组样本细菌群落组成差异显著(图 3).原始芦苇湿地土壤的差异优势物种包括4目2科, 退化土壤的差异优势物种包括2纲5目8科4属, 修复土壤的差异优势物种包括3纲8目8科9属.在这些差异优势物种中, Rhizobiales(根瘤菌目)和Bacillus(芽孢杆菌属)属于根际促生菌, 是一类能直接或间接促进植物生长及吸收营养物质的根际有益微生物[37]; 而Burkholderiaceae(伯克氏菌科)中包含众多动植物病原菌如Burkholderia(伯克氏菌属)和Ralstonia(青枯菌属)[38].原始芦苇湿地土壤中含有相对丰度较高的Rhizobiales, 退化湿地土壤中含有相对丰度较高的Burkholderiaceae, 修复湿地土壤中则含有相对丰度较高的Bacillus, 表明湿地退化导致部分根际促生菌相对丰度下降, 致病菌相对丰度上升.此外, 退化湿地土壤中含有相对丰度较高的Sphingomonas(鞘氨醇单胞菌属)和Rubrobacter(红色杆菌属), Sphingomonas可以利用多环芳烃等有机污染物作为碳源和能源[39], Rubrobacter(红色杆菌属)具有耐受高强度辐射的特性, 可以生长在极端高辐射环境中[40].由于人类活动造成的松嫩平原湿地退化以及土壤污染、植被覆盖度下降和土壤贫瘠化[10], 而使这些抗逆性强的微生物的相对丰度增加.

图 3 Fig. 3 (a)差异物种进化分支, (b)差异物种LDA值图 3 基于LEfSe分析的细菌群落差异物种进化分支和LDA值 Fig. 3 Cladogram and LDA scores of different bacterial communities based on LEfSe analysis

各样地产甲烷菌科水平群落组成如图 4所示.3组样本共同的产甲烷菌优势菌科为Methanosarcinaceae(甲烷八叠球菌科)和Methanobacteriaceae(甲烷杆菌科), 在其他关于湿地产甲烷菌群落的研究中, 这两种产甲烷菌科也占据优势地位[41].但3组样本产甲烷菌优势菌科的相对丰度存在差异(图 5).原始湿地土壤中Methanobacteriaceae相对丰度最高, 约为63%, 其次是Methanosarcinaceae, 约为30%.退化的和正在修复的芦苇湿地土壤中Methanosarcinaceae相对丰度最高, 分别为48%和46%, 而Methanobacteriaceae相对丰度差异较大, 分别为24%和12%.此外, 修复土壤中还出现了相对丰度较高的Methanomicrobiaceae(甲烷微菌科), 该菌科在修复土壤中相对丰度为16%, 而在原始土壤和退化土壤中的相对丰度很低, 均为1%左右.在这些优势菌科中, Methanomicrobiaceae属于氢型产甲烷菌, Methanobacteriaceae同时具有氢型和甲基型产甲烷的活性, Methanosarcinaceae则可以利用氢型、乙酸型和甲基型这3种主要途径产甲烷[42].目前的研究发现只有Methanosarcinaceae和Methanosaetaceae(鬃毛甲烷菌科)能够利用乙酸途径产甲烷[6], 表明湿地退化导致乙酸型产甲烷菌的相对丰度上升.根据上述3类产甲烷菌在3个样地中的相对丰度情况, 氢型营养代谢路径为湿地产甲烷菌基本的甲烷产生代谢路径[43].Methanosarcinaceae由于具有较高合成效率和低能耗的乙酸型代谢路径能够在自然环境中占据优势生态位[44], 是松嫩平原芦苇湿地产甲烷菌群落的重要组成部分; Methanobacteriaceae相对丰度的变化则可能归因于其以水作为代谢底物氧化甲基化合物的甲基型甲烷产生代谢路径[45], 原始的、退化的和正在修复的芦苇湿地土壤的含水率依次降低, Methanobacteriaceae的相对丰度也由此依次下降; Methanomicrobiaceae可能因与Methanobacteriaceae存在竞争关系而只在修复的湿地土壤中相对丰度有所上升.

图 4 Fig. 4 图 4 产甲烷菌科水平群落组成 Fig. 4 Community composition of methanogens at the family level 图 5 Fig. 5 图 5 产甲烷菌优势菌科相对丰度 Fig. 5 Relative abundances of dominant methanogens at the family level

各样地甲烷好氧氧化菌属水平群落组成如图 6所示, 优势菌属相对丰度差异如图 7所示.原始湿地土壤中优势菌属为Methylocystis(甲基孢囊菌属)和 Methylosarcina(甲基八叠球菌属), 相对丰度分别约为43%和28%.退化湿地土壤中优势菌属为Methylocystis、Methylobacter(甲基杆菌属)、Methylomonas(甲基单胞菌属)和Methylococcus(甲基球菌属), 相对丰度分别约为27%、27%、24%和12%.修复湿地土壤中优势菌属为Methylocystis(甲基孢囊菌属), 相对丰度约为45%.Methylocystis属于Type Ⅱ型甲烷好氧氧化菌, 普遍存在于湿地等多种生态系统中[22, 46, 47], 该菌属是3组湿地土壤样本的共同优势菌属.Methylocystis能够在非常低的甲烷浓度下生长, 是贫营养型甲烷氧化菌, 某些Methylocystis菌株甚至可以在大气中氧化甲烷长达3个月[48], 并且Methylocystis在周期性排水的生态系统中普遍存在.因此, Type Ⅱ型甲烷好氧氧化菌在松嫩平原芦苇湿地中发挥着重要作用.Methylosarcina属于Type Ⅰ型甲烷好氧氧化菌, 最初分离自垃圾填埋场土壤[49], 在已有的关于松嫩平原湿地甲烷氧化菌的研究中, 该菌属不占据主要地位[50], 然而在本研究中, Methylosarcina在原始芦苇湿地土壤中占优势地位, 在退化湿地土壤和修复湿地土壤中的相对丰度较低.原始湿地土壤中出现高丰度的Methylosarcina可能与样地土壤理化性质及土壤性状有关, 这种甲烷好氧氧化指示菌可能是该湿地区域潜在的发挥碳汇作用的微生物资源.

图 6 Fig. 6 图 6 甲烷好氧氧化菌属水平群落组成 Fig. 6 Community composition of methanotrophs at the genus level 图 7 Fig. 7 图 7 甲烷好氧氧化菌优势菌属相对丰度 Fig. 7 Relative abundances of dominant methanotrophs at the genus level

此外, 在退化湿地样品中占优势地位的Methylobacter、Methylomonas和Methylococcus其相对丰度均显著高于原始湿地和修复湿地.这3种菌属均属于Type Ⅰ型甲烷好氧氧化菌.其中, Methylobacter和Methylomonas可以适应极端低温环境[51], 研究者在西伯利亚永久冻土和高寒湿地中检测到了大量的Methylobacter[52, 53], Methylomonas分类下也已经分离到了能在深层地下水中生存的嗜冷菌株[54].Methylococcus则可以耐受极端高温环境, 有研究者从热泉中发现了该属下的嗜热菌株[55].可见, 芦苇湿地退化会导致耐受极端环境的甲烷氧化菌属的相对丰度增加.

3 结论

(1) 与原始芦苇湿地相比, 松嫩平原退化的芦苇湿地土壤理化性质发生了变化.具体表现为: 退化湿地土壤总碳、总氮和土壤含水率均显著降低.

(2) 芦苇湿地退化导致土壤细菌和产甲烷菌的α多样性降低, 甲烷氧化菌的α多样性升高, 而细菌和产甲烷菌的α多样性与土壤含水率呈显著正相关关系, 含水率越高的湿地土壤产甲烷菌的多样性也越高.

(3) 湿地退化导致根际促生菌Rhizobiales的相对丰度下降, 致病菌Burkholderiaceae、耐污染细菌Sphingomonas、抗辐射细菌Rubrobacter以及Type Ⅰ型耐受极端环境的甲烷好氧氧化菌Methylobacter、Methylomonas和Methylococcus的相对丰度上升.而在原始及修复的芦苇湿地中优势菌属为Type Ⅱ型甲烷好氧氧化菌Methylocystis.因此, 不同的湿地状态可以间接改变土壤性状进而改变湿地甲烷代谢菌群落结构.

致谢: 中国科学院沈阳应用生态研究所及中国科学院东北地理与农业生态研究所在野外考察及采样方面提供帮助, 在此一并致谢！

【本文地址】

松嫩平原芦苇湿地退化与修复过程中土壤细菌和甲烷代谢微生物的群落结构

松嫩平原芦苇湿地退化与修复过程中土壤细菌和甲烷代谢微生物的群落结构

今日新闻

推荐新闻