环境科学  2022, Vol. 43 Issue (1): 314-328   PDF    
岗南水库沉积物好氧反硝化菌群落时空分布特征
张紫薇, 陈召莹, 张甜娜, 周石磊, 崔建升, 罗晓     
河北科技大学环境科学与工程学院, 河北省污染防治生物技术实验室, 石家庄 050018
摘要: 为研究岗南水库沉积物中好氧反硝化菌群的时空分布特征,运用好氧反硝化功能基因(napA)高通量测序技术,进行了α-多样性分析、微生物群落组成分析、微生物群落差异分析以及微生物网络分析.结果表明,α-多样性分析显示Chao1指数、Observed species指数和Shannon指数秋冬季均高于春夏季,即秋冬季采样点丰富度和多样性均较高.空间分布上,过渡区采样点的菌群的Chao1指数、Observed species指数、Shannon指数和Simpson指数均高于入库区和主库区采样点,即过渡区采样点菌群丰富度和多样性最高.岗南水库沉积物的好氧反硝化菌群主要属于变形菌门(Protebacteria),第一优势纲为β-Proteobacteria纲,第一优势属为陶厄氏菌属(Thauera).韦恩图分析(Venn)表明岗南水库沉积物各采样点好氧反硝化菌群分布均存在明显的时空差异.从空间水平基于LEfSe分析得到的差异物种比在季节水平的多.随机森林分析筛选得到标志物种中,Comamonas在所有空间分组采样点贡献程度最高,Pectobacterium在所有季节分组采样点贡献程度最高.网络分析表明存在9个主要模块包含297个节点565个边,关键物种分析得到47个关键物种OTUs.主坐标分析(principal co-ordinates analysis)显示关键物种的空间分布差异显著(Adonis,P < 0.001).综上,通过对岗南水库沉积物好氧反硝化菌群组成结构及分布规律的研究,为后期高效好氧反硝化菌分离鉴定提供了科学依据.
关键词: 岗南水库      沉积物      好氧反硝化      网络分析      生物多样性分析     
Spatial and Temporal Distribution of Aerobic Denitrification Bacterial Community in Sediments of Gangnan Reservoir
ZHANG Zi-wei , CHEN Zhao-ying , ZHANG Tian-na , ZHOU Shi-lei , CUI Jian-sheng , LUO Xiao     
Pollution Prevention Biotechnology Laboratory of Hebei Province, College of Environmental Sciences and Engineering, Hebei University of Science and Technology, Shijiazhuang 050018, China
Abstract: In order to study the spatial and temporal distribution characteristics of aerobic denitrification bacteria in the sediments of Gangnan Reservoir, the aerobic denitrification bacterial community was analyzed using a MiSeq high-throughput sequencing technique based on the napA gene. Moreover, the composition, diversity, difference, and network analysis of the aerobic denitrification bacterial community were investigated. The results showed higher α-diversity (Chao1 index, Observed species index, and Shannon index) in autumn and winter than that in spring and summer; that is, the richness and diversity of sampling sites in autumn and winter were higher. In terms of spatial distribution, the Chao1 index, Observed species index, Shannon index, and Simpson index of the bacterial community in the sampling sites in the transition zone were higher than those in the sampling sites in the inflow area and the main reservoir area, indicating that the richness and diversity of the bacterial community in the sampling sites in the transition zone were the highest. The aerobic denitrifying bacteria in the sediments of Gangnan Reservoir mainly belong to Protebacteria. The first dominant class was β-Proteobacteria, and the first dominant genus was Thauera. A Venn diagram analysis indicated that the community of aerobic denitrification bacteria exhibited significant spatial and temporal differences. At the spatial level, there were higher numbers of different species based on LEfSe analysis than that at the seasonal level. Among the marker species screened by random forest analysis, Comamonas had the highest contribution in all spatial sampling sites, whereas Pectobacterium had the highest contribution in all seasonal sampling sites. Network analysis showed that there were nine modules, including 565 edges of 297 nodes; 47 keystone species were obtained based on the indicator OTUs analysis and network analysis. The principal co-ordinates analysis revealed that keystone species in the spatial distribution exhibited significant differences (Adonis, P < 0.001). The results of this study provide scientific basis for the separation and identification of efficient aerobic denitrifying bacteria.
Key words: Gangnan Reservoir      sediment      aerobic denitrification      network analysis      biodiversity analysis     

随着我国对城市地下水开采的限制, 水库成为城市的主要供水水源[1], 到目前为止, 水源水库都存在一定程度的水质污染或富营养化问题, 例如: 缓丘地带黄河镇水库[2]、大连碧流河水库[3]、山东周村水库[4]和三峡水库[5], 严重的已经丧失作为供给水源的功能.从而导致水厂处理难度加大、制水成本与水质超标风险大幅提高, 严重威胁到城市的饮用水安全.水库中的沉积物是水生态系统中氮素物质的汇集地[6~8], 对污染物迁移转化发挥重要作用.沉积物中氮素的去除对水库内源污染控制有十分重要的意义.

岗南水库作为石家庄的饮用水源地, 其水库沉积物中富含有大量微生物, 且微生物在水-沉积物的物质循环中起重要作用[9], 因此沉积物中的微生物与内源污染负荷关系紧密.好氧反硝化菌群的出现, 对水源水库氮素的消减和水体功能恢复有指导意义.Robertson等[10]的研究在废水脱硫和反硝化系统中首次分离出1株好氧反硝化菌Thiosphaera pantotropha(脱氮副球菌, Paracoccus denitrificans), 该菌群不仅可以降低氮素含量和去除有机污染物[11], 还可以同时进行硝化和反硝化, 在水体内源污染的控制过程中发挥重要作用.近年来, 国内外学者对好氧反硝化菌群已经进行了大量的研究[11, 12], 其中针对水源水库沉积物中好氧反硝化菌群的研究主要集中在高效菌株的筛选分离[13]和脱氮效率[14], 但实验室培养基驯化的高效菌株不适于实际应用, 很大程度是因为不清楚实际环境中好氧反硝化菌的分布特征以及驱动自然界好氧反硝化群落演变的影响因子.而沉积物是好氧反硝化菌群的聚集地, 目前对水源水库沉积物好氧反硝化微生物时空分布特征分析较少, 因此探明岗南水库沉积物好氧反硝化菌群实际分布特征, 可以为高效好氧反硝化菌的定向筛选提供数据支持, 对污染负荷环境的脱氮以及富营养化控制有十分重要的意义.

本课题组前期基于好氧反硝化功能基因(napA)的高通量测序[15, 16], 研究了好氧反硝化微生物群落特征.本文以岗南水库沉积物为研究对象, 通过好氧反硝化功能基因(napA)的测序技术探究从入库河口到坝前主库区沉积物中好氧反硝化菌群的时空分布特征, 本文应用α-多样性分析、韦恩分析、LEfSe分析和随机森林分析阐明岗南水库沉积物好氧反硝化菌群结构组成、多样性和差异性, 结合微生物网络分析解析好氧反硝化菌群落演变的关键物种, 以期为进一步研究水源水库氮素消除、水质改善和饮用水水质保障提供参考.

1 材料与方法 1.1 研究区概况

岗南水库位于河北省平山县滹沱河干流, 是一座防洪灌溉、发电及库区养鱼等综合利用的大型水利枢纽工程.水库于1962年竣工, 总库容15.71亿m3, 控制流域面积15 900 km2, 流域属大陆性气候, 多年平均降雨量在地区分布上相差较大, 且80%以上集中在汛期6~9月, 多年平均气温为12.9~4.1℃, 多年平均蒸发量1 239 mm, 多年平均风速为1.7~9.0 m ·s-1; 多年平均径流量为16.6亿m3.

1.2 样品的采集与测定

本研究自入库河口到坝前主库区共选取6个代表性采样点采集沉积物样品, 根据水库空间位置, 将岗南水库划分为3个区域, 分别为入库区、过渡区和主库区.本文选取1号采样点代表入库区, 2号和3号采样点代表过渡区, 4号、5号和6号采样点代表主库区, 采样点位区域如图 1所示.分别于2018年9月(秋季)、2018年12月(冬季)、2019年3月(春季)和2019年6月(夏季)用彼得森采泥器采集表层沉积物(0~2 cm), 来研究岗南水库沉积物好氧反硝化菌群时空分布特征.

图 1 岗南水库沉积物采样点分布示意 Fig. 1 Distribution of sediment sampling sites in Gangnan Reservoir

1.3 DNA提取和测序

沉积物在派森诺生物公司进行微生物群落多样性分析, 基于napA功能基因对24个沉积物样品测序, 用V66F: 5′-TAYTTYYTNHSNAARATHATGTA YGG-3′和V67R: 5′-DATNGGRTGCATYTCNGCCAT RTT-3′对水库的napA基因进行扩增[16].进而分析岗南水库沉积物基于功能基因(napA)的好氧反硝化细菌的群落组成.采用QIIME2对原始数据进行拼接、过滤和去除嵌合体, 得到有效数据, 并将有效序列按照97%的一致性聚类成为OTU, 并采样Nt数据库对OTU进行注释, 用Vsearch进行数据分析.

1.4 测序分析

应用Chao1指数、Observed species指数、香农指数(Shannon index)、辛普森指数(Simpson index)、覆盖度(Goods coverage)指数和Pielou指数进行α-多样性分析.在派森诺基因云平台, 根据OTU丰度绘制Venn图[17], 进行物种差异性分析.基于LEfSe分析得到差异物种, 基于随机森林分析筛选标志物种.

1.5 网络分析

应用微生物网络可视化分析水源水库沉积物在不同时期好氧反硝化群落演变特征以及筛选出关键物种[18, 19].应用PCoA(Principal co-ordinates analysis)分析关键物种, 研究关键物种组成的相似性和差异性[20].利用置换多因素方差分析(Adonis analysis)研究关键物种在季节和空间分组在统计学上的差异[21].

1.6 数据分析方法

使用QIIME2软件对原始数据进行处理.使用ArcGIS软件完成采样图的绘制.数据处理由派森诺基因云平台完成.用Excel软件完成指数相关计算.通过R分析进行微生物网络的构建, 研究物种群落间相互协作关系.

2 结果与分析 2.1 微生物多样性分析 2.1.1 不同季节好氧反硝化细菌群生物多样性分析

本研究采用α-多样性来表征好氧反硝化细菌生物多样性.其中Chao1指数和Observed species指数均表征菌群丰富度, 指数越大说明菌群丰度越高.从季节的角度来看(图 2), Chao1指数呈现了一直上升过程, 从春季的2 616.03±164.55上升到夏季的2 623.28±361.52, 上升到秋季的2 690.25±449.72, 继续上升到冬季的2 775.43±211.69. Observed species指数呈现了先下降后上升再下降的过程, 从春季的1 730.05±84.23, 下降到夏季的1 722.18±238.52, 再上升到秋季的1 930.37±259.91, 再下降到冬季的1 885.48±181.76.秋冬季要比春夏季指数高, 表明秋冬季菌群丰富度更高.两个指数均没有显著差异(P>0.05).表征菌群多样性的Shannon指数从季节分布上看, 秋冬季指数高于春夏季, 其中秋季指数最高, 夏季最低, 表明秋季好氧反硝化菌群多样性最多, 夏季最低.Simpson指数在4个季节没有显著差异(P>0.05).说明岗南水库沉积物好氧反硝化菌群在秋季有稳定的生存环境.Goods coverage指数用来评估测序对群落中物种的覆盖度, 4个季节的指数均维持在0.98以上, 表明测序可以很好地覆盖物种信息; Pielou指数表征好氧反硝化菌群的均匀性, 4个季节指数季节差异为: 秋季>冬季>春季>夏季, 说明秋季物种分布最均匀, 反之, 夏季最低.由此可以看出岗南水库好氧反硝化细菌多样性季节演替显著, 且不同库区菌群多样性季节演替存在差异.

图 2 岗南水库沉积物不同季节好氧反硝化细菌群生物多样性指数 Fig. 2 Microbial community diversity of aerobic denitrification bacteria in different seasons in sediments of Gangnan Reservoir

水库主库区采样点的Chao1指数变化趋势为: 秋季(2 643.04)>春季(2 624.45)>冬季(2 617.55)>夏季(2 523.19).Observed species指数变化趋势为: 秋季(1 955.27)>冬季(1 742.90)>春季(1 674.47)>夏季(1 640.77).Shannon指数的变化趋势为秋季(7.73)>冬季(7.12)>春季(7.10)>夏季(6.93).Simpson指数变化趋势为秋季>春季>冬季>夏季. Pielou指数显示秋季最高, 夏季最低.说明主库区好氧反硝化菌群在秋季丰富度和多样性均最高, 夏季均最低.主库区的季节演替规律与岗南水库整体季节演替规律有重合点, 都是秋季多样性最高, 夏季最低.

2.1.2 不同空间好氧反硝化细菌群生物多样性分析

从空间的角度看(图 3), 过渡区采样点的菌群的Chao1指数和Observed species指数均高于入库区和主库区采样点, 入库区的2 709.66±465.68和1 878.20±291.72上升到过渡区的2 770.83±299.11和1 881.94±187.01, 下降到主库区的2 602.06±252.98和1 753.35±197.72, 说明过渡区菌群丰富度最高, 主库区最低.表征菌群多样性的Shannon指数和Simpson指数均呈现先上升后下降的趋势, 过渡区采样点指数均为最高, 主库区最低.说明岗南水库沉积物过渡区好氧反硝化菌群多样性最高, 主库区最低.Goods coverage指数均维持在0.98以上.表征好氧反硝化菌群的均匀性的Pielou指数, 在过渡区最高, 入库区最低.水库的不同空间位置, 沉积物微生物丰富度和多样性也不同.

图 3 岗南水库沉积物不同空间好氧反硝化细菌群生物多样性指数 Fig. 3 Microbial community diversity of aerobic denitrification bacteria in different areas in sediments of Gangnan Reservoir

2.2 微生物物种组成分析

通过物种注释, 水库水体4个季节24个样品共有9 582个OTU, 主要的门类有10个, 分别为变形菌门(Proteobacteria)、异常球菌-栖热菌门(Deinococcus-Thermus)、脊索动物门(Chordata)、放线菌(Actinobacteria)、厚壁菌门(Firmicutes)、硝化螺旋菌门(Nitrospirae)、浮霉菌门(Planctomycetes)、拟杆菌门(Bacteroidetes)、蓝细菌门(Cyanobacteria)和疣微菌门(Verrucomicrobia).变形菌门(Proteobacteria)和异常球菌-栖热菌门(Deinococcus-Thermus)在不同季节的24个采样点都有出现.其中变形菌门(Proteobacteria)均为优势菌门, 所占比例为86.29% ~99.57%(图 4), 与黄小兰等[22]研究的鄱阳湖河交错带水体细菌群落和王欢等[23]、冯亮等[24]研究的碧流河水库优势菌门的结果一致.从空间分布上看, 变形菌门含量从入库区(93.73%)增长到过渡区(94.75%)再增长到主库区(97.66%), 整体呈递增的趋势.从季节分布上来看, 夏季占比最高为98.12%, 春季最低占94.82%.在纲水平, 主要以β-Proteobacteria、γ-Proteobacteria、α-Proteobacteria和δ-Proteobacteria为主. β-Proteobacteria是最大的纲, 所占比例为39.42% ~51.75%, 其中在入库区的点位, 分别在秋季达到最高值, 在冬季达最低值.从入库区到主库区含量变化呈现先降低后增高的趋势, 整体上含量减少, 入库区含量最高, 过渡区最低.在季节分布上β-Proteobacteria春季含量最高, 冬季最低.第二优势纲为γ-Proteobacteria, 所占比例为7.87% ~20.89%. γ-Proteobacteria从入库区到主库区含量变化呈现先增高后降低的趋势, 整体上水流进入主库区含量降低; 在季节分布上, 春季含量最高, 夏季最低.第三优势纲为α-Proteobacteria, 从入库区到主库区含量呈下降趋势, 水流进入主库区含量呈现最低值; 在季节分布上, 秋季含量最高, 夏季含量最低. 不同季节不同空间微生物的组成不同, 这与环境因素密切相关[25].

图 4 好氧反硝化菌群门、纲水平群落结构及分布 Fig. 4 Community structure and distribution of aerobic denitrifying bacteria at the phylum and class level

在属水平(图 5), 陶厄氏菌属(Thauera)为第一优势属, 所占比例为7.02% ~18.84%, 水流在进入水库主库区的过程中, 含量逐渐降低; 季节分布上含量变化趋势为: 春季>夏季>冬季>秋季.第二优势属为贪铜菌属(Cupriavidus), 空间上, 从入库区到主库区, 含量呈递减的趋势, 入库区的1号点位在秋季含量偏高, 占23.65%; 季节分布上的趋势为: 秋季>夏季>春季>冬季, 秋季较其他季节含量明显偏高, 占9.41%.第一、二优势属均从属于β-Proteobacteria.第三优势属为Sutterella, 季节分布上含量变化为: 夏季>冬季>春季>秋季.岗南水库沉积物好氧反硝化菌群结构组成在不同季节不同空间菌群含量均存在差异.

图 5 好氧反硝化菌群属水平群落结构及分布 Fig. 5 Community structure and distribution of aerobic denitrification bacteria at the genus level

2.3 好氧反硝化菌群组成差异性分析

图 6显示同一空间位置不同季节采样点的OTU种类的差异性.入库区所有采样点中, 4个季节共有的OTU种类为482个, 春季、夏季、秋季和冬季独有的OTU种类为544、700、460和812个.过渡区所有采样点中, 4个季节共有的OTU种类为841个, 春季、夏季、秋季和冬季独有的OTU种类为709、621、1 211和817个.主库区所有采样点中, 4个季节共有的OTU种类为1 162个, 春季、夏季、秋季和冬季独有的OTU种类为683、740、1 202和707个.同一空间位置不同季节采样点共有OTU个数从入库区到过渡区再到主库区逐渐增多, 主库区最多.过渡区4个季节独有的OTU个数偏高.秋季独有的OTU种类在过渡区和主库区采样点均为最高, 在入库区为最低.这与多样性分析结果主库区多样性在秋季最高相一致.夏季和秋季独有的OTU个数在过渡区最高, 在入库区最低.说明夏秋季热分层现象对过渡区好氧反硝化菌群的组成影响较大.综上, 说明岗南水库沉积物好氧反硝化菌群同一位置不同季节分布存在显著差异.

图 6 岗南水库沉积物在同一位置不同季节OTU分布的Venn分析 Fig. 6 Venn analysis in sediment OTU distribution in the same location in different seasons in Gangnan Reservoir

图 7显示同一季节不同空间位置采样点的OTU种类的差异性, 春季所有采样点中, 3个位置共有的OTU个数为853个, 入库区、过渡区和主库区独有的OTU个数为494、1 002和1 488个.夏季3个位置共有的OTU个数为884个, 入库区、过渡区和主库区独有的OTU种类为606、846和1 480个; 秋季3个位置共有的OTU个数为776个, 入库区、过渡区和主库区独有的OTU种类为414、1 116和1 694个.冬季3个位置共有的OTU个数为944个, 入库区、过渡区和主库区独有的OTU种类为648、1 010和1 416个.每个季节主库区采样点独有的OTU种类个数均为最多, 入库区均为最少.每个季节不同分区共有的OTU数目, 在冬季最多, 秋季最少.综上, 说明岗南水库沉积物好氧反硝化菌群在同一季节不同空间分布存在显著差异.

图 7 岗南水库沉积物在同一季节不同位置OTU分布的Venn分析 Fig. 7 Venn analysis in sediment OTU distribution at different locations in the same season in Gangnan Reservoir

用LEfSe分析岗南水库沉积物中不同分组间好氧反硝化菌群差异性, 研究与物种层级关系相一致的差异物种.其中图 8显示样本群落中从门到属(从内圈到外圈)主要分类单元的分类等级关系.表 1从空间角度显示了岗南水库沉积物好氧反硝化细菌的LEfSe分析LDA得分值, LDA值越高说明该分类单元的差异越显著, 即该菌群对组间差异性贡献越大.

LDA临界值为2 图 8 岗南水库沉积物好氧反硝化细菌的LEfSe分析 Fig. 8 LEfSe analysis of aerobic denitrifying bacteria in sediments of Gangnan Reservoir

表 1 岗南水库沉积物好氧反硝化细菌的LEfSe分析LDA得分值 Table 1 LDA score value of LEfSe analysis of aerobic denitrifying bacteria in sediments of Gangnan Reservoir

在空间分组中, 入库区采样点分组中差异物种在纲水平的菌群为ε-Proteobacteria; 在目水平的菌群有Alteromonadales、Oceanospirillales和Campylo-bacterales; 在科水平的菌群有Chromatiaceae、Camp-ylobacteraceae、Ferrimonadaceae、Thiobacillaceae、Sphingomonadaceae和Rhizobiaceae; 在属水平的菌群有OceanithermusPhotobacteriumThiocystisSulfurospirillumGeobacterCampylobacterSulfuritortusNovosphingobiumAgrobacteriumFerrimonas.过渡区采样点分组差异物种为γ-Proteobacteria、Gallionellaceae和Sulfuricella. 主库区采样点分组差异物种为门水平的差异性菌群为Proteobacteria; 目水平的菌群有Nitrosomonadales; 科水平的菌群有Budviciaceae、Comamonadaceae、Rhodocyclaceae、Sutterellaceae和Sterolibacteriaceae; 属水平的菌群有LysobacterSerratiaYangiaComamonasAzospiraSutterella.所有差异物种对组间差异性贡献最大的是主库区中的Proteobacteria, 贡献最低的为入库区的Oceanithermus.其中入库区对组间贡献最大的菌群为Alteromonadales、最小的为Oceanithermus.过渡区中对组间贡献最大的菌群为γ-Proteobacteria, 最低的为Sulfuricella.且γ-Proteobacteria在过渡区对数得分值最高.主库区对组间贡献最大的菌群为Proteobacteria, 最小的为Lysobacter.

在季节分组中, 差异物种出现在夏季和秋季.在夏季分组中在科水平的差异物种为Pectobacteriaceae, 在属水平的差异物种为Pectobacterium.在秋季分组中在科水平的差异物种为Yersiniaceae、Burkholderiaceae, 在属水平的差异物种为CupriavidusSalmonellaRhodopseudomonasTardiphaga.所有差异物种对组间差异性贡献最大的为秋季分组中Burkholderiaceae和Cupriavidus, 贡献最低的为秋季分组Tardiphaga.其中贡献最大的科为Burkholderiaceae, 最低的科为Yersiniaceae.贡献最大的属为Cupriavidus, 最低的属为Tardiphaga.从空间分布的水平的LEfSe分析的结果比在季节分布水平的差异物种多, 不同季节、空间分析的差异物种不同.

用随机森林分析寻找不同分组下的标志物种, 从目、科和属这3个分类水平列举了重要性排名靠前的菌群.图 9显示在属水平标志物种的重要性排名.从上到下, 物种对分类模型的重要性依次递减, 其中重要性靠前的物种是组间差异的标志物种.

柱状图横坐标表示物种对分类模型的重要性打分值, 热图显示了物种在各分组的丰度分布 图 9 岗南水库沉积物好氧反硝化细菌的随机森林分析(属水平) Fig. 9 Random forest analysis of aerobic denitrifying bacteria in sediments of Gangnan Reservoir

在空间分组中, 在目水平, Cellvibrionales贡献程度最高, 在过渡区含量较高.在科水平Chromatiaceae贡献最高, 在入库区含量最多.在属水平, Comamonas贡献最高, 在主库区含量较多.其中ThiocystisSulfuricellaAgrobacteriumSaccharo-phagusOceanithermusCampylobacterGeobacterSulfuritortusAzospirillumSalmonellaLaribacterCupriavidusPhotobacteriumFerrimonas在入库区和过渡区含量较高, 主库区含量低. Comamonas在所有空间分组采样点贡献程度最高.

在季节分组中, 在目水平Enterobacterales贡献程度最大, 在春夏季含量较高.在科水平中, Shewanellaceae、Pectobacteriaceae、Yersiniaceae和Methylobacteriaceae贡献程度较高, 其中Shewanellaceae贡献程度最大, 在春季丰度较高. Pectobacteriaceae在春夏季丰度较高, Methylobacteriaceae和Yersiniaceae在秋季含量较高.在属水平PectobacteriumShewanellaSinorhizobiumSalmonellEdwardsiellaBurkholderiaPluralibacter贡献程度较高, 即差异解释率较大.其中Pectobacterium贡献程度最高, 其在春季和夏季丰度均较高.在空间分组中和季节分组中, 丰度值的变化均比较显著.结果表明随机森林分析的标志物种与LEfSe分析结果的差异物种有交叉相同的物种, 但差异贡献程度最高物种不同.

2.4 微生物网络分析

选取物种(OTU)间的斯皮尔曼相关系数r >0.9, 显著性P < 0.001的物种来构建微生物OTU-OTU网络, 研究岗南水库沉积物好氧反硝化菌物种间的关系, 网络分析结果见表 2, 图 10(a)表示物种间网络模块结构, 图 10(b)表示有关节点的物种分类, 图中节点大小与网络节点的度成比例, 节点与节点间边的颜色显示了物种间的作用关系: 红色为正向关系, 绿色为负向关系, 边的宽度与相关系数r成比例.微生物OTU-OTU网络结果表明显示, 正向边占比大于负向边, 说明岗南水库沉积物好氧反硝化菌群物种间合作关系大于竞争关系[16].模块化值高于0.4认为网络具有模块结构[26], 如图 10(a)划分为9个模块, 即有9个相似的生态位[26].其中模块1占比13.47%, 模块2占比11.78%, 模块3占比11.11%, 模块4占比10.10%, 模块5占比8.75%, 模块6占比6.06%, 模块7占比3.03%, 模块8占比2.69%, 模块9占比33.01%.图 10(b)中未分类的菌群占42.76%、β-Proteobacteria占32.66%、γ-Proteobacteria占13.13%、α-Proteobacteria占10.77%、ε-Proteobacteria占0.34%和δ-Proteobacteria占0.34%.

表 2 微生物OTO-OTU网络分析结果 Table 2 Microbial OTO-OTU network analysis results

图 10 岗南水库沉积物好氧反硝化菌群微生物网络 Fig. 10 Network of aerobic denitrification bacteria in sediments of Gangnan Reservoir

由节点数据可知占比较高的模块中模块9包括未分类菌群(44.90%)、β-Proteobacteria(29.60%)、α-Proteobacteria(15.30%)和γ-Proteobacteria(10.20%); 模块1包括未分类菌群(47.50%)、β-Proteobacteria(27.50%)、α-proteobacteria(12.50%)和γ-Proteobacteria(12.50%).

基于节点度大于5和中介中心值小于1 000的标准[11], 得到47个关键物种OTU, 分布在模块1中的有7个OTU, 分布在模块2中的有13个OTU, 分布在模块3中的有7个OTU, 分布在模块4中的有7个OTU, 分布在模块5中的有6个OTU, 分布在模块6中的有6个OTU, 分布在模块9中的有1个OTU.关键物种OTU如下: 9、13、14、18、27、30、31、41、43、59、82、91、116、131、167、200、209、215、238、306、310、315、335、351、355、427、459、460、509、518、727、732、749、750、764、864、1 012、1 035、1 073、1 141、1 366、1 512、1 564、1 612、2 331、2 473和3 015.其中分为未分类菌群(44.68%)、β-Proteobacteria(36.17%)、γ-Proteobacteria(12.77%)和α-Proteobacteria(6.38%), 主要有13个属, 分别为ThaueraAzoarcusAzospiraAzospirillumCitrobacterComamonasPseudomonasShewanellaSinorhizobiumSulfuricellaSulfuritaleaThioalkalivibrio和unclassified.除去未分类菌群, Thauera含量最多.

应用PCoA分析对关键物种OTU进一步分析, 根据图 11可知岗南水库沉积物好氧反硝化菌群中47个关键物种在时空分布存在一定差异性, PCoA1和PCoA2分别解释了整体46.62%和9.60%, 总共解释整体的56.22%, PCoA1占主体.不同季节采样点基于Adonis分析(R2=0.120 19, P>0.05)显示季节差异性不明显.不同空间采样点基于Adonis分析(R2=0.559 22, P < 0.001)显示空间分布差异性显著.其中入库区采样点主要分布在第一、三象限, 其中春季采样点主要分布在第一象限, 夏秋冬季主要分布在第三象限.过渡区的采样点主要分布在第一、二和三象限, 其中秋季主要分布在第二象限, 春季分布在第一象限, 夏冬季主要分布在第三象限.主库区采样点主要分布在第二四象限, 其中秋季主要分布在第二象限, 春夏冬季主要分布在第四象限.入库区采样点与过渡区采样点较聚集, 说明关键物种在入库区和过渡区菌群结构有相似性, 与主库区的菌群组成差异性明显.春夏冬季的采样点较聚集, 说明春夏冬季菌群结构相似, 秋季采样点较分散表明秋季采样点菌群结构差异性明显.说明不同空间和不同时间对好氧反硝化菌群中的关键物种有一定的选择性[27].

图 11 关键物种不同时空主成分分析 Fig. 11 Principal component analysis of key species in different time and space

3 讨论 3.1 好氧反硝化菌群多样性的时空分布特征

从岗南水库入库河口到坝前主库区沉积物中好氧反硝化菌群丰富度和多样性呈现先上升后下降的趋势, 过渡区菌群丰富度和多样性最高.黄小兰等[22]的研究发现鄱阳湖河交汇处细菌比湖区内丰富度高.刘晓丹等[28]的研究发现尼洋河沉积物上游、中游和下游不同河段沉积物微生物差异性明显.陈星等[29]的研究发现三峡水库库中生物多样性指数高于库首和库尾.说明微生物的多样性和丰富度与湖库所处的地理环境和空间位置都有关系, 岗南水库不同的空间位置沉积物好氧反硝化菌群的丰富度和多样性均不同.在季节分布上, 好氧反硝化菌群丰富度和多样性在秋季最高, 其中秋季水库存在热分层现象, 上下层水体物质交换受阻, 会进一步影响到沉积物中微生物的分布及多样性[30].后续可深入探究热分层条件下好氧反硝化菌群垂向分布特征, 为热分层时期水库水质改善提供科学依据.谭八梅等[31]的研究指出海刺参养殖池塘水体菌群秋季丰富度和多样性均为最低, 与本研究结果相反.王欢等[23]的研究指出全年碧流河水库细菌多样性秋季最高, 与本研究的结果一致.菌群的多样性与所处的环境密切相关[31], 所以研究结果各不相同.本文单独对主库区的采样点进行研究, 发现主库区采样点水样好氧反硝化菌群演替规律和水库整体演替规律有重合点, 都是秋季多样性最高, 夏季最低.这与秦宇等[32]在三峡库区中段水体中微生物多样性夏季较低的研究结果一致.说明秋季的水力环境可能更有利于好氧反硝化菌群的繁殖[33].综上, 不同季节和空间均会对好氧反硝化菌群的多样性产生影响.

3.2 好氧反硝化菌群组成的时空分布特征

岗南水库的好氧反硝化菌群主要以变形菌门为主, 与白洋淀好氧反硝化菌群组成[34]和密云水库入库河流水体反硝化菌群[35]的研究结果一致, 说明可能淡水湖泊中主要好氧反硝化菌群均为变形菌门[24], 且变形菌门可能与沉积物中含氮有机物等营养物质有紧密关系[36~38]. 变形菌门在氮素消除和有机物的降解中发挥重要作用, 且大部分在厌氧条件下可以进行反硝化[39~42], 是水生态系统污染恢复的重要参与者.Romírez等[43]的研究发现墨西哥湾西北部的表层沉积物中变形菌门受深度、距海岸线的距离和温度溶解氧等影响, 且不同细菌在不同位置占优势, δ-Proteobacteria在内陆架、γ-Proteobacteria在上坡沉积物占优势, 菌群分布受环境影响显著.Tang等[44]的研究发现温度和变形菌门的协同作用可能是厌氧/缺氧/好氧(A2O)污水处理中出水有机物季节变化的机制. Ke等[45]的研究发现长期风化的铝土矿渣处置区中pH和电导率的改变对变形菌门的多样性丰度均有影响, 人类可以干预处置区进行生态重建. 已有研究均表明沉积物中好氧反硝化菌组成结构是与环境相适应的结果.

从空间分布上看, 变形菌门(优势菌门)在主库区含量最高, 从季节分布上来看, 夏季占比最高.与多样性的分布特征有明显差异.在纲水平, 主要以β-Proteobacteria、γ-Proteobacteria和α-Proteobacteria为主.研究表明变形菌纲包含大多数反硝化菌群[46, 47]. β-Proteobacteria是最大的纲, 入库区含量最高.在季节分布上, 春季含量最高.第二和三优势纲为γ-Proteobacteria和α-Proteobacteria, 3种优势纲中存在可以进行异养或自养反硝化的反硝化菌群[48], 将NO3-最终转化为N2, 在氮循环水质净化过程中扮演着很重要的角色.岗南水库沉积物中存在的主要纲类与国内外研究沉积物反硝化菌群结果一致[49~51].优势属为陶厄氏菌属(Thauera), 其反硝化活性较高[52], 能有效地去除氮素[53, 54], 对有机物有较强的降解能力.第二优势属贪铜菌属(Cupriavidus), 空间上入库区含量最高, 季节分布上秋季较其他季节含量明显偏高.第一、二优势属均从属于β-Proteobacteria.第三优势属为Sutterella, 空间上含量变化与贪铜菌属(Cupriavidus)一致, 季节分布上含量变化为夏季含量最高. Aeromonas(气单胞菌属)能引起鱼类得急性传染病, 对水产养殖业有潜在风险[55].贪铜菌属(Cupriavidus)有异养硝化和好氧反硝化的能力, 可高效利用硝态氮和亚硝态氮, 有一定的耐金属特性, 在处理焦化废水中发挥着重要作用[56].其中Pseudomonas(假单胞菌属)春季含量偏高, 有研究表明, 其对硝酸盐氮有较好的去除效果, 是比较大的反硝化细菌的类群[47]. Aeromonas(气单胞菌属)还被鉴定为具有理想性能的关键反硝化聚磷菌[57].岗南水库沉积物好氧反硝化菌群结构组成在不同季节不同空间菌群含量均存在差异, 与OTU水平的韦恩分析相互印证.

3.3 差异物种与关键物种分析

不同季节和不同空间基于LEfSe分析的差异物种不同.随机森林分析得到的标志物种与LEfSe分析的差异物种相互印证, 除差异解释率最大的物种外大部分差异物种和标志物种相同. 本研究结果发现AzospiraAzospirillumComamonasShewanellaSinorhizobiumSulfuricella既是差异物种和标志物种又是关键物种.因此这些菌群在岗南水库好氧反硝化菌群中占重要地位, 在今后应加强对这些菌群的研究.其中Azospira sp. KJ可还原高氯酸盐, 有高效脱氯的效果[58]. Azospirillum lipoferum AL-3可以帮助马铃薯植物对抗早疫病, 从而降低环境污染[59]. Sinorhizobium的存在可以提高畜牧饲料的产量和有效地低温固氮[60]. Comamonas可以降解多环芳烃和杂环芳烃等污染物[61], 对邻-甲醛和对-甲酚也有降解特性[62]. 不同的差异物种、标志物种和关键物种对不同的环境因素适应能力不同, 比如一些菌群与光照和温度相关. Saejung等[63]的研究发现, 不同光照条件可以提高Rhodopseudomonas faecalis PA2的生物量和类胡萝卜素等有价值物质产生. Rhodopseudomonas有较高温度的适应性, 生物产氢速率也得到了提高[64].Zhang等[65]的研究发现Pseudomonas sp. Y39-6在低温条件可去除废水中的硝酸盐. Adyasari等[66]的研究发现海底地下水中的Burkholderiaceae主要存在于低盐度样品中.不同环境因素会造成差异物种、标志物种和关键物种的不同. 综上, 不同空间、不同季节和各种环境因素都会引起水库沉积物好氧反硝化菌群组成的改变, 进而对水库水质演变产生影响.基于差异物种和关键物种的解析, 可以为将来研究水库微生物驱动的氮素循环提供必要的数据支持.

4 结论

(1) α-多样性分析显示从季节分布上看, 秋冬季菌群丰富度和多样性均较高, 其中秋季好氧反硝化菌群多样性最高; 从空间上看, 过渡区的好氧反硝化菌群丰富度和多样性均为最高.岗南水库好氧反硝化细菌的优势菌门为变形菌门(Proteobacteria), 所占比例为86.29% ~99.57%.在纲水平, 主要以β-Proteobacteria、γ-Proteobacteria和α-Proteobacteria为主.第一优势属为陶厄氏菌属(Thauera).

(2) Venn图显示岗南水库沉积物好氧反硝化菌群时空分布存在显著差异.LEfSe分析显示在季节分组中, 差异物种出现在夏季和秋季.所有差异物种对组间差异性贡献最大的为秋季分组中Burkholderiaceae和Cupriavidus, 贡献最低的为秋季分组的Tardiphaga.在空间分组中, 所有差异物种对组间差异性贡献最大的为主库区的Proteobacteria, 贡献最低的为入库区的Oceanithermus.随机森林分析显示在季节分组中, Pectobacterium贡献程度最高, 差异解释率最大, 其在春季和夏季丰度均较高.在空间分组中, Comamonas贡献程度最高, 在主库区丰度较高.

(3) 网络分析显示好氧反硝化菌群物种间合作关系大于竞争关系.得到47个关键物种OTU, 已分类的菌群中Thauera含量最多.应用主成分分析对关键物种进一步分析表明, 关键物种中秋季采样点菌群组成与其他季节差异性明显, 而春季、夏季和冬季菌群组成有一定的相似性.主库区关键物种采样点菌群组成与入库区和过渡区差异性明显.

参考文献
[1] 黄廷林. 水源水库水质污染原位控制与改善是饮用水水质安全保障的首要前提[J]. 给水排水, 2017, 43(1): 1-3, 69. DOI:10.3969/j.issn.1002-8471.2017.01.001
[2] 文新宇, 谌书, 杨远坤, 等. 缓丘地带黄河镇水库微生物群落多样性特征[J]. 应用与环境生物学报, 2021, 27(6): 1-10.
Wen X Y, Chen S, Yang Y K, et al. Characteristics of microbial community diversity in in the gentle hilly area-Huanghe Town Reservoir[J]. Chinese Journal of Applied & Environmental Biology, 2021, 27(6): 1-10.
[3] 原晖. 大连市大型水库富营养化状况及影响因素研究[J]. 水资源开发与管理, 2021(2): 12-16.
Yuan H. Study on eutrophication status and influencing factors of large reservoirs in Dalian[J]. The Global Seabuckthorn Research and Development, 2021(2): 12-16.
[4] 黄廷林, 曾明正, 邱晓鹏, 等. 温带季节性分层水库浮游植物功能类群的时空演替[J]. 中国环境科学, 2016, 36(4): 1157-1166.
Huang T L, Zeng M Z, Qiu X P, et al. Phytoplankton functional groups and their spatial and temporal distribution characteristics in a temperate seasonally stratified reservoir[J]. China Environmental Science, 2016, 36(4): 1157-1166. DOI:10.3969/j.issn.1000-6923.2016.04.029
[5] 王顺天, 雷俊山, 贾海燕, 等. 三峡水库2003~2017年水质变化特征及成因分析[J]. 人民长江, 2020, 51(10): 47-53, 127.
Wang S T, Lei J S, Jia H Y, et al. Water quality characteristics and causes analysis of Three Gorges Reservoir from 2003 to 2017[J]. Yangtze River, 2020, 51(10): 47-53, 127.
[6] 康丽娟, 许海, 朱广伟, 等. 太湖主要环湖河道沉积物反硝化潜力及其控制因子[J]. 环境科学学报, 2021, 41(4): 1393-1400.
Kang L J, Xu H, Zhu G W, et al. Sediment denitrification potential and its influencing factors in the main rivers of Lake Taihu[J]. Acta Scientiae Circumstantiae, 2021, 41(4): 1393-1400.
[7] Gardner W S, Yang L Y, Cotner J B, et al. Nitrogen dynamics in sandy freshwater sediments (Saginaw Bay, Lake Huron)[J]. Journal of Great Lakes Research, 2001, 27(1): 84-97. DOI:10.1016/S0380-1330(01)70624-7
[8] Gross A, Boyd C E, Wood C W. Nitrogen transformations and balance in channel catfish ponds[J]. Aquacultural Engineering, 2000, 24(1): 1-14. DOI:10.1016/S0144-8609(00)00062-5
[9] 屈建航, 李宝珍, 袁红莉. 沉积物中微生物资源的研究方法及其进展[J]. 生态学报, 2007, 27(6): 2636-2641.
Qu J H, Li B Z, Yuan H L. Methodology in studies of microbial resource in sediment[J]. Acta Ecologica Sinica, 2007, 27(6): 2636-2641. DOI:10.3321/j.issn:1000-0933.2007.06.058
[10] Robertson L A, Kuenen J G. Thiosphaera pantotropha gen. nov. sp. nov., a facultatively anaerobic, facultatively autotrophic sulphur bacterium[J]. Microbiology, 1983, 129(9): 2847-2855. DOI:10.1099/00221287-129-9-2847
[11] 汤茵琪, 李阳, 常素云, 等. 好氧反硝化菌强化生态浮床对水体氮与有机物净化机理[J]. 生态学杂志, 2017, 36(2): 569-576.
Tang Y Q, Li Y, Chang S Y, et al. Mechanisms of nitrogen and organic matter removal from water bodies by aerobic denitrifiers enhanced ecological floating bed[J]. Chinese Journal of Ecology, 2017, 36(2): 569-576.
[12] 李炳堂, 周志勤, Naidu R, 等. 锁磷剂联合好氧反硝化菌修复富营养化水体[J]. 环境科学, 2021, 42(4): 1861-1869.
Li B T, Zhou Z Q, Naidu R, et al. Combined remediation of eutrophic water by phoslock(R) and aerobic denitrifying bacteria[J]. Environmental Science, 2021, 42(4): 1861-1869.
[13] 王慎, 张思思, 许尤, 等. 水库沉积物贫营养型好氧反硝化菌(Pseudomonas sp.)分离及其对微污染水体脱氮效率分析[J]. 环境科学学报, 2020, 40(3): 915-921.
Wang S, Zhang S S, Xu Y, et al. Isolation of aerobic denitrifying bacteria (Pseudomonas sp.) from reservoir sediments and analysis of its denitrification efficiency in micro-polluted water[J]. Acta Scientiae Circumstantiae, 2020, 40(3): 915-921.
[14] 蒋云霞. 水源水库沉积物中好氧反硝化菌群筛选及脱氮特性研究[D]. 西安: 西安建筑科技大学, 2020.
[15] 周石磊, 张艺冉, 孙悦, 等. 异养硝化-好氧反硝化菌富集驯化过程中微生物种群演变特征——典型城市景观水系[J]. 中国环境科学, 2019, 39(11): 4831-4839.
Zhou S L, Zhang Y R, Sun Y, et al. Characteristics of bacterial community structure during the enrichment and domestication of heterotrophic nitrification-aerobic denitrification bacteria based on the typical city landscape water[J]. China Environmental Science, 2019, 39(11): 4831-4839. DOI:10.3969/j.issn.1000-6923.2019.11.043
[16] 周石磊, 孙悦, 岳哿丞, 等. 雄安新区-白洋淀冬季冰封期水体好氧反硝化菌群落空间分布特征及驱动因素[J]. 环境科学, 2020, 41(5): 2177-2187.
Zhou S L, Sun Y, Yue G C, et al. Spatial distribution characteristics and driving factors of aerobic denitrification bacterial community structure from Baiyangdian Lake in Xiong'an New Area during the winter freezing period[J]. Environmental Science, 2020, 41(5): 2177-2187.
[17] 寄博华, 李玮, 常军军, 等. 滇池湖滨湿地不同挺水植物区沉积物细菌群落结构特征[J]. 生态与农村环境学报, 2020, 36(3): 390-398.
Ji B H, Li W, Chang J J, et al. Characteristics of bacterial community structure in sediments with different emergent plants in Dianchi Lake-side Wetland[J]. Journal of Ecology and Rural Environment, 2020, 36(3): 390-398.
[18] Li Y, Wu H N, Shen Y, et al. Statistical determination of crucial taxa indicative of pollution gradients in sediments of Lake Taihu, China[J]. Environmental Pollution, 2019, 246: 753-762. DOI:10.1016/j.envpol.2018.12.087
[19] Zhang B G, Zhang J, Liu Y, et al. Co-occurrence patterns of soybean rhizosphere microbiome at a continental scale[J]. Soil Biology and Biochemistry, 2018, 118: 178-186. DOI:10.1016/j.soilbio.2017.12.011
[20] 周石磊, 孙悦, 岳哿丞, 等. 白洋淀养殖区春季沉积物垂向微生物群落特征及驱动因素[J]. 环境科学学报, 2020, 40(5): 1722-1733.
Zhou S L, Sun Y, Yue G C, et al. Vertical distribution characteristics and driving factors of microbial community of spring culture area sediments in Baiyangdian Lake[J]. Acta Scientiae Circumstantiae, 2020, 40(5): 1722-1733.
[21] Anderson M J, Ellingsen K E, Mcardle B H. Multivariate dispersion as a measure of beta diversity[J]. Ecology Letters, 2006, 9(6): 683-693. DOI:10.1111/j.1461-0248.2006.00926.x
[22] 黄小兰, 简正军, 叶长盛, 等. 鄱阳湖河湖交错带枯水期水体细菌群落空间分布[J]. 微生物学通报, 2021, 48(6): 1862-1873.
Huang X L, Jian Z J, Ye C S, et al. Spatial distribution of bacterial community in River-lake ecotone of Poyang Lake during dry season[J]. Microbiology China, 2021, 48(6): 1862-1873.
[23] 王欢, 赵文, 谢在刚, 等. 碧流河水库细菌群落结构特征及其关键驱动因子[J]. 环境科学, 2018, 39(8): 3660-3669.
Wang H, Zhao W, Xie Z G, et al. Bacterial community structure characteristics in the Biliuhe Reservoir and its key driving factors[J]. Environmental Science, 2018, 39(8): 3660-3669.
[24] 冯亮, 袁春燕, 杨超, 等. 好氧反硝化生物脱氮技术的研究进展[J]. 微生物学通报, 2020, 47(10): 3342-3354.
Feng L, Yuan C Y, Yang C, et al. Research progress in nitrogen removal by aerobic denitrification[J]. Microbiology China, 2020, 47(10): 3342-3354.
[25] 朱婷婷, 田从魁. 水库底泥中微生物多样性及其与环境因子相关性分析[J]. 北京大学学报(自然科学版), 2018, 54(3): 625-632.
Zhu T T, Tian C K. Analysis on Microbial diversity in the sediments and its relationship with environmental factors in a reservoir[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2018, 54(3): 625-632.
[26] 沈烽, 黄睿, 曾巾, 等. 分子生态网络分析研究进展[J]. 环境科学与技术, 2016, 39(S1): 94-98, 173.
Shen F, Huang R, Zeng J, et al. Progress of the molecular ecology network analysis[J]. Environmental Science & Technology, 2016, 39(S1): 94-98, 173.
[27] 韩文杰, 吴迪, 周家中, 等. 长三角地区MBBR泥膜复合污水厂低温季节微生物多样性分析[J]. 环境科学, 2020, 41(11): 5037-5049.
Han W J, Wu D, Zhou J Z, et al. Microbial diversity analysis of WWTPs based on Hybrid-MBBR process in a low temperature season in the Yangtze River Delta[J]. Environmental Science, 2020, 41(11): 5037-5049.
[28] 刘晓丹, 黄毅, 王永花, 等. 西藏尼洋河沉积物中微生物群落结构特征分析[J]. 环境科学, 2020, 41(7): 3249-3256.
Liu X D, Huang Y, Wang Y H, et al. Structural characteristics of microbial communities in the sediments of the Niyang River in Tibet[J]. Environmental Science, 2020, 41(7): 3249-3256.
[29] 陈星, 叶麟, 谭路, 等. 三峡水库蓄水后期香溪河库湾浮游甲壳动物群落结构[J]. 水生生物学报, 2020, 44(5): 1070-1079.
Chen X, Ye L, Tan L, et al. Crustacean zooplankton community in Xiangxi bay during the later impoundment period of The Three Gorges Reservoir[J]. Acta Hydrobiologica Sinica, 2020, 44(5): 1070-1079.
[30] 王丽娜, 陈辉煌, 刘乐冕, 等. 亚热带分层水库固氮微生物时空分布格局[J]. 生态学报, 2016, 36(18): 5827-5837.
Wang L N, Chen H H, Liu L M, et al. Spatiotemporal patterns of diazotrophic microorganisms in a subtropical stratified reservoir[J]. Acta Ecologica Sinica, 2016, 36(18): 5827-5837.
[31] 谭八梅, 王荦, 裴泓霖, 等. 不同季节刺参养殖池塘水体菌群结构与功能特征研究[J]. 渔业科学进展, 2021, 42(3): 77-88.
Tan B M, Wang L, Pei H L, et al. The seasonal structural and functional characteristics of bacterial community in Apostichopus japonicus culture pond water[J]. Progress in Fishery Sciences, 2021, 42(3): 77-88.
[32] 秦宇, 张曦, 王紫薇, 等. 三峡水库中段水体微生物群落结构季节变化[J/OL]. 水生态学杂志, 2020, https://kns.cnki.net/kcms/detail/42.1785.X.20201104.1212.002.html, 2020-11-06.
Qin Y, Zhang X, Wang Z W, et al. Community structure of microorganisms and seasonal variation characteristics in the middle of the Three Gorges Reservoir[J/OL]. Journal of Hydroecology, 2020, https://kns.cnki.net/kcms/detail/42.1785.X.20201104.1212.002.html, 2020-11-06.
[33] 圣平, 于一尊, 王东升, 等. 通过PCR-DGGE技术分析鄱阳湖不同河口底泥中微生物多样性[J]. 江西科学, 2015, 33(3): 312-317, 329.
Sheng P, Yu Y Z, Wang D S, et al. Analysis of bacterial diversities in different estuarine sediments of Poyang Lake by PCR-DGGE[J]. Jiangxi Science, 2015, 33(3): 312-317, 329.
[34] 周石磊, 张甜娜, 陈召莹, 等. 输水情景下白洋淀好氧反硝化菌群落对溶解性有机物的响应[J]. 环境科学, 2021, 42(9): 4319-4331.
Zhou S L, Zhang T N, Chen Z Y, et al. Response of aerobic denitrification bacterial community structure to dissolved organic matter in Baiyangdian lake during the water delivery period[J]. Environmental Science, 2021, 42(9): 4319-4331.
[35] 李垒, 孟庆义, 叶飞, 等. 密云水库入库河流水体总细菌和反硝化菌群落组成与结构[J]. 水生态学杂志, 2018, 39(6): 44-51.
Li L, Meng Q Y, Ye F, et al. Microbial community structures of total bacteria and denitrifying bacteria in the tributaries of Miyun reservoir[J]. Journal of Hydroecology, 2018, 39(6): 44-51.
[36] 刘鹄, 孟婷, 王敏, 等. 水源地水库沉积物中营养盐与微生物群落的空间分布特征[J]. 水土保持通报, 2016, 36(3): 24-29, 37.
Liu H, Meng T, Wang M, et al. Spatial distribution of nutrients and bacterial communities in sediments of water resource reservoir[J]. Bulletin of Soil and Water Conservation, 2016, 36(3): 24-29, 37.
[37] 邱莹, 靳燕, 苏振华, 等. 城市再生水河道沉积物细菌群落空间变化分析: 以京津冀北运河为例[J]. 环境科学, 2021, 42(5): 2287-2295.
Qiu Y, Jin Y, Su Z H, et al. Analysis of the spatial changes in bacterial communities in urban reclaimed water channel sediments: a case study of the North Canal River[J]. Environmental Science, 2021, 42(5): 2287-2295.
[38] 于妍, 王悦悦, 方社贤, 等. 白洋淀表层沉积物细菌多样性及影响因素[J]. 环境工程学报, 2021, 15(3): 1121-1130.
Yu Y, Wang Y Y, Fang S X, et al. Bacterial diversity in surface sediments of Baiyangdian Lake and its influencing factors[J]. Chinese Journal of Environmental Engineering, 2021, 15(3): 1121-1130.
[39] McLellan S L, Newton R J, Vandewalle J L, et al. Sewage reflects the distribution of human faecal Lachnospiraceae[J]. Environmental Microbiology, 2013, 15(8): 2213-2227. DOI:10.1111/1462-2920.12092
[40] Wu W Z, Yang L H, Wang J L. Denitrification performance and microbial diversity in a packed-bed bioreactor using PCL as carbon source and biofilm carrier[J]. Applied Microbiology and Biotechnology, 2013, 97(6): 2725-2733. DOI:10.1007/s00253-012-4110-4
[41] Hou L F, Zhou Q, Wu Q P, et al. Spatiotemporal changes in bacterial community and microbial activity in a full-scale drinking water treatment plant[J]. Science of the Total Environment, 2018, 625: 449-459. DOI:10.1016/j.scitotenv.2017.12.301
[42] 陈倩茹, 夏雪, 王川, 等. 沙颍河下游城市黑臭内河沉积物微生物群落季节变化特征[J]. 水生生物学报, 2021, 45(1): 182-189.
Chen Q R, Xia X, Wang C, et al. Seasonal variation characteristics of microbial communities in the sediments of urban black-odorous river in lower reaches of Shaying River[J]. Acta Hydrobiologica Sinica, 2021, 45(1): 182-189.
[43] Ramírez D, Vega-Alvarado L, Taboada B, et al. Bacterial diversity in surface sediments from the continental shelf and slope of the North West gulf of Mexico and the presence of hydrocarbon degrading bacteria[J]. Marine Pollution Bulletin, 2020, 150. DOI:10.1016/j.marpolbul.2019.110590
[44] Tang G, Zheng X, Li X L, et al. Variation of effluent organic matter (EfOM) during anaerobic/anoxic/oxic (A2O) wastewater treatment processes[J]. Water Research, 2020, 178. DOI:10.1016/j.watres.2020.115830
[45] Ke W S, Zhang X C, Zhu F, et al. Appropriate human intervention stimulates the development of microbial communities and soil formation at a long-term weathered bauxite residue disposal area[J]. Journal of Hazardous Materials, 2021, 405. DOI:10.1016/j.jhazmat.2020.124689
[46] 唐婧, 徐小蓉, 商传禹, 等. 南明河城区河段细菌多样性与环境因子的关系[J]. 微生物学报, 2015, 55(8): 1050-1059.
Tang J, Xu X R, Shang C Y, et al. Association of bacterial diversity in city area of Nanming River with environmental factors[J]. Acta Microbiologica Sinica, 2015, 55(8): 1050-1059.
[47] 杨文焕, 石大钧, 张元, 等. 高原湖泊沉积物中反硝化微生物的群落特征——以包头南海湖为例[J]. 中国环境科学, 2020, 40(1): 431-438.
Yang W H, Shi D J, Zhang Y, et al. Community characteristics of denitrifying microorganisms in plateau lake sediments[J]. China Environmental Science, 2020, 40(1): 431-438. DOI:10.3969/j.issn.1000-6923.2020.01.048
[48] Schreier H J, Mirzoyan N, Saito K. Microbial diversity of biological filters in recirculating aquaculture systems[J]. Current Opinion in Biotechnology, 2010, 21(3): 318-325. DOI:10.1016/j.copbio.2010.03.011
[49] Jiang X L, Yao L, Guo L D, et al. Multi-scale factors affecting composition, diversity, and abundance of sediment denitrifying microorganisms in Yangtze lakes[J]. Applied Microbiology and Biotechnology, 2017, 101(21): 8015-8027. DOI:10.1007/s00253-017-8537-5
[50] 康鹏亮, 黄廷林, 张海涵, 等. 西安市典型景观水体水质及反硝化细菌种群结构[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.
[51] 吴玲. 富营养化太湖沉积物中微生物群落及对环境因子的响应[J]. 应用与环境生物学报, 2019, 25(6): 1470-1476.
Wu L. Response of microbial community structure in eutrophic Taihu sediments to environmental factors[J]. Chinese Journal of Applied & Environmental Biology, 2019, 25(6): 1470-1476.
[52] 徐富强, 桂梦瑶, 杜俊逸, 等. 典型工业污染物对好氧反硝化菌群脱氮性能及群落结构的影响[J]. 环境工程学报, 2019, 13(10): 2442-2450.
Xu F Q, Gui M Y, Du J Y, et al. Effects of typical industrial pollutants on the denitrification performance and the community structure of aerobic denitrifying bacteria[J]. Chinese Journal of Environmental Engineering, 2019, 13(10): 2442-2450. DOI:10.12030/j.cjee.201812007
[53] 鲍林林, 王晓燕, 陈永娟, 等. 北运河沉积物中主要脱氮功能微生物的群落特征[J]. 中国环境科学, 2016, 36(5): 1520-1529.
Bao L L, Wang X Y, Chen Y J, et al. Diversity, abundance and distribution of nirS-type denitrifiers and Anammox bacteria in sediments of Beiyun River[J]. China Environmental Science, 2016, 36(5): 1520-1529. DOI:10.3969/j.issn.1000-6923.2016.05.035
[54] 杜瑞芳, 李靖宇, 赵吉. 乌梁素海湖滨湿地细菌群落结构多样性[J]. 微生物学报, 2014, 54(10): 1116-1128.
Du R F, Li J Y, Zhao J. Bacterial diversity in littoral wetland of Wuliangsuhai Lake[J]. Acta Microbiologica Sinica, 2014, 54(10): 1116-1128.
[55] 马富平. 主养草鱼池塘微生物群落结构分析及嗜水气单胞菌拮抗菌筛选[D]. 重庆: 西南大学, 2017.
[56] 孙智毅. 一株异养硝化—好氧反硝化Cupriavidus sp. S1的筛选及降解特性研究[D]. 太原: 太原理工大学, 2017.
[57] Zhang X X, Bahtiar F, Wang Y, et al. Achieving simultaneous biological nutrient removal and sludge minimization from marine ship sewage based on an innovative landscape integrated ecological treatment system (LIETS)[J]. Ecological Engineering, 2020, 156. DOI:10.1016/j.ecoleng.2020.105989
[58] Puigserver D, Herrero J, Parker B L, et al. Natural attenuation of pools and plumes of carbon tetrachloride and chloroform in the transition zone to bottom aquitards and the microorganisms involved in their degradation[J]. Science of the Total Environment, 2020, 712. DOI:10.1016/j.scitotenv.2019.135679
[59] Mehmood T, Li G H, Anjum T, et al. Azospirillum lipoferum strain AL-3 reduces early blight disease of potato and enhance yield[J]. Crop Protection, 2021, 139. DOI:10.1016/j.cropro.2020.105349
[60] Robin C, Sultan-Tubeileh K, Obaton M, et al. Nitrogen fixation and growth of annual Medicago-Sinorhizobium associations at low temperature[J]. European Journal of Agronomy, 2005, 22(3): 267-275. DOI:10.1016/j.eja.2004.04.001
[61] 王乔, 郑瑞, 孙学婷, 等. 睾丸酮丛毛单胞菌对羊草根际土壤PAHs降解及细菌群落结构的影响[J]. 生物工程学报, 2020, 36(12): 2657-2673.
Wang Q, Zheng R, Sun X T, et al. Effects of Comamonas testosteroni on PAHs degradation and bacterial community structure in Leymus chinensis rhizosphere soil[J]. Chinese Journal of Biotechnology, 2020, 36(12): 2657-2673.
[62] 张立志, 余思彤, 袁欣, 等. 丛毛单胞菌对邻甲酚及对甲酚的降解特性[J]. 环境污染与防治, 2020, 42(7): 820-825, 832.
Zhang L Z, Yu S T, Yuan X, et al. Degradation characteristics of o-cresol and p-cresol by Comamonas sp.[J]. Environmental Pollution and Control, 2020, 42(7): 820-825, 832.
[63] Saejung C, Chanthakhot T. Single-phase and two-phase cultivations using different light regimes to improve production of valuable substances in the anoxygenic photosynthetic bacterium Rhodopseudomonas faecalis PA2[J]. Bioresource Technology, 2021, 328. DOI:10.1016/j.biortech.2021.124855
[64] Du Toit J P, Pott R W M. Heat-acclimatised strains of Rhodopseudomonas palustris reveal higher temperature optima with concomitantly enhanced biohydrogen production rates[J]. International Journal of Hydrogen Energy, 2021, 46(21): 11564-11572. DOI:10.1016/j.ijhydene.2021.01.068
[65] Zhang D Y, Liu Y, Han Y X, et al. Nitrate removal from low C/N wastewater at low temperature by immobilized Pseudomonas sp. Y39-6 with versatile nitrate metabolism pathways[J]. Bioresource Technology, 2021, 326. DOI:10.1016/j.biortech.2021.124794
[66] Adyasari D, Hassenrück C, Oehler T, et al. Microbial community structure associated with submarine groundwater discharge in northern Java (Indonesia)[J]. Science of the Total Environment, 2019, 689: 590-601. DOI:10.1016/j.scitotenv.2019.06.193