环境科学  2025, Vol. 46 Issue (5): 2805-2815   PDF    
金沙江下游梯级水库浮游细菌群落结构及生态过程
陈露欣1,2, 李哲1, 汤琼1,2, 鲁伦慧1     
1. 中国科学院重庆绿色智能技术研究院, 重庆 400714;
2. 中国科学院大学重庆学院, 重庆 400714
摘要: 浮游细菌是淡水生态系统的重要组成部分, 对碳氮磷等生源要素的生物地球化学循环过程具有重要作用. 大坝运行会导致水库的水文条件发生显著变化, 影响淡水生态系统. 探究不同水文条件下浮游细菌群落结构和生态过程是了解大坝对河流生态系统影响的关键. 以金沙江下游溪洛渡-向家坝梯级水库为研究对象, 利用16S rRNA扩增测序分析了水库不同调度运行下水体中浮游细菌群落结构和群落构建生态过程的差异. 结果表明:Simpson指数、Chao指数、水温、溶解性有机碳(DOC)浓度和总磷(TP)浓度在低水位时期更高, 水力停留时间、pH、水位在高水位时期更高. 高、低水位时期浮游细菌群落组成差异显著. 绿弯菌门(Chloroflexi)、芽单胞菌门(Gemmatimonadota)、蓝细菌门(Cyanobacteria)、硝化螺旋菌门(Nitrospirota)和噬氢菌门(Hydrogenedentes)对水位变化环境敏感. 蓝细菌门(Cyanobacteria)和硝化螺旋菌门是浮游细菌网络的关键物种. 低水位时期浮游细菌网络关键物种类别和丰度高于高水位时期, 低水位时期浮游细菌网络复杂性和稳定性更高. 高、低水位时期浮游细菌群落构建均以随机性过程为主导. 因此, 水库不同调度运行影响了浮游细菌群落结构和生态过程, 进而影响了碳氮元素的循环过程. 研究结果对于认识筑坝蓄水下河库生境水生态系统结构及功能, 以及指导高坝大库管理具有重要的学术价值与科学意义.
关键词: 水库调度运行      梯级水库      浮游细菌      网络特性      群落构建     
Structure and Ecological Processes of the Bacterioplankton Community in the Cascade Reservoirs of the Lower Reaches of the Jinsha River
CHEN Lu-xin1,2 , LI Zhe1 , TANG Qiong1,2 , LU Lun-hui1     
1. Chongqing Institute of Green and Intelligence Technology, Chinese Academy of Sciences, Chongqing 400714, China;
2. Chongqing School, University of Chinese Academy of Sciences, Chongqing 400714, China
Abstract: Bacterioplankton are an important component of freshwater ecosystems and play a crucial role in the biogeochemical cycling of biogenic elements, such as carbon, nitrogen, and phosphorus. The operation of dams can cause notable changes in the hydrological conditions of reservoirs, affecting freshwater ecosystems. Exploring the community structure and ecological processes of bacterioplankton under different hydrological conditions is key to understanding the impact of dams on river ecosystems. Taking the Xiluodu-Xiangjiaba cascade reservoirs in the lower reaches of the Jinsha River as the research object, 16S rRNA amplification sequencing was used to analyze the differences in the community structure and ecological process of bacterioplankton in the water under different scheduling operations of the reservoir. The results indicated that the Simpson index, Chao index, water temperature, dissolved organic carbon (DOC) concentration, and total phosphorus (TP) concentration were higher during the low water level period, whereas the hydraulic retention time, pH, and water level were higher during the high water level period. Significant differences were observed in the composition of bacterioplankton communities between the high and low water level periods. Chloroflexi, Gemmatimonadota, Cyanobacteria, Nitrospirota, and Hydrogenedentes were sensitive to changes in water level environments. Cyanobacteria and Nitrospirota were key species in the bacterioplankton networks. During the low water level period, the diversity and abundance of key species in the bacterioplankton network were higher than those during the high water level period, resulting in a more complex and stable network. During both high and low water level periods, the assembly of bacterioplankton communities were dominated by stochastic processes. Therefore, different scheduling operations of reservoirs affected the community structure and ecological processes of bacterioplankton, thereby affecting the cycling process of carbon and nitrogen elements. The research results have important academic value and scientific importance for understanding the structure and function of aquatic ecosystems in river and reservoir habitats under dam construction and water storage, as well as guiding the management of high dams and large reservoirs.
Key words: reservoir operation scheduling      cascade reservoir      bacterioplankton      network characteristics      community assembly     

大坝是对河流最大的人为干扰之一[1], 在协调流域内水电、供水、防洪、航运和灌溉[2]的同时, 导致河流主干及其相关水系的水位变化, 影响微生物群落和营养物的组成[3], 导致物种多样性下降和基本生态系统服务功能的丧失[4~7].

浮游细菌是淡水生态系统的重要组成部分, 在参与水生生态系统的关键生物地球化学循环和水生微生物食物网中扮演重要角色[8].有研究表明, 长江湖库区域的物种丰富度和细菌的代谢活动低于河流区域[9], 表层的浮游细菌更易受到环境扰动的影响[10].在水库调度运行影响下, 水位、水力停留时间、水量、悬浮泥沙流量和流速等水文因素通过改变生源物质的输移进而改变浮游细菌群落[11, 12].在生态系统中, 时间和空间上共存的物种通过物质、能量和信息的交换产生联系, 形成复杂的相互作用, 从根本上影响着种群、群落和生态系统过程, 对淡水生态系统的生物多样性、生态系统服务功能以及生态系统对环境变化的响应有着深刻的影响[13].生态网络是对生态系统中物质、能量流动进行模拟的结构模型, 在群落尺度上, 物种间可以同时以多种方式发生相互作用, 共同组成了生物共现网络(co-occurrence networks).而通过解析生物共现网络的拓扑结构可以作为水生态系统中长期动态指示和水生态监测指标.理解群落构建过程可以解释生物多样性格局的形成和维持, 对于生物多样性保护起到关键作用[14].浮游细菌群落构建过程受到环境过滤、竞争和捕食等相互作用(确定性因素)和物种死亡、出生(随机性因素)等因素的影响[15].例如, 随着盐度的增加, 确定性过程逐渐主导浮游细菌群落的聚集, 从而影响微生物群落结构[16].同样的, 大坝的运行会影响河流环境的物理化学性质, 变化的水温[17~19]、盐度[20]、pH[21, 22]和养分[23]等因素会影响水生生态系统中的浮游细菌群落构建过程.通过对浮游细菌群落构建机制的研究, 可以了解到不同生态过程对群落内物种组成和多样性维持贡献的大小[24].

金沙江属长江上游河段, 金沙江下游河段渡口至宜宾段长782 km, 落差729 m, 水力资源丰富.溪洛渡-向家坝梯级电站是金沙江下游最后两级电站, 溪洛渡水电站位于四川省雷波县和云南省永善县交界的金沙江下游江段溪洛渡峡谷.向家坝水电站坝址位于向家坝峡谷出口处, 左岸是四川省宜宾县, 右岸为云南省水富县, 电站的开发任务以发电为主, 同时改善通航条件, 结合防洪和拦沙, 兼顾灌溉, 并且具有为上游梯级电站起到反调节的作用.根据环评报告信息, 溪洛渡水库正常蓄水位600 m时水库回水长199 km, 相应库容115.7亿m3, 调节库容64.6亿m3, 具有不完全年调节性.向家坝水电站水库正常蓄水位380 m, 相应库容49.77亿m3, 调节库容9.03亿m3, 电站总装机容量6 400 MW.梯级水库的建设对金沙江下游地区的流域连通性和区域生态环境造成负面影响[25], 不仅会影响库区水环境和水生生态系统的演替, 还影响了氮、磷等营养盐输送[26], 使库区水体自净能力减弱[27].然而, 对梯级水库不同调度运行下, 浮游细菌群落结构和生态过程的研究较少.因此, 以金沙江下游溪洛渡-向家坝梯级水库为研究对象, 基于16S rRNA高通量测序技术, 探究在大坝不同调度运行下, 不同时期浮游细菌群落结构和相互作用, 分析不同水库运行时期浮游细菌群落构建机制, 通过揭示不同水位时期浮游细菌的生态模式, 旨在为梯级水库水生态研究及管理提供科学依据.

1 材料与方法 1.1 研究区域与样品采集

本研究以金沙江下游溪洛渡(XLD)、向家坝(XJB)梯级水库为研究对象, 共布设了15个采样点[图 1(a)].分别于2018年高水位(1月、3月、11月)和低水位(5月、7月、9月)时期开展了6次野外采样.采样点位信息如表 1所示.采样点位坐标和与大坝距离均利用Google Earth量取获得.使用10 L的采水桶采集河道中泓线表层0.5 m水样, 并用无菌瓶将采集的水样充分混合.将其中的2 L水样使用0.22 µm混合纤维滤膜(Millipore, Bedford, USA)过滤.滤膜保存至液氮罐中运输回实验室, 放置于-86℃超低温冰箱保存至DNA提取.原水低温运输回实验室进行理化指标测试.

图 1 溪洛渡-向家坝梯级水库采样点及水位变化 Fig. 1 Sampling points and water level changes in Xiluodu-Xiangjiaba cascade reservoirs

表 1 金沙江下游溪洛渡-向家坝梯级水库采样点位信息1) Table 1 XLD-XJB cascade reservoirs sampling point distribution information

溪洛渡-向家坝梯级水库因同时需要满足防洪和发电的功能, 根据环评报告中溪洛渡-向家坝水库工程运行规程, 溪洛渡水库正常蓄水位600 m, 汛期(6月1日至9月10日)为满足汛期防洪及发电需求, 限制水位560 m运行, 死水位540 m. 9月中旬开始蓄水, 9月底蓄至正常蓄水位, 12月下旬至翌年5月底为供水期, 5月底水库降至死水位.此外, 根据2018年采样时期溪洛渡-向家坝水位变化图[图 1(b)]可知, 2018年1月、3月和11月处于水库调度运行的高水位运行期, 采样时5月、7月和9月处于水库调度运行的低水位时期.

1.2 理化指标分析

在采样现场利用多参数水质分析仪(YSI EXO2, YSI Incorporated, USA)测定pH、水温(water temperature, WT)和溶解氧(dissolved oxygen, DO).水样经过0.45 μm玻璃纤维膜(Whatman, Membrane Solutions, USA)过滤后, 使用总有机碳分析仪(Vario TOC select, Elementar, GER)测定溶解性有机碳(dissolved organic carbon, DOC).原水采用过硫酸钾消解-钼酸铵分光光度法测定总磷(total phosphorus, TP), 采用过硫酸钾消解-紫外分光光度法测定总氮(total nitrogen, TN).水位(water level, WL)和水力停留时间(hydraulic retention time, HRT)数据来自三峡集团.

1.3 水体eDNA提取、PCR扩增和高通量测序

使用FastDNA SPIN kit(Mo Bio Laboratories, USA)进行微生物DNA提取, 1%琼脂糖凝胶电泳和NanoDrop2000对DNA提取效果进行质检和浓度测定, 提取的DNA用于质量检测、扩增、纯化和测序.采用引物338F(5′-ACTCCTACGGGAGGCAGCAGC-3′)和806R(5′-GGACTACHVGGGTWTCTAAT-3′)进行PCR扩增[28].生物测序在上海美吉生物医药科技有限公司的Illumina MiSeq PE250平台进行.

利用QIIME软件对扩增产物的原始序列进行质量过滤, UPARSE(版本7.1)以97%的相似性对OTU进行聚类.嵌合序列被UCHIME识别, 随后被剔除.使用SILVA 138/16S_bacteria数据库的RDP分类器对序列进行分类.利用QIIME2进行质量控制、物种注释和统计计算.

1.4 共现网络(co⁃occurrence network)分析

微生物共现网络分析是研究微生物群落结构及其内部相互作用的常用工具[29].为提高数据集的可靠性, 首先保留出现频率超过20%的样本中存在的OTU来构建浮游细菌共现网络.然后利用“psych”包计算这些OTU之间所有可能的成对Spearman秩相关(r), 随后将具有显著统计学意义和可靠(P<0.05, |r|>0.6)的OTU筛选出来进行共现网络构建[30].拓扑参数利用“igraph”包分析, 包括平均度(average degree)、聚类系数(clustering coefficient)、平均路径长度(average path length)、网络直径(diameter)、图密度(density)和模块化(modularity).较高的平均度、聚类系数、图密度和模块化以及较低的平均路径长度和直径表明网络连接性更强, 反映了网络的潜在复杂性更高[31, 32].通过平均拓扑参数的标准化分数得到一个指数来反映网络的复杂性.平均路径长度和网络直径被计算为变量(x-1)的倒数.微生物群落的稳定性由模块化和关键物种决定[33].因此, 利用模块化和关键物种衡量浮游细菌群落稳定性.关键物种对生态系统稳定性具有重要作用[34].基于ZiPi分析以识别不同水位时期浮游细菌网络节点的拓扑作用.在微生物生态网络中, 节点表示群落中的物种, 不同节点的拓扑学角色可以作为确定关键物种的依据.模块枢纽(module hubs)、网络枢纽(network hubs)和连接器(connectors)区域内的节点在模块间和各自的模块内起重要作用, 因此可作为判断其是否为关键物种的指标[35].所有的分析在R(版本4.2.2)中进行, 在Gephi(版本0.10.1)中实现共现网络可视化.

1.5 群落构建机制(community assembly)

应用Stegen等 [36, 37]开发的零建模方法分析确定不同水位时期浮游细菌群落构建机制.在该方法中, β平均最近分类单元距离(β mean nearest taxon distance, βMNTD)表示群落之间的成对系统发育更替, β最近分类单元指数(β nearest taxon index, βNTI)表示由观察到的βMNTD之间的标准差计算出的环境影响. βNTI<-2表示同质选择, βNTI>2表示异质选择, 这两者都是确定性过程. |βNTI|<2表示随机性过程.通过生成999个局部群落构建排列来获得Bray-Curtis相异性分布, 以及通过比较经验观察到的Bray-Curtis和模拟的零分布来计算修正的Raup-Crick指数(RCBray).根据修正的RCBray, 随机过程可分为3个生态过程:扩散限制(RCBray>0.95)、同质扩散(RCBray<-0.95)和生态漂变|RCBray|<0.95[36, 37].

1.6 其他生物信息统计分析

使用“ggplot2”、“ggpubr”和“cowplot”包进行主成分分析(principal component analysis, PCA)和绘图.应用随机森林算法(random forests)来帮助识别对不同水位运行时期的敏感性物种.使用R包“randomForest”构建了100个随机森林重复, 并使用函数“combine”将它们合并, 以获得稳健的重要性排序.最后通过R包“strucchange”实现函数断点, 计算基于线性回归模型的最优断点, 识别对高低水位条件变化响应最强烈的敏感门.分析和绘图在R(版本4.2.2)、SPSS和Origin(2022)中实现.

2 结果与分析 2.1 不同水库运行时期浮游细菌群落α-多样性和环境要素

在金沙江下游溪洛渡-向家坝梯级水库, 高水位和低水位时期浮游细菌α-多样性指数和水体理化指标如图 2所示, 通过比较各个指标的均值, 发现低水位时期浮游细菌群落Simpson指数(高水位时期平均值:0.95, 低水位时期平均值:0.96)和Chao指数(高水位:947.60, 低水位:1 036.03)显著高于高水位时期.河流水体pH和WT在高低水位时期差异显著. ρ(DOC)(高水位:2.27 mg·L-1, 低水位:2.99 mg·L-1)和ρ(TP)(高水位:0.07 mg·L-1, 低水位:0.09 mg·L-1)在低水位时期显著高于高水位时期.水力停留时间(高水位:12.23 d, 低水位:6.33 d)和水位(高水位:468.88 m, 低水位:439.42 m)在高水位时期显著高于低水位时期.Shannon指数、Pd、DO和TN浓度无显著差异.

Shannon表示香农指数, Simpson表示辛普森指数, Chao表示物种丰富度指数, PD表示系统发育多样性指数, HRT表示水力停留时间, DO表示溶解氧, WT表示水温, WL表示水位, TP表示总磷, DOC表示溶解性有机碳, TN表示总氮, NS表示P>0.05, *表示0.01<P<0.05, **表示0.01<P<0.001, ***表示P<0.001 图 2 高水位和低水位时期表层水体理化指标和浮游细菌α-多样性指数 Fig. 2 Physicochemical parameters of surface water during high and low water levels, bacterioplankton α-diversity index

将不同水库运行时期(高水位、低水位)的水位变化和水力停留时间分别与浮游细菌α⁃多样性进行线性拟合分析(图 3).研究发现, 不同水位时期浮游细菌α-多样性指数与水位变化均正相关, 而浮游细菌α-多样性指数与水力停留时间均负相关.

Shannon表示香农指数, Simpson表示辛普森指数, Chao表示物种丰富度指数, PD表示系统发育多样性指数, ΔWL表示水位变化, HRT表示水力停留时间, 阴影区域表示95%的置信区间;NS表示P>0.05, *表示0.01<P<0.05, **表示0.01<P<0.001, ***表示P<0.001 图 3 高水位和低水位时期浮游细菌α-多样性指数与水位变化和水力停留时间线性拟合 Fig. 3 Linear fitting of the α-diversity index of bacterioplankton during high and low water levels with water and hydraulic retention time

2.2 溪洛渡-向家坝梯级水库浮游细菌群落组成

高水位和低水位时期浮游细菌群落由32个相同的门组成, 但是在丰度上存在差异。高水位时期浮游细菌群落以放线菌门(Actinobacteriota)、变形菌门(Proteobacteria)、拟杆菌门(Bacteroidota)、酸杆菌门(Acidobacteriota)和浮霉菌门(Planctomycetota)为主, 分别占46.91%、23.53%、12.13%、4.96%和3.34%.低水位时期浮游细菌群落以放线菌门、变形菌门、拟杆菌门、蓝细菌门(Cyanobacteria)和浮霉菌门为主, 分别占44.96%、23.36%、12.48%、5.43%和3.95%[图 4(a)].

(a)不同时期浮游细菌群落组成;(b)基于主成分分析(PCA)分析浮游细菌群落差异;(c)基于随机森林分类的不同时期敏感门鉴定, 阴影部分表示高、低水位时期敏感的物种 图 4 溪洛渡-向家坝梯级水库高水位和低水位时期浮游细菌群落变化 Fig. 4 Changes in bacterioplankton communities during high and low water levels in XLD-XJB cascade reservoirs

通过PCA分析发现PC1轴对高水位和低水位时期的浮游细菌群落有明显的区分[图 4(b)].基于随机森林进一步比较了高水位和低水位两个时期浮游细菌群落的敏感性物种[图 4(c)], 发现分别是绿弯菌门(Chloroflexi)、芽单胞菌门(Gemmatimonadota)、蓝细菌门、硝化螺旋菌门(Nitrospirota)和噬氢菌门(Hydrogenedentes).利用Mantel检验将浮游细菌群落和环境变量进行了相关性分析, 如图 5所示.总体而言, 浮游细菌群落与pH、DO和DOC等理化指标有很强的相关性(P<0.01).

*表示0.01<P<0.05, **表示0.01<P<0.001, ***表示P<0.001 图 5 环境因子与浮游细菌群落的相关性分析 Fig. 5 Correlation analysis between environmental factors and planktonic bacterial communities

2.3 溪洛渡-向家坝梯级水库浮游细菌共现网络及关键物种

本研究通过构建共现网络来评估浮游细菌群落共现模式[图 6(a)].为了确定高水位和低水位对浮游细菌网络复杂性和群落稳定性的影响, 笔者计算了高低水位时期浮游细菌网络的拓扑参数, 如表 2所示, 再对拓扑参数进行计算, 得到浮游细菌网络的复杂性指数, 高水位时期浮游细菌网络复杂性指数为0.17, 低水位时期浮游细菌复杂性指数为0.83, 即低水位时期浮游细菌网络更复杂.基于ZiPi识别关键物种, 结果表明, 高、低水位时期浮游细菌网络都有模块中心点和连接节点, 但在两个时期的浮游细菌网络中都没有观察到网络枢纽[图 6(b)].低水位时期浮游细菌网络关键物种类别显著多于高水位时期.低水位时期浮游细菌网络有34个关键物种, 其中有20个被归类为细菌网络中的模块枢纽, 14个被归类为高度连接到多个模块的连接器.关键物种在门水平上属于浮霉菌门、放线菌门、变形菌门、蓝细菌门、酸杆菌门、拟杆菌门、厚壁菌门(Firmicutes)和蛭弧菌门(Bdellovibrionota).高水位时期浮游细菌网络有19个关键物种, 其中有6个被归类为细菌网络中的模块枢纽, 13个被归类为高度连接到多个模块的连接器.关键物种在门水平上属于蓝细菌门、拟杆菌门、放线菌门、变形菌门、髌骨菌门(Patescibacteria)、厚壁菌门、硝化螺旋菌门和未分类的门.不同水库调度运行时期浮游细菌网络没有相同的关键物种.低水位时期浮游细菌网络模块化指数(0.29)高于高水位时期(0.26), 关键物种类别和丰度也高于高水位时期, 即低水位时期浮游细菌群落稳定性和浮游细菌网络复杂性更高.

(a)高水位时期浮游细菌共现网络和关键物种分析;(b)低水位时期浮游细菌共现网络和关键物种分析 图 6 高水位和低水位时期浮游细菌共现网络和关键物种 Fig. 6 Co-occurrence network and keystone taxa of bacterioplankton during different water level regulation periods

表 2 溪洛渡-向家坝梯级水库浮游细菌网络拓扑参数 Table 2 Topological parameters of bacterioplankton network in XLD-XJB cascade reservoirs

2.4 不同水库运行时期浮游细菌群落构建机制

微生物群落的组成受到确定性过程(同质和异质选择)和随机性过程(均质扩散、扩散限制和生态漂移)的控制.每个生态过程的相对贡献用零模型方法进行量化(图 7).高水位时期, 有22.22%的浮游细菌|βNTI|小于2, 而低水位时期有77.78%的浮游细菌群落|βNTI|小于2.因此, 高水位和低水位时期浮游细菌群落主要由随机性过程构建[图 7(a)].通过量化不同水位时期各个生态过程对浮游细菌群落构建贡献, 发现高水位时期浮游细菌群落构建主要受均质扩散(60.52%)、异质选择(31.23%)和生态漂变(8.06%)驱动, 而同质选择和扩散限制的作用不明显.低水位时期浮游细菌群落主要受到均质扩散(68.98%)、异质选择(22.22%)和漂变(8.64%)的驱动, 同质选择和扩散限制的影响较小[图 7(b)].

(a)高水位和低水位时期浮游细菌的βNTI, 绿色阴影部分表示-2<βNTI<2, 红色阴影部分表示βNTI>2和βNTI<-2;(b)高水位和低水位时期浮游细菌群落生态过程量化 图 7 高水位和低水位时期浮游细菌群群落构建过程 Fig. 7 Construction process of bacterioplankton communities during high and low water levels

3 讨论 3.1 不同时期浮游细菌群落的分布模式

大坝不同运行时期浮游细菌群落结构和多样性以及河流水环境存在显著差异[12, 38, 39].大坝蓄水和泄水不同运行方式差异导致高水位时期水力停留时间更长且水位更高[40].通过对高低水位时期浮游细菌群落组成进行分析, 发现在低水位时期蓝细菌门的相对丰度(5.43%)高于高水位时期(1.51%).蓝藻会为浮游细菌提供营养物质, 促进浮游细菌的生长[41].同时, 高水位时期水温显著低于低水位时期.高水位时期水温低, 浮游细菌的新陈代谢和生长速率减慢, 低水位时期水温高, 细菌新陈代谢旺盛生长速率加快[42, 43].因此在低水位时期浮游细菌Simpson指数和Chao指数高于高水位时期.高水位时期水力停留时间更长, 增加了微生物和营养物的滞留时间, 从而促进了营养物质的生物地球化学循环和物理转化[44~46], 使水体中的营养物质浓度降低, 导致高水位时期浮游细菌营养物DOC和TP浓度低于低水位时期.而低水位时期水温更高, 这可能是由于季节导致的差异.

本研究还观察到蓝细菌门、硝化螺旋菌门对高低水位两种水文环境敏感.蓝细菌门的相对丰度(高水位时期:1.51%;低水位时期:5.43%)在低水位时期更高, 硝化螺旋菌门相对丰度(高水位时期:1.67%;低水位时期:1.12%)在高水位时期更高.蓝细菌能进行产氧光合作用[47], 部分具有固氮能力, 对氮磷浓度变化有敏感的响应[48].硝化螺旋菌可将亚硝酸盐氧化成硝酸盐, 含有丰富的硝化功能[49], 是氮循环中的关键微生物群落[50].硝化螺旋菌对pH和养分水平的变化较为敏感[51].TP浓度在低水位时期更高, pH在高水位时期更高, 因此蓝细菌门和硝化螺旋菌门在高水位和低水位时期的差异可能是由变化的环境导致的, 并且这种环境变化会直接或间接[52, 53]地影响河流生态系统功能.但是, TN浓度在高、低水位时期无显著差异, 这可能是由于硝化螺旋菌可将亚硝酸盐氧化成硝酸盐[54], 但水体中TN的浓度不会改变, 而蓝细菌固氮释放到环境中的有机氮化合物被其他微生物利用[55].

3.2 不同时期浮游细菌群落生态过程

有研究发现, 低水位时期浮游细菌网络表现出更大的复杂性.这可能是由于营养物浓度会影响微生物间的相互作用[56, 57], 低水位时期DOC和TP浓度更高, 导致物种之间的相互作用强度增加[56], 使低水位时期浮游细菌网络更复杂.浮游细菌群落的稳定性可用模块化和关键物种来表征[33].结果表明低水位时期浮游细菌群落稳定性高于高水位时期.大型河流水库中的浮游细菌群落会直接或间接地受到蓄水和泄水的影响[2, 58], 水位变化和水力停留时间是其两个显著的特征.水位和水力停留时间的变化影响了水体的营养物浓度[3, 59]和水温[60~62]等物理指标.关键物种在微生物群落结构和功能中起着重要的作用[63, 64], 高、低水位时期浮游细菌网络关键物种差异显著, 没有相同的关键物种.之前的研究表明关键物种容易受到环境变化的影响[65, 66], 通过比较高低水位时期河流水环境理化指标, 发现高水位时期水体理化指标pH、水温、DOC和TP浓度与低水位时期差异显著, 并且理化指标与浮游细菌群落显著相关.因此环境差异导致高低水位时期关键物种差异显著.而微生物群落的模块化程度与水温和pH值显著相关[65], 即高低水位通过影响营养物浓度、水温等理化指标进而影响了浮游细菌群落的稳定性和物种间的相互作用.

从研究结果可以观察到高水位时期和低水位时期浮游细菌群落构建都以随机性过程为主导, 具体是由均质扩散驱动.河流具有高度的连通性和流动性[67~69], 导致微生物在水体中具有强的扩散能力[70]、诸多研究也表明河流生境具有同质化现象[71], 这些结果表明河流的流动性、连通性和大坝运行都会影响浮游细菌群落.此外, 研究发现在低水位时期均质扩散对浮游细菌群落构建的贡献增加, 异质选择贡献降低, 高水位时期均质扩散减少.这是因为水文条件的变化会影响微生物群落[72], 高水位时期水力停留时间较长, 阻碍了水体的流动[73], 导致浮游细菌扩散率降低, 因此均质扩散在高水位时期作用较弱.综上所述, 大坝通过泄水和蓄水影响了河流水文条件, 这样的运行模式导致不同水位时期两个浮游细菌群落构建过程的差异.

4 结论

(1)大坝不同运行时期浮游细菌群落α-多样性和河流理化指标差异显著, 大多数指标在低水位时期更高.

(2)蓝细菌门、硝化螺旋菌门是对高低水位敏感的物种, 参与了不同元素的循环过程.

(3)低水位时期浮游细菌网络稳定性和网络复杂性更高, 关键物种的种类更多.

(4)高、低水位时期浮游细菌群落构建均以随机性过程为主导, 低水位时期均质扩散对浮游细菌群落构建的贡献增加, 异质选择贡献降低.

参考文献
[1] Baxter A M, Johnson L, Edgerton J, et al. Structure and function of denitrifying bacterial assemblages in low-order Indiana Streams[J]. Freshwater Science, 2012, 31(2): 304-317. DOI:10.1899/11-066.1
[2] Ruiz-González C, Proia L, Ferrera I, et al. Effects of large river dam regulation on bacterioplankton community structure[J]. FEMS Microbiology Ecology, 2013, 84(2): 316-331. DOI:10.1111/1574-6941.12063
[3] Smith M A, Kominoski J S, Price R M, et al. Linking seasonal changes in organic matter composition and nutrients to shifting hydraulic gradients in coastal urban canals[J]. Water Resources Research, 2023, 59(2). DOI:10.1029/2022WR033334
[4] 刘红涛, 胡天龙, 王慧, 等. 典型水稻土细菌亚类群—泛化种、特化种的群落构建及功能潜力[J]. 土壤学报, 2023, 60(2): 546-557.
Liu H T, Hu T L, Wang H, et al. Community assembly and functional potential of habitat generalists and specialists in typical paddy soils[J]. Acta Pedologica Sinica, 2023, 60(2): 546-557.
[5] Reid A J, Carlson A K, Creed I F, et al. Emerging threats and persistent conservation challenges for freshwater biodiversity[J]. Biological Reviews, 2019, 94(3): 849-873. DOI:10.1111/brv.12480
[6] Vörösmarty C J, McIntyre P B, Gessner M O, et al. Global threats to human water security and river biodiversity[J]. Nature, 2010, 467(7315): 555-561. DOI:10.1038/nature09440
[7] Li F L, Guo F, Gao W, et al. Environmental DNA biomonitoring reveals the interactive effects of dams and nutrient enrichment on aquatic multitrophic communities[J]. Environmental Science & Technology, 2022, 56(23): 16952-16963.
[8] Wang Y, Lu L H, Hong Y G, et al. Divergent responses of taxonomic and predicted functional profiles of bacterioplankton to reservoir impoundment[J]. Environmental Research, 2020, 182. DOI:10.1016/j.envres.2019.109083
[9] 胡愈炘, 张静, 黄杰, 等. 长江流域河流和湖库的浮游细菌群落差异[J]. 环境科学, 2022, 43(3): 1414-1423.
Hu Y X, Zhang J, Huang J, et al. Characteristics of bacterioplankton community between river and lake/reservoir in the Yangtze River basin[J]. Environmental Science, 2022, 43(3): 1414-1423.
[10] 王洵, 张佳佳, 袁秋生, 等. 筑坝蓄水对不同水深浮游微生物群落结构和种间互作模式的影响[J]. 环境科学, 2022, 43(7): 3623-3634.
Wang X, Zhang J J, Yuan Q S, et al. Effects of damming and impoundment on planktonic microbial community structure and interspecific interaction patterns in different water depths[J]. Environmental Science, 2022, 43(7): 3623-3634.
[11] 王雪, 刘晋仙, 柴宝峰, 等. 宁武亚高山湖泊细菌群落的时空格局及驱动机制[J]. 环境科学, 2019, 40(7): 3285-3294.
Wang X, Liu J X, Chai B F, et al. Spatio-temporal patterns of microbial communities and their driving mechanisms in subalpine lakes, Ningwu, Shanxi[J]. Environmental Science, 2019, 40(7): 3285-3294.
[12] Geng M M, Wang K L, Yang N, et al. Spatiotemporal water quality variations and their relationship with hydrological conditions in Dongting Lake after the operation of the Three Gorges Dam, China[J]. Journal of Cleaner Production, 2021, 283. DOI:10.1016/j.jclepro.2020.124644
[13] 万凌琳, 陈芷凡, 郭佳, 等. 生物共现网络原理及其在淡水生态系统评估中的应用[J]. 湖泊科学, 2022, 34(6): 1765-1787.
Wan L L, Chen Z F, Guo J, et al. Principle and application of co-occurrence networks for freshwater ecosystem assessment[J]. Journal of Lake Sciences, 2022, 34(6): 1765-1787.
[14] Rosindell J, Hubbell S P, Etienne R S. The unified neutral theory of biodiversity and biogeography at age ten[J]. Trends in Ecology & Evolution, 2011, 26(7): 340-348.
[15] Zhou J Z, Ning D L. Stochastic community assembly: does it matter in microbial ecology?[J]. Microbiology and Molecular Biology Reviews, 2017, 81(4). DOI:10.1128/mmbr.00002-17
[16] Wang B L, Liu N, Yang M L, et al. Co-occurrence of planktonic bacteria and archaea affects their biogeographic patterns in China's coastal wetlands[J]. Environmental Microbiome, 2021, 16(1). DOI:10.1186/S40793-021-00388-9
[17] Shapiro R S, Cowen L E. Thermal control of microbial development and virulence: Molecular mechanisms of microbial temperature sensing[J]. mBio, 2012, 3(5). DOI:10.1128/mBio.00238-12
[18] Luo Z M, Liu J X, Zhao P Y, et al. Biogeographic patterns and assembly mechanisms of bacterial communities differ between habitat generalists and specialists across elevational gradients[J]. Frontiers in Microbiology, 2019, 10. DOI:10.3389/fmicb.2019.00169
[19] Yuan M M, Guo X, Wu L W, et al. Climate warming enhances microbial network complexity and stability[J]. Nature Climate Change, 2021, 11(4): 343-348. DOI:10.1038/s41558-021-00989-9
[20] Mo Y Y, Zhang W J, Wilkinson D M, et al. Biogeography and co-occurrence patterns of bacterial generalists and specialists in three subtropical marine bays[J]. Limnology and Oceanography, 2021, 66(3): 793-806. DOI:10.1002/lno.11643
[21] Liu S, Ren H X, Shen L D, et al. pH levels drive bacterial community structure in sediments of the Qiantang River as determined by 454 pyrosequencing[J]. Frontiers in Microbiology, 2015, 6. DOI:10.3389/fmicb.2015.00285
[22] 程豹, 望雪, 徐雅倩, 等. 澜沧江流域浮游细菌群落结构特征及驱动因子分析[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.
[23] Hu A Y, Yang X Y, Chen N W, et al. Response of bacterial communities to environmental changes in a mesoscale subtropical watershed, Southeast China[J]. Science of the Total Environment, 2014, 472: 746-756. DOI:10.1016/j.scitotenv.2013.11.097
[24] Prach K, Walker L R. Four opportunities for studies of ecological succession[J]. Trends in Ecology & Evolution, 2011, 26(3): 119-123.
[25] 侯轶群, 邹曦, 陈小娟, 等. 金沙江下游大坝建设对纵向连通性的影响研究[J]. 水生态学杂志, 2024, 45(1): 26-31.
Hou Y Q, Zou X, Chen X J, et al. Impact of dam construction on the longitudinal connectivity of the lower Jinsha River[J]. Journal of Hydroecology, 2024, 45(1): 26-31.
[26] 刘恒, 刘九夫, 唐海行. 澜沧江流域(云南段)水资源开发利用现状及趋势分析[J]. 水科学进展, 1998, 9(1): 70-76.
Liu H, Liu J F, Tang H X. Present and future of water resources development in Lancang River basin in Yunnan Province[J]. Advances in Water Science, 1998, 9(1): 70-76. DOI:10.3321/j.issn:1001-6791.1998.01.012
[27] 纪道斌, 龙良红, 徐慧, 等. 梯级水库建设对水环境的累积影响研究进展[J]. 水利水电科技进展, 2017, 37(3): 7-14.
Ji D B, Long L H, Xu H, et al. Advances in study on cumulative effects of construction of cascaded reservoirs on water environment[J]. Advances in Science and Technology of Water Resources, 2017, 37(3): 7-14.
[28] Tang X M, Xie G J, Shao K Q, et al. Contrast diversity patterns and processes of microbial community assembly in a river-lake continuum across a catchment scale in northwestern China[J]. Environmental Microbiome, 2020, 15(1). DOI:10.1186/s40793-020-00356-9
[29] Mo Y Y, Peng F, Gao X F, et al. Low shifts in salinity determined assembly processes and network stability of microeukaryotic plankton communities in a subtropical urban reservoir[J]. Microbiome, 2021, 9(1). DOI:10.1186/s40168-021-01079-w
[30] Kallmeyer J, Pockalny R, Adhikari R R, et al. Global distribution of microbial abundance and biomass in subseafloor sediment[J]. Proceedings of the National Academy of Sciences of the United States of America, 2012, 109(40): 16213-16216.
[31] Ma B, Wang H Z, Dsouza M, et al. Geographic patterns of co-occurrence network topological features for soil microbiota at continental scale in eastern China[J]. The ISME Journal, 2016, 10(8): 1891-1901. DOI:10.1038/ismej.2015.261
[32] Barberán A, Bates S T, Casamayor E O, et al. Using network analysis to explore co-occurrence patterns in soil microbial communities[J]. The ISME Journal, 2012, 6(2): 343-351. DOI:10.1038/ismej.2011.119
[33] Zhao D Y, Shen F, Zeng J, et al. Network analysis reveals seasonal variation of co-occurrence correlations between Cyanobacteria and other bacterioplankton[J]. Science of the Total Environment, 2016, 573: 817-825. DOI:10.1016/j.scitotenv.2016.08.150
[34] Amit G, Bashan A. Top-down identification of keystone taxa in the microbiome[J]. Nature Communications, 2023, 14(1). DOI:10.1038/s41467-023-39459-5
[35] Zhou J Z, Deng Y, Luo F, et al. Functional molecular ecological networks[J]. mBio, 2010, 1(4). DOI:10.1128/mBio.00169-10
[36] Stegen J C, Lin X J, Fredrickson J K, et al. Quantifying community assembly processes and identifying features that impose them[J]. The ISME Journal, 2013, 7(11): 2069-2079. DOI:10.1038/ismej.2013.93
[37] Stegen J C, Lin X J, Konopka A E, et al. Stochastic and deterministic assembly processes in subsurface microbial communities[J]. The ISME Journal, 2012, 6(9): 1653-1664. DOI:10.1038/ismej.2012.22
[38] 丁宁, 杨莹莹, 万年新, 等. 调蓄水库细菌群落季节变化特征与影响因素[J]. 环境科学, 2023, 44(3): 1484-1496.
Ding N, Yang Y Y, Wan N X, et al. Seasonal variation and influencing factors of bacterial communities in storage reservoirs[J]. Environmental Science, 2023, 44(3): 1484-1496.
[39] Spinti R A, Condon L E, Zhang J. The evolution of dam induced river fragmentation in the United States[J]. Nature Communications, 2023, 14(1). DOI:10.1038/s41467-023-39194-x
[40] 冉祥滨, 于志刚, 姚庆祯, 等. 水库对河流营养盐滞留效应研究进展[J]. 湖泊科学, 2009, 21(5): 614-622.
Ran X B, Yu Z G, Yao Q Z, et al. Advances in nutrient retention of dams on river[J]. Journal of Lake Sciences, 2009, 21(5): 614-622. DOI:10.3321/j.issn:1003-5427.2009.05.002
[41] 陈姝, 郭照冰, 方华, 等. CO2浓度升高对湖泊浮游藻类与浮游细菌耦合关系的影响[J]. 中国环境科学, 2015, 35(7): 2209-2216.
Chen S, Guo Z B, Fang H, et al. Effects of elevated CO2 concentration on the coupling relationship between planktonic algae-planktonic bacteria[J]. China Environmental Science, 2015, 35(7): 2209-2216. DOI:10.3969/j.issn.1000-6923.2015.07.049
[42] 张雅洁, 李珂, 朱浩然, 等. 北海湖微生物群落结构随季节变化特征[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.
[43] 沈烽, 赵大勇, 黄睿, 等. 南京玄武湖浮游细菌群落结构的季节变化及其与环境因子的关系[J]. 湖泊科学, 2017, 29(3): 662-669.
Shen F, Zhao D Y, Huang R, et al. Seasonal variation of bacterioplankton community structure in Xuanwu Lake (Nanjing) and its relationship with environmental factors[J]. Journal of Lake Sciences, 2017, 29(3): 662-669.
[44] Fearnside P M. Hydroelectric dams in the Brazilian Amazon as sources of 'greenhouse' gases[J]. Environmental Conservation, 1995, 22(1): 7-19. DOI:10.1017/S0376892900034020
[45] Friedl G, Wüest A. Disrupting biogeochemical cycles-Consequences of damming[J]. Aquatic Sciences, 2002, 64(1): 55-65. DOI:10.1007/s00027-002-8054-0
[46] Van Cappellen P, Maavara T. Rivers in the Anthropocene: global scale modifications of riverine nutrient fluxes by damming[J]. Ecohydrology & Hydrobiology, 2016, 16(2): 106-111.
[47] 王梦梦, 龚正文, 陈成, 等. 湖泊蓝藻多聚磷酸盐代谢机制及功能[J]. 湖泊科学, 2023, 35(3): 780-794.
Wang M M, Gong Z W, Chen C, et al. Metabolic mechanism and function of cyanobacterial polyphosphate in lakes[J]. Journal of Lake Sciences, 2023, 35(3): 780-794.
[48] Chen M Y, Teng W K, Zhao L, et al. Phylogenomics uncovers evolutionary trajectory of nitrogen fixation in Cyanobacteria[J]. Molecular Biology and Evolution, 2022, 39(9). DOI:10.1093/molbev/msac171
[49] Miao Y, Liao R H, Zhang X X, et al. Metagenomic insights into Cr(VI) effect on microbial communities and functional genes of an expanded granular sludge bed reactor treating high-nitrate wastewater[J]. Water Research, 2015, 76: 43-52. DOI:10.1016/j.watres.2015.02.042
[50] Daims H, Lebedeva E V, Pjevac P, et al. Complete nitrification by Nitrospira bacteria[J]. Nature, 2015, 528(7583): 504-509. DOI:10.1038/nature16461
[51] Li Y T, Gao Y, Chen W, et al. Shifts in bacterial diversity, interactions and microbial elemental cycling genes under cadmium contamination in paddy soil: Implications for altered ecological function[J]. Journal of Hazardous Materials, 2024, 461. DOI:10.1016/j.jhazmat.2023.132544
[52] Giling D P, Beaumelle L, Phillips H R P, et al. A niche for ecosystem multifunctionality in global change research[J]. Global Change Biology, 2019, 25(3): 763-774. DOI:10.1111/gcb.14528
[53] Hautier Y, Tilman D, Isbell F, et al. Anthropogenic environmental changes affect ecosystem stability via biodiversity[J]. Science, 2015, 348(6232): 336-340. DOI:10.1126/science.aaa1788
[54] Li Y, Gao L, Niu L H, et al. Developing a statistical-weighted index of biotic integrity for large-river ecological evaluations[J]. Journal of Environmental Management, 2021, 277. DOI:10.1016/j.jenvman.2020.111382
[55] 谢平. 论蓝藻水华的发生机制: 从生物进化、生物地球化学和生态学视点[M]. 北京: 科学出版社, 2007.
[56] Ratzke C, Barrere J, Gore J. Strength of species interactions determines biodiversity and stability in microbial communities[J]. Nature Ecology & Evolution, 2020, 4(3): 376-383.
[57] Hu J L, Amor D R, Barbier M, et al. Emergent phases of ecological diversity and dynamics mapped in microcosms[J]. Science, 2022, 378(6615): 85-89. DOI:10.1126/science.abm7841
[58] Li Z, Lu L H, Guo J S, et al. Responses of spatial-temporal dynamics of bacterioplankton community to large-scale reservoir operation: a case study in the Three Gorges Reservoir, China[J]. Scientific Reports, 2017, 7. DOI:10.1038/srep42469
[59] 刘霞, 刘宝贵, 陈宇炜, 等. 鄱阳湖浮游植物叶绿素a及营养盐浓度对水位波动的响应[J]. 环境科学, 2016, 37(6): 2141-2148.
Liu X, Liu B G, Chen Y W, et al. Responses of nutrients and chlorophyll a to water level fluctuations in Poyang Lake[J]. Environmental Science, 2016, 37(6): 2141-2148.
[60] Olsson F, Mackay E, Barker P, et al. Can reductions in water residence time be used to disrupt seasonal stratification and control internal loading in a eutrophic monomictic lake?[J]. Journal of Environmental Management, 2021, 304. DOI:10.1016/j.jenvman.2021.114169
[61] Wang S Y, Gao Y, Jia J J, et al. Water level as the key controlling regulator associated with nutrient and gross primary productivity changes in a large floodplain-lake system (Lake Poyang), China[J]. Journal of Hydrology, 2021, 599. DOI:10.1016/j.jhydrol.2021.126414
[62] Guarch-Ribot A, Butturini A. Hydrological conditions regulate dissolved organic matter quality in an intermittent headwater stream. From drought to storm analysis[J]. Science of the Total Environment, 2016, 571: 1358-1369. DOI:10.1016/j.scitotenv.2016.07.060
[63] Banerjee S, Schlaeppi K, van der Heijden M G A. Keystone taxa as drivers of microbiome structure and functioning[J]. Nature Reviews Microbiology, 2018, 16(9): 567-576. DOI:10.1038/s41579-018-0024-1
[64] Xue Y Y, Chen H H, Yang J R, et al. Distinct patterns and processes of abundant and rare eukaryotic plankton communities following a reservoir cyanobacterial bloom[J]. The ISME Journal, 2018, 12(9): 2263-2277. DOI:10.1038/s41396-018-0159-0
[65] Liu S W, Yu H, Yu Y H, et al. Ecological stability of microbial communities in Lake Donghu regulated by keystone taxa[J]. Ecological Indicators, 2022, 136. DOI:10.1016/j.ecolind.2022.108695
[66] D'Amen M, Mod H K, Gotelli N J, et al. Disentangling biotic interactions, environmental filters, and dispersal limitation as drivers of species co-occurrence[J]. Ecography, 2018, 41(8): 1233-1244. DOI:10.1111/ecog.03148
[67] Isabwe A, Yang J R, Wang Y M, et al. Community assembly processes underlying phytoplankton and bacterioplankton across a hydrologic change in a human-impacted river[J]. Science of the Total Environment, 2018, 630: 658-667. DOI:10.1016/j.scitotenv.2018.02.210
[68] Wang Y M, Liu L M, Chen H H, et al. Spatiotemporal dynamics and determinants of planktonic bacterial and microeukaryotic communities in a Chinese subtropical river[J]. Applied Microbiology and Biotechnology, 2015, 99(21): 9255-9266. DOI:10.1007/s00253-015-6773-0
[69] Grill G, Lehner B, Thieme M, et al. Author correction: mapping the world's free-flowing rivers[J]. Nature, 2019, 572(7768). DOI:10.1038/s41586-019-1379-9
[70] Chen W D, Ren K X, Isabwe A, et al. Stochastic processes shape microeukaryotic community assembly in a subtropical river across wet and dry seasons[J]. Microbiome, 2019, 7(1). DOI:10.1186/s40168-019-0749-8
[71] Logares R, Lindström E S, Langenheder S, et al. Biogeography of bacterial communities exposed to progressive long-term environmental change[J]. The ISME Journal, 2013, 7(5): 937-948.
[72] Niño-García J P, Ruiz-González C, del Giorgio P A. Interactions between hydrology and water chemistry shape bacterioplankton biogeography across boreal freshwater networks[J]. The ISME Journal, 2016, 10(7): 1755-1766.
[73] Spaak J W, Baert J M, Baird D J, et al. Shifts of community composition and population density substantially affect ecosystem function despite invariant richness[J]. Ecology Letters, 2017, 20(10): 1315-1324.