环境科学  2023, Vol. 44 Issue (2): 847-856   PDF    
西安市景观湖泊水体细菌群落结构分析与代谢功能预测
王敏1,2, 张雨桐1,2, 黄晨1,2, 任杰辉1,2, 万甜1,2, 程文1,2     
1. 西安理工大学省部共建西北旱区生态水利国家重点实验室, 西安 710048;
2. 西安理工大学水利水电学院, 西安 710048
摘要: 湖泊水体中细菌群落多样性及其代谢功能研究对城市景观湖泊水环境健康发展有着重要意义.以汉城湖、南湖、兴庆湖和桃花潭为研究区域,使用Illumina Nova高通量测序技术研究4个不同类型湖泊中群落组成和结构,使用冗余分析方法分析细菌群落与环境因子的相关性,采用Tax4Fun预测细菌群落的代谢功能.结果表明,不同湖泊水体细菌群落多样性各不相同,其中桃花潭α多样性最高,兴庆湖最低,同一湖泊水体的细菌群落结构有聚集趋势;4个湖泊水体的细菌优势门类相似,主要是变形菌门、拟杆菌门和放线菌门,不同之处仅在于细菌的比例;4个湖泊水体属水平细菌群落分布有较大差异;水体理化性质与细菌群落存在显著相关性,其中最关键因子为氮磷营养盐(P < 0.05);Tax4Fun功能预测表明,4个湖泊水体细菌群落代谢功能相似,氨基酸代谢和碳水化合物代谢相对丰度最高,此外,其潜在的耐药性细菌污染较大以及致人类患传染性疾病可能性较大.研究结果有助于了解西安市城市湖泊水环境细菌现状,为城市湖泊治理及可持续发展提供理论依据.
关键词: 城市湖泊      细菌群落      多样性      环境因子      代谢功能     
Bacterial Community Structure and the Prediction of Metabolic Function in Landscape Lake Water in Xi'an
WANG Min1,2 , ZHANG Yu-tong1,2 , HUANG Chen1,2 , REN Jie-hui1,2 , WAN Tian1,2 , CHENG Wen1,2     
1. State Key Laboratory of Eco-hydraulics in Northwest Arid Region, Xi'an University of Technology, Xi'an 710048, China;
2. Institute of Water Resources and Hydro-Electric Engineering, Xi'an University of Technology, Xi'an 710048, China
Abstract: The study of bacterial community diversity and metabolic function in lake water is of great significance to the healthy development of the urban landscape lake water environment. In this study, Hancheng Lake, Nanhu Lake, Xingqing Lake, and Taohuatan Lake were selected as the study area. Illumina Nova high-throughput sequencing technology was used to study the composition and structure of bacterial communities in four different types of lakes. The correlation between bacterial communities and environmental factors was analyzed using the RDA method. Tax4Fun was used to predict the metabolic function of bacterial communities. The results showed that the diversity of bacterial communities was different in different lakes. The α diversity of Taohuatan Lake was the highest, and that of Xingqing Lake was the lowest. The bacterial community structure of the same lake had a trend of aggregation. The dominant phyla of bacteria in the four lakes were similar, mainly Proteobacteria, Bacteroidetes, and Actinobacteria. The difference was only in the proportion of bacteria. Additionally, the distribution of bacterial communities at the genus level in the four lakes was quite different. The physical and chemical properties of lakes were significantly correlated with bacterial community, and the most critical factor was nitrogen and phosphorus nutrients(P < 0.05). Tax4Fun function prediction showed that the metabolic function of bacterial communities in the four lakes was similar, and the relative abundance of amino acid metabolism and carbohydrate metabolism was the highest. In addition, the potential drug-resistant bacterial pollution and the possibility of human infectious diseases were higher. These results are helpful for understanding the current situation of bacteria in urban lakes in Xi'an and providing a theoretical basis for the management and sustainable development of urban lakes.
Key words: urban lake      bacterial community      diversity      environmental factors      metabolic function     

细菌是水生生态系统的重要组成部分, 在生物地球化学循环中起着不可替代的作用[1].比如, 细菌通过利用溶解的有机化合物并将其转化为对自养生物有用的无机化学物质来维持水生生态系统的生产力[2].细菌通过同化和异化等代谢过程完成碳、氮和磷等元素的降解, 也可以促进有机物的降解和污染物的转化, 从而维持水生生态系统的平衡[3].例如, 蓝藻门会导致水体富营养化严重, 被污染水体可能出现食物链污染, 对生态系统造成重大的环境威胁[4]; 当在废水和废弃物厌氧处理系统中起重要作用的厚壁菌门存在时, 说明水体已经面临污染胁迫[5], 与氮和甲烷代谢相关的浮游细菌在水生生态系统中发挥着不成比例的重要作用, 尽管它们在细菌群落中的比例很小[6].此前, 对细菌的研究多集中在群落结构和分布特征上, 对其代谢功能研究展开较少.因此, 本研究系统分析水生生态系统中细菌群落的结构、分布和其代谢功能特征, 对维护水体环境有着极为重要的意义.

西安市作为北方缺水城市, 城市湖泊补水较少, 水体交换速率慢, 且受周边居民生产生活的影响, 城市湖泊水体环境较差[7].张锋等[8]对汉城湖的水环境进行了全面调查, 发现部分水质指标已经超出景观娱乐用水水质A类标准的最低要求.城市景观水体大多处于静止、封闭或半封闭状态, 自净能力差, 水生生物单一, 且生态系统不稳定[9].

因此, 本文以西安市景观湖泊(汉城湖、南湖、兴庆湖和桃花潭)为研究对象, 对4个湖泊进行实地调查以及理化指标测定; 并使用Illumina Nova高通量测序技术获得湖泊水体细菌群落结构特征, 分析水体中细菌群落分布规律; 应用RDA分析得出细菌群落与环境因子的相关性; 利用Tax4Fun预测细菌群落代谢功能, 以期为城市景观湖泊水环境健康发展提供依据.

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

本研究涉及西安市区4个典型的城市湖泊, 分别位于西安市主城区的东、西、南、北这4个方向, 每个湖泊周围环境有明显差异, 可以代表西安市不同类型城市湖泊水体的特征.

南湖位于西安市东南郊, 来水水源为曲江水厂第一道水, 全长1 088 m, 宽窄不等, 最宽处552 m, 平均水深1.68 m, 水力停留时间约38 d.南湖因地势开凿, 湖形曲折, 南北长, 东西窄[10].

汉城湖位于西安市西北方向, 补水来自沣河, 湖体自南向北流入渭河, 全长6 270 m, 宽20~100 m不等, 水力停留时间约14 d.在2009年前, 汉城湖是城市污水主要排污口, 城市河道改造后, 承接兴庆湖、护城河、老城区和西北郊区部分城市排污、雨洪排泄任务[11].

桃花潭位于浐河下游, 补水来自临近河流和西安市第三污水处理厂出水, 全长1 200 m, 宽约660 m, 平均水深1.3 m[12].

兴庆湖位于西安市城区东南角, 补水为南郊水厂和雨水, 面积约1.1×105 m2, 平均水深1.6m, 水力停留时间约310 d, 该湖是西安市东郊地区的防洪缓冲库[13].

1.2 样品采集及处理

分别在西安市汉城湖、兴庆湖、南湖、桃花潭的进水口、湖中心和出水口进行采样点布设, 具体位置如图 1所示.

图 1 采样点布设示意 Fig. 1 Schematic diagram of sampling point layout

水样于2020年8月采集, 使用简易采水器在水面以下0.5 m处取3个水样, 各1 L.装入无菌聚乙烯瓶内, 临时保存在4℃以下的保温箱, 直至运送到实验室.在超净工作台, 使用0.22 μm醋酸纤维滤膜过滤截留水样中的细菌, 将滤膜保存在-80℃以下, 直至进行细菌群落测序.

1.3 水体理化指标测定

对南湖(N11~13)、兴庆湖(N21~23)、桃花潭(N31~33)和汉城湖(N41~43)共12个采样点的水体理化性质进行测定, 测定指标及检测方法, 如表 1所示.

表 1 水质指标及其测定方法 Table 1 Water quality index and its determination method

1.4 DNA提取

采用标准化CTAB法提取样本DNA.即将样品加入含有溶菌酶和1 000 μL CTAB缓冲液的裂解液中.用酚(pH8.0)∶氯仿∶异戊醇(25∶24∶1)提取DNA, 在-20℃条件下用异丙醇沉淀, 然后经75%乙醇洗涤2次, 用ddH2O溶解.最后加入1 μL RNase A在37℃下消化RNA 15 min, 用1%琼脂糖凝胶电泳检测DNA浓度和纯度, 根据DNA浓度, 无菌水稀释到1 ng·μL-1[14].

1.5 PCR扩增

16S V4区扩增采用带Barcode的特异引物, 带正向引物515F(5′-GTGCCAGCMGCCGCGGTAA-3′)和反向引物806 R (5′-GGACTACHVGGGTWTCTA AT-3′)的通用引物[15].PCR反应采用Phusion® High-Fidelity PCR Master Mix with GC Buffer(New England Biolabs, Ipswich, MA, USA), 总体积为30 μL(Phusion Master Mix 15 μL, 6μmol·L-1引物, 10 μL gDNA, 2 μL H2O).

PCR反应程序为:在98℃进行1 min, 随后30个循环(98℃, 10 s; 50℃, 30 s; 72℃, 30 s), 然后在72℃下进行最后延伸步骤5 min.PCR产物使用2%浓度的琼脂糖凝胶进行电泳检测, 根据PCR产物浓度进行等量混样, 充分混匀后使用2%的琼脂糖凝胶电泳检测PCR产物, 对目的条带使用Qiagen公司提供的胶回收试剂盒回收产物.

1.6 测序数据处理

使用Illumina公司TruSeq DNA PCR-Free Library Preparation Kit建库试剂盒进行文库的构建, 构建好的文库经过Qubit定量和文库检测, 合格后, 使用NovaSeq 6000进行上机测序.

使用FLASH将读取的成对的数据合并, 根据特异性Barcode分配给每个样品, 利用QIIME软件包对FLASH生成的原始标签进行筛选和分析.

根据Barcode序列和PCR扩增引物序列从下机数据中拆分出各样本数据, 截去Barcode和引物序列后使用FLASH(V1.2.7, http://ccb.jhu.edu/software/FLASH/), 拼接得到的原始Tags数据(raw tags)[16].利用Qiime(V1.9.1, http://qiime.org/scripts/split_libraries_fastq.html)[17]对FLASH生成的数据进行筛选和分析. 经以上处理, Tags序列通过(https://github.com/torognes/vsearch/)[18]进行去除嵌合体序列的处理, 得到最终的有效数据(effective tags).

1.7 Tax4Fun功能预测分析

使用Tax4Fun预测了浮游细菌群落的潜在功能概况, Tax4Fun是一个开源的R包, 它基于16S rRNA标记基因概况预测微生物群落的功能能力.根据SILVA 123数据库使用封闭参考获得OTU表, 使用京都基因和基因组百科全书(KEGG)途径参考使用标准化OTU表计算代谢能力[19].

1.8 数据分析

利用Origin 2018软件对水体理化性质、细菌丰度进行分析处理及绘图; 利用Canoco5.0软件对细菌群落物种信息进行去趋势对应分析(detrended correspondence analysis, DCA), 后基于冗余分析(redundency analysis, RDA)探究细菌群落与环境因子间的响应关系; 最后利用Tax4Fun进行细菌群落代谢功能预测.

2 结果与讨论 2.1 西安市4个湖泊水质现状

图 2为西安市4个湖泊基本水质指标分布状况, 水体pH在6.98~7.61之间, 水体基本呈中性或弱碱性.溶解氧在5.26~10.69 mg·L-1之间, 其中汉城湖溶解氧最高, 桃花潭溶解氧最低.桃花潭COD、TP、TN和NH4+-N均为最高, 分别达到地表劣Ⅴ类、Ⅴ类、劣Ⅴ类和Ⅲ类水质标准; 南湖、兴庆湖和汉城湖COD达到Ⅳ类, TP达到Ⅳ类, NH4+-N达到Ⅱ类, 南湖和兴庆湖TN达到Ⅴ类, 汉城湖仅达劣Ⅴ类.总体来说, 4个湖泊均存在氮营养盐污染, 其中汉城湖和桃花潭氮超标严重.

图 2 4个湖泊水质现状 Fig. 2 Status of water quality in the four lakes

2.2 细菌群落多样性分析

α多样性用于分析样本内的细菌群落多样性, 通过单样本的多样性分析可以反映样本内的细菌群落的丰富度和多样性, 并用一系列统计学指数分析评估各样本中细菌群落的物种丰富度和多样性的差异[20].选用Simpson、Chao1、ACE和Shannon等指数来反映细菌群落的相对丰度.Chao1指数和ACE指数用来表征物种的丰富度[21], Chao1指数和ACE指数数值越大, 表明群落丰富度越高; Simpson指数和Shannon指数用来表征物种的多样性[22], Simpson指数和Shannon指数数值越大, 表明群落多样性越高.

图 3所示, 西安市4个湖泊水体Chao1指数平均值分别为:南湖(1 019.25)、兴庆湖(793.16)、桃花潭(1 030.92)和汉城湖(932.25); ACE指数平均值分别为:南湖(1 038.66)、兴庆湖(817.50)、桃花潭(1 072.22)和汉城湖(939.77); Shannon指数平均数分别为:南湖(5.99)、兴庆湖(5.34)、桃花潭(6.84)和汉城湖(6.44); Simpson指数平均数分别为:南湖(0.95)、兴庆湖(0.94)、桃花潭(0.98)和汉城湖(0.97).结果表明, 西安市4个湖泊水体中细菌群落丰富度高低为:桃花潭>南湖>汉城湖>兴庆湖; 细菌群落多样性高低为:桃花潭>汉城湖>南湖>兴庆湖.综合α多样性各项指标, 桃花潭均高于其他湖泊, 说明桃花潭水体细菌多样性较高, 水体环境适合细菌的生长繁殖[23], 这可能与桃花潭存在一定程度富营养化有关.

图 3 湖泊水体细菌多样性指数 Fig. 3 Bacterial diversity index of lake water

2.3 湖泊细菌群落结构的空间分布规律 2.3.1 门水平细菌群落结构分布

图 4为西安市4个城市湖泊水样中细菌在门水平的相对丰度分布结果.每个样本选取在门水平上最大丰度前10的物种, 生成物种相对丰度柱形累加图进行分析.可以看出, 湖泊水体中变形菌门(Proteobacteria)相对丰度(18.3%~73.8%)显著高于其他门类, 在地表水细菌群落中是常见的优势菌群.其他主要菌群门类依次为:放线菌门(Actinobacteria, 2.7%~49.7%)、拟杆菌门(Bacteroidetes, 3.0%~24.8%)、蓝藻门(Cyanobacteria, 0.9%~16.3%)、疣微菌门(Verrucomicrobia, 0.9%~12.7%)、厚壁菌门(Firmicutes, 0.8%~7.1%)和绿弯菌门(Chloroflexi, 0.8%~4.4%), 这与其他水环境中细菌群落组成基本相似[24].在4个湖泊中相对丰度均较少的有:酸杆菌门(Acidobacteria, 0.1%~1.5%)、Parcubacteria菌门(0%~1.2%)和单胞菌门(Gemmatimonadetes, 0.1%~1.1%).

图 4 门水平细菌结构分布 Fig. 4 Bacterial structure distribution at the phylum level

在门分类水平上, 4个湖泊水体的细菌群落分布有所差异.兴庆湖中优势菌群(平均丰度>5%)为变形菌门和拟杆菌门; 汉城湖中优势菌群为放线菌门、变形菌门、拟杆菌门、蓝藻门和疣微菌门; 南湖中优势菌群为变形菌门、放线菌门和蓝藻门; 桃花潭中优势菌群为放线菌门、变形菌门、蓝藻门和拟杆菌门.其中, 汉城湖细菌丰富度最高, 包含细菌种群数量最多.从各湖泊3个采样点来看, 湖泊进水、湖心和出水中细菌群落分布和相对丰度变化不大.

水体中的细菌对水质有明显的指示意义[25].变形菌门在有机物分解和循环中有重要作用, 可以去除水体中的磷, 并将亚硝酸盐转化成硝酸盐, 其高丰度与其高营养可用性有关[26]; 拟杆菌门与许多哺乳动物的肠道微生物群相关, 其中一些物种被认为是有效的粪便指示剂, 其丰度较高可能由于住宅区废水直接排放造成, 这些废水可能含有硝酸盐、氨和粪便等污染物[27]; 放线菌门是与水体净化过程相关的重要细菌, 可以产生多种抗生素、细菌素和酶, 对水环境有消毒作用, 可抑制并杀灭病原细菌[28]; 蓝藻门是水体富营养化主要指示生物, 相对较高的蓝藻门的分布, 表明湖泊存在潜在的富营养化的风险[29]; 疣微菌门在富营养环境下可利用硫酸盐降解多糖, 在寡营养条件下可以通过化能自养代谢获取能量[30].这些细菌的分布主要与湖泊水体水质和环境状况有关.

2.3.2 属水平细菌群落结构分布

在属分类水平上, 不同湖泊水体中的优势菌属及其数量存在一定差异, 属水平细菌相对丰度如图 5所示.南湖的优势菌属(平均丰度>5%)为:unidentified_α-变形菌属(unidentified_α-Proteobacteria, 15.9%~20.9%)和黑草属(Buchnera, 3.6%~22.6%), 兴庆湖的优势菌属为:Pseudarcicella(18.0%~20.0%)、栖湖菌属(Limnohabitans, 10.4%~12.0%)、多核杆菌属(Polynucleobacter, 8.9%~10.6%)、红育菌属(Rhodoferax, 7.6%~7.8%)、Malikia(3.7%~7.6%)和Pseudorhodobacter(3.8%~6.6%), 汉城湖的优势菌属为:unidentified_微酸菌属(unidentified_Acidimicrobiia, 2.0%~8.7%), 桃花潭的优势菌属为:多核杆菌属(Polynucleobacter, 3.7%~5.3%).

图 5 属水平细菌结构分布 Fig. 5 Bacterial structure distribution at the genus level

4个湖泊水体中优势属主要属于变形菌门、放线菌门和拟杆菌门, 这与门水平优势菌分析结果一致.在这些属水平的细菌中, 属于γ-变形菌纲的布赫纳氏菌属、LimnohabitansMalikia、多核杆菌属和红育菌属分别在固氮[31]、硝酸盐还原[32]、反硝化作用[33]、溶解性有机物降解[34]和硝化作用[35]等氮的循环过程中发挥作用; 属于α-变形菌纲的unidentified_α-变形菌属能够利用有机物去除污染[36], 桃花潭中α-变形菌属相对含量较低, 说明桃花潭中能起到脱氮除磷效果的菌群较少, 而通过对比分析水质检测结果, 发现桃花潭中氮磷化合物浓度高于其他采样点; 属于微酸菌纲的unidentified_微酸菌属是水质净化重要细菌, 其在汉城湖、南湖和桃花潭中较高的相对丰度表明其水环境生态修复及净化效果较好; 属于拟杆菌纲的Pseudarcicella为指示水体所受的粪便污染提供了途径, 其在兴庆湖中较高的相对丰度表明兴庆湖粪源污染较严重, 可能存在未经处理的生活废水排入[37].

2.4 环境因子对细菌群落差异的影响

由于细菌对水环境的变化十分敏感, 细菌功能和多样性与环境因子之间的关系对于促进水环境的健康至关重要[38].本研究选取6个湖泊水体基本理化性质指标(pH、DO、COD、NH4+-N、TN和TP)作为环境因子, 运用线性模型RDA分析属水平前8种细菌群落与环境因子间的响应关系, 结果如图 6所示.检验结果显示, TN(r2=0.765, P=0.002)、TP(r2=0.571, P=0.004)、NH4+-N(r2=0.540, P=0.004)和pH(r2=0.316, P=0.046)对湖泊水体细菌群落的影响达到显著水平(P < 0.05).以上相关性分析进一步解析了西安市4个城市湖泊水体细菌群落结构与环境因子的响应关系.

1.pH, 2.COD, 3.TP, 4.NH4+-N, 5.DO, 6.TN; A. Buchnera, B. unidentified_α-Proteobacteria, C. Pseudarcicella, D. Limbohabitans, E. Polynucleobacter, F. unidentified_Acidimicrobiia, G. Rhodoferax, H. Malikia; 红色剪头表示环境因子, 蓝色箭头表示细菌群落 图 6 水体细菌群落与关键环境因子的RDA分析 Fig. 6 Redundancy analysis of bacterial community and key environmental factors in water

从细菌与环境因子的关联性上看, 第一象限中, BuchneraLimbohabitansMalikiaPseudarcicella与NH4+-N、TN、TP和COD呈正相关关系, 与pH、DO呈负相关关系; 第二象限中, unidentified_α-变形菌属和unidentified_微酸菌属与pH、DO呈正相关关系, 与NH4+-N、TN、TP和COD呈负相关关系; 第四象限中, RhodoferaxPolynucleobacter与NH4+-N、TN、TP和COD呈正相关关系, 与pH和DO呈负相关关系.本研究发现影响细菌群落的最大影响因素为氮磷营养盐(TN、TP和NH4+-N), 这与其他学者对影响淡水细菌群落关键因子的研究结果一致[39].

2.5 Tax4Fun功能预测细菌群落代谢功能

Tax4Fun是用来分析预测细菌群落的代谢功能的R程序包, 根据数据库注释结果, 选取每个样本或分组在各注释层级上最大丰度排名前10的功能信息, 生成功能相对丰度柱状堆积图如图 7所示.西安市4个湖泊细菌群落代谢功能依次包括:代谢(metabolism)、遗传信息处理(genetic-information-processing)、环境信息处理(environmental-information-processing)、细胞进程(cellular-process)、人类疾病(human-diseases)和有机体系统(organismal-system).为更好地理解不同细菌群落的功能作用, 使用基于KEGG第2层级代谢途径数据分析库分析的Tax4Fun功能预测来获得细菌代谢途径.主要的代谢功能有:氨基酸代谢(amino acid metabolism)、碳水化合物代谢(carbohydrate metabolism)、能量代谢(energy metabolism)、核苷酸代谢(nucleotide metabolism)以及辅因子和维生素的代谢(metabolism of cofactors and vitamins)等.其中, 氨基酸代谢和碳水化合物代谢丰度最高, 营养来源是影响水环境细菌群落结构的潜在动力[40].氨基酸是有机氮的主要赋存形态, 其降解对细菌繁殖有促进作用, 代谢有益于细菌对氨基酸的利用等[41].氨基酸代谢相关细菌能促进氨基酸产生和腐殖质合成[42].动植物腐败物使湖泊中有机物积累, 异养型细菌利用碳水化合物代谢途径分解有机物[43], 湖泊中优势菌门(变形菌门和拟杆菌门)可以降解有机物, 进一步证明该区域有机物浓度较高.细菌基因功能一方面调节环境以维持自身生命活动, 另一方面细菌细胞间的群体感应抵抗了环境波动以维持细菌群落的相对稳定.由此可见, 细菌群落是水生生态系统中至关重要的一部分[44].

图 7 细菌群落代谢功能相对丰度 Fig. 7 Bacterial community functional relative abundance

景观湖泊分布于居民区, 西安市的城市景观湖泊大都具有游览观光和休闲娱乐, 甚至防洪保安和农业灌溉等功能, 对市民的生产生活及身体健康有重要影响[45].Sandifer等[46]研究表明环境中的细菌菌群与人体健康密切相关.西安市4个湖泊中存在的细菌种类较多, 其中不乏和人类疾病相关的细菌, 如图 8所示, 包括:耐药性(drug resistance)、传染性疾病(ifectious disease)、癌症(cancers)、内分泌及代谢性疾病(endocrine and metabolic)、神经退行性疾病(neurodegenerative disease)、心血管疾病(cardiorascular)、物质依赖(substance dependence)和免疫性疾病(immune disease). 其中以耐药性与传染性疾病相对丰度最高, 免疫性疾病最低, 这表明西安市4个湖泊存在潜在的耐药性细菌污染以及传染性疾病风险. 耐药细菌以及致病细菌可通过食物链向人类传播, 也会通过各种途径传播至周围环境[47]. 由于抗生素在水环境中的大量残留, 携带耐药性基因的细菌在环境中急剧增长, 因此, 对水环境潜在的耐药性细菌的监测和深入研究十分必要[48]; 水环境污染导致溶血性霍乱弧菌已成为我国细菌性食源性疾病的首要病因, 应加强传染性疾病的预防, 以确保水环境及人类身体健康[49].

图 8 细菌群落致病程度预测 Fig. 8 Bacterial community pathogenicity prediction

3 结论

(1) 从西安市4个湖泊水体细菌群落α多样性来看, 细菌多样性高低为:桃花潭>汉城湖>南湖>兴庆湖, 桃花潭细菌多样性最高, 兴庆湖最低.

(2) 变形菌门是4个湖泊水体中细菌群落占比最高的门类, 变形菌门在有机物分解和循环中有重要作用, 是地表水中常见的优势菌群.不同湖泊水体中的优势菌属及其数量存在一定差异, 但大都属于变形菌门、放线菌门和拟杆菌门.

(3) 西安市城市湖泊水体的理化性质与细菌群落结构存在显著相关性, 这些相关性也解释了湖泊水质分布规律, 氮磷营养盐为其细菌群落的关键环境因子(P < 0.05).

(4) Tax4Fun功能预测分析显示, 西安市4个湖泊中氨基酸代谢和碳水化合物代谢这2种代谢功能丰度最高, 细菌代谢较为活跃; 4个湖泊存在潜在的耐药性细菌污染以及传染性疾病风险较大, 因此应关注西安市景观湖泊水体细菌的潜在风险, 以保障水环境健康及可持续发展.

参考文献
[1] 程豹, 望雪, 徐雅倩, 等. 澜沧江流域浮游细菌群落结构特征及驱动因子分析[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.
[2] Zhu J Y, Zhang J X, Li Q, et al. Phylogenetic analysis of bacterial community composition in sediment contaminated with multiple heavy metals from the Xiangjiang River in China[J]. Marine Pollution Bulletin, 2013, 70(1-2): 134-139. DOI:10.1016/j.marpolbul.2013.02.023
[3] Wang L, Zhang J, Li H L, et al. Shift in the microbial community composition of surface water and sediment along an urban river[J]. Science of the Total Environment, 2018, 627: 600-612. DOI:10.1016/j.scitotenv.2018.01.203
[4] Pham T L, Utsumi M. An overview of the accumulation of microcystins in aquatic ecosystems[J]. Journal of Environmental Management, 2018, 213: 520-529.
[5] 张旭. 淡水水体底泥微生物群落特征及微生物指标评价研究[D]. 上海: 上海大学, 2016.
Zhang X. Microbial community characteristics and microbial indicators evaluation study in fresh water sediment[D]. Shanghai: Shanghai University, 2016.
[6] 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
[7] 黄晨. 西安市典型景观湖泊中抗生素分布特征及其生态风险评价[D]. 西安: 西安理工大学, 2021.
Huang C. Characteristics of antibiotic distribution and ecological risk assessment in typical landscape lakes in Xi'an[D]. Xi'an: Xi'an University of Technology, 2021.
[8] 张锋, 卢雪, 杨晓慧, 等. 西安市汉城湖水环境的调研报告[J]. 西安文理学院学报(自然科学版), 2017, 20(5): 109-112, 117.
Zhang F, Lu X, Yang X H, et al. Investigation report on water environment of Hancheng Lake in Xi'an City[J]. Journal of Xi'an University (Natural Science Edition), 2017, 20(5): 109-112, 117. DOI:10.3969/j.issn.1008-5564.2017.05.024
[9] 王敏, 黄晨, 隋蕾, 等. 兴庆湖中典型抗生素污染特征及其生态风险评价[J]. 西安理工大学学报, 2020, 36(1): 17-24.
Wang M, Huang C, Sui L, et al. Status of antibiotics contamination and ecological risks assessment on Xingqing Lake[J]. Journal of Xi'an University of Technology, 2020, 36(1): 17-24.
[10] 张彦波. 模糊综合评价法在地表水环境质量评价中的应用——以西安市曲江池水环境为例[J]. 西安石油大学学报(社会科学版), 2016, 25(4): 1-6.
Zhang Y B. Application of fuzzy comprehensive evaluation method to the assessment of surface water environment quality—with the example of surface water environment in Xi'an Qujiang Pool[J]. Journal of Xi'an Shiyou University (Social Science Edition), 2016, 25(4): 1-6. DOI:10.3969/j.issn.1008-5645.2016.04.001
[11] 刘清, 王松, 许江娟, 等. 西安市汉城湖浮游植物群落结构与水质状况分析[J]. 安全与环境工程, 2017, 24(3): 48-56.
Liu Q, Wang S, Xu J J, et al. Phytoplankton community structure and water quality condition of Hancheng Lake in Xi'an city[J]. Safety and Environmental Engineering, 2017, 24(3): 48-56.
[12] 李智, 王怡, 王文怀. 不同水源补给对景观水体水质及浮游动物的影响[J]. 中国给水排水, 2021, 37(5): 91-96, 105.
Li Z, Wang Y, Wang W H. Effects of different water sources supply on water quality and zooplankton in landscape waters[J]. China Water & Wastewater, 2021, 37(5): 91-96, 105.
[13] 崔芳. 利用水平潜流人工湿地净化城市湖泊污水——以西安市兴庆湖为例[J]. 湿地科学, 2015, 13(2): 207-210.
Cui F. Treatment of urban lake water by horizontal subsurface flow constructed wetlands: a case of Xingqinghu Lake in Xi'an city[J]. Wetland Science, 2015, 13(2): 207-210.
[14] Lutz K A, Wang W Q, Zdepski A, et al. Isolation and analysis of high quality nuclear DNA with reduced organellar DNA for plant genome sequencing and resequencing[J]. BMC Biotechnology, 2011, 11(1). DOI:10.1186/1472-6750-11-54
[15] Caporaso J G, Lauber C L, Walters W A, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample[J]. Proceedings of the National Academy of Sciences of the United States of America, 2011, 108(S1): 4516-4522.
[16] Mago Dč T, Salzberg S L. FLASH: fast length adjustment of short reads to improve genome assemblies[J]. Bioinformatics, 2011, 27(21): 2957-2963. DOI:10.1093/bioinformatics/btr507
[17] Caporaso J G, Kuczynski J, Stombaugh J, et al. QⅡME allows analysis of high-throughput community sequencing data[J]. Nature Methods, 2010, 7(5): 335-336. DOI:10.1038/nmeth.f.303
[18] Rognes T, Flouri T, Nichols B, et al. VSEARCH: a versatile open source tool for metagenomics[J]. PeerJ, 2016, 4. DOI:10.7717/peerj.2584
[19] Aßhauer K P, Wemheuer B, Daniel R, et al. Tax4Fun: predicting functional profiles from metagenomic 16S rRNA data[J]. Bioinformatics, 2015, 31(17): 2882-2884. DOI:10.1093/bioinformatics/btv287
[20] 马凤娟, 孙杰, 徐培智, 等. 生防菌XP1对香蕉枯萎病防效及土壤细菌群落多样性的影响[J]. 南方农业学报, 2019, 50(9): 1981-1989.
Ma F J, Sun J, Xu P Z, et al. Effects of biocontrol bacteria XP1 on banana fusarium wilt and diversity of bacterial community in soil[J]. Journal of Southern Agriculture, 2019, 50(9): 1981-1989.
[21] 蒋旭瑶, 吉喜燕, 黄德英, 等. 不同植物类型复合垂直流人工湿地根系微生物群落结构的研究[J]. 农业环境科学学报, 2019, 38(1): 176-183.
Jiang X Y, Ji X Y, Huang D Y, et al. Microbial community structure in the roots of three kinds of plants in integrated vertical flow constructed wetlands[J]. Journal of Agro-Environment Science, 2019, 38(1): 176-183.
[22] 王佳, 彭剑峰, 宋永会, 等. 浑河底泥微生物群落的季节性变化特征[J]. 环境科学研究, 2016, 29(2): 202-210.
Wang J, Peng J F, Song Y H, et al. Seasonal changes of microbial community distribution in sediments of Hun River[J]. Research of Environmental Sciences, 2016, 29(2): 202-210.
[23] 王晓丽, 其勒格尔. 黄河内蒙古段表层沉积物细菌多样性及群落结构类型[J]. 生态学报, 2020, 40(2): 578-589.
Wang X L, Qi Legeer. Bacterial diversity and community structure in surface sediments of Yellow River from inner Mongolia section[J]. Acta Ecologica Sinica, 2020, 40(2): 578-589.
[24] 薛银刚, 刘菲, 孙萌, 等. 太湖竺山湾春季浮游细菌群落结构及影响因素[J]. 环境科学, 2018, 39(3): 1151-1158.
Xue Y G, Liu F, Sun M, et al. Community structure and influencing factors of bacterioplankton in Spring in Zhushan Bay, Lake Taihu[J]. Environmental Science, 2018, 39(3): 1151-1158.
[25] 胡愈炘, 张静, 黄杰, 等. 长江流域河流和湖库的浮游细菌群落差异[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.
[26] 赵志瑞, 吴会清, 毕玉方, 等. 补给水质与社会活动对白洋淀湿地微生物的影响特征[J]. 环境科学, 2021, 42(7): 3272-3280.
Zhao Z R, Wu H Q, Bi Y F, et al. Effects of water supply quality and social activity on the microbial community in Baiyangdian wetland[J]. Environmental Science, 2021, 42(7): 3272-3280.
[27] Chakravarthy S S, Pande S, Kapoor A, et al. Comparison of denitrification between Paracoccus sp. and Diaphorobacter sp.[J]. Applied Biochemistry and Biotechnology, 2011, 165(1): 260-269.
[28] Koskinen K, Hultman J, Paulin L, et al. Spatially differing bacterial communities in water columns of the northern Baltic Sea[J]. FEMS Microbiology Ecology, 2011, 75(1): 99-110.
[29] 徐超, 张军毅, 朱冰川, 等. 夏季太湖梅梁湾水体中细菌的群落结构[J]. 环境监控与预警, 2015, 7(1): 37-40.
Xu C, Zhang J Y, Zhu B C, et al. Extensive profiling of a bacterial community during summer in Meiliang Bay, Lake Taihu[J]. Environmental Monitoring and Forewarning, 2015, 7(1): 37-40.
[30] 王泽铭, 李传虹, 马巧丽, 等. 湿度盐度pH协同驱动锡林河景观疣微菌群空间异质性[J]. 微生物学报, 2021, 61(6): 1728-1742.
Wang Z M, Li C H, Ma Q L, et al. Moisture, salinity and pH co-driving spatial heterogeneity of Verrucomicrobial populations in Xilin River landscape[J]. Acta Microbiologica Sinica, 2021, 61(6): 1728-1742.
[31] 韩雪梅, 龚子乐, 杨晓明, 等. 汛期前后老鹳河干流人类干扰下浮游细菌多样性及功能预测[J]. 环境科学, 2021, 42(2): 831-841.
Han X M, Gong Z L, Yang X M, et al. Diversity and function prediction of bacterioplankton under human disturbance in the main stream of the Laoguan River before and after the flood season[J]. Environmental Science, 2021, 42(2): 831-841.
[32] 朱金山, 秦海兰, 孙启耀, 等. 冬季小流域水体微生物多样性及影响因素[J]. 环境科学, 2020, 41(11): 5016-5026.
Zhu J S, Qin H L, Sun Q Y, et al. Microbial diversity and influencing factors in a small watershed in winter[J]. Environmental Science, 2020, 41(11): 5016-5026.
[33] 周石磊, 张艺冉, 孙悦, 等. 异养硝化-好氧反硝化菌富集驯化过程中微生物种群演变特征——典型城市景观水系[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.
[34] 张菲, 田伟, 孙峰, 等. 丹江口库区表层浮游细菌群落组成与PICRUSt功能预测分析[J]. 环境科学, 2019, 40(3): 1252-1260.
Zhang F, Tian W, Sun F, et al. Community structure and predictive functional analysis of surface water bacterioplankton in the Danjiangkou reservoir[J]. Environmental Science, 2019, 40(3): 1252-1260.
[35] 梁沪莲, 郭小雅, 刘洋, 等. 基于高通量测序的4种硝化细菌富集培养物微生物群落结构分析[J]. 微生物学通报, 2017, 44(9): 2112-2119.
Liang H L, Guo X Y, Liu Y, et al. Microbial community of four nitrobacteria enrichments by high-throughput sequencing[J]. Microbiology China, 2017, 44(9): 2112-2119.
[36] Zhang L, Zhong M M, Li X C, et al. River bacterial community structure and co-occurrence patterns under the influence of different domestic sewage types[J]. Journal of Environmental Management, 2020, 266. DOI:10.1016/j.jenvman.2020.110590
[37] Zhang J X, Yang Y Y, Zhao L, et al. Distribution of sediment bacterial and archaeal communities in plateau freshwater lakes[J]. Applied Microbiology and Biotechnology, 2015, 99(7): 3291-3302.
[38] 姜磊. 会仙湿地生态恢复中的细菌群落结构及其与环境因子间的关联研究[D]. 桂林: 桂林理工大学, 2019.
Jiang L. Bacterial community structure and its relationship with environmental factors in ecological restoration Huixian wetland[D]. Guilin: Guilin University of Technology, 2019.
[39] Wang P, Chen B, Yuan R Q, et al. Characteristics of aquatic bacterial community and the influencing factors in an urban river[J]. Science of the Total Environment, 2016, 569-570: 382-389.
[40] Wilhelm L, Besemer K, Fragner L, et al. Altitudinal patterns of diversity and functional traits of metabolically active microorganisms in stream biofilms[J]. The ISME Journal, 2015, 9(11): 2454-2464.
[41] 万甜, 何梦夏, 任杰辉, 等. 渭河流域水体细菌群落的环境响应及生态功能预测[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.
[42] 丁钰珮, 杜宇佳, 高广磊, 等. 呼伦贝尔沙地樟子松人工林土壤细菌群落结构与功能预测[J]. 生态学报, 2021, 41(10): 4131-4139.
Ding Y P, Du Y J, Gao G L, et al. Soil bacterial community structure and functional prediction of Pinus sylvestris var. mongolica plantations in the Hulun Buir sandy land[J]. Acta Ecologica Sinica, 2021, 41(10): 4131-4139.
[43] Liu L L, Sun F F, Zhao H B, et al. Compositional changes of sedimentary microbes in the Yangtze River Estuary and their roles in the biochemical cycle[J]. Science of the Total Environment, 2021, 760. DOI:10.1016/j.scitotenv.2020.143383
[44] 许郑超, 邓燕飞, 周军, 等. 中华绒螯蟹养殖池塘菌群结构和功能特征[J]. 生态学杂志, 2021, 40(7): 2223-2233.
Xu Z C, Deng Y F, Zhou J, et al. Structural and functional characteristics of microbial communities in aquaculture ponds of Chinese mitten crab Eriocheir sisensis L.[J]. Chinese Journal of Ecology, 2021, 40(7): 2223-2233.
[45] 蒋子银. 城市公园生态景观美景度评价研究——以西安市为例[D]. 西安: 西安建筑科技大学, 2017.
Jiang Z Y. Beauty evaluation and research of ecological landscape in Xi'an city[D]. Xi'an: Xi'an University of Architecture and Technology, 2017.
[46] Sandifer P A, Sutton-Grier A E, Ward B P. Exploring connections among nature, biodiversity, ecosystem services, and human health and well-being: opportunities to enhance health and biodiversity conservation[J]. Ecosystem Services, 2015, 12: 1-15.
[47] Zhou H, Zhang D G, Jiang Z H, et al. Changes in the soil microbial communities of alpine steppe at Qinghai-Tibetan Plateau under different degradation levels[J]. Science of the Total Environment, 2019, 651: 2281-2291.
[48] 黎丽雯, 邹立扣, 周杨, 等. 成都府河水体细菌动态变化及耐药性[J]. 湖泊科学, 2012, 24(1): 96-103.
Li L W, Zou L K, Zhou Y, et al. Dynamics of bacterial diversity and antibiotic resistance of Escherichia coli in Fu River, Chengdu[J]. Journal of Lake Sciences, 2012, 24(1): 96-103.
[49] 王江权, 康敉, 郑祥, 等. 海河流域典型河流粪源性指示微生物的污染特征及其时空分布[J]. 环境科学学报, 2017, 37(1): 138-145.
Wang J Q, Kang M, Zheng X, et al. Occurrence and temporal-spatial distribution of fecal indicator microorganisms in three rivers of the Haihe River basin[J]. Acta Scientiae Circumstantiae, 2017, 37(1): 138-145.