裴泓霖 王 荦 谭八梅 李伟嘉 甄 昊 丁 君 常亚青

(大连海洋大学 农业农村部北方海水增养殖重点实验室 大连 116023)

刺参(Apostichopus j aponicus)因具有丰富的营养和药用价值，深受消费者青睐。随着刺参养殖产业的不断发展，池塘养殖已成为刺参增养殖的重要方式(常亚青等, 2006)。在集约化养殖过程中，残饵及养殖生物粪便的堆积易引起底质恶化，对养殖生态系统健康存在潜在危害(董双林, 2015; 王姣, 2016)。

养殖池塘沉积物菌群作为养殖生态系统的重要组成部分，不仅在水体物质循环及能量流动过程中发挥巨大的作用(Bentzon-Tilia et al, 2016)，对其养殖生物的生长、消化、免疫、疾病等均存在相关性(Becker et al, 2004; Wang et al, 2005; 高菲等, 2010)。研究表明，蓝细菌(Cyanobacteria)能将溶解有机碳转化为自身的颗粒有机碳，从而完成碳转换与储存(周伟华等,2007)；假单胞菌属(Pseudomonas)菌群能利用反硝化作用将硝酸盐还原为氮气(肖晶晶等, 2009)；链霉菌属(Streptomyces)、芽孢杆菌属(Bacillus)、节杆菌属(Arthrobacter)中某些菌种对有机磷和无机磷的分解具有明显作用(孙亚凯, 2006)；脱氮硫杆菌(Thiobacillus denitrificans ) 、 反 硝 化 硫 微 螺 菌(Thiomicrospira deni trificans) 、 多 样 硫 杆 菌(Thiobacillusd versutus )等菌群能从多种还原态无机硫的氧化中获得能量，同时，还原硝酸盐或亚硝酸盐产生氮气(Pokorna et al, 2015)；嘉利翁氏菌属(Gallionella)、海杆菌属(Mariprofundus)、Sideroxydans 能够将Fe2+氧化为Fe3+，而Fe3+可以作为电子受体进行能量代谢，从而完成铁循环(Edwards et al, 2003; Emerson et al,2004)。因此，有效利用某些微生物的生理功能来保持养殖池塘生态系统的微生态平衡，对防治池塘养殖环境恶化以及促进养殖生物健康生长具有重要意义(王鑫毅等, 2019)。然而，养殖池塘生态系统菌群结构极易受到外部环境理化因子的影响，季节变动、营养物质输入等变化均可造成池塘养殖环境菌群种类和丰度的改变，间接影响养殖生物的生长及健康(Sunde et al, 2001; Suzuki et al, 2019)。

为探明不同季节刺参养殖池塘沉积物菌群结构及其影响因素，本研究选取辽宁省营口地区刺参养殖池塘为研究对象，基于高通量测序技术，揭示不同季节刺参养殖池塘沉积物菌群结构的季节性差异和共性，查明影响刺参养殖池塘沉积物菌群结构的主导环境因子，以期为刺参养殖池塘微生态调控和绿色健康养殖提供一定的理论依据。

1 材料与方法

1.1 样品采集

选取辽宁营口(40°37′46″ N；122°8′59″ E)刺参养殖池塘为研究对象，于2019 年1 月(冬季)、4 月(春季)、8 月(夏季)、10 月(秋季)，利用柱状采泥器，在同一刺参养殖池塘对角线两侧离岸1.5 m 及正中间处采集表层(5 cm)沉积物样品，装入无菌离心管，-80℃冷冻保存。将每次采集的5 个沉积物样品等量混匀送检，每个样品进行3 次重复。

采样池塘仅养殖刺参单一品种，放养刺参密度约为8～10 头/m2，刺参规格约为100 g/头。池塘面积约为0.1 km2，水深2～3 m，底质多泥质，水源来自近海海域，利用潮汐带来的潮差进水和排水。全年不投喂饲料及益生菌，保证刺参自然生长。

1.2 环境因子测定

环境因子监测点位同沉积物采样点，每个采样点测定环境指标均进行3 次重复。采用HACH(HQd)水质分析仪器，在养殖池塘现场测定温度、盐度和pH。沉积物总有机碳(TOC)的测定采用重铬酸钾氧化-还原容量法，总氮(TN)采用过硫酸钾氧化法，总磷(TP)采用钼锑抗分光光度法，氨氮(-N)采用次溴酸盐氧化法，硝酸盐(-N)采用镉柱还原法，亚硝酸盐(N)采用萘乙二胺分光光度法，磷酸盐(-P)采用磷钼蓝分光光度法，硫化物采用亚甲基蓝分光光度法，氧化还原电位采用碱性高锰酸钾法。

1.3 DNA 提取

选用OMEGA Soil DNA 试剂盒提取刺参养殖池塘环境样品总DNA，为确保DNA 无降解无污染，DNA 浓度要求≥50 ng/μl，OD260/280=1.8～2.0。利用1%琼脂糖凝胶电泳检测DNA，条带清晰完整、无拖尾，质量合格。

1.4 PCR 扩增与高通量测序

细菌16S rRNA 基因V3～V4片段的扩增引物为341F(5′-CCTACGGGNGGCWGCAG-3′)和805R(5′-G ACTACHVGGGTATCTAATCC-3′)，应用Illumina MiSeq测序平台完成序列测序。

1.5 数据分析

采用Usearch 软件，以97%阈值开展OTU(Optical Transform Unit)操作分类单元聚类(Blaxter et al ,2005)。选用RDP classifier 软件(Wang et al, 2007)进行物种分类操作，阈值设置为0.8，低于该阈值的分类结果被划归为Unclassified 一类。利用Qiime 软件(Caporaso et al, 2010)进行多样性分析。基于OTU 聚类和环境因子信息，利用R 语言的Vegan package(Oksanen et al, 2010)进行环境因子相关性分析。通过最大似然法构建进化树。根据物种分类结果得到不同水平上各Rank 的丰度值，采用STAMP 软件(版本2.1.3)，比较样本之间丰度差异，P＜0.05 为差异显著，比较样本组间丰度差异。

2 结果与分析

2.1 环境因子监测

辽宁营口刺参养殖池塘不同季节沉积物环境因子监测数据如表1 所示。根据中华人民共和国国家标准《海洋沉积物质量(GB 18668-2002)》，所有样品均符合国家Ⅰ类沉积物标准。所监测的环境因子理化指标，均呈显著的季节性差异。其中，总氮、氨氮、硝酸盐、总磷和硫化物含量均在夏季最高；TOC 含量在秋季最高，春季最低。

2.2 高通量测序结果

不同季节刺参养殖池塘沉积物样品高通量测序结果见表 2，从表 2 可以看出，共得到原始序列37690～87024 条，有效序列41562～79288 条，聚类于36980～3913 个 OTUs。样品测序覆盖率为 95.9%～99.6%，表明测序结果可以真实反映样本信息。稀释曲线趋于平坦，表示测序数据深度足够反映样品的多样性。

2.3 α 多样性分析

Chao1 和Ace 指数常用于估算样品中菌群的丰度；Shannon 和Simpson 指数常用于估算样品中菌群的多样性。

表1 不同季节刺参养殖池塘沉积物环境因子Tab.1 Environmental factors of sediment in Apostichopus japonicus culture pond in different seasons

表2 高通量测序结果Tab.2 Results of high-throughput sequencing

比较不同季节α 多样性指数得出(图1)，营口地区刺参养殖池塘沉积物菌群丰度和多样性具有显著性季节差异。其中，夏季菌群丰度和多样性显著升高；其他季节菌群丰度和多样性均有不同程度的下降，具体表现为秋冬两季菌群丰度极显著性下降(P＜0.01)，且冬季菌群物种丰度达四季最低。

2.4 β 多样性分析

在OTU 水平上，根据加权Unifrac 距离对不同季节样本进行三维主坐标分析(图 2)。其中，PCoA1(Principal Co-ordinates analysis)的贡献率为55%，PCoA2的贡献率为35%，PCoA3 的贡献率为4%，总贡献率为94%。同一季节样品均聚集在一起，显示出较好的生物学重复性。不同季节相比较，秋季与其他季节沉积物菌群结构具有显著差异(P＜0.05)。

图1 α 多样性指数Fig.1 Alpha diversity index

图2 PCoA 三维分析Fig.2 PCoA three-dimensional analysis

2.5 菌群结构特征

基于GraPhlAn 分类和系统发育信息可视化图(图3)可以看出，不同季节刺参养殖池塘沉积物第一优势菌门均为变形菌门(Proteobacteria)(相对丰度＞43.19%)，次优势菌门存在显著差异。其中，绿弯菌门(Chloroflexi)(相对丰度＞7%)为春季次优势菌门；拟杆菌门(Bacteroidetes)(相对丰度＞8.13%)为夏季次优势菌门，厚壁菌门(Firmicutes)(相对丰度＞17.96%)为秋、冬季次优势菌门。

菌属水平下，优势菌属呈现出显著的季节性差异(P＜0.05)。其中，春季优势菌主要隶属于Sulfurovum、Thioprofundum、除硫单胞菌属(Desulfuromonas)、Desulfosarcina 和Desulfopila。夏季优势菌属主要隶属于Sulfurovum、Halioglobus、Haliea、Thioprofundum、Desulfuromonas、Desulfopila 和Desulfobulbus。秋季优势菌主要隶属于Sulfurovum、Desulfuromonas、Desulfuromusa、Desulfosarcina、Desulfopila 和海杆菌属(Marinobacter)。冬季优势菌主要隶属于柠檬酸杆菌属(Citrobacter)、假单胞菌属(Pseudomonas)、Desulfuromonas、Pelobacter、Desulfosarcina、Desulfopila、Desulfobulbus、不动杆菌属(Acinetobacter)和Exiguobacterium。春、夏、秋、冬 4 季沉积物优势菌属均包含脱硫八叠菌Desulfosarcina、Desulfuromonas 和Desulfopila。

LEfSe (LDA Effect Size)差异性分析结果显示(图4)，共有不同分类水平的98 个细菌群类在4 个季节中存在显著差异。其中，变形菌门在4 个季节中都存在差异菌群，如柠檬酸杆菌属、不动杆菌属、假单胞菌属是隶属于春季的差异菌属；外硫红螺旋菌科(Ectothiorhodospiraceae)、交替单胞菌目(Alteromonadales)是隶属于夏季的差异菌目； 弯曲菌目(Campylobacterales)、Sulfurovum 是隶属于秋季的差异菌；着色菌目Chromatiales、Thiohalomonas 是隶属于冬季的差异菌目。此外，春季刺参养殖池塘沉积物差异菌门主要隶属于拟杆菌门(Bacteroides)和酸杆菌门(Acidobacteria)，其代表性差异菌属分别为Cytophagla和Gp21。夏季沉积物差异菌门主要隶属于酸杆菌门、疣微菌门(V e r r u c o m i c r o b i a)和浮霉菌门(P l a n c t o m y c e t e s)，分别以 G p 2 3、疣微菌科(Verrucomicrobiaceae)和浮霉菌科(Planctomycetaceae)为代表。秋季沉积物差异菌门隶属于拟杆菌门和厚壁菌门，其代表性差异菌属为 Mangrovibacterium 和Natranaerovirga。冬季沉积物差异菌门则隶属于放线菌门(Actinobacteria)以Actinomycetales 为代表。

2.6 菌群结构与环境因子相关性分析

RDA 是将对应分析与多元回归分析混合的一种分析方法。营口地区刺参养殖池塘沉积物菌群与环境因子RDA(Redundancy Analysis)分析结果显示(图5)，春季和夏季在第2 象限中基本处于同一位置，显示影响春夏两季沉积物菌群结构的环境因子相似。根据环境因子与原点连线的长度可知，温度、TN、pH 和TOC 是营口刺参养殖池塘沉积物菌群的主导环境因子。其中，温度、TN、pH 与春、夏季刺参养殖池塘沉积物菌群相对丰度呈正相关，与秋、冬季菌群相对丰度呈负相关；TOC 与秋、冬季刺参养殖池塘菌群相对丰度呈正相关，与春、夏季菌群相对丰度呈负相关。

图4 不同季节沉积物菌群LEfSe 分析Fig.4 LEfSe analysis of sediment bacterial community in different seasons

图5 沉积物菌群与环境因子RDA 分析Fig.5 RDA analysis of sediment bacterial community and environmental factors

不同季节刺参养殖池塘沉积物差异菌群与主导环境因子相关性分析结果见表3。其中，拟杆菌门与TOC 含量呈极显著负相关(P＜0.01)，表明拟杆菌门受沉积物有机物质含量影响显著。酸杆菌门与pH 和TN 呈极显著正相关(P＜0.01)。疣微菌门和浮霉菌门与TN 含量呈极显著正相关(P＜0.01)。厚壁菌门与TN 呈极显著负相关(P＜0.01)，与TOC 含量呈极显著正相关(P＜0.01)。放线菌门则与pH 和TOC 呈显著相关性(P＜0.05)。

表3 沉积物优势菌群与主导环境因子相关性分析Tab.3 Correlation analysis between dominant bacterial community and environment factors in sediment

3 讨论

3.1 不同季节刺参养殖池塘沉积物菌群多样性

微生物是养殖池塘生态系统的重要组成部分，在生态系统物质循环和能量流动过程中发挥着重要作用，对维持生态系统平衡具有重要意义(Rungrassamee et al, 2016; Xiong et al, 2016; 金笑等, 2017)。本研究结果显示，营口地区不同季节刺参养殖池塘沉积物菌群丰度和多样性呈现出显著的季节演替特征。其中，夏季刺参养殖池塘沉积物菌群丰度和多样性最高，冬季沉积物菌群丰度和多样性则最低。已有的研究表明，温度是影响环境微生物生命活动的重要因素之一，高温通过加速细菌体内酶活性，促进微生物代谢，加速有机物分解，为微生物生长繁殖提供有利条件(李彬等, 2010; Cardona et al, 2016)。同时，氮、磷等营养物质的增多，也为细菌生长创造了适宜的条件，使菌群丰度和多样性更丰富(董学兴等, 2019)。营口地区夏季刺参养殖池塘水环境温度和氮磷含量相对较高，为沉积物中微生物的繁殖提供了良好的生存环境，因此，夏季沉积物菌群丰度和多样性显著高于其他季节。

3.2 不同季节刺参养殖池塘沉积物菌群结构及其影响因素

本研究中，不同季节营口地区刺参养殖池塘沉积物第1 优势菌门均为变形菌门，与已有的研究报道一致(窦妍等, 2014、2016; 丁斯予等, 2019)。变形菌门作为养殖池塘环境中的第1 优势菌门，主要是由其在原核生物分子生物学分类中的绝对优势所决定的(Gupta et al, 2000)。变形菌门下多个细菌群类均为池塘底栖环境中常见的固氮菌，参与池塘生态系统氮循环，在池塘养殖生态系统中占据着重要的生态学地位(刘思亮等, 2002; 龚骏等, 2010)。刺参养殖池塘沉积物除变形菌门为第1 优势菌门外，其他优势菌门为拟杆菌门、酸杆菌门、疣微菌门、厚壁菌门等，与营底栖生活的缢蛏(Sinonovacula const ricta)相近(王毅鑫等, 2019)，与杂食性的罗非鱼池塘沉积物菌群差别较大(史丽娜等, 2015)。养殖品种食性的差异会导致饲料投入不同，造成池塘沉积物中有机物和营养盐浓度产生差异。除饲料不同之外，刺参的摄食方式为吞食沉积物，造成沉积物菌群结构发生一定改变(于东祥,2010)。因此，刺参养殖池塘沉积物菌群结构具有特异性。

比较不同季节刺参养殖池塘沉积物菌群结构特征的结果显示，刺参养殖池塘沉积物菌群存在显著的季节性差异，特异性菌群相对丰度与其生境因子之间存在显著性相关，影响沉积物菌群的主导环境因子为温度、TOC、TN 和pH，且所有主导环境因子均呈现出显著的地域差异。已有的研究表明，微生物生长繁殖所需的营养物质主要依靠环境中的无机盐和有机物，其群落组成对环境因子变化有着较强的敏感性(王轶南等, 2010; Bell et al, 2013; 李建光等, 2014;Yang et al, 2015; Amabebe et al, 2020)。微生物菌群特征是其生境的反映，与环境理化性质密切相关且相互限制影响(闫法军等, 2014)。其中，春季刺参养殖池塘沉积物特异性菌群主要以拟杆菌门为代表。拟杆菌门细菌与微生物蛋白质和脂类等有机物的转化密切相关(Ren et al, 2015)，是维持环境碳循环的主要动力，且与春季沉积物环境中TOC 含量显著性相关。夏季沉积物特异性菌群以酸杆菌门、疣微菌门和浮霉菌门为代表。夏季随着温度的升高，养殖生物代谢加快，排泄物增多，易导致硫、铁等元素在池底大量堆积，使养殖池塘底部沉积物pH 下降，为具有嗜酸特性的酸杆菌门繁殖提供了有利条件(Barns et al, 2007;刘彩霞等, 2016)。此外，夏季沉积物中含氮物质的增高会造成水体富营养化，导致蓝藻暴发，也会对水生生物的生长产生不良影响(付莹等, 2018)。然而，疣微菌门细菌的存在可以有效利用夏季沉积物环境中的氮元素，加速氮循环进程(Shen et al, 2017)。夏季沉积环境中高含量的氮和硫，则有利于浮霉菌门细菌的生长繁殖(黄佩蓓等, 2014)。厚壁菌门细菌能够有效利用环境中硝酸盐和亚硝酸盐，进而降低沉积物中的氮含量，作为秋季特异性菌群，与其沉积物中较低的氮含量呈现出一致性。冬季随着气温的下降，沉积物中微生物丰度和多样性均受到不同程度的抑制，而放线菌门细菌因能形成抗逆性极强的芽孢，从而具有极强的环境适应性(王蓉等, 2013)，成为冬季沉积物特异性菌群。此外，放线菌门具有降解淀粉、蛋白等大分子的功能，产生的抗菌素类物质能抑制环境中其他革兰氏阳性细菌发挥作用，有成为益生菌的潜力(李文均等, 2002)。

值得注意的是，夏季刺参养殖池塘沉积物中厚壁菌门相对丰度极显著低于其他季节，相对丰度仅为0.5%。而隶属于厚壁菌门、芽孢杆菌目(Bacillales)的芽孢杆菌属是典型的益生菌，被广泛用于水产养殖池塘环境微生态调控(杜佗等, 2017)。芽孢杆菌可以产生多种消化酶和营养物质，供养殖动物吸收利用，还可以通过改善水产动物体内微生态环境进而提高动物机体免疫力，同时，具有净化水质和清除自由基的功效(苏艳莉等, 2016)。因此，建议在夏季池塘中添加适当含量的芽孢杆菌，用于协助刺参安全度夏。

4 结论

本研究结果表明，营口地区刺参养殖池塘沉积物菌群丰度和多样性存在显著性季节变化特征，温度、总有机碳、总氮和pH 是影响沉积物菌群的主要环境因子。本研究初步探明沉积物菌群与主导环境因子间的相关性，而菌群与环境因子间的互作机制还有待于进一步深入探究。