Spatio-temporal changes of bacterioplankton communities in Litopenaeus vannamei desalinated ponds and their responses to physicochemical factors
-
摘要: 浮游细菌调控是对虾养殖水体环境控制策略的核心内容,探究浮游细菌群落构建的一般规律,可进一步推动对虾养殖水体水质调控技术的研究。运用16S rRNA高通量测序技术,对凡纳滨对虾 (Litopenaeus vannamei) 淡化养殖池塘进行了12次周际调查。结果显示,48个样品共获得2 854个操作分类单元 (Operational taxonomic units, OTU,97%相似性),序列比对发现古细菌2门1纲1科1属,细菌30门59纲98目199科433属,其中优势菌群25属。优势菌群在组成上有较高的相似性,但各池优势菌群在分布和相对丰度变动上有较大差异。各池系统发育多样性指数总平均值为77.57,变幅为24.39~111.65;香农多样性指数总平均值为3.96,变幅为2.64~5.06;物种丰富度指数总平均值为716,变幅为229~1 054。非度量多维标度分析 (Non-metric multidimensional scaling, NMDS) 表明各池塘浮游细菌群落在养殖初期差异较大,中、后期差异减小,冗余分析 (Redundancy analysis, RDA) 显示活性磷、碱度、溶解氧和硫化物可显著影响浮游细菌的群落结构。Abstract: The regulation of bacterioplankton is the core content of the environmental control strategy in shrimp aquaculture. Exploring the general rules of the construction of bacterioplankton community can further promote the research on water quality regulation for shrimp ponds. Using 16s rRNA high-throughput sequencing technology, we conducted 12 weekly surveys in Litopenaeus vannamei desalinated ponds.The results indicate that a total of 2 854 OTUs (97% similarity) were obtained from 48 samples (Archaea belonged to 2 phyla, 1 class, 1 family and 1 genus; and bacteria belonged to 30 phyla, 59 classes, 98 orders, 199 families and 433 genera, among which 25 genera were dominant flora). The dominant flora had high similarity in the composition, but with great differences in the distribution and relative abundance in each pond. The total average phylogenetic diversity index was 77.57, ranging from 24.39 to 111.65; the total average Shannon diversity index was 3.96, ranging from 2.64 to 5.06; the total average species richness index was 716, ranging from 229 to 1 054. NMDS analysis shows that the community structure of bacterioplankton in each pond varied greatly at the early stage of aquaculture, but not so greatly at the middle and late stages. The results of redundancy analysis shows that the labile phosphorus, total alkalinity, dissolved oxygen and sulfide were the main environmental factors affecting the distribution characteristics of community structure of bacterioplankton.
-
斑节对虾 (Penaeus monodon) 适盐度广、生长迅速、出肉率高,是我国和世界三大主要养殖对虾之一[1]。目前我国斑节对虾的养殖模式主要有池塘高密度养殖、虾蟹混养和虾鱼混养等方式[2]。虾蟹混养是一种高空间利用率、高收益的虾类和蟹类2种不同甲壳动物混合养殖的复合生态养殖模式[3]。蟹类主要是锯缘青蟹 (Scylla serrata),由于蟹类为杂食性,能够摄食病弱或死亡的鱼虾、小型底栖动物等,对于对虾残饵的利用和水质净化有一定作用,可促进养殖塘底层物质循环和能量流动,提高营养物质的利用率[4]。虾蟹混养是一种复合生态养殖模式,因其持续稳定的养殖产量而备受欢迎。
浮游植物是指在水体中生活的微小植物,大小在几微米至几百微米之间,是水生生态系统中最主要的初级生产者,是食物网的基础[5]。浮游植物通过光合作用将水体中无机物转换为有机物,是水体中动物的主要食物来源,并且能够提高养殖水体中溶氧 (DO) 含量、吸收水体中的有毒有害物质、降低氨氮含量,保持养殖水体的稳定,其群落结构的变化对水生生态系统有着重要影响,主要包括能量流动、物质循环和信息传递[6]。浮游植物的生长繁殖极易受到外界环境的影响,其对环境条件的变化能够迅速做出响应,因此被视作水生生态系统的重要指示物[7-8]。
浮游植物作为饵料生物的基础,在池塘养殖过程中发挥了重要作用,与其相关的研究也日益受到重视。李俊伟等[9]研究对虾池塘混养鲻鱼 (Mugil cephalus) 和罗非鱼 (Oreochroms mossambcus) 对水环境及对虾生长的影响,探讨了浮游植物的变化特征;孔谦[10]也对凡纳滨对虾 (Litopenaeus vannamei) 与鲻鱼混养中精养池的理化因子进行了研究。目前有关虾蟹混养的研究较少,特别是针对混养池塘浮游植物群落结构的研究更是鲜见报道,本研究对广州市番禺区3个虾蟹混养池塘进行逐月调查,分析和探讨了养殖过程中浮游植物种类组成及其数量变化状况,旨在摸清虾蟹混养池塘浮游植物种群结构特征及其与环境因子的关系,以期为虾蟹混养的健康发展提供理论依据。
1. 材料与方法
1.1 养殖池塘
在广州市番禺区十六涌选取3个虾蟹混养池塘 (1#、2#、3#),水深2 m,3月初开始投放虾苗,每20 d投放1次;4月初开始投放蟹苗,每2个月投放1次。虾蟹分别投放4次,6月初开始捕获虾蟹。2018年3—11月对3个混养池塘进行逐月调查,采集池塘水质样品和浮游植物样品,带回实验室测定和鉴定。
1.2 样品的采集、测定和鉴定
1.2.1 水质样品
1) 样品采集和保存。采用容积为5 L的有机玻璃采水器,在每个池塘四角和中央采集表层水,用滤膜现场过滤,用玻璃瓶收集适量过滤水保存,带回实验室测定。2) 样品测定。营养盐的分析测定参照国家标准《海洋调查规范第2部分:海洋水文观测》 (GB/T 12763.2—2007)、《海洋调查规范第4部分:海洋化学要素观测》 (GB/T12763.4—2007) 和《海洋监测规范第4部分:海水分析》 (GB 17378.4—2007);叶绿素a的分析测定参照《海洋监测规范第7部分:近海污染生态监测和生物监测》 (GB 17378.7—2007) 中规定的方法。
1.2.2 浮游植物样品
1) 样品采集和保存。在每个池塘四角和中央用5 L采水器各采水样一份,经浮游生物网过滤后混合装入1 L采样瓶,并现场加入适量的5%甲醛溶液固定保存,带回实验室[11]。在实验室静置24 h,弃去上清液,浓缩后于高倍镜下鉴定浮游植物种类,并进行数量分析。2018年9月由于突发事件未采集到样品。2) 样品鉴定。采用显微镜鉴定。浮游植物样品定性分析和种类鉴定方法依据《中国淡水藻类》[12]、《中国海藻志》[13]、《中国淡水藻类》[14]、《中国淡水生物图谱》[15]、《环境微生物图谱》[16]。在10×40倍光学显微镜下进行,每次取0.1 mL样品放入浮游植物方格计数板计数,观察计数格,每个样品计数2片,2次求平均值作为最终结果,2片计数结果差异不超过15%,若相差较大,则再取1次计数,然后取数值相差较小的2片取平均值,即为浮游植物细胞密度[17]。
1.3 分析和评价方法
浮游植物物优势度 (Y) 应用以下公式计算:
$$ Y=\dfrac{{{n_i}}}{N}\times{f_i} $$ 式中ni为第i种的个体数;fi是该种在各站中出现的频率;N为所有站每个种出现的总个体数。
Shannon-Weaver多样性指数的计算公式为:
$$H' = - \sum\limits_{i = 1}^S {{P_i}{{\log }_2}{P_i}} $$ 式中H′为种类多样性指数,S为样品中的种类总数,Pi为第i种的个体数与总个体数的比值。
数据分析、处理以及图表的制作使用Excel 2016、Origin 2017、Canoco 5等软件。
2. 结果
2.1 种类组成
浮游植物调查结果见附录1 (详见http://dx.doi.org/10.12131/20190233的资源附件)。共鉴定出浮游植物5门29种,其中硅藻种类最多 (19种),占浮游植物种类数的65.5%;其次是蓝藻,共5种,占浮游植物种类数17.2%;甲藻3种,绿藻和裸藻各1种。每次采样中各塘浮游植物出现频次最多的有2门5种,其中硅藻门4种,分别为尖布纹藻 (Gyrosigma acuminatum)、新月菱形藻 (Nitzschia closterium)、柔弱菱形藻 (N. delicatissima) 和直舟形藻 (Navicula directa),出现的频率分别为10、10、22和9;蓝藻门1种,为小颤藻 (Oscillatoria tenuis),出现的频率为10 (图1-a)。
各池塘的浮游植物种类在时间上的差异不明显,基本呈先波动下降后小幅度增加的趋势,1#塘浮游植物种类数波动较大。各池塘出现种类数最多的月份各不相同 (图2-a)。其中1#塘出现种类数最多的为5月 (10种);2#和3#塘种类数最大值均出现在养殖初期 (分别为6种和8种)。各月份均以硅藻门种类数最多,其次为蓝藻门和甲藻门,绿藻门和裸藻门的种类仅在养殖后期出现。
2.2 细胞密度
各池塘浮游植物细胞密度的变化趋势见图2-b。浮游植物细胞密度介于2.56×103~189.2×103个·m−3,最大值出现在养殖前期的3—4月,最小值在7—8月。1#塘4月细胞密度最大 (124.8×103个·m−3),其次是5月 (52.92×103个·m−3),7月最低 (6.0×103个·m−3);2#塘3月细胞密度最大 (189.2×103个·m−3),其次是4月 (119.6×103个·m−3),8月最低 (2.56×103个·m−3);3#塘3月细胞密度最大 (118.8×103个·m−3),其次是4月 (82.13×103个·m−3),10月最低 (7.2×103个·m−3)。整个调查期间,3个池塘浮游植物细胞密度的最高值或次高值均出现在3—4月,而较低值均出现在7—8月 (图2-b)。
各次调查浮游植物细胞密度的组成中均以硅藻门占较高的比例,其他4门藻类所占比例相对较低;各养殖池塘各门浮游植物细胞密度所占百分比变化波动不大,同一月份不同养殖塘浮游植物细胞密度存在较大差异 (图 1-b)。
2.3 优势种组成
各次调查中浮游植物出现的优势种见表1。其中以7月和11月2次调查出现的优势种最多 (均为6种);3月和5月出现的优势种最少 (均为3种)。调查期间各混养池塘出现的优势种有硅藻门的奇异菱形藻、直舟形藻、柔弱菱形藻和尖布纹藻等;甲藻门的大角角藻;蓝藻门的小颤藻。其中柔弱菱形藻在8次调查中均为优势种,尖布纹藻、新月菱形藻和小颤藻在4次调查中为优势种,奇直舟形藻均在3次调查中为优势种,其他优势种在1~2次调查中为当次调查的优势种。由此可知,养殖初期优势种主要为硅藻门的种类,养殖中后期蓝藻门的种类优势度逐渐增加,成为优势种。
表 1 混养池塘浮游植物优势种Table 1. Dominant species of phytoplankton in polyculture pond优势种
Dominant species3月
Mar.4月
Apr.5月
May6月
Jun.7月
Jul.8月
Aug.10月
Oct.11月
Nov.奇异菱形藻 Nitzschia paradoxa 0.411 0.079 直舟形藻 Navicula directa 0.336 0.216 0.06 柔弱菱形藻 Nitzschia delicatissima 0.180 0.528 0.628 0.056 0.221 0.108 0.400 0.085 近缘斜纹藻 Pleurosigma affine 0.037 大角角藻 Ceratium macroceros 0.145 尖布纹藻 Gyrosigma acuminatum 0.445 0.197 0.028 0.026 小颤藻 Oscillatoria tenuis 0.060 0.134 0.386 0.307 寒带菱形藻 Nitzschia frigida 0.032 新月菱形藻 Nitzschia closterium 0.021 0.081 0.101 0.053 二形栅藻 Scenedesmus dimorphus 0.081 0.127 梭形裸藻 Euglena acus 0.020 银灰平裂藻 Merismopedia glauca 0.020 琼氏圆筛藻 Coscinodiscus jonesianus 0.035 平滑双眉藻 Amphora lavis 0.028 伏氏海毛藻 Thalassiothrix frauenfeldii 0.222 大颤藻 Oscillatoria maxima 0.047 2.4 Shannon-Weaver多样性指数
各混养池塘浮游植物多样性指数的变化趋势见图2-c。调查期间1#塘多样性指数介于0.802~2.140,在5月达到最高值,其次是3月,最低出现在6月;2#塘多样性指数介于0.998~2.247,在11月达到最高值,其次是7月,最低出现在8月;3#塘多样性指数介于0.468~2.597,在4月达到最高值,其次是5月,最低出现在6月。
各池塘多样性指数的变化趋势各不相同,1#塘4—6月的多样性指数呈先上升后下降的趋势;2#塘3—6月的多样性指数变化幅度较小,7—11月则呈先上升后下降再上升的趋势;3#塘4—11月的多样性指数呈先下降后上升的趋势。
各养殖塘内浮游植物多样性指数总体上表现为养殖前期低养殖后期高。根据马建新等[18]的浮游植物多样性指数评价标准,1#和2#塘全年的生物多样性都处于较好的范围内,3#塘养殖后期生物多样性水平较差 (图2-c)。
2.5 浮游植物与主要环境因子的关系
通过对浮游植物数据去趋势分析 (DCA),发现排序轴长为2.561,小于3,浮游植物群落分布更接近线性模型,根据之前学者的研究内容[19-20],本研究选取水温、pH、DO、总氮 (TN)、盐度、亚硝酸盐 (NO2 −-N)、硝酸盐 (NO3 −-N)、磷酸盐 (PO4 3−)和铵盐 (NH4 +) 9种环境因子和主要浮游植物优势种的细胞密度进行冗余分析 (RDA),浮游植物代码和环境因子分见表2和表3。养殖池塘浮游植物和环境因子关系的RDA排序图结果见图3。轴1和轴2特征值分别为0.263和0.18,前2个环境因子与物种轴之间的相关系数分别为0.939和0.922,拟合度较好,排序轴可以较好地反映浮游植物优势种和环境因子之间的关系;物种-环境累计百分比分别为36.1%和60.9% (表4)。养殖塘中浮游植物细胞密度主要受到营养盐、盐度和pH的影响,不同藻类受到环境因子的影响差异显著 (图3)。本次调查中,B1 (小颤藻) 作为蓝藻门的主要优势种,其细胞密度与水温、pH和DO成显著正相关,水温和pH的影响较大;A2 (柔弱菱形藻)、A4 (奇异菱形藻)、A3 (直舟形藻) 和A5 (近缘斜纹藻) 的细胞密度和NO2 −-N、NH4 +、PO4 3−、NO3 −-N、TN和盐度呈显著正相关,NO3 −-N和盐度对硅藻门影响较大C1 (大角角藻) 细胞密度与水温和PO4 3−呈显著正相关。
表 2 冗余分析中浮游植物代码种类Table 2. Phytoplankton species codes for redundancy analysis代码
Code种类
Species代码
Code种类
SpeciesA1 尖布纹藻
Gyrosigma acuminatumA6 中肋骨条藻
Skeletonema costatumA2 柔弱菱形藻
Nitzschia delicatissimaA7 洛伦菱形藻
Nitzschia lorenzianaA3 直舟形藻
Navicula directaB1 小颤藻
Oscillatoria tenuisA4 奇异菱形藻
Nitzschia paradoxaC1 大角角藻
Ceratium macrocerosA5 近缘斜纹藻
Pleurosigma affine表 3 养殖池塘浮游植物优势种冗余统计分析信息Table 3. Summary statistic for axes of RDA performed on dominant species in pond指标 Factors 轴 Axe 总方差
Total variance1 2 3 4 特征值 Eigenvalues 0.263 0.180 0.141 0.063 1.000 种类环境相关性 Species-environment correlations 0.939 0.922 0.963 0.632 累计变量百分比 Cumulative percentage variance 种类数据 Species data 26.3 44.3 58.4 64.6 种类环境相关 Species-environment relation 36.1 60.9 80.3 88.9 总特征值 Sum of all eigenvalues 1.000 所有典范特征值 Sum of all canonical eigenvalues 0.727 3. 讨论
3.1 浮游植物细胞密度的变化
本研究表明,养殖塘水体浮游植物细胞密度整体变化规律呈养殖前期较高,中后期较低的趋势,这与养殖塘水体营养盐含量的变化趋势一致 (图4)。营养盐,尤其是N、P含量是影响浮游植物密度的重要因子[21],摄食行为[22]和恶劣天气[23]也会对浮游植物细胞密度产生影响。Teeling等[24]研究表明,浮游植物细胞密度与水体内N、P含量成正相关关系,在一定范围内且其他影响因子一致的情况下,N、P含量越高,浮游植物细胞密度越高,这与本研究结果相一致。Burford和Williams[25]研究表明,水体中有机氮的主要来源是过剩的饵料和虾类排泄物,大量的有机氮在养殖水体内富集,影响水质情况进而影响对虾正常生长。养殖中后期,浮游植物细胞密度的变化可能与蟹类摄食病弱或死亡的对虾、粪便和残饵有关,蟹类的摄食行为会降低水体内有机氮的含量。李志斐等[26]研究发现鲮 (Cirrhinus molitorella) 和凡纳滨对虾混养后水体内的浮游植物细胞密度也呈降低趋势,富营养化程度得到缓解。在虾苗投放前,养殖池塘水体环境相对稳定,浮游植物群落自然演替,藻类群落结构具有相对稳定性。虾苗投放后,塘内的理化因子和浮游生物种群均受到影响,从而导致浮游植物细胞密度改变。对虾通过摄食一定数量的浮游动物,间接地影响了浮游植物的群落结构,幼虾还会滤食部分浮游植物,直接影响浮游植物群落[23]。养殖前期,对虾主要以天然饵料如硅藻浮游动物等为主要食物来源,养殖后期人工投放饵料占据虾食的98%[27]。本研究结果显示,养殖后期3个养殖塘浮游植物数量变化基本趋于平缓。查广才和周昌清[28]研究发现台风、暴雨、持续高温等恶劣天气会导致养殖水环境不稳定,导致浮游植物细胞密度发生较大波动。本研究中7—8月浮游植物细胞密度处于较低水平,这可能与华南地区7—8月强降水和台风等恶劣天气有关。
3.2 浮游植物种类的变化
本次研究发现,斑节对虾与锯缘青蟹混养养殖塘浮游植物种类共鉴定出5门29种,硅藻和蓝藻所占比重较大,其他藻类较少。郭永坚等[29]在凡纳滨对虾-鲻混养池塘的研究中鉴定出浮游植物7门65种,李志斐等[26]在鲮和凡纳滨对虾混养养殖塘研究中鉴定出浮游植物6门44种,不同水体浮游植物种类数的差异可能跟研究水域及采样时间差异有关[30-32]。本次调查结果显示,不同季节浮游植物种类呈现春秋季>夏季的规律,不同养殖阶段浮游植物优势种不同,养殖初期主要优势种以硅藻为主,养殖中后期优势种逐渐出现了蓝藻门。这可能与养殖后期盐度下降和pH升高有关。盐度是调节浮游植物群落结构的重要因子,养殖初期盐度较高,优势种以硅藻门为主,如柔弱菱形藻和直舟行藻等;养殖后期盐度下降,蓝藻门种类较多并成为优势种。这与吴斌和廖思明[33]在广西北海凡纳滨对虾养殖池塘浮游植物群落结构的研究结果一致,盐度较高的池塘优势种往往是绿藻门和硅藻门,盐度较低的池塘优势种是蓝藻门和甲藻门。谢立民等[34]也指出盐度是影响浮游植物群落结构的重要因子,蓝藻门在盐度较低的虾塘中为主要优势种,当盐度较高时,舟行藻等硅藻门为主要优势种。本次调查结果显示,pH介于7.78~9.21,水体一直处于弱碱性,而浮游植物的最适生长范围为7.75~8.75,过高或过低都会导致浮游植物生物量的减少。pH是由水体内理化因子、浮游生物和养殖物种综合作用的结果,周晴等[35]研究表明,pH高的水体中,蓝藻大量生长繁殖。养殖后期水体内pH上升,蓝藻成为虾塘优势种。浮游植物Shannon-Weaver多样指数介于0.802~2.140,不同季节多样性指数变化差异较大,春秋两季处于较高水平,夏季处于较低水平。养殖中后期,蓝藻大量繁殖成为优势种,因此降低了夏季浮游植物的多样性水平,武秀国等[36]和宋庆洋等[37]的研究中也揭示了这一规律。
4. 结论
在虾蟹混养养殖塘调查中,共发现浮游植物29种,其中硅藻门种类占绝对优势,其次是蓝藻门,甲藻门、绿藻门和裸藻门的种类数相对较少,浮游植物多样性存在较大的季节差异。浮游植物种类数总体为养殖初期较多,养殖中后期较少。浮游植物细胞密度在养殖初期 (3—5月) 较高,养殖中后期 (6—11月) 较低。
浮游植物主要的优势种主要为直舟形藻、奇异菱形藻、柔弱菱形藻、尖布纹藻、大角角藻和小颤藻等。浮游植物多样性指数均处于较高水平,由此可见其浮游植物群落结构较为稳定。
营养盐、盐度和pH是影响池塘浮游植物群落结构的主要环境因子。浮游植物的细胞密度、种类及分布都与环境因子的动态变化密切相关,环境因子的变动可能会改变养殖塘内浮游植物的群落结构。通过调控环境因子来改善浮游植物的群落结构,从而达到优化养殖环境,保持虾塘内生态系统平衡,增强抗干扰能力,进而提高产量和增加收益的目标。
-
图 3 3口池塘的α多样性系数分析
注:箱体上中下线分别为75、50 (中位数) 和25分位数,轴须线最长不超过1.5倍箱体范围,黑色空心圆表示平均数;差异显著性用* (P<0.05)、** (P<0.01)以及*** (P<0.001) 表示;图中的样本量:A:n=22、B:n=12、C:n=12。
Figure 3. α diversity index analysis of bacterioplankton in three ponds
Note: The upper, middle and lower lines of the box are 75, 50 (Median) and 25 quantiles, respectively. The maximum length of whiskers shall not exceed 1.5 times of the box range. The black hollow circles represent the average values. The significant differences were represented by * (P<0.05), ** (P<0.01) and *** (P<0.001). The numbers of replicated samples in this figure are: A: n=22; B: n=12; C: n=12.
图 7 优势菌群分布和群落多样性 (H') 与主要理化因子之间的关系 (II型标尺)
注:Aci. 不动杆菌属;Aer. 气单胞菌属;Fla. 黄杆菌属;Gem. 芽殖杆菌属;GpI. GpIIa;Ilu. 微酸菌属;Lim. 湖栖菌属;Pol. 多核杆菌;Rho. 红杆菌属;Sed. 沉积物杆状菌属;Spa. 发光细菌属;Sph. 鞘脂菌属。
Figure 7. Species associations of dominan flora and diversity (H') with environmental factors (Scaling II)
Note: Aci. Acinetobacter sp.; Aer. Aeromonas sp.; Fla. Flavobacterium sp.; Gem. Gemmobacter sp.; GpI. GpIIa; Ilu. Ilumatobacter sp.; Lim. Limnohabitans sp.; Pol. Polynucleobacter sp.; Rho. Rhodobacter sp.; Sed. Sediminibacterium sp.; Spa. Spartobacteria genera incertae sedis; Sph. Sphingomonas sp..
表 1 凡纳滨对虾养殖池塘水体主要理化因子
Table 1 Environmental factors in L. vannamei ponds
环境因子
Environmental factor池塘A
Pond A池塘B
Pond B池塘C
Pond C水温 Temperature/℃ 26.19±2.01 26.15±1.57 26.37±2.00 pH 8.10±0.50 8.29±0.37 8.30±0.38 溶解氧质量浓度 DO/(mg·L−1) 8.00±0.68 8.45±1.17 8.11±1.15 铵态氮质量浓度 NH4-N/(mg·L−1) 0.45±0.73 0.44±0.39 0.33±0.27 亚硝酸氮质量浓度 NO2-N/(mg·L−1) 0.01±0.01 0.01±0.01 0.01±0.01 硝酸态氮质量浓度 NO3-N/(mg·L−1) 0.34±0.44 0.41±0.42 0.48±0.53 活性磷质量浓度 PO4-P/(mg·L−1) 0.12±0.14 0.37±0.24 0.22±0.28 活性硅酸盐质量浓度 SiO3-Si/(mg·L−1) 11.2±3.55 12.38±7.38 4.03±2.98 高锰酸盐指数 CODMn/(mg·L−1) 8.34±7.19 6.12±3.18 5.34±2.91 叶绿素 a 质量浓度 Chl-a/(mg·L−1) 75.78±131.01 56.13±42.81 58.36±61.48 硫化物质量浓度 Sul/(mg·L−1) 0.03±0.04 0.02±0.01 0.02±0.01 矿化度质量浓度 Mineralization degree/(mg·L−1) 1209.44±2445.58 813.39±171.69 814.35±41.91 碱度 ALK/(mg·L−1) 138.18±12.06 115.32±33.82 84.46±14.73 总硬度 Total hardness/(mg·L−1) 362.76±395.28 476.5±227.34 349.36±60.18 表 2 各样本有效序列数据统计
Table 2 Valid sequences of each sample
样品
Sample条形码
Barcode有效序列
Valid sequence/条碱基数
Base number平均长度
Mean length/bp最短序列长度
Min. length/bp最长序列长度
Max. length/bpS1 GTAACA 82 736 34 604 966 418.26 363 469 S2 CCAGAC 81 905 34 619 470 422.68 361 476 S3 GGTGAA 56 859 23 683 959 416.54 367 475 S4 TGCATC 85 212 35 667 198 418.57 353 469 S5 TCGACC 83 246 34 547 936 415.01 365 470 S6 GTCGCG 73 829 30 547 465 413.76 362 453 S7 CGGATG 83 662 34 668 763 414.39 356 465 S8 GTGAAA 84 384 34 811 366 412.54 350 466 S9 ATCTTG 100 945 41 590 813 412.01 350 452 S10 TATGCA 73 801 30 659 241 415.43 352 459 S11 GTAACA 86 471 35 924 249 415.45 354 471 S12 GCGAGG 90 432 37 722 526 417.14 351 465 S13 CACGAT 54 439 22 735 194 417.63 373 465 S14 GCGGTA 44 227 18 478 585 417.81 352 471 S15 TATCGA 61 431 25 638 627 417.36 359 473 S16 ATCACG 52 384 21 792 829 416.02 376 471 S17 CGGATG 99 211 41 215 038 415.43 360 473 S18 CGCATA 100 001 41 427 908 414.27 372 476 S19 TGCATC 65 712 27 246 907 414.64 359 476 S20 TCAGTA 76 029 31 690 841 416.83 362 471 S21 CGGCAC 75 202 31 405 455 417.61 364 469 S22 ATCACG 71 605 30 328 221 423.55 352 470 S23 CGGATG 63 358 26 635 233 420.39 367 470 S24 GTGAAA 56 429 23 698 925 419.98 350 471 S25 TCAGTA 92 400 38 487 742 416.53 355 435 S26 GAAGTG 87 305 37 043 958 424.31 368 448 S27 TCGACC 94 075 39 472 386 419.58 360 464 S28 CTTGTA 52 273 22 054 548 421.91 373 472 S29 GTTTCG 44 697 18 684 447 418.02 360 462 S30 ATCTTG 59 737 25 072 871 419.72 365 465 S32 GCCATC 78 462 32 993 239 420.5 357 474 S33 TGTGTT 69 613 29 167 734 419.00 351 474 S34 CTTGTA 56 083 23 551 548 419.94 356 471 S35 GTTTCG 45 701 18 993 334 415.60 359 464 S36 TTCGTA 46 337 19 212 573 414.63 356 468 S37 CCAGAC 51 978 21 598 326 415.53 372 436 S38 AGCAGT 77 233 31 835 991 412.21 350 470 S39 GAGGAA 75 226 31 060 291 412.89 370 468 S40 AAGGTA 46 920 19 633 227 418.44 352 450 S41 ATCACG 43 983 18 242 559 414.76 352 469 S42 TAGGAC 66 639 27 702 230 415.71 356 470 S43 TGGACG 49 206 20 369 621 413.97 357 472 S44 AGAACA 50 322 20 765 128 412.65 356 469 S45 GGTGTG 41 169 17 011 770 413.22 350 471 S46 AACTAT 67 691 28 070 176 414.68 357 468 S47 ACTGCG 60 220 25 559 117 424.43 359 474 S48 TGTGTT 94 261 39 467 506 418.7 355 474 S49 TAGGAC 86 904 36 530 365 420.35 353 467 -
[1] CABELLO F C. Heavy use of prophylactic antibiotics in aquaculture: a growing problem for human and animal healthy and for the environment[J]. Environ Microbiol, 2006, 8: 1137-1144. doi: 10.1111/j.1462-2920.2006.01054.x
[2] MA Z, SONG X, WAN R, et al. A modified water quality index for intensive shrimp ponds of Litopenaeus vannamei[J]. Ecol Indic, 2013, 24: 287-293. doi: 10.1016/j.ecolind.2012.06.024
[3] GONZALEZ J M, SIMO R, MASSANA R, et al. Bacterial community structure associated with a dimethylsulfo niopropionate-producing North Atlantic algal bloom[J]. Appl Environ Microb, 2000, 66: 4237-4246. doi: 10.1128/AEM.66.10.4237-4246.2000
[4] WEI C L, ZENG Y H, TANG K X, et al. Comparison of bacterioplankton communities in three mariculture ponds farming different commercial animals in subtropical Chinese coast[J]. Hydrobiologia, 2009, 632: 107-126. doi: 10.1007/s10750-009-9831-8
[5] CARON D A. Inorganic nutrients, bacteria, and the microbial loop[J]. Microb Ecol, 1994, 28: 295-298. doi: 10.1007/BF00166820
[6] COTNER J B, BIDDANDA B. Small players, large role: microbial influence on biogeochemical processes in pelagic aquatic ecosystems[J]. Ecosystems, 2002, 5: 105-121. doi: 10.1007/s10021-001-0059-3
[7] LEMONNIER H, COURTIES C, MUGNIER C, et al. Nutrient and microbial dynamics in eutrophying shrimp ponds affected or unaffected by vibriosis[J]. Mar Pollut Bull 2010, 60: 402-411.
[8] AMANN R I. Phylogenetic identification and in situ detection of individual microbial cells without cultivation[J]. Microbiol Rev, 1995, 59: 143-169. doi: 10.1128/mr.59.1.143-169.1995
[9] LI C, LIU J, CHEN X, et al. Determinism governs the succession of disturbed bacterioplankton communities in a coastal maricultural ecosystem[J]. Sci Total Environ, 2022, 828: 154457. doi: 10.1016/j.scitotenv.2022.154457
[10] LIU C H, CHANG F H, YANG J W, et al. Free-living marine bacterioplankton composition and diversity along the Kuroshio region[J]. Deep-Sea Res I, 2022, 183: 103741. doi: 10.1016/j.dsr.2022.103741
[11] LU Z B, REN Z M, LIN W H, et al. Succession, sources, and assembly of bacterial community in the developing crab larval microbiome[J]. Aquaculture, 2021, 548: 737600.
[12] LIU Y P, PAN B Z, ZHU X Z, et al. Patterns of microbial communities and their relationships with water quality in a large-scale water transfer system[J]. J Environ Manage, 2022, 319: 115678. doi: 10.1016/j.jenvman.2022.115678
[13] 靳燕, 邱莹, 董志, 等. 北运河浮游细菌集合群落空间变化的环境解释[J]. 中国环境科学, 2021, 41(3): 1378-1386. doi: 10.3969/j.issn.1000-6923.2021.03.041 [14] 杨淑芳. 凡纳滨对虾养殖池塘水质动态变化及细菌多样性研究[D]. 厦门: 集美大学, 2017: 10-19. [15] 张皓. 养殖环境微生物群落的动态变化及与水环境的互作[D]. 苏州: 苏州大学, 2015: 16-29. [16] YANG W, ZHENG C, ZHENG Z M, et al. Nutrient enrichment during shrimp cultivation alters bacterioplankton assemblies and destroys community stability[J]. Ecotox Environ Safe, 2018, 156: 366-374. doi: 10.1016/j.ecoenv.2018.03.043
[17] YANG W Z, ZHENG J Y, CHEN Q, et al. Succession of bacterioplankton community in intensive shrimp (Litopenaeus vannamei) aquaculture systems[J]. Aquaculture, 2018, 497: 200-213. doi: 10.1016/j.aquaculture.2018.07.053
[18] 国家环境保护总局水和废水检测分析方法编委会. 水和废水检测分析方法(第四版)[M]. 北京: 中国环境科学出版社, 2002: 30-152. [19] EDGAR R C. Search and clustering orders of magnitude faster than BLAST[J]. Bioinformatics, 2010, 26(19): 2460-2461. doi: 10.1093/bioinformatics/btq461
[20] AKSNES D L, WASSMANN P. Modeling the significance of zooplankton grazing for export production[J]. Limnol Oceanogr, 1993, 38(5): 978-985. doi: 10.4319/lo.1993.38.5.0978
[21] FAITH D P. Conservation evaluation and phylogenetic diversity[J]. Biol Conserv, 1992, 61: 1-10. doi: 10.1016/0006-3207(92)91201-3
[22] SHANNON C E, WEAVER W. The mathematical theory of communication[M]. Urbana, IL: University of Illinois Press, 1949: 55-187.
[23] 杨家新. 微生物生态学[M]. 北京: 化学工业出版社, 2004: 45-134. [24] HAN J H, BAEK K, LEE M H. Lacibacter nakdongensis sp. nov., isolated from river sediment[J]. Int J Syst Evol Microbiol, 2017, 67: 352-356. doi: 10.1099/ijsem.0.001631
[25] WITCZAK P, CYDZIK-KWIATKOWSKA A, ZIELIŃSKA M. Treatment of the liquid phase of digestate from a biogas plant for water reuse[J]. Bio Tech, 2019, 276: 226-235. doi: 10.1016/j.biortech.2018.12.077
[26] XUAN L, SHENG Z, LU J, et al. Bacterioplankton community responses and the potential ecological thresholds along disturbance gradients[J]. Sci Total Environ, 2019, 696: 134015. doi: 10.1016/j.scitotenv.2019.134015
[27] YANG Y, GAO Y, HUANG X, et al. Adaptive shifts of bacterioplankton communities in response to nitrogen enrichment in a highly polluted river[J]. Environ Pollut, 2019, 245: 290-299. doi: 10.1016/j.envpol.2018.11.002
[28] COMTE J, del GIORGIO P A. Composition influences the pathway but not the outcome of the metabolic response of bacterioplankton to resource shifts[J]. PLoS One, 2011, 6: e25266. doi: 10.1371/journal.pone.0025266
[29] ADAMOVICH B V, MEDVINSKY A B, NIKITINA L V, et al. Relations between variations in the lake bacterioplankton abundance and the lake trophic state: evidence from the 20-year monitoring[J]. Ecol Indic, 2019, 7: 120-129.
[30] 吕锡斌, 吴云成, 陈良强, 等. 赤水河流域浮游细菌群落特征及其与水质的关系[J]. 环境科学学报, 2021, 41(11): 4596-4605. doi: 10.13671/j.hjkxxb.2021.0125 [31] JING X Y, GOU H, GONG Y, et al. Seasonal dynamics of the coastal bacterioplankton at intensive fish-farming areas of the Yellow Sea, China revealed by high-throughput sequencing[J]. Mar Pollut Bull, 2019, 139: 366-375. doi: 10.1016/j.marpolbul.2018.12.052
[32] SHADE A, HANNES P, STEVEN D A, et al. Fundamentals of microbial community resistance and resilience[J]. Front Microbiol, 2012, 19(3): 417.
[33] 郭远涛. 吸附固定化对藻菌共生体在畜禽养殖废水处理中的强化效果研究[D]. 南昌: 南昌大学, 2021: 12-36. [34] CAVENDER-BARES J, KOZAK K H, FINE P V A, et al. The merging of community ecology and phylogenetic biology[J]. Ecol Lett, 2009, 12: 693-715. doi: 10.1111/j.1461-0248.2009.01314.x
[35] CHESSON P. Mechanisms of maintenance of species diversity[J]. Annu Rev Ecol Syst, 2000, 31: 343-366. doi: 10.1146/annurev.ecolsys.31.1.343
[36] FAN L M, HU G D, QIU L P, et al. Variations in bacterioplankton communities in aquaculture ponds and the influencing factors during the peak period of culture[J]. Environ Pollut, 2020, 258: 113656. doi: 10.1016/j.envpol.2019.113656
[37] IVES A R, CARPENTER S R. Stability and diversity of ecosystems[J]. Science, 2007, 317: 58-62. doi: 10.1126/science.1133258
[38] YACHI S, LOREAU M. Biodiversity and ecosystem productivity in a fluctuating environment: the insurance hypothesis[J]. Proc Natl Acad Sci USA, 1999, 96: 1463-1468. doi: 10.1073/pnas.96.4.1463
[39] 胡晓娟, 李卓佳, 曹煜成, 等. 强天气干扰条件下粤西凡纳滨对虾养殖池塘细菌群落动态特征[J]. 南方水产科学, 2012, 8(5): 52-59. doi: 10.3969/j.issn.2095-0780.2012.05.008 [40] CONNELL J H. Diversity in tropical rain forests and coral reefs: high diversity of trees and corals is maintained only in a non-equilibrium state[J]. Science, 1978, 199: 1302-1310. doi: 10.1126/science.199.4335.1302
-
期刊类型引用(8)
1. 傅建军,安睿,朱文彬,王兰梅,罗明坤,董在杰. 套养池塘主要养殖阶段浮游生物的群落结构特征. 水产学杂志. 2024(01): 104-112 . 百度学术
2. 宋光同,王芬,徐笑娜,朱成骏,陈祝,李翔,周翔,叶圣陶,蒋业林. 中华鳖养殖过程中浮游生物群落结构的演变. 中国农学通报. 2024(08): 157-164 . 百度学术
3. 季雨来,吴召仕,张京,刘鹏,裴鹏娜,张帅,黄佳聪. 太湖流域平原农业区典型圩塘浮游植物群落特征及其影响因子. 湖泊科学. 2024(05): 1380-1392 . 百度学术
4. 闵文武,王龙燕,陈飞雄,周其椿. 集装箱养殖生态净化池塘中浮游植物群落结构多样性分析. 水产科技情报. 2024(05): 302-310 . 百度学术
5. 蒲炜佳,董世鹏,张东旭,于力业,谢宜成,徐羡,王芳,李由明. 三疣梭子蟹池塘综合养殖系统浮游植物群落结构及其与环境因子的关系. 中国水产科学. 2022(04): 549-561 . 百度学术
6. 侯德昌,张莹莹,魏文志. 不同中华鳖养殖模式浮游植物功能群特征及水环境评价. 安徽农业科学. 2022(09): 96-99+170 . 百度学术
7. 马景雪,张培玉,王宗兴,郑明刚,高萍,曲凌云,王波,郑风荣. 青岛崂山湾近海扇贝养殖区细菌多样性及环境因子分析. 海洋科学进展. 2022(02): 307-319 . 百度学术
8. 张晓蕾,王强,张国奇,周陆,李廷发,张玉,赵思雅. 池塘循环流水养殖模式中浮游植物群落结构的空间变化研究. 南方水产科学. 2021(03): 36-45 . 本站查看
其他类型引用(2)