环境科学  2024, Vol. 45 Issue (10): 6086-6095   PDF    
湖泊沉积物的氮代谢微生物、功能基因及代谢途径分析:以贵州省红枫湖为例
刘倩1, 陈晓1, 李彦澄1,2, 何雨含1, 李江1,2     
1. 贵州大学资源与环境工程学院, 喀斯特地质资源与环境教育部重点实验室, 贵阳 550025;
2. 贵州喀斯特环境生态系统教育部野外科学观测研究站, 贵阳 550025
摘要: 氮循环是湖泊生态系统中物质循环和能量流动的重要过程, 在湖泊沉积物中, 由微生物驱动的氮循环在维持湖泊生态系统平衡中发挥了重要作用.选择贵州省红枫湖作为研究对象, 采集湖泊沉积物进行理化性质分析, 并利用宏基因组学技术分析沉积物中的氮代谢微生物、功能基因、代谢途径及其相互关系.结果表明, 红枫湖沉积物中的氮代谢微生物主要为细菌, 其相对丰度前3的优势菌属分别为硫杆菌属(Thiobacillus, 16.64%)、红长命菌属(Rubrivivax, 9.43%)和硝化螺旋菌属(Nitrospira, 7.09%);共检测到6条氮代谢途径, 分别为固氮、硝化、反硝化、同化硝酸盐还原、异化硝酸盐还原和完全硝化, 主要以反硝化和异化硝酸盐还原途径为主, 未检测到厌氧氨氧化途径;细菌和古菌参与了6条氮代谢途径, 真核生物参与了3条途径, 分别是异化硝酸盐还原、反硝化和完全硝化;氨氮、硝酸盐氮和总磷是氮代谢功能基因分布的主要环境影响因子, 但在不同环境中存在差异;氮代谢功能基因与微生物之间绝大多数呈正相关关系(95.04%), 其中, 丰度最高且宿主最多的氮代谢基因分别为narGnarZnxrA. 研究结果进一步阐释了贵州喀斯特地区典型湖泊的氮循环过程.
关键词: 红枫湖      宏基因组学      氮循环      功能基因      代谢途径     
Functional Genes and Metabolic Pathways of Nitrogen Metabolism Microorganisms in Lake Sediments: A Case Study of Hongfeng Lake, Guizhou Province
LIU Qian1 , CHEN Xiao1 , LI Yan-cheng1,2 , HE Yu-han1 , LI Jiang1,2     
1. Key Laboratory of Karst Georesources and Environment, Ministry of Education, College of Resources and Environmental Engineering, Guizhou University, Guiyang 550025, China;
2. Guizhou Karst Environmental Ecosystem Observation and Research Station, Ministry of Education, Guiyang 550025, China
Abstract: The nitrogen cycle is of great importance for material circulation and energy flow in lake ecosystems. It is driven by microorganisms in lake sediments and can contribute to balancing lake ecosystems. In this study, physical and chemical properties of the sediments sampled from Hongfeng Lake in Guizhou Province were assayed and analyzed using metagenomics to reveal relevant microorganisms, functional genes, metabolic pathways, and their relationships throughout nitrogen metabolism. The results showed that bacteria were dominant, and the top three relative abundant genera were Thiobacillus (16.64%), Rubrivivax (9.43%), and Nitrospira (7.09%). Only six pathways, including nitrogen fixation, nitrification, denitrification, assimilatory nitrate reduction, dissimilatory nitrate reduction, and complete nitrification, were detected in total, of which denitrification and dissimilatory nitrate reduction were the primary processes, but anaerobic ammonia oxidation was not detected. Bacteria and archaea participated in these six pathways, while eukaryotes only functioned in dissimilatory nitrate reduction, denitrification, and complete nitrification. Ammonia nitrogen, nitrate nitrogen, and total phosphorus, as main environmental factors affecting the distribution of functional genes for nitrogen metabolism, differentiated with each other in their respective real-world conditions. A positive correlation (95.04%) was observed between the functional genes and microorganisms, and narG, narZ, and nxrA possessed the highest abundance and the highest host genes. On this basis, these findings are expected to further elucidate the nitrogen cycle of typical karst lakes in Guizhou Province.
Key words: Hongfeng Lake      metagenomics      nitrogen cycle      functional genes      metabolism pathway     

氮循环是湖泊生态系统的物质循环和能量流动的重要过程, 是微生物和化学过程介导的生物地球化学循环[1], 湖泊沉积物中的大量微生物是其主要驱动力[2].对于自然环境中的氮循环, 国内外学者开展了大量研究, 包括森林土壤、海洋和湖泊等自然环境.如徐婷婷等[3]研究了长白山阔叶红松林土壤的硝化、反硝化过程及其功能微生物, 证明了降雨减少在一定程度上抑制了长白山阔叶红松林土壤净硝化作用, 氮的交互作用促进了土壤净硝化作用.Yu等[4]研究了南海浅湾和深海沉积物的氮循环微生物群落, 发现存在固氮、硝化、反硝化和异化硝酸盐还原为铵和厌氧氨氧化5个过程, 在浅湾和深海测得的丰度较高的功能微生物种类存在较大差异.Li等[5]探究了苏必利尔湖不同深度沉积物的氮循环, 估算了沉积物-水界面的硝酸盐和氨氮通量, 计算了沉积物的硝化和反硝化速率, 并发现大型淡水湖沉积物氮的循环与海洋系统相似.此外, 盐碱湖泊的氮循环也是研究热点, 如Feng等[6]通过探究温暖季节滨海富营养化湖泊沉积物中微生物群落的动态变化及其多种环境因子对氮循环的影响, 证实了海水入侵引起盐度升高, 进而增加了反硝化、异化硝酸盐还原和氨氧化作用相关基因的丰度, 降低了固氮和同化硝酸盐还原相关基因的表达.孙小溪[7]从反硝化、厌氧氨氧化和硝化这3条氮代谢途径定量研究青藏高原不同盐度湖泊沉积物氮素转化、相关功能微生物群落结构及环境因子的影响规律, 证明盐度是氮移除和氮荷载的决定性因素.

活性氮作为农业肥料被过度施用, 会对下游生态环境产生不利影响[8, 9], 淡水湖是活性氮的有效汇, 可将高达90%的活性氮转化为氮气[10].在湖泊生态系统中, 不同形态的氮素营养物质通过各种途径进入水体和沉积物并参与湖泊氮素生物地球化学循环, 使湖泊生态系统中富含氮物质[11].但是, 过量的营养物质(如氮、磷和有机物等)输入也会导致湖泊富营养化、水华和赤潮等问题, 进而影响沉积物中微生物群落特征和功能[12].因此, 研究湖泊沉积物中的氮代谢微生物、功能基因和代谢途径, 对了解湖泊的生物地球化学功能至关重要.

本研究选择贵州省喀斯特地区的典型湖泊(红枫湖)作为研究对象, 采集湖泊沉积物进行理化性质分析, 然后利用宏基因组学技术分析沉积物中的氮代谢微生物、功能基因和代谢途径, 识别主要的氮代谢功能微生物和功能基因, 进一步分析氮代谢功能基因与环境因子和微生物之间的关系, 以期进一步阐释贵州喀斯特地区典型湖泊的氮循环过程, 并为贵州省贵阳市的饮用水源地保护提供理论依据.

1 材料与方法 1.1 样品采集

贵州省贵阳市红枫湖位于106°20′19″~106°27'53″E, 26°24′39″~26°34'24″N之间.红枫湖的控制集流水域面积为1 551 km2, 正常库容为6.01亿m3, 相应水面面积为57.2 km2, 设计正常水位高程为1 240 m, 死水位高程为1 227.5 m, 现为贵阳市的重要集中式饮用水源地.根据红枫湖的水域分区特点, 于2023年2月, 在红枫湖选取5个典型干流水域, 3个典型支流水域进行样品采集, 共设置8个采样点, 如图 1所示.每个样点选取3.0 m×3.0 m作为1个样方, 利用沉积物采样器抓取湖底0~10 cm的表层沉积物, 每个样方内随机采集3次样品, 并进行混合, 然后分为2份, 一份装入50 mL灭菌离心管后立即放入干冰桶(-80℃)中保存, 用于宏基因组学分析, 另一份装入无菌采样袋中, 冷藏带回实验室, 用于常规理化指标分析.

图 1 红枫湖采样点分布示意 Fig. 1 Distribution of sampling sites in Hongfeng Lake

1.2 常规理化指标的测定

取过筛后的沉积物进行理化指标的测定.沉积物pH使用pH计(Pocket Pro+ Multi2)测定, 电导率(EC)采用电极法(HJ 802-2016)测定, 硫酸根(SO42-)采用硫酸钡沉淀法测定, 有机质(SOM)和有机碳(TOC)采用重铬酸钾容量法-外加热法测定, 全氮(TN)采用半微量凯氏定氮法测定, 全磷(TP)采用碱熔-钼锑抗分光光度法测定, 氨氮(NH4+-N)采用KCl浸提-靛酚蓝比色法测定[13], 硝酸盐氮(NO3--N)采用紫外分光光度法(GB/T 32737-2016)测定, 亚硝酸盐氮(NO2--N)采用氯化钾溶液提取-分光光度法(HJ 634-2012)测定[14].

1.3 宏基因组学分析

将采集的沉积物样品在低温条件下送至上海某生物医药科技有限公司进行宏基因组学分析.主要的测试步骤包括:采用Fast DNA SPIN Kit for Soil试剂盒(MP bio, CA, USA)进行样品DNA抽提;完成DNA抽提后, 检测DNA浓度和纯度, 利用1%琼脂糖凝胶电泳检测DNA完整性;使用Covaris M220将检测合格的DNA剪切为约400 bp的片段, 用于构建PE文库;文库构建完成后进行桥式PCR, 最后利用Illumina NovaSeq 6000测序平台进行高通量测序.

宏基因测序数据分析参照曲疆奇等的方法[15], 高通量测序后, 首先使用fastp进行质量剪切, 保留高质量的pair-end reads和single-end reads;然后使用MEGAHIT对优化序列进行拼接组装;使用MetaGene对拼接结果中的contigs进行ORFs预测;使用CD-HIT软件将所有样品中预测出来的基因序列进行聚类(默认参数为:90%一致性, 90%覆盖度), 构建非冗余基因集;使用SOAPaligner软件分别将每个样品的高质量reads与非冗余基因集进行比对(95%一致性), 统计基因在对应样品中的丰度信息;使用Diamond将氮代谢基因集的氨基酸序列与KEGG数据库进行比对(BLASTP比对参数设置期望值e-value为1e-5), 获得氮代谢物种和功能注释, 然后基于物种对应的基因丰度总和计算该物种的相对丰度.

1.4 数据处理

本文基于代谢途径建立了氮代谢基因集, 该基因集包含所有氮代谢微生物.利用Excel 2019软件进行数据处理, 使用Origin 2021软件绘制功能基因丰度图, 通过美吉生物云平台(https://cloud.majorbio.com)计算氮代谢功能基因与环境因子间的Pearson相关关系并进行可视化分析, 采用Tbtools软件绘制门水平和属水平氮代谢微生物群落丰度图, 采用R语言psych包生成微生物与氮代谢功能基因相关性矩阵, 再利用Gephi软件绘制二者网络图, 进行可视化分析.

2 结果与讨论 2.1 沉积物理化指标分析

红枫湖沉积物的理化性质如表 1所示, 沉积物的pH在8.08~8.36之间, 平均值为8.24, 总体上呈弱碱性.EC为454.0~1 324.0 μS·cm-1, ω(TP)为457.96~1 186.09 mg·kg-1, ω(TN)为1 742.83~4 134.17 mg·kg-1, ω(SOM)为44.35~153.37 g·kg-1.与贵阳市的另外一座典型湖泊阿哈水库[16]相比, 红枫湖沉积物有机质和总氮均偏低, 总磷略高.采样点H1的有机质含量最高, 湖泊沉积物有机质来源分为内源和外源, 此处人为活动少, 因此可能主要是内源输入, 如水生植物腐烂和鱼类、底栖动物、浮游动物的代谢产物长期积累[17].除采样点H5以外, 各点亚硝酸盐氮明显高于氨氮和硝酸盐氮, 表明亚硝酸盐氮是红枫湖沉积物中主要的无机氮存在形式.且pH值越高的位置亚硝酸盐浓度普遍更高, 如采样点H3和H6, 这一规律也与徐亚同[18]的研究结果相似, 即当pH偏离最适反硝化的值时, 反硝化速率降低, 亚硝酸盐逐渐积累.采样点H5附近有大量农田、种植园、农场和一家农业科技企业, 人类的大量活动使得此处沉积物TN、TP和氨氮含量相对较高.

表 1 沉积物基本理化性质 Table 1 Basic physical and chemical properties of sediments

2.2 沉积物氮代谢功能基因

在红枫湖沉积物中共检测出6条氮代谢途径, 分别为固氮、硝化、反硝化、同化硝酸盐还原、异化硝酸盐还原和完全硝化, 如图 2所示.氮代谢功能基因的丰度大小表现为:反硝化 > 异化硝酸盐还原 > 完全硝化 > 硝化 > 同化硝酸盐还原 > 固氮, 且反硝化基因丰度远高于固氮基因丰度, 可能是由于硝酸盐氮的去除主要依赖于反硝化和异化硝酸盐还原过程, 且以反硝化过程为主.丰度前3的氮代谢功能基因分别是narGnrfAnapA, 相关研究表明较高丰度的narGnapA可以促进硝酸盐还原为亚硝酸盐, 且这两种基因广泛存在于反硝化菌中, 在反硝化过程中占主导地位[19].当河流和湖泊等汇入海洋系统之前, 接近一半的活性氮在反硝化的作用下, 以N2的形式释放到大气中[20].反硝化过程作为去除水体中硝态氮的主要途径之一, 能够将硝酸盐氮永久地转化为氮气[21], 这一作用对于饮用水源地来说是非常有益的, 有利于湖泊水体中氮素指标维持在良好的水平.

图 2 红枫湖沉积物中氮代谢功能基因丰度 Fig. 2 Abundance of functional genes of nitrogen metabolism in sediments of Hongfeng Lake

此外, 在红枫湖沉积物中未检测到厌氧氨氧化途径的功能基因, 韩雪恪等[22]总结了厌氧氨氧化过程的影响因素主要有底物浓度、有机物、溶解氧、温度和pH值等.结合当时的季节和红枫湖的实际情况, 表明这可能与高浓度亚硝酸盐和高浓度TOC的作用有关.红枫湖沉积物ω(NO2--N)平均值为235.88 mg·kg-1, Strous等[23]的研究表明在亚硝酸盐浓度高于100 mg·L-1时, 就会完全抑制厌氧氨氧化过程.红枫湖沉积物ω(TOC)平均值为48.14 g·kg-1, 与其他淡水湖泊相比, 如平寨水库[24]、太湖[25]和鄱阳湖[26, 27]等, 沉积物有机质的TOC含量较高.根据高美云等[28]研究的结果可知, 随着有机碳的含量增高, 对厌氧氨氧化过程具有抑制作用, 同时减少了厌氧氨氧化菌的活性.由于多个因素的综合作用, 使得厌氧氨氧化过程受到抑制, 因而检测不出相关的功能基因.

2.3 氮代谢微生物群落结构

红枫湖沉积物中细菌、古菌和真核生物的群落结构组成如图 3所示.细菌共获得41门、90纲、158目、245科和413属, 细菌是介导红枫湖沉积物氮代谢的主要微生物类群.相对丰度前3的优势菌门分别是变形菌门(Proteobacteria, 72.77%)、硝化螺旋菌门(Nitrospirae, 5.73%)和未分类的细菌(unclassified_d__Bacteria, 4.16%).相对丰度前3的优势菌属分别是硫杆菌属(Thiobacillus, 16.64%)、红长命菌属(Rubrivivax, 9.43%)和硝化螺旋菌属(Nitrospira, 7.09%).而岱海沉积物中的优势菌属为硫杆菌属(Thiobacillus)、丙酸杆菌属(Gammaproteobacteria)和暖绳菌属(Caldilinea[29], 二者的优势菌属具有一定差异. 岱海属于咸水湖, 且位于中高纬度地区, 而红枫湖属于淡水湖泊, 且位于喀斯特低纬度地区, 因此湖泊的生境会使沉积物中的微生物存在差异.

(a)门水平细菌丰度, (b)门水平古菌和真核生物丰度, (c)属水平细菌丰度, (d)属水平古菌和真核生物丰度;仅展示细菌、古菌和真核生物在各氮代谢途径丰度前5的门和丰度前10的属, 其余归为others;门属水平物种丰度进行lg(X+1)转化;A.完全硝化路径, B.硝化路径, C.反硝化路径, D.异化硝酸盐还原路径, E.固氮, F.同化硝酸盐还原路径;1、2、3、4、5、6、7和8对应采样点H1、H2、H3、H4、H5、H6、H7和H8 图 3 氮代谢微生物群落丰度 Fig. 3 Abundance of nitrogen-metabolizing microbial community

固氮路径的细菌群落多样性最低, 其中, 甲基杆菌属(Methylobacter)为固氮途径优势菌属, 在H5样点占比最高, 达到45.7%, 该菌属是湖泊中通常存在的好氧甲烷氧化菌, 可以在水生系统的食物网中发挥关键作用, 促进甲烷衍生碳向更高营养水平的转移[30].该结果与草海沉积物的优势甲烷氧化菌菌属[31]和三峡库区夏季万州段底泥中优势甲烷氧化菌菌属[32]一致.

反硝化和异化硝酸盐还原路径具有丰富的细菌群落组成, 且两个路径的细菌群落组成相似, 主要优势菌属均为硫杆菌属(Thiobacillus)、红长命菌属(Rubrivivax)和硝化螺旋菌属(Nitrospira).Thiobacillus中的细菌具有自养反硝化功能, 在厌氧或缺氧的条件下能以硫元素或硫化物为电子供体, 以NO3-为电子受体, 通过自养反硝化作用, 可以起到良好的硫、氮脱除效果[33]. Rubrivivax在氮代谢路径中作为优势菌属的研究报道较少, 能在光照条件下具有脱氮的功能[34], 可能对红枫湖水体的脱氮起到了至关重要的作用. Nitrospira具有硝化作用, 可将亚硝酸盐氧化成硝酸盐.以上结果表明氮、硫代谢在红枫湖沉积物生态系统中联系紧密, 并为氮硫的去除做出了较大贡献.此外, 嗜氢菌属(Hydrogenophilus)和硫单胞菌属(Thiomonas)分别是反硝化和异化硝酸盐还原路径的特有菌属, 且以上菌属只存在于部分采样点.Hydrogenophilus多在石油、煤矿区等有机物丰富的环境以优势菌属形式存在, 属于有机污染物降解菌[35].同化硝酸盐还原路径优势菌属主要是Parasulfuritortus和硫杆菌属(Thiobacillus), 且在湖区分布均匀.

古菌共获得7门、10纲、12目、14科和18属, 其中, 异化硝酸盐还原路径相比其他5条路径拥有更丰富的古菌群落, Candidatus_Thermoplasmatota是该途径主要的优势古菌门.Candidatus_Prometheoarchaeum是异化硝酸盐还原路径的特有古菌属.

除了细菌和古菌外, 真核生物也参与了红枫湖沉积物的氮循环, 分别是异化硝酸盐还原路径、反硝化路径和完全硝化路径. 根霉菌属(Rhizopus)是完全硝化路径的优势真核生物属, 在此路径中占比为28.6%, 属于毛霉菌门(Mucoromycota).蜉蝣属(Ephemera)是异化硝酸盐还原路径和反硝化路径共同的的优势真核生物属, 在以上路径中占比分别为4.7%和58.1%, 属于节肢动物门(Arthropoda).根据已有的证据, 真核生物可能有助于氮循环中的还原部分, 以上证据来自原生动物属(Loxodes)的纤毛虫, 它们生活在淡水湖中缺氧的深度以下[36, 37].Kamp等[38]认为在反硝化和异化硝酸盐还原途径中, 真核生物的作用很大, 在某些特殊的环境中, 它们的贡献甚至超过了原核生物.

2.4 沉积物环境因子与氮代谢功能基因的关联分析

为了解红枫湖沉积物样本中的氮代谢功能基因与主要环境因子(SOM、NO2--N、pH、TN、NH4+-N、TP、NO3--N和EC)的关系, 进行Pearson相关性分析, 结果如图 4所示.SOM、NO2--N、pH和TN与红枫湖沉积物中氮代谢功能基因之间相关性不显著.在氮代谢基因中, NO3-同化还原成NO2-过程中的narB基因仅与NO3--N呈正相关(P < 0.05);NH4+氧化为NH2OH过程中的pmoA-amoA基因与NH4+-N和TP正相关(P < 0.05), 与NO3--N负相关(P < 0.05);NO2-同化还原成NH4+过程中的nirA基因仅与NH4+-N正相关(P < 0.05);NH4+氧化为NH2OH过程中的pmoB-amoB基因与NH4+-N正相关(P < 0.05);N2氧化为NH4+过程中的nifK基因仅与NH4+-N正相关(P < 0.05), 与NO3--N负相关(P < 0.05);N2氧化为NH4+过程中的nifD基因与NH4+-N显著正相关(P < 0.01), 与NO3--N显著负相关(P < 0.01), 与TP正相关(P < 0.05);N2氧化为NH4+过程中的nifH基因仅与TP正相关(P < 0.05);NO2-还原成NO过程中的nirS基因仅与NH4+-N正相关(P < 0.05);NH2OH氧化为NO2-过程的hao基因仅与NO3--N正相关(P < 0.05);NO2-异化还原成NH4+过程中的nrfA基因仅与EC负相关(P < 0.05).说明显著影响红枫湖沉积物中氮代谢功能基因分布的主要环境因子为NH4+-N、NO3--N和TP, 这与谭秋阳等[39]的研究结果相似.

∗表示P<0.05, ∗∗表示P<0.01 图 4 红枫湖沉积物中氮代谢功能基因与环境因子的关联分析 Fig. 4 Relationships between functional genes and factors of nitrogen metabolism in sediments of Hongfeng Lake

nirSnirK是反硝化路径中编码亚硝酸盐还原酶(Nir)的功能基因, Nir能将亚硝酸盐还原为一氧化氮, 是反硝化过程的标志性反应[40].本研究中, 仅nirS与NH4+-N呈正相关, NO3--N和TN与nirSnirK基因并无显著相关性, 而在俞琳等[41]的研究中得出NO3--N含量是影响河流沉积物反硝化作用的重要因素, 反硝化速率与NO3--N含量呈显著正相关关系.在李佳霖等[42]的研究中得出弥河的反硝化作用与沉积物的TN呈正相关.以上差异说明NO3--N和TN对反硝化过程及其功能基因的影响在不同的环境条件下响应不同.NO3--N和NH4+-N对固氮路径功能基因(nifDnifKnifH)的影响表现出一致的相关性, 且影响最大.TP对固氮基因也表现出一致的相关性, 杨文焕等[43]研究发现, 较高的C/N和TP是固氮菌生存的必要营养条件, 可以促进固氮反应的顺利进行.NH4+-N的增加会抑制固氮微生物的生长, 因此NH4+-N和TP会显著影响固氮作用相关基因丰度.在生物脱氮过程中, 由氨单加氧酶将铵根(NH4+)氧化为羟胺(NH2OH), 氨单加氧酶是由amoAamoBamoC这3个基因编码的αβγ这3个亚基组成[44], 且amoA基因被广泛用于生态系统中氨氧化菌AOA和AOB的研究[45].羟胺是完全硝化过程的中间产物, hao是编码羟胺氧化还原酶的基因, 可以将羟胺氧化为亚硝酸盐[46].相关研究表明, 沉积物NH4+-N和TP对硝化作用影响较大, 是硝化过程的关键影响因素[47].而磷是微生物生长所必须的营养元素, 沉积物中磷的存在有利于氨氧化细菌的生存, 并且磷酸盐的存在也能够形成磷酸盐缓冲体系, 可暂时性延缓因硝化作用造成的沉积物pH值下降, 有利于硝化作用的进行[48].

2.5 微生物与氮代谢功能基因共现性网络分析

在属水平下, 基于Spearman相关性系数(P)计算, 选择0.8作为相关系数阈值, 构建总丰度前40的氮代谢微生物(不包括未分类的细菌)与氮代谢功能基因的共现网络, 分析氮代谢功能基因与微生物的共现性关系[49]图 5), 并绘制出红枫湖沉积物中氮循环过程示意图(图 6).共现网络由60个节点和141条边组成, 平均加权度为11.267, 密度为0.08, 模块化指数为0.481(大于0.4), 证明网络图具有较强的模块化结构, 对环境变化具有较高的抵抗力[50].模块性分析发现氮代谢微生物与功能基因共现性网络存在6个模块, 前3个模块占比高达75%, 主要归属于变形菌门, 变形菌门是水环境中主要的原核生物[51], 具有丰富的生态功能, 能够通过介导硫、甲烷氧化, 硫酸盐还原和反硝化等作用驱动碳、氢、硫、氮等物质循环[52].功能基因与微生物之间正相关性(95.04%)的占比远大于负相关性(4.96%), 样本中丰度较高的功能基因如norBnorCnarInarVnasAnapASulfurisoma呈显著正相关, nirSZoogloeaLysobacter呈显著正相关, narGnarZnxrAPelomonas呈显著正相关, 以上结果意味着在该系统中氮代谢微生物的主要生存方式为互利共生.vnfK与厌氧粘菌属(Anaeromyxobacter)呈显著负相关, 虽然二者丰度不高, 但为网络中负相关最强的边.相关研究表明, 大多数不利连接被认为是由不相容的生态位或直接的食物竞争引起的[53]. narGnarZnxrAnarHnarYnxrBnarInarV连接程度较高, 为该网络的枢纽基因[54]. narGnarZnxrA这三者属于同源基因[1], 在该网络中宿主最多, 共有11个属.与narG相关性最强的前3种微生物分别是Pelomonas(0.95)、Sulfuritalea(0.90)和Sulfurisoma(0.88), 其中嗜糖假单胞菌属(Pelomonas)参与了同化硝酸盐还原途径、反硝化途径和厌氧氨氧化途径, 在亚硝酸盐向一氧化氮的转化中发挥了重要作用.

每个节点大小与连接数成正比, 红色连接表示正相关, 绿色连接表示负相关, 括号内的数值表示每个模块的占比 图 5 属水平上氮代谢功能基因与微生物群落关系的网络分析 Fig. 5 Network analysis of the relationship between nitrogen metabolism functional genes and microbial community at the genus level

图 6 红枫湖沉积物氮代谢过程示意 Fig. 6 Schematic of nitrogen metabolism in sediments of Hongfeng Lake

3 结论

(1)红枫湖沉积物中氮代谢微生物主要为细菌, 其中相对丰度前3的优势菌属分别为硫杆菌属(Thiobacillus, 16.64%)、红长命菌属(Rubrivivax, 9.43%)和硝化螺旋菌属(Nitrospira, 7.09%).古菌在异化硝酸盐还原途径丰度最高, 真核生物参与了异化硝酸盐还原途径、反硝化途径和完全硝化途径.

(2)红枫湖沉积物中存在6条氮代谢途径, 分别为固氮、硝化、反硝化、同化硝酸盐还原、异化硝酸盐还原和完全硝化.丰度前3的氮代谢功能基因分别是narGnrfAnapA.反硝化和异化硝酸盐还原是红枫湖沉积物中氮去除的主要途径.

(3)氮代谢功能基因主要的环境影响因子为NH4+-N、NO3--N和TP, 与其他环境的驱动因子存在差异.微生物与氮代谢功能基因之间主要呈正相关关系, 占比为95.04%.氮代谢枢纽基因分别为narGnarZnxrAnarHnarYnxrBnarInarV.

参考文献
[1] 张博雅, 余珂. 微生物基因数据库在氮循环功能基因注释中的应用[J]. 微生物学通报, 2020, 47(9): 3021-3038.
Zhang B Y, Yu K. Application of microbial gene databases in the annotation of nitrogen cycle functional genes[J]. Microbiology China, 2020, 47(9): 3021-3038.
[2] 张亚平, 万宇, 聂青, 等. 湖泊水体中氮的生物地球化学过程及其生态学意义[J]. 南京大学学报(自然科学), 2016, 52(1): 5-15.
Zhang Y P, Wan Y, Nie Q, et al. Biogeochemical processes of nitrogen in lake systems and their ecological importances[J]. Journal of Nanjing University (Natural Science), 2016, 52(1): 5-15.
[3] 徐婷婷, 郑俊强, 韩士杰, 等. 长白山阔叶红松林土壤氮转化过程对长期施氮和降水变化的响应[J]. 应用生态学报, 2018, 29(9): 2797-2807.
Xu T T, Zheng J Q, Han S J, et al. Responses of soil nitrogen transformation to long-term nitrogen fertilization and precipitation changes in a broad-leaved Korean pine forest in Changbai Mountains, China[J]. Chinese Journal of Applied Ecology, 2018, 29(9): 2797-2807.
[4] Yu T T, Li M, Niu M Y, et al. Difference of nitrogen-cycling microbes between shallow bay and deep-sea sediments in the South China Sea[J]. Applied Microbiology and Biotechnology, 2018, 102(1): 447-459. DOI:10.1007/s00253-017-8594-9
[5] Li J Y, Katsev S. Nitrogen cycling in deeply oxygenated sediments: Results in Lake Superior and implications for marine sediments[J]. Limnology and Oceanography, 2014, 59(2): 465-481. DOI:10.4319/lo.2014.59.2.0465
[6] Feng L J, Zhang Z L, Yang G F, et al. Microbial communities and sediment nitrogen cycle in a coastal eutrophic lake with salinity and nutrients shifted by seawater intrusion[J]. Environmental Research, 2023, 225. DOI:10.1016/j.envres.2023.115590
[7] 孙小溪. 青藏高原湖泊微生物介导的氮移除过程及其环境效应[D]. 武汉: 中国地质大学, 2022.
Sun X X. Microbially mediated nitrogen removal process and its environmental effect in the Qing-Tibetan lakes[D]. Wuhan: China University of Geosciences, 2022.
[8] Galloway J N, Townsend A R, Erisman J W, et al. Transformation of the nitrogen cycle: Recent trends, questions, and potential solutions[J]. Science, 2008, 320(5878): 889-892. DOI:10.1126/science.1136674
[9] Canfield D E, Glazer A N, Falkowski P G. The evolution and future of Earth's nitrogen cycle[J]. Science, 2010, 330(6001): 192-196. DOI:10.1126/science.1186120
[10] Bernhardt E S. Cleaner lakes are dirtier lakes[J]. Science, 2013, 342(6155): 205-206. DOI:10.1126/science.1245279
[11] 赵宝刚, 张夏彬, 昝逢宇, 等. 不同湖泊表层沉积物氮形态的分布特征与影响因素[J]. 中国环境科学, 2021, 41(2): 837-847.
Zhao B G, Zhang X B, Zan F Y, et al. Distribution characteristics and influencing factors of nitrogen forms in surface sediments of different lakes[J]. China Environmental Science, 2021, 41(2): 837-847. DOI:10.3969/j.issn.1000-6923.2021.02.038
[12] 鲍林林, 陈永娟, 王晓燕. 河流沉积物氮循环主要微生物的生态特征[J]. 微生物学通报, 2015, 42(6): 1141-1150.
Bao L L, Chen Y J, Wang X Y. Ecological characteristics of nitrogen cycling microbes in river sediments[J]. Microbiology China, 2015, 42(6): 1141-1150.
[13] 鲍士旦. 土壤农化分析[M]. ((第三版)). 北京: 中国农业出版社, 2000.
[14] Chen X F, Wang K, Li X, et al. Microcystis blooms aggravate the diurnal alternation of nitrification and nitrate reduction in the water column in Lake Taihu[J]. Science of the Total Environment, 2021, 767. DOI:10.1016/j.scitotenv.2020.144884
[15] 曲疆奇, 张清靖, 杨浩辰, 等. 养殖尾水处理系统填料生物膜氮循环微生物功能特征[J]. 水生态学杂志, 2023, 44(2): 104-113.
Qu J Q, Zhang Q J, Yang H C, et al. Functional characteristics of nitrogen cycling microbes in the biofilm carrier of an aquaculture tailwater treatment system[J]. Journal of Hydroecology, 2023, 44(2): 104-113.
[16] 严莲英, 范成五, 刘桂华, 等. 喀斯特山区水库底泥养分与重金属的含量特征-以阿哈水库为例[J]. 安全与环境学报, 2017, 17(5): 2031-2037.
Yan L Y, Fan C W, Liu G H, et al. Case study of the characteristic features of the nutritious and heavy metal contents in sediments of the Karst rock mountain reservoir in Guizhou[J]. Journal of Safety and Environment, 2017, 17(5): 2031-2037.
[17] 张小春, 代丹, 秦成, 等. 梅溪湖表层沉积物营养盐和有机质分布特征及污染评价[J]. 水生态学杂志, 2023, 44(5): 107-113.
Zhang X C, Dai D, Qin C, et al. Distribution and pollution assessment of nutrients and organic matter in surface sediments of Meixi Lake[J]. Journal of Hydroecology, 2023, 44(5): 107-113.
[18] 徐亚同. pH值、温度对反硝化的影响[J]. 中国环境科学, 1994, 14(4): 308-313.
Xu Y T. The influence of pH values and temperature on denitrification[J]. China Environmental Science, 1994, 14(4): 308-313. DOI:10.3321/j.issn:1000-6923.1994.04.010
[19] Zeng W, Zhang J, Wang A Q, et al. Denitrifying phosphorus removal from municipal wastewater and dynamics of "Candidatus Accumulibacter" and denitrifying bacteria based on genes of ppk1, narG, nirS and nirK[J]. Bioresource Technology, 2016, 207: 322-331. DOI:10.1016/j.biortech.2016.02.016
[20] Galloway J N, Aber J D, Erisman J W, et al. The nitrogen cascade[J]. BioScience, 2003, 53(4): 341-356. DOI:10.1641/0006-3568(2003)053[0341:TNC]2.0.CO;2
[21] Seitzinger S, Harrison J A, Böhlke J K, et al. Denitrification across landscapes and waterscapes: A synthesis[J]. Ecological Applications, 2006, 16(6): 2064-2090. DOI:10.1890/1051-0761(2006)016[2064:DALAWA]2.0.CO;2
[22] 韩雪恪, 王峥嵘, 彭永臻, 等. 几种重要因素对厌氧氨氧化过程的影响综述[J]. 中国环境科学, 2023, 43(5): 2220-2227.
Han X K, Wang Z R, Peng Y Z, et al. A review of several important factors influencing the Anammox process[J]. China Environmental Science, 2023, 43(5): 2220-2227. DOI:10.3969/j.issn.1000-6923.2023.05.013
[23] Strous M, Kuenen J G, Jetten M S M. Key physiology of anaerobic ammonium oxidation[J]. Applied and Environmental Microbiology, 1999, 65(7): 3248-3250. DOI:10.1128/AEM.65.7.3248-3250.1999
[24] 邹艳, 周忠发, 孔杰, 等. 喀斯特水库及入库河流TOC、TN的变化特征及影响机制[J]. 地球与环境, 2023, 51(3): 255-265.
Zou Y, Zhou Z F, Kong J, et al. Variation characteristics and influence mechanisms of TOC and TN in karst reservoirs and inlet rivers[J]. Earth and Environment, 2023, 51(3): 255-265.
[25] Xue Q J, Steinman A D, Xie L Q, et al. Seasonal variation and potential risk assessment of microcystins in the sediments of Lake Taihu, China[J]. Environmental Pollution, 2020, 259. DOI:10.1016/j.envpol.2019.113884
[26] 王毛兰, 赖建平, 胡珂图, 等. 鄱阳湖表层沉积物有机碳、氮同位素特征及其来源分析[J]. 中国环境科学, 2014, 34(4): 1019-1025.
Wang M L, Lai J P, Hu K T, et al. Compositions and sources of stable organic carbon and nitrogen isotopes in surface sediments of Poyang Lake[J]. China Environmental Science, 2014, 34(4): 1019-1025.
[27] 金笑, 寇文伯, 于昊天, 等. 鄱阳湖不同区域沉积物细菌群落结构、功能变化及其与环境因子的关系[J]. 环境科学研究, 2017, 30(4): 529-536.
Jin X, Kou W B, Yu H T, et al. Environmental factors influencing the spatial distribution of sediment bacterial community structure and function in Poyang Lake[J]. Research of Environmental Sciences, 2017, 30(4): 529-536.
[28] 高美云, 刘兴国, 曾宪磊, 等. 有机碳对养殖池塘沉积物中反硝化、厌氧氨氧化的影响[J]. 环境工程学报, 2018, 12(1): 49-56.
Gao M Y, Liu X G, Zeng X L, et al. Effects of organic carbon for denitrification and anaerobic ammonium oxidation in sediments of aquaculture pond[J]. Chinese Journal of Environmental Engineering, 2018, 12(1): 49-56.
[29] 尹强. 岱海沉积物中微生物群落结构及功能特征研究[D]. 包头: 内蒙古科技大学, 2021.
Yin Q. Microbial community structure and function in Daihai sediments feature research[D]. Baotou: Inner Mongolia University of Science & Technology, 2021.
[30] Van Grinsven S, Damsté J S S, Harrison J, et al. Nitrate promotes the transfer of methane-derived carbon from the methanotroph Methylobacter sp. to the methylotroph Methylotenera sp. in eutrophic lake water[J]. Limnology and Oceanography, 2020, 66(3): 878-891.
[31] 艾佳, 吕杨, 钟雄, 等. 贵州草海喀斯特高原湖泊湿地甲烷氧化菌群落特征及功能探析[J]. 湖泊科学, 2022, 34(3): 906-918.
Ai J, Lv Y, Zhong X, et al. The community characteristics and functions of methanotrophs in karst Lake Caohai, Guizhou Plateau[J]. Journal of Lake Sciences, 2022, 34(3): 906-918.
[32] 彭兴意, 秦宇, 舒钰清, 等. 三峡库区夏季万州段底泥甲烷功能菌群落对甲烷排放的影响[J]. 环境工程学报, 2022, 16(3): 1028-1038.
Peng X Y, Qin Y, Shu Y Q, et al. Effects of functional methane community in sediments of Wanzhou section of the Three Gorges Reservoir on methane emissions in summer[J]. Chinese Journal of Environmental Engineering, 2022, 16(3): 1028-1038.
[33] 武鑫, 赵晋忠. 脱氮硫杆菌对废水中硫氮的脱除[J]. 应用与环境生物学报, 2013, 19(1): 175-178.
Wu X, Zhao J Z. Removal of nitrate and sulfide from wastewater by Thiobacillusdenitrificans[J]. Chinese Journal of Applied and Environmental Biology, 2013, 19(1): 175-178.
[34] 康鹏亮, 黄廷林, 张海涵, 等. 西安市典型景观水体水质及反硝化细菌种群结构[J]. 环境科学, 2017, 38(12): 5174-5183.
Kang P L, Huang T L, Zhang H H, et al. Water quality and diversity of denitrifier community structure of typical scenic water bodies in Xi'an[J]. Environmental Science, 2017, 38(12): 5174-5183.
[35] 秦胜红, 李新宇, 李旭, 等. 不同土壤细菌种群结构对氯嘧磺隆胁迫的响应及降解菌系的获得[J]. 生态学杂志, 2012, 31(7): 1724-1732.
Qin S H, Li X Y, Li X, et al. Responses of bacterial community structure in different soils to chlorimuron-ethyl and acquisition of chlorimuron-ethyl degrading bacterial population[J]. Chinese Journal of Ecology, 2012, 31(7): 1724-1732.
[36] Finlay B J. Nitrate respiration by Protozoa (Loxodes spp.) in the hypolimnetic nitrite maximum of a productive freshwater pond[J]. Freshwater Biology, 1985, 15(3): 333-346. DOI:10.1111/j.1365-2427.1985.tb00205.x
[37] Finlay B J, Span A S W, Harman J M P. Nitrate respiration in primitive eukaryotes[J]. Nature, 1983, 303(5915): 333-336. DOI:10.1038/303333a0
[38] Kamp A, Høgslund S, Risgaard-Petersen N, et al. Nitrate storage and dissimilatory nitrate reduction by eukaryotic microbes[J]. Frontiers in Microbiology, 2015, 6. DOI:10.3389/fmicb.2015.01492
[39] 谭秋阳, 吕金燕, 白献宇, 等. 典型城市内河菌群组成与氮循环功能垂向分布及溯源分析[J]. 环境科学, 2023, 44(3): 1497-1507.
Tan Q Y, Lv J Y, Bai X Y, et al. Vertical distribution and source tracking analysis of bacteria composition and nitrogen metabolism function of a typical urban inland river[J]. Environmental Science, 2023, 44(3): 1497-1507.
[40] 王杨. 不同酸度土壤硝化和反硝化活性的差异[D]. 大连: 大连交通大学, 2014.
Wang Y. Difference in the activity of nitrification and denitrification with different soil acidity[D]. Dalian: Dalian Jiaotong University, 2014.
[41] 俞琳, 王东启, 邓焕广, 等. 上海市苏州河与崇明岛河流沉积物反硝化速率及其影响因素[J]. 生态与农村环境学报, 2015, 31(1): 100-105.
Yu L, Wang D Q, Deng H G, et al. Denitrification rates in sediments of Suzhou River and rivers in Chongming Island of Shanghai and their affecting factors[J]. Journal of Ecology and Rural Environment, 2015, 31(1): 100-105.
[42] 李佳霖, 秦松. 弥河沉积物的反硝化和厌氧氨氧化过程[J]. 中国环境科学, 2021, 41(4): 1588-1596.
Li J L, Qin S. Denitrification and anammox processes in sediment of Mihe River, China[J]. China Environmental Science, 2021, 41(4): 1588-1596. DOI:10.3969/j.issn.1000-6923.2021.04.011
[43] 杨文焕, 张元, 王志超, 等. 寒旱区湖泊沉积物中固氮微生物群落特征-以包头南海湖为例[J]. 中国环境科学, 2020, 40(6): 2674-2682.
Yang W H, Zhang Y, Wang Z C, et al. Community characteristics of nitrogen-fixing microorganisms in lake sediment-taking Nanhaihu lake as example[J]. China Environmental Science, 2020, 40(6): 2674-2682. DOI:10.3969/j.issn.1000-6923.2020.06.038
[44] 潘彦羽, 代嫣然, 王飞华, 等. 东湖表层沉积物中氨氧化古菌和氨氧化细菌丰度及多样性研究[J]. 水生生物学报, 2018, 42(2): 406-415.
Pan Y Y, Dai Y R, Wang F H, et al. Abundance and diversity of ammonia-oxidizing archaea and ammonia-oxidizing bacteria in the surface sediments of Donghu Lake[J]. Acta Hydrobiologica Sinica, 2018, 42(2): 406-415.
[45] Petersen D G, Blazewicz S J, Firestone M, et al. Abundance of microbial genes associated with nitrogen cycling as indices of biogeochemical process rates across a vegetation gradient in Alaska[J]. Environmental Microbiology, 2012, 14(4): 993-1008. DOI:10.1111/j.1462-2920.2011.02679.x
[46] 吴启凤, 丁晨雨, 何腾霞, 等. 羟胺氧化还原酶在生物脱氮中的研究进展[J]. 微生物学通报, 2022, 49(3): 1120-1134.
Wu Q F, Ding C Y, He T X, et al. Research progress of hydroxylamine oxidoreductase in biological denitrification[J]. Microbiology China, 2022, 49(3): 1120-1134.
[47] 白洁, 陈春涛, 赵阳国, 等. 辽河口湿地沉积物硝化细菌及硝化作用研究[J]. 环境科学, 2010, 31(12): 3011-3017.
Bai J, Chen C T, Zhao Y G, et al. Studies on nitrobacteria and nitrification in Liaohe estuary wetland sediments[J]. Environmental Science, 2010, 31(12): 3011-3017.
[48] 代惠萍, 顾斌, 屈继旗, 等. 铵态氮和水溶性磷酸盐配比对土壤硝化的影响[J]. 西北农业学报, 2009, 18(5): 189-193.
Dai H P, Gu B, Qu J Q, et al. Effect of match ratio of ammonium nitrogen and water solubility phosphas on soil nitrification[J]. Acta Agriculturae Boreali-Occidentalis Sinica, 2009, 18(5): 189-193. DOI:10.3969/j.issn.1004-1389.2009.05.041
[49] Liu D P, Zhang Y D, Yang Q L, et al. Fate of ofloxacin in rural wastewater treatment facility: Removal performance, pathways and microbial characteristics[J]. Bioresource Technology, 2023, 371. DOI:10.1016/j.biortech.2023.128611
[50] Newman M E J. Modularity and community structure in networks[J]. Proceedings of the National Academy of Sciences of the United States of America, 2006, 103(23): 8577-8582.
[51] 张雅洁, 李珂, 朱浩然, 等. 北海湖微生物群落结构随季节变化特征[J]. 环境科学, 2017, 38(8): 3319-3329.
Zhang Y J, Li K, Zhu H R, et al. Community structure of microorganisms and its seasonal variation in Beihai Lake[J]. Environmental Science, 2017, 38(8): 3319-3329.
[52] Zhou Z C, Tran P Q, Kieft K, et al. Genome diversification in globally distributed novel marine Proteobacteria is linked to environmental adaptation[J]. The ISME Journal, 2020, 14(8): 2060-2077. DOI:10.1038/s41396-020-0669-4
[53] Zhu M H, Zhang M L, Yuan Y B, et al. Responses of microbial communities and their interactions to ibuprofen in a bio-electrochemical system[J]. Journal of Environmental Management, 2021, 289. DOI:10.1016/j.jenvman.2021.112473
[54] Zhang F, Wang L Y, Bai P X, et al. Identification of regulatory networks and hub genes controlling nitrogen uptake in tea plants [Camellia sinensis (L.) O. Kuntze][J]. Journal of Agricultural and Food Chemistry, 2020, 68(8): 2445-2456. DOI:10.1021/acs.jafc.9b06427