环境科学  2022, Vol. 43 Issue (7): 3635-3644   PDF    
天津独流减河流域不同等级河流沉积物细菌区系及功能辨识
刘嘉元1, 丰玥1, 杨雪纯1, 张彦2, 李德生1, 刘福德1     
1. 天津理工大学环境科学与安全工程学院, 天津 300384;
2. 天津市生态环境科学研究院, 天津 300191
摘要: 细菌群落对环境变化的敏感性较高, 但在相似等级的生态条件下其代谢功能或具有趋同性.为探讨微生物群落的这种环境属性并验证其用于流域功能单元划分方面的可行性, 本文以天津市独流减河流域上游汇水段河流表层沉积物为研究对象, 从河道干流(MS)、左岸支流(LT)和右岸支流(RT)这3个连续的空间单元入手, 研究沉积物中细菌的多样性、结构和功能代谢丰度, 并探讨其作为流域功能单元划分关键因子的可行性.结果表明, 不同等级河流沉积物之间细菌的Shannon指数和Simpson指数无显著差异(P>0.05), 但Chao1、Ace、Observed_species和PD_whole_tree指数均呈现支流显著高于干流的现象(P < 0.05).细菌优势菌门的相对丰度在不同等级河流沉积物中差异不显著, 但其功能代谢丰度则表现为干流显著高于支流(P<0.05), 而不同支流之间差异不显著(P>0.05).本研究发现, 除可交换态磷和硝态氮之外, 其他理化性质在不同河流空间单元之间差异不显著, 但pH和铁铝结合态磷显著影响了细菌群落的结构组成, 有机质和总氮则与细菌功能代谢密切相关.研究同时发现, 基于细菌代谢丰度可将河道干流和支流清晰地分成了Level-1和Level-2两个不同的功能单元, 表明沉积物中细菌群落的功能代谢丰度可以作为流域功能单元划分的依据.
关键词: 独流减河流域      表层沉积物      细菌多样性      代谢功能      功能分区     
Identification of Bacterial Flora and Metabolic Function of Sediments in Different Channels of Duliujian River Basin, Tianjin
LIU Jia-yuan1 , FENG Yue1 , YANG Xue-chun1 , ZHANG Yan2 , LI De-sheng1 , LIU Fu-de1     
1. School of Environmental Science and Safety Engineering, Tianjin University of Technology, Tianjin 300384, China;
2. Tianjin Academy of Eco-Environmental Sciences, Tianjin 300191, China
Abstract: Bacterial communities are highly sensitive to environmental changes, but their metabolic functions may be convergent under similar ecological conditions. In order to test this environmental attribute of the bacterial community and verify the feasibility in using the bacterial metabolic data to divide functional units at the river basin scale, the surface sediments in three continuous spatial units of the main stream (MS), left tributary (LT), and right tributary (RT) of the Duliujian River basin were selected as the research objects. Therefore, 16S rDNA high-throughput sequencing and the PICRUSt platform were used to analyze the bacterial diversity and metabolic function of sediments in different units. The results showed that there were no significant differences for the Shannon and Simpson indices between the different river channels (P>0.05); however, the bacterial diversity indices of Chao1, ACE, Observed_species, and PD_whole_tree of LT and RT were significantly higher than those of the MS (P < 0.05). Moreover, the dominant bacterial phyla were not significantly different among the different river channels. In the present study, a total of 41 metabolic pathways were predicted based on the KEGG database, among which 34 metabolic pathways such as membrane transport, amino acid metabolism, and carbohydrate metabolism were significantly different. The functional metabolic abundance was higher in MS than that in RT and LT (P < 0.05); however, there was no significant difference between those of RT and LT (P>0.05). Except for exchangeable phosphorus and nitrate nitrogen, there were no significant differences in other physicochemical properties among the different river channels (P>0.05). Furthermore, pH and iron/aluminum-phosphorus were found to significantly affect the bacterial structure, and SOM and TN were found to significantly affect the bacterial metabolic function. Overall, the sediments of MS and tributaries (LT and RT) of the river were clearly divided into level-1 and level-2 groups, indicating that the metabolic abundance of the bacterial community can be used as the key indicator in the division of functional units at the river basin scale. This study provided a theoretical basis and technical support for ecological functional division and management in coastal river basins.
Key words: Duliujian River basin      surface sediment      bacterial diversity      metabolism function      functional division     

沉积物在水生态系统中发挥着重要作用, 并可作为地球化学元素循环和细菌群落生存的场所[1~4].与上覆水不同, 沉积物的特征不会在短时间内发生改变, 能够储存并“记录”一段时间内生态系统的改变, 有研究表明, 沉积物作为细菌群落的生存场所是生态系统多样性的基础[5~7].因此研究沉积物细菌多样性和功能代谢丰度的变化有助于深入了解水生态系统的功能特征及其影响因子.

细菌在生物地球化学循环过程中起着至关重要的作用[8~13].细菌群落的多样性、结构和功能代谢易受到环境变化和人为活动的影响, 进而影响到生态系统的平衡状态.例如, 盐胁迫会使得细菌群落的TAC循环、脂肪酸生物合成以及乙醛酸和二羟酸代谢等代谢途径进行重新匹配[14].此外, 不同生境条件对于细菌功能代谢的影响程度也不同.Lee等[15]研究了潮滩芦苇和碱蓬的根际和非根际的细菌群落的变化, 发现沉积物环境条件的改变显著影响了细菌群落的多样性和结构; Su等[16]对杭州湾沉积物细菌的氮转化途径进行了研究, 发现在人为干扰条件下磷素是影响细菌群落结构的关键因素, 并且盐分和无机氮分别作为自然和人为因素也会对细菌的氮代谢功能产生显著影响.

近年来, 细菌群落的某些代谢途径被用作生态系统健康评价和预测的手段之一[17~20].有研究发现细菌可以通过脂质代谢和异生物素代谢等代谢途径, 对持久性有机污染物进行降解, 并以此来确保自身的结构和功能恢复到原始水平[21, 22].目前针对细菌群落结构和功能代谢途径的研究很多, 但在水生态功能分区领域尚无研究先例, 沉积物中细菌群落的多样性、结构和功能代谢丰度能否反映研究区域生态因子和特色功能表达的区域单元划分尚不得而知.同时有研究表明, 细菌对环境变化的敏感性较高, 当所处环境发生改变时, 会通过自身的代谢途径做出应激性表达[23], 但细菌群落的代谢功能在相似等级的生态条件下, 或具有趋同性[24].因此, 本文试图验证细菌群落的多样性、结构和代谢丰度在相同等级河流沉积物中是否具有趋同性, 从而探讨细菌群落特征用于流域功能单元划分的可行性.

本文选择了天津独流减河流域上游汇水段作为研究区域, 通过16S rDNA测序技术研究细菌群落多样性、结构和功能在不同级别河流沉积物中的变化特征和影响因子, 并探讨以沉积物中细菌群落功能代谢丰度为依据进行功能单元划分的可行性, 以期为我国水生态功能区划及流域空间管理提供理论基础和技术支撑.

1 材料与方法 1.1 研究地点概况

独流减河属于海河水系, 是海河南系下游最大的支流.在上游汇水区, 属于天津西青区和静海区的交接地带, 由不同支流(子牙河、大清河和南运河)汇入组成, 属于人工开挖泄洪河道, 设计流量3 200 m3·s-1, 地处中纬度, 受亚欧大陆和太平洋的影响, 春季多风, 干旱少雨; 夏季炎热, 雨水集中; 秋季气爽, 冷暖适中; 冬季寒冷, 干燥少雪.区内无地带性特殊植被, 主要以狗尾草、芦苇、芦竹、马唐和牛筋草等禾本科植物为主.

1.2 样品采集和处理

本研究沉积物样品采集于2020年9月, 在独流减河汇水区进行布点采样(图 1), 根据干流流向划分为左岸支流(left tributary, LT)、河道干流(main stream, MS)和右岸支流(right tributary, RT)3个不同的空间单元.根据河网密度和支流分布情况在每个空间单元各设置5~7个采样点, 共设置了20个采样点, 在每个采样点随机采集3~4个沉积物样品剔除杂物后混合均匀.选择代表性样点(图 1中彩色填充的样点, 确保每条支流均有采样点布设)的部分样品装入离心管后立即放入干冰壶中保存, 随后在上海中科新生命生物科技有限公司利用Illumina MiSeq测序平台进行16S rDNA测序.剩余样品装入聚乙烯袋中带回实验室, 一部分置于-20℃冰箱中保存, 用于测定沉积物的含水率和氮的形态指标; 另一部分在自然状态下风干, 用于测定沉积物的pH、盐含量、有机质、总氮、总磷和磷形态指标.

彩色填充符号采样点沉积物进行理化性质测定和16S rDNA测序, 黑色填充符号采样点沉积物仅进行理化性质相关测定; 蓝色箭头表示河流流向 图 1 研究区域沉积物采样点位示意 Fig. 1 Sampling sites of sediments in this study

1.3 沉积物理化性质研究方法

本研究中, 沉积物的pH和盐含量(salt)分别采用pH计(Delta320, 梅特勒-托利多, 瑞士)和电导率仪(DDSJ-308A, 雷兹, 中国)测定; 有机质(SOM)采用重铬酸钾-外加热法测定; 硝态氮(NO3--N)、亚硝态氮(NO2--N)和铵态氮(NH4+-N)采用紫外分光光度法测定; 总氮(TN)采用元素分析仪(Vario EL Ⅲ, 艾力蒙塔, 德国)测定; 总磷(TP)通过电感耦合等离子体原子发射光谱(VISTA-MPX, 瓦里安, 美国)测定; 含水率(MC)采用烘干重法测定[24]; 本研究采用欧盟计划的标准测量测试方法(SMT法)提取沉积物中的可交换态磷(Ex-P)、铁铝结合态磷(Fe/Al-P)和钙结合态磷(Ca-P), 并用紫外分光光度法测定[25].

1.4 沉积物DNA、提取及测序技术

采用CTAB方法对样本的基因组DNA进行提取, 在1%琼脂糖凝胶上监测DNA浓度和纯度.根据测序区域的选择, 使用带有Barcode的特异性引物在V3~V4区域进行PCR扩增; 使用0.2μmol·L-1的正向引物ACTCCTACGGGAGGCAGCA和反向引物GGACTACHVGGGTWTCTAAT, 以及约10 ng的模板DNA.将PCR产物以等密度比混合, 用AxyPrepDNA凝胶提取试剂盒(AXYGEN)纯化混合PCR产物.最后, 使用Illumina Novaseq6000进行上机测序.

1.5 统计方法

本研究采用Surfer 15绘制研究区域采样点; 采用“ggsignif”、“ggcor”、“factoextra”和“psych”R包的统计方法对理化指标和微生物群落进行ANOVA检验、T检验、Mantel检验、Pearson相关性分析、主成分分析和Spearman相关性分析, 并通过软件Gephi 0.9.2、“ggcor”和“ggplot2”R包进行可视化展示.

2 结果与分析 2.1 沉积物理化性质变化

在左岸支流、干流和右岸支流这3个连续空间单元中, 沉积物理化性质的变化范围大小不同(图 2), pH变化范围为6.66~7.26, 含水率变化范围为22.07%~42.71%, 含盐量变化范围为0.010%~0.047%, ω(NH4+-N)的变化范围为5.52~24.12 mg·kg-1, ω(NO2--N)的变化范围为0.16~1.10 mg·kg-1, ω(Fe/Al-P)的变化范围为20.74~48.82 mg·kg-1, ω(Ca-P)的变化范围为111.86~319.24 mg·kg-1, ω(TN)的变化范围为515~2 580 mg·kg-1, ω(TP)的变化范围为181.35~1 116.514 mg·kg-1, ω(SOM)的变化范围为20.56~42.83 mg·g-1.而ω(NO3--N)呈显著下降趋势(P < 0.01), 变化范围为0.012~0.071 mg·kg-1; ω(Ex-P)呈显著上升趋势(P < 0.01), 变化范围为3.27~13.29 mg·kg-1.

空白三角形表示各空间单元沉积物中平均含量, 黑色圆点表示各空间单元所存在的异常值 图 2 独流减河流域不同等级河流沉积物中理化性质的变化 Fig. 2 Changes in physicochemical properties of sediments in different channels of Duliujian River basin

2.2 细菌群落α多样性变化

细菌群落α多样性的分布特征如图 3所示.方差分析显示, 支流沉积物中细菌的Ace、Chao1、Observed_species和PD_whole_tree指数均显著大于干流(P < 0.05).然而, 干流和支流细菌的Shannon指数和Simpson指数无显著差异(P>0.05).

空白三角形表示各空间单元沉积物中平均含量, 黑色圆点表示各空间单元所存在的异常值 图 3 独流减河流域不同等级河流沉积物细菌群落α多样性的变化 Fig. 3 Changes in α diversity of bacterial communities of sediments in different channels of Duliujian River basin

2.3 沉积物细菌群落优势门丰度变化特征

图 4显示了独流减河干流和支流沉积物细菌群落优势门的相对丰度情况, 根据其相对丰度的降序排列情况, 这些细菌门被鉴定为Proteobacteria(36.39%)、Actinobacteria(19.51%)、Chloroflexi(14.65%)、Acidobacteria(9.93%)、Bacteroidetes(6.09%)、Gemmatimonadetes(4.97%)、Firmicutes(3.12%)、Patescibacteria(1.37%)、Rokubacteria(0.84%)、Nitrospirae(0.80%)和others(2.34%), 共占整个细菌群落的97.66%.

图 4 独流减河流域不同等级河流沉积物中细菌群落在门水平上优势门的相对丰度 Fig. 4 Relative abundance of dominant phyla of bacterial communities of sediments in different channels of Duliujian River basin

除others以外, 这些优势细菌门在不同空间单元之间均没有显著差异(P>0.05), 但变化范围各不相同(图 4), Proteobacteria的变化范围为26.71%~57.90%, Actinobacteria的变化范围为3.64%~28.12%, Chloroflexi的变化范围为4.03%~32.20%, Acidobacteria的变化范围为2.33%~15.86%, Bacteroidetes的变化范围为1.86%~10.44%, Gemmatimonadetes的变化范围为0.96%~13.32%, Firmicutes的变化范围为0.64%~14.32%, Patescibacteria的变化范围为0.33%~3.70%, Rokubacteria的变化范围为0.20%~1.61%, Nitrospirae的变化范围为0.13%~4.45%.本研究发现others相对丰度表现为左支流显著高于右支流(P < 0.05), 变化范围为0.85%~3.77%.

2.4 沉积物中细菌群落KEGG数据库功能预测

基于KEGG数据库的预测结果, 图 5分别显示了干流与支流代谢丰度的差异, 共包含了膜运输、氨基酸代谢、碳水化合物代谢、复制与修复、能量代谢、不良特征、转运、辅助因子和维生素的代谢、脂质代谢和外源生物降解与代谢等41条功能代谢途径.经T检验, 共有34条代谢途径差异显著(P < 0.05), 其中, 干流沉积物细菌的代谢丰度显著高于左岸支流[图 5(a), P < 0.05]和右岸支流[图 5(b), P < 0.05], 但其在左右岸支流之间并无显著差异(P>0.05).

图 5 独流减河流域不同等级河流沉积物中基于KEGG数据库预测的细菌群落功能代谢丰度 Fig. 5 Functional metabolic abundance of bacterial community predicted by KEGG database in different channels of Duliujian River basin

2.5 沉积物理化性质与细菌群落相互关系

相关性热图显示了沉积物中理化性质和细菌参数的内部关系(图 6).所有理化性质按元素循环可以分为两大类, 例如, 在氮形态转化的过程中, NO3--N和含盐量呈极显著正相关(P < 0.001), 含水率和TN呈显著正相关(P < 0.05); 而在磷循环的过程中, Fe/Al-P与含水率和TN呈显著负相关性(P < 0.05), Ca-P和NH4+-N呈显著正相关性(P < 0.05).此外, SOM和TN呈显著正相关性(P < 0.01).当把细菌群落作为整体考虑其与沉积物理化性质的关系时, 细菌群落的多样性参数几乎不受沉积物理化性质的影响, 而细菌群落的优势门受pH和Fe/Al-P影响显著, 细菌群落的功能代谢则受TN和SOM影响显著.

圆圈大小表示Spearman相关系数大小; 线条宽度表示mantel test相关系数大小; ***、**和*分别代表在0.001、0.01、和0.05水平上的显著性 图 6 沉积物理化性质与细菌群落多样性、细菌群落优势门和细菌代谢功能之间的相关热图 Fig. 6 Correlation heatmap of physicochemical properties, bacterial diversity, bacterial dominant phyla, and bacterial functional metabolism of sediments in different channels of Duliujian River basin

2.6 细菌群落功能代谢丰度的网络联系

基于KEGG数据库预测的细菌功能代谢丰度进行主成分分析后发现, 前两个主成分共解释了不同等级河流92.3%的差异变化, 可以清晰地将干流和支流沉积物划分为Level-1和Level-2两个功能单元[图 7(a)].

(a) Level-1表示河道干流采样点, Level-2表示左右岸支流各采样点; (b)和(c)表示网络模型之间通过Spearman相关
(R>0.6, P < 0.05)进行相互关联, 红色和绿色圆圈分别表示河道干流和左右岸支流沉积物中细菌群落代谢途径
图 7 基于KEGG数据库预测的细菌群落功能代谢丰度的主成分分析和网络共现分析 Fig. 7 Principal component analysis and network co-occurrence analysis based on the functional metabolic abundance of bacterial community predicted by KEGG database

通过共现网络分析比较细菌群落代谢途径连接的复杂性, 不同功能单元呈现出不同的网络组成和拓扑特征.基于KEGG数据库所预测的细菌代谢丰度之间的Spearman相关性导致了更复杂的代谢途径结构(R>0.6; P < 0.05), 即Level-1[41个节点, 57条边, 平均度为6.268, 图密度0.157, 图 7(c)]和Level-2[41个节点, 1079条边, 平均度为26.317, 图密度0.658, 图 7(b)].总的来说, 支流相比于干流拥有更为复杂的网络结构, 细菌代谢关联更为紧密.

3 讨论

沉积物作为流域生态系统的重要组成部分, 在一定程度上可以反映河流上覆水环境的时空变化特征[26~29], 而沉积物主要通过理化性质和微生物组成的改变影响上覆水环境.在流域污染防治过程中, 统筹上下游、左右岸和干支流的关系, 协同治理是关键.在本研究中, 左岸支流沉积物中NO3--N的含量显著高于干流与右岸支流, 而右岸支流沉积物中Ex-P的含量显著高于干流和左岸支流.河流一般具有“左凸右凹”(左岸堆积, 右岸侵蚀)的环境属性, 即受到地转偏向力等因素的影响时河流右岸容易呈现出河床加深、堆积的现象, 而这一过程促使沉积物发生“再悬浮”, 间隙水中更多的溶解氧被释放到上覆水中, 氧的释放有利于硝化反应过程中NH4+-N向NO3--N的转变, 从而导致右岸支流沉积物中NO3--N的含量较高.磷属于典型的沉积型循环物质, 由外源输入或内源生物分解所产生的磷素, 在左岸冲刷作用下最终会在底泥中沉积, 增加了左岸支流沉积物中可交换磷的数量.除NO3--N和Ex-P之外, 其他沉积物的理化性质在流域不同空间单元之间均无显著差异, 说明干流和支流的差异并非主要由沉积物理化性质的改变引起, 也可能与沉积物中细菌区系的变化有关.有研究表明, 细菌群落可以反映流域生态系统的完整性特征, 并通过多样性和结构的改变适应外部环境的变迁[30], 因此, 不同生境条件下细菌群落的差异是可预测的.本研究发现, 除细菌群落的Simpson指数和Shannon指数在不同空间单元之间无显著差异性变化外, Chao1、Ace、Observed_species和PD_whole_tree指数等均表现为支流显著高于干流的特征, 可能和支流采样点的空间异质性更高有关, 空间异质性越高越有利于细菌种群结构的复杂性和多样性.

近年来, 细菌及其代谢途径被广泛用作评判生态系统污染修复状况的关键指标[31].因此, 探索细菌群落组成、结构和代谢功能对于了解区域生态系统的健康和潜在生态修复能力具有重要意义.在本研究中, 优势细菌门相对丰度之间并无显著差异.细菌群落在应对环境变化的同时, 能够利用环境中不同类型的营养物质[32~34], 应对环境变化所带来的扰动[35], 较为丰富的细菌群落会表现出相似的恢复力或抗性, 并具有相似的特征丰度[36].此外, 细菌群落在应对某种环境变化时, 通常会采取持久性休眠策略, 休眠中的细菌可以承受环境变化的扰动[37].因此, 河流沉积物中可能存在大量休眠细菌, 这增加了不同空间单元中优势菌群的相似程度.

尽管细菌群落优势门的相对丰度差异不大, 但细菌群落的代谢丰度具有显著差异.有研究表明, 外界环境的变化会影响细菌群落的功能代谢[38, 39].在本研究中, 基于KEGG数据库预测的细菌不同代谢途径的丰度值均表现出河道干流显著高于支流的趋势, 但支流之间并无显著差异, 该发现也证实了细菌群落的代谢水平在相似等级的生态条件下具有趋同性[24].经Mantel test分析表明, 沉积物中TN和SOM的含量对细菌群落的功能代谢影响显著, 说明在河道干流中细菌群落可能通过更快的生命周转过程不断为自身提供营养物质, 从而加快了该生境中细菌的代谢过程.另外, 氨基酸代谢在细菌群落中一直占据主导地位, 而氨基酸代谢可通过脱氮、转氨和联合脱氮作用, 分解成细菌所需的能量和CO2[40], 河道干流沉积物中氨基酸代谢的丰度较高, 从而提高了该生境中其他代谢途径的丰度值.细菌群落的功能代谢除了受到环境条件制约外, 还会受到空间变量的限制, 而且空间变量对细菌群落结构以及功能代谢的影响大于环境变量[41].流域汇水段不同河流空间单元的差异可能和细菌群落的生存习性密切相关, 在相对狭小或流速较慢的空间单元内, 细菌群落之间存在竞争模式, 或与其他生物相互作用, 如营养汲取、捕食等, 不利于细菌群落的生长繁衍和功能性状的表达.此外, 当生存空间的等级发生改变时, 细菌群落可以通过代谢途径准确地对周围环境做出适应性和应激性表达, 这些功能性状的表达同时会决定上覆水的健康状况[42].本研究发现, 河道干流相比较左右岸支流, 具有横截面宽、流速快和水深等特点(独流减河河道干流平均宽度为720.02 m; 左岸支流平均宽度为57.56 m; 右岸支流平均宽度为25.70 m), 更适合细菌群落的功能表达.因此, 细菌群落功能代谢丰度的变化同时预示着不同空间单元的变迁, 揭示了新的功能分区单元的形成.

Wood等[43]曾对湿地生态系统中细菌群落的代谢途径进行过研究, 结果显示, 在深度小于75 cm的水环境中进行有氧呼吸的细菌群落的相对丰度显著高于深度大于75 cm的水环境, 细菌代谢途径的变化与垂直尺度上环境因子的改变是一致的.在本研究中, 基于细菌群落代谢丰度数据的主成分分析可以清晰地将河道干流和左右岸支流划分为Level-1(河道干流)和Level-2(左右岸支流)两个不同的功能单元.其中, Level-1反映了自然特征总体趋势和生态系统综合服务功能的格局; Level-2是在Level-1的划分基础上进一步细化的结果, 能够更好地体现水系的功能特征[44].网格共现分析显示左右岸支流细菌群落的代谢途径之间关联性更紧密, 呈现出比河道干流更复杂的网络特征, 验证了基于细菌功能代谢丰度进行水生态功能区划分的合理性.此外, 为消除样本数量对网络共现分析结果产生的影响, 分别对Level-2的左右岸支流进行了共现分析, 结果表明左岸和右岸支流的网络结构均比干流复杂(左岸支流包含41个节点, 344条边, 平均度为8.39, 图密度为0.21; 右岸支流包含41个节点, 615条边, 平均度为15, 图密度为0.375).因此, 可以选择细菌群落的功能代谢丰度作为流域水生态功能分区的关键指标.但也应清醒地认识到本研究结论是基于小流域尺度上的探索性研究, 未来需在更大尺度和更多流域空间单元中做进一步的验证和评估.总之, 本研究是为提高流域水生态系统管理水平而进行的理论探索, 结果反映流域尺度上特色功能单元划分的微生物学基础, 同时可为我国水生态功能分区及流域空间管理提供理论基础和技术支持.

4 结论

(1) 本研究发现支流沉积物中细菌群落的Chao1、Ace、Observed_species和PD_whole_tree指数均显著高于干流, 但Simpson指数和Shannon指数在河道干流和左右岸支流之间无显著性差异.本研究同时发现, 沉积物理化性质的变化对细菌群落的多样性无显著影响, 故该类参数在流域功能单元划分中存在不确定性.

(2) 沉积物中优势菌门的相对丰度在不同等级河流中差异并不显著, 但pH和Fe/Al-P显著影响了细菌群落优势菌门的分布, 故河流沉积物中磷的异质性在其细菌群落的结构组成上占据一定的主导地位.

(3) 干流沉积物中细菌群落的功能代谢丰度显著高于支流, 但不同支流之间差异并不显著, 表明细菌群落的代谢丰度在同一等级的河道中具有趋同性, 适合用作流域功能单元划分的关键参数.研究同时发现, 细菌代谢功能受SOM和TN的影响显著, 故流域管理中需关注氮的输入情况.

参考文献
[1] Long Y C, Jiang J, Hu X J, et al. The response of microbial community structure and sediment properties to anthropogenic activities in Caohai wetland sediments[J]. Ecotoxicology and Environmental Safety, 2021, 211. DOI:10.1016/j.ecoenv.2021.111936
[2] 周石磊, 孙悦, 岳哿丞, 等. 白洋淀养殖区春季沉积物垂向微生物群落特征及驱动因素[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.
[3] 李明月, 杨雨虹, 米铁柱, 等. 海洋沉积物中细菌DNA和RNA水平群落差异[J]. 环境科学, 2020, 41(5): 2485-2495.
Li M Y, Yang Y H, Mi T Z, et al. Differences between DNA- and RNA-based bacterial communities in marine sediments[J]. Environmental Science, 2020, 41(5): 2485-2495.
[4] 于小彦, 张平究, 张经纬, 等. 城市河流沉积物微生物量分布和群落结构特征[J]. 环境科学学报, 2020, 40(2): 585-596.
Yu X Y, Zhang P J, Zhang J W, et al. Characteristics of distribution patterns of microbial biomass and community structures in the sediments from urban river[J]. Acta Scientiae Circumstantiae, 2020, 40(2): 585-596.
[5] Zhang M Z, Wu Z J, Sun Q Y, et al. Response of chemical properties, microbial community structure and functional genes abundance to seasonal variations and human disturbance in Nanfei River sediments[J]. Ecotoxicology and Environmental Safety, 2019, 183. DOI:10.1016/j.ecoenv.2019.109601
[6] Zhang L, Li X C, Fang W K, et al. Impact of different types of anthropogenic pollution on bacterial community and metabolic genes in urban river sediments[J]. Science of the Total Environment, 2021, 793. DOI:10.1016/j.scitotenv.2021.148475
[7] 雷远达, 燕鸿宇, 钟彩英, 等. 沉积物记录揭示的深圳湾红树林生态系统稳态转变[J]. 生态学报, 2020, 40(23): 8539-8548.
Lei Y D, Yan H Y, Zhong C Y, et al. Regime shift recorded by sediments from the Futian mangrove ecosystem in the Shenzhen Bay[J]. Acta Ecologica Sinica, 2020, 40(23): 8539-8548.
[8] Sanyal S K, Brugger J, Etschmann B, et al. Metal resistant bacteria on gold particles: Implications of how anthropogenic contaminants could affect natural gold biogeochemical cycling[J]. Science of the Total Environment, 2020, 727. DOI:10.1016/j.scitotenv.2020.138698
[9] 余炎炎, 李梦莎, 刘啸林, 等. 大兴安岭典型永久冻土土壤细菌群落组成和多样性[J]. 微生物学通报, 2020, 47(9): 2759-2770.
Yu Y Y, Li M S, Liu X L, et al. Soil bacterial community composition and diversity of typical permafrost in Greater Khingan Mountains[J]. Microbiology China, 2020, 47(9): 2759-2770.
[10] 孙峰, 田伟, 张菲, 等. 丹江口库区库滨带植被土壤细菌群落多样性及PICRUSt功能预测分析[J]. 环境科学, 2019, 40(1): 421-429.
Sun F, Tian W, Zhang F, et al. Composition and predictive functional analysis of rhizosphere bacterial communities in riparian buffer strips in the Danjiangkou Reservoir, China[J]. Environmental Science, 2019, 40(1): 421-429.
[11] 李萍, 史荣久, 赵峰, 等. 大兴安岭落叶松林不同演替阶段土壤细菌群落结构与功能潜势[J]. 应用生态学报, 2019, 30(1): 95-107.
Li P, Shi R J, Zhao F, et al. Soil bacterial community structure and predicted functions in the larch forest during succession at the Greater Khingan Mountains of Northeast China[J]. Chinese Journal of Applied Ecology, 2019, 30(1): 95-107.
[12] 王蕊, 吴宪, 李刚, 等. 林地转型耕地对东北丘陵区白浆土cbbL细菌群落丰度和结构的影响[J]. 环境科学, 2019, 40(12): 5561-5569.
Wang R, Wu X, Li G, et al. Effects of conversion of forest to arable land on the abundance and structure of the cbbL-harboring bacterial community in albic soil of the hilly region of Northeast China[J]. Environmental Science, 2019, 40(12): 5561-5569.
[13] 张双双, 靳振江, 贾远航, 等. 岩溶地区不同土地利用方式土壤固碳细菌群落结构特征[J]. 环境科学, 2019, 40(1): 412-420.
Zhang S S, Jin Z J, Jia Y H, et al. Community structure of CO2-fixing soil bacteria from different land use types in karst areas[J]. Environmental Science, 2019, 40(1): 412-420. DOI:10.3969/j.issn.1000-6923.2019.01.049
[14] Wang M, Zhao S W, Wang L F, et al. Salt stress-induced changes in microbial community structures and metabolic processes result in increased soil cadmium availability[J]. Science of the Total Environment, 2021, 782. DOI:10.1016/j.scitotenv.2021.147125
[15] Lee H, Heo Y M, Kwon S L, et al. Environmental drivers affecting the bacterial community of intertidal sediments in the Yellow Sea[J]. Science of The Total Environment, 2021, 755. DOI:10.1016/j.scitotenv.2020.142726
[16] Su Z G, Dai T J, Tang Y S, et al. Sediment bacterial community structures and their predicted functions implied the impacts from natural processes and anthropogenic activities in coastal area[J]. Marine Pollution Bulletin, 2018, 131: 481-495. DOI:10.1016/j.marpolbul.2018.04.052
[17] 程豹, 望雪, 徐雅倩, 等. 澜沧江流域浮游细菌群落结构特征及驱动因子分析[J]. 环境科学, 2018, 39(8): 3649-3659.
Cheng B, Wang X, Xu Y Q, et al. Bacterioplankton community structure in the Lancang River basin and the analysis of its driving environmental factors[J]. Environmental Science, 2018, 39(8): 3649-3659.
[18] 万甜, 何梦夏, 任杰辉, 等. 渭河流域水体细菌群落的环境响应及生态功能预测[J]. 环境科学, 2019, 40(8): 3588-3595.
Wan T, He M X, Ren J H, et al. Environmental response and ecological function prediction of aquatic bacterial communities in the Weihe River Basin[J]. Environmental Science, 2019, 40(8): 3588-3595.
[19] Hanif M A, Guo Z M, Moniruzzaman M, et al. Plant taxonomic diversity better explains soil fungal and bacterial diversity than functional diversity in restored forest ecosystems[J]. Plants, 2019, 8(11). DOI:10.3390/plants8110479
[20] 张拓, 徐飞, 怀宝东, 等. 松花江下游沿江湿地土地利用变化对土壤细菌群落多样性的影响[J]. 环境科学, 2020, 41(9): 4273-4283.
Zhang T, Xu F, Huai B D, et al. Effects of land use changes on soil bacterial community diversity in the riparian wetland along the downstream of Songhua River[J]. Environmental Science, 2020, 41(9): 4273-4283.
[21] Miao Y, Heintz M B, Bell C H, et al. Profiling microbial community structures and functions in bioremediation strategies for treating 1, 4-dioxane-contaminated groundwater[J]. Journal of Hazardous Materials, 2021, 408. DOI:10.1016/j.jhazmat.2020.124457
[22] Zhao R X, Liu J, Feng J, et al. Microbial community composition and metabolic functions in landfill leachate from different landfills of China[J]. Science of the Total Environment, 2021, 767. DOI:10.1016/j.scitotenv.2020.144861
[23] Zhao L J, Huang Y X, Paglia K, et al. Metabolomics reveals the molecular mechanisms of copper induced cucumber leaf (Cucumis sativus) senescence[J]. Environmental Science & Technology, 2018, 52(12): 7092-7100.
[24] Liu F D, Mo X, Kong W J, et al. Soil bacterial diversity, structure, and function of Suaeda salsa in rhizosphere and non-rhizosphere soils in various habitats in the Yellow River Delta, China[J]. Science of the Total Environment, 2020, 740. DOI:10.1016/j.scitotenv.2020.140144
[25] 黄宏, 罗娜, 张硕. 海州湾海洋牧场表层沉积物中磷形态及释放通量[J]. 水产学报, 2020, 44(6): 959-968.
Huang H, Luo N, Zhang S. Forms and releasing flux of phosphorus in the surface sediments in Haizhou Bay marine ranching[J]. Journal of Fisheries of China, 2020, 44(6): 959-968.
[26] Li M M, Kong F L, Li Y, et al. Impact of sub-watershed characteristics and changes on sediment DOM quantity and quality over Jiaozhou Bay[J]. Ecological Indicators, 2020, 118. DOI:10.1016/j.ecolind.2020.106777
[27] 郝文超, 王从锋, 杨正健, 等. 氧化还原循环过程中沉积物磷的形态及迁移转化规律[J]. 环境科学, 2019, 40(2): 640-648.
Hao W C, Wang C F, Yang Z J, et al. Speciation and transformation of phosphorus in sediments during the redox cycle[J]. Environmental Science, 2019, 40(2): 640-648.
[28] 杜先, 荀凡, 王亚蕊, 等. 蓝藻碎屑堆积对湖泊沉积物矿化特征的影响[J]. 湖泊科学, 2020, 32(6): 1671-1682.
Du X, Xun F, Wang Y R, et al. Effects of algal-detritus accumulation on the mineralization characteristics of lake sediment[J]. Journal of Lake Sciences, 2020, 32(6): 1671-1682.
[29] 汤宁, 李如忠, 王聿庆, 等. 污水厂尾水受纳河段沉积物磷形态及释放风险效应[J]. 环境科学, 2020, 41(2): 801-808.
Tang N, Li R Z, Wang Y Q, et al. Phosphorus forms and release risk of sediments in urban sewage treatment plant effluent and receiving stream reach[J]. Environmental Science, 2020, 41(2): 801-808.
[30] Liao H H, Yen J Y, Guan Y J, et al. Differential responses of stream water and bed sediment microbial communities to watershed degradation[J]. Environment International, 2020, 134. DOI:10.1016/j.envint.2019.105198
[31] Obieze C C, Chikere C B, Selvarajan R, et al. Functional attributes and response of bacterial communities to nature-based fertilization during hydrocarbon remediation[J]. International Biodeterioration & Biodegradation, 2020, 154. DOI:10.1016/j.ibiod.2020.105084
[32] Lee S A, Kim J M, Kim Y, et al. Different types of agricultural land use drive distinct soil bacterial communities[J]. Scientific Reports, 2020, 10. DOI:10.1038/s41598-020-74193-8
[33] 吴玲. 富营养化太湖沉积物中微生物群落及对环境因子的响应[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 and Environmental Biology, 2019, 25(6): 1470-1476.
[34] 胡越航, 宣丽霞, 裘琼芬. 海水养殖尾水直排对河道沉积物细菌群落的影响[J]. 环境科学, 2019, 40(9): 4185-4194.
Hu Y H, Xuan L X, Qiu Q F. Bacterial community composition in channel sediment in response to mariculture wastewater drainage[J]. Environmental Science, 2019, 40(9): 4185-4194.
[35] Giudice A L, Caruso G, Rizzo C, et al. Bacterial communities versus anthropogenic disturbances in the Antarctic coastal marine environment[J]. Environmental Sustainability, 2019, 2(3): 297-310. DOI:10.1007/s42398-019-00064-2
[36] Jiao S, Wang J M, Wei G H, et al. Dominant role of abundant rather than rare bacterial taxa in maintaining agro-soil microbiomes under environmental disturbances[J]. Chemosphere, 2019, 235: 248-259. DOI:10.1016/j.chemosphere.2019.06.174
[37] Xie G X, Kong X L, Kang J L, et al. Community-level dormancy potential regulates bacterial beta-diversity succession during the co-composting of manure and crop residues[J]. Science of the Total Environment, 2021, 772. DOI:10.1016/j.scitotenv.2021.145506
[38] Reddy B, Pandey J, Dubey S K. Assessment of environmental gene tags linked with carbohydrate metabolism and chemolithotrophy associated microbial community in River Ganga[J]. Gene, 2019, 704: 31-41. DOI:10.1016/j.gene.2019.04.004
[39] Cao H Y, Du Y J, Gao G L, et al. Afforestation with Pinus sylvestris var. mongolica remodelled soil bacterial community and potential metabolic function in the Horqin Desert[J]. Global Ecology and Conservation, 2021, 29. DOI:10.1016/j.gecco.2021.e01716
[40] Mohsenpour S F, Hennige S, Willoughby N, et al. Integrating micro-algae into wastewater treatment: a review[J]. Science of the Total Environment, 2021, 752. DOI:10.1016/j.scitotenv.2020.142168
[41] Wang H H, Tian T, Gong Y, et al. Both environmental and spatial variables affect bacterial functional diversity in mangrove sediments at an island scale[J]. Science of the Total Environment, 2021, 753. DOI:10.1016/j.scitotenv.2020.142054
[42] Zhang L, Wang Z Y, Cai H, et al. Long-term agricultural contamination shaped diversity response of sediment microbiome[J]. Journal of Environmental Sciences, 2021, 99: 90-99. DOI:10.1016/j.jes.2020.06.013
[43] Wood J R, Burge O R, Bolstridge N, et al. Vertical distribution of prokaryotes communities and predicted metabolic pathways in New Zealand wetlands, and potential for environmental DNA indicators of wetland condition[J]. PLoS One, 2021, 16(1). DOI:10.1371/journal.pone.0243363
[44] 李法云, 范志平, 张博, 等. 辽河流域水生态功能一级分区指标体系与技术方法[J]. 气象与环境学报, 2012, 28(5): 83-89.
Li F Y, Fan Z P, Zhang B, et al. An index system and a technical method for the first level aquatic ecoregion in the Liaohe River basin[J]. Journal of Meteorology and Environment, 2012, 28(5): 83-89.