文章信息
- 王丽娜, 陈辉煌, 刘乐冕, 余正, 杨军
- WANG Lina, CHEN Huihuang, LIU Lemian, YU Zheng, YANG Jun.
- 亚热带分层水库固氮微生物时空分布格局
- Spatiotemporal patterns of diazotrophic microorganisms in a subtropical stratified reservoir
- 生态学报[J]. 2016, 36(18): 5827-5837
- Acta Ecologica Sinica[J]. 2016, 36(18): 5827-5837
- http://dx.doi.org/10.5846/stxb201504250854
-
文章历史
- 收稿日期: 2015-04-25
- 网络出版日期: 2016-01-05
2. 中国科学院大学, 北京 100049
2. University of Chinese Academy of Sciences, Beijing 100049, China
生物固氮是指通过生物作用将大气中的氮气还原成氨的生态过程, 该过程通常由特殊且多样的固氮微生物来实现[1-2]。随着分子技术与基因组分析方法被引入到固氮微生物自然群落研究中, 人们对水环境固氮微生物群落有了一定的认识[3-5]。目前, 利用固氮基因研究固氮微生物的分布及多样性逐渐得到国际学者的认同并被广泛应用[1-3, 6]。例如, 利用nifH基因检测到的固氮微生物包括:蓝藻(蓝细菌)、变形菌、厚壁菌及放线菌等[3, 7-12]。固氮作用的活性和调节机制在不同的微生物类群中存在很大差异, 其生态学或生物学过程和机制并没有完全了解清楚[6, 13]。通常影响微生物固氮作用的因素不是一个因子能够解释的, 而是由一系列的生态因子和固氮微生物群落组成共同决定[14]。Dang[15]等人研究中国南海12个站位的固氮微生物多样性、丰度、生物地理分布及群落结构, 结果显示温度是影响沉积物固氮微生物分布的最显著因子, 暗示在沉积物中固氮作用是很活跃的, 其他一些因子也会对固氮微生物群落组成造成一定的影响, 如PO43--P等。Farnelid[16]等人于2007、2008和2011年在波罗的海从表层到底层取样, 基于454高通量测序发现固氮微生物多样性很高, 并且具有明显的垂直分层规律, 认为氮固定过程不仅仅限于氮素缺乏的表层水体, 而在底层缺氧环境同样很重要。陈因[17]等人以蓝藻Anabaena 7120为例研究光照对固氮活性的影响, 发现一系列光合抑制剂对固氮蓝藻活性有抑制作用, 而且弱光下此种抑制会加强。Howarth[18]等人总结了固氮作用在淡水、河口及海洋系统中生物地球化学控制机制, 认为氧气和溶解性无机氮对固氮作用有抑制作用, 低水平的氮磷比、以及微量元素铁和钼对固氮作用有促进作用。总的来看, 已有的水生固氮微生物研究多限于海洋环境, 固氮作用与水质, 特别是氮的含量密切相关[6, 9-12]。然而, 人们对淡水固氮微生物与环境因子之间的联系及其时空动态变化的研究较少[8]。
在淡水生态系统中, 生物固氮受到氮磷浓度与光照强度的影响[19], 然而不同固氮微生物多样性和活性的时空分布规律还不是很明确[13]。Short[20]等人在切萨皮克湾生态系统中发现两种nifH基因型有着不同的季节和空间分布格局, 暗示这些生物群落受到不同环境过程的影响。Moisander[9]等人研究切萨皮克湾水柱和沉积物固氮微生物组成的时空变化, 发现多样性最高的是nifH Cluster 3, 但是Cluster 1(α-, β-, γ-变形菌, 蓝藻)具有较高的丰度;底层水体固氮微生物群落具有较高的相似性, 而表层水体群落变异较大;盐度、溶解性无机氮、溶解性有机碳和溶解性有机磷等环境因子与固氮微生物群落组成有显著的关系。Yu[12]等人研究东圳水库中固氮微生物秋季垂直分布特征, 发现温跃层的出现使水温与溶解氧发生明显垂直变化, 从而对固氮微生物群落组成造成了显著影响。以上研究报道显示, 水生固氮微生物的垂直分布和时空动态过程的相关研究是最新发展趋势和研究热点。
我国亚热带深水水库多属于典型的暖单次混合型水库, 除冬季有短暂的混合期外, 在一年当中会持续出现热分层现象[21]。在这些水库中, 明显的热分层是春季、夏季与秋季的典型特征, 在温跃层上下水体中明显差异的环境对于微生物群落的组成、多样性及丰度可能有着显著的影响[22]。显然, 亚热带分层水库同样非常适合研究固氮微生物垂直分布和季节动态变化规律。本研究以厦门汀溪水库为例, 综合应用克隆文库、定量PCR、定量RT-PCR技术分析固氮微生物在4个季节不同水层的组成、多样性、活性及丰度。本研究的目的是:1)揭示分层水库固氮微生物群落的季节和空间分布格局;2)探讨固氮微生物优势类群的特征与成因。
1 材料与方法 1.1 样品采集及理化分析汀溪水库(24°48′N, 118°08′E)位于福建省厦门市, 是一座典型的亚热带深水水库。本研究选取汀溪水库湖泊区的深水区设置站位进行采样, 此站位水深可达25 m。水体样品采集周期为1a:2012年的7月(夏季), 10月(秋季), 以及2013年的1月(冬季), 4月(春季)。每次采样均采集5个水层的样品, 分别是:0、6、12、18、24 m(夏), 0、10、15、20、25 m(秋), 0、5、10、15、21 m(冬), 0、5、10、15、20 m(春)。每份样品均采集10 L水样, 水样采集后立即送往实验室进行过滤, 采用孔径为0.22 μm的滤膜(直径47 mm, Millipore, 美国)过滤1 L水样, 并保留4份重复样品。将载有微生物的滤膜储存在-80℃的冰箱, 用于后续DNA和RNA的提取。
现场利用美国哈希HACH多参数水质分析仪(Hydrolab DS5)监测水温、溶解氧、pH、叶绿素a等参数, 每隔1 m测定1次。水体营养盐分析在实验室进行, 总氮、氨氮、硝态氮、总磷、正磷酸盐均按照国家标准进行测定。
1.2 DNA和RNA提取根据说明书, 利用试剂盒Fast DNA spin kit(MP Biomedicals, 美国)和E.Z.N.A. total RNA kit(Omega Bio-Tek, 美国)分别提取固氮微生物DNA和RNA。提取的RNA立即反转录成cDNA, 使用的是试剂盒OneStep RT-PCR kit Version 2.0(Takara, 日本)。反转录程序是:37 ℃, 15 min;95 ℃, 5 s。提纯后的DNA和cDNA均置于冰箱-40 ℃保存备用。
1.3 PCR扩增及克隆文库分析本研究扩增nifH基因选择的引物是nifH-F (5′-CGTAGGTTGCGACCCTAAGGCTGA-3′)和nifH-R (5′-GCATACATCGCCATCATTTCACC-3′)[23]。PCR的扩增体系是50 μL, 包含:5 μL 10 × PCR缓冲液(10 nmol/L), 4 μL dNTPs(10 mmol/L), 0.5 μL Taq聚合酶(5 U/μL), 前后引物各1 μL(10 μmol/L), 100 ng DNA模板, 然后使用ddH2O补充到50 μL。PCR反应体系如下:94 ℃预变性5 min, 然后40个循环(94 ℃变性30 s, 58 ℃退火30 s, 72 ℃延伸40 s), 最终72 ℃延伸10 min。PCR产物利用1%的琼脂糖凝胶进行电泳检测。利用胶回收试剂盒提纯PCR产物, 然后根据操作指南将其连接到pMD18-T载体(Takara, 日本)。将连接后的载体转移到大肠杆菌DH5α(Takara, 日本)中, 置于液体培养基中37 ℃培养过夜。利用添加氨苄的固体培养基进行蓝白斑筛选, 挑取白色克隆进行测序分析。本文选用四个季节的表层和底层样品, 共构建8个克隆文库, 每个克隆文库挑选70个克隆子进行测序。
1.4 定量PCR定量PCR与克隆文库所用的引物是相同的。利用Applied Biosystems 7500 Real-Time PCR System(ABI, 美国)进行定量PCR, 每个样品设置3个重复[12]。本研究采用的是20 μL体系:0.5 μL前后引物(10 μmol/L), 10 μL SYBR Premix Ex TaqTM, 2 μL DNA或cDNA模板和7 μL RNase-free water。本研究的标准曲线的R2值为0.994(DNA)和0.991(RNA), 扩增效率分别是97.45%(DNA)和98.19%(RNA)。RNA/DNA比值可以作为反映固氮微生物基因表达活性的参数, 因此也进行了分析。
1.5 数据分析克隆文库的操作分类单元(Operational taxonomic unit, OTU)划分采用97%的相似度, 利用软件MOTHUR v.1.33.0分析8个文库的稀释曲线, 然后计算OTU个数、覆盖率、Shannon-Wiener指数、Simpson倒数指数、Pielou均匀度指数、以及反映OTU丰富度的Chao1指数和ACE指数[24]。将测序获得的序列在GenBank数据库中进行比对, 查找最相似序列并确定其分类归属, 统计分析主要类群的OTUs比例和序列比例。利用PRIMER 6.0计算样品间的Bray-Curtis相似系数, 进而对样品进行聚类分析[25]。利用SPSS 19.0进行单因素方差方法, 确定混合期及分层水体中不同水层之间水体理化及生物指标的差异是否显著(P < 0.05), 应用Spearman相关分析确定DNA拷贝数、RNA拷贝数及RNA/DNA比值与环境因子之间的相关性是否显著(P < 0.05)。
2 结果 2.1 水库环境因子的时空变化汀溪水库是一座典型的亚热带分层水库, 在冬季呈现短时期的上下混合, 在春季、夏季和秋季出现稳定的热分层(图 1)。本研究中, 春季4月份、夏季7月份、秋季10月份水库的温跃层分别出现在3—14 m、18—22 m、21—24 m。冬季1月份处于上下混合的状态, 水体物理化学参数从表层到底层无明显差异。春季分层开始出现时, 温跃层跨度较大, 为11 m。夏季和秋季的温跃层厚度相对狭窄(厚3—4 m)。显然, 随着热分层状态的变化, 水体物理化学参数也呈现梯度变化, 使得不同季节同一水层或同一季节不同水层的生境特征差异显著(表 1)。通常, 分层时期的水温、pH、溶解氧、叶绿素a在湖上层均显著高于湖下层(P < 0.05)。营养盐在不同水层差异也比较显著, 湖上层总氮、氨氮、总磷、正磷酸盐均低于湖下层。
变量 Variable |
混合期 Mixing layer (n=5) |
湖上层 Epilimnon (n=8) |
温跃层 Metalimnion (n=3) |
湖下层 Hypolimnion (n=4) |
水温Water temperature /℃ | 16.00±0.073a | 25.55±0.56b | 20.05±0.97c | 16.54±0.13a |
pH | 8.06±0.17a | 7.50±0.25a | 7.39±0.24a | 6.59±0.08b |
溶解氧Dissolved oxygen/(mg/L) | 8.82±0.40a | 6.79±0.56b | 4.95±0.75c | 0.34±0.09d |
叶绿素a Chlorophll a/(μg/L) | 5.17±1.02a | 4.39±0.65a | 3.23±0.76ab | 0.98±0.07b |
总氮Total nitrogen/(mg/L) | 1.35±0.02a | 1.06±0.03b | 1.32±0.06a | 1.30±0.06a |
氨氮NH 4-N/(mg/L) | 0.06±0.004a | 0.22±0.04ab | 0.10±0.03a | 0.35±0.09b |
硝态氮NO x-N/(mg/L) | 1.03±0.007a | 1.51±0.07b | 1.27±0.07c | 1.29±0.04c |
总磷Total phosphorus/(mg/L) | 0.05±0.001a | 0.02±0.005b | 0.02±0.003b | 0.06±0.006a |
正磷酸盐PO 4-P/(mg/L) | 0.02±0.001a | 0.004±0.001b | 0.01±0.005ab | 0.015±0.003a |
DNA拷贝DNA copies/L | (3.36±0.34) ×107a | (4.12±0.49) ×107a | (4.56±0.93)×107a | (5.10±0.51)×107a |
RNA拷贝RNA copies/L | (1.47±0.32)×105a | (1.10±0.23) ×105a | (6.30±1.93)×104a | (1.64±0.59)×105a |
RNA/DNA | 0.0043±0.00095a | 0.0023±0.00029ab | 0.0013±0.00022b | 0.0027±0.00079ab |
不同小写字母表示显著性差异(P < 0.05), 数据为平均值±标准误 |
固氮微生物的丰度和活性随着季节和水层的变化发生明显的变化(图 2)。冬季固氮微生物的nifH DNA丰度随着深度的增加有一个缓慢上升的趋势, 在底层达到最大值4.52×107拷贝/L;春季固氮微生物nifH DNA丰度在温跃层(5, 10 m)要低于湖上层和湖下层, 在湖下层15 m水深处为最大值4.95×107拷贝/L;夏季nifH DNA丰度在6—12 m急剧上升, 在温跃层开始点18 m达最高值7.66×107拷贝/L, 24 m稍微下降, 但并不明显;秋季nifH DNA丰度随着深度的增加而逐渐减少, 水深20 m水层丰度最低, 有趣的是nifH DNA丰度最大值位于秋季表层。固氮微生物nifHRNA水平变化趋势与DNA存在一定的差异性, 而且RNA丰度显著低于DNA丰度, 相差2个数量级。冬季水深10 m和15 m处的nifH RNA丰度要明显高于其他水层;春季随水深增加先降低后上升, 在10 m处为最低值;夏季, 随着深度的增加呈现明显上升的趋势, 在底层24 m达到最高值;秋季呈现总体下降的趋势, 在15—20 m急剧下降, 与nifH DNA变化趋势一致。RNA/DNA是一个能够表征基因表达活性的指标, 而4个季节表现出不同的垂直分布格局。例如, 4个季节活性最高出现的位置不同:10 m(冬季)、20 m(春季)、24 m(夏季)、15 m(秋季)。方差分析显示, nifH DNA和RNA在不同水层之间差异不显著, 而RNA/DNA在温跃层显著低于混合期水层(表 1)。显然, 固氮微生物群落的DNA、RNA丰度变化趋势呈现出一定差异性;基因表达活性的显著特征是在春季和夏季的底层活性最高, 而在秋季湖上层的活性显著高于湖下层。据此推测, 温跃层的出现是造成湖上层与湖下层差异的主要原因。
2.3 固氮微生物的多样性4个季节的表层和底层8个样品克隆文库共成功获得451条序列, 尽管稀释曲线没有平缓, 在97%相似度水平共鉴定106个固氮微生物OTUs(图 3)。主要类群包括:蓝藻(Cyanobacteria)46个OTUs、α-变形菌(Alphaproteobacteria)21个OTUs、β-变形菌(Betaproteobacteria)3个OTUs、γ-变形菌(Gammaproteobacteria)9个OTUs、厚壁菌(Firmicutes)1个OTU, 以及其他未知的固氮细菌18个OTUs、未知序列8个OTUs(表S1)。其中98个OTUs的序列与GenBank中各种生态坏境中的固氮微生物序列相似, 相似度水平范围是80%—99%。然而, 在50%相似性水平下, 20条序列(8个OTUs)在数据库中没有找到相似的序列, 因此属于未知序列。8个克隆文库中, 冬季混合期水库表层和底层OTUs数目并没有差异, 然而分层期水库表层和底层显示出明显的差异。例如, 春季和夏季底层OTUs数目高于表层, 而秋季表层OTUs数目稍高于底层(表 2)。在本研究测序范围内, 固氮微生物群落的多样性指数, 包括Shannon-Wiener指数、Simpson倒数指数、Pielou均匀度指数均呈现一定时空变化, 变化趋势与OTUs相似(表 2)。根据稀释曲线计算的OTU丰富度指数Chao1和ACE均表现出分层差异, 冬、春、夏3个季节水库底层大于表层, 而秋季却呈现为表层高于底层的特征(表 2)。
季节 Season |
水深/m Depth |
序列数 Sequences number |
OTU 0.03个数 OTU 0.03 number |
覆盖率/% Coverage |
香农威 纳指数 Shannon- Wiener index |
Simpson 倒数指数 Simpson inverse index |
Pielou均匀 度指数 Pielou evenness |
Chao1 指数 Chao1 index |
ACE指数 ACE index |
冬季Winter | 0 | 44 | 23 | 68.18 | 2.87 | 18.5 | 0.63 | 38 | 69 |
冬季 | 21 | 37 | 25 | 48.65 | 2.76 | 12.8 | 0.61 | 80 | 491 |
春季Spring | 0 | 67 | 12 | 91.04 | 1.57 | 3.03 | 0.44 | 17 | 36 |
春季 | 20 | 69 | 33 | 68.12 | 2.96 | 14.7 | 0.60 | 108 | 182 |
夏季Summer | 0 | 64 | 16 | 85.94 | 1.94 | 5.00 | 0.50 | 27 | 95 |
夏季 | 24 | 61 | 36 | 54.10 | 3.22 | 21.3 | 0.62 | 99 | 445 |
秋季Autumn | 0 | 66 | 17 | 87.88 | 2.19 | 6.67 | 0.55 | 25 | 40 |
秋季 | 25 | 43 | 12 | 88.37 | 1.70 | 3.13 | 0.47 | 14 | 16 |
除秋季底层蓝藻OTUs的比例为16.67%外, 在所有样品中蓝藻物种所占比例是最大的一类, OTUs占据38.89%—56.52%, 序列占据34.21%—76.92%(图 4)。α-变形菌OTUs为8.33%—27.78%, 其中秋季表层多样性最高, 序列丰度为2.33—10.29%;β-变形菌在冬季底层、春季表层、夏季表层, 以及秋季表层和底层均没有检测到, 而且OTUs和序列所占比例均小于5%;γ-变形菌在4个季节的表层与底层都被检测到, OTUs为4.00%—18.75%, 序列为2.70%—13.24%;厚壁菌仅仅在秋季表层被检测到, OTUs为5.55%, 序列为1.52%;未知固氮菌在4个季节中均被发现, 而且底层OTUs和序列数均高于表层;未知序列在夏季表层和秋季表层之外其他水层都有检测到, 而且底层均高于表层(图 4)。
聚类分析结果显示, 4个季节表层与底层的固氮微生物群落组成明显不同(图 5)。冬季表层与底层的群落组成比较相似, 差异相对较小。存在分层现象的春、夏、秋3个季节, 表层固氮微生物群落组成最相似;春、夏底层的群落组成较相似, 秋季底层显示出其独特性。据此推测, 水体分层对固氮微生物群落组成的影响比季节差异的影响更明显。
2.5 固氮微生物群落与环境因子的关系汀溪水库nifH DNA丰度与pH、溶解氧、叶绿素a显著负相关, 与硝态氮显著正相关;nifH RNA丰度则与氨氮存在极显著的负相关关系;RNA/DNA与水温显著负相关(表 3)。显然, 热分层发生时温跃层的出现及其导致的水体理化分层与固氮微生物群落组成、丰度和活性密切相关。
相关性Correlation | DNA | RNA | RNA/DNA |
水温Water temperature | NS | NS | -0.320* |
溶解氧Dissolved oxygen | -0.316* | NS | NS |
总氮Total nitrogen | NS | NS | NS |
硝态氮NO x -N | 0.382* | NS | NS |
正磷酸盐PO4-P | NS | NS | NS |
pH | -0.433** | NS | NS |
叶绿素a Chlorophyll a | -0.315* | NS | NS |
氨氮NH4-N | NS | -0.362** | NS |
总磷Total phosphorus | NS | NS | NS |
* P < 0.05; ** P < 0.01; NS:不显著 not significant |
汀溪水库在水体分层期间(春、夏、秋), 固氮微生物群落在上下水层空间上的差异要明显高于季节上的差异。在冬季处于混合期, 水体上下混合导致表层到底层物理化学因子比较接近, 同样表层与底层的固氮微生物群落组成也比较接近, 然而nifH DNA和RNA在底层的丰度却高于表层。利用nifH基因可以检测到广泛的自养和异养固氮微生物:蓝藻、古菌、厚壁菌、螺旋菌、α-变形菌、β-变形菌、γ-变形菌和δ-变形菌[1]。本研究同样检测到蓝藻、α-变形菌、β-变形菌、γ-变形菌、厚壁菌、未知固氮菌, 还有少量不能确定归属的新序列。在所有的已知类群中, 蓝藻序列占据最优势的地位。水库在春季开始出现分层, 固氮微生物群落组成在表层和底层就呈现出明显的差异, 例如春季底层OTU数目, 以及多样性指数均明显高于表层, 而且表层没有检测到β-变形菌。夏季相较于春季而言温跃层下降, 有趣的是夏季底层检测到的nifH DNA丰度、RNA丰度、RNA/DNA、OTU数目都明显高于表层, 说明底层固氮微生物丰度、活性、多样性高于表层。秋季温跃层位置继续下降, 表层OTU数目高于底层;而且表层nifH DNA丰度、RNA丰度、RNA/DNA明显高于底层, 呈现出与夏季相反的垂直分布格局。Moisander[9]等人发现在4月、10月切萨皮克湾的底层水体固氮微生物群落组成比较相似, 而表层样品变化更加明显, 这与本文的结果存在差异。汀溪水库4月、7月和10月表层水体的固氮微生物组成更加相似, 而底层样品多样性高, 群落组成变化较大。有趣的是, 在秋季表层检测到厚壁菌(Firmicutes)。厚壁菌微生物一般大量出现在土壤环境[26], 本文的检测结果暗示具有固氮作用的厚壁菌也能生存在水体环境或者来源于库区流域的土壤。由于在春夏表层没有发现β-变形菌, 在底层却检测到, 因此认为具有固氮功能的β-变形菌在分层水库中可能更加偏好底层的低温、缺氧的环境。
随着大量相关研究的进行, 固氮微生物的多样性和复杂性逐步被人们所认知[1], 但是影响水生态系统固氮微生物的分布及多样性的因子还不是很清楚[20]。在贫营养海洋中, 固氮微生物的丰度低于湖泊或者近海岸水体, 但是多样性却更高[1]。固氮蓝藻是湖泊和水库中常见的水华种类[27], 而且通常占据优势地位, 如纽约乔治湖[28]。有趣的是, 氮是莫诺湖的限制因子, 尽管在垂直水柱中没有检测到固氮作用的发生, 表层和底层水体却发现多样性很高且群落差异显著的固氮微生物[1]。河口具有多样的固氮微生物, 如切萨皮克湾与纽斯河具有丰富的α-变形菌、β-变形菌、γ-变形菌[1, 29]。同样, 在河口沉积物和沼泽中, 固氮微生物是非常多样的, 而且不论氮含量高低或者环境变化, 总是会保持相对稳定的状态[30-31]。综上所述, 固氮微生物种类在水环境中具有广泛的分布, 包括微生物垫、湖泊、盐沼、河口等多种生境, 而在不同环境下存在着各自的分布特点[20]。本文的研究结果在一定程度和前人的结果一致, 例如, 蓝藻是淡水生境优势类群;但也存在一定的特殊性, 比如检测到厚壁菌。值得注意的是, 本研究中8个OTUs在GenBank上没有找到比较接近的序列, 暗示这些OTUs可能是新的未知序列, 迫切需要后续的实验分析确定其分类地位信息。此外, 这些未知固氮菌和未知序列通常在4个季节底层的出现频率高于表层, 主要原因可能是过去对湖泊水库的研究往往忽视底层水体, 对于底层水体微生物的研究远远少于表层。
3.2 固氮蓝藻占优势的原因在汀溪水库的8个克隆文库中, 固氮微生物的多样性非常高, 更重要的特征是蓝藻占据绝对优势。湖泊水库蓝藻优势的存在有其物理、化学、生物因素等多重原因[32]。研究表明, 水库水体上下混合运动会影响物种间的相互作用, 在小范围内存在水体运动的环境中, 有竞争作用的物种也可以共同存在[33]。Zani[28]等人在乔治湖发现固氮蓝藻与其他固氮菌共同存在, 这与本文在同一样品中检测到多样的固氮微生物是一致的, 虽然有竞争作用, 但是可以共同生存。通常, 氧气对固氮酶活性具有抑制作用, 蓝藻占据优势地位也与自身生理结构有关[1]。例如蓝藻有两种固氮机制:空间上区分开光合作用和固氮作用;时间上错开光合作用和固氮作用[34]。因此, 蓝藻在水体表层和底层都可以进行固氮作用, 而且氧气不会对其固氮作用造成不利影响。有研究报道, 海洋蓝藻束毛藻(eg. Trichodesmium spp.)是热带贫营养盐海域优势的固氮微生物[35], 而具有异形胞的蓝藻在淡水和半咸水环境中都可以成为占优势的类群[36]。具有异形胞的蓝藻相对于非异形胞蓝藻更加适合固氮作用[35], 因为异形胞可以空间上分开光合作用与固氮作用, 这样不仅固氮效率更高, 而且也有效避免了氧气的抑制作用, 因此可以解释在氮磷比低的淡水湖泊和水库中固氮蓝藻容易形成水华[37-38]。变形菌门的固氮微生物多数是异养型固氮微生物, 相对喜欢底层低氧或厌氧环境[39]。Toepel[40]等人研究了Cyanothece sp. ATCC 51142株在白天黑夜及持续性光照条件下固氮基因表达差异性, 实验证明其在白天会储存糖原, 连续光照能够诱导晚上进行固氮作用。另一方面, 水库稳定的分层现象也为蓝藻提供了有利的条件, 蓝藻可以在分层水体上下自由移动, 这是因为蓝藻体内有特殊的气泡结构, 可以通过自身调节提供浮力[41]。蓝藻可以借助气泡在水层中找到光、营养盐等合适生境的水层, 进行生长繁殖[42]。由此可以解释, 在亚热带深水水库蓝藻不论季节和空间上的变化都可以占据优势地位。
Langois[43]等人发现蓝藻更偏好表层水体, 这与我们研究结果一致, 因为汀溪水库4个季节表层蓝藻的丰度比例和物种多样性均高于底层。蓝藻喜欢表层并不是一个因素决定的, 可能是蓝藻属于自养光合生物, 偏好较高水温、更适宜光照的表层水体, 或者是同其他生态因子的综合效应[27, 35]。Bentzon-Tilia[44]等人研究了两个河口湾的固氮蓝藻, 发现两个河口湾表层水体比海洋束毛藻(Trichodesmium)水华爆发期的固氮效率还要高。由此推测, 水库生态系统中蓝藻在固氮作用中扮演着很重要的角色。Voss[45]等人发现随着海水深度变化的总溶解铁与固氮作用是显著相关的。这需要后续的野外控制实验验证水库中微量元素溶解铁与固氮过程的关系, 进而可以更加深入的了解蓝藻固氮的生态过程与机制。另外一种原因, 可能是蓝藻体内有气泡、具有很强的漂浮能力, 也更有利于在表层占据优势地位[41-43]。本文研究结果提示, 在亚热带深水水库全年具有蓝藻水华爆发的条件和可能, 因此在蓝藻水华预防和控制时要充分考虑水体分层的变化规律。
4 结论本研究鉴定固氮微生物OTUs超过100个, 表明汀溪水库中具有多样性很高的固氮微生物, 主要类群包括蓝藻、α-变形菌、β-变形菌、γ-变形菌、厚壁菌, 以及少量未知的固氮细菌和序列。固氮微生物的丰度、活性、多样性及群落组成具有明显的季节和垂直分布格局, 在春夏秋水库热分层时期, 群落组成在水体表层和底层的空间变化明显大于时间变化, 暗示热分层对固氮微生物群落具有更为显著的影响。春夏秋冬4个季节里, 蓝藻在水库中始终占据最优势地位, 特别是在水库分层时期(春、夏、秋), 蓝藻的优势更加突出。高温和热分层环境, 以及蓝藻的独特生理结构有利于表层水体固氮蓝藻占优势。有趣的是, 汀溪水库夏季和秋季固氮微生物丰度和活性呈现相反的垂直变化趋势。在后续研究中, 有必要更加深入系统研究并揭示热分层条件下固氮微生物群落变化的生态机制, 进而为预防和控制蓝藻水华发生提供对策。
[1] | Zehr J P, Jenkins B D, Short S M, Steward G F. Nitrogenase gene diversity and microbial community structure:a cross-system comparison. Environmental Microbiology , 2003, 5 (7) : 539–554. DOI:10.1046/j.1462-2920.2003.00451.x |
[2] | Zehr J P, Paerl H W. Molecular ecological aspects of nitrogen fixation in the marine environment//Kirchman D L, ed. Microbial Ecology of the Oceans, Second Edition. New York:John Wiley & Sons, Inc., 2008:481-525. |
[3] | Gaby J C, Buckley D H. A comprehensive evaluation of PCR primers to amplify the nifH gene of nitrogenase. PLoS One , 2012, 7 (7) : e42149. DOI:10.1371/journal.pone.0042149 |
[4] | DeLong E F, Preston C M, Mincer T, Rich V, Hallam S J, Frigaard N-U, Martinez A, Sullivan M B, Edwards R, Brito B R, Chisholm S W, Karl D M. Community genomics among stratified microbial assemblages in the ocean's interior. Science , 2006, 311 (5760) : 496–503. DOI:10.1126/science.1120250 |
[5] | Short S M, Zehr J P. Nitrogenase gene expression in the Chesapeake Bay Estuary. Environmental Microbiology , 2007, 9 (6) : 1591–1596. DOI:10.1111/emi.2007.9.issue-6 |
[6] | Farnelid H, Öberg T, Riemann L. Identity and dynamics of putative N2-fixing picoplankton in the Baltic Sea proper suggest complex patterns of regulation. Environmental Microbiology Reports , 2009, 1 (2) : 145–154. DOI:10.1111/emi4.2009.1.issue-2 |
[7] | Young J P W. Phylogenetic classification of nitrogen-fixing organisms//Stacey G, Evans H J, Burris R H, eds. Biological Nitrogen Fixation. New York, USA:Chapman and Hall, 1992:43-86. |
[8] | Gaby J C, Buckley D H. A global census of nitrogenase diversity. Environmental Microbiology , 2011, 13 (7) : 1790–1799. DOI:10.1111/j.1462-2920.2011.02488.x |
[9] | Moisander P H, Morrison A E, Ward B B, Jenkins B D, Zehr J P. Spatial-temporal variability in diazotroph assemblages in Chesapeake Bay using an oligonucleotide nifH microarray. Environmental Microbiology , 2007, 9 (7) : 1823–1835. DOI:10.1111/j.1462-2920.2007.01304.x |
[10] | Church M J, Björkman K M, Karl D M, Saito M A, Zehr J P. Regional distributions of nitrogen-fixing bacteria in the Pacific Ocean. Limnology and Oceanography , 2008, 53 (1) : 63–77. DOI:10.4319/lo.2008.53.1.0063 |
[11] | Langlois R J, Hümmer D, LaRoche J. Abundances and distributions of the dominant nifH phylotypes in the Northern Atlantic Ocean. Applied and Environmental Microbiology , 2008, 74 (6) : 1922–1931. DOI:10.1128/AEM.01720-07 |
[12] | Yu Z, Zhou J, Yang J, Yu X Q, Liu L M. Vertical distribution of diazotrophic bacterial community associated with temperature and oxygen gradients in a subtropical reservoir. Hydrobiologia , 2014, 741 (1) : 69–77. DOI:10.1007/s10750-014-1832-6 |
[13] | Zehr J P, Church M J, Moisander P H. Diversity, distribution and biogeochemical significance of nitrogen-fixing microorganisms in anoxic and suboxic ocean environments//Neretin L N, ed. Past and Present Water Column Anoxia. Netherlands:Springer, 2006:337-369. |
[14] | Mahaffey C, Michaels A F, Capone D G. The conundrum of marine N2 fixation. American Journal of Science , 2005, 305 (6-8) : 546–595. DOI:10.2475/ajs.305.6-8.546 |
[15] | Dang H Y, Yang J Y, Li J, Luan X W, Zhang Y B, Gu G Z, Xue R R, Zong M Y, Klotz M G. Environment-dependent distribution of the sediment nifH-harboring microbiota in the northern South China Sea. Applied and Environmental Microbiology , 2013, 79 (1) : 121–132. DOI:10.1128/AEM.01889-12 |
[16] | Farnelid H, Bentzon-Tilia M, Andersson A F, Bertisson S, Jost G, Labrenz M, Jürgens K, Riemann L. Active nitrogen-fixing heterotrophic bacteria at and below the chemocline of the central Baltic Sea. The ISME Journal , 2013, 7 (7) : 1413–1423. DOI:10.1038/ismej.2013.26 |
[17] | 陈因, 方大惟. 蓝藻Anabaena 7120固氮的光调节. 植物生理学报 , 1983, 9 (1) : 51–59. |
[18] | Howarth R W, Marino R, Cole J J. Nitrogen fixation in freshwater, estuarine, and marine ecosystems. 2. Biogeochemical controls. Limnology and Oceanography , 1988, 33 : 688–701. |
[19] | Havens K E, James R T, East T L, Smith V H. N:P ratios, light limitation, and cyanobacterial dominance in a subtropical lake impacted by non-point source nutrient pollution. Environmental Pollution , 2003, 122 (3) : 379–390. DOI:10.1016/S0269-7491(02)00304-4 |
[20] | Short S M, Jenkins B D, Zehr J P. Spatial and temporal distribution of two diazotrophic bacteria in the Chesapeake Bay. Applied and Environmental Microbiology , 2004, 70 (4) : 2186–2192. DOI:10.1128/AEM.70.4.2186-2192.2004 |
[21] | Wang S, Qian X, Han B P, Wang Q H, Ding Z F. Physical limnology of a typical subtropical reservoir in south China. Lake and Reservoir Management , 2011, 27 (2) : 149–161. DOI:10.1080/07438141.2011.573613 |
[22] | Yu Z, Yang J, Zhou J, Yu X Q, Liu L M, Lv H. Water stratification affects the microeukaryotic community in a subtropical deep reservoir. Journal of Eukaryotic Microbiology , 2014, 61 (2) : 126–133. DOI:10.1111/jeu.2014.61.issue-2 |
[23] | Olson J B, Steppe T F, Litaker R W, Paerl H W. N2-fixing microbial consortia associated with the ice cover of Lake Bonney, Antarctica. Microbial Ecology , 1998, 36 (3-4) : 231–238. |
[24] | Schloss P D, Westcott S L, Ryabin T, Hall J R, Hartmann M, Hollister E B, Lesniewski R A, Oakley B B, Parks D H, Robinson C J, Sahl J W, Stres B, Thallinger G G, Van Horn D J, Weber C F. Introducing mothur:open-source, platform-independent, community-supported software for describing and comparing microbial communities. Applied and Environmental Microbiology , 2009, 75 (23) : 7537–7541. DOI:10.1128/AEM.01541-09 |
[25] | Clarke K R, Gorley R N. PRIMER v6:User Manual/Tutorial. Plymouth:PRIMER-E, 2006. http://www.scirp.org/reference/ReferencesPapers.aspx?ReferenceID=1328480 |
[26] | Kizilova A K, Titova L V, Kravchenko I K, Iutinskaya G A. Evaluation of the diversity of nitrogen-fixing bacteria in soybean rhizosphere by nifH gene analysis. Microbiology , 2012, 81 (5) : 621–629. DOI:10.1134/S0026261712050116 |
[27] | Paerl H W. Nuisance phytoplankton blooms in coastal, estuarine, and inland waters. Limnology and Oceanography , 1988, 33 (4part2) : 823–843. DOI:10.4319/lo.1988.33.4part2.0823 |
[28] | Zani S, Mellon M T, Collier J L, Zehr J P. Expression of nifH genes in natural microbial assemblages in Lake George, New York, detected by reverse transcriptase PCR. Applied and Environmental Microbiology , 2000, 66 (7) : 3119–3124. DOI:10.1128/AEM.66.7.3119-3124.2000 |
[29] | Affourtit J, Zehr J P, Paerl H W. Distribution of nitrogen-fixing microorganisms along the Neuse River Estuary, North Carolina. Microbial Ecology , 2001, 41 (2) : 114–123. |
[30] | Piceno Y M, Lovell C R. Stability in natural bacterial communities:Ⅰ. Nutrient addition effects on rhizosphere diazotroph assemblage composition. Microbial Ecology , 2000, 39 (1) : 32–40. DOI:10.1007/s002489900192 |
[31] | Bagwell C E, Lovell C R. Persistence of selected Spartina alterniflora rhizoplane diazotrophs exposed to natural and manipulated environmental variability. Applied and Environmental Microbiology , 2000, 66 (11) : 4625–4633. |
[32] | Shade A, Jones S E, McMahon K D. The influence of habitat heterogeneity on freshwater bacterial community composition and dynamics. Environmental Microbiology , 2008, 10 (4) : 1057–1067. DOI:10.1111/j.1462-2920.2007.01527.x |
[33] | Kerr B, Riley M A, Feldman M W, Bohannan B J. Local dispersal promotes biodiversity in a real-life game of rock-paper-scissors. Nature , 2002, 418 (6894) : 171–174. DOI:10.1038/nature00823 |
[34] | Kumar K, Mella-Herrera R A, Golden J W. Cyanobacterial heterocysts. Cold Spring Harbor Perspectives in Biology , 2010, 2 (4) : a000315. |
[35] | Karl D, Michaels A, Bergman B, Capone D, Carpenter E, Letelier R, Lipschultz F, Paerl H, Sigman D, Stal L. Dinitrogen fixation in the world's oceans. Biogeochemistry , 2002, 57/58 (1) : 47–98. |
[36] | Laamanen M J. Environmental factors affecting the occurrence of different morphological forms of cyanoprokaryotes in the northern Baltic Sea. Journal of Plankton Research , 1997, 19 (10) : 1385–1403. DOI:10.1093/plankt/19.10.1385 |
[37] | Walsby A E. The permeability of heterocysts to the gases nitrogen and oxygen. Proceedings of the Royal society of London. Series B. Biological Sciences , 1985, 226 (1244) : 345–366. DOI:10.1098/rspb.1985.0099 |
[38] | Staal M, Meysman F J R, Stal L J. Temperature excludes N2-fixing heterocystous cyanobacteria in the tropical oceans. Nature , 2003, 425 (6957) : 504–507. DOI:10.1038/nature01999 |
[39] | Jayakumar A, Al-Rshaidat M M D, Ward B B, Mulholland M R. Diversity, distribution, and expression of diazotroph nifH genes in oxygen-deficient waters of the Arabian Sea. FEMS Microbiology Ecology , 2012, 82 (3) : 597–606. DOI:10.1111/j.1574-6941.2012.01430.x |
[40] | Toepel J, Welsh E, Summerfield T C, Pakrasi H, Sherman L A. Differential transcriptional analysis of the cyanobacterium Cyanothece sp. strain ATCC 51142 during light-dark and continuous-light growth. Journal of Bacteriology , 2008, 190 (11) : 3904–3913. DOI:10.1128/JB.00206-08 |
[41] | Walsby A E, Hayes P K, Boje R, Stal L J. The selective advantage of buoyancy provided by gas vesicles for planktonic cyanobacteria in the Baltic Sea. New Phytologist , 1997, 136 (3) : 407–417. DOI:10.1046/j.1469-8137.1997.00754.x |
[42] | Walsby A E. Mechanisms of buoyancy regulation by planktonic cyanobacteria with gas vesicles//Fay P, Van Baalen C, eds. The Cyanobacteria. Amsterdam:Elsevier Science Publishers, 1987:377-414. |
[43] | Langlois R J, LaRoche J, Raab P A. Diazotrophic diversity and distribution in the tropical and subtropical Atlantic Ocean. Applied and Environmental Microbiology , 2005, 71 (12) : 7910–7919. DOI:10.1128/AEM.71.12.7910-7919.2005 |
[44] | Bentzon-Tilia M, Traving S J, Mantikci M, Knudsen-Leerbeck H, Hansen J L S, Markager S, Riemann L. Significant N2 fixation by heterotrophs, photoheterotrophs and heterocystous cyanobacteria in two temperate estuaries. The ISME Journal , 2015, 9 (2) : 273–285. DOI:10.1038/ismej.2014.119 |
[45] | Voss M, Croot P, Lochte K, Mills M, Peeken I. Patterns of nitrogen fixation along 10° N in the tropical Atlantic. Geophysical Research Letters , 2004, 31 (23) : L23S09. |