ACTA Scientiarum Naturalium Universitatis Pekinensis
水库底泥中微生物多样性及其与环境因子相关性分析
1. 国家环境保护饮用水水源地管理技术重点实验室, 深圳市饮用水水源地安全保障重点实验室, 深圳市水环境中新型污染物检测与控制重点实验室, 深圳市环境科学研究院, 深圳 518001; 2. 生物资源保护与利用湖北省重点实验室, 湖北民族学院,恩施 445000; † 通信作者, E-mail: pkutianck@163.com
摘要 为评价广东某水库底泥中微生物多样性及其与环境因子的相关性, 为该区域微生物的生态功能研究提供理论基础, 选择水库库中、外部引水入口及主要支流入口共 6 个采样点, 进行不同季节的水质特征分析,结合 PCR 扩增技术和 T-RFLP 技术, 对不同时期各点水体中微生物进行多样性分析, 并分析微生物群落结构多样性与水库环境条件的关系。从总体上看, 水库底泥中各采样点间隙水的理化性质具有较高的时空差异性, 库中不同形式氮的浓度和有机碳浓度低于外部供水入口及支流, 各指标浓度大体上呈现丰水期>枯水期>平水期的现象。枯水期水库底泥的微生物比平水期和丰水期丰富, 不同时期不同采样点的细菌多样性分布也呈现不同规律。平水期的优势菌属主要有Family1_uncultured, Bacteroidetes_uncultured, Sphingobacteriales_ uncultured和anaerolineaceae, 丰水期的优势菌属主要有Peptostreptococcaceae_incertae_sedis, Anaerolineaceae, SCI-84, Xanthomonadales和 Clostridium, 枯水期的优势菌属主要有Rhodocyclaceae, Fusobacteriales, BSV26_ norank, Comamonadaceae和anaerolineaceae。水库底泥中细菌具有丰富的多样性, 水体微生物群落结构在不同点位具有一定的差异性, 受环境因子的综合影响, 不同形式氮浓度的影响尤为明显。关键词 水库; 底泥; 微生物多样性; 环境因子; 相关性分析中图分类号 X172
behaved the higher level of bacterial community diversity and richer degree in dry season than that in normal season and wet season. The dominant bacteria were Family1_uncultured, Bacteroidetes_uncultured, Sphingobacteriales_uncultured, Anaerolineaceae in normal season, Peptostreptococcaceae_incertae_sedis, Anaerolineaceae, SC-I-84, Xanthomonadales, and Clostridium in wet season, Rhodocyclaceae, Fusobacteriales, Bsv26_norank, Comamonadaceae, Anaerolineaceae in dry season. Bacterial diversity was rich in the sediments in the reservoir. Correlation analysis showed that the microbial community structure was different at different sites, effected by environmental factors, among which the effects of different forms of nitrogen concentration were most obvious. Key words reservoir; sediments; microbial diversity; environmental factors; correlation analysis
随着工业化和城市化进程的加快, 我国很多自然水体受到严重污染, 大量污染物沉积到底泥中,使得底泥成为微生物种类丰富、物质交换频繁的载体。水体和底泥中微生物是水域生态系统的重要组成部分, 微生物群落状况与水质状况有密切的联系,其结构的变化对水质污染负荷积累有很好的响应。微生物在生源要素循环、污染物降解和净化等方面
[1–4]起着非常重要的作用 。因此, 沉积物中微生物的群落结构是衡量水体生态系统稳定性的重要指标之一。认识底泥中微生物群落分布特征, 对于阐明微生物群落与其生境的关系, 有针对性地开展污染系统的生物修复具有重要意义。环境条件的改变会引起微生物群落结构及其组成的变化[5–7], 不同季节水库底泥中的微生物群落结构亦存在季节性差异。目前, 关于河流和湖泊沉积物中微生物多样性的研究已有一些报道, 但多针对单一季节, 关于水库底泥微生物群落结构季节性变化特征的研究比较少见[8]。对不同环境条件下微生物群落结构的多样性进行研究, 对于揭示微生物与环境之间的关系以及环境综合评价分析等方面都具有重要的作用[9–10]。本文分别考察平水期、枯水期和丰水期广东某水库的水质理化指标、底泥的微生物群落结构情况以及响应关系, 旨在更深入地掌握水库底泥微生物的时空分布及群落结构的季节性变化特征, 并通过测序, 了解不同时期底泥中的主要微生物, 以期为水库生态系统的保护和管理提供参考。
1材料和方法1.1样品的采集
分别于 2014 年 4 月(平水期)、2014 年 8 月(丰水期)和2015 年 1 月(枯水期), 在广东省某水库采集底泥样品。根据水库及其主要入库支流的环境与基本特征, 本研究布设 6 个采样点, 包括外部引水入口 S7、水库库中 3 个点(S8, S9, S10)、外部引水入口 S11 和主要支流入库口 S12。用抓斗式采样器
采集底泥样品, 样品采集、保存与运输符合标准HJ/T 166—2004。
1.2 理化指标的测定
测定泥样间隙水中TOC, NH4+-N, NO2−-N, NO3−N, TN, ph 和 ORP。NH4+-N 测定采用钠氏试剂分光光度法(GB 7479—87), NO2−-N 测定采用 N-(1-萘基)乙二胺分光光度法(GB 13580.7—92), NO3−-N 测定采用酚二磺酸分光光度法(GB/T 7480—1987), TN 测定采用过硫酸钾氧化–紫外分光光度法(GB/T 11894— 1989), TOC 测定采用燃烧氧化–非分散红外吸收法(GB 13193— 91), ph 和 ORP 测定采用瑞士 Mettler公司 S210 Seven Compact ph 计。
1.3 样品基因组总DNA 提取和 PCR 扩增
滤纸上的微生物利用液氮冷却, 研磨成粉后用Mpbio 公司的 FASTDNA® Spin Kit for Soil 提取基因组 DNA。底泥中微生物基因组 DNA 直接用 Mpbio公司的 FASTDNA® Spin Kit for Soil 提取。
®23 μl PCR 反应体系,采用 使用 PCR 仪(ABI, Geneamp 9700 型), 对细菌16S V4-V5区间基因进行扩增, 细菌16S V4-V5区间基因 PCR 扩增的反应体系和 PCR 程序如下所示。
扩增引物: 515F: GTGCCAGCMGCCGCGG, 907R: CCGTCAATTCMTTTRAGTTT。
反应体系: 5×Fastpfu Buffer 4 μl, 2.5 mm DNTPS 2 μl, Forward Primer (5 μm) 0.6 μl, Reverse Primer (5 μm) 0.6 μl, Fastpfu Polymerase 0.4 μl, Template DNA 10 ng, 补 ddh2o 至 20 μl。
PCR 程序: 1) 1 cycle × (95°C 3分钟); 2) 25 cycles × (95oc 30秒, 55oc 30秒, 72oc 45秒); 3) 1 cycle × (72oc 10分钟); 4) 10oc 直至结束。
1.4 T-RFLP 分析
利用正向引物27F (序列: 5’-AGAGTTTGATCC TGGCTCAG-3’, 5’-AFM 标记)与反向引物1492R (序列: 5’-GGTTACCTTGTTACGACTT-3’)进行16S
RDNA 扩增。扩增反应采用 25 μl 体系, 模板量为2 μl, 阴性对照用 ddh2o 代替模板 DNA。PCR 扩增条件为: 94oc 5 分钟, 1 cycle; 94oc 30 秒, 54oc 30秒, 72oc 1 分钟 30 秒, 35 cycles; 延伸 72oc 10分钟, 1 cycle。PCR 产物纯化后, 取适量用 Alu I 进行酶切, 反应采用 10 μl 体系, 37oc 3 小时酶切, 完成后65oc 20分钟, 终止反应。取1 μl 酶切产物+12 μl甲酰胺+0.2 μl 内标混匀后, 在 ABI310 遗传分析仪上进行电泳, 电泳时间为25分钟。酶切产物通过毛细管电泳 ELITE300进行荧光扫描。结果分析: 每一个荧光峰至少代表一个细菌或几个近缘的细菌,可粗略地认为峰面积与细菌在某一部位的微生态群落中的数量有关。
1.5 微生物群落结构以及多样性分析
将滤纸上的微生物利用液氮冷却, 研磨成粉后用 Mpbio 公司的 FASTDNA® Spin Kit for Soil 提取基因组DNA。底泥中微生物基因组DNA直接用Mpbio公司的 FASTDNA® Spin Kit for Soil 提取。聚合酶链式扩增引物采用细菌通用引物515F: GTGCCAGC MGCCGCGG, 907R: CCGTCAATTCMTTTRAGTTT。每个样品的扩增均做 3 次重复, PCR反应体系为20 μl, 其中包括4 μl的5 × Fastpfu Buffer, 2 μl 的2.5 mm DNTPS, 引物 (5 μm) 各 0.8 μl, 0.4 μl 的Fastpfu Polymerase 以及 10 ng 模板 Dna。illumina Miseq 测序从2%琼脂糖凝胶中回收扩增子, 依照Axyprep DNA Gel Extraction Kit 的说明书进行纯化(Axygen Biosciences, Union City, CA, US), 再用Quantifluor™-st (Promega, US)进行定量分析。将纯化后的扩增子等量混合, 然后根据 Illumina Miseq 测序平台的标准流程进行双端(2×250 bp 或2×300 bp)测序。原始数据提交到 NCBI Sequence Read Archive (SRA)数据库。
1.6 冗余分析
结合高通量测序所得微生物群落结构组成结果与所测定的环境因子数据, 应用 Canoco 软件完成冗余分析(redundancy analysis, RDA)。根据去趋势波动分析(detrended correspondence analysis, DCA)结果, 选用 RDA 进行群落结构与环境因子的影响分析。在 RDA 二维图中, 用带箭头的线段表示环境因子, 从样方点到数量型环境因子的箭头投影长度可表示样方内环境因子的影响值, 进而反映该环境因子对样方点群落组成的影响。
2 结果与讨论2.1 底泥样品的理化性质分析
底泥样品间隙水的理化性质分析结果如表 1 所示。库中不同点位的底泥(S8, S9, S10)中, 各环境因子的检测结果在平水期、丰水期和枯水期具有一定的差异性, 库内底泥中不同形式氮的浓度和有机碳浓度均低于外部引水入口和入库支流(S7, S11, S12), 且外部引水入口和支流入口处的底泥环境同样存在一定的时空差异性。各指标浓度大体上呈现丰水期>枯水期>平水期的现象。在外部引水入口S11的底泥中, 不同形式氮的浓度和有机碳浓度明显高于库区以及外部引水入口S7 和支流入口 S12。
2.2 水库沉积物中微生物多样性分析
底泥微生物多样性的测序分析结果如表 2 所示。除平水期、丰水期和枯水期检测的序列条数在25000~30000之间, 覆盖率达97%以上, 按97%的相似度进一步归类后, 得到一定数量的 Operational Ttaxonomic Unit (OTU)。不同季节沉积物中微生物多样性的 Shannon 指数差异如图 1 所示。底泥样品的 Shannon 指数大部分处于6~7之间, 枯水期多样性同样高于平水期和丰水期。这主要由于枯水期水流量小, 流速较慢, 底泥环境比较稳定, 利于其中的微生物生存繁殖; 在丰水期和平水期, 底泥受到水流的冲刷, 环境存在波动, 不利于其中的微生物生存繁殖。
图2~4显示水库沉积物菌群多样性的 T-RFLP分析结果, 可以看出, 底泥样品中平水期的细菌多样性低于丰水期和枯水期, 不同时期细菌多样性的差异性可能由于丰水期和枯水期降水或者地下水的补给造成底泥环境的波动, 使得底泥中细菌群落结构发生变化。
2.3 沉积物微生物群落结构组成
沉积物微生物群落结构组成如图5 所示。平水期的优势菌属主要有 Family1_uncultured, Bacteroidetes, Sphingobacteriales_uncultured 和 Anaerolineaceae,丰水期的优势菌属主要有 Peptostreptococcaceae_ incertae_sedis, Anaerolineaceae, SC-I-84, Xanthomonadales 和 Clostridium, 枯水期的优势菌属主要有Rhodocyclaceae, Fusobacteriales, Bsv26_norank, Comamonadaceae 和 Anaerolineaceae。丰水期, 库中采样点 S10的主要菌属与 S8, S9不同, 主要是Vadinha17_norank 和 Anaerolineaceae, 而采样点 S8, S9的主要菌属是 Peptostreptococcaceae_incertae_
sedis。说明即便在库内不同采样点, 底泥的菌属组成也存在差异。外部引水和支流入口中, 采样点S11的主要菌属是 Peptostreptococcaceae_incertae_ sedis, 采样点 S7的主要菌属是 Sc-i-84_norank 和Anaerolineaceae, 采样点 S12的主要菌属是 Xanthomonadales和 Caldidate, 底泥的微生物组成与丰水期的水流补给结果并不一致。枯水期, 库中采样点S9的主要菌属与 S8, S10不同, 主要是 BSV26_ norank, 比例达到10%, 而采样点 S8, S10的主要菌属除 Bsv26_norank 外, 还包括 Rhodocyclaceae 和Fusobacte-riales。外部引水入口 S11的主要菌属与
外部引水入口 S7和支流入口 S12不同, 主要是 Rhodocyclaceae,外部引水入口 S7的主要菌属是 Rhodocyclaceae 和 Bsv26_norank, 支流入口 S12的主要菌属是 Comamonadaceae。在不同时期底泥中微生物的菌属组成存在较大差异, 主要原因是底泥的环境随着丰水期水量的增大波动较剧烈, 对微生物组成影响较大。
2.4 沉积物的微生物群落结构基本特征及其对环境因子响应分析
沉积物微生物群落结构与水质因子响应 RDA分析结果如图 6 所示。平水期, 外部引水入口 S7 受亚硝氮影响, 外部引水入口 S11 的微生物群落结构主要受亚硝氮和总氮的影响, 氨氮、硝氮和总有机碳浓度的影响也很大, 支流入口 S12 主要受有机碳和硝氮浓度影响, 氨氮、总氮和亚硝氮浓度的影响也比较明显, 外部引水入口 S7 主要受亚硝氮浓度的影响, 库中 3个采样点的微生物群落结构与各环境因子没有呈现正相关的关系。丰水期, 外部引水入口 S7和支流入口 S12的微生物群落结构主要受亚硝氮和有机碳浓度的影响, 库中 S8的微生物群落结构除了受亚硝氮、有机碳浓度的影响, 还与总氮、氨氮和硝氮都呈正相关的关系, 库中 S9 和外部引水入口 S11 的微生物群落结构主要受总氮、氨氮和硝氮的影响。枯水期, 外部引水入口 S7 的微生物群落结构主要受各种形式氮浓度的影响, 其中影响最大的是亚硝氮, 库中 S8, S9 和 S10 这 3个点的微生物群落结构与各环境因子之间没有呈现正相 关的关系, 外部引水入口 S11的微生物群落结构与各种环境因子之间都有正相关性, 其中影响最大的是氨氮、总氮和硝氮。以上结果表明在水库的底泥中, 不同时期、不同采样地点, 环境因子对微生物群落组成的影响不同, 环境因子对微生物群落组成存在综合影响, 其中不同形式氮的浓度是主要影响因素。
3 讨论
各指标浓度大体上呈现丰水期>枯水期>平水期的现象。库中底泥环境具有较高的时空差异性,且底泥中的有机物浓度普遍偏高, 尤其是在外部引水入口处, 需要对污染源做进一步的调查和防控。底泥中的微生物多样性十分丰富, 种群结构复杂,且底泥在枯水期的细菌丰度和多样性高于平水期和丰水期, 枯水期水库的水文水利条件可能更有利于微生物的繁殖。大量研究表明, 生境中环境因子以及其他营养元素等显著影响微生物的种类组成, 某种营养元素含量增加, 相对应的功能性微生物数量也会随之改变[11]。我们也得到了类似的研究结果。相关性分析结果表明, 同一时期不同采样点的沉积物中微生物组成存在一定的差异性, 在不同时期的微生物组成也不同, 但同一时期不同点位的微生物组成相差不大。平水期的优势菌属主要有 Family1_ uncultured, Bacteroidetes_uncultured, Sphingobacteriales_uncultured和 Anaerolineaceae, 丰水期的优势菌属主要有 Peptostreptococcaceae_incertae_sedis, Ana-
erolineaceae, SC-I-84, Xanthomonadales 和 Clostridium,枯水期的优势菌属主要有 Rhodocyclaceae, Fusobacteriales, Bsv26_norank, Comamonadaceae 和Anaerolineaceae。不同时期的优势菌属呈现一定的相似性, 同时存在一定的差异性。环境因子响应分析结果表明, 在不同时期、不同采样点, 影响微生 物群落结构组成的主要环境因子不同, 但是都受到氮浓度的影响, 这是由于自然界中的固氮作用、氨化作用、硝化作用和反硝化作用都离不开微生物的参与。
参考文献
[1] 柴晓娟, 骆大伟, 吴春笃, 等. 水体中微生物分布及与环境要素的相关性研究. 人民长江, 2008, 39 (3): 45–48 [2] 任丽娟, 何聃, 邢鹏, 等. 湖泊水体细菌多样性及其生态功能研究进展. 生物多样性, 2013, 21(4): 421– 432 [3] Newton R J, Jones S E, Eiler A, et al. A guide to the natural history of freshwater lake bacteria. Microbiology and Molecular Biology Reviews, 2011, 75(1): 14–49 [4] Zwart G, Crump B C, Agteveld M P K V, et al. Typical freshwater bacteria: an analysis of available 16S RRNA gene sequences from plankton of lakes and rivers. Aquatic Microbical Ecology, 2002, 28(2): 141– 155 [5] 全向春, 王育来, 何孟常, 等. FISH 方法解析大辽河入海口沉积物中微生物的垂直分布特征. 环境科学学报, 2009, 29(7): 1502–1509 [6] Chen J, Xie P. Tissue distributions and seasonal dynamics of the hepatotoxic microcystins-lr and -RR in two freshwater shrimps, Palaemon modestus and Macrobrachium nipponensis, from a large shallow, eutrophic lake of the subtropical China. Toxicon, 2005, 45(5): 615–625 [7] 周霜艳, 陈小娜, 崔磊, 等. 秦皇岛新开河河口邻近海域细菌多样性及其影响因子. 微生物学通报, 2016, 43(12): 2578–2593 [8] 宋巍, 安德荣, 刘雪, 等. DGGE 分析东江流域农村饮用水源中微生物多样性及其与环境因子相关性.微生物学通报, 2009, 32(9): 1311–1317 [9] Wang Z H, Yang J Q, Zhou J, et al. Composition and structure of bacterial communities in waste water of aquatic products processing factories. Research Journal of Biotechnology, 2014, 9(2): 65–70 [10] 王佳, 彭剑峰, 宋永会, 等. 浑河底泥微生物群落的季节性变化特征. 环境科学研究, 2016, 29(2): 202–210 [11] 圣平, 于一尊, 王东升, 等. 通过 PCR-DGGE 技术分析鄱阳湖不同河口底泥中微生物多样性. 江西科学, 2015, 33(3): 312–317