南京地区斑点叉尾鮰养殖池塘水体微生物群落结构研究

钟立强，王明华，张世勇，姜虎成，陈校辉，朱广伟，边文冀*

（1.江苏省淡水水产研究所，南京210017；2.中国科学院南京地理与湖泊研究所 湖泊与环境国家重点实验室 太湖湖泊生态系统研究站，南京210008；3.中国科学院大学，北京100049）

中国是全球最大的渔业养殖和水产品消费国家，水产养殖产量占世界水产养殖总产量的60%以上[1]。而内陆淡水池塘养殖又是其中最主要的养殖方式，2018 年养殖面积达252.7 万hm2，养殖产量2 212.2 万t，占全国水产品总产量的32.9%[2]。池塘养殖时，强调“养鱼先养水”，即鱼苗在下塘养殖前，要提前对养殖池塘的水体进行施肥等管理工作，使得养殖水体保持“肥、活、嫩、爽”。从生态学角度讲，这个过程就是通过施肥加速养殖水体富营养化，促进藻类及浮游动物的生长，从而提高水体生产力，为鱼类提供充足的生物饵料，形成相对稳定并适合鱼类生活的微生态环境。

与自然水域生态系统不同，淡水池塘养殖是特殊的人为操纵生态系统。充足的营养盐通过上行效应提高鱼类产量的同时，大量投入的有机质也为池塘生态系统食物网的腐屑链提供了能量和物质基础。细菌在腐屑链中占据着重要地位，一方面细菌起着降解作用，是物质和能量的贮存者，另一方面，它们是重要的生产者，能够将有机质转化为细菌菌体蛋白，直接被浮游动物和鱼类摄食[3]。已有研究表明，细菌群落与水体的营养水平及浮游植物的生物量存在明显关系[4-5]。养殖池塘中，由于富营养化和鱼类的牧食压力，小型浮游植物和浮游动物占据优势[6]，进而可能对池塘生态系统中的细菌群落结构造成影响。因此，对养殖池塘水体细菌群落结构开展研究十分必要。本文采用16S rRNA 基因片段高通量测序技术对南京斑点叉尾鮰（Ictalurus punctatus）淡水养殖池塘水体细菌群落的周年变化开展了研究，目的主要是：（1）了解斑点叉尾鮰养殖池塘水体细菌群落组成、多样性的年度变化特点；（2）明确斑点叉尾鮰养殖池塘水体细菌群落结构的主要影响因子。这些研究结果将对斑点叉尾鮰养殖池塘水体的精细化管理和健康养殖提供帮助。

1 材料与方法

1.1 养殖池塘及样品采集

试验在江苏省淡水水产研究所国家斑点叉尾鮰遗传育种中心养殖基地开展。养殖池塘面积0.21 hm2（50 m×42 m），位于南京市江宁区禄口镇（31°77′03″N，118°88′29″E）。养殖开始前一个月抽干池水，漂白粉清塘后，加水至1.5 m左右，池塘中央设置一台增氧机。4—10 月，每晚23：00 开启增氧机至次日凌晨4：00，晴天中午11：00 点开启2 h，保证养殖水体溶氧充足。夏季每月月中补水，春秋季则每2 月补水一次。池塘内放入2 747 尾斑点叉尾鮰鱼种（平均18.91 g·尾-1），养殖从11月底至次年12月初。每日早晚两次投喂漂浮性颗粒配合饲料（蛋白含量35%），饲料投喂量为鱼体质量的3%左右。

养殖周期内，每月的月末开展水样采集。每次上午10：00，利用Schindler 采水器（5 L）采集3处池塘表层50 cm 的水样，放入灭菌塑料瓶。水样放入4 ℃保温箱带回实验室，分析时样品为3 个点水样的混合样。取300 mL 混合水样经0.22 µm 的聚碳酸酯膜（Millipore，Cork，Ireland）负压过滤，得到滤膜样品，放入2 mL 的无菌离心管，-80 ℃超低温冰箱冷冻保存，留待DNA 提取。另取200 mL 混合水样4 ℃低温保存，进行水质的化学因子分析。

1.2 池塘水质因子测定

利用多功能水质参数仪（YSI6920v2，USA）现场测定温度（Temp）、酸碱度（pH）、总溶解性悬浮物（TDS）、浊度（Tur）以及溶解氧（DO）等环境因子。塞氏盘测定水体透明度（SD）。其他水质指标在实验室内按以下方法测定：总氮（TN）采用碱性过硫酸钾消解紫外分光光度法（GB 11894—1989）；铵态氮（NH-N）采用纳氏试剂分光光度法（GB 7479—1987）；硝态氮（NO-N）采用酚二磺酸分光光度法（GB 7480—1987）；亚硝态氮（NO-N）采用分光光度法（GB 7493—1987）；总磷（TP）和正磷酸盐（PO-P）采用磷钼蓝比色法（HJ 670-2013）；总有机碳（TOC）采用燃烧-非分散红外吸收法（GB 13193—1991）；化学需氧量（CODMn）采用高锰酸钾滴定法（GB 11892—1989）。水样经Whatman GF/C 0.45µm滤膜过滤后，分光光度法（HJ 897—2017）测定叶绿素a（Chl a）含量。

1.3 DNA提取和16S rRNA基因扩增

水样过滤收集滤膜，用预先经75%酒精擦拭过的医用剪刀剪碎后，利用Qiagen DNeasy Blood & Tissue Kit试剂盒（Qiagen，CA，USA）提取微生物总DNA。提取的DNA 经紫外分光光度计和1%的琼脂糖凝胶电泳检测浓度及纯度后，-20 ℃保存备用。以通用引物338F 和806R 扩增16S rRNA 基因V3-V4 高变区[7]，建库及高通量测序在上海美吉生物医药科技有限公司完成。

1.4 数据优化与分析

生物信息学分析前，先用Trimmomatic软件[8]去除测序原始序列的引物和barcode 标签，再用FLASH 软件[9]对序列进行拼接，使用VSEARCH（1.9.6）将序列以97%相似性聚类成操作分类单元（Operational taxonomic units，OTUs），比对的参考数据库是Silva_132 16SrRNA database（http：//www.arb-sliva.de）[10]。

所有统计分析及可视化分析均利用R 语言中的vegan和ggplot2程序包（Version 3.2.2）以及上海美吉生物医药科技有限公司I-sanger 云平台完成。多样性分析前，根据每个样本中最少的序列数进行抽平。利用mothur 软件计算α 多样性[11]。计算β 多样性前，对抽平后的物种数据进行Hellinger 转化[12]，再进行非度量多维尺度分析（Nonmetric multi-dimensional scaling，NMDS）[13]。对环境因子进行999次Monte Carlo置换检验的前向选择后，采用冗余分析（Redundancy analysis，RDA）对环境因子与微生物群落进行典范对应分析[14]。采用相似性分析（Analysis of similarities，ANOSIM）检验季节间差异。最后，通过多物种线性差异判别分析（Linear discriminant analysis effect size，LEfSe），根据分类学组成对样本按照不同分组条件进行线性判别分析（LDA），找出具有季节显著性差异的菌群[15]。

2 结果与分析

2.1 水体理化特征

年度养殖试验过程中，水温（Temperature）呈现出明显的季节性变化，从冬季的1月末逐步升高，夏季7月末达到峰值33.78 ℃，之后缓慢下降，除春季和秋季外，其他季节间水温都存在显著性差异（P＜0.05）。总溶解性悬浮物（TDS）和水体透明度（SD）都呈现出持续下降的趋势，且都是冬春季显著高于夏秋季（P＜0.05）。全年水样pH 介于7.68～9.50，始终保持偏碱性，冬季水体的碱性显著强于夏季（P＜0.05）。水体溶氧（DO）含量为6.12～16.39 mg·L-1，适合鱼类养殖，冬季水体的溶氧含量也显著高于其他3 个季节（P＜0.05）。水体总氮（TN）和总磷（TP）浓度始终维持较高水平，一直处于富营养化状态，四季间TN 变化不大，但冬季TP 浓度显著低于夏秋季节（P＜0.05）。氨氮（NH+4-N）、亚硝酸盐（NO-2-N）、硝酸盐（NO-3-N）随季节的变化表现出一定的差异，但正磷酸盐（PO3-4-P）则全年无显著性变化。总有机碳（TOC）也表现为冬春季显著高于夏秋季（P＜0.05）。叶绿素a（Chl a）夏秋季要高于冬春季，其中夏季显著高于其他3 个季节，并在7 月达到年度峰值，冬季的叶绿素a 浓度最低。水体理化特征见表1。

表1 斑点叉尾鮰养殖池塘水体理化性质Table 1 Summary of physicochemical parameters for water samples of channel catfish pond

2.2 微生物群落多样性特征

经优化处理后，共获得356 877条序列，平均每个样品有29 741 条高质量序列，序列平均长度为434.55 bp。这些序列被聚类为1 023 个OTUs。稀释曲线分析中所有样本的都接近平缓，表明本次试验抽样充分，测序深度足以反映每个样本细菌群落的生物信息。

选取Ace指数、Chao1指数、Shannon指数和Simp-son指数用来反映养殖池塘水体细菌群落的物种丰富度和多样性。结果显示（图1），养殖年度内，Ace 指数和Chao1指数变化趋势较为一致：冬春季呈S形变化，最小值出现在冬季2 月，并与相邻的1 月以及初春的3 月形成细菌群落物种丰富度的小低谷，之后的3 个月又持续上升；7 月Ace 指数和Chao1 指数大幅下降，之后锯齿状震荡上升，并在10 月达到年度最大值。使用单因素方差分析（ANOVA）对4 个季度的Ace 指数和Chao1 指数进行多重比较，结果表明，两个指数只有秋季和冬季之间存在显著性差异（P＜0.05）。年度内养殖池塘水体细菌群落的Shannon 指数和Simpson指数变化趋势呈相反的趋势，表明细菌群落多样性变化相对稳定，预测结果较为一致。Shannon 指数和Simpson指数为对应的“山”字状，7月是细菌群落多样性的最低值，且与其他月份差别较大，3月则是倒数第二低值，其他月份变化都不大（图1）。单因素方差分析与多重比较（ANOVA）结果显示，4 个季度Shannon指数和Simpson指数间无显著性差异（P＞0.05）。

图1 斑点叉尾鮰养殖池塘水体微生物群落α多样性Figure 1 The alpha diversity of bacterial community in channel catfish pond

使用非度量多维尺度分析（NMDS）展示养殖池塘水样细菌群落的β 多样性。NMDS 分布图（图2A）显示，所有样本呈现出较好的季节性排列趋势，即4个季度的样本可以被明显地区分，同一季度内的样品间距离更近，其微生物群落也更加相似。同时，整个年度水体细菌群落也呈现出连续的演替特征，从排序图的左侧连续往右排列，只有夏季样品的群落稍有例外，与秋季略有重叠，但整体排在右侧。NMDS 分析的Stress值为0.042，也表明排序结果可靠。细菌群落层次聚类分析的结果与NMDS类似，所有群落聚为两簇（图2B）：夏季和秋季样品（11 月样品除外）聚为一簇，冬季和春季的样品聚为另一簇。各分簇内部，每个季节的样品也聚类为一支。

2.3 微生物群落结构

1 023 个OUT 被注释为27 个门（Phyla）、60 个纲（Classes）和347个属（Genera）。为了减少群落的复杂性和冗余性，相对丰度低于1%的群落都归为稀有种，并为“Others”[16]。

养殖池塘水体样本中OTU 数量从高到低呈现为秋季＞春季＞夏季＞冬季。4个季节样本共享的OUT有324个，占全部OTU 的31.67%，3个季节和2个季节样品间共享OUT分别为274个和267个，占比26.78%和26.09%。而单季样品独有的OUT 数量最少，只有158个，占比15.44%。其中，冬季和秋季样品独有OUT 最多，分别为48个和47个，占比4.69%和4.59%，夏季样品独有OUT最少，仅29个，占比2.83%。

图2 斑点叉尾鮰养殖池塘水体微生物群落NMDS排序图和层级聚类分析Figure 2 NMDS plot and hierarchical clustering analysis of the compositional dissimilarities of the bacterial communities

ANOSIM 检验表明4个季度养殖池塘水体样品间细菌群落组成存在显著性差异（P＜0.01）。然而，分析细菌群落结构组成，所有样品在门水平上组成较为相似。相对丰度从高到低依次为：蓝藻门Cyanobacteria（29%）、放线菌门Actinobacteria（25%）、变形菌门Proteobacteria（20%）和 拟 杆 菌 门 Bacteroidetes（12%），这4 个门占全部群落的86%。而疣微菌门Verrucomicrobia（3.9%）、浮 霉 菌 门Planctomycetes（3.2%）、厚壁菌门Firmicutes（2.5%）、绿菌门Chlorobi（1.8%）和绿弯菌门Chloroflexi（1.4%）相对丰度都低于4%。

门水平上，冬季水体样品中放线菌门相对丰度最高，之后依次为变形菌门、蓝藻门以及拟杆菌门。而春、夏、秋3 季的水体样品中蓝藻门相对丰度最高，之后分别为放线菌门、变形菌门和拟杆菌门。各个季度样品内其他菌门的相对丰度都低于5%。从时间角度上看，养殖池塘水体中4 个主要菌门的相对丰度呈现出截然不同的变化趋势（图3）。蓝藻门的相对丰度呈单峰形变化，夏季达到最高峰（37%），春秋季略有下降（30%和32%），冬季最低（17.45%）。其中，3 月和7 月水体中蓝藻门相对丰度出现了两次明显的升高，分别为49.97%和56.77%。放线菌门相对含量则表现为增长的趋势，从春季的20%一直增加至冬季的29%。变形菌门1 年内的相对丰度保持相对稳定，维持在20%左右。拟杆菌门则呈U字形变化，即冬春季高（15%左右），夏秋季低（8.5%）。

4 个季度样品间的细菌群落结构比较显示，厚壁菌门、绿菌门和Saccharibacteria门存在显著性差异（P＜0.05）。其中，冬春季的厚壁菌门相对数量显著高于夏秋季，绿菌门相对数量则是夏秋季显著高于冬春季，而秋季Saccharibacteria 门的相对数量显著高于其他3个季节。

图3 斑点叉尾鮰养殖池塘水体微生门级水平细菌群落结构的年度变动Figure 3 Annual dynamic of bacterial community composition at the phylum level

属水平上，养殖池塘水体的优势菌属为hgcI_clade、蓝细菌属norank_c__Cyanobacteria、聚球藻属Syn⁃echococcus、CL500-29、CL500-3和多核杆菌属Polynu⁃cleobacterium。聚类热图表明，属水平上，所有样品聚为冬春和夏秋两簇（图4）。其中，肠球菌属Enterococ⁃cus、芽孢杆菌属Bacillus、乳球菌属Lactococcus、单胞菌属Polaromonas、红育菌属Rhodoferax和黄杆菌属Flavobacterium的相对丰度在冬春季明显高于夏秋季。而甲基嗜酸菌属Candidatus_Methylacidiphilum、球衣菌属Sphaerotilus、玫瑰单胞菌属Roseomonas、不动杆菌属Acidibacter和湖沉积杆菌属Limnobacter的相对丰度则是夏秋季明显高于冬春季。

2.4 季节差异性菌群分析

图4 养殖池塘水体微生物属级水平的分层聚类热图Figure 4 Bacterial community heatmap analysis in genus level

4个季度样品间的细菌群落结构比较在3个菌门发现了显著性差异，在时间维度上进一步采用LefSe软件对养殖水体细菌群落样本进行线性判别分析（LDA），以找到各个季节显著性对应的群落或物种。对LDA 判别值大于3.5 的细菌类群进行分析，共有不同分类水平的53个细菌类群，在4个季节间具有显著差异（图5）。春、夏和秋3 个季节的显著优势种相对集中，主要分布于一两个门。其中，春季的水体中，优势性菌群主要集中于厚壁菌门中的肠球菌属Entero⁃coccus，以及蓝藻门中的SubsectionIV。夏季水样中，优势性菌群集中于变形菌门中的δ-变形菌纲Deltaproteobacteria，寡弯菌属Oligoflexaceae，以及变形菌门中的球衣菌属Sphaerotilus和蓝藻门中的颤藻属Oscillatoria。秋季的样品优势菌群则主要集中于绿菌门中的OPB56 属、Saccharibacteria 门，以 及 蓝 藻 门 中 的 属FamilyI和变形菌门中的红螺菌目Rhodospirillales 与玫瑰单胞菌属Roseomonas。冬季具有显著优势菌群最多，共有2 个纲、6 个目、9 个科和10 个属。其中，疣微菌门中的疣微菌纲Verrucomicrobiae 疣微菌科Verrucomicrobiaceae、丰佑菌纲OpitutaeVadin中的HA64属、放线菌门中的PeM15 纲、从棒状杆菌纲Corynebacteriales中的分枝杆菌属Mycobaterium、厚壁菌门中的芽孢杆菌目Bacillales 芽孢杆菌属Bacillus、链球菌科Streptococcaceae 乳球菌属Lactococcus、变形菌门中的红细菌目Rhodobacterales 到红杆菌属Rhodobacte⁃raceae、MNG 科、以及其他零星分布的普雷沃氏菌科Prevotellaceae、球形红杆菌属Candidatus_Aquirestis、嗜甲基菌属Methylotenera、单胞菌属Polaromona和Limnohabitans属。

图5 斑点叉尾鮰养殖池塘四季间差异性菌群多物种线性差异判别分析Figure 5 LEfSe identified the most differentially abundant taxa in the channel catfish pond among four seasons

2.5 细菌群落变化影响因素

与样本菌群组成相关的环境因子很多，但其中部分因子存在共线性，所以先进行筛选，保留相互作用较小的环境因子。方差膨胀因子（Variance inflation factor，VIF）分析表明自变量间的多重共线性关系严重，过滤掉VIF 大于20的环境因子NO-3-N 和Chl a，多次筛选后，剩下的环境因子对应的VIF 值全部小于20。Monte Carlo 排序检验表明Temp、TDS、pH、NO2-N、TOC与细菌群落之间存在显著的相关性（P＜0.05）。这5个因子对细菌群落格局的影响大小（r2）依次顺序是：TDS＞Temp＞pH＞NO-2-N＞TOC（表2）。

利用前置检验挑选出来的5 个与细菌群落显著相关的环境因子进行RDA 排序分析，结果表明，RDA1 轴和RDA2 轴分别解释了细菌群落变异的44.53%和10.98%，主要菌群被聚为2 簇（图6）。第1簇由冬季和春季样品组成，第2 簇由夏季和秋季样品组成，两簇被RDA2轴分隔。其中，第1簇冬春季样本分布比较紧密，主要受pH、TOC和TDS的影响；第2簇夏秋季样本分布较为分散，且夏秋季样品被RDA1 轴分 隔，主 要 受NO-2-N 和Temp 的 影 响。RDA 结 果 显示，斑点叉尾鮰养殖池塘水体微生物具有比较明显的季节性，4个季节的样品沿着RDA1轴从左往右排列。

3 讨论

本文通过分析斑点叉尾鮰养殖池塘水体细菌群落的年度变化过程，阐明斑点叉尾鮰养殖池塘水体细菌群落的变化规律。研究表明，养殖池塘的微生物群落与营养盐、水质、养殖对象及日常管理都具有重要的关系[3]。本研究结果表明，斑点叉尾鮰养殖池塘水体优势菌群为蓝藻门、放线菌门、变形菌门和拟杆菌门，与三疣梭子蟹和吉富罗非鱼养殖塘水体较为相似，但是主要优势菌群的相对丰度排序略有差别[17-18]。不同养殖品种，甚至同一养殖品种的不同养殖阶段，其养殖水体的微生物群落都有明显的差别，这与养殖对象的栖息方式、食性特征以及养殖周期等差异有关[18]。斑点叉尾鮰和吉富罗非鱼养殖过程中投饵量大，产生大量的残饵和粪便，往往造成过量的营养物输入至养殖池塘，同时，斑点叉尾鮰和罗非鱼都是底栖性鱼类，因此优势细菌群落较为相似。

ANOSIM 检验表明斑点叉尾鮰养殖池塘水体4个季度样品间细菌群落组成存在显著性差异（P＜0.01），NMDS 分析（图2A）也验证了所有样本呈现出较好的季节性排列趋势。斑点叉尾鮰养殖池塘水体优势菌中，拟杆菌门相对丰度冬春季显著高于夏秋季节，LefSe 分析进一步表明，主要差异菌群为冬季的普雷沃氏菌科Prevotellaceae、球形红杆菌属Candida⁃tus_Aquirestis。普雷沃氏菌科Prevotellaceae 能够利用和降解蛋白质及各类有机物[19]，而斑点叉尾鮰养殖池塘的总有机碳显著高于夏秋季节，为普雷沃氏菌科细菌的繁殖提供了适合的环境。球形红杆菌属Candi⁃datus_Aquirestis广泛分布于内陆高硬度水体，冬春季也会在自然水体中形成明显的生物量高峰，秋冬季则维持较低的丰度[20]，与斑点叉尾鮰养殖池塘水体的明显季节性变化一致。

表2 环境因子对排序结果的相关性系数及显著性检验值Table 2 Coefficient correlations between physicochemical parameters and RDA axis

图6 细菌群落与环境因子之间的RDA分析Figure 6 RDA bioplot of the distribution of bacterial community and environmental variables

蓝藻门作为广泛分布于各类水体的细菌门类，其丰度状况与水体富营养化程度密切相关。水产养殖环境中，因为营养盐和有机物的密集投入，水体长期维持富营养化状态，蓝藻水华爆发已经变得越来越频繁[21]。斑点叉尾鮰池塘的蓝藻含量虽然没有发现显著的季节性差异，但是其相对丰度依然存在季节性的变化规律。春季3月和夏季7月斑点叉尾鮰养殖池塘出现了两次蓝藻生物量的高峰，蓝藻生物量占总细菌生物量的49.98%和56.77%，同时，3 月和7 月蓝藻生物量的上升导致池塘水体细菌群落α 多样性显著降低。LefSe分析表明，春、夏和秋季蓝藻门的显著优势菌群各不相同。春季优势蓝藻主要是念珠藻目的g__norank_f__FamilyI_o__SubsectionIV属，大 都 呈 丝状，形成厚壁孢子，在长换水周期、高pH、高温和高光照强度等条件下，更具有优势[22]。斑点叉尾鮰养殖池塘水体没有剧烈波动，冬春季pH 一直较高，3 月末温度已经接近15 ℃，营养盐也一直维持较高浓度，因此在3月末，蓝藻迅速生长，并导致了叶绿素a浓度的上升。之后4 月蓝藻大量死亡并氧化分解，消耗水体溶氧并产生乳酸等酸性物质，导致4 月池塘溶氧含量降低，pH 也略有下降。夏季优势的颤藻属Oscillatoria是常见的淡水蓝藻群体，呈丝状不分枝，会节律性颤动。斑点叉尾鮰养殖池塘水体夏季颤藻属快速生长，成为优势类群，并在7 月形成蓝藻水华的最高峰，此时池塘水体叶绿素a 浓度也是全年最高，这个过程与开放湖库的夏季蓝藻水华演替较为相似[23]。秋季g__norank_f__FamilyI_o__SubsectionI属的蓝藻替代颤藻属成为优势蓝藻类群，该属归属于蓝藻纲Subsection I，皆为单细胞蓝藻，有空泡，并合成藻胆蛋白，是湖泊、水库等天然水体中秋季常见的水华蓝藻种类[24]。

Monte Carlo 检验和RDA 排序分析表明斑点叉尾鮰养殖池塘水体细菌群落与环境因子间具有显著的相关性，主要的环境因子包括TDS、Temp、pH、NO-N和TOC。5 个环境因子都表现出明显的季节性差异，冬春季的TDS、pH 和TOC 显著高于夏秋季，Temp 和NO-N 则是夏秋季显著高于冬春季。其中，TDS 和Temp 对细菌群落的影响最大。TDS 能干扰水体的氧传递和降低细菌活性，从而影响硝酸盐、有机物和COD 的去除[25]。大量的研究都表明温度对细菌群落组成有重要影响。温度会影响细菌降解有机物、再生营养以及参与微食物链等方面的生物学功能[26-27]，同时还对水体浮游植物的演替产生影响，改变水体微生态系统[28]，从而驱动细菌群落组成的改变。水温季节性的波动造成栖息环境的季节性变化，进一步推动水体细菌群落组成呈现出季节性的响应[29]。此外，不同菌群对温度的敏感差异，也会导致不同温度下细菌群落结构的差异[30]。斑点叉尾鮰养殖池塘春季和夏季的水温升高提高了蓝藻门的生长速率，并给蓝藻门带来了竞争优势[31]，最终与其他因素共同导致斑点叉尾鮰养殖池塘3月和7月的两次蓝藻爆发。

pH 被认为是生物化学地球转换过程中重要的调控因子，它能调节水体中离子和微量金属的利用率，从而对细菌的生长起到抑制或者促进的作用[32]。养殖水体光合作用的强弱对pH影响较大，因此pH会随着光照强度和初级生产力的改变呈现出明显的季节性变化。本研究中pH 冬春明显高于夏秋，这可能与本试验池塘水体初级生产力相对较低有关，叶绿素a浓度明显低于草鱼、乌鳢和罗非鱼等品种的养殖池塘[33-34]。同时，夏季池塘补水量大，可能也会引起pH的降低。N 和P 一直被认为是细菌和浮游植物生长的限制性因子[35]，TOC 也被认为能够直接影响细菌的生物量和组成，从而改变细菌群落结构[36]。然而，斑点叉尾鮰养殖水体一直呈人为的富营养化状态，大部分的营养盐（TN、NH-N、NO-N、TP 和PO-P）都没有对细菌群落结构的塑造起到显著性作用。而NO-N是水体中氮素存在的一个中间态，本研究池塘溶氧持续充足，整个养殖年度内NO-N 浓度一直维持在较低的状态。因此，NO-N的低浓度加上夏秋高冬春低的波动，会对细菌群落如亚硝酸盐硝化细菌等产生显著的影响。斑点叉尾鮰养殖池塘水体的C/N[可由TOC（/NH-N+NO-N+NO-N）粗略计算]只有7～14，远低于细菌生长的最适范围20～30[37]，因此C（主要是TOC）是斑点叉尾鮰养殖池塘细菌生长的限制因子。拟杆菌门的菌群能够高效利用低碳物质[38]，冬春季的TOC 浓度显著高于夏秋季节（P＜0.05，表1），因此冬春季的拟杆菌门Bacteroidetes 丰度也显著高于夏秋季。

4 结论

本研究运用16S rDNA 扩增技术研究斑点叉尾鮰养殖池塘微生物群落，结果表明斑点叉尾鮰养殖池塘水体细菌群落结构呈现出季节性变化，优势菌群为蓝藻门、放线菌门、变形菌门和拟杆菌门。环境因子中，总溶解性悬浮物、水温、酸碱度、亚硝酸盐和总有机碳对水体细菌群落结构有显著的影响。但是本研究局限于一个地方的斑点叉尾鮰池塘主养模式，而斑点叉尾鮰养殖已经遍布全国，因此，后续研究仍需进一步针对斑点叉尾鮰的不同养殖模式在不同地域开展大尺度的养殖池塘微生物群落结构比较，从而为不同地域不同模式的斑点叉尾鮰养殖管理提供精确理论指导。